mod_JHUGen.F90

Go to the documentation of this file.

MODULE modjhugen
   use modparameters
   use modjhugenmela
   use modkinematics
#if compiler==1
   use ifport
#endif
   implicit none
 
   public :: initfirsttime,resetpdfs
 
   CONTAINS
 
 
 
SUBROUTINE initfirsttime(pdftable,pdfstrlength,pdfmember,collider_sqrts)
   use modparameters
   use modkinematics
   use modmisc
   implicit none
   integer :: pdfstrlength
   character(len=pdfstrlength) pdftable
   integer :: pdfmember
   real(8), intent(in) :: collider_sqrts
 
   call setjhugendefaults()
 
   includeinterference=.true.
   includegammastar=.true.
   mphotoncutoff=4d0*gev
   widthscheme=0
   pdfset=3      ! 1: CTEQ6L1   2: MRSW with best fit, 2xx: MSTW with eigenvector set xx=01..40
   lhapdfstring = pdftable
   lhapdfmember = pdfmember
   lenlhapdfstring = pdfstrlength
 
!---------------------------
   call printlogo(io_stdout, "JHUGen MELA")
!---------------------------
   call resetmubarhgabarh()
!---------------------------
   call computeewvariables()
   call computeckmelements(vckm_ud, vckm_us, vckm_cd, vckm_cs, vckm_ts, vckm_tb)
   print *, "JHUGen CKM initialization"
   print *, "Vud = ",vckm_ud
   print *, "Vus = ",vckm_us
   print *, "Vub = ",vckm_ub
   print *, "Vcd = ",vckm_cd
   print *, "Vcs = ",vckm_cs
   print *, "Vcb = ",vckm_cb
   print *, "Vtd = ",vckm_td
   print *, "Vts = ",vckm_ts
   print *, "Vtb = ",vckm_tb
!---------------------------
 
#if useLHAPDF==1
   if( lhapdfstring.eq."" ) then
      print *, "Need to specify pdf file for LHAPDF"
      stop
   endif
#else
   if( lhapdfstring.eq."" .and. pdfset.eq.3) then
      print *, "Need to specify pdf file for NNPDF"
      stop
   endif
#endif
 
   if( isazdecay(decaymode1) ) then
      m_v = m_z
      ga_v= ga_z
   elseif( isawdecay(decaymode1) ) then
      m_v = m_w
      ga_v= ga_w
   elseif( isaphoton(decaymode1) ) then
      m_v = 0d0
      ga_v= 0d0
   endif
 
   call initpdfs()
   IF( collider.EQ.1) THEN
      collider_energy  = collider_sqrts
   ELSE
      print *,"Collider not implemented."
      stop
   ENDIF
 
#if useCollier==1
   collier_maxnloopprops = -1
   collier_maxrank = -1
#endif
   call initcollier(4,3) ! Arguments for ggZH
 
   return
 
END SUBROUTINE
 
SUBROUTINE setjhugendefaults()
   use modparameters
   use modhashcollection
   implicit none
 
   call setdefaultckm()
 
   collider=1
   vegasit1=-1
   vegasnc0=-1
   vegasnc1=-1
   vegasnc2=-1
   process = 0   ! select 0, 1 or 2 to represent the spin of the resonance
   pchannel=2
   decaymode1=0  ! Z/W+
   decaymode2=0  ! Z/W-
! !       DecayMode=0:  Z --> l+ l- (l=e,mu)
! !       DecayMode=1:  Z --> q qbar (q=u,d,c,s,b)
! !       DecayMode=2:  Z --> tau+ tau-
! !       DecayMode=3:  Z --> nu nubar (nu=nu_e,nu_mu,nu_tau)
! !       DecayMode=4:  W --> l nu_l (l=e,mu)
! !       DecayMode=5:  W --> q qbar' (q=u,c, qbar'=d,s)
! !       DecayMode=6:  W --> tau nu_tau
! !       DecayMode=7:  photon
! !       DecayMode=8:  Z --> l+ l- (l=e,mu,tau)
! !       DecayMode=9:  Z --> anything
! !       DecayMode=10: W --> l nu_l (l=e,mu,tau)
! !       DecayMode=11: W --> anything
   topdecays=-1
   taudecays=-1
   h_dk = .false.
   unweighted =.true.
   mufacmultiplier = 1d0
   murenmultiplier = 1d0
   facscheme = krenfacscheme_default
   renscheme = krenfacscheme_default
   offshellreson=.true.
   offshellv1=.true.
   offshellv2=.true.
 
   lheprodfile=""
   readlhefile=.false.
   convertlhefile=.false.
   readcsmax=.false.
   readpmzz = .false.
   pmzzfile="PMZZdistribution.out"
   pmzzevals=-1
   doprintpmzz = .false.
   printpmzzintervals = 20
   generateevents=.false.
   requestnleptons = -1
   requestos=-1
   requestossf=-1
   counttauasany = .true.
   writefailedevents=0
 
   call setuphashes()
 
END SUBROUTINE setjhugendefaults
 
SUBROUTINE initpdfvalues()
   use modparameters
   use modkinematics
   implicit none
 
#if useLHAPDF==0
   IF (alphas_mz .LE. 0d0) THEN
      WRITE(6,*) .le.'alphas_mz  0:',alphas_mz
      WRITE(6,*) 'continuing with alphas_mz=0.118'
      alphas_mz=0.118d0
   ENDIF
#endif
 
   mu_fact = m_reso ! Set pdf scale to resonance mass by default, later changed as necessary in the EvalWeighted/EvalUnweighted subroutines
   mu_ren = m_reso ! Set renorm. scale to resonance mass by default, later changed as necessary in the EvalWeighted/EvalUnweighted subroutines
   call evalalphas() ! Set alphas at default Mu_Ren. Notice ModParameters::ComputeQCDVariables is automatically called!
   return
END SUBROUTINE
 
SUBROUTINE initpdfs()
 
   use modparameters
   use modkinematics
   implicit none
 
#if useLHAPDF==1
 
   implicit none
   DOUBLE PRECISION alphasPDF
 
     call initpdfsetbyname(trim(lhapdfstring,lenlhapdfstring)) ! Let LHAPDF handle everything
     call initpdf(lhapdfmember)
 
     alphas_mz=alphaspdf(zmass_pdf)
     ! Dummy initialization, just in case. These values are not used.
     !nloops_pdf = 1
     zmass_pdf = m_z
 
#else
 
     zmass_pdf = m_z ! Take zmass_pdf=M_Z in pdfs that do not specify this value
 
     if( pdfset.eq.1 ) then ! CTEQ6L1
        call setctq6(4)  ! 4    CTEQ6L1  Leading Order cteq6l1.tbl
 
        alphas_mz=0.130d0
        !nloops_pdf=1
     elseif( pdfset.eq.3 ) then  ! NNPDF 3.0 LO with a_s=0.13
        call nnpdfdriver(lhapdfstring,lenlhapdfstring)
        call nninitpdf(lhapdfmember)
 
        alphas_mz=0.130d0
        !nloops_pdf=1
        zmass_pdf=91.199996948242188d0*gev
     elseif( (pdfset.eq.2) .or. (pdfset.ge.201 .and. pdfset.le.240) ) then ! MSTW2008 and variations
        alphas_mz=0.13939d0
        !nloops_pdf=1
     else ! Everything else
        write(6,*) "mod_JHUGen.F90::InitPDFs: PDFSet",pdfset,"QCD parameters are unknown. Please double-check! Stopping JHUGen..."
        stop
        ! Could also have used these instead of the stop statement, but why introduce arbitrary number?
        !alphas_mz = 0.13229060d0
        !nloops_pdf = 1
     endif
 
#endif
 
     call initpdfvalues() ! Call this only once
   return
END SUBROUTINE
 
subroutine resetpdfs(pdftable,pdfstrlength,pdfmember,pdfid)
implicit none
integer :: pdfstrlength
character(len=pdfstrlength) pdftable
integer :: pdfmember
integer, intent(in) :: pdfid
   pdfset=pdfid      ! 1: CTEQ6L1   2: MRSW with best fit, 2xx: MSTW with eigenvector set xx=01..40    3: NNPDF3.0
   lhapdfstring = pdftable
   lhapdfmember = pdfmember
   lenlhapdfstring = pdfstrlength
#if useLHAPDF==1
   if( lhapdfstring.eq."" ) then
      print *, "Need to specify pdf file for LHAPDF"
      stop
   endif
#else
   if( lhapdfstring.eq."" .and. pdfset.eq.3) then
      print *, "Need to specify pdf file for NNPDF"
      stop
   endif
#endif
   call initpdfs()
end subroutine
 
 
END MODULE modjhugen