JHUGen MELA  v2.4.1
Matrix element calculations as used in JHUGen. MELA is an important tool that was used for the Higgs boson discovery and for precise measurements of its structure and interactions. Please see the website https://spin.pha.jhu.edu/ and papers cited there for more details, and kindly cite those papers when using this code.
mod_Parameters.F90
Go to the documentation of this file.
1 MODULE modparameters
2 implicit none
3 save
4 !
5 !
6 character(len=*),parameter :: jhugen_version="v7.5.4"
7 !
8 !
9 !=====================================================
10 !internal
11 integer, public, parameter :: dp = selected_real_kind(15)
12 real(8), public, parameter :: tol = 0.0000001d0
13 integer, public, parameter :: invalidmode=-1,wwmode=00,zzmode=01,zgsmode=02,gszmode=03,gsgsmode=04,zgmode=05,gsgmode=06,ggmode=07
14 integer, public, parameter :: wwpmode=10,wpwmode=11,wpwpmode=12
15 integer, public, parameter :: zzpmode=20,zpzmode=21,zpzpmode=22
16 integer, public, parameter :: gszpmode=30,zpgsmode=31,zpgmode=32
17 integer, public :: collider,pchannel,process,decaymode1,decaymode2,topdecays,taudecays
18 integer, public :: vegasit1,vegasnc0,vegasnc1,vegasnc2,pmzzevals
19 real(8), public :: collider_energy
20 integer, public :: facscheme,renscheme
21 real(8), public :: mufacmultiplier,murenmultiplier,customupperscalebound, customlowerscalebound
22 integer, public :: vegasit1_default,vegasnc0_default,vegasnc1_default,vegasnc2_default
23 integer, public :: numhistograms
24 integer, public :: requestnleptons(1:2) = -1
25 integer, public :: requestnjets(1:2) = -1
26 integer, public :: requestos(1:2) = -1
27 integer, public :: requestossf(1:2) = -1
28 logical, public :: unweighted,offshellreson,offshellv1,offshellv2,readlhefile,convertlhefile,doprintpmzz
29 logical, public :: readcsmax,generateevents,counttauasany,hasleptonfilter, foundhiggsmass, foundhiggswidth
30 integer, public :: writefailedevents
31 logical, public :: filesopened = .false.
32 integer, public, parameter :: krenfacscheme_default=0
33 integer, public, parameter :: krenfacscheme_mhstar=1
34 integer, public, parameter :: krenfacscheme_mjjhstar=2
35 integer, public, parameter :: krenfacscheme_mjj_mhstar=3
36 integer, public, parameter :: krenfacscheme_mj_mj_mhstar=4
37 integer, public, parameter :: krenfacscheme_mjj=5
38 integer, public, parameter :: krenfacscheme_mj_mj=6
39 integer, public, parameter :: krenfacscheme_mjhstar=7
40 integer, public, parameter :: krenfacscheme_mj_mhstar=8
41 integer, public, parameter :: krenfacscheme_mj=9
42 integer, public, parameter :: krenfacscheme_maxptj=10
43 integer, public, parameter :: krenfacscheme_minptj=11
44 integer, public, parameter :: krenfacscheme_custom_scale=12
45 integer, public, parameter :: nrenfacschemes=13
46 integer, public, parameter :: maxpart = 30
47 integer, public, parameter :: nmaxchannels=200 ! Maximum 250 from vegas_common xi dimensions
48 integer(8), public :: evalcounter=0
49 integer(8), public :: rejecounter=0
50 integer(8), public :: accepcounter=0
51 integer(8), public :: alertcounter=0
52 integer(8), public :: accepcounter_part(-6:6,-6:6)=0,rejecounter_part(-6:6,-6:6)=0,requevents(-6:+6,-6:+6)=0
53 real(8), public :: crosssecmax(-6:+6,-6:+6)=0,crosssec(-6:+6,-6:+6)=0
54 integer(8), public :: requevents2(1:nmaxchannels)=0
55 integer(8), public :: accepcounter_part2(1:nmaxchannels)=0
56 real(8), public :: crosssec2(1:nmaxchannels)=0
57 real(8), public :: crosssecmax2(1:nmaxchannels)=0
58 integer, public :: ipart_sel, jpart_sel, ichann_sel
59 real(8) :: time_start,time_end,time_int
60 logical, public :: warmup
61 character(len=500) :: datafile
62 character(len=500) :: csmaxfile
63 character(len=500) :: logfile
64 character(len=500) :: lheprodfile
65 ! PDFset variables, present regardless of useLHAPDF value due to MELA
66 character(len=100) :: lhapdfstring
67 character(len=500) :: lhapdf_data_path
68 integer, public :: lhapdfmember, lenlhapdfstring ! lenLHAPDFString is needed in MELA
69 integer, public :: pdfset
70 ! End PDFset variables
71 #if useCollier==1
72 ! COLLIER initialization variables
73 integer, public :: collier_maxnloopprops = -1
74 integer, public :: collier_maxrank = -1
75 ! End COLLIER initialization variables
76 #endif
77 logical, public :: includeinterference, reweightinterference, writegit
78 real(8), public :: m_v,ga_v, m_vprime,ga_vprime, m_v_ps,ga_v_ps, m_z_ps,ga_z_ps, m_w_ps,ga_w_ps
79 real(dp), public, parameter :: pi =3.141592653589793238462643383279502884197_dp
80 real(dp), public, parameter :: sqrt2 = 1.4142135623730950488016887242096980786_dp
81 real(dp), public, parameter :: gamma_0 = 0.5772156649015328606065120900824024310421_dp !Euler–Mascheroni constant
82 real(dp), public, parameter :: pisq = pi**2
83 real(8), public, parameter :: one = 1.0d0, mone = -1.0d0
84 real(8), public, parameter :: half = 0.5d0,two = 2.0d0
85 real(8), public, parameter :: zero = 0d0
86 complex(8), parameter, public :: czero = (0d0,0d0)
87 complex(8), parameter, public :: cone = 1.0d0
88 complex(8), parameter, public :: ci=(0d0,1.0d0)
89 complex(8), parameter, public :: ne=(0d0,1.0d0)
90 real(8), public, parameter :: percent=1d0/100d0
91 real(8), public, parameter :: symmfac=1d0/2d0, spinavg=1d0/4d0, quarkcolavg=1d0/3d0, gluoncolavg=1d0/8d0
92 ! Units
93 real(8), public, parameter :: gev=1d0/100d0 ! we are using units of 100GeV, i.e. Lambda=10 is 1TeV
94 real(8), public, parameter :: tev=1d+3*gev
95 real(8), public, parameter :: mev=1d-3*gev
96 real(8), public, parameter :: kev=1d-6*gev
97 real(8), public, parameter :: ev=1d-9*gev
98 real(8), public, parameter :: meter=1d6
99 real(8), public, parameter :: milimeter=1d-3*meter
100 real(8), public, parameter :: ctauunit=milimeter
101 real(8), public, parameter :: barn=1d-28*(meter**2)
102 real(8), public, parameter :: picobarn=1e-12*barn
103 real(8), public, parameter :: femtobarn=1e-15*barn
104 real(8), public, parameter :: attobarn=1e-18*barn
105 real(8), public, parameter :: second=1d15
106 real(8), public, parameter :: milisecond=1d-3*second
107 real(8), public, parameter :: speedoflight=299792458*meter/second
108 real(8), public, parameter :: planckconstant=4.13567d-15*ev*second
109 real(8), public, parameter :: hbar=planckconstant/(2d0*pi)
110 real(8), public, parameter :: hbarc=hbar*speedoflight
111 real(8), public, parameter :: hbarc2=hbarc**2
112 real(8), public, parameter :: hbarc2_fbgev2 = hbarc2/(femtobarn*(gev**2))
113 real(8), public, parameter :: hbarc2_pbgev2 = hbarc2/(picobarn*(gev**2))
114 real(8), public, parameter :: hbarc2xsecunit = hbarc2_fbgev2*(gev**2) ! Units of xsec are reported in fb. Mutliplication by GeV**2 is to cancel out the units of ME*PS.
115 ! real(8),public :: GlobalMax=-1d99
116 ! real(8),public :: GlobalMin=+1d99
117 ! integer,parameter :: NPart=200
118 ! real(8),public :: PartitionMax(0:NPart,0:NPart)=-1d99
119 real(8),public :: mincs=1d10,maxcs=0d0,avgcs=0d0
120 integer, public :: br_z_ll_counter=0
121 integer, public :: br_z_inv_counter=0
122 integer, public :: br_z_uu_counter=0
123 integer, public :: br_z_dd_counter=0
124 integer, public :: br_w_ll_counter=0
125 integer, public :: br_w_ud_counter=0
126 integer, public :: br_counter(1:5,1:5)=0
127 integer, public :: leptinevent(0:8) = 0
128 integer, public :: jetsinevent(0:8) = 0
129 logical, public :: reweightdecay = .false.
130 integer, public :: userseed = 0
131 integer, public :: widthscheme = 0 ! 1=running BW-width, 2=fixed BW-width (default), 3=Passarino's CPS
132 integer, public :: widthschemein = 0 ! 1=running BW-width, 2=fixed BW-width (default), 3=Passarino's CPS
133 real(8), public :: mubarh = -999d0 !for CPS
134 real(8), public :: gabarh = -999d0 !for CPS
135 real(8), public :: maxinputmhstar = -999d0, mininputmhstar = 1d15, mhstarforphasespace
136 logical, public :: readpmzz
137 character(len=500), public :: pmzzfile = "PMZZdistribution.out"
138 real(8), public :: pmzz_mreso = -1d0
139 integer, public, parameter :: pmzzsize = 10000
140 real(8), public :: pmzzdistribution(1:pmzzsize,1:2) !huge array, in normal cases will never get near the edge
141 integer, public :: pmzzminindex=-1, pmzzmaxindex=-1 !store the largest and smallest values currently used
142 complex(8), public :: printpmzz !real part is the minimum, imaginary part is the maximum
143 integer, public :: printpmzzintervals
144 integer, public :: vbfoffsh_run=-1
145 logical, public :: findcrosssectionwithweights = .false.
146 real(8), public :: vprimedecaylengthmasscutofffactor = -1d0 ! x Ga_Vprime
147 real(8), public :: crosssectionwithweights = 0d0, crosssectionwithweightserrorsquared = 0d0
148 !=====================================================
149 
150 
151 
152 !=====================================================
153 !these fixed random seeds are NOT actually used to generate events
154 !the seed provided via the command line (or, if none is provided, one generated randomly using the system time)
155 ! is put in place of the first seed that makes a difference in generating a random number
156 ! then the seeds are used to generate the (compiler dependent) required number of seeds
157 ! and THOSE are used for event generation
158 !note that even if the same seed is provided, results are compiler dependent
159 integer, public, parameter :: nmaxseeds = 33
160 integer, public, parameter :: defaultseeds(1:nmaxseeds) = (/847362834,470115596,392845769,191039475,372910496,192049687,695820194,218930493,902943834,471748302,123958674,390534012,938576849,386472918,938576483,891928354,593857698,938576432,948576849,192847564 ,218940493,902943854,471748304,123958670,390534013,938576846,386472917,938576484,891928359,593857698,938576434,948576848,192847565/)
161 integer, public :: theseeds(1:nmaxseeds) = defaultseeds
162 !changing the default seeds is not advised, since then results from before the change
163 ! will not be reproducible
164 !=====================================================
165 
166 !=====================================================
167 !switches - should be set on the command line
168 logical, public :: fix_channels_ratio = .false.
169 
170 real(8), public :: channels_ratio_fix = 0.25d0 ! desired ratio of N_qq/(N_qq+N_gg)
171 
172 logical, public :: writeweightedlhe = .false.
173 
174 logical, public :: randomizevvdecays = .true. ! randomize DecayMode1 and DecayMode2 in H-->VV and TTBAR decays whenever appropriate
175 
176 logical, public :: useunformattedread = .false. !Set this to true if the regular reading fails for whatever reason
177 
178 logical, public :: hbbdecays =.false. ! default to false so H in V* > VH (Process = 51) does not decay to bbbar
179 
180 logical, public :: h_dk =.false. ! default to false so H in V* > VH (Process = 50) does not decay to bbbar
181 !=====================================================
182 
183 ! new VH
184 character(len=2), public :: vh_pc = "lo" ! VH partonic channel and mode selection, in development.
185  ! "ee" ( = e+ e- @LO)
186  ! "gg" ( = triangles + boxes of gg)
187  ! "qq" ( = q q~ @LO)
188  ! "lo" ( = q q~ @LO)
189  ! "tr" ( = triangles of gg)
190  ! "bo" ( = boxes of gg)
191  ! "in" ( = interference = 2*dble(box*dconjg(triangle)) of gg)
192  ! "qg" ( = real - dipoles, for g q/q~ > VH + q/q~, for development only)
193  ! "gq" ( = virtual + I + K + P, for g q/q~ > VH + q/q~, for development only)
194  ! "sb" ( = real - dipoles, for q q~ @NLO, for development only)
195  ! "sp" ( = virtual + I + K + P, for q q~ @NLO, for development only)
196  ! "nl" ( = full oneloop = q ~ @LO + NLO + gg + gq)
197  ! VH_PC overrides Pchannel.
198 
199 real(8), public :: alpha_dip = 1d0 !extra non physical degree of freedom for dipoles. Vary to check indepedence (of alpha_dip).
200 ! new VH
201 !=====================================================
202 
203 !=====================================================
204 !cuts - should be set on the command line
205 real(8), public :: ptjetcut = -1d0*gev ! jet min pt
206 real(8), public :: etajetcut = -1d0 ! jet max |eta|
207 real(8), public :: detajetcut = -1d0 ! min difference in eta between jets
208 real(8), public :: rjet = -1d0 ! jet deltaR
209 real(8), public :: mjjcut = 0d0*gev ! minimum mJJ
210 real(8), public :: m4l_minmax(1:2) = (/ -1d0,-1d0 /)*gev ! min and max for m_4l. default is (-1,-1): m_4l ~ Higgs resonance (on-shell)
211 real(8), public :: m2l_minmax(1:2) = (/ 0d0,14000d0 /)*gev ! min and max for m_V in VH production;
212 real(8), public :: mvh_minmax(1:2) = (/ 0d0,14000d0 /)*gev ! min and max for m_VH in VH production;
213 logical, public :: includegammastar = .false. ! include offshell photons?
214 logical, public :: includevprime = .false.
215 real(8), public :: mphotoncutoff = -1d0*gev ! minimum |mass_ll| for offshell photons when includeGammaStar = .true. or in VBF bkg
216 real(8), public :: ptlepcut = -1d0*gev ! lepton min pt
217 real(8), public :: etalepcut = 999d0 ! lepton max |eta|
218 real(8), public :: pthcut = 0d0*gev
219 logical, public :: jetsoppositeeta = .false. ! Ensures associated jets are in opposite hemispheres. Useful to obtain VBF/VBS topology from inclusive production
220 !=====================================================
221 
222 !=====================================================
223 !constants
224 real(8), public :: m_top = 173.2d0 *gev ! top quark mass
225 real(8), public :: ga_top = 2.0d0 *gev ! top quark width
226 real(8), public :: m_z = 91.1876d0 *gev ! Z boson mass (PDG-2011)
227 real(8), public :: ga_z = 2.4952d0 *gev ! Z boson width(PDG-2011)
228 real(8), public :: m_w = 80.399d0 *gev ! W boson mass (PDG-2011)
229 real(8), public :: ga_w = 2.085d0 *gev ! W boson width(PDG-2011)
230 real(8), public :: m_reso = 125.0d0 *gev ! X resonance mass (spin 0, spin 1, spin 2, can be overwritten by command line argument)
231 real(8), public :: ga_reso = 0.00407d0 *gev ! X resonance width
232 real(8), public :: m_reso2 = -1d0 *gev ! second resonance mass (spin 0 in off-shell MCFM, can be overwritten by command line argument)
233 real(8), public :: ga_reso2= 0d0 *gev ! second resonance width
234 
235 real(8), public :: higgsdecaylengthmm = 0d0 ! Higgs decay length in [mm]
236 real(8), public :: vprimedecaylengthmm = 0d0 ! Vprime decay length in [mm]
237 
238 real(8), public :: m_bot = 4.75d0 *gev ! bottom quark mass
239 real(8), public :: m_charm = 1.275d0 *gev ! charm quark mass
240 real(8), public :: m_el = 0.00051100d0 *gev ! electron mass
241 real(8), public :: m_mu = 0.10566d0 *gev ! muon mass
242 real(8), public :: m_tau = 1.7768d0 *gev ! tau mass
243 real(8), public :: ga_tau = 2.267d-12 *gev ! tau width
244 real(8), public :: m_bot_4gen = 100000d0*gev ! b mass for fourth-generation loop in gg MCFM
245 real(8), public :: m_top_4gen = 100000d0*gev ! t mass for fourth-generation loop in gg MCFM
246 
247 real(8), public :: gf = 1.16639d-5/gev**2 ! Fermi constant
248 real(8), public :: alpha_qed = 1d0/128d0 ! el.magn. coupling
249 real(8), public :: vev ! = 1.0d0/sqrt(Gf*sqrt(2.0d0))
250 real(8), public :: gwsq ! = 4.0d0 * M_W**2/vev**2 ! weak constant squared
251 real(dp), public :: esq ! = 4.0d0 * pi * alpha_QED ! Fundamental charge
252 
253 real(8), public :: xw = 0.23119d0 ! sin**2(Theta_Weinberg) (PDG-2008)
254 real(8), public :: sitw ! = dsqrt(xw) ! sin(Theta_Weinberg) (PDG-2008)
255 real(8), public :: twosc ! = sqrt(4.0_dp*xw*(1.0_dp-xw))
256 real(8), public, parameter :: lhc_energy=13d0 *tev ! LHC hadronic center of mass energy
257 real(8), public, parameter :: tev_energy=1.96d0 *tev ! Tevatron hadronic center of mass energy
258 real(8), public, parameter :: ilc_energy=250d0 *gev ! Linear collider center of mass energy
259 !command line: epPolarization, emPolarization
260 real(8), public :: pol_a = 0d0 ! e+ polarization. 0: no polarization, 100: helicity = 1, -100: helicity = -1
261 real(8), public :: pol_b = 0d0 ! e- polarization. 0: no polarization, 100: helicity = 1, -100: helicity = -1
262 
263 ! PDF and QCD scale variables, set in main::InitPDFNonConstVals if not a parameter
264 integer, public, parameter :: nqflavors_pdf = 5 ! Number of flavors enforced to the PDF, used in ModParameters::EvalAlphaS()
265 integer, public, parameter :: nloops_pdf = 1 ! alpha_s order
266 real(8), public :: zmass_pdf ! Z mass used in pdf toward the QCD scale, reset later in main per PDF if needed
267 real(8), public :: mu_fact ! pdf factorization scale (set to M_Reso in main.F90)
268 real(8), public :: mu_ren ! QCD renormalization (alpha_s) scale (set to M_Reso in main.F90)
269 real(dp), public :: alphas ! strong coupling per event, set to some reasonable value
270 real(dp), public :: alphas_mz ! strong coupling at M_Z, reset later in main per PDF
271 real(dp), public :: gs ! = sqrt(alphas*4.0_dp*pi)
272 real(8), public :: breitwignercutoff = 20d0 ! This is the cutoff to the phase space calculation in mod_CrossSection.F90
273 
274 
275 ! CKM squared matrix entries
276 real(8), public :: vckm_ud! = 0.974285d0
277 real(8), public :: vckm_us! = 0.225290d0
278 real(8), public :: vckm_cs! = 0.9734244d0
279 real(8), public :: vckm_cd! =-0.225182d0
280 real(8), public :: vckm_tb! = 0.99912367d0
281 real(8), public :: vckm_ts! =-0.040920069d0
282 real(8), public :: vckm_cb! = dsqrt(1d0-VCKM_cd**2-VCKM_cs**2)
283 real(8), public :: vckm_ub! = dsqrt(1d0-VCKM_ud**2-VCKM_us**2)
284 real(8), public :: vckm_td! = dsqrt(1d0-VCKM_tb**2-VCKM_ts**2)
285 
286 
287 ! absolute branching fraction (taken from PDG-2014)
288 real(8), public, parameter :: br_z_ll = 10.10d0*percent ! leptonic Z branching
289 real(8), public, parameter :: br_z_hadr = 69.91d0*percent ! hadronic Z branching
290 real(8), public, parameter :: br_z_inv = 100d0*percent - br_z_ll - br_z_hadr ! invisible Z branching
291 real(8), public, parameter :: br_z_uu = 11.6d0*percent ! up upbar Z branching
292 real(8), public, parameter :: br_z_cc = 11.6d0*percent ! chm chmbar Z branching
293 real(8), public, parameter :: br_z_dd = 15.6d0*percent ! dn dnbar Z branching
294 real(8), public, parameter :: br_z_ss = 15.6d0*percent ! str strbar Z branching
295 real(8), public, parameter :: br_z_bb = br_z_hadr - br_z_uu - br_z_dd - br_z_cc - br_z_ss
296 real(8), public, parameter :: br_w_ll = 32.72d0*percent ! leptonic W branching
297 real(8), public, parameter :: br_w_hadr = 100d0*percent - br_w_ll ! hadronic W branching
298 
299 
300 ! derived branching fractions
301 real(8), public, parameter :: br_z_ee = 1d0/3d0*br_z_ll ! electron Z branching
302 real(8), public, parameter :: br_z_mm = 1d0/3d0*br_z_ll ! muon Z branching
303 real(8), public, parameter :: br_z_tt = 1d0/3d0*br_z_ll ! tau Z branching
304 real(8), public, parameter :: br_z_nn = 1d0/3d0*br_z_inv ! neutrino Z branching
305 real(8), public, parameter :: br_w_en = 1d0/3d0*br_w_ll ! electron W branching
306 real(8), public, parameter :: br_w_mn = 1d0/3d0*br_w_ll ! muon W branching
307 real(8), public, parameter :: br_w_tn = 1d0/3d0*br_w_ll ! electron W branching
308 real(8), public, parameter :: br_w_ud = 1d0/2d0*br_w_hadr ! u-d W branching
309 real(8), public, parameter :: br_w_cs = 1d0/2d0*br_w_hadr ! c-s W branching
310 
311 real(8), public, parameter :: brlept_z_ee = br_z_ee/br_z_ll ! Z branching fraction Ga(el)/Ga(leptonic)
312 real(8), public, parameter :: brlept_z_mm = br_z_mm/br_z_ll ! Z branching fraction Ga(mu)/Ga(leptonic)
313 real(8), public, parameter :: brlept_z_tt = br_z_tt/br_z_ll ! Z branching fraction Ga(tau)/Ga(leptonic)
314 real(8), public, parameter :: brlept_z_nn = br_z_nn/br_z_inv ! Z branching fraction Ga(neu)/Ga(invisible)
315 real(8), public, parameter :: brlept_w_en = br_w_en/br_w_ll ! W branching fraction Ga(el)/Ga(leptonic)
316 real(8), public, parameter :: brlept_w_mn = br_w_mn/br_w_ll ! W branching fraction Ga(mu)/Ga(leptonic)
317 real(8), public, parameter :: brlept_w_tn = br_w_tn/br_w_ll ! W branching fraction Ga(tau)/Ga(leptonic)
318 
319 real(8), public, parameter :: brhadr_z_uu = br_z_uu/br_z_hadr ! Z branching fraction Ga(up)/Ga(hadronic)
320 real(8), public, parameter :: brhadr_z_cc = br_z_cc/br_z_hadr ! Z branching fraction Ga(chm)/Ga(hadronic)
321 real(8), public, parameter :: brhadr_z_dd = br_z_dd/br_z_hadr ! Z branching fraction Ga(don)/Ga(hadronic)
322 real(8), public, parameter :: brhadr_z_ss = br_z_ss/br_z_hadr ! Z branching fraction Ga(str)/Ga(hadronic)
323 real(8), public, parameter :: brhadr_z_bb = br_z_bb/br_z_hadr ! Z branching fraction Ga(bot)/Ga(hadronic)
324 real(8), public, parameter :: brhadr_w_ud = br_w_ud/br_w_hadr ! W branching fraction Ga(up)/Ga(hadronic)
325 real(8), public, parameter :: brhadr_w_cs = br_w_cs/br_w_hadr ! W branching fraction Ga(chm)/Ga(hadronic)
326 
327 real(8), public :: scale_alpha_z_uu = 1.037560d0 ! scaling factor of alpha (~partial width) for Z > u u~, c c~
328 real(8), public :: scale_alpha_z_dd = 1.037560d0 ! scaling factor of alpha (~partial width) for Z > d d~, s s~, b b~
329 real(8), public :: scale_alpha_z_ll = 1d0 ! scaling factor of alpha (~partial width) for Z > l+ l- (l=e,mu)
330 real(8), public :: scale_alpha_z_tt = 1d0 ! scaling factor of alpha (~partial width) for Z > tau+ tau-
331 real(8), public :: scale_alpha_z_nn = 1d0 ! scaling factor of alpha (~partial width) for Z > nu nu~
332 real(8), public :: scale_alpha_w_ud = 1.038200d0 ! scaling factor of alpha (~partial width) for W > u d
333 real(8), public :: scale_alpha_w_cs = 1.038200d0 ! scaling factor of alpha (~partial width) for W > c s
334 real(8), public :: scale_alpha_w_ln = 1d0 ! scaling factor of alpha (~partial width) for W > l nu (l=e,mu)
335 real(8), public :: scale_alpha_w_tn = 1d0 ! scaling factor of alpha (~partial width) for W > tau nu
336 ! real(8), public, parameter :: scale_alpha_Z_ll = (1d0-(Br_Z_uu+Br_Z_cc)*scale_alpha_Z_uu-(Br_Z_dd+Br_Z_ss+Br_Z_bb)*scale_alpha_Z_dd)/(3d0*Br_Z_nn+3d0*Br_Z_ee) ! scaling factor of alpha (~partial width) for Z > l+ l- which restores the total width
337 ! real(8), public, parameter :: scale_alpha_Z_nn = scale_alpha_Z_ll ! scaling factor of alpha (~partial width) for Z > nu nu~ which restores total width
338 ! real(8), public, parameter :: scale_alpha_W_ln = (1d0-Br_W_ud*scale_alpha_W_ud-Br_W_cs*scale_alpha_W_cs)/(3d0*Br_W_en) ! scaling factor of alpha (~partial width) for W > l nu
339 !
340 ! sum rules
341 ! 1 = 3*Br_Z_nn + 3*Br_Z_ee + 2*Br_Z_uu + 3*Br_Z_dd
342 !
343 ! sum rule with scaling factors
344 ! 1 = 3*(Br_Z_nn*scale_alpha_Z_nn) + 3*(Br_Z_ee*scale_alpha_Z_ll) + 2*(Br_Z_uu*scale_alpha_Z_uu) + 3*(Br_Z_dd*scale_alpha_Z_dd)
345 !=====================================================
346 
347 
348 !=====================================================
349 !resonance couplings
350 
351 ! Lambda scale enters in two places
352 ! overall scale for x-section and in power suppressed
353 ! operators/formfactors (former r).
354 real(8), public, parameter :: lambda = 1000d0 *gev
355 real(8), public, parameter :: lambda2 = 1000d0 *gev ! for second resonance
356 
357 !--------------------!
358 !-----! Spin-0 !-----!
359 !--------------------!
360 !-- parameters that define on-shell spin 0 coupling to SM fields, see note
361  logical, public, parameter :: generate_as = .false. ! .true. uses ah* instead of gh*
362  complex(8), public, parameter :: ahg1 = (1.0d0,0d0)
363  complex(8), public, parameter :: ahg2 = (0d0,0d0)
364  complex(8), public, parameter :: ahg3 = (0d0,0d0) ! pseudoscalar
365  complex(8), public, parameter :: ahz1 = (1.0d0,0d0)
366  complex(8), public, parameter :: ahz2 = (0d0,0d0) ! this coupling does not contribute for gamma+gamma final states
367  complex(8), public, parameter :: ahz3 = (0d0,0d0) ! pseudoscalar
368 
369 !-- parameters that define off-shell spin 0 coupling to SM fields, see note
370 !-- Hgg couplings to gluons for point-like vertices
371  complex(8), public :: ghg2 = (1.0d0,0d0)
372  complex(8), public :: ghg3 = (0d0,0d0)
373  complex(8), public :: ghg4 = (0d0,0d0) ! pseudoscalar
374 
375 !-- Hgg couplings to gluons for point-like vertices (4th generation vertex for MCFM)
376  complex(8), public :: ghg2_4gen = (0d0,0d0)
377  complex(8), public :: ghg3_4gen = (0d0,0d0)
378  complex(8), public :: ghg4_4gen = (0d0,0d0)
379 
380 !-- HVV' couplings to ZZ/ZA/AA and WW
381  complex(8), public :: ghz1 = (2.0d0,0d0) ! SM=2 (MCFM => =1)
382  complex(8), public :: ghz2 = (0d0,0d0)
383  complex(8), public :: ghz3 = (0d0,0d0)
384  complex(8), public :: ghz4 = (0d0,0d0) ! pseudoscalar
385 
386 !-- parameters that define q^2 dependent form factors
387  complex(8), public :: ghz1_prime = (0d0,0d0)
388  complex(8), public :: ghz1_prime2= (0d0,0d0)
389  complex(8), public :: ghz1_prime3= (0d0,0d0)
390  complex(8), public :: ghz1_prime4= (0d0,0d0)
391  complex(8), public :: ghz1_prime5= (0d0,0d0)
392  complex(8), public :: ghz1_prime6= (0d0,0d0)
393  complex(8), public :: ghz1_prime7= (0d0,0d0)
394 
395  complex(8), public :: ghz2_prime = (0d0,0d0)
396  complex(8), public :: ghz2_prime2= (0d0,0d0)
397  complex(8), public :: ghz2_prime3= (0d0,0d0)
398  complex(8), public :: ghz2_prime4= (0d0,0d0)
399  complex(8), public :: ghz2_prime5= (0d0,0d0)
400  complex(8), public :: ghz2_prime6= (0d0,0d0)
401  complex(8), public :: ghz2_prime7= (0d0,0d0)
402 
403  complex(8), public :: ghz3_prime = (0d0,0d0)
404  complex(8), public :: ghz3_prime2= (0d0,0d0)
405  complex(8), public :: ghz3_prime3= (0d0,0d0)
406  complex(8), public :: ghz3_prime4= (0d0,0d0)
407  complex(8), public :: ghz3_prime5= (0d0,0d0)
408  complex(8), public :: ghz3_prime6= (0d0,0d0)
409  complex(8), public :: ghz3_prime7= (0d0,0d0)
410 
411  complex(8), public :: ghz4_prime = (0d0,0d0)
412  complex(8), public :: ghz4_prime2= (0d0,0d0)
413  complex(8), public :: ghz4_prime3= (0d0,0d0)
414  complex(8), public :: ghz4_prime4= (0d0,0d0)
415  complex(8), public :: ghz4_prime5= (0d0,0d0)
416  complex(8), public :: ghz4_prime6= (0d0,0d0)
417  complex(8), public :: ghz4_prime7= (0d0,0d0)
418 
419  complex(8), public :: ghzgs1_prime2= (0d0,0d0)
420  complex(8), public :: ghzgs2 = (0d0,0d0)
421  complex(8), public :: ghzgs3 = (0d0,0d0)
422  complex(8), public :: ghzgs4 = (0d0,0d0)
423  complex(8), public :: ghgsgs2 = (0d0,0d0)
424  complex(8), public :: ghgsgs3 = (0d0,0d0)
425  complex(8), public :: ghgsgs4 = (0d0,0d0)
426 
427  real(8), public, parameter :: lambda_z1 = 10000d0*gev
428  real(8), public, parameter :: lambda_z2 = 10000d0*gev
429  real(8), public, parameter :: lambda_z3 = 10000d0*gev
430  real(8), public, parameter :: lambda_z4 = 10000d0*gev
431  real(8), public, parameter :: lambda_zgs1 = 10000d0*gev
432  real(8), public, parameter :: lambda_q = 10000d0*gev
433 
434  integer, public :: cz_q1sq = 0 ! Sign of q1,2,12**2 for the following Lambda's, set to 1 or -1 to get q**2-dependence from these form factor Lambdas
435  integer, public :: cz_q2sq = 0
436  integer, public :: cz_q12sq = 0
437  ! These Lambdas all have a numerical value of 1d0
438  real(8), public :: lambda_z11 = 100d0*gev ! For Z1
439  real(8), public :: lambda_z21 = 100d0*gev
440  real(8), public :: lambda_z31 = 100d0*gev
441  real(8), public :: lambda_z41 = 100d0*gev
442  real(8), public :: lambda_z12 = 100d0*gev ! For Z2
443  real(8), public :: lambda_z22 = 100d0*gev
444  real(8), public :: lambda_z32 = 100d0*gev
445  real(8), public :: lambda_z42 = 100d0*gev
446  real(8), public :: lambda_z10 = 100d0*gev ! For the Higgs
447  real(8), public :: lambda_z20 = 100d0*gev
448  real(8), public :: lambda_z30 = 100d0*gev
449  real(8), public :: lambda_z40 = 100d0*gev
450 
451 !-- extra HWW couplings for weak boson fusion when WW-spin-0 couplings are required to be different from ZZ-spin-0
452 !-- note: ZZ-spin-0 couplings are used in processes other than VBF, and WW is distinguished from ZZ only in case distinguish_HWWcouplings=.true.
453  logical, public :: distinguish_hwwcouplings=.false.
454  complex(8), public :: ghw1 = (0d0,0d0)
455  complex(8), public :: ghw2 = (0d0,0d0)
456  complex(8), public :: ghw3 = (0d0,0d0)
457  complex(8), public :: ghw4 = (0d0,0d0)
458 
459 !-- parameters that define q^2 dependent form factors in WBF WW-spin-0 case described above
460  complex(8), public :: ghw1_prime = (0d0,0d0)
461  complex(8), public :: ghw1_prime2= (0d0,0d0)
462  complex(8), public :: ghw1_prime3= (0d0,0d0)
463  complex(8), public :: ghw1_prime4= (0d0,0d0)
464  complex(8), public :: ghw1_prime5= (0d0,0d0)
465  complex(8), public :: ghw1_prime6= (0d0,0d0)
466  complex(8), public :: ghw1_prime7= (0d0,0d0)
467 
468  complex(8), public :: ghw2_prime = (0d0,0d0)
469  complex(8), public :: ghw2_prime2= (0d0,0d0)
470  complex(8), public :: ghw2_prime3= (0d0,0d0)
471  complex(8), public :: ghw2_prime4= (0d0,0d0)
472  complex(8), public :: ghw2_prime5= (0d0,0d0)
473  complex(8), public :: ghw2_prime6= (0d0,0d0)
474  complex(8), public :: ghw2_prime7= (0d0,0d0)
475 
476  complex(8), public :: ghw3_prime = (0d0,0d0)
477  complex(8), public :: ghw3_prime2= (0d0,0d0)
478  complex(8), public :: ghw3_prime3= (0d0,0d0)
479  complex(8), public :: ghw3_prime4= (0d0,0d0)
480  complex(8), public :: ghw3_prime5= (0d0,0d0)
481  complex(8), public :: ghw3_prime6= (0d0,0d0)
482  complex(8), public :: ghw3_prime7= (0d0,0d0)
483 
484  complex(8), public :: ghw4_prime = (0d0,0d0)
485  complex(8), public :: ghw4_prime2= (0d0,0d0)
486  complex(8), public :: ghw4_prime3= (0d0,0d0)
487  complex(8), public :: ghw4_prime4= (0d0,0d0)
488  complex(8), public :: ghw4_prime5= (0d0,0d0)
489  complex(8), public :: ghw4_prime6= (0d0,0d0)
490  complex(8), public :: ghw4_prime7= (0d0,0d0)
491 
492  real(8), public, parameter :: lambda_w1 = 10000d0*gev
493  real(8), public, parameter :: lambda_w2 = 10000d0*gev
494  real(8), public, parameter :: lambda_w3 = 10000d0*gev
495  real(8), public, parameter :: lambda_w4 = 10000d0*gev
496 
497  integer, public :: cw_q1sq = 0 ! Sign of q1,2,12**2 for the following Lambda's, set to 1 or -1 to get q**2-dependence from these form factor Lambdas
498  integer, public :: cw_q2sq = 0
499  integer, public :: cw_q12sq = 0
500  real(8), public :: lambda_w11 = 100d0*gev ! For W+
501  real(8), public :: lambda_w21 = 100d0*gev
502  real(8), public :: lambda_w31 = 100d0*gev
503  real(8), public :: lambda_w41 = 100d0*gev
504  real(8), public :: lambda_w12 = 100d0*gev ! For W-
505  real(8), public :: lambda_w22 = 100d0*gev
506  real(8), public :: lambda_w32 = 100d0*gev
507  real(8), public :: lambda_w42 = 100d0*gev
508  real(8), public :: lambda_w10 = 100d0*gev ! For the Higgs
509  real(8), public :: lambda_w20 = 100d0*gev
510  real(8), public :: lambda_w30 = 100d0*gev
511  real(8), public :: lambda_w40 = 100d0*gev
512 
513 !-- HVV contact terms
514  complex(8), public :: ghzzp1 = (0d0,0d0)
515  complex(8), public :: ghzzp2 = (0d0,0d0)
516  complex(8), public :: ghzzp3 = (0d0,0d0)
517  complex(8), public :: ghzzp4 = (0d0,0d0)
518 
519 !-- parameters that define q^2 dependent form factors
520  complex(8), public :: ghzzp1_prime = (0d0,0d0)
521  complex(8), public :: ghzzp1_prime2= (0d0,0d0)
522  complex(8), public :: ghzzp1_prime3= (0d0,0d0)
523  complex(8), public :: ghzzp1_prime4= (0d0,0d0)
524  complex(8), public :: ghzzp1_prime5= (0d0,0d0)
525  complex(8), public :: ghzzp1_prime6= (0d0,0d0)
526  complex(8), public :: ghzzp1_prime7= (0d0,0d0)
527 
528  complex(8), public :: ghzzp2_prime = (0d0,0d0)
529  complex(8), public :: ghzzp2_prime2= (0d0,0d0)
530  complex(8), public :: ghzzp2_prime3= (0d0,0d0)
531  complex(8), public :: ghzzp2_prime4= (0d0,0d0)
532  complex(8), public :: ghzzp2_prime5= (0d0,0d0)
533  complex(8), public :: ghzzp2_prime6= (0d0,0d0)
534  complex(8), public :: ghzzp2_prime7= (0d0,0d0)
535 
536  complex(8), public :: ghzzp3_prime = (0d0,0d0)
537  complex(8), public :: ghzzp3_prime2= (0d0,0d0)
538  complex(8), public :: ghzzp3_prime3= (0d0,0d0)
539  complex(8), public :: ghzzp3_prime4= (0d0,0d0)
540  complex(8), public :: ghzzp3_prime5= (0d0,0d0)
541  complex(8), public :: ghzzp3_prime6= (0d0,0d0)
542  complex(8), public :: ghzzp3_prime7= (0d0,0d0)
543 
544  complex(8), public :: ghzzp4_prime = (0d0,0d0)
545  complex(8), public :: ghzzp4_prime2= (0d0,0d0)
546  complex(8), public :: ghzzp4_prime3= (0d0,0d0)
547  complex(8), public :: ghzzp4_prime4= (0d0,0d0)
548  complex(8), public :: ghzzp4_prime5= (0d0,0d0)
549  complex(8), public :: ghzzp4_prime6= (0d0,0d0)
550  complex(8), public :: ghzzp4_prime7= (0d0,0d0)
551 
552 !-- Zpgs
553  complex(8), public :: ghzpgs1_prime2= (0d0,0d0)
554  complex(8), public :: ghzpgs2 = (0d0,0d0)
555  complex(8), public :: ghzpgs3 = (0d0,0d0)
556  complex(8), public :: ghzpgs4 = (0d0,0d0)
557 
558 !-- ZpZp
559  complex(8), public :: ghzpzp1 = (0d0,0d0)
560  complex(8), public :: ghzpzp2 = (0d0,0d0)
561  complex(8), public :: ghzpzp3 = (0d0,0d0)
562  complex(8), public :: ghzpzp4 = (0d0,0d0)
563 
564 !-- parameters that define q^2 dependent form factors
565  complex(8), public :: ghzpzp1_prime = (0d0,0d0)
566  complex(8), public :: ghzpzp1_prime2= (0d0,0d0)
567  complex(8), public :: ghzpzp1_prime3= (0d0,0d0)
568  complex(8), public :: ghzpzp1_prime4= (0d0,0d0)
569  complex(8), public :: ghzpzp1_prime5= (0d0,0d0)
570  complex(8), public :: ghzpzp1_prime6= (0d0,0d0)
571  complex(8), public :: ghzpzp1_prime7= (0d0,0d0)
572 
573  complex(8), public :: ghzpzp2_prime = (0d0,0d0)
574  complex(8), public :: ghzpzp2_prime2= (0d0,0d0)
575  complex(8), public :: ghzpzp2_prime3= (0d0,0d0)
576  complex(8), public :: ghzpzp2_prime4= (0d0,0d0)
577  complex(8), public :: ghzpzp2_prime5= (0d0,0d0)
578  complex(8), public :: ghzpzp2_prime6= (0d0,0d0)
579  complex(8), public :: ghzpzp2_prime7= (0d0,0d0)
580 
581  complex(8), public :: ghzpzp3_prime = (0d0,0d0)
582  complex(8), public :: ghzpzp3_prime2= (0d0,0d0)
583  complex(8), public :: ghzpzp3_prime3= (0d0,0d0)
584  complex(8), public :: ghzpzp3_prime4= (0d0,0d0)
585  complex(8), public :: ghzpzp3_prime5= (0d0,0d0)
586  complex(8), public :: ghzpzp3_prime6= (0d0,0d0)
587  complex(8), public :: ghzpzp3_prime7= (0d0,0d0)
588 
589  complex(8), public :: ghzpzp4_prime = (0d0,0d0)
590  complex(8), public :: ghzpzp4_prime2= (0d0,0d0)
591  complex(8), public :: ghzpzp4_prime3= (0d0,0d0)
592  complex(8), public :: ghzpzp4_prime4= (0d0,0d0)
593  complex(8), public :: ghzpzp4_prime5= (0d0,0d0)
594  complex(8), public :: ghzpzp4_prime6= (0d0,0d0)
595  complex(8), public :: ghzpzp4_prime7= (0d0,0d0)
596 
597 
598  complex(8), public :: ezp_el_left = (0d0,0d0)
599  complex(8), public :: ezp_el_right = (0d0,0d0)
600  complex(8), public :: ezp_mu_left = (0d0,0d0)
601  complex(8), public :: ezp_mu_right = (0d0,0d0)
602  complex(8), public :: ezp_ta_left = (0d0,0d0)
603  complex(8), public :: ezp_ta_right = (0d0,0d0)
604  complex(8), public :: ezp_nue_left = (0d0,0d0) !same for NuMu and NuTau
605  complex(8), public :: ezp_nue_right = (0d0,0d0) !same for NuMu and NuTau
606  complex(8), public :: ezp_up_left = (0d0,0d0)
607  complex(8), public :: ezp_up_right = (0d0,0d0)
608  complex(8), public :: ezp_chm_left = (0d0,0d0)
609  complex(8), public :: ezp_chm_right = (0d0,0d0)
610  complex(8), public :: ezp_dn_left = (0d0,0d0)
611  complex(8), public :: ezp_dn_right = (0d0,0d0)
612  complex(8), public :: ezp_str_left = (0d0,0d0)
613  complex(8), public :: ezp_str_right = (0d0,0d0)
614  complex(8), public :: ezp_bot_left = (0d0,0d0)
615  complex(8), public :: ezp_bot_right = (0d0,0d0)
616  complex(8), public :: ezp_top_left = (0d0,0d0)
617  complex(8), public :: ezp_top_right = (0d0,0d0)
618 
619  real(8), public :: m_zprime = -1d0 ! <0: CT interaction, >=0: Heavy Zprime propagator
620  real(8), public :: ga_zprime = 0d0
621 !--
622 !-- HVV contact terms
623  complex(8), public :: ghwwp1 = (0d0,0d0)
624  complex(8), public :: ghwwp2 = (0d0,0d0)
625  complex(8), public :: ghwwp3 = (0d0,0d0)
626  complex(8), public :: ghwwp4 = (0d0,0d0)
627 
628 !-- parameters that define q^2 dependent form factors
629  complex(8), public :: ghwwp1_prime = (0d0,0d0)
630  complex(8), public :: ghwwp1_prime2= (0d0,0d0)
631  complex(8), public :: ghwwp1_prime3= (0d0,0d0)
632  complex(8), public :: ghwwp1_prime4= (0d0,0d0)
633  complex(8), public :: ghwwp1_prime5= (0d0,0d0)
634  complex(8), public :: ghwwp1_prime6= (0d0,0d0)
635  complex(8), public :: ghwwp1_prime7= (0d0,0d0)
636 
637  complex(8), public :: ghwwp2_prime = (0d0,0d0)
638  complex(8), public :: ghwwp2_prime2= (0d0,0d0)
639  complex(8), public :: ghwwp2_prime3= (0d0,0d0)
640  complex(8), public :: ghwwp2_prime4= (0d0,0d0)
641  complex(8), public :: ghwwp2_prime5= (0d0,0d0)
642  complex(8), public :: ghwwp2_prime6= (0d0,0d0)
643  complex(8), public :: ghwwp2_prime7= (0d0,0d0)
644 
645  complex(8), public :: ghwwp3_prime = (0d0,0d0)
646  complex(8), public :: ghwwp3_prime2= (0d0,0d0)
647  complex(8), public :: ghwwp3_prime3= (0d0,0d0)
648  complex(8), public :: ghwwp3_prime4= (0d0,0d0)
649  complex(8), public :: ghwwp3_prime5= (0d0,0d0)
650  complex(8), public :: ghwwp3_prime6= (0d0,0d0)
651  complex(8), public :: ghwwp3_prime7= (0d0,0d0)
652 
653  complex(8), public :: ghwwp4_prime = (0d0,0d0)
654  complex(8), public :: ghwwp4_prime2= (0d0,0d0)
655  complex(8), public :: ghwwp4_prime3= (0d0,0d0)
656  complex(8), public :: ghwwp4_prime4= (0d0,0d0)
657  complex(8), public :: ghwwp4_prime5= (0d0,0d0)
658  complex(8), public :: ghwwp4_prime6= (0d0,0d0)
659  complex(8), public :: ghwwp4_prime7= (0d0,0d0)
660 
661  complex(8), public :: ghwpwp1 = (0d0,0d0)
662  complex(8), public :: ghwpwp2 = (0d0,0d0)
663  complex(8), public :: ghwpwp3 = (0d0,0d0)
664  complex(8), public :: ghwpwp4 = (0d0,0d0)
665 
666 !-- parameters that define q^2 dependent form factors
667  complex(8), public :: ghwpwp1_prime = (0d0,0d0)
668  complex(8), public :: ghwpwp1_prime2= (0d0,0d0)
669  complex(8), public :: ghwpwp1_prime3= (0d0,0d0)
670  complex(8), public :: ghwpwp1_prime4= (0d0,0d0)
671  complex(8), public :: ghwpwp1_prime5= (0d0,0d0)
672  complex(8), public :: ghwpwp1_prime6= (0d0,0d0)
673  complex(8), public :: ghwpwp1_prime7= (0d0,0d0)
674 
675  complex(8), public :: ghwpwp2_prime = (0d0,0d0)
676  complex(8), public :: ghwpwp2_prime2= (0d0,0d0)
677  complex(8), public :: ghwpwp2_prime3= (0d0,0d0)
678  complex(8), public :: ghwpwp2_prime4= (0d0,0d0)
679  complex(8), public :: ghwpwp2_prime5= (0d0,0d0)
680  complex(8), public :: ghwpwp2_prime6= (0d0,0d0)
681  complex(8), public :: ghwpwp2_prime7= (0d0,0d0)
682 
683  complex(8), public :: ghwpwp3_prime = (0d0,0d0)
684  complex(8), public :: ghwpwp3_prime2= (0d0,0d0)
685  complex(8), public :: ghwpwp3_prime3= (0d0,0d0)
686  complex(8), public :: ghwpwp3_prime4= (0d0,0d0)
687  complex(8), public :: ghwpwp3_prime5= (0d0,0d0)
688  complex(8), public :: ghwpwp3_prime6= (0d0,0d0)
689  complex(8), public :: ghwpwp3_prime7= (0d0,0d0)
690 
691  complex(8), public :: ghwpwp4_prime = (0d0,0d0)
692  complex(8), public :: ghwpwp4_prime2= (0d0,0d0)
693  complex(8), public :: ghwpwp4_prime3= (0d0,0d0)
694  complex(8), public :: ghwpwp4_prime4= (0d0,0d0)
695  complex(8), public :: ghwpwp4_prime5= (0d0,0d0)
696  complex(8), public :: ghwpwp4_prime6= (0d0,0d0)
697  complex(8), public :: ghwpwp4_prime7= (0d0,0d0)
698 
699  complex(8), public :: ewp_el_left = (0d0,0d0)
700  complex(8), public :: ewp_el_right = (0d0,0d0)
701  complex(8), public :: ewp_mu_left = (0d0,0d0)
702  complex(8), public :: ewp_mu_right = (0d0,0d0)
703  complex(8), public :: ewp_ta_left = (0d0,0d0)
704  complex(8), public :: ewp_ta_right = (0d0,0d0)
705  complex(8), public :: ewp_up_left = (0d0,0d0)
706  complex(8), public :: ewp_up_right = (0d0,0d0)
707  complex(8), public :: ewp_chm_left = (0d0,0d0)
708  complex(8), public :: ewp_chm_right = (0d0,0d0)
709  complex(8), public :: ewp_top_left = (0d0,0d0)
710  complex(8), public :: ewp_top_right = (0d0,0d0)
711 
712  real(8), public :: m_wprime = -1d0 ! <0: CT interaction, >=0: Heavy Zprime propagator
713  real(8), public :: ga_wprime = 0d0
714 !--
715 
716 
717 
718 
719 !-- second resonance (H2) couplings for MCFM interface
720 !-- Hgg couplings to gluons for point-like vertices
721  complex(8), public :: gh2g2 = (0d0,0d0)
722  complex(8), public :: gh2g3 = (0d0,0d0)
723  complex(8), public :: gh2g4 = (0d0,0d0) ! pseudoscalar
724 
725 !-- Hgg couplings to gluons for point-like vertices (4th generation vertex for MCFM)
726  complex(8), public :: gh2g2_4gen = (0d0,0d0)
727  complex(8), public :: gh2g3_4gen = (0d0,0d0)
728  complex(8), public :: gh2g4_4gen = (0d0,0d0)
729 
730 !-- HVV' couplings to ZZ/ZA/AA and WW
731  complex(8), public :: gh2z1 = (0.0d0,0d0)
732  complex(8), public :: gh2z2 = (0.0d0,0d0)
733  complex(8), public :: gh2z3 = (0.0d0,0d0)
734  complex(8), public :: gh2z4 = (0.0d0,0d0) ! pseudoscalar
735 
736  complex(8), public :: gh2zgs2 = (0.00d0,0d0)
737  complex(8), public :: gh2zgs3 = (0.00d0,0d0)
738  complex(8), public :: gh2zgs4 = (0.00d0,0d0)
739  complex(8), public :: gh2gsgs2 = (0.00d0,0d0)
740  complex(8), public :: gh2gsgs3 = (0.00d0,0d0)
741  complex(8), public :: gh2gsgs4 = (0.00d0,0d0)
742 
743 !-- parameters that define q^2 dependent form factors
744  complex(8), public :: gh2z1_prime = (0.0d0,0d0)
745  complex(8), public :: gh2z1_prime2= (0.0d0,0d0)
746  complex(8), public :: gh2z1_prime3= (0.0d0,0d0)
747  complex(8), public :: gh2z1_prime4= (0.0d0,0d0)
748  complex(8), public :: gh2z1_prime5= (0.0d0,0d0)
749  complex(8), public :: gh2z1_prime6= (0.0d0,0d0)
750  complex(8), public :: gh2z1_prime7= (0.0d0,0d0)
751 
752  complex(8), public :: gh2z2_prime = (0.0d0,0d0)
753  complex(8), public :: gh2z2_prime2= (0.0d0,0d0)
754  complex(8), public :: gh2z2_prime3= (0.0d0,0d0)
755  complex(8), public :: gh2z2_prime4= (0.0d0,0d0)
756  complex(8), public :: gh2z2_prime5= (0.0d0,0d0)
757  complex(8), public :: gh2z2_prime6= (0.0d0,0d0)
758  complex(8), public :: gh2z2_prime7= (0.0d0,0d0)
759 
760  complex(8), public :: gh2z3_prime = (0.0d0,0d0)
761  complex(8), public :: gh2z3_prime2= (0.0d0,0d0)
762  complex(8), public :: gh2z3_prime3= (0.0d0,0d0)
763  complex(8), public :: gh2z3_prime4= (0.0d0,0d0)
764  complex(8), public :: gh2z3_prime5= (0.0d0,0d0)
765  complex(8), public :: gh2z3_prime6= (0.0d0,0d0)
766  complex(8), public :: gh2z3_prime7= (0.0d0,0d0)
767 
768  complex(8), public :: gh2z4_prime = (0.0d0,0d0)
769  complex(8), public :: gh2z4_prime2= (0.0d0,0d0)
770  complex(8), public :: gh2z4_prime3= (0.0d0,0d0)
771  complex(8), public :: gh2z4_prime4= (0.0d0,0d0)
772  complex(8), public :: gh2z4_prime5= (0.0d0,0d0)
773  complex(8), public :: gh2z4_prime6= (0.0d0,0d0)
774  complex(8), public :: gh2z4_prime7= (0.0d0,0d0)
775 
776  complex(8), public :: gh2zgs1_prime2= (0.0d0,0d0)
777 
778  real(8), public, parameter :: lambda2_z1 = 10000d0*gev
779  real(8), public, parameter :: lambda2_z2 = 10000d0*gev
780  real(8), public, parameter :: lambda2_z3 = 10000d0*gev
781  real(8), public, parameter :: lambda2_z4 = 10000d0*gev
782  real(8), public, parameter :: lambda2_zgs1 = 10000d0*gev
783  real(8), public, parameter :: lambda2_q = 10000d0*gev
784 
785  integer, public :: c2z_q1sq = 0 ! Sign of q1,2,12**2 for the following Lambda2's, set to 1 or -1 to get q**2-dependence from these form factor Lambda2s
786  integer, public :: c2z_q2sq = 0
787  integer, public :: c2z_q12sq = 0
788  real(8), public :: lambda2_z11 = 100d0*gev ! For Z1
789  real(8), public :: lambda2_z21 = 100d0*gev
790  real(8), public :: lambda2_z31 = 100d0*gev
791  real(8), public :: lambda2_z41 = 100d0*gev
792  real(8), public :: lambda2_z12 = 100d0*gev ! For Z2
793  real(8), public :: lambda2_z22 = 100d0*gev
794  real(8), public :: lambda2_z32 = 100d0*gev
795  real(8), public :: lambda2_z42 = 100d0*gev
796  real(8), public :: lambda2_z10 = 100d0*gev ! For the Higgs
797  real(8), public :: lambda2_z20 = 100d0*gev
798  real(8), public :: lambda2_z30 = 100d0*gev
799  real(8), public :: lambda2_z40 = 100d0*gev
800 
801 
802 !-- extra HWW couplings for weak boson fusion when WW-spin-0 couplings are required to be different from ZZ-spin-0
803 !-- note: ZZ-spin-0 couplings are used in processes other than VBF, and WW is distinguished from ZZ only in case distinguish_HWWcouplings=.true.
804 ! logical, public :: distinguish_HWWcouplings=.false.
805  complex(8), public :: gh2w1 = (0.0d0,0d0)
806  complex(8), public :: gh2w2 = (0.0d0,0d0)
807  complex(8), public :: gh2w3 = (0.0d0,0d0)
808  complex(8), public :: gh2w4 = (0.0d0,0d0)
809 
810 !-- parameters that define q^2 dependent form factors in WBF WW-spin-0 case described above
811  complex(8), public :: gh2w1_prime = (0.0d0,0d0)
812  complex(8), public :: gh2w1_prime2= (0.0d0,0d0)
813  complex(8), public :: gh2w1_prime3= (0.0d0,0d0)
814  complex(8), public :: gh2w1_prime4= (0.0d0,0d0)
815  complex(8), public :: gh2w1_prime5= (0.0d0,0d0)
816  complex(8), public :: gh2w1_prime6= (0.0d0,0d0)
817  complex(8), public :: gh2w1_prime7= (0.0d0,0d0)
818 
819  complex(8), public :: gh2w2_prime = (0.0d0,0d0)
820  complex(8), public :: gh2w2_prime2= (0.0d0,0d0)
821  complex(8), public :: gh2w2_prime3= (0.0d0,0d0)
822  complex(8), public :: gh2w2_prime4= (0.0d0,0d0)
823  complex(8), public :: gh2w2_prime5= (0.0d0,0d0)
824  complex(8), public :: gh2w2_prime6= (0.0d0,0d0)
825  complex(8), public :: gh2w2_prime7= (0.0d0,0d0)
826 
827  complex(8), public :: gh2w3_prime = (0.0d0,0d0)
828  complex(8), public :: gh2w3_prime2= (0.0d0,0d0)
829  complex(8), public :: gh2w3_prime3= (0.0d0,0d0)
830  complex(8), public :: gh2w3_prime4= (0.0d0,0d0)
831  complex(8), public :: gh2w3_prime5= (0.0d0,0d0)
832  complex(8), public :: gh2w3_prime6= (0.0d0,0d0)
833  complex(8), public :: gh2w3_prime7= (0.0d0,0d0)
834 
835  complex(8), public :: gh2w4_prime = (0.0d0,0d0)
836  complex(8), public :: gh2w4_prime2= (0.0d0,0d0)
837  complex(8), public :: gh2w4_prime3= (0.0d0,0d0)
838  complex(8), public :: gh2w4_prime4= (0.0d0,0d0)
839  complex(8), public :: gh2w4_prime5= (0.0d0,0d0)
840  complex(8), public :: gh2w4_prime6= (0.0d0,0d0)
841  complex(8), public :: gh2w4_prime7= (0.0d0,0d0)
842 
843  real(8), public, parameter :: lambda2_w1 = 10000d0*gev
844  real(8), public, parameter :: lambda2_w2 = 10000d0*gev
845  real(8), public, parameter :: lambda2_w3 = 10000d0*gev
846  real(8), public, parameter :: lambda2_w4 = 10000d0*gev
847 
848  integer, public :: c2w_q1sq = 0 ! Sign of q1,2,12**2 for the following Lambda2's, set to 1 or -1 to get q**2-dependence from these form factor Lambda2s
849  integer, public :: c2w_q2sq = 0
850  integer, public :: c2w_q12sq = 0
851  real(8), public :: lambda2_w11 = 100d0*gev ! For W+
852  real(8), public :: lambda2_w21 = 100d0*gev
853  real(8), public :: lambda2_w31 = 100d0*gev
854  real(8), public :: lambda2_w41 = 100d0*gev
855  real(8), public :: lambda2_w12 = 100d0*gev ! For W-
856  real(8), public :: lambda2_w22 = 100d0*gev
857  real(8), public :: lambda2_w32 = 100d0*gev
858  real(8), public :: lambda2_w42 = 100d0*gev
859  real(8), public :: lambda2_w10 = 100d0*gev ! For the Higgs
860  real(8), public :: lambda2_w20 = 100d0*gev
861  real(8), public :: lambda2_w30 = 100d0*gev
862  real(8), public :: lambda2_w40 = 100d0*gev
863 
864 
865 !-- Anomalous coupling parameters in triple and quartic gauge couplings
866 ! notation wrt paper: dV_N=d^{NWW}*d_1^N, dP_N=d^{NWW}*d_2^N, dM_N=d^{NWW}*d_3^N, dFour_N=d^{NWW}*d_4^N, for N=Z,\gamma
867 ! dZZWpWm = d^{ZZWW}*sw^2/cw^2, dZAWpWm = d^{Z\gammaWW}*sw/cw, dAAWpWm = d^{\gamma\gammaWW}
868  complex(8), public :: dv_a= (1.0d0,0d0)
869  complex(8), public :: dp_a= (1.0d0,0d0)
870  complex(8), public :: dm_a= (1.0d0,0d0)
871  complex(8), public :: dv_z= (1.0d0,0d0)
872  complex(8), public :: dp_z= (1.0d0,0d0)
873  complex(8), public :: dm_z= (1.0d0,0d0)
874  complex(8), public :: daawpwm= (1.0d0,0d0)
875  complex(8), public :: dzawpwm= (1.0d0,0d0)
876  complex(8), public :: dzzwpwm= (1.0d0,0d0)
877  complex(8), public :: dfour_a= (0.0d0,0d0)
878  complex(8), public :: dfour_z= (0.0d0,0d0)
879 
880 
881 !-- Hff couplings for ttbar+H and bbar+H
882  complex(8), public :: kappa = (1d0,0d0)
883  complex(8), public :: kappa_tilde = (0d0,0d0)
884 !-- first resonance (H) couplings
885  ! Usable in JHUGen-only processes
886  complex(8), public :: kappa_top = (1d0,0d0)
887  complex(8), public :: kappa_tilde_top = (0d0,0d0)
888  complex(8), public :: kappa_bot = (1d0,0d0)
889  complex(8), public :: kappa_tilde_bot = (0d0,0d0)
890  ! Extra MCFM couplings to fourth-generation quarks
891  complex(8), public :: kappa_4gen_top = (0d0,0d0)
892  complex(8), public :: kappa_tilde_4gen_top = (0d0,0d0)
893  complex(8), public :: kappa_4gen_bot = (0d0,0d0)
894  complex(8), public :: kappa_tilde_4gen_bot = (0d0,0d0)
895 
896 !-- seconds resonance (H2) couplings for MCFM interface
897  complex(8), public :: kappa2_top = (0d0,0d0)
898  complex(8), public :: kappa2_tilde_top = (0d0,0d0)
899  complex(8), public :: kappa2_bot = (0d0,0d0)
900  complex(8), public :: kappa2_tilde_bot = (0d0,0d0)
901 
902  complex(8), public :: kappa2_4gen_top = (0d0,0d0)
903  complex(8), public :: kappa2_tilde_4gen_top = (0d0,0d0)
904  complex(8), public :: kappa2_4gen_bot = (0d0,0d0)
905  complex(8), public :: kappa2_tilde_4gen_bot = (0d0,0d0)
906 
907 
908 !--------------------!
909 !-----! Spin-1 !-----!
910 !--------------------!
911 !---parameters that define spin 1 coupling to SM fields, see note
912  complex(8), public :: zprime_qq_left = (1.0d0,0d0)
913  complex(8), public :: zprime_qq_right = (1.0d0,0d0)
914  complex(8), public :: zprime_zz_1 = (0d0,0d0)! =1 for JP=1- vector
915  complex(8), public :: zprime_zz_2 = (0d0,0d0)! =1 for JP=1+ pseudovector
916 
917 
918 !--------------------!
919 !-----! Spin-2 !-----!
920 !--------------------!
921 !-- parameters that define spin 2 coupling to SM fields, see note
922 ! minimal coupling corresponds to a1 = b1 = b5 = 1 everything else 0
923  complex(8), public :: a1 = (0d0,0d0) ! g1 -- c.f. draft
924  complex(8), public :: a2 = (0d0,0d0) ! g2
925  complex(8), public :: a3 = (0d0,0d0) ! g3
926  complex(8), public :: a4 = (0d0,0d0) ! g4
927  complex(8), public :: a5 = (0d0,0d0) ! pseudoscalar, g8
928  complex(8), public :: graviton_qq_left = (1.0d0,0d0)! graviton coupling to quarks
929  complex(8), public :: graviton_qq_right = (1.0d0,0d0)
930 
931 !-- see mod_Graviton for these two parameters
932  logical, public, parameter :: generate_bis = .true.
933  logical, public, parameter :: use_dynamic_mg = .true.
934 
935  ! all b' below are g's in the draft
936  complex(8), public :: b1 = (0d0,0d0)
937  complex(8), public :: b2 = (0d0,0d0)
938  complex(8), public :: b3 = (0d0,0d0)
939  complex(8), public :: b4 = (0d0,0d0)
940  complex(8), public :: b5 = (0d0,0d0) ! this coupling does not contribute to V+gamma final states
941  complex(8), public :: b6 = (0d0,0d0) ! this coupling does not contribute to V+gamma final states
942  complex(8), public :: b7 = (0d0,0d0) ! this coupling does not contribute to V+gamma final states
943  complex(8), public :: b8 = (0d0,0d0)
944  complex(8), public :: b9 = (0d0,0d0) ! this coupling does not contribute to V+gamma final states
945  complex(8), public :: b10 =(0d0,0d0) ! this coupling does not contribute to V+gamma final states
946 
947  complex(8), public :: bzzp1 = (0d0,0d0)
948  complex(8), public :: bzzp2 = (0d0,0d0)
949  complex(8), public :: bzzp3 = (0d0,0d0)
950  complex(8), public :: bzzp4 = (0d0,0d0)
951  complex(8), public :: bzzp5 = (0d0,0d0)
952  complex(8), public :: bzzp6 = (0d0,0d0)
953  complex(8), public :: bzzp7 = (0d0,0d0)
954  complex(8), public :: bzzp8 = (0d0,0d0)
955  complex(8), public :: bzzp9 = (0d0,0d0)
956  complex(8), public :: bzzp10 =(0d0,0d0)
957 
958  complex(8), public :: bzpzp1 = (0d0,0d0)
959  complex(8), public :: bzpzp2 = (0d0,0d0)
960  complex(8), public :: bzpzp3 = (0d0,0d0)
961  complex(8), public :: bzpzp4 = (0d0,0d0)
962  complex(8), public :: bzpzp5 = (0d0,0d0)
963  complex(8), public :: bzpzp6 = (0d0,0d0)
964  complex(8), public :: bzpzp7 = (0d0,0d0)
965  complex(8), public :: bzpzp8 = (0d0,0d0)
966  complex(8), public :: bzpzp9 = (0d0,0d0)
967  complex(8), public :: bzpzp10 =(0d0,0d0)
968 
969  complex(8), public :: bzgs1 = (0d0,0d0)
970  complex(8), public :: bzgs2 = (0d0,0d0)
971  complex(8), public :: bzgs3 = (0d0,0d0)
972  complex(8), public :: bzgs4 = (0d0,0d0)
973  complex(8), public :: bzgs8 = (0d0,0d0)
974 
975  complex(8), public :: bgsgs1 = (0d0,0d0)
976  complex(8), public :: bgsgs2 = (0d0,0d0)
977  complex(8), public :: bgsgs3 = (0d0,0d0)
978  complex(8), public :: bgsgs4 = (0d0,0d0)
979  complex(8), public :: bgsgs8 = (0d0,0d0)
980 
981  complex(8), public :: bzpgs1 = (0d0,0d0)
982  complex(8), public :: bzpgs2 = (0d0,0d0)
983  complex(8), public :: bzpgs3 = (0d0,0d0)
984  complex(8), public :: bzpgs4 = (0d0,0d0)
985  complex(8), public :: bzpgs8 = (0d0,0d0)
986 
987 
988  complex(8), public, parameter :: c1 = (1.0d0,0d0)
989  complex(8), public, parameter :: c2 = (0d0,0d0)
990  complex(8), public, parameter :: c3 = (0d0,0d0)
991  complex(8), public, parameter :: c41= (0d0,0d0)
992  complex(8), public, parameter :: c42= (0d0,0d0)
993  complex(8), public, parameter :: c5 = (0d0,0d0)
994  complex(8), public, parameter :: c6 = (0d0,0d0) ! this coupling does not contribute to gamma+gamma final states
995  complex(8), public, parameter :: c7 = (0d0,0d0) ! this coupling does not contribute to gamma+gamma final states
996 
997 
998 ! user Zff couplings
999  real(8), public :: daz_lep_left = (0d0)
1000  real(8), public :: daz_lep_right = (0d0)
1001  real(8), public :: daz_neu_left = (0d0)
1002  real(8), public :: daz_neu_right = (0d0)
1003  real(8), public :: daz_qup_left = (0d0)
1004  real(8), public :: daz_qup_right = (0d0)
1005  real(8), public :: daz_qdn_left = (0d0)
1006  real(8), public :: daz_qdn_right = (0d0)
1007 
1008 !=====================================================
1009 
1010 !=====================================================
1011 !internal
1012 !--- B0_PDF=(11.-2.*NF/3.)/4./PI
1013 real(dp), public, parameter :: b0_pdf(0:6) = (/ 0.8753521870054244d0,0.822300539308126d0,0.7692488916108274d0,0.716197243913529d0,0.6631455962162306d0,0.6100939485189321d0,0.5570423008216338d0 /)
1014 real(dp), public, parameter :: nf = 5.0_dp
1015 real(dp), public, parameter :: xn = 3.0_dp
1016 real(dp), public, parameter :: ca = 3.0_dp
1017 real(dp), public, parameter :: cf = 4.0_dp/3.0_dp
1018 real(dp), public, parameter :: tr = 0.5_dp
1019 real(dp), public, parameter :: avegg = 1.0_dp/4.0_dp/64.0_dp
1020 real(dp), public, parameter :: aveqg = 1.0_dp/4.0_dp/24.0_dp
1021 real(dp), public, parameter :: aveqq = 1.0_dp/4.0_dp/9.0_dp
1022 
1023 ! Particle isospin and charges
1024 real(8), public, parameter :: t3ll= -0.5d0
1025 real(8), public, parameter :: t3lr= 0d0
1026 real(8), public, parameter :: t3nl= 0.5d0
1027 real(8), public, parameter :: t3nr= 0d0
1028 real(8), public, parameter :: t3ul= 0.5d0
1029 real(8), public, parameter :: t3ur= 0d0
1030 real(8), public, parameter :: t3dl= -0.5d0
1031 real(8), public, parameter :: t3dr= 0d0
1032 real(8), public, parameter :: qll = -1d0
1033 real(8), public, parameter :: qlr = -1d0
1034 real(8), public, parameter :: qnl = 0d0
1035 real(8), public, parameter :: qnr = 0d0
1036 real(8), public, parameter :: qul = 2d0/3d0
1037 real(8), public, parameter :: qur = 2d0/3d0
1038 real(8), public, parameter :: qdl = -1d0/3d0
1039 real(8), public, parameter :: qdr = -1d0/3d0
1040 
1041 
1042 ! V-f-fbar couplings:
1043 ! g_R(f) = -e*sw/cw*Q(f) = e/2/sw/cw * a(b,c)R,
1044 ! g_L(f) = -e*sw/cw*Q(f) + e/sw/cw*T3(f) = e/2/sw/cw * a(b,c)L
1045 ! with
1046 ! aR(f) = -2*sw**2*Q(f),
1047 ! aL(f) = -2*sw**2*Q(f) + 2*T3(f).
1048 ! for V = Z-boson,
1049 ! bR = 0
1050 ! bL = dsqrt(2)*cw
1051 ! for V = W-boson,
1052 ! and
1053 ! cR = -2*sw*cw*Q(f)
1054 ! cL = -2*sw*cw*Q(f)
1055 ! for V = photon*
1056 !
1057 real(8), public :: overallcouplvffsq ! Overall coupling squared that goes along with the ones below
1058 real(8), public :: ar_lep ! = -2d0*sitW**2*(-1d0)
1059 real(8), public :: al_lep ! = -2d0*sitW**2*(-1d0)-1d0
1060 real(8), public :: ar_neu ! = -2d0*sitW**2*(0d0)
1061 real(8), public :: al_neu ! = -2d0*sitW**2*(0d0)+1d0
1062 real(8), public :: ar_qup ! = -2d0*sitW**2*(2d0/3d0)
1063 real(8), public :: al_qup ! = -2d0*sitW**2*(2d0/3d0)+1d0
1064 real(8), public :: ar_qdn ! = -2d0*sitW**2*(-1d0/3d0)
1065 real(8), public :: al_qdn ! = -2d0*sitW**2*(-1d0/3d0)-1d0
1066 real(8), public :: bl ! = dsqrt(2d0)*dsqrt(1d0-sitW**2)
1067 real(8), public :: br ! = 0d0
1068 real(8), public :: cr_lep ! = -2d0*sitW*dsqrt(1d0-sitW**2)*(-1d0)
1069 real(8), public :: cl_lep ! = -2d0*sitW*dsqrt(1d0-sitW**2)*(-1d0)
1070 real(8), public :: cr_neu ! = -2d0*sitW*dsqrt(1d0-sitW**2)*(0d0)
1071 real(8), public :: cl_neu ! = -2d0*sitW*dsqrt(1d0-sitW**2)*(0d0)
1072 real(8), public :: cr_qup ! = -2d0*sitW*dsqrt(1d0-sitW**2)*(2d0/3d0)
1073 real(8), public :: cl_qup ! = -2d0*sitW*dsqrt(1d0-sitW**2)*(2d0/3d0)
1074 real(8), public :: cr_qdn ! = -2d0*sitW*dsqrt(1d0-sitW**2)*(-1d0/3d0)
1075 real(8), public :: cl_qdn ! = -2d0*sitW*dsqrt(1d0-sitW**2)*(-1d0/3d0)
1076 
1077 ! Coupling normalizations based on the common factor to be gwsq
1078 real(dp), public :: couplwffsq
1079 real(dp), public :: couplzffsq
1080 real(dp), public :: couplazff
1081 real(dp), public :: couplaffsq
1082 !----------------------------------------------------------------------------------
1083 
1084 integer, public, target :: up_ = 1
1085 integer, public, target :: dn_ = 2
1086 integer, public, target :: chm_ = 3
1087 integer, public, target :: str_ = 4
1088 integer, public, target :: top_ = 5
1089 integer, public, target :: bot_ = 6
1090 integer, public, target :: elp_ = 7
1091 integer, public, target :: mup_ = 8
1092 integer, public, target :: tap_ = 9
1093 integer, public, target :: glu_ = 10
1094 integer, public, target :: pho_ = 11
1095 integer, public, target :: z0_ = 12
1096 integer, public, target :: wp_ = 13
1097 integer, public, target :: nue_ = 14
1098 integer, public, target :: num_ = 15
1099 integer, public, target :: nut_ = 16
1100 integer, public, target :: hig_ = 25
1101 integer, public, target :: zpr_ = 32
1102 integer, public, target :: zpr2_ = 33
1103 integer, public, target :: wppr_ = 34
1104 integer, public, target :: gra_ = 39
1105 
1106 integer, public, target :: aup_ = -1
1107 integer, public, target :: adn_ = -2
1108 integer, public, target :: achm_ = -3
1109 integer, public, target :: astr_ = -4
1110 integer, public, target :: atop_ = -5
1111 integer, public, target :: abot_ = -6
1112 integer, public, target :: elm_ = -7
1113 integer, public, target :: mum_ = -8
1114 integer, public, target :: tam_ = -9
1115 integer, public, target :: wm_ = -13
1116 integer, public, target :: anue_ = -14
1117 integer, public, target :: anum_ = -15
1118 integer, public, target :: anut_ = -16
1119 integer, public, target :: wmpr_ = -34
1120 
1121 integer, public, parameter :: not_a_particle_ = -9000
1122 real(8), public, parameter :: mom_not_a_particle(1:4) = (/0d0,0d0,0d0,0d0/)
1123 
1124 integer, public, parameter :: pdfglu_ = 0
1125 integer, public, parameter :: pdfdn_ = 1
1126 integer, public, parameter :: pdfup_ = 2
1127 integer, public, parameter :: pdfstr_ = 3
1128 integer, public, parameter :: pdfchm_ = 4
1129 integer, public, parameter :: pdfbot_ = 5
1130 integer, public, parameter :: pdftop_ = 6 ! Dummy
1131 integer, public, parameter :: pdfadn_ = -1
1132 integer, public, parameter :: pdfaup_ = -2
1133 integer, public, parameter :: pdfastr_ = -3
1134 integer, public, parameter :: pdfachm_ = -4
1135 integer, public, parameter :: pdfabot_ = -5
1136 integer, public, parameter :: pdfatop_ = -6 ! Dummy
1137 
1138 integer,parameter :: io_stdout=6
1139 integer,parameter :: io_lheoutfile=14
1140 integer,parameter :: io_histofile=15
1141 integer,parameter :: io_lheinfile=16
1142 integer,parameter :: io_logfile=17
1143 integer,parameter :: io_csmaxfile=18
1144 integer,parameter :: io_lheoutfile2=19
1145 integer,parameter :: io_lheoutfile3=20
1146 integer,parameter :: io_tmpfile=21 !to use for whatever purpose is needed, but close afterwards
1147 
1148 integer, public :: debugcounter(0:10) = 0
1149 real(8), public :: debugvar(0:10) = 0d0
1150 
1151 
1152 integer, public :: ijpartons(1:2)=0
1153 
1154 !=====================================================
1155 
1156 
1157 CONTAINS
1158 
1159 
1160 function hvvspinzerodynamiccoupling (index,sWplus,sWminus,sWW,tryWWcoupl)
1161 integer, intent(in) :: index
1162 real(8), intent(in) :: swplus, swminus, sww
1163 real(8) :: swplus_signed, swminus_signed, sww_signed, qsqcompoundfactor
1164 logical,optional :: trywwcoupl
1165 complex(8) :: hvvspinzerodynamiccoupling
1166 complex(8) :: vvcoupl(1:8)
1167 real(8) :: lambda_v
1168 real(8) :: lambda_v120(1:3)
1169 logical :: forcezzcoupl
1170 logical :: computeqsqcompundcoupl
1171 
1172  if(present(trywwcoupl)) then
1173  forcezzcoupl = (.not.trywwcoupl .or. .not.distinguish_hwwcouplings .or. (index.gt.4 .and. index.lt.12))
1174  else
1175  forcezzcoupl = .true.
1176  endif
1177  computeqsqcompundcoupl = .false.
1178  swplus_signed=0d0
1179  swminus_signed=0d0
1180  sww_signed=0d0
1181  vvcoupl(:)=czero
1182  hvvspinzerodynamiccoupling=czero
1183  if( forcezzcoupl ) then
1184  if(cz_q1sq.ne.0) swplus_signed=abs(swplus)*dble(sign(1,cz_q1sq))
1185  if(cz_q2sq.ne.0) swminus_signed=abs(swminus)*dble(sign(1,cz_q2sq))
1186  if(cz_q12sq.ne.0) sww_signed=abs(sww)*dble(sign(1,cz_q12sq))
1187  if(cz_q1sq.ne.0 .or. cz_q2sq.ne.0 .or. cz_q12sq.ne.0) computeqsqcompundcoupl=.true.
1188  if(index.eq.1) then ! ZZ 1-4
1189  vvcoupl = (/ ghz1, ghz1_prime, ghz1_prime2, ghz1_prime3, ghz1_prime4, ghz1_prime5, ghz1_prime6, ghz1_prime7 /)
1190  lambda_v = lambda_z1
1191  lambda_v120 = (/ lambda_z11, lambda_z12, lambda_z10 /)
1192  elseif(index.eq.2) then
1193  vvcoupl = (/ ghz2, ghz2_prime, ghz2_prime2, ghz2_prime3, ghz2_prime4, ghz2_prime5, ghz2_prime6, ghz2_prime7 /)
1194  lambda_v = lambda_z2
1195  lambda_v120 = (/ lambda_z21, lambda_z22, lambda_z20 /)
1196  elseif(index.eq.3) then
1197  vvcoupl = (/ ghz3, ghz3_prime, ghz3_prime2, ghz3_prime3, ghz3_prime4, ghz3_prime5, ghz3_prime6, ghz3_prime7 /)
1198  lambda_v = lambda_z3
1199  lambda_v120 = (/ lambda_z31, lambda_z32, lambda_z30 /)
1200  elseif(index.eq.4) then
1201  vvcoupl = (/ ghz4, ghz4_prime, ghz4_prime2, ghz4_prime3, ghz4_prime4, ghz4_prime5, ghz4_prime6, ghz4_prime7 /)
1202  lambda_v = lambda_z4
1203  lambda_v120 = (/ lambda_z41, lambda_z42, lambda_z40 /)
1204  elseif(index.eq.5) then ! Zgs 1
1205  vvcoupl(3) = ghzgs1_prime2
1206  lambda_v = lambda_zgs1
1207  lambda_v120 = (/ lambda_z11, lambda_z12, lambda_z10 /)
1208  elseif(index.eq.6) then ! Zgs 2-4
1209  vvcoupl(1) = ghzgs2
1210  lambda_v = 1d0 ! Not present
1211  lambda_v120 = (/ lambda_z21, lambda_z22, lambda_z20 /)
1212  elseif(index.eq.7) then
1213  vvcoupl(1) = ghzgs3
1214  lambda_v = 1d0 ! Not present
1215  lambda_v120 = (/ lambda_z31, lambda_z32, lambda_z30 /)
1216  elseif(index.eq.8) then
1217  vvcoupl(1) = ghzgs4
1218  lambda_v = 1d0 ! Not present
1219  lambda_v120 = (/ lambda_z41, lambda_z42, lambda_z40 /)
1220  elseif(index.eq.9) then ! gsgs 2-4
1221  vvcoupl(1) = ghgsgs2
1222  lambda_v = 1d0 ! Not present
1223  lambda_v120 = (/ lambda_z21, lambda_z22, lambda_z20 /)
1224  elseif(index.eq.10) then
1225  vvcoupl(1) = ghgsgs3
1226  lambda_v = 1d0 ! Not present
1227  lambda_v120 = (/ lambda_z31, lambda_z32, lambda_z30 /)
1228  elseif(index.eq.11) then
1229  vvcoupl(1) = ghgsgs4
1230  lambda_v = 1d0 ! Not present
1231  lambda_v120 = (/ lambda_z41, lambda_z42, lambda_z40 /)
1232  elseif(index.eq.12) then ! ZpZ 1-4
1233  vvcoupl = (/ ghzzp1, ghzzp1_prime, ghzzp1_prime2, ghzzp1_prime3, ghzzp1_prime4, ghzzp1_prime5, ghzzp1_prime6, ghzzp1_prime7 /)
1234  lambda_v = lambda_z1
1235  lambda_v120 = (/ lambda_z11, lambda_z12, lambda_z10 /)
1236  elseif(index.eq.13) then
1237  vvcoupl = (/ ghzzp2, ghzzp2_prime, ghzzp2_prime2, ghzzp2_prime3, ghzzp2_prime4, ghzzp2_prime5, ghzzp2_prime6, ghzzp2_prime7 /)
1238  lambda_v = lambda_z2
1239  lambda_v120 = (/ lambda_z11, lambda_z12, lambda_z10 /)
1240  elseif(index.eq.14) then
1241  vvcoupl = (/ ghzzp3, ghzzp3_prime, ghzzp3_prime2, ghzzp3_prime3, ghzzp3_prime4, ghzzp3_prime5, ghzzp3_prime6, ghzzp3_prime7 /)
1242  lambda_v = lambda_z3
1243  lambda_v120 = (/ lambda_z11, lambda_z12, lambda_z10 /)
1244  elseif(index.eq.15) then
1245  vvcoupl = (/ ghzzp4, ghzzp4_prime, ghzzp4_prime2, ghzzp4_prime3, ghzzp4_prime4, ghzzp4_prime5, ghzzp4_prime6, ghzzp4_prime7 /)
1246  lambda_v = lambda_z4
1247  lambda_v120 = (/ lambda_z11, lambda_z12, lambda_z10 /)
1248  elseif(index.eq.16) then ! ZpZp 1-4
1249  vvcoupl = (/ ghzpzp1, ghzpzp1_prime, ghzpzp1_prime2, ghzpzp1_prime3, ghzpzp1_prime4, ghzpzp1_prime5, ghzpzp1_prime6, ghzpzp1_prime7 /)
1250  lambda_v = lambda_z1
1251  lambda_v120 = (/ lambda_z11, lambda_z12, lambda_z10 /)
1252  elseif(index.eq.17) then
1253  vvcoupl = (/ ghzpzp2, ghzpzp2_prime, ghzpzp2_prime2, ghzpzp2_prime3, ghzpzp2_prime4, ghzpzp2_prime5, ghzpzp2_prime6, ghzpzp2_prime7 /)
1254  lambda_v = lambda_z2
1255  lambda_v120 = (/ lambda_z11, lambda_z12, lambda_z10 /)
1256  elseif(index.eq.18) then
1257  vvcoupl = (/ ghzpzp3, ghzpzp3_prime, ghzpzp3_prime2, ghzpzp3_prime3, ghzpzp3_prime4, ghzpzp3_prime5, ghzpzp3_prime6, ghzpzp3_prime7 /)
1258  lambda_v = lambda_z3
1259  lambda_v120 = (/ lambda_z11, lambda_z12, lambda_z10 /)
1260  elseif(index.eq.19) then
1261  vvcoupl = (/ ghzpzp4, ghzpzp4_prime, ghzpzp4_prime2, ghzpzp4_prime3, ghzpzp4_prime4, ghzpzp4_prime5, ghzpzp4_prime6, ghzpzp4_prime7 /)
1262  lambda_v = lambda_z4
1263  lambda_v120 = (/ lambda_z11, lambda_z12, lambda_z10 /)
1264  elseif(index.eq.20) then ! Zpgs 1
1265  vvcoupl(3) = ghzpgs1_prime2
1266  lambda_v = lambda_zgs1
1267  lambda_v120 = (/ lambda_z11, lambda_z12, lambda_z10 /)
1268  elseif(index.eq.21) then ! Zpgs 2-4
1269  vvcoupl(1) = ghzpgs2
1270  lambda_v = 1d0 ! Not present
1271  lambda_v120 = (/ lambda_z21, lambda_z22, lambda_z20 /)
1272  elseif(index.eq.22) then
1273  vvcoupl(1) = ghzpgs3
1274  lambda_v = 1d0 ! Not present
1275  lambda_v120 = (/ lambda_z31, lambda_z32, lambda_z30 /)
1276  elseif(index.eq.23) then
1277  vvcoupl(1) = ghzpgs4
1278  lambda_v = 1d0 ! Not present
1279  lambda_v120 = (/ lambda_z41, lambda_z42, lambda_z40 /)
1280  endif
1281  else
1282  if(cw_q1sq.ne.0) swplus_signed=abs(swplus)*dble(sign(1,cw_q1sq))
1283  if(cw_q2sq.ne.0) swminus_signed=abs(swminus)*dble(sign(1,cw_q2sq))
1284  if(cw_q12sq.ne.0) sww_signed=abs(sww)*dble(sign(1,cw_q12sq))
1285  if(cw_q1sq.ne.0 .or. cw_q2sq.ne.0 .or. cw_q12sq.ne.0) computeqsqcompundcoupl=.true.
1286  if(index.eq.1) then
1287  vvcoupl = (/ ghw1, ghw1_prime, ghw1_prime2, ghw1_prime3, ghw1_prime4, ghw1_prime5, ghw1_prime6, ghw1_prime7 /)
1288  lambda_v = lambda_w1
1289  lambda_v120 = (/ lambda_w11, lambda_w12, lambda_w10 /)
1290  elseif(index.eq.2) then
1291  vvcoupl = (/ ghw2, ghw2_prime, ghw2_prime2, ghw2_prime3, ghw2_prime4, ghw2_prime5, ghw2_prime6, ghw2_prime7 /)
1292  lambda_v = lambda_w2
1293  lambda_v120 = (/ lambda_w21, lambda_w22, lambda_w20 /)
1294  elseif(index.eq.3) then
1295  vvcoupl = (/ ghw3, ghw3_prime, ghw3_prime2, ghw3_prime3, ghw3_prime4, ghw3_prime5, ghw3_prime6, ghw3_prime7 /)
1296  lambda_v = lambda_w3
1297  lambda_v120 = (/ lambda_w31, lambda_w32, lambda_w30 /)
1298  elseif(index.eq.4) then
1299  vvcoupl = (/ ghw4, ghw4_prime, ghw4_prime2, ghw4_prime3, ghw4_prime4, ghw4_prime5, ghw4_prime6, ghw4_prime7 /)
1300  lambda_v = lambda_w4
1301  lambda_v120 = (/ lambda_w41, lambda_w42, lambda_w40 /)
1302  elseif(index.eq.12) then
1303  vvcoupl = (/ ghwwp1, ghwwp1_prime, ghwwp1_prime2, ghwwp1_prime3, ghwwp1_prime4, ghwwp1_prime5, ghwwp1_prime6, ghwwp1_prime7 /)
1304  lambda_v = lambda_w1
1305  lambda_v120 = (/ lambda_w11, lambda_w12, lambda_w10 /)
1306  elseif(index.eq.13) then
1307  vvcoupl = (/ ghwwp2, ghwwp2_prime, ghwwp2_prime2, ghwwp2_prime3, ghwwp2_prime4, ghwwp2_prime5, ghwwp2_prime6, ghwwp2_prime7 /)
1308  lambda_v = lambda_w2
1309  lambda_v120 = (/ lambda_w11, lambda_w12, lambda_w10 /)
1310  elseif(index.eq.14) then
1311  vvcoupl = (/ ghwwp3, ghwwp3_prime, ghwwp3_prime2, ghwwp3_prime3, ghwwp3_prime4, ghwwp3_prime5, ghwwp3_prime6, ghwwp3_prime7 /)
1312  lambda_v = lambda_w3
1313  lambda_v120 = (/ lambda_w11, lambda_w12, lambda_w10 /)
1314  elseif(index.eq.15) then
1315  vvcoupl = (/ ghwwp4, ghwwp4_prime, ghwwp4_prime2, ghwwp4_prime3, ghwwp4_prime4, ghwwp4_prime5, ghwwp4_prime6, ghwwp4_prime7 /)
1316  lambda_v = lambda_w4
1317  lambda_v120 = (/ lambda_w11, lambda_w12, lambda_w10 /)
1318  elseif(index.eq.16) then
1319  vvcoupl = (/ ghwpwp1, ghwpwp1_prime, ghwpwp1_prime2, ghwpwp1_prime3, ghwpwp1_prime4, ghwpwp1_prime5, ghwpwp1_prime6, ghwpwp1_prime7 /)
1320  lambda_v = lambda_w1
1321  lambda_v120 = (/ lambda_w11, lambda_w12, lambda_w10 /)
1322  elseif(index.eq.17) then
1323  vvcoupl = (/ ghwpwp2, ghwpwp2_prime, ghwpwp2_prime2, ghwpwp2_prime3, ghwpwp2_prime4, ghwpwp2_prime5, ghwpwp2_prime6, ghwpwp2_prime7 /)
1324  lambda_v = lambda_w2
1325  lambda_v120 = (/ lambda_w11, lambda_w12, lambda_w10 /)
1326  elseif(index.eq.18) then
1327  vvcoupl = (/ ghwpwp3, ghwpwp3_prime, ghwpwp3_prime2, ghwpwp3_prime3, ghwpwp3_prime4, ghwpwp3_prime5, ghwpwp3_prime6, ghwpwp3_prime7 /)
1328  lambda_v = lambda_w3
1329  lambda_v120 = (/ lambda_w11, lambda_w12, lambda_w10 /)
1330  elseif(index.eq.19) then
1331  vvcoupl = (/ ghwpwp4, ghwpwp4_prime, ghwpwp4_prime2, ghwpwp4_prime3, ghwpwp4_prime4, ghwpwp4_prime5, ghwpwp4_prime6, ghwpwp4_prime7 /)
1332  lambda_v = lambda_w4
1333  lambda_v120 = (/ lambda_w11, lambda_w12, lambda_w10 /)
1334  endif
1335  endif
1336 
1337  if(vvcoupl(2).ne.czero) hvvspinzerodynamiccoupling = hvvspinzerodynamiccoupling + vvcoupl(2) * lambda_v**4/(lambda_v**2 + abs(swplus))/(lambda_v**2 + abs(swminus))
1338  if(vvcoupl(3).ne.czero) hvvspinzerodynamiccoupling = hvvspinzerodynamiccoupling + vvcoupl(3) * ( swplus + swminus )/lambda_v**2
1339  if(vvcoupl(4).ne.czero) hvvspinzerodynamiccoupling = hvvspinzerodynamiccoupling + vvcoupl(4) * ( swplus - swminus )/lambda_v**2
1340  if(vvcoupl(5).ne.czero) hvvspinzerodynamiccoupling = hvvspinzerodynamiccoupling + vvcoupl(5) * ( sww )/lambda_q**2
1341  if(vvcoupl(6).ne.czero) hvvspinzerodynamiccoupling = hvvspinzerodynamiccoupling + vvcoupl(6) * ( swplus**2 + swminus**2 )/lambda_v**4
1342  if(vvcoupl(7).ne.czero) hvvspinzerodynamiccoupling = hvvspinzerodynamiccoupling + vvcoupl(7) * ( swplus**2 - swminus**2 )/lambda_v**4
1343  if(vvcoupl(8).ne.czero) hvvspinzerodynamiccoupling = hvvspinzerodynamiccoupling + vvcoupl(8) * ( swplus * swminus )/lambda_v**4
1344 
1345  if(index.eq.1) then
1346  if(computeqsqcompundcoupl) then
1347  qsqcompoundfactor = (lambda_v120(1)**2 + swplus_signed)*(lambda_v120(2)**2 + swminus_signed)*(lambda_v120(3)**2 + sww_signed)
1348  if(qsqcompoundfactor.ne.0d0) then
1349  qsqcompoundfactor = (lambda_v120(1)*lambda_v120(2)*lambda_v120(3))**2/qsqcompoundfactor
1350  endif
1351  hvvspinzerodynamiccoupling = hvvspinzerodynamiccoupling * qsqcompoundfactor
1352  endif
1353  if(vvcoupl(1).ne.czero) hvvspinzerodynamiccoupling = hvvspinzerodynamiccoupling + vvcoupl(1)
1354  else
1355  if(vvcoupl(1).ne.czero) hvvspinzerodynamiccoupling = hvvspinzerodynamiccoupling + vvcoupl(1)
1356  if(computeqsqcompundcoupl) then
1357  qsqcompoundfactor = (lambda_v120(1)**2 + swplus_signed)*(lambda_v120(2)**2 + swminus_signed)*(lambda_v120(3)**2 + sww_signed)
1358  if(qsqcompoundfactor.ne.0d0) then
1359  qsqcompoundfactor = (lambda_v120(1)*lambda_v120(2)*lambda_v120(3))**2/qsqcompoundfactor
1360  endif
1361  hvvspinzerodynamiccoupling = hvvspinzerodynamiccoupling * qsqcompoundfactor
1362  endif
1363  endif
1364 
1365 return
1366 end function
1367 
1368 
1369 function vpffcoupling(jhuid, hel, useWp)
1370 integer, intent(in) :: jhuid
1371 integer, intent(in) :: hel
1372 logical, intent(in) :: usewp
1373 complex(8) :: vpffcoupling
1374 integer :: absid
1375 
1376  vpffcoupling=czero
1377  if (abs(hel).ne.1) then
1378  return
1379  endif
1380  absid=abs(jhuid)
1381 
1382  if(usewp) then
1383  if (absid.eq.abs(up_) .or. absid.eq.abs(dn_)) then
1384  if(hel.eq.-1) then
1385  vpffcoupling=ewp_up_left
1386  else
1387  vpffcoupling=ewp_up_right
1388  endif
1389  elseif (absid.eq.abs(chm_) .or. absid.eq.abs(str_)) then
1390  if(hel.eq.-1) then
1391  vpffcoupling=ewp_chm_left
1392  else
1393  vpffcoupling=ewp_chm_right
1394  endif
1395  elseif (absid.eq.abs(top_) .or. absid.eq.abs(bot_)) then
1396  if(hel.eq.-1) then
1397  vpffcoupling=ewp_top_left
1398  else
1399  vpffcoupling=ewp_top_right
1400  endif
1401  elseif (absid.eq.abs(elp_) .or. absid.eq.abs(nue_)) then
1402  if(hel.eq.-1) then
1403  vpffcoupling=ewp_el_left
1404  else
1405  vpffcoupling=ewp_el_right
1406  endif
1407  elseif (absid.eq.abs(mup_) .or. absid.eq.abs(num_)) then
1408  if(hel.eq.-1) then
1409  vpffcoupling=ewp_mu_left
1410  else
1411  vpffcoupling=ewp_mu_right
1412  endif
1413  elseif (absid.eq.abs(tap_) .or. absid.eq.abs(nut_)) then
1414  if(hel.eq.-1) then
1415  vpffcoupling=ewp_ta_left
1416  else
1417  vpffcoupling=ewp_ta_right
1418  endif
1419  endif
1420  else
1421  if (absid.eq.abs(up_)) then
1422  if(hel.eq.-1) then
1423  vpffcoupling=ezp_up_left
1424  else
1425  vpffcoupling=ezp_up_right
1426  endif
1427  elseif (absid.eq.abs(dn_)) then
1428  if(hel.eq.-1) then
1429  vpffcoupling=ezp_dn_left
1430  else
1431  vpffcoupling=ezp_dn_right
1432  endif
1433  elseif (absid.eq.abs(chm_)) then
1434  if(hel.eq.-1) then
1435  vpffcoupling=ezp_chm_left
1436  else
1437  vpffcoupling=ezp_chm_right
1438  endif
1439  elseif (absid.eq.abs(str_)) then
1440  if(hel.eq.-1) then
1441  vpffcoupling=ezp_str_left
1442  else
1443  vpffcoupling=ezp_str_right
1444  endif
1445  elseif (absid.eq.abs(top_)) then
1446  if(hel.eq.-1) then
1447  vpffcoupling=ezp_top_left
1448  else
1449  vpffcoupling=ezp_top_right
1450  endif
1451  elseif (absid.eq.abs(bot_)) then
1452  if(hel.eq.-1) then
1453  vpffcoupling=ezp_bot_left
1454  else
1455  vpffcoupling=ezp_bot_right
1456  endif
1457  elseif (absid.eq.abs(elp_)) then
1458  if(hel.eq.-1) then
1459  vpffcoupling=ezp_el_left
1460  else
1461  vpffcoupling=ezp_el_right
1462  endif
1463  elseif (absid.eq.abs(nue_)) then
1464  if(hel.eq.-1) then
1465  vpffcoupling=ezp_nue_left
1466  else
1467  vpffcoupling=ezp_nue_right
1468  endif
1469  elseif (absid.eq.abs(mup_)) then
1470  if(hel.eq.-1) then
1471  vpffcoupling=ezp_mu_left
1472  else
1473  vpffcoupling=ezp_mu_right
1474  endif
1475  elseif (absid.eq.abs(num_)) then
1476  if(hel.eq.-1) then
1477  vpffcoupling=ezp_nue_left
1478  else
1479  vpffcoupling=ezp_nue_right
1480  endif
1481  elseif (absid.eq.abs(tap_)) then
1482  if(hel.eq.-1) then
1483  vpffcoupling=ezp_ta_left
1484  else
1485  vpffcoupling=ezp_ta_right
1486  endif
1487  elseif (absid.eq.abs(nut_)) then
1488  if(hel.eq.-1) then
1489  vpffcoupling=ezp_nue_left
1490  else
1491  vpffcoupling=ezp_nue_right
1492  endif
1493  endif
1494  endif
1495 
1496 end function
1497 function vpffcoupling_pdg(pdgid, hel, useWp)
1498 integer, intent(in) :: pdgid
1499 integer, intent(in) :: hel
1500 logical, intent(in) :: usewp
1501 complex(8) :: vpffcoupling_pdg
1502 integer :: jhuid
1503  jhuid=convertlhereverse(pdgid)
1504  vpffcoupling_pdg=vpffcoupling(jhuid,hel,usewp)
1505 end function
1506 
1507 
1508 FUNCTION scalefactor(id1in,id2in)
1509 implicit none
1510 real(8) :: scalefactor
1511 integer :: id1, id2, id1in, id2in
1512 id1 = abs(id1in)
1513 id2 = abs(id2in)
1514 
1515 ! W->qq
1516 if((id1.eq.convertlhe(up_) .and. id2.eq.convertlhe(dn_)) .or. (id1.eq.convertlhe(dn_) .and. id2.eq.convertlhe(up_)))then
1517  scalefactor = scale_alpha_w_ud
1518 elseif((id1.eq.convertlhe(up_) .and. id2.eq.convertlhe(str_)) .or. (id1.eq.convertlhe(str_) .and. id2.eq.convertlhe(up_)))then
1519  scalefactor = scale_alpha_w_ud
1520 elseif((id1.eq.convertlhe(up_) .and. id2.eq.convertlhe(bot_)) .or. (id1.eq.convertlhe(bot_) .and. id2.eq.convertlhe(up_)))then
1521  scalefactor = scale_alpha_w_ud
1522 elseif((id1.eq.convertlhe(chm_) .and. id2.eq.convertlhe(dn_)) .or. (id1.eq.convertlhe(dn_) .and. id2.eq.convertlhe(chm_)))then
1523  scalefactor = scale_alpha_w_cs
1524 elseif((id1.eq.convertlhe(chm_) .and. id2.eq.convertlhe(str_)) .or. (id1.eq.convertlhe(str_) .and. id2.eq.convertlhe(chm_)))then
1525  scalefactor = scale_alpha_w_cs
1526 elseif((id1.eq.convertlhe(chm_) .and. id2.eq.convertlhe(bot_)) .or. (id1.eq.convertlhe(bot_) .and. id2.eq.convertlhe(chm_)))then
1527  scalefactor = scale_alpha_w_cs
1528 ! W-> td
1529 elseif((id1.eq.convertlhe(top_) .and. id2.eq.convertlhe(dn_)) .or. (id1.eq.convertlhe(dn_) .and. id2.eq.convertlhe(top_)))then
1530  scalefactor = 1d0
1531 elseif((id1.eq.convertlhe(top_) .and. id2.eq.convertlhe(str_)) .or. (id1.eq.convertlhe(str_) .and. id2.eq.convertlhe(top_)))then
1532  scalefactor = 1d0
1533 elseif((id1.eq.convertlhe(top_) .and. id2.eq.convertlhe(bot_)) .or. (id1.eq.convertlhe(bot_) .and. id2.eq.convertlhe(top_)))then
1534  scalefactor = 1d0
1535 ! W->lnu
1536 elseif((abs(id1).eq.abs(convertlhe(nut_)) .and. abs(id2).eq.abs(convertlhe(tap_))) .or. (abs(id1).eq.abs(convertlhe(tap_)) .and. abs(id2).eq.abs(convertlhe(nut_))))then
1537  scalefactor = scale_alpha_w_tn
1538 elseif((abs(id1).eq.abs(convertlhe(num_)) .and. abs(id2).eq.abs(convertlhe(mup_))) .or. (abs(id1).eq.abs(convertlhe(mup_)) .and. abs(id2).eq.abs(convertlhe(num_))))then
1539  scalefactor = scale_alpha_w_ln
1540 elseif((abs(id1).eq.abs(convertlhe(nue_)) .and. abs(id2).eq.abs(convertlhe(elp_))) .or. (abs(id1).eq.abs(convertlhe(elp_)) .and. abs(id2).eq.abs(convertlhe(nue_))))then
1541  scalefactor = scale_alpha_w_ln
1542 ! Z->qq
1543 elseif((id1.eq.convertlhe(up_) .and. id2.eq.convertlhe(up_)) .or. (id1.eq.convertlhe(chm_) .and. id2.eq.convertlhe(chm_)))then
1544  scalefactor = scale_alpha_z_uu
1545 elseif((id1.eq.convertlhe(dn_) .and. id2.eq.convertlhe(dn_)) .or. (id1.eq.convertlhe(str_) .and. id2.eq.convertlhe(str_)) .or. (id1.eq.convertlhe(bot_) .and. id2.eq.convertlhe(bot_)))then
1546  scalefactor = scale_alpha_z_dd
1547 ! Z-> ll, nunu
1548 elseif(id1.eq.convertlhe(tap_) .and. id2.eq.convertlhe(tap_))then
1549  scalefactor = scale_alpha_z_tt
1550 elseif((id1.eq.convertlhe(mup_) .and. id2.eq.convertlhe(mup_)) .or. (id1.eq.convertlhe(elp_) .and. id2.eq.convertlhe(elp_)))then
1551  scalefactor = scale_alpha_z_ll
1552 elseif((id1.eq.convertlhe(nut_) .and. id2.eq.convertlhe(nut_)) .or. (id1.eq.convertlhe(num_) .and. id2.eq.convertlhe(num_)) .or. (id1.eq.convertlhe(nue_) .and. id2.eq.convertlhe(nue_)))then
1553  scalefactor = scale_alpha_z_nn
1554 ! Everything else
1555 else
1556  scalefactor = 1d0
1557 endif
1558 
1559 return
1560 END FUNCTION
1561 
1562 FUNCTION ckmbare(id1in,id2in)
1563 implicit none
1564 real(8) :: ckmbare
1565 integer :: id1, id2, id1in, id2in
1566 id1 = abs(id1in)
1567 id2 = abs(id2in)
1568 if((id1.eq.convertlhe(up_) .and. id2.eq.convertlhe(dn_)) .or. (id1.eq.convertlhe(dn_) .and. id2.eq.convertlhe(up_)))then
1569  ckmbare= vckm_ud
1570 elseif((id1.eq.convertlhe(up_) .and. id2.eq.convertlhe(str_)) .or. (id1.eq.convertlhe(str_) .and. id2.eq.convertlhe(up_)))then
1571  ckmbare= vckm_us
1572 elseif((id1.eq.convertlhe(up_) .and. id2.eq.convertlhe(bot_)) .or. (id1.eq.convertlhe(bot_) .and. id2.eq.convertlhe(up_)))then
1573  ckmbare= vckm_ub
1574 elseif((id1.eq.convertlhe(chm_) .and. id2.eq.convertlhe(dn_)) .or. (id1.eq.convertlhe(dn_) .and. id2.eq.convertlhe(chm_)))then
1575  ckmbare= vckm_cd
1576 elseif((id1.eq.convertlhe(chm_) .and. id2.eq.convertlhe(str_)) .or. (id1.eq.convertlhe(str_) .and. id2.eq.convertlhe(chm_)))then
1577  ckmbare= vckm_cs
1578 elseif((id1.eq.convertlhe(chm_) .and. id2.eq.convertlhe(bot_)) .or. (id1.eq.convertlhe(bot_) .and. id2.eq.convertlhe(chm_)))then
1579  ckmbare= vckm_cb
1580 elseif((id1.eq.convertlhe(top_) .and. id2.eq.convertlhe(dn_)) .or. (id1.eq.convertlhe(dn_) .and. id2.eq.convertlhe(top_)))then
1581  ckmbare= vckm_td
1582 elseif((id1.eq.convertlhe(top_) .and. id2.eq.convertlhe(str_)) .or. (id1.eq.convertlhe(str_) .and. id2.eq.convertlhe(top_)))then
1583  ckmbare= vckm_ts
1584 elseif((id1.eq.convertlhe(top_) .and. id2.eq.convertlhe(bot_)) .or. (id1.eq.convertlhe(bot_) .and. id2.eq.convertlhe(top_)))then
1585  ckmbare= vckm_tb
1586 else
1587  ckmbare= 1d0
1588 endif
1589 
1590 END FUNCTION
1591 
1592 FUNCTION ckm(id1in,id2in)
1593 implicit none
1594 real(8) :: ckm
1595 integer :: id1in, id2in
1596  ckm=ckmbare(id1in, id2in)*sqrt(scalefactor(id1in,id2in))
1597 END FUNCTION
1598 
1599 ! for VH use
1600 FUNCTION zffbare(id1in,id2in,h1,h2)
1601 implicit none
1602 real(8) :: zffbare
1603 integer :: id1, id2, id1in, id2in
1604 real(8) :: h1,h2
1605 
1606 id1 = abs(id1in)
1607 id2 = abs(id2in)
1608 
1609  if(id1.ne.id2)then
1610  print*,"Not a valid Z to fermion pair vertex!",id1in,id2in
1611  stop
1612  !Zll
1613  else if(id1.eq.convertlhe(elm_).or.id1.eq.convertlhe(mum_).or.id1.eq.convertlhe(tam_))then
1614  if((id1in*h1).gt.0d0)then
1615  zffbare = ar_lep
1616  else
1617  zffbare = al_lep
1618  endif
1619  !Zuu
1620  else if(id1.eq.convertlhe(up_).or.id1.eq.convertlhe(chm_).or.id1.eq.convertlhe(top_))then
1621  if((id1in*h1).gt.0d0)then
1622  zffbare = ar_qup
1623  else
1624  zffbare = al_qup
1625  endif
1626  !Zdd
1627  else if(id1.eq.convertlhe(dn_).or.id1.eq.convertlhe(str_).or.id1.eq.convertlhe(bot_))then
1628  if((id1in*h1).gt.0d0)then
1629  zffbare = ar_qdn
1630  else
1631  zffbare = al_qdn
1632  endif
1633  !Z nu nu
1634  else if(id1.eq.convertlhe(nue_).or.id1.eq.convertlhe(num_).or.id1.eq.convertlhe(nut_))then
1635  if((id1in*h1).lt.0d0)then
1636  zffbare = al_neu
1637  else
1638  print*,"No right-handed neutrino here!",id1in,h1
1639  stop
1640  endif
1641  else
1642  print*,"Not a valid Z to fermion pair vertex!",id1in,id2in
1643  stop
1644  endif
1645 
1646 return
1647 END FUNCTION
1648 
1649 FUNCTION zff(id1in,id2in,h1,h2)
1650 implicit none
1651 real(8) :: zff
1652 integer :: id1in, id2in
1653 real(8) :: h1,h2
1654  zff = zffbare(id1in,id2in,h1,h2) * sqrt(scalefactor(id1in,id2in))
1655 return
1656 END FUNCTION
1657 
1658 FUNCTION aff(id1in,id2in,h1,h2)
1659 implicit none
1660 real(8) :: aff
1661 integer :: id1, id2, id1in, id2in
1662 real(8) :: h1,h2
1663 logical, save :: printedneutrinowarning = .false.
1664 
1665 id1 = abs(id1in)
1666 id2 = abs(id2in)
1667 
1668  if(id1.ne.id2)then
1669  print*,"Not a valid Z to fermion pair vertex!",id1in,id2in
1670  stop
1671  !photon ll
1672  else if(id1.eq.convertlhe(elm_).or.id1.eq.convertlhe(mum_).or.id1.eq.convertlhe(tam_))then
1673  if((id1in*h1).gt.0d0)then
1674  aff = qlr
1675  else
1676  aff = qll
1677  endif
1678  !photon uu
1679  else if(id1.eq.convertlhe(up_).or.id1.eq.convertlhe(chm_).or.id1.eq.convertlhe(top_))then
1680  if((id1in*h1).gt.0d0)then
1681  aff = qur
1682  else
1683  aff = qul
1684  endif
1685  !photon dd
1686  else if(id1.eq.convertlhe(dn_).or.id1.eq.convertlhe(str_).or.id1.eq.convertlhe(bot_))then
1687  if((id1in*h1).gt.0d0)then
1688  aff = qdr
1689  else
1690  aff = qdl
1691  endif
1692  !photon nu nu
1693  else if(id1.eq.convertlhe(nue_).or.id1.eq.convertlhe(num_).or.id1.eq.convertlhe(nut_))then
1694  aff = 0d0
1695  if (.not. printedneutrinowarning) then
1696  print*,"Warning, gamma > nu nu~ gives 0.",id1in,id2in
1697  printedneutrinowarning = .true.
1698  endif
1699  else
1700  print*,"Not a valid photon to fermion pair vertex!",id1in,id2in
1701  stop
1702  endif
1703 
1704 return
1705 END FUNCTION
1706 
1707 ! for VH use
1708 
1709 
1710 FUNCTION convertlhereverse(Part) ! PDG/PDF->JHU
1711 implicit none
1712 integer :: convertlhereverse
1713 integer :: part
1714 
1715  if( part.eq.0 ) then ! 0=Glu_ is not the official LHE convention
1716  convertlhereverse = glu_
1717  elseif( part.eq.1 ) then
1718  convertlhereverse = dn_
1719  elseif( part.eq.2 ) then
1720  convertlhereverse = up_
1721  elseif( part.eq.3 ) then
1722  convertlhereverse = str_
1723  elseif( part.eq.4 ) then
1724  convertlhereverse = chm_
1725  elseif( part.eq.5 ) then
1726  convertlhereverse = bot_
1727  elseif( part.eq.6 ) then
1728  convertlhereverse = top_
1729  elseif( part.eq.-1 ) then
1730  convertlhereverse = adn_
1731  elseif( part.eq.-2 ) then
1732  convertlhereverse = aup_
1733  elseif( part.eq.-3 ) then
1734  convertlhereverse = astr_
1735  elseif( part.eq.-4 ) then
1736  convertlhereverse = achm_
1737  elseif( part.eq.-5 ) then
1738  convertlhereverse = abot_
1739  elseif( part.eq.-6 ) then
1740  convertlhereverse = atop_
1741  elseif( part.eq.21 ) then
1742  convertlhereverse = glu_
1743  elseif( part.eq.11 ) then
1744  convertlhereverse = elm_
1745  elseif( part.eq.22 ) then
1746  convertlhereverse = pho_
1747  elseif( part.eq.23 ) then
1748  convertlhereverse = z0_
1749  elseif( part.eq.24 ) then
1750  convertlhereverse = wp_
1751  elseif( part.eq.-24 ) then
1752  convertlhereverse = wm_
1753  elseif( part.eq.-11 ) then
1754  convertlhereverse = elp_
1755  elseif( part.eq.13 ) then
1756  convertlhereverse = mum_
1757  elseif( part.eq.-13 ) then
1758  convertlhereverse = mup_
1759  elseif( part.eq.15 ) then
1760  convertlhereverse = tam_
1761  elseif( part.eq.-15 ) then
1762  convertlhereverse = tap_
1763  elseif( part.eq.12 ) then
1764  convertlhereverse = nue_
1765  elseif( part.eq.-12) then
1766  convertlhereverse = anue_
1767  elseif( part.eq.14) then
1768  convertlhereverse = num_
1769  elseif( part.eq.-14) then
1770  convertlhereverse = anum_
1771  elseif( part.eq.16 ) then
1772  convertlhereverse = nut_
1773  elseif( part.eq.-16) then
1774  convertlhereverse = anut_
1775  elseif( part.eq.+25) then
1776  convertlhereverse = hig_
1777  elseif( part.eq.-25) then
1778  convertlhereverse = hig_
1779  elseif( part.eq.32) then
1780  convertlhereverse = zpr_
1781  elseif( part.eq.33) then
1782  convertlhereverse = zpr2_
1783  elseif( part.eq.34) then
1784  convertlhereverse = wppr_
1785  elseif( part.eq.-34) then
1786  convertlhereverse = wmpr_
1787  else
1788  print *, "MYLHE format not implemented for ",part
1789  stop
1790  endif
1791 
1792 
1793 END FUNCTION
1794 
1795 
1796 
1797 FUNCTION convertlhe(Part) ! JHU->PDG/PDF
1798 implicit none
1799 integer :: convertlhe
1800 integer :: part
1801 
1802 
1803  if( part.eq.0 ) then
1804  convertlhe = 0
1805  elseif( part.eq.glu_ ) then
1806  convertlhe = 21
1807  elseif( part.eq.elm_ ) then
1808  convertlhe = 11
1809  elseif( part.eq.elp_ ) then
1810  convertlhe =-11
1811  elseif( part.eq.mum_ ) then
1812  convertlhe = 13
1813  elseif( part.eq.mup_ ) then
1814  convertlhe =-13
1815  elseif( part.eq.tam_ ) then
1816  convertlhe = 15
1817  elseif( part.eq.tap_ ) then
1818  convertlhe =-15
1819  elseif( part.eq.nue_ ) then
1820  convertlhe = 12
1821  elseif( part.eq.anue_) then
1822  convertlhe =-12
1823  elseif( part.eq.num_) then
1824  convertlhe = 14
1825  elseif( part.eq.anum_) then
1826  convertlhe =-14
1827  elseif( part.eq.nut_ ) then
1828  convertlhe = 16
1829  elseif( part.eq.anut_) then
1830  convertlhe =-16
1831  elseif( part.eq.up_ ) then
1832  convertlhe = 2
1833  elseif( part.eq.aup_ ) then
1834  convertlhe =-2
1835  elseif( part.eq.dn_ ) then
1836  convertlhe = 1
1837  elseif( part.eq.adn_ ) then
1838  convertlhe =-1
1839  elseif( part.eq.chm_ ) then
1840  convertlhe = 4
1841  elseif( part.eq.achm_) then
1842  convertlhe =-4
1843  elseif( part.eq.str_ ) then
1844  convertlhe = 3
1845  elseif( part.eq.astr_) then
1846  convertlhe =-3
1847  elseif( part.eq.bot_ ) then
1848  convertlhe = 5
1849  elseif( part.eq.abot_) then
1850  convertlhe =-5
1851  elseif( part.eq.top_ ) then
1852  convertlhe = 6
1853  elseif( part.eq.atop_) then
1854  convertlhe =-6
1855  elseif( part.eq.z0_) then
1856  convertlhe =23
1857  elseif( part.eq.wp_) then
1858  convertlhe =24
1859  elseif( part.eq.wm_) then
1860  convertlhe =-24
1861  elseif( part.eq.pho_) then
1862  convertlhe =22
1863  elseif( part.eq.hig_) then
1864  convertlhe =25
1865  elseif( part.eq.zpr_) then
1866  convertlhe =32
1867  elseif( part.eq.zpr2_) then
1868  convertlhe =33
1869  elseif( part.eq.wppr_) then
1870  convertlhe =34
1871  elseif( part.eq.wmpr_) then
1872  convertlhe =-34
1873  elseif( part.eq.gra_) then
1874  convertlhe =39
1875  elseif( part.le.not_a_particle_) then
1876  convertlhe = not_a_particle_
1877  else
1878  print *, "LHE format not implemented for ",part
1879  stop
1880  endif
1881 
1882 END FUNCTION
1883 
1884 
1885 FUNCTION converttopartindex(Part) ! JHU->PDF
1886 implicit none
1887 integer :: converttopartindex
1888 integer :: part
1889 
1890 
1891  if( part.eq.glu_ ) then
1892  converttopartindex = pdfglu_
1893  elseif( part.eq.up_ ) then
1894  converttopartindex = pdfup_
1895  elseif( part.eq.aup_ ) then
1896  converttopartindex = pdfaup_
1897  elseif( part.eq.dn_ ) then
1898  converttopartindex = pdfdn_
1899  elseif( part.eq.adn_ ) then
1900  converttopartindex = pdfadn_
1901  elseif( part.eq.chm_ ) then
1902  converttopartindex = pdfchm_
1903  elseif( part.eq.achm_) then
1904  converttopartindex = pdfachm_
1905  elseif( part.eq.str_ ) then
1906  converttopartindex = pdfstr_
1907  elseif( part.eq.astr_) then
1908  converttopartindex = pdfastr_
1909  elseif( part.eq.bot_ ) then
1910  converttopartindex = pdfbot_
1911  elseif( part.eq.abot_) then
1912  converttopartindex = pdfabot_
1913  else
1914  print *, "Unsuccessful conversion to a parton ME array index from ",part
1915  stop
1916  endif
1917 
1918 END FUNCTION
1919 
1920 
1921 FUNCTION convertfrompartindex(Part) ! PDF->JHU
1922 implicit none
1923 integer :: convertfrompartindex
1924 integer :: part
1925 
1926 
1927  if( part.eq.pdfglu_ ) then
1928  convertfrompartindex = glu_
1929  elseif( part.eq.pdfup_ ) then
1930  convertfrompartindex = up_
1931  elseif( part.eq.pdfaup_ ) then
1932  convertfrompartindex = aup_
1933  elseif( part.eq.pdfdn_ ) then
1934  convertfrompartindex = dn_
1935  elseif( part.eq.pdfadn_ ) then
1936  convertfrompartindex = adn_
1937  elseif( part.eq.pdfchm_ ) then
1938  convertfrompartindex = chm_
1939  elseif( part.eq.pdfachm_) then
1940  convertfrompartindex = achm_
1941  elseif( part.eq.pdfstr_ ) then
1942  convertfrompartindex = str_
1943  elseif( part.eq.pdfastr_) then
1944  convertfrompartindex = astr_
1945  elseif( part.eq.pdfbot_ ) then
1946  convertfrompartindex = bot_
1947  elseif( part.eq.pdfabot_) then
1948  convertfrompartindex = abot_
1949  elseif( part.eq.pdftop_ ) then
1950  convertfrompartindex = top_
1951  elseif( part.eq.pdfatop_) then
1952  convertfrompartindex = atop_
1953  else
1954  print *, "Unsuccessful conversion to a parton id from the ME array index ",part
1955  stop
1956  endif
1957 
1958 END FUNCTION
1959 
1960 
1961 subroutine setmass(mass, ipart)
1962 implicit none
1963 real(8), intent(in) :: mass
1964 integer, intent(in) :: ipart
1965 integer :: Part
1966 
1967  part=abs(ipart)
1968  if( part.eq.abs(elm_) ) then
1969  m_el = mass
1970  elseif( part.eq.abs(mum_) ) then
1971  m_mu = mass
1972  elseif( part.eq.abs(tam_) ) then
1973  m_tau = mass
1974  elseif( part.eq.abs(chm_) ) then
1975  m_charm = mass
1976  elseif( part.eq.abs(bot_) ) then
1977  m_bot = mass
1978  elseif( part.eq.abs(top_) ) then
1979  m_top = mass
1980  elseif( part.eq.abs(z0_) ) then
1981  m_z = mass
1982  elseif( part.eq.abs(zpr_) ) then
1983  m_zprime = mass
1984  elseif( part.eq.abs(wp_) ) then
1985  m_w = mass
1986  elseif( part.eq.abs(wppr_) ) then
1987  m_wprime = mass
1988  elseif( part.eq.abs(hig_) ) then
1989  m_reso = mass
1990  endif
1991 
1992 end subroutine setmass
1993 
1994 
1995 subroutine setdecaywidth(width, ipart)
1996 implicit none
1997 real(8), intent(in) :: width
1998 integer, intent(in) :: ipart
1999 integer :: Part
2000 
2001  part=abs(ipart)
2002  if( part.eq.abs(tam_) ) then
2003  ga_tau = width
2004  elseif( part.eq.abs(top_) ) then
2005  ga_top = width
2006  elseif( part.eq.abs(z0_) ) then
2007  ga_z = width
2008  elseif( part.eq.abs(zpr_) ) then
2009  ga_zprime = width
2010  elseif( part.eq.abs(wp_) ) then
2011  ga_w = width
2012  elseif( part.eq.abs(wppr_) ) then
2013  ga_wprime = width
2014  elseif( part.eq.abs(hig_) ) then
2015  ga_reso = width
2016  endif
2017 
2018 END subroutine setdecaywidth
2019 
2020 FUNCTION getmass(Part)
2021 implicit none
2022 real(8) :: getmass
2023 integer :: part
2024 
2025 
2026  if( part.eq.glu_ ) then
2027  getmass = 0d0
2028  elseif( abs(part).eq.abs(elm_) ) then
2029  getmass = m_el
2030  elseif( abs(part).eq.abs(mum_) ) then
2031  getmass = m_mu
2032  elseif( abs(part).eq.abs(tam_) ) then
2033  getmass = m_tau
2034  elseif( abs(part).eq.abs(nue_) ) then
2035  getmass = 0d0
2036  elseif( abs(part).eq.abs(num_) ) then
2037  getmass = 0d0
2038  elseif( abs(part).eq.abs(nut_) ) then
2039  getmass = 0d0
2040  elseif( abs(part).eq.abs(up_) ) then
2041  getmass = 0d0
2042  elseif( abs(part).eq.abs(dn_) ) then
2043  getmass = 0d0
2044  elseif( abs(part).eq.abs(chm_) ) then
2045  getmass = m_charm
2046  elseif( abs(part).eq.abs(str_) ) then
2047  getmass = 0d0
2048  elseif( abs(part).eq.abs(bot_) ) then
2049  getmass = m_bot
2050  elseif( abs(part).eq.abs(top_) ) then
2051  getmass = m_top
2052  elseif( abs(part).eq.abs(z0_) ) then
2053  getmass = m_z
2054  elseif( abs(part).eq.abs(zpr_) ) then
2055  getmass = m_zprime
2056  elseif( abs(part).eq.abs(wp_) ) then
2057  getmass = m_w
2058  elseif( abs(part).eq.abs(wppr_) ) then
2059  getmass = m_wprime
2060  elseif( abs(part).eq.abs(pho_) ) then
2061  getmass = 0d0
2062  elseif( abs(part).eq.abs(hig_) ) then
2063  getmass = m_reso
2064  elseif( part.eq.not_a_particle_) then
2065  getmass = 0d0
2066  else
2067  print *, "Error in getMass",part
2068  stop
2069  endif
2070 
2071 END FUNCTION
2072 
2073 FUNCTION getdecaywidth(Part)
2074 implicit none
2075 real(8) :: getdecaywidth
2076 integer :: part
2077 
2078  getdecaywidth = 0d0
2079  if( abs(part).eq.abs(top_) ) then
2080  getdecaywidth = ga_top
2081  elseif( abs(part).eq.abs(z0_) ) then
2082  getdecaywidth = ga_z
2083  elseif( abs(part).eq.abs(zpr_) ) then
2084  getdecaywidth = ga_zprime
2085  elseif( abs(part).eq.abs(wp_) ) then
2086  getdecaywidth = ga_w
2087  elseif( abs(part).eq.abs(wppr_) ) then
2088  getdecaywidth = ga_wprime
2089  elseif( abs(part).eq.abs(hig_) ) then
2090  getdecaywidth = ga_reso
2091  elseif( abs(part).eq.abs(tam_) ) then
2092  getdecaywidth = ga_tau
2093  endif
2094 
2095 END FUNCTION
2096 
2097 
2098 subroutine sethiggsmass(jH,mass)
2099 implicit none
2100 real(8) :: mass
2101 integer :: jH
2102  if (jh .eq. 1) then
2103  m_reso=mass
2104  elseif (jh .eq.2) then
2105  m_reso2=mass
2106  endif
2107 END subroutine
2108 
2109 subroutine sethiggsdecaywidth(jH,width)
2110 implicit none
2111 real(8) :: width
2112 integer :: jH
2113  if (jh .eq. 1) then
2114  ga_reso=width
2115  elseif (jh .eq.2) then
2116  ga_reso2=width
2117  endif
2118 END subroutine
2119 
2120 FUNCTION gethiggsmass(jH)
2121 implicit none
2122 real(8) :: gethiggsmass
2123 integer :: jh
2124  gethiggsmass=0d0
2125  if (jh .eq. 1) then
2126  gethiggsmass = m_reso
2127  elseif (jh .eq.2) then
2128  gethiggsmass = m_reso2
2129  endif
2130 END FUNCTION
2131 
2132 FUNCTION gethiggsdecaywidth(jH)
2133 implicit none
2134 real(8) :: gethiggsdecaywidth
2135 integer :: jh
2136  gethiggsdecaywidth=0d0
2137  if (jh .eq. 1) then
2138  gethiggsdecaywidth = ga_reso
2139  elseif (jh .eq.2) then
2140  gethiggsdecaywidth = ga_reso2
2141  endif
2142 END FUNCTION
2143 
2144 
2145 FUNCTION getparticle(Part)
2146 implicit none
2147 character(len=3) :: getparticle
2148 integer :: part
2149 
2150 
2151  if( part.eq.glu_ ) then
2152  getparticle = "glu"
2153  elseif( part.eq.0 ) then
2154  getparticle = "glu"
2155  elseif( part.eq.elm_ ) then
2156  getparticle = "el-"
2157  elseif( part.eq.elp_ ) then
2158  getparticle = "el+"
2159  elseif( part.eq.mum_ ) then
2160  getparticle = "mu-"
2161  elseif( part.eq.mup_ ) then
2162  getparticle = "mu+"
2163  elseif( part.eq.tam_ ) then
2164  getparticle = "ta-"
2165  elseif( part.eq.tap_ ) then
2166  getparticle = "t+-"
2167  elseif( part.eq.nue_ ) then
2168  getparticle = "nuE"
2169  elseif( part.eq.anue_ ) then
2170  getparticle = "AnE"
2171  elseif( part.eq.num_ ) then
2172  getparticle = "nuM"
2173  elseif( part.eq.anum_ ) then
2174  getparticle = "AnM"
2175  elseif( part.eq.nut_ ) then
2176  getparticle = "nuT"
2177  elseif( part.eq.anut_ ) then
2178  getparticle = "AnT"
2179  elseif( part.eq.up_ ) then
2180  getparticle = " up"
2181  elseif( part.eq.aup_ ) then
2182  getparticle = "Aup"
2183  elseif( part.eq.dn_ ) then
2184  getparticle = " dn"
2185  elseif( part.eq.adn_ ) then
2186  getparticle = "Adn"
2187  elseif( part.eq.chm_ ) then
2188  getparticle = "chm"
2189  elseif( part.eq.achm_ ) then
2190  getparticle = "Achm"
2191  elseif( part.eq.str_ ) then
2192  getparticle = "str"
2193  elseif( part.eq.astr_ ) then
2194  getparticle = "Astr"
2195  elseif( part.eq.bot_ ) then
2196  getparticle = "bot"
2197  elseif( part.eq.abot_ ) then
2198  getparticle = "Abot"
2199  elseif( part.eq.top_ ) then
2200  getparticle = "top"
2201  elseif( part.eq.atop_ ) then
2202  getparticle = "Atop"
2203  elseif( part.eq.z0_ ) then
2204  getparticle = " Z0"
2205  elseif( part.eq.zpr_ ) then
2206  getparticle = " Zprime"
2207  elseif( part.eq.wp_ ) then
2208  getparticle = " W+"
2209  elseif( part.eq.wm_ ) then
2210  getparticle = " W-"
2211  elseif( part.eq.wppr_ ) then
2212  getparticle = " W+prime"
2213  elseif( part.eq.wmpr_ ) then
2214  getparticle = " W-prime"
2215  elseif( part.eq.pho_ ) then
2216  getparticle = "pho"
2217  elseif( part.eq.hig_ ) then
2218  getparticle = "Hig"
2219  else
2220  print *, "Error in getParticle",part
2221  stop
2222  endif
2223 
2224 
2225 END FUNCTION
2226 
2227 
2228 FUNCTION getlheparticle(PartLHE)
2229 implicit none
2230 character(len=3) :: getlheparticle
2231 integer :: partlhe,part
2232 
2233 
2234  part = convertlhereverse(partlhe)
2235  getlheparticle = getparticle(part)
2236 
2237 END FUNCTION
2238 
2239 
2240 
2241 
2242 FUNCTION isazdecay(DKMode)
2243 implicit none
2244 logical :: isazdecay
2245 integer :: dkmode
2246 
2247 
2248  if( dkmode.eq.0 ) then
2249  isazdecay = .true.
2250  elseif( dkmode.eq.1 ) then
2251  isazdecay = .true.
2252  elseif( dkmode.eq.2 ) then
2253  isazdecay = .true.
2254  elseif( dkmode.eq.3 ) then
2255  isazdecay = .true.
2256  elseif( dkmode.eq.8 ) then
2257  isazdecay = .true.
2258  elseif( dkmode.eq.9 ) then
2259  isazdecay = .true.
2260  elseif( dkmode.lt.0 ) then
2261  isazdecay = ( &
2262  dkmode.eq.-2*2 .or. & ! Z->ee
2263  dkmode.eq.-3*3 .or. & ! Z->mumu
2264  dkmode.eq.-5*5 .or. & ! Z->dd
2265  dkmode.eq.-7*7 .or. & ! Z->uu
2266  dkmode.eq.-11*11 .or. & ! Z->ss
2267  dkmode.eq.-13*13 .or. & ! Z->cc
2268  dkmode.eq.-17*17 & ! Z->bb
2269  )
2270  else
2271  isazdecay=.false.
2272  endif
2273 
2274 END FUNCTION
2275 
2276 
2277 FUNCTION isawdecay(DKMode)
2278 implicit none
2279 logical :: isawdecay
2280 integer :: dkmode
2281 
2282 
2283  if( dkmode.eq.4 ) then
2284  isawdecay = .true.
2285  elseif( dkmode.eq.5 ) then
2286  isawdecay = .true.
2287  elseif( dkmode.eq.6 ) then
2288  isawdecay = .true.
2289  elseif( dkmode.eq.10 ) then
2290  isawdecay = .true.
2291  elseif( dkmode.eq.11 ) then
2292  isawdecay = .true.
2293  elseif( dkmode.lt.0 ) then
2294  isawdecay = ( &
2295  dkmode.eq.-2*1 .or. & ! W->enu
2296  dkmode.eq.-3*1 .or. & ! W->munu
2297  dkmode.eq.-5*7 .or. & ! W->du
2298  dkmode.eq.-5*13 .or. & ! W->dc
2299  dkmode.eq.-11*7 .or. & ! W->su
2300  dkmode.eq.-11*13 .or. & ! W->sc
2301  dkmode.eq.-17*7 .or. & ! W->bu
2302  dkmode.eq.-17*13 & ! W->bc
2303  )
2304  else
2305  isawdecay=.false.
2306  endif
2307 
2308 
2309 END FUNCTION
2310 
2311 
2312 FUNCTION isaphoton(DKMode)
2313 implicit none
2314 logical :: isaphoton
2315 integer :: dkmode
2316 
2317 
2318  if( dkmode.eq.7 ) then
2319  isaphoton = .true.
2320  else
2321  isaphoton=.false.
2322  endif
2323 
2324 
2325 END FUNCTION
2326 
2327 
2328 
2329 FUNCTION isdowntypequark(PartType)
2330 implicit none
2331 logical :: isdowntypequark
2332 integer :: parttype
2333  isdowntypequark = ( abs(parttype).eq.abs(dn_) .or. abs(parttype).eq.abs(str_) .or. abs(parttype).eq.abs(bot_) )
2334 END FUNCTION
2335 FUNCTION islhedowntypequark(PartType)
2336 implicit none
2337 logical :: islhedowntypequark
2338 integer :: parttype
2339  islhedowntypequark = ( abs(parttype).eq.1 .or. abs(parttype).eq.3 .or. abs(parttype).eq.5 )
2340 END FUNCTION
2341 
2342 FUNCTION isuptypequark(PartType)
2343 implicit none
2344 logical :: isuptypequark
2345 integer :: parttype
2346  isuptypequark = ( abs(parttype).eq.abs(up_) .or. abs(parttype).eq.abs(chm_) .or. abs(parttype).eq.abs(top_) )
2347 END FUNCTION
2348 FUNCTION islheuptypequark(PartType)
2349 implicit none
2350 logical :: islheuptypequark
2351 integer :: parttype
2352  islheuptypequark = ( abs(parttype).eq.2 .or. abs(parttype).eq.4 .or. abs(parttype).eq.6 )
2353 END FUNCTION
2354 
2355 FUNCTION isuptypelightquark(PartType)
2356 implicit none
2357 logical :: isuptypelightquark
2358 integer :: parttype
2359  isuptypelightquark = ( abs(parttype).eq.abs(up_) .or. abs(parttype).eq.abs(chm_))
2360 END FUNCTION
2361 FUNCTION islheuptypelightquark(PartType)
2362 implicit none
2363 logical :: islheuptypelightquark
2364 integer :: parttype
2365  islheuptypelightquark = ( abs(parttype).eq.2 .or. abs(parttype).eq.4 )
2366 END FUNCTION
2367 
2368 FUNCTION isaquark(PartType)
2369 implicit none
2370 logical :: isaquark
2371 integer :: parttype
2372  isaquark=isuptypequark(parttype) .or. isdowntypequark(parttype)
2373 END FUNCTION
2374 FUNCTION isalhequark(PartType)
2375 implicit none
2376 logical :: isalhequark
2377 integer :: parttype
2378  isalhequark=islheuptypequark(parttype) .or. islhedowntypequark(parttype)
2379 END FUNCTION
2380 
2381 FUNCTION isalightquark(PartType)
2382 implicit none
2383 logical :: isalightquark
2384 integer :: parttype
2385  isalightquark=isuptypelightquark(parttype) .or. isdowntypequark(parttype)
2386 END FUNCTION
2387 FUNCTION isalhelightquark(PartType)
2388 implicit none
2389 logical :: isalhelightquark
2390 integer :: parttype
2391  isalhelightquark=islheuptypelightquark(parttype) .or. islhedowntypequark(parttype)
2392 END FUNCTION
2393 
2394 
2395 FUNCTION isaneutrino(PartType)
2396 implicit none
2397 logical :: isaneutrino
2398 integer :: parttype
2399  isaneutrino = ( abs(parttype).eq.abs(nue_) .or. abs(parttype).eq.abs(num_) .or. abs(parttype).eq.abs(nut_) )
2400 END FUNCTION
2401 FUNCTION isalheneutrino(PartType)
2402 implicit none
2403 logical :: isalheneutrino
2404 integer :: parttype
2405  isalheneutrino = ( abs(parttype).eq.12 .or. abs(parttype).eq.14 .or. abs(parttype).eq.16 )
2406 END FUNCTION
2407 
2408 FUNCTION isalepton(PartType)! note that lepton means charged lepton here
2409 implicit none
2410 logical :: isalepton
2411 integer :: parttype
2412  isalepton = ( abs(parttype).eq.abs(elp_) .or. abs(parttype).eq.abs(mup_) .or. abs(parttype).eq.abs(tap_) )
2413 END FUNCTION
2414 FUNCTION isalhelepton(PartType)! note that lepton means charged lepton here
2415 implicit none
2416 logical :: isalhelepton
2417 integer :: parttype
2418  isalhelepton = ( abs(parttype).eq.11 .or. abs(parttype).eq.13 .or. abs(parttype).eq.15 )
2419 END FUNCTION
2420 
2421 
2422 FUNCTION isagluon(PartType)
2423 implicit none
2424 logical :: isagluon
2425 integer :: parttype
2426  isagluon = (abs(parttype).eq.glu_)
2427 END FUNCTION
2428 
2429 FUNCTION isalhegluon(PartType)
2430 implicit none
2431 logical :: isalhegluon
2432 integer :: parttype
2433  isalhegluon = (abs(parttype).eq.convertlhe(glu_))
2434 END FUNCTION
2435 
2436 
2437 FUNCTION isajet(PartType)
2438 implicit none
2439 logical :: isajet
2440 integer :: parttype
2441  isajet = (isalightquark(parttype) .or. isagluon(parttype))
2442 END FUNCTION
2443 
2444 FUNCTION isalhejet(PartType)
2445 implicit none
2446 logical :: isalhejet
2447 integer :: parttype
2448  isalhejet = (isalhelightquark(parttype) .or. isalhegluon(parttype))
2449 END FUNCTION
2450 
2451 
2452 FUNCTION isaboson(PartType)
2453 implicit none
2454 logical :: isaboson
2455 integer :: parttype
2456  isaboson = ( abs(parttype).eq.abs(pho_) .or. abs(parttype).eq.abs(z0_) .or. abs(parttype).eq.abs(wp_) .or. abs(parttype).eq.abs(hig_) )
2457 END FUNCTION
2458 
2459 FUNCTION isalheboson(PartType)
2460 implicit none
2461 logical :: isalheboson
2462 integer :: parttype
2463  isalheboson = ( abs(parttype).ge.22 .and. abs(parttype).le.25 )
2464 END FUNCTION
2465 
2466 function coupledvertex(id,hel,useAHcoupl)
2467  implicit none
2468  integer, optional :: useahcoupl
2469  integer, intent(in) :: id(1:2),hel
2470  integer :: testahcoupl
2471  integer :: coupledvertex
2472 
2473  testahcoupl = 0
2474  if(present(useahcoupl)) then
2475  testahcoupl = useahcoupl
2476  endif
2477  if( (&
2478  (id(1).eq.elp_ .and. id(2).eq.nue_) .or. (id(2).eq.elp_ .and. id(1).eq.nue_) .or. &
2479  (id(1).eq.mup_ .and. id(2).eq.num_) .or. (id(2).eq.mup_ .and. id(1).eq.num_) .or. &
2480  (id(1).eq.tap_ .and. id(2).eq.nut_) .or. (id(2).eq.tap_ .and. id(1).eq.nut_) .or. &
2481  (id(1).eq.up_ .and. (id(2).eq.adn_ .or. id(2).eq.astr_ .or. id(2).eq.abot_)) .or. (id(2).eq.up_ .and. (id(1).eq.adn_ .or. id(1).eq.astr_ .or. id(1).eq.abot_)) .or. &
2482  (id(1).eq.chm_ .and. (id(2).eq.adn_ .or. id(2).eq.astr_ .or. id(2).eq.abot_)) .or. (id(2).eq.chm_ .and. (id(1).eq.adn_ .or. id(1).eq.astr_ .or. id(1).eq.abot_)) .or. &
2483  (id(1).eq.top_ .and. (id(2).eq.adn_ .or. id(2).eq.astr_ .or. id(2).eq.abot_)) .or. (id(2).eq.top_ .and. (id(1).eq.adn_ .or. id(1).eq.astr_ .or. id(1).eq.abot_)) &
2484  ) .and. hel.lt.0) then
2485  coupledvertex=wp_
2486  elseif( (&
2487  (id(1).eq.elm_ .and. id(2).eq.anue_) .or. (id(2).eq.elm_ .and. id(1).eq.anue_) .or. &
2488  (id(1).eq.mum_ .and. id(2).eq.anum_) .or. (id(2).eq.mum_ .and. id(1).eq.anum_) .or. &
2489  (id(1).eq.tam_ .and. id(2).eq.anut_) .or. (id(2).eq.tam_ .and. id(1).eq.anut_) .or. &
2490  (id(1).eq.aup_ .and. (id(2).eq.dn_ .or. id(2).eq.str_ .or. id(2).eq.bot_)) .or. (id(2).eq.aup_ .and. (id(1).eq.dn_ .or. id(1).eq.str_ .or. id(1).eq.bot_)) .or. &
2491  (id(1).eq.achm_ .and. (id(2).eq.dn_ .or. id(2).eq.str_ .or. id(2).eq.bot_)) .or. (id(2).eq.achm_ .and. (id(1).eq.dn_ .or. id(1).eq.str_ .or. id(1).eq.bot_)) .or. &
2492  (id(1).eq.atop_ .and. (id(2).eq.dn_ .or. id(2).eq.str_ .or. id(2).eq.bot_)) .or. (id(2).eq.atop_ .and. (id(1).eq.dn_ .or. id(1).eq.str_ .or. id(1).eq.bot_)) &
2493  ) .and. hel.lt.0) then
2494  coupledvertex=wm_
2495  elseif( (&
2496  (id(1).eq.elm_ .and. id(2).eq.elp_) .or. (id(2).eq.elm_ .and. id(1).eq.elp_) .or. &
2497  (id(1).eq.mum_ .and. id(2).eq.mup_) .or. (id(2).eq.mum_ .and. id(1).eq.mup_) .or. &
2498  (id(1).eq.tam_ .and. id(2).eq.tap_) .or. (id(2).eq.tam_ .and. id(1).eq.tap_) .or. &
2499  (id(1).eq.up_ .and. id(2).eq.aup_) .or. (id(2).eq.up_ .and. id(1).eq.aup_) .or. &
2500  (id(1).eq.dn_ .and. id(2).eq.adn_) .or. (id(2).eq.dn_ .and. id(1).eq.adn_) .or. &
2501  (id(1).eq.chm_ .and. id(2).eq.achm_) .or. (id(2).eq.chm_ .and. id(1).eq.achm_) .or. &
2502  (id(1).eq.str_ .and. id(2).eq.astr_) .or. (id(2).eq.str_ .and. id(1).eq.astr_) .or. &
2503  (id(1).eq.top_ .and. id(2).eq.atop_) .or. (id(2).eq.top_ .and. id(1).eq.atop_) .or. &
2504  (id(1).eq.bot_ .and. id(2).eq.abot_) .or. (id(2).eq.bot_ .and. id(1).eq.abot_) &
2505  ) .and. hel.ne.0) then
2506  if(testahcoupl.eq.1) then
2507  coupledvertex=pho_
2508  elseif(testahcoupl.eq.2) then
2509  coupledvertex=hig_
2510  else
2511  coupledvertex=z0_
2512  endif
2513  elseif( (&
2514  (id(1).eq.nue_ .and. id(2).eq.anue_) .or. (id(2).eq.nue_ .and. id(1).eq.anue_) .or. &
2515  (id(1).eq.num_ .and. id(2).eq.anum_) .or. (id(2).eq.num_ .and. id(1).eq.anum_) .or. &
2516  (id(1).eq.nut_ .and. id(2).eq.anut_) .or. (id(2).eq.nut_ .and. id(1).eq.anut_) &
2517  ) .and. hel.lt.0) then
2518  coupledvertex=z0_ ! Only Z coupling to nuL-nubR
2519  else
2520  coupledvertex=not_a_particle_
2521  endif
2522 
2523  return
2524 end function coupledvertex
2525 
2526 function coupledvertexisdiagonal(id,hel)
2527  implicit none
2528  integer, intent(in) :: id(1:2),hel
2529  logical :: coupledvertexisdiagonal
2530 
2531  if( (&
2532  (id(1).eq.elp_ .and. id(2).eq.nue_) .or. (id(2).eq.elp_ .and. id(1).eq.nue_) .or. & ! W+ -> l+ nu
2533  (id(1).eq.mup_ .and. id(2).eq.num_) .or. (id(2).eq.mup_ .and. id(1).eq.num_) .or. &
2534  (id(1).eq.tap_ .and. id(2).eq.nut_) .or. (id(2).eq.tap_ .and. id(1).eq.nut_) .or. &
2535  (id(1).eq.up_ .and. id(2).eq.adn_) .or. (id(2).eq.up_ .and. id(1).eq.adn_) .or. & ! W+ -> qu qbd
2536  (id(1).eq.chm_ .and. id(2).eq.astr_) .or. (id(2).eq.chm_ .and. id(1).eq.astr_) .or. &
2537  (id(1).eq.top_ .and. id(2).eq.abot_) .or. (id(2).eq.top_ .and. id(1).eq.abot_) .or. &
2538  (id(1).eq.elm_ .and. id(2).eq.anue_) .or. (id(2).eq.elm_ .and. id(1).eq.anue_) .or. & ! W- -> l- nub
2539  (id(1).eq.mum_ .and. id(2).eq.anum_) .or. (id(2).eq.mum_ .and. id(1).eq.anum_) .or. &
2540  (id(1).eq.tam_ .and. id(2).eq.anut_) .or. (id(2).eq.tam_ .and. id(1).eq.anut_) .or. &
2541  (id(1).eq.aup_ .and. id(2).eq.dn_ ) .or. (id(2).eq.aup_ .and. id(1).eq.dn_) .or. & ! W- -> qd qbu
2542  (id(1).eq.achm_ .and. id(2).eq.str_) .or. (id(2).eq.achm_ .and. id(1).eq.str_) .or. &
2543  (id(1).eq.atop_ .and. id(2).eq.bot_) .or. (id(2).eq.atop_ .and. id(1).eq.bot_) .or. &
2544  (id(1).eq.nue_ .and. id(2).eq.anue_) .or. (id(2).eq.nue_ .and. id(1).eq.anue_) .or. & ! Z -> nu nu
2545  (id(1).eq.num_ .and. id(2).eq.anum_) .or. (id(2).eq.num_ .and. id(1).eq.anum_) .or. &
2546  (id(1).eq.nut_ .and. id(2).eq.anut_) .or. (id(2).eq.nut_ .and. id(1).eq.anut_) &
2547  ) .and. hel.lt.0) then
2548  coupledvertexisdiagonal=.true.
2549  elseif( (&
2550  (id(1).eq.elm_ .and. id(2).eq.elp_) .or. (id(2).eq.elm_ .and. id(1).eq.elp_) .or. & ! Z -> ll
2551  (id(1).eq.mum_ .and. id(2).eq.mup_) .or. (id(2).eq.mum_ .and. id(1).eq.mup_) .or. &
2552  (id(1).eq.tam_ .and. id(2).eq.tap_) .or. (id(2).eq.tam_ .and. id(1).eq.tap_) .or. &
2553  (id(1).eq.up_ .and. id(2).eq.aup_) .or. (id(2).eq.up_ .and. id(1).eq.aup_) .or. & ! Z -> qq
2554  (id(1).eq.dn_ .and. id(2).eq.adn_) .or. (id(2).eq.dn_ .and. id(1).eq.adn_) .or. &
2555  (id(1).eq.chm_ .and. id(2).eq.achm_) .or. (id(2).eq.chm_ .and. id(1).eq.achm_) .or. &
2556  (id(1).eq.str_ .and. id(2).eq.astr_) .or. (id(2).eq.str_ .and. id(1).eq.astr_) .or. &
2557  (id(1).eq.top_ .and. id(2).eq.atop_) .or. (id(2).eq.top_ .and. id(1).eq.atop_) .or. &
2558  (id(1).eq.bot_ .and. id(2).eq.abot_) .or. (id(2).eq.bot_ .and. id(1).eq.abot_) &
2559  ) .and. hel.ne.0) then
2560  coupledvertexisdiagonal=.true.
2561  else
2562  coupledvertexisdiagonal=.false.
2563  endif
2564 
2565  return
2566 end function coupledvertexisdiagonal
2567 
2568 
2569 
2570 ! Counts *charged* leptons to be specific
2571 FUNCTION countleptons( MY_IDUP )
2572 implicit none
2573 integer :: my_idup(:),countleptons
2574 integer :: i
2575 
2576  countleptons = 0
2577  do i = 1,size(my_idup)
2578  if( isalepton( my_idup(i) ) ) countleptons=countleptons+1
2579  enddo
2580 
2581 
2582 RETURN
2583 END FUNCTION
2584 
2585 ! Counts light jets
2586 FUNCTION countjets( MY_IDUP )
2587 implicit none
2588 integer :: my_idup(:),countjets
2589 integer :: i
2590 
2591  countjets = 0
2592  do i = 1,size(my_idup)
2593  if( isajet( my_idup(i) ) ) countjets=countjets+1
2594  enddo
2595 
2596 
2597 RETURN
2598 END FUNCTION
2599 
2600 
2601 
2602 FUNCTION su2flip(Part)
2603 implicit none
2604 integer :: su2flip
2605 integer :: part
2606 
2607  if( abs(part).eq.up_ ) then
2608  su2flip = sign(1,part)*dn_
2609  elseif( abs(part).eq.dn_ ) then
2610  su2flip = sign(1,part)*up_
2611  elseif( abs(part).eq.chm_ ) then
2612  su2flip = sign(1,part)*str_
2613  elseif( abs(part).eq.str_ ) then
2614  su2flip = sign(1,part)*chm_
2615  elseif( abs(part).eq.bot_ ) then
2616  su2flip = sign(1,part)*top_
2617  elseif( abs(part).eq.top_ ) then
2618  su2flip = sign(1,part)*bot_
2619  elseif( abs(part).eq.elp_ ) then
2620  su2flip = sign(1,part)*nue_
2621  elseif( abs(part).eq.mup_ ) then
2622  su2flip = sign(1,part)*num_
2623  elseif( abs(part).eq.tap_ ) then
2624  su2flip = sign(1,part)*nut_
2625  elseif( abs(part).eq.nue_ ) then
2626  su2flip = sign(1,part)*elp_
2627  elseif( abs(part).eq.num_ ) then
2628  su2flip = sign(1,part)*mup_
2629  elseif( abs(part).eq.nut_ ) then
2630  su2flip = sign(1,part)*tap_
2631  else
2632  print *, "Error: Invalid flavor in SU2flip ",part
2633  stop
2634  endif
2635 
2636 END FUNCTION
2637 
2638 
2639 Function coupltolhewp(Part)
2640  implicit none
2641  logical :: coupltolhewp
2642  integer, intent(in) :: part(1:2)
2643 
2644  coupltolhewp = &
2645  ((part(1).eq.2 .and.part(2).eq.-1) &
2646  .or.(part(1).eq.2 .and.part(2).eq.-3) &
2647  .or.(part(1).eq.2 .and.part(2).eq.-5) &
2648  .or.(part(1).eq.4 .and.part(2).eq.-1) &
2649  .or.(part(1).eq.4 .and.part(2).eq.-3) &
2650  .or.(part(1).eq.4 .and.part(2).eq.-5) &
2651  .or.(part(1).eq.-1.and.part(2).eq. 2) &
2652  .or.(part(1).eq.-3.and.part(2).eq. 2) &
2653  .or.(part(1).eq.-5.and.part(2).eq. 2) &
2654  .or.(part(1).eq.-1.and.part(2).eq. 4) &
2655  .or.(part(1).eq.-3.and.part(2).eq. 4) &
2656  .or.(part(1).eq.-5.and.part(2).eq. 4))
2657 
2658  return
2659 end function coupltolhewp
2660 
2661 
2662 Function coupltolhewm(Part)
2663  implicit none
2664  logical :: coupltolhewm
2665  integer, intent(in) :: part(1:2)
2666 
2667  coupltolhewm = &
2668  ((part(1).eq.-2.and.part(2).eq. 1) &
2669  .or.(part(1).eq.-2.and.part(2).eq. 3) &
2670  .or.(part(1).eq.-2.and.part(2).eq. 5) &
2671  .or.(part(1).eq.-4.and.part(2).eq. 1) &
2672  .or.(part(1).eq.-4.and.part(2).eq. 3) &
2673  .or.(part(1).eq.-4.and.part(2).eq. 5) &
2674  .or.(part(1).eq. 1.and.part(2).eq.-2) &
2675  .or.(part(1).eq. 3.and.part(2).eq.-2) &
2676  .or.(part(1).eq. 5.and.part(2).eq.-2) &
2677  .or.(part(1).eq. 1.and.part(2).eq.-4) &
2678  .or.(part(1).eq. 3.and.part(2).eq.-4) &
2679  .or.(part(1).eq. 5.and.part(2).eq.-4))
2680 
2681  return
2682 end function coupltolhewm
2683 
2684 
2685 
2686 
2687 
2688 FUNCTION chargeflip(Part)
2689 implicit none
2690 integer :: chargeflip
2691 integer :: part
2692 
2693  if( (abs(part).ge.1 .and. abs(part).le.9) .or. (abs(part).ge.13 .and. abs(part).le.16)) then! quarks, leptons, W's, neutrinos
2694  chargeflip = -part
2695  elseif( (abs(part).ge.10 .and. abs(part).le.12) ) then! glu,pho,Z0
2696  chargeflip = +part
2697  elseif( part.eq.25 .or. part.eq.32 .or. part.eq.39 ) then ! Higgs,Zprime,Graviton
2698  chargeflip = +part
2699  else
2700  print *, "Error: Invalid flavor in ChargeFlip"
2701  stop
2702  endif
2703 
2704 END FUNCTION
2705 
2706 
2707 subroutine computeckmelements(inVCKM_ud, inVCKM_us, inVCKM_cd, inVCKM_cs, inVCKM_ts, inVCKM_tb, inVCKM_ub, inVCKM_cb, inVCKM_td)
2708 implicit none
2709 real(8) :: inVCKM_ud
2710 real(8) :: inVCKM_us
2711 real(8) :: inVCKM_cd
2712 real(8) :: inVCKM_cs
2713 real(8) :: inVCKM_ts
2714 real(8) :: inVCKM_tb
2715 real(8), optional :: inVCKM_ub
2716 real(8), optional :: inVCKM_cb
2717 real(8), optional :: inVCKM_td
2718 real(8) :: sumVsq(1:3),diffVsq
2719 
2720  vckm_ud=invckm_ud
2721  vckm_us=invckm_us
2722  vckm_cd=invckm_cd
2723  vckm_cs=invckm_cs
2724  vckm_ts=invckm_ts
2725  vckm_tb=invckm_tb
2726 
2727  if(present(invckm_ub)) then
2728  vckm_ub = invckm_ub
2729  else
2730  diffvsq = 1d0-vckm_ud**2-vckm_us**2
2731  if (diffvsq.ge.0d0) then
2732  vckm_ub = sqrt(diffvsq)
2733  else
2734  vckm_ub = 0d0
2735  endif
2736  endif
2737  if(present(invckm_cb)) then
2738  vckm_cb = invckm_cb
2739  else
2740  diffvsq = 1d0-vckm_cd**2-vckm_cs**2
2741  if (diffvsq.ge.0d0) then
2742  vckm_cb = sqrt(diffvsq)
2743  else
2744  vckm_cb = 0d0
2745  endif
2746  endif
2747  if(present(invckm_td)) then
2748  vckm_td = invckm_td
2749  else
2750  diffvsq = 1d0-vckm_td**2-vckm_ts**2
2751  if (diffvsq.ge.0d0) then
2752  vckm_tb = sqrt(diffvsq)
2753  else
2754  vckm_tb = 0d0
2755  endif
2756  endif
2757 
2758  sumvsq(1) = vckm_ub**2+vckm_us**2+vckm_ud**2
2759  vckm_ud=vckm_ud/sqrt(sumvsq(1))
2760  vckm_us=vckm_us/sqrt(sumvsq(1))
2761  vckm_ub=vckm_ub/sqrt(sumvsq(1))
2762  sumvsq(2) = vckm_cb**2+vckm_cs**2+vckm_cd**2
2763  vckm_cd=vckm_cd/sqrt(sumvsq(2))
2764  vckm_cs=vckm_cs/sqrt(sumvsq(2))
2765  vckm_cb=vckm_cb/sqrt(sumvsq(2))
2766  sumvsq(3) = vckm_tb**2+vckm_ts**2+vckm_td**2
2767  vckm_td=vckm_td/sqrt(sumvsq(3))
2768  vckm_ts=vckm_ts/sqrt(sumvsq(3))
2769  vckm_tb=vckm_tb/sqrt(sumvsq(3))
2770 
2771 end subroutine computeckmelements
2772 
2773 subroutine setdefaultckm()
2774 implicit none
2775  vckm_ud = 0.974285d0
2776  vckm_us = 0.225290d0
2777  vckm_cs = 0.9734244d0
2778  vckm_cd =-0.225182d0
2779  vckm_tb = 0.99912367d0
2780  vckm_ts =-0.040920069d0
2781  call computeckmelements(vckm_ud, vckm_us, vckm_cd, vckm_cs, vckm_ts, vckm_tb)
2782 end subroutine setdefaultckm
2783 
2784 subroutine computeewvariables()
2785 implicit none
2786 
2787  ! Calculate fundamental couplings
2788  vev = 1.0d0/sqrt(gf*sqrt(2.0d0))
2789  gwsq = 4.0d0 * m_w**2/vev**2
2790  sitw = sqrt(xw)
2791  twosc = sqrt(4d0*xw*(1d0-xw))
2792 
2793  esq = 4.0_dp*pi*alpha_qed
2794  overallcouplvffsq = esq/(twosc**2) ! ~= gwsq/4.0_dp/(1.0_dp-xw)
2795  ! Z couplings
2796  ar_lep =2d0*(t3lr-qlr*xw)
2797  al_lep =2d0*(t3ll-qll*xw)
2798  ar_neu =2d0*(t3nr-qnr*xw)
2799  al_neu =2d0*(t3nl-qnl*xw)
2800  ar_qup =2d0*(t3ur-qur*xw)
2801  al_qup =2d0*(t3ul-qul*xw)
2802  ar_qdn =2d0*(t3dr-qdr*xw)
2803  al_qdn =2d0*(t3dl-qdl*xw)
2804 
2805  ! Zff couplings: adding user specified values (default 0.0) to SM values
2806  al_lep = al_lep + daz_lep_left
2807  ar_lep = ar_lep + daz_lep_right
2808  al_neu = al_neu + daz_neu_left
2809  ar_neu = ar_neu + daz_neu_right
2810  al_qup = al_qup + daz_qup_left
2811  ar_qup = ar_qup + daz_qup_right
2812  al_qdn = al_qdn + daz_qdn_left
2813  ar_qdn = ar_qdn + daz_qdn_right
2814 
2815  ! W couplings
2816  bl = sqrt(2d0*(1d0-xw))
2817  br = 0d0
2818  ! A couplings
2819  cr_lep = -twosc*qlr
2820  cl_lep = -twosc*qll
2821  cr_neu = -twosc*qnr
2822  cl_neu = -twosc*qnl
2823  cr_qup = -twosc*qur
2824  cl_qup = -twosc*qul
2825  cr_qdn = -twosc*qdr
2826  cl_qdn = -twosc*qdl
2827 
2828  ! Normalizations used in VH and VBF
2829  couplwffsq = gwsq/2.0_dp
2830  couplzffsq = gwsq/4.0_dp/(1.0_dp-xw)
2831  !couplZffsq = alpha_QED*4d0*pi/xw/4.0_dp/(1.0_dp-xw)
2832  couplazff = -gwsq*sitw/2.0_dp/sqrt(1.0_dp-xw)
2833  couplaffsq = gwsq*xw
2834 
2835 
2836 
2837 
2838 end subroutine computeewvariables
2839 
2840 subroutine computeqcdvariables()
2841 implicit none
2842  gs = sqrt(alphas*4.0_dp*pi)
2843 end subroutine computeqcdvariables
2844 
2845 
2846 
2847 ! QCD scale from MCFM
2848 ! Implementation into JHUGen by Ulascan Sarica, Dec. 2015
2849 subroutine evalalphas()
2850 ! use ModParameters
2851  IMPLICIT NONE
2852 #if useLHAPDF==1
2853 !--- This is simply a wrapper to the LHAPDF implementation of the running coupling alphas, in the style of the native MCFM routine
2854  DOUBLE PRECISION alphasPDF
2855  REAL(DP) :: Q
2856  q = mu_ren/gev
2857  alphas=alphaspdf(q)
2858 #else
2859 ! Evaluation of strong coupling constant alphas
2860 ! Original Author: R.K. Ellis
2861 ! q -- Scale at which alpha_s is to be evaluated
2862 ! alphas_mz -- ModParameters value of alpha_s at the mass of the Z-boson
2863 ! nloops_pdf -- ModParameters value of the number of loops (1,2, or 3) at which the beta function is evaluated to determine running.
2864 ! If you somehow need a more complete implementation, check everything at or before commit 28472c5bfee128dde458fd4929b4d3ece9519ab8
2865  INTEGER, PARAMETER :: NF6=6
2866  INTEGER, PARAMETER :: NF5=5
2867  INTEGER, PARAMETER :: NF4=4
2868  INTEGER, PARAMETER :: NF3=3
2869  INTEGER, PARAMETER :: NF2=2
2870  INTEGER, PARAMETER :: NF1=1
2871  INTEGER, PARAMETER :: NF0=0
2872 
2873  IF (mu_ren .LE. 0d0) THEN
2874  WRITE(6,*) .le.'ModParameters::EvalAlphaS: Mu_Ren 0, Mu_Ren (GeV) = ',(mu_ren*gev)
2875  stop
2876  ENDIF
2877  IF (nqflavors_pdf .GT. nf6 .OR. nqflavors_pdf .LT. nf0) THEN
2878  WRITE(6,*) 'ModParameters::EvalAlphaS: nQflavors_pdf has to be between 0 and 6. nQflavors_pdf = ',nqflavors_pdf
2879  stop
2880  ENDIF
2881  IF (nqflavors_pdf .NE. nf5) THEN
2882  WRITE(6,*) 'ModParameters::EvalAlphaS: nQflavors_pdf invalid, nQflavors_pdf = ',nqflavors_pdf
2883  WRITE(6,*) 'ModParameters::EvalAlphaS: Check 28472c5bfee128dde458fd4929b4d3ece9519ab8'
2884  stop
2885  ENDIF
2886  IF (nloops_pdf .NE. 1) THEN
2887  WRITE(6,*) 'ModParameters::EvalAlphaS: nloops_pdf invalid, nloops_pdf = ',nloops_pdf
2888  WRITE(6,*) 'ModParameters::EvalAlphaS: Check 28472c5bfee128dde458fd4929b4d3ece9519ab8'
2889  stop
2890  ENDIF
2891 
2892  alphas=alphas_mz/(1.0_dp+alphas_mz*b0_pdf(nqflavors_pdf)*2.0_dp*dlog((mu_ren/zmass_pdf)))
2893 #endif
2894  ! Calculate the derived couplings
2895  call computeqcdvariables()
2896  RETURN
2897 end subroutine evalalphas
2898 
2899 
2900 
2901 !-----------------------------------------------------------------------------
2902 !
2903  real(8) function massfrun(mf,scale)
2904 !
2905 !-----------------------------------------------------------------------------
2906 !
2907 ! This function returns the 'nloop' value of a MSbar fermion mass
2908 ! at a given scale.
2909 !
2910 ! INPUT: mf = MSbar mass of fermion at MSbar fermion mass scale
2911 ! scale = scale at which the running mass is evaluated
2912 ! asmz = AS(MZ) : this is passed to alphas(scale,asmz,2)
2913 ! nloop = # of loops in the evolutionC
2914 !
2915 ! COMMON BLOCKS: COMMON/QMASS/CMASS,BMASS,TMASS
2916 ! contains the MS-bar masses of the heavy quarks.
2917 !
2918 ! EXTERNAL: double precision alphas(scale,asmz,2)
2919 !
2920 !-----------------------------------------------------------------------------
2921 !
2922 ! use ModParameters
2923  implicit none
2924 !
2925 ! ARGUMENTS
2926 !
2927  real(8), intent(in) :: mf, scale
2928  real(8) scale_temp_ren
2929  !integer , intent(in) :: nloop
2930 !
2931 ! LOCAL
2932 !
2933  real(8) beta0, beta1,gamma0,gamma1
2934  real(8) as,asmf,l2
2935  integer nfrun
2936 
2937  scale_temp_ren=mu_ren
2938 !
2939 ! EXTERNAL
2940 !
2941 ! double precision alphas
2942 ! external alphas
2943 !
2944 ! COMMON
2945 !
2946 ! real *8 cmass,bmass,tmass
2947 ! COMMON/QMASS/CMASS,BMASS,TMASS
2948 !
2949 ! CONSTANTS
2950 !
2951 ! double precision One, Two, Three, Pi
2952  !parameter( One = 1d0, Two = 2d0, Three = 3d0 )
2953  !parameter( Pi = 3.14159265358979323846d0)
2954 
2955  if ( mf.gt.m_top ) then
2956  nfrun = 6
2957  else
2958  nfrun = 5
2959  end if
2960 
2961  beta0 = ( 11d0 - 2d0/3d0 *nfrun )/4d0
2962  !beta1 = ( 102d0 - 38d0/3d0*nf )/16d0
2963  gamma0= 1d0
2964  !gamma1= ( 202d0/3d0 - 20d0/9d0*nf )/16d0
2965  !A1 = -beta1*gamma0/beta0**2+gamma1/beta0
2966  mu_ren=scale
2967  call evalalphas()
2968  as=alphas
2969 
2970  mu_ren=mf
2971  call evalalphas()
2972  asmf=alphas
2973  !l2 = (1d0+A1*as/Pi)/(one+A1*asmf/Pi)
2974 
2975  massfrun = mf * (as/asmf)**(gamma0/beta0)
2976 
2977  !if(nloop.eq.2) massfrun=massfrun*l2
2978 
2979  mu_ren=scale_temp_ren
2980  call evalalphas()
2981 
2982  return
2983  end function massfrun
2984 
2985 
2986 
2987 
2988 
2989 
2990 
2991 
2992 
2993 !========================================================================
2994 !---- THESE ARE POLARIZATION ROUTINES
2995  ! -- massless vector polarization subroutine
2996  function pol_mless(p,i,outgoing)
2997  implicit none
2998  complex(dp), intent(in) :: p(4)
2999  integer, intent(in) :: i
3000  logical, intent(in),optional :: outgoing
3001  ! -------------------------------
3002  integer :: pol
3003  real(dp) :: p0,px,py,pz
3004  real(dp) :: pv,ct,st,cphi,sphi
3005  complex(dp) :: pol_mless(4)
3006 
3007 !^^^IFmp
3008 ! p0=(p(1)+conjg(p(1)))/two
3009 ! px=(p(2)+conjg(p(2)))/two
3010 ! py=(p(3)+conjg(p(3)))/two
3011 ! pz=(p(4)+conjg(p(4)))/two
3012 !^^^ELSE
3013  p0=real(p(1),dp)
3014  px=real(p(2),dp)
3015  py=real(p(3),dp)
3016  pz=real(p(4),dp)
3017 !^^^END
3018 
3019 
3020  pv=sqrt(abs(p0**2))
3021  ct=pz/pv
3022  st=sqrt(abs(1.0_dp-ct**2))
3023 
3024  if (st < tol) then
3025  cphi=1.0_dp
3026  sphi=0.0_dp
3027  else
3028  cphi= px/pv/st
3029  sphi= py/pv/st
3030  endif
3031 
3032 
3033  ! -- distinguish between positive and negative energies
3034  if ( p0 > 0.0_dp) then
3035  pol=i
3036  else
3037  pol=-i
3038  endif
3039 
3040  ! -- take complex conjugate for outgoing
3041  if (present(outgoing)) then
3042  if (outgoing) pol = -pol
3043  endif
3044 
3045  pol_mless(1)=czero
3046  pol_mless(2)=ct*cphi/sqrt2 - ci*pol*sphi/sqrt2
3047  pol_mless(3)=ct*sphi/sqrt2 + ci*pol*cphi/sqrt2
3048  pol_mless(4)=-st/sqrt2
3049 
3050  end function pol_mless
3051 
3052  function pol_mless2(p,i,out)
3053  implicit none
3054  integer, intent(in) :: i
3055  complex(dp), intent(in) :: p(4)
3056  character(len=*), intent(in):: out
3057  complex(dp) :: pol_mless2(4)
3058  ! -------------------------------------
3059 
3060  if (out == 'out') then
3061  pol_mless2 = pol_mless(p,i,outgoing=.true.)
3062  else
3063  pol_mless2 = pol_mless(p,i,outgoing=.false.)
3064  endif
3065  end function pol_mless2
3066 
3067  function pol_dk2mom(plepton,antilepton,i,outgoing)
3068  implicit none
3069  integer, intent(in) :: i
3070  integer :: j
3071  complex(dp), intent(in) :: plepton(1:4),antilepton(1:4)
3072  logical, intent(in),optional :: outgoing
3073  complex(dp) :: pol_dk2mom(4),ub(4),v(4),q(4),qsq
3074 
3075 
3076  q=plepton+antilepton
3077  qsq=q(1)**2-q(2)**2-q(3)**2-q(4)**2
3078 
3079  ub(:)=ubar0(plepton,i)
3080  v(:)=v0(antilepton,-i)
3081  !---Now return in Kirill's notation 1=E,2=px,3=py,4=pz
3082  ! This is an expression for (-i)/qsq* (-i) Ub(+/-)) Gamma^\mu V(-/+)
3083  pol_dk2mom(1)=-(ub(2)*v(4)+v(2)*ub(4)+ub(1)*v(3)+v(1)*ub(3))
3084  pol_dk2mom(2)=-(-ub(1)*v(4)+v(1)*ub(4)-ub(2)*v(3)+v(2)*ub(3))
3085  pol_dk2mom(3)=-ci*(ub(1)*v(4)+v(1)*ub(4)-ub(2)*v(3)-v(2)*ub(3))
3086  pol_dk2mom(4)=-(ub(2)*v(4)-v(2)*ub(4)-ub(1)*v(3)+v(1)*ub(3))
3087 
3088 
3089  do j=1,4
3090  pol_dk2mom(j)=pol_dk2mom(j)/qsq
3091  enddo
3092 
3093  ! -- do nothing in this case
3094  if (present(outgoing)) then
3095  !if (outgoing) pol_dk2mom = conjg(pol_dk2mom)
3096  endif
3097 
3098  end function pol_dk2mom
3099 
3100 ! -- massive vector polarization subroutine
3101  function pol_mass(p,i,outgoing)
3102  implicit none
3103  integer, intent(in) :: i
3104  integer :: pol
3105  complex(8), intent(in) :: p(4)
3106  logical, intent(in),optional :: outgoing
3107  complex(8) :: pol_mass(4)
3108  complex(8) :: msq,m
3109  real(8) :: p0,px,py,pz, pv,pvsq
3110  real(8) :: ct,st,cphi,sphi
3111 
3112  p0=dreal(p(1))
3113  px=dreal(p(2))
3114  py=dreal(p(3))
3115  pz=dreal(p(4))
3116 
3117  pv=px**2 + py**2 + pz**2
3118  m=p0**2-pv
3119  m=sqrt(m)
3120  pv=sqrt(pv)
3121 
3122  if(cdabs(pv/m).lt.1d-8) then
3123  if(i.eq.0) then
3124  pol_mass(1:3)=czero
3125  pol_mass( 4 )=cone
3126  return
3127  endif
3128  ct = 1d0; st=0d0
3129  else
3130  ct= pz/pv
3131  st= dsqrt(dabs(1.0d0-ct**2))
3132  endif
3133 
3134 
3135  if (st .lt. 1d-15) then
3136  cphi=1.0d0
3137  sphi=0d0
3138  else
3139  cphi= px/pv/st
3140  sphi= py/pv/st
3141  endif
3142 
3143 
3144 ! i=0 is longitudinal polarization
3145 ! the following ifstatement distinguishes between
3146 ! positive and negative energies
3147  if ( p0 .gt. 0.0d0) then
3148  pol=i
3149  else
3150  pol=-i
3151  endif
3152 
3153  ! -- take complex conjugate for outgoing
3154  if (present(outgoing)) then
3155  if (outgoing) pol = -pol
3156  endif
3157 
3158  if(pol.eq.-1 .or. pol.eq.1) then
3159  pol_mass(1)=czero
3160  pol_mass(2)=(ct*cphi-pol*ci*sphi)/sqrt2
3161  pol_mass(3)=(ct*sphi+pol*ci*cphi)/sqrt2
3162  pol_mass(4)=-st/sqrt2
3163  else if(pol.eq.0) then
3164  pol_mass(1)= pv/m
3165  pol_mass(2)= p0/m/pv*px
3166  pol_mass(3)= p0/m/pv*py
3167  pol_mass(4)= p0/m/pv*pz
3168  else
3169  print *,"wrong helicity setting in pol_mass"
3170  stop
3171  endif
3172 
3173  end function pol_mass
3174 
3175  function pol_mass2(p,i,out)
3176  implicit none
3177  integer, intent(in) :: i
3178  complex(dp), intent(in) :: p(4)
3179  character(len=*), intent(in):: out
3180  complex(dp) :: pol_mass2(4)
3181  ! -------------------------------------
3182 
3183  if (out == 'out') then
3184  pol_mass2 = pol_mass(p,i,outgoing=.true.)
3185  else
3186  pol_mass2 = pol_mass(p,i,outgoing=.false.)
3187  endif
3188  end function pol_mass2
3189 
3190 !---- THESE ARE SPINOR ROUTINES
3191 ! ubar spinor, massless
3192  function ubar0(p,i)
3193  implicit none
3194  complex(dp), intent(in) :: p(4)
3195  integer, intent(in) :: i
3196  complex(dp) :: ubar0(4)
3197  complex(dp) :: fc, fc2
3198  real(dp) :: p0,px,py,pz,mass
3199 
3200 
3201  p0=real(p(1),dp)
3202  px=real(p(2),dp)
3203  py=real(p(3),dp)
3204  pz=real(p(4),dp)
3205  mass=dsqrt(dabs(p0**2-px**2-py**2-pz**2))
3206  if( mass.lt.1d-4 ) mass=0d0
3207 
3208 
3209  fc2 = p0 + pz
3210  fc=sqrt(fc2)
3211 
3212  if (abs(fc2).gt. tol) then
3213  if (i.eq.1) then
3214  ubar0(1)=czero
3215  ubar0(2)=czero
3216  ubar0(3)=fc
3217  ubar0(4)=(px-ci*py)/fc
3218  elseif (i.eq.-1) then
3219  ubar0(1)=(px+ci*py)/fc
3220  ubar0(2)=-fc
3221  ubar0(3)=czero
3222  ubar0(4)=czero
3223  else
3224  stop 'ubar0: i out of range'
3225  endif
3226  else
3227  if (i.eq.1) then
3228  ubar0(1) = czero
3229  ubar0(2) = czero
3230  ubar0(3) = czero
3231  ubar0(4) = sqrt(cone*two*p0)
3232  elseif (i.eq.-1) then
3233  ubar0(1) = sqrt(cone*(two*p0))
3234  ubar0(2) = czero
3235  ubar0(3) = czero
3236  ubar0(4) = czero
3237  else
3238  stop 'ubar0: i out of range'
3239  endif
3240  endif
3241 
3242 
3243 ! if (i.eq.1) then
3244 ! ubar0(1)=dcmplx(mass,0d0)/fc
3245 ! ubar0(2)=czero
3246 ! ubar0(3)=fc
3247 ! ubar0(4)=dcmplx(px,-py)/fc
3248 ! elseif (i.eq.-1) then
3249 ! ubar0(1)=dcmplx(px,py)/fc
3250 ! ubar0(2)=-fc
3251 ! ubar0(3)=czero
3252 ! ubar0(4)=-dcmplx(mass,0d0)/fc
3253 ! else
3254 ! stop 'ubar0: i out of range'
3255 ! endif
3256 
3257 
3258 
3259  end function ubar0
3260 
3261  ! -- v0 spinor, massless
3262  function v0(p,i)
3263  implicit none
3264  complex(dp), intent(in) :: p(4)
3265  integer, intent(in) :: i
3266  complex(dp) :: v0(4)
3267  complex(dp) :: fc2, fc
3268  real(dp) :: p0,px,py,pz,mass
3269 
3270  p0=real(p(1),dp)
3271  px=real(p(2),dp)
3272  py=real(p(3),dp)
3273  pz=real(p(4),dp)
3274  mass=dsqrt(dabs(p0**2-px**2-py**2-pz**2))
3275  if( mass.lt.1d-4 ) mass=0d0
3276 
3277 
3278  fc2 = p0 + pz
3279  fc=sqrt(fc2)
3280 
3281  if (abs(fc2).gt. tol) then
3282  if (i.eq.1) then
3283  v0(1)=czero
3284  v0(2)=czero
3285  v0(3)=(px-ci*py)/fc
3286  v0(4)=-fc
3287  elseif (i.eq.-1) then
3288  v0(1)=fc
3289  v0(2)=(px+ci*py)/fc
3290  v0(3)=czero
3291  v0(4)=czero
3292  else
3293  stop 'v0: i out of range'
3294  endif
3295  else
3296  if (i.eq.1) then
3297  v0(1)=czero
3298  v0(2)=czero
3299  v0(3)=sqrt(cone*two*p0)
3300  v0(4)=czero
3301  elseif (i.eq.-1) then
3302  v0(1)=czero
3303  v0(2)=sqrt(cone*two*p0)
3304  v0(3)=czero
3305  v0(4)=czero
3306  else
3307  stop 'v0: i out of range'
3308  endif
3309  endif
3310 
3311 
3312 ! if (i.eq.+1) then
3313 ! v0(1)=czero
3314 ! v0(2)=dcmplx(mass,0d0)/fc
3315 ! v0(3)=dcmplx(px,-py)/fc
3316 ! v0(4)=-fc
3317 ! elseif (i.eq.-1) then
3318 ! v0(1)=fc
3319 ! v0(2)=dcmplx(px,py)/fc
3320 ! v0(3)=dcmplx(-mass,0d0)/fc
3321 ! v0(4)=czero
3322 ! else
3323 ! stop 'v0: i out of range'
3324 ! endif
3325 
3326 
3327 
3328  end function v0
3329 
3330 subroutine spinoru(p,za,zb,s)
3331 !---Calculate spinor products
3332 !---taken from MCFM & modified by R. Rontsch, May 2015
3333 !---extended to deal with negative energies ie with all momenta outgoing
3334 !---Arbitrary conventions of Bern, Dixon, Kosower, Weinzierl,
3335 !---za(i,j)*zb(j,i)=s(i,j)
3336  implicit none
3337  real(dp) :: p(:,:),two
3338  integer, parameter :: mxpart=14
3339  complex(dp):: c23(mxpart),f(mxpart),rt(mxpart),za(:,:),zb(:,:),czero,cone,ci
3340  real(dp) :: s(:,:)
3341  integer i,j,N
3342 
3343  n = size(p,1)
3344 ! if (size(p,1) .ne. N) then
3345 ! call Error("spinorz: momentum mismatch",size(p,1))
3346 ! endif
3347  two=2d0
3348  czero=dcmplx(0d0,0d0)
3349  cone=dcmplx(1d0,0d0)
3350  ci=dcmplx(0d0,1d0)
3351 
3352 
3353 !---if one of the vectors happens to be zero this routine fails.
3354  do j=1,n
3355  za(j,j)=czero
3356  zb(j,j)=za(j,j)
3357 
3358 !-----positive energy case
3359  if (p(j,4) .gt. 0d0) then
3360  rt(j)=dsqrt(p(j,4)+p(j,1))
3361  c23(j)=dcmplx(p(j,3),-p(j,2))
3362  f(j)=cone
3363  else
3364 !-----negative energy case
3365  rt(j)=dsqrt(-p(j,4)-p(j,1))
3366  c23(j)=dcmplx(-p(j,3),p(j,2))
3367  f(j)=ci
3368  endif
3369  enddo
3370  do i=2,n
3371  do j=1,i-1
3372  s(i,j)=two*(p(i,4)*p(j,4)-p(i,1)*p(j,1)-p(i,2)*p(j,2)-p(i,3)*p(j,3))
3373  za(i,j)=f(i)*f(j)*(c23(i)*dcmplx(rt(j)/rt(i))-c23(j)*dcmplx(rt(i)/rt(j)))
3374 
3375  if (abs(s(i,j)).lt.1d-5) then
3376  zb(i,j)=-(f(i)*f(j))**2*dconjg(za(i,j))
3377  else
3378  zb(i,j)=-dcmplx(s(i,j))/za(i,j)
3379  endif
3380  za(j,i)=-za(i,j)
3381  zb(j,i)=-zb(i,j)
3382  s(j,i)=s(i,j)
3383  enddo
3384  enddo
3385 
3386  end subroutine spinoru
3387 
3388 
3389 
3390 
3391 !========================================================================
3392 
3393 !========================================================================
3394 ! Common initialization functions that may be called multiple times if needed
3395 ! Check arXiv:1604.06792 for the parameters
3396 subroutine initcollier(Nmax, Rmax)
3397 #if useCollier==1
3398 use collier
3399 implicit none
3400 integer, intent(in) :: Nmax, Rmax
3401 integer :: supNmax, supRmax
3402  supnmax = max(nmax, collier_maxnloopprops)
3403  suprmax = max(rmax, collier_maxrank)
3404  if ((supnmax .gt. collier_maxnloopprops) .or. (suprmax .gt. collier_maxrank)) then
3405  call init_cll(supnmax,suprmax,'')
3406  call setmode_cll(1)
3407  endif
3408 #else
3409 implicit none
3410 integer, intent(in) :: Nmax, Rmax
3411  return
3412 #endif
3413 end subroutine
3414 
3415 
3416 
3417 END MODULE
modparameters::lambda2_w42
real(8), public lambda2_w42
Definition: mod_Parameters.F90:858
modparameters::b5
complex(8), public b5
Definition: mod_Parameters.F90:940
modparameters::vbfoffsh_run
integer, public vbfoffsh_run
Definition: mod_Parameters.F90:144
modparameters::abot_
integer, target, public abot_
Definition: mod_Parameters.F90:1111
modparameters::b1
complex(8), public b1
Definition: mod_Parameters.F90:936
modparameters::b0_pdf
real(dp), dimension(0:6), parameter, public b0_pdf
Definition: mod_Parameters.F90:1013
modparameters::gh2w4_prime
complex(8), public gh2w4_prime
Definition: mod_Parameters.F90:835
modparameters::ghzpzp1_prime7
complex(8), public ghzpzp1_prime7
Definition: mod_Parameters.F90:571
modparameters::gamma_0
real(dp), parameter, public gamma_0
Definition: mod_Parameters.F90:81
modparameters::kappa_tilde_top
complex(8), public kappa_tilde_top
Definition: mod_Parameters.F90:887
modparameters::ezp_top_right
complex(8), public ezp_top_right
Definition: mod_Parameters.F90:617
modparameters::lambda2_w3
real(8), parameter, public lambda2_w3
Definition: mod_Parameters.F90:845
modparameters::alphas_mz
real(dp), public alphas_mz
Definition: mod_Parameters.F90:270
modparameters::lambda_w31
real(8), public lambda_w31
Definition: mod_Parameters.F90:502
modparameters::bzpgs8
complex(8), public bzpgs8
Definition: mod_Parameters.F90:985
modparameters::zpr2_
integer, target, public zpr2_
Definition: mod_Parameters.F90:1102
modparameters::ghzpzp4_prime4
complex(8), public ghzpzp4_prime4
Definition: mod_Parameters.F90:592
modparameters::lambda2_z11
real(8), public lambda2_z11
Definition: mod_Parameters.F90:788
modparameters::ghwwp4_prime5
complex(8), public ghwwp4_prime5
Definition: mod_Parameters.F90:657
modparameters::gh2w1_prime2
complex(8), public gh2w1_prime2
Definition: mod_Parameters.F90:812
modparameters::m_vprime
real(8), public m_vprime
Definition: mod_Parameters.F90:78
modparameters::achm_
integer, target, public achm_
Definition: mod_Parameters.F90:1108
modparameters::includegammastar
logical, public includegammastar
Definition: mod_Parameters.F90:213
modparameters::ghz4_prime6
complex(8), public ghz4_prime6
Definition: mod_Parameters.F90:416
modparameters::debugcounter
integer, dimension(0:10), public debugcounter
Definition: mod_Parameters.F90:1148
modparameters::ijpartons
integer, dimension(1:2), public ijpartons
Definition: mod_Parameters.F90:1152
modparameters::ghzzp1_prime
complex(8), public ghzzp1_prime
Definition: mod_Parameters.F90:520
modparameters::nloops_pdf
integer, parameter, public nloops_pdf
Definition: mod_Parameters.F90:265
modparameters::daz_neu_right
real(8), public daz_neu_right
Definition: mod_Parameters.F90:1002
modparameters::glu_
integer, target, public glu_
Definition: mod_Parameters.F90:1093
modparameters::bzpgs3
complex(8), public bzpgs3
Definition: mod_Parameters.F90:983
modparameters::a5
complex(8), public a5
Definition: mod_Parameters.F90:927
modparameters::pdfstr_
integer, parameter, public pdfstr_
Definition: mod_Parameters.F90:1127
modparameters::br_z_nn
real(8), parameter, public br_z_nn
Definition: mod_Parameters.F90:304
modparameters::hasleptonfilter
logical, public hasleptonfilter
Definition: mod_Parameters.F90:29
modparameters::vegasnc2
integer, public vegasnc2
Definition: mod_Parameters.F90:18
modparameters::ezp_top_left
complex(8), public ezp_top_left
Definition: mod_Parameters.F90:616
modparameters::vev
real(8), public vev
Definition: mod_Parameters.F90:249
modparameters::gh2z1_prime7
complex(8), public gh2z1_prime7
Definition: mod_Parameters.F90:750
modparameters::elm_
integer, target, public elm_
Definition: mod_Parameters.F90:1112
modparameters::gh2zgs2
complex(8), public gh2zgs2
Definition: mod_Parameters.F90:736
modparameters::lambda_w1
real(8), parameter, public lambda_w1
Definition: mod_Parameters.F90:492
modparameters::ghzpzp2_prime7
complex(8), public ghzpzp2_prime7
Definition: mod_Parameters.F90:579
modparameters::gh2w2_prime5
complex(8), public gh2w2_prime5
Definition: mod_Parameters.F90:823
modparameters::twosc
real(8), public twosc
Definition: mod_Parameters.F90:255
modparameters::ghz2_prime
complex(8), public ghz2_prime
Definition: mod_Parameters.F90:395
modparameters::couplazff
real(dp), public couplazff
Definition: mod_Parameters.F90:1080
modparameters::ghw1
complex(8), public ghw1
Definition: mod_Parameters.F90:454
modparameters::krenfacscheme_mj
integer, parameter, public krenfacscheme_mj
Definition: mod_Parameters.F90:41
modparameters::widthschemein
integer, public widthschemein
Definition: mod_Parameters.F90:132
modparameters::daz_qdn_left
real(8), public daz_qdn_left
Definition: mod_Parameters.F90:1005
modparameters::pdfastr_
integer, parameter, public pdfastr_
Definition: mod_Parameters.F90:1133
modparameters::ghzzp1_prime7
complex(8), public ghzzp1_prime7
Definition: mod_Parameters.F90:526
modparameters::ewp_chm_right
complex(8), public ewp_chm_right
Definition: mod_Parameters.F90:708
modparameters::ghz4_prime4
complex(8), public ghz4_prime4
Definition: mod_Parameters.F90:414
modparameters::zffbare
real(8) function zffbare(id1in, id2in, h1, h2)
Definition: mod_Parameters.F90:1601
modparameters::ghwpwp3
complex(8), public ghwpwp3
Definition: mod_Parameters.F90:663
modparameters::kappa_tilde_4gen_bot
complex(8), public kappa_tilde_4gen_bot
Definition: mod_Parameters.F90:894
modparameters::zpgsmode
integer, parameter, public zpgsmode
Definition: mod_Parameters.F90:16
modparameters::ghzpzp1
complex(8), public ghzpzp1
Definition: mod_Parameters.F90:559
modparameters::ghz3_prime5
complex(8), public ghz3_prime5
Definition: mod_Parameters.F90:407
modparameters::pdfset
integer, public pdfset
Definition: mod_Parameters.F90:69
modparameters::dm_z
complex(8), public dm_z
Definition: mod_Parameters.F90:873
modparameters::krenfacscheme_mhstar
integer, parameter, public krenfacscheme_mhstar
Definition: mod_Parameters.F90:33
modparameters::c2w_q1sq
integer, public c2w_q1sq
Definition: mod_Parameters.F90:848
modparameters::hvvspinzerodynamiccoupling
complex(8) function hvvspinzerodynamiccoupling(index, sWplus, sWminus, sWW, tryWWcoupl)
Definition: mod_Parameters.F90:1161
modparameters::ghwwp4_prime3
complex(8), public ghwwp4_prime3
Definition: mod_Parameters.F90:655
modparameters::ghg2
complex(8), public ghg2
Definition: mod_Parameters.F90:371
modparameters::gs
real(dp), public gs
Definition: mod_Parameters.F90:271
modparameters::brlept_z_mm
real(8), parameter, public brlept_z_mm
Definition: mod_Parameters.F90:312
modparameters::mphotoncutoff
real(8), public mphotoncutoff
Definition: mod_Parameters.F90:215
modparameters::distinguish_hwwcouplings
logical, public distinguish_hwwcouplings
Definition: mod_Parameters.F90:453
modparameters::lhapdfmember
integer, public lhapdfmember
Definition: mod_Parameters.F90:68
modparameters::defaultseeds
integer, dimension(1:nmaxseeds), parameter, public defaultseeds
Definition: mod_Parameters.F90:160
modparameters::nmaxchannels
integer, parameter, public nmaxchannels
Definition: mod_Parameters.F90:47
modparameters::vprimedecaylengthmm
real(8), public vprimedecaylengthmm
Definition: mod_Parameters.F90:236
modparameters::pdfatop_
integer, parameter, public pdfatop_
Definition: mod_Parameters.F90:1136
modparameters::lambda_z10
real(8), public lambda_z10
Definition: mod_Parameters.F90:446
modparameters::daz_neu_left
real(8), public daz_neu_left
Definition: mod_Parameters.F90:1001
modparameters::bzpzp1
complex(8), public bzpzp1
Definition: mod_Parameters.F90:958
modparameters::tev
real(8), parameter, public tev
Definition: mod_Parameters.F90:94
modparameters::ghzzp3_prime
complex(8), public ghzzp3_prime
Definition: mod_Parameters.F90:536
modparameters::isagluon
logical function isagluon(PartType)
Definition: mod_Parameters.F90:2423
modparameters::ghw2_prime6
complex(8), public ghw2_prime6
Definition: mod_Parameters.F90:473
modparameters::pdfglu_
integer, parameter, public pdfglu_
Definition: mod_Parameters.F90:1124
modparameters::coupledvertex
integer function coupledvertex(id, hel, useAHcoupl)
Definition: mod_Parameters.F90:2467
modparameters::pdftop_
integer, parameter, public pdftop_
Definition: mod_Parameters.F90:1130
modparameters::c7
complex(8), parameter, public c7
Definition: mod_Parameters.F90:995
modparameters::bzpzp6
complex(8), public bzpzp6
Definition: mod_Parameters.F90:963
modparameters::tev_energy
real(8), parameter, public tev_energy
Definition: mod_Parameters.F90:257
modparameters::ghzzp4_prime2
complex(8), public ghzzp4_prime2
Definition: mod_Parameters.F90:545
modparameters::gh2z4_prime
complex(8), public gh2z4_prime
Definition: mod_Parameters.F90:768
modparameters::ghg3_4gen
complex(8), public ghg3_4gen
Definition: mod_Parameters.F90:377
modparameters::zgmode
integer, parameter, public zgmode
Definition: mod_Parameters.F90:13
modparameters::ghwwp2_prime7
complex(8), public ghwwp2_prime7
Definition: mod_Parameters.F90:643
modparameters::ghwwp4_prime4
complex(8), public ghwwp4_prime4
Definition: mod_Parameters.F90:656
modparameters::krenfacscheme_mj_mj
integer, parameter, public krenfacscheme_mj_mj
Definition: mod_Parameters.F90:38
modparameters::m_tau
real(8), public m_tau
Definition: mod_Parameters.F90:242
modparameters::lambda_z41
real(8), public lambda_z41
Definition: mod_Parameters.F90:441
modparameters::mvh_minmax
real(8), dimension(1:2), public mvh_minmax
Definition: mod_Parameters.F90:212
modparameters::lambda2_z10
real(8), public lambda2_z10
Definition: mod_Parameters.F90:796
modparameters::ptjetcut
real(8), public ptjetcut
Definition: mod_Parameters.F90:205
modparameters::mubarh
real(8), public mubarh
Definition: mod_Parameters.F90:133
modparameters::vckm_cb
real(8), public vckm_cb
Definition: mod_Parameters.F90:282
modparameters::bzgs1
complex(8), public bzgs1
Definition: mod_Parameters.F90:969
modparameters::lambda2_w4
real(8), parameter, public lambda2_w4
Definition: mod_Parameters.F90:846
modparameters::ghwpwp1
complex(8), public ghwpwp1
Definition: mod_Parameters.F90:661
modparameters::cl_neu
real(8), public cl_neu
Definition: mod_Parameters.F90:1071
modparameters::io_csmaxfile
integer, parameter io_csmaxfile
Definition: mod_Parameters.F90:1143
modparameters::jetsinevent
integer, dimension(0:8), public jetsinevent
Definition: mod_Parameters.F90:128
modparameters::kappa2_tilde_4gen_top
complex(8), public kappa2_tilde_4gen_top
Definition: mod_Parameters.F90:903
modparameters::pmzzevals
integer, public pmzzevals
Definition: mod_Parameters.F90:18
modparameters::doprintpmzz
logical, public doprintpmzz
Definition: mod_Parameters.F90:28
modparameters::b10
complex(8), public b10
Definition: mod_Parameters.F90:945
modparameters::gh2w3_prime4
complex(8), public gh2w3_prime4
Definition: mod_Parameters.F90:830
modparameters::reweightdecay
logical, public reweightdecay
Definition: mod_Parameters.F90:129
modparameters::ezp_up_left
complex(8), public ezp_up_left
Definition: mod_Parameters.F90:606
modparameters::gh2z3_prime2
complex(8), public gh2z3_prime2
Definition: mod_Parameters.F90:761
modparameters::krenfacscheme_custom_scale
integer, parameter, public krenfacscheme_custom_scale
Definition: mod_Parameters.F90:44
modparameters::gh2w2
complex(8), public gh2w2
Definition: mod_Parameters.F90:806
modparameters::ghzgs4
complex(8), public ghzgs4
Definition: mod_Parameters.F90:422
modparameters::pmzzfile
character(len=500), public pmzzfile
Definition: mod_Parameters.F90:137
modparameters::b8
complex(8), public b8
Definition: mod_Parameters.F90:943
modparameters::zprime_zz_1
complex(8), public zprime_zz_1
Definition: mod_Parameters.F90:914
modparameters::vckm_cs
real(8), public vckm_cs
Definition: mod_Parameters.F90:278
modparameters::vegasit1
integer, public vegasit1
Definition: mod_Parameters.F90:18
modparameters::lambda_z12
real(8), public lambda_z12
Definition: mod_Parameters.F90:442
modparameters::getlheparticle
character(len=3) function getlheparticle(PartLHE)
Definition: mod_Parameters.F90:2229
modparameters::ghg3
complex(8), public ghg3
Definition: mod_Parameters.F90:372
modparameters::aveqq
real(dp), parameter, public aveqq
Definition: mod_Parameters.F90:1021
modparameters::br_counter
integer, dimension(1:5, 1:5), public br_counter
Definition: mod_Parameters.F90:126
modparameters::bzpgs2
complex(8), public bzpgs2
Definition: mod_Parameters.F90:982
modparameters::gh2z2_prime
complex(8), public gh2z2_prime
Definition: mod_Parameters.F90:752
modparameters::mufacmultiplier
real(8), public mufacmultiplier
Definition: mod_Parameters.F90:21
modparameters::offshellreson
logical, public offshellreson
Definition: mod_Parameters.F90:28
modparameters::ckm
real(8) function ckm(id1in, id2in)
Definition: mod_Parameters.F90:1593
modparameters::ghzpzp2_prime6
complex(8), public ghzpzp2_prime6
Definition: mod_Parameters.F90:578
modparameters::ghz3_prime2
complex(8), public ghz3_prime2
Definition: mod_Parameters.F90:404
modparameters::dzzwpwm
complex(8), public dzzwpwm
Definition: mod_Parameters.F90:876
collier
Definition: COLLIER.F90:30
modparameters::ghz2
complex(8), public ghz2
Definition: mod_Parameters.F90:382
modparameters::ghw2_prime7
complex(8), public ghw2_prime7
Definition: mod_Parameters.F90:474
modparameters::ghzzp2_prime5
complex(8), public ghzzp2_prime5
Definition: mod_Parameters.F90:532
modparameters::a2
complex(8), public a2
Definition: mod_Parameters.F90:924
modparameters::ga_wprime
real(8), public ga_wprime
Definition: mod_Parameters.F90:713
modparameters::ghwwp4_prime7
complex(8), public ghwwp4_prime7
Definition: mod_Parameters.F90:659
modparameters::maxpart
integer, parameter, public maxpart
Definition: mod_Parameters.F90:46
modparameters::brhadr_z_ss
real(8), parameter, public brhadr_z_ss
Definition: mod_Parameters.F90:322
modparameters::lambda_z21
real(8), public lambda_z21
Definition: mod_Parameters.F90:439
modparameters::counttauasany
logical, public counttauasany
Definition: mod_Parameters.F90:29
modparameters::mu_fact
real(8), public mu_fact
Definition: mod_Parameters.F90:267
modparameters::ghzpzp2
complex(8), public ghzpzp2
Definition: mod_Parameters.F90:560
modparameters::theseeds
integer, dimension(1:nmaxseeds), public theseeds
Definition: mod_Parameters.F90:161
modparameters::ghzzp4_prime3
complex(8), public ghzzp4_prime3
Definition: mod_Parameters.F90:546
modparameters::ghwpwp3_prime6
complex(8), public ghwpwp3_prime6
Definition: mod_Parameters.F90:688
modparameters::customlowerscalebound
real(8), public customlowerscalebound
Definition: mod_Parameters.F90:21
modparameters::kappa_4gen_bot
complex(8), public kappa_4gen_bot
Definition: mod_Parameters.F90:893
modparameters::bzpzp8
complex(8), public bzpzp8
Definition: mod_Parameters.F90:965
modparameters::dn_
integer, target, public dn_
Definition: mod_Parameters.F90:1085
modparameters::detajetcut
real(8), public detajetcut
Definition: mod_Parameters.F90:207
modparameters::murenmultiplier
real(8), public murenmultiplier
Definition: mod_Parameters.F90:21
modparameters::brlept_w_mn
real(8), parameter, public brlept_w_mn
Definition: mod_Parameters.F90:316
modparameters::ahz2
complex(8), parameter, public ahz2
Definition: mod_Parameters.F90:366
modparameters::ghw2_prime3
complex(8), public ghw2_prime3
Definition: mod_Parameters.F90:470
modparameters::ahz3
complex(8), parameter, public ahz3
Definition: mod_Parameters.F90:367
modparameters::logfile
character(len=500) logfile
Definition: mod_Parameters.F90:63
modparameters::bzgs4
complex(8), public bzgs4
Definition: mod_Parameters.F90:972
modparameters::ghzzp4_prime4
complex(8), public ghzzp4_prime4
Definition: mod_Parameters.F90:547
modparameters::writegit
logical, public writegit
Definition: mod_Parameters.F90:77
modparameters::ga_v_ps
real(8), public ga_v_ps
Definition: mod_Parameters.F90:78
modparameters::m2l_minmax
real(8), dimension(1:2), public m2l_minmax
Definition: mod_Parameters.F90:211
modparameters::kappa_bot
complex(8), public kappa_bot
Definition: mod_Parameters.F90:888
modparameters::ghzpzp3_prime5
complex(8), public ghzpzp3_prime5
Definition: mod_Parameters.F90:585
modparameters::bzzp2
complex(8), public bzzp2
Definition: mod_Parameters.F90:948
modparameters::lambda_z31
real(8), public lambda_z31
Definition: mod_Parameters.F90:440
modparameters::ghzpzp3_prime2
complex(8), public ghzpzp3_prime2
Definition: mod_Parameters.F90:582
modparameters::pol_mless2
complex(dp) function, dimension(4) pol_mless2(p, i, out)
Definition: mod_Parameters.F90:3053
modparameters::c2z_q12sq
integer, public c2z_q12sq
Definition: mod_Parameters.F90:787
modparameters::m_reso2
real(8), public m_reso2
Definition: mod_Parameters.F90:232
modparameters::spinoru
subroutine spinoru(p, za, zb, s)
Definition: mod_Parameters.F90:3331
modparameters::qul
real(8), parameter, public qul
Definition: mod_Parameters.F90:1036
modparameters::ghz3_prime6
complex(8), public ghz3_prime6
Definition: mod_Parameters.F90:408
modparameters::ga_w
real(8), public ga_w
Definition: mod_Parameters.F90:229
modparameters::initcollier
subroutine initcollier(Nmax, Rmax)
Definition: mod_Parameters.F90:3397
modparameters::hbarc2_pbgev2
real(8), parameter, public hbarc2_pbgev2
Definition: mod_Parameters.F90:113
modparameters::gh2z2
complex(8), public gh2z2
Definition: mod_Parameters.F90:732
modparameters::ghwpwp2_prime2
complex(8), public ghwpwp2_prime2
Definition: mod_Parameters.F90:676
modparameters::ghzgs2
complex(8), public ghzgs2
Definition: mod_Parameters.F90:420
modparameters::pdfdn_
integer, parameter, public pdfdn_
Definition: mod_Parameters.F90:1125
modparameters::scale_alpha_w_tn
real(8), public scale_alpha_w_tn
Definition: mod_Parameters.F90:335
modparameters::t3ll
real(8), parameter, public t3ll
Definition: mod_Parameters.F90:1024
modparameters::crosssec
real(8), dimension(-6:+6,-6:+6), public crosssec
Definition: mod_Parameters.F90:53
modparameters::pmzzsize
integer, parameter, public pmzzsize
Definition: mod_Parameters.F90:139
modparameters::lambda2_w10
real(8), public lambda2_w10
Definition: mod_Parameters.F90:859
modparameters::kappa2_4gen_bot
complex(8), public kappa2_4gen_bot
Definition: mod_Parameters.F90:904
modparameters::topdecays
integer, public topdecays
Definition: mod_Parameters.F90:17
modparameters::krenfacscheme_mjjhstar
integer, parameter, public krenfacscheme_mjjhstar
Definition: mod_Parameters.F90:34
modparameters::zpzpmode
integer, parameter, public zpzpmode
Definition: mod_Parameters.F90:15
modparameters::ghwpwp3_prime
complex(8), public ghwpwp3_prime
Definition: mod_Parameters.F90:683
modparameters::gh2z2_prime6
complex(8), public gh2z2_prime6
Definition: mod_Parameters.F90:757
modparameters::gh2w1_prime7
complex(8), public gh2w1_prime7
Definition: mod_Parameters.F90:817
modparameters::ezp_el_left
complex(8), public ezp_el_left
Definition: mod_Parameters.F90:598
modparameters::ghwwp2_prime5
complex(8), public ghwwp2_prime5
Definition: mod_Parameters.F90:641
modparameters::gh2g3_4gen
complex(8), public gh2g3_4gen
Definition: mod_Parameters.F90:727
modparameters::ghzpzp1_prime5
complex(8), public ghzpzp1_prime5
Definition: mod_Parameters.F90:569
modparameters::cz_q12sq
integer, public cz_q12sq
Definition: mod_Parameters.F90:436
modparameters::ahg1
complex(8), parameter, public ahg1
Definition: mod_Parameters.F90:362
modparameters::gh2z2_prime2
complex(8), public gh2z2_prime2
Definition: mod_Parameters.F90:753
modparameters::ghz1_prime
complex(8), public ghz1_prime
Definition: mod_Parameters.F90:387
modparameters::cf
real(dp), parameter, public cf
Definition: mod_Parameters.F90:1017
modparameters::convertlhereverse
integer function convertlhereverse(Part)
Definition: mod_Parameters.F90:1711
modparameters::gh2z4_prime6
complex(8), public gh2z4_prime6
Definition: mod_Parameters.F90:773
modparameters::crosssecmax
real(8), dimension(-6:+6,-6:+6), public crosssecmax
Definition: mod_Parameters.F90:53
modparameters::ahg2
complex(8), parameter, public ahg2
Definition: mod_Parameters.F90:363
modparameters::krenfacscheme_default
integer, parameter, public krenfacscheme_default
Definition: mod_Parameters.F90:32
modparameters::brlept_w_en
real(8), parameter, public brlept_w_en
Definition: mod_Parameters.F90:315
modparameters::gh2gsgs2
complex(8), public gh2gsgs2
Definition: mod_Parameters.F90:739
modparameters::gh2z1
complex(8), public gh2z1
Definition: mod_Parameters.F90:731
modparameters::sitw
real(8), public sitw
Definition: mod_Parameters.F90:254
modparameters::ghwwp2
complex(8), public ghwwp2
Definition: mod_Parameters.F90:624
modparameters::one
real(8), parameter, public one
Definition: mod_Parameters.F90:83
modparameters::lambda_z2
real(8), parameter, public lambda_z2
Definition: mod_Parameters.F90:428
modparameters::gh2zgs4
complex(8), public gh2zgs4
Definition: mod_Parameters.F90:738
modparameters::fix_channels_ratio
logical, public fix_channels_ratio
Definition: mod_Parameters.F90:168
modparameters::zgsmode
integer, parameter, public zgsmode
Definition: mod_Parameters.F90:13
modparameters::convertfrompartindex
integer function convertfrompartindex(Part)
Definition: mod_Parameters.F90:1922
modparameters::sqrt2
real(dp), parameter, public sqrt2
Definition: mod_Parameters.F90:80
modparameters::m_zprime
real(8), public m_zprime
Definition: mod_Parameters.F90:619
modparameters::time_end
real(8) time_end
Definition: mod_Parameters.F90:59
modparameters::gh2z1_prime5
complex(8), public gh2z1_prime5
Definition: mod_Parameters.F90:748
modparameters::ghzpzp3_prime
complex(8), public ghzpzp3_prime
Definition: mod_Parameters.F90:581
modparameters::ezp_bot_right
complex(8), public ezp_bot_right
Definition: mod_Parameters.F90:615
modparameters::isajet
logical function isajet(PartType)
Definition: mod_Parameters.F90:2438
modparameters::scale_alpha_w_ln
real(8), public scale_alpha_w_ln
Definition: mod_Parameters.F90:334
modparameters::lambda2_z1
real(8), parameter, public lambda2_z1
Definition: mod_Parameters.F90:778
modparameters::br_w_en
real(8), parameter, public br_w_en
Definition: mod_Parameters.F90:305
modparameters::ewp_top_right
complex(8), public ewp_top_right
Definition: mod_Parameters.F90:710
modparameters::bzpzp3
complex(8), public bzpzp3
Definition: mod_Parameters.F90:960
modparameters::two
real(8), parameter, public two
Definition: mod_Parameters.F90:84
modparameters::couplaffsq
real(dp), public couplaffsq
Definition: mod_Parameters.F90:1081
modparameters::ghz1_prime7
complex(8), public ghz1_prime7
Definition: mod_Parameters.F90:393
modparameters::getmass
real(8) function getmass(Part)
Definition: mod_Parameters.F90:2021
modparameters::ghwpwp2_prime3
complex(8), public ghwpwp2_prime3
Definition: mod_Parameters.F90:677
modparameters::includevprime
logical, public includevprime
Definition: mod_Parameters.F90:214
modparameters::pmzzdistribution
real(8), dimension(1:pmzzsize, 1:2), public pmzzdistribution
Definition: mod_Parameters.F90:140
modparameters::time_int
real(8) time_int
Definition: mod_Parameters.F90:59
modparameters::hbarc2xsecunit
real(8), parameter, public hbarc2xsecunit
Definition: mod_Parameters.F90:114
modparameters::c6
complex(8), parameter, public c6
Definition: mod_Parameters.F90:994
modparameters::ghwpwp3_prime7
complex(8), public ghwpwp3_prime7
Definition: mod_Parameters.F90:689
modparameters::ci
complex(8), parameter, public ci
Definition: mod_Parameters.F90:88
modparameters::writeweightedlhe
logical, public writeweightedlhe
Definition: mod_Parameters.F90:172
modparameters::crosssec2
real(8), dimension(1:nmaxchannels), public crosssec2
Definition: mod_Parameters.F90:56
modparameters::avegg
real(dp), parameter, public avegg
Definition: mod_Parameters.F90:1019
modparameters::gh2z4
complex(8), public gh2z4
Definition: mod_Parameters.F90:734
modparameters::lambda2_z2
real(8), parameter, public lambda2_z2
Definition: mod_Parameters.F90:779
modparameters::gethiggsdecaywidth
real(8) function gethiggsdecaywidth(jH)
Definition: mod_Parameters.F90:2133
modparameters::zff
real(8) function zff(id1in, id2in, h1, h2)
Definition: mod_Parameters.F90:1650
modparameters::gh2w4
complex(8), public gh2w4
Definition: mod_Parameters.F90:808
modparameters::lambda
real(8), parameter, public lambda
Definition: mod_Parameters.F90:354
modparameters::bgsgs1
complex(8), public bgsgs1
Definition: mod_Parameters.F90:975
modparameters::bgsgs8
complex(8), public bgsgs8
Definition: mod_Parameters.F90:979
modparameters::ghwpwp1_prime
complex(8), public ghwpwp1_prime
Definition: mod_Parameters.F90:667
modparameters::graviton_qq_left
complex(8), public graviton_qq_left
Definition: mod_Parameters.F90:928
modparameters::ghzzp3_prime3
complex(8), public ghzzp3_prime3
Definition: mod_Parameters.F90:538
modparameters::ghzpzp1_prime4
complex(8), public ghzpzp1_prime4
Definition: mod_Parameters.F90:568
modparameters::cr_qdn
real(8), public cr_qdn
Definition: mod_Parameters.F90:1074
modparameters::dv_a
complex(8), public dv_a
Definition: mod_Parameters.F90:868
modparameters::pol_b
real(8), public pol_b
Definition: mod_Parameters.F90:261
modparameters::brhadr_w_cs
real(8), parameter, public brhadr_w_cs
Definition: mod_Parameters.F90:325
modparameters::ghzzp1_prime3
complex(8), public ghzzp1_prime3
Definition: mod_Parameters.F90:522
modparameters::ghzpzp1_prime6
complex(8), public ghzpzp1_prime6
Definition: mod_Parameters.F90:570
modparameters::bzzp10
complex(8), public bzzp10
Definition: mod_Parameters.F90:956
modparameters::bgsgs4
complex(8), public bgsgs4
Definition: mod_Parameters.F90:978
modparameters::io_logfile
integer, parameter io_logfile
Definition: mod_Parameters.F90:1142
modparameters::c42
complex(8), parameter, public c42
Definition: mod_Parameters.F90:992
modparameters::ghz2_prime5
complex(8), public ghz2_prime5
Definition: mod_Parameters.F90:399
modparameters::gh2z1_prime2
complex(8), public gh2z1_prime2
Definition: mod_Parameters.F90:745
modparameters::pol_mass2
complex(dp) function, dimension(4) pol_mass2(p, i, out)
Definition: mod_Parameters.F90:3176
modparameters::cw_q2sq
integer, public cw_q2sq
Definition: mod_Parameters.F90:498
modparameters::ipart_sel
integer, public ipart_sel
Definition: mod_Parameters.F90:58
modparameters::lambda_w2
real(8), parameter, public lambda_w2
Definition: mod_Parameters.F90:493
modparameters::ghgsgs4
complex(8), public ghgsgs4
Definition: mod_Parameters.F90:425
modparameters::cr_neu
real(8), public cr_neu
Definition: mod_Parameters.F90:1070
modparameters::m_top_4gen
real(8), public m_top_4gen
Definition: mod_Parameters.F90:245
modparameters::ghw3_prime5
complex(8), public ghw3_prime5
Definition: mod_Parameters.F90:480
modparameters::t3dl
real(8), parameter, public t3dl
Definition: mod_Parameters.F90:1030
modparameters::wwpmode
integer, parameter, public wwpmode
Definition: mod_Parameters.F90:14
modparameters::rejecounter
integer(8), public rejecounter
Definition: mod_Parameters.F90:49
modparameters::gh2w2_prime6
complex(8), public gh2w2_prime6
Definition: mod_Parameters.F90:824
modparameters::adn_
integer, target, public adn_
Definition: mod_Parameters.F90:1107
modparameters::ghw4_prime4
complex(8), public ghw4_prime4
Definition: mod_Parameters.F90:487
modparameters::gh2w4_prime2
complex(8), public gh2w4_prime2
Definition: mod_Parameters.F90:836
modparameters::m_bot_4gen
real(8), public m_bot_4gen
Definition: mod_Parameters.F90:244
modparameters::wm_
integer, target, public wm_
Definition: mod_Parameters.F90:1115
modparameters::kappa2_tilde_bot
complex(8), public kappa2_tilde_bot
Definition: mod_Parameters.F90:900
modparameters::alpha_dip
real(8), public alpha_dip
Definition: mod_Parameters.F90:199
modparameters::facscheme
integer, public facscheme
Definition: mod_Parameters.F90:20
modparameters::zero
real(8), parameter, public zero
Definition: mod_Parameters.F90:85
modparameters::ghzpzp4_prime
complex(8), public ghzpzp4_prime
Definition: mod_Parameters.F90:589
modparameters::ghzpzp4_prime5
complex(8), public ghzpzp4_prime5
Definition: mod_Parameters.F90:593
modparameters::br_z_ee
real(8), parameter, public br_z_ee
Definition: mod_Parameters.F90:301
modparameters::chm_
integer, target, public chm_
Definition: mod_Parameters.F90:1086
modparameters::etalepcut
real(8), public etalepcut
Definition: mod_Parameters.F90:217
modparameters::kappa2_top
complex(8), public kappa2_top
Definition: mod_Parameters.F90:897
modparameters::m_v_ps
real(8), public m_v_ps
Definition: mod_Parameters.F90:78
modparameters::up_
integer, target, public up_
Definition: mod_Parameters.F90:1084
modparameters::lambda_z3
real(8), parameter, public lambda_z3
Definition: mod_Parameters.F90:429
modparameters::bot_
integer, target, public bot_
Definition: mod_Parameters.F90:1089
modparameters::cl_lep
real(8), public cl_lep
Definition: mod_Parameters.F90:1069
modparameters::krenfacscheme_minptj
integer, parameter, public krenfacscheme_minptj
Definition: mod_Parameters.F90:43
modparameters::gwsq
real(8), public gwsq
Definition: mod_Parameters.F90:250
modparameters::gh2w4_prime5
complex(8), public gh2w4_prime5
Definition: mod_Parameters.F90:839
modparameters::ghzzp1
complex(8), public ghzzp1
Definition: mod_Parameters.F90:514
modparameters::milimeter
real(8), parameter, public milimeter
Definition: mod_Parameters.F90:99
modparameters::krenfacscheme_mjj
integer, parameter, public krenfacscheme_mjj
Definition: mod_Parameters.F90:37
modparameters::ghwwp4_prime2
complex(8), public ghwwp4_prime2
Definition: mod_Parameters.F90:654
modparameters::brlept_z_tt
real(8), parameter, public brlept_z_tt
Definition: mod_Parameters.F90:313
modparameters::ghzzp1_prime6
complex(8), public ghzzp1_prime6
Definition: mod_Parameters.F90:525
modparameters::massfrun
real(8) function massfrun(mf, scale)
Definition: mod_Parameters.F90:2904
modparameters::half
real(8), parameter, public half
Definition: mod_Parameters.F90:84
modparameters::br_w_mn
real(8), parameter, public br_w_mn
Definition: mod_Parameters.F90:306
modparameters::ghzzp4
complex(8), public ghzzp4
Definition: mod_Parameters.F90:517
modparameters::cr_lep
real(8), public cr_lep
Definition: mod_Parameters.F90:1068
modparameters::io_lheoutfile
integer, parameter io_lheoutfile
Definition: mod_Parameters.F90:1139
modparameters::lhapdf_data_path
character(len=500) lhapdf_data_path
Definition: mod_Parameters.F90:67
modparameters::lambda_w30
real(8), public lambda_w30
Definition: mod_Parameters.F90:510
modparameters::lheprodfile
character(len=500) lheprodfile
Definition: mod_Parameters.F90:64
modparameters::gh2w2_prime4
complex(8), public gh2w2_prime4
Definition: mod_Parameters.F90:822
modparameters::ghz3_prime3
complex(8), public ghz3_prime3
Definition: mod_Parameters.F90:405
modparameters::lambda_w22
real(8), public lambda_w22
Definition: mod_Parameters.F90:505
modparameters::br_z_mm
real(8), parameter, public br_z_mm
Definition: mod_Parameters.F90:302
modparameters::v0
complex(dp) function, dimension(4) v0(p, i)
Definition: mod_Parameters.F90:3263
modparameters::tol
real(8), parameter, public tol
Definition: mod_Parameters.F90:12
modparameters::ghwwp1_prime5
complex(8), public ghwwp1_prime5
Definition: mod_Parameters.F90:633
modparameters::a3
complex(8), public a3
Definition: mod_Parameters.F90:925
modparameters::kev
real(8), parameter, public kev
Definition: mod_Parameters.F90:96
modparameters::ghw3_prime6
complex(8), public ghw3_prime6
Definition: mod_Parameters.F90:481
modparameters::ezp_chm_right
complex(8), public ezp_chm_right
Definition: mod_Parameters.F90:609
modparameters::ga_reso
real(8), public ga_reso
Definition: mod_Parameters.F90:231
modparameters::h_dk
logical, public h_dk
Definition: mod_Parameters.F90:180
modparameters::ghwwp1_prime
complex(8), public ghwwp1_prime
Definition: mod_Parameters.F90:629
modparameters::second
real(8), parameter, public second
Definition: mod_Parameters.F90:105
modparameters::gh2g4_4gen
complex(8), public gh2g4_4gen
Definition: mod_Parameters.F90:728
modparameters::gh2z4_prime5
complex(8), public gh2z4_prime5
Definition: mod_Parameters.F90:772
modparameters::isaphoton
logical function isaphoton(DKMode)
Definition: mod_Parameters.F90:2313
modparameters::io_lheoutfile3
integer, parameter io_lheoutfile3
Definition: mod_Parameters.F90:1145
modparameters::islhedowntypequark
logical function islhedowntypequark(PartType)
Definition: mod_Parameters.F90:2336
modparameters::br_z_ll_counter
integer, public br_z_ll_counter
Definition: mod_Parameters.F90:120
modparameters::gra_
integer, target, public gra_
Definition: mod_Parameters.F90:1104
modparameters::lambda_z11
real(8), public lambda_z11
Definition: mod_Parameters.F90:438
modparameters::bzzp1
complex(8), public bzzp1
Definition: mod_Parameters.F90:947
modparameters::zprime_qq_right
complex(8), public zprime_qq_right
Definition: mod_Parameters.F90:913
modparameters::br_z_ss
real(8), parameter, public br_z_ss
Definition: mod_Parameters.F90:294
modparameters::zzmode
integer, parameter, public zzmode
Definition: mod_Parameters.F90:13
modparameters::m_top
real(8), public m_top
Definition: mod_Parameters.F90:224
modparameters::ghzpgs3
complex(8), public ghzpgs3
Definition: mod_Parameters.F90:555
modparameters::nf
real(dp), parameter, public nf
Definition: mod_Parameters.F90:1014
modparameters::br_w_ud
real(8), parameter, public br_w_ud
Definition: mod_Parameters.F90:308
modparameters::gh2g4
complex(8), public gh2g4
Definition: mod_Parameters.F90:723
modparameters::m_wprime
real(8), public m_wprime
Definition: mod_Parameters.F90:712
modparameters::b9
complex(8), public b9
Definition: mod_Parameters.F90:944
modparameters::lambda2_w32
real(8), public lambda2_w32
Definition: mod_Parameters.F90:857
modparameters::ghz3_prime4
complex(8), public ghz3_prime4
Definition: mod_Parameters.F90:406
modparameters::lambda_w4
real(8), parameter, public lambda_w4
Definition: mod_Parameters.F90:495
modparameters::isalhelightquark
logical function isalhelightquark(PartType)
Definition: mod_Parameters.F90:2388
modparameters::convertlhefile
logical, public convertlhefile
Definition: mod_Parameters.F90:28
modparameters::ar_qdn
real(8), public ar_qdn
Definition: mod_Parameters.F90:1064
modparameters::lambda2_w2
real(8), parameter, public lambda2_w2
Definition: mod_Parameters.F90:844
modparameters::dzawpwm
complex(8), public dzawpwm
Definition: mod_Parameters.F90:875
modparameters::bzpzp9
complex(8), public bzpzp9
Definition: mod_Parameters.F90:966
modparameters::gh2gsgs3
complex(8), public gh2gsgs3
Definition: mod_Parameters.F90:740
modparameters::c41
complex(8), parameter, public c41
Definition: mod_Parameters.F90:991
modparameters::ckmbare
real(8) function ckmbare(id1in, id2in)
Definition: mod_Parameters.F90:1563
modparameters::ghw3_prime4
complex(8), public ghw3_prime4
Definition: mod_Parameters.F90:479
modparameters::ghwwp3_prime
complex(8), public ghwwp3_prime
Definition: mod_Parameters.F90:645
modparameters::requevents
integer(8), dimension(-6:+6,-6:+6), public requevents
Definition: mod_Parameters.F90:52
modparameters::ghz3
complex(8), public ghz3
Definition: mod_Parameters.F90:383
modparameters::coupledvertexisdiagonal
logical function coupledvertexisdiagonal(id, hel)
Definition: mod_Parameters.F90:2527
modparameters::krenfacscheme_mjhstar
integer, parameter, public krenfacscheme_mjhstar
Definition: mod_Parameters.F90:39
modparameters::ghwpwp2
complex(8), public ghwpwp2
Definition: mod_Parameters.F90:662
modparameters::gh2z3_prime7
complex(8), public gh2z3_prime7
Definition: mod_Parameters.F90:766
modparameters::gh2w3_prime6
complex(8), public gh2w3_prime6
Definition: mod_Parameters.F90:832
modparameters::gh2zgs1_prime2
complex(8), public gh2zgs1_prime2
Definition: mod_Parameters.F90:776
modparameters::isalepton
logical function isalepton(PartType)
Definition: mod_Parameters.F90:2409
modparameters::bzgs2
complex(8), public bzgs2
Definition: mod_Parameters.F90:970
modparameters::vckm_tb
real(8), public vckm_tb
Definition: mod_Parameters.F90:280
modparameters::computeqcdvariables
subroutine computeqcdvariables()
Definition: mod_Parameters.F90:2841
modparameters::krenfacscheme_mj_mhstar
integer, parameter, public krenfacscheme_mj_mhstar
Definition: mod_Parameters.F90:40
modparameters::ghwpwp3_prime3
complex(8), public ghwpwp3_prime3
Definition: mod_Parameters.F90:685
modparameters::pisq
real(dp), parameter, public pisq
Definition: mod_Parameters.F90:82
modparameters::ghzzp2_prime6
complex(8), public ghzzp2_prime6
Definition: mod_Parameters.F90:533
modparameters::writefailedevents
integer, public writefailedevents
Definition: mod_Parameters.F90:30
modparameters::computeewvariables
subroutine computeewvariables()
Definition: mod_Parameters.F90:2785
modparameters::ghg2_4gen
complex(8), public ghg2_4gen
Definition: mod_Parameters.F90:376
modparameters::t3nr
real(8), parameter, public t3nr
Definition: mod_Parameters.F90:1027
modparameters::ghwpwp1_prime5
complex(8), public ghwpwp1_prime5
Definition: mod_Parameters.F90:671
modparameters::mone
real(8), parameter, public mone
Definition: mod_Parameters.F90:83
modparameters::dm_a
complex(8), public dm_a
Definition: mod_Parameters.F90:870
modparameters::attobarn
real(8), parameter, public attobarn
Definition: mod_Parameters.F90:104
modparameters::isuptypelightquark
logical function isuptypelightquark(PartType)
Definition: mod_Parameters.F90:2356
modparameters::brhadr_z_bb
real(8), parameter, public brhadr_z_bb
Definition: mod_Parameters.F90:323
modparameters::aff
real(8) function aff(id1in, id2in, h1, h2)
Definition: mod_Parameters.F90:1659
modparameters::bzgs8
complex(8), public bzgs8
Definition: mod_Parameters.F90:973
modparameters::debugvar
real(8), dimension(0:10), public debugvar
Definition: mod_Parameters.F90:1149
modparameters::ghz1_prime4
complex(8), public ghz1_prime4
Definition: mod_Parameters.F90:390
modparameters::gsgmode
integer, parameter, public gsgmode
Definition: mod_Parameters.F90:13
modparameters::vegasnc1
integer, public vegasnc1
Definition: mod_Parameters.F90:18
modparameters::unweighted
logical, public unweighted
Definition: mod_Parameters.F90:28
modparameters::gh2w2_prime2
complex(8), public gh2w2_prime2
Definition: mod_Parameters.F90:820
modparameters::ghwwp3_prime5
complex(8), public ghwwp3_prime5
Definition: mod_Parameters.F90:649
modparameters::nrenfacschemes
integer, parameter, public nrenfacschemes
Definition: mod_Parameters.F90:45
modparameters::gh2z2_prime7
complex(8), public gh2z2_prime7
Definition: mod_Parameters.F90:758
modparameters::scale_alpha_z_ll
real(8), public scale_alpha_z_ll
Definition: mod_Parameters.F90:329
modparameters::crosssectionwithweightserrorsquared
real(8), public crosssectionwithweightserrorsquared
Definition: mod_Parameters.F90:147
modparameters::lambda2_zgs1
real(8), parameter, public lambda2_zgs1
Definition: mod_Parameters.F90:782
modparameters::ghz1_prime2
complex(8), public ghz1_prime2
Definition: mod_Parameters.F90:388
modparameters::dp_z
complex(8), public dp_z
Definition: mod_Parameters.F90:872
modparameters::ghzzp4_prime5
complex(8), public ghzzp4_prime5
Definition: mod_Parameters.F90:548
modparameters::bzpzp5
complex(8), public bzpzp5
Definition: mod_Parameters.F90:962
modparameters::gh2g3
complex(8), public gh2g3
Definition: mod_Parameters.F90:722
modparameters::bl
real(8), public bl
Definition: mod_Parameters.F90:1066
modparameters::lambda2_z12
real(8), public lambda2_z12
Definition: mod_Parameters.F90:792
modparameters::ewp_mu_left
complex(8), public ewp_mu_left
Definition: mod_Parameters.F90:701
modparameters::ghz4_prime5
complex(8), public ghz4_prime5
Definition: mod_Parameters.F90:415
modparameters::offshellv1
logical, public offshellv1
Definition: mod_Parameters.F90:28
modparameters::ghw4_prime3
complex(8), public ghw4_prime3
Definition: mod_Parameters.F90:486
modparameters::ezp_mu_right
complex(8), public ezp_mu_right
Definition: mod_Parameters.F90:601
modparameters::ghzzp2_prime
complex(8), public ghzzp2_prime
Definition: mod_Parameters.F90:528
modparameters::ezp_chm_left
complex(8), public ezp_chm_left
Definition: mod_Parameters.F90:608
modparameters::vegasnc2_default
integer, public vegasnc2_default
Definition: mod_Parameters.F90:22
modparameters::getparticle
character(len=3) function getparticle(Part)
Definition: mod_Parameters.F90:2146
modparameters::pmzzmaxindex
integer, public pmzzmaxindex
Definition: mod_Parameters.F90:141
modparameters::lambda_zgs1
real(8), parameter, public lambda_zgs1
Definition: mod_Parameters.F90:431
modparameters::cz_q2sq
integer, public cz_q2sq
Definition: mod_Parameters.F90:435
modparameters::vckm_ts
real(8), public vckm_ts
Definition: mod_Parameters.F90:281
modparameters::gh2z3_prime4
complex(8), public gh2z3_prime4
Definition: mod_Parameters.F90:763
modparameters::pchannel
integer, public pchannel
Definition: mod_Parameters.F90:17
modparameters::csmaxfile
character(len=500) csmaxfile
Definition: mod_Parameters.F90:62
modparameters::ghwwp3_prime2
complex(8), public ghwwp3_prime2
Definition: mod_Parameters.F90:646
modparameters::ewp_ta_left
complex(8), public ewp_ta_left
Definition: mod_Parameters.F90:703
modparameters::includeinterference
logical, public includeinterference
Definition: mod_Parameters.F90:77
modparameters::ghw1_prime6
complex(8), public ghw1_prime6
Definition: mod_Parameters.F90:465
modparameters::ca
real(dp), parameter, public ca
Definition: mod_Parameters.F90:1016
modparameters::gh2z3_prime3
complex(8), public gh2z3_prime3
Definition: mod_Parameters.F90:762
modparameters::ghzzp2
complex(8), public ghzzp2
Definition: mod_Parameters.F90:515
modparameters::mincs
real(8), public mincs
Definition: mod_Parameters.F90:119
modparameters::pmzzminindex
integer, public pmzzminindex
Definition: mod_Parameters.F90:141
modparameters::m_el
real(8), public m_el
Definition: mod_Parameters.F90:240
modparameters::ghwpwp2_prime
complex(8), public ghwpwp2_prime
Definition: mod_Parameters.F90:675
modparameters::requestnleptons
integer, dimension(1:2), public requestnleptons
Definition: mod_Parameters.F90:24
modparameters::generate_bis
logical, parameter, public generate_bis
Definition: mod_Parameters.F90:932
modparameters::kappa
complex(8), public kappa
Definition: mod_Parameters.F90:882
modparameters::jetsoppositeeta
logical, public jetsoppositeeta
Definition: mod_Parameters.F90:219
modparameters::ghzzp3_prime7
complex(8), public ghzzp3_prime7
Definition: mod_Parameters.F90:542
modparameters::kappa_4gen_top
complex(8), public kappa_4gen_top
Definition: mod_Parameters.F90:891
modparameters::requevents2
integer(8), dimension(1:nmaxchannels), public requevents2
Definition: mod_Parameters.F90:54
modparameters::lambda2_z40
real(8), public lambda2_z40
Definition: mod_Parameters.F90:799
modparameters::ghwwp1
complex(8), public ghwwp1
Definition: mod_Parameters.F90:623
modparameters::gh2z3
complex(8), public gh2z3
Definition: mod_Parameters.F90:733
modparameters::ahz1
complex(8), parameter, public ahz1
Definition: mod_Parameters.F90:365
modparameters::ga_reso2
real(8), public ga_reso2
Definition: mod_Parameters.F90:233
modparameters::isalhequark
logical function isalhequark(PartType)
Definition: mod_Parameters.F90:2375
modparameters::ghzzp1_prime4
complex(8), public ghzzp1_prime4
Definition: mod_Parameters.F90:523
modparameters::lambda2_w12
real(8), public lambda2_w12
Definition: mod_Parameters.F90:855
modparameters::printpmzzintervals
integer, public printpmzzintervals
Definition: mod_Parameters.F90:143
modparameters::scale_alpha_w_cs
real(8), public scale_alpha_w_cs
Definition: mod_Parameters.F90:333
modparameters::lambda_z32
real(8), public lambda_z32
Definition: mod_Parameters.F90:444
modparameters::ghw2_prime5
complex(8), public ghw2_prime5
Definition: mod_Parameters.F90:472
modparameters::alertcounter
integer(8), public alertcounter
Definition: mod_Parameters.F90:51
modparameters::ev
real(8), parameter, public ev
Definition: mod_Parameters.F90:97
modparameters::islheuptypelightquark
logical function islheuptypelightquark(PartType)
Definition: mod_Parameters.F90:2362
modparameters::brlept_w_tn
real(8), parameter, public brlept_w_tn
Definition: mod_Parameters.F90:317
modparameters::lambda2_w11
real(8), public lambda2_w11
Definition: mod_Parameters.F90:851
modparameters::xn
real(dp), parameter, public xn
Definition: mod_Parameters.F90:1015
modparameters::lambda_w10
real(8), public lambda_w10
Definition: mod_Parameters.F90:508
modparameters::gh2z2_prime5
complex(8), public gh2z2_prime5
Definition: mod_Parameters.F90:756
modparameters::lambda2_w41
real(8), public lambda2_w41
Definition: mod_Parameters.F90:854
modparameters::ghw2
complex(8), public ghw2
Definition: mod_Parameters.F90:455
modparameters::crosssectionwithweights
real(8), public crosssectionwithweights
Definition: mod_Parameters.F90:147
modparameters::ghw2_prime4
complex(8), public ghw2_prime4
Definition: mod_Parameters.F90:471
modparameters::ar_neu
real(8), public ar_neu
Definition: mod_Parameters.F90:1060
modparameters::ghzpgs4
complex(8), public ghzpgs4
Definition: mod_Parameters.F90:556
modparameters::tam_
integer, target, public tam_
Definition: mod_Parameters.F90:1114
modparameters::br_w_ll_counter
integer, public br_w_ll_counter
Definition: mod_Parameters.F90:124
modparameters::ghzpzp2_prime3
complex(8), public ghzpzp2_prime3
Definition: mod_Parameters.F90:575
modparameters::scale_alpha_w_ud
real(8), public scale_alpha_w_ud
Definition: mod_Parameters.F90:332
modparameters::couplwffsq
real(dp), public couplwffsq
Definition: mod_Parameters.F90:1078
modparameters::c5
complex(8), parameter, public c5
Definition: mod_Parameters.F90:993
modparameters::isaneutrino
logical function isaneutrino(PartType)
Definition: mod_Parameters.F90:2396
modparameters::bzzp5
complex(8), public bzzp5
Definition: mod_Parameters.F90:951
modparameters::xw
real(8), public xw
Definition: mod_Parameters.F90:253
modparameters::m_mu
real(8), public m_mu
Definition: mod_Parameters.F90:241
modparameters::zprime_qq_left
complex(8), public zprime_qq_left
Definition: mod_Parameters.F90:912
modparameters::brlept_z_ee
real(8), parameter, public brlept_z_ee
Definition: mod_Parameters.F90:311
modparameters::isalhelepton
logical function isalhelepton(PartType)
Definition: mod_Parameters.F90:2415
modparameters::gh2z3_prime6
complex(8), public gh2z3_prime6
Definition: mod_Parameters.F90:765
modparameters::gh2g2
complex(8), public gh2g2
Definition: mod_Parameters.F90:721
modparameters::bzpgs4
complex(8), public bzpgs4
Definition: mod_Parameters.F90:984
modparameters::picobarn
real(8), parameter, public picobarn
Definition: mod_Parameters.F90:102
modparameters::ghwpwp2_prime7
complex(8), public ghwpwp2_prime7
Definition: mod_Parameters.F90:681
modparameters::ghw1_prime3
complex(8), public ghw1_prime3
Definition: mod_Parameters.F90:462
modparameters::al_qdn
real(8), public al_qdn
Definition: mod_Parameters.F90:1065
modparameters::ubar0
complex(dp) function, dimension(4) ubar0(p, i)
Definition: mod_Parameters.F90:3193
modparameters::ghwpwp4_prime2
complex(8), public ghwpwp4_prime2
Definition: mod_Parameters.F90:692
modparameters::gh2g2_4gen
complex(8), public gh2g2_4gen
Definition: mod_Parameters.F90:726
modparameters::ghwwp2_prime2
complex(8), public ghwwp2_prime2
Definition: mod_Parameters.F90:638
modparameters::aup_
integer, target, public aup_
Definition: mod_Parameters.F90:1106
modparameters::io_lheoutfile2
integer, parameter io_lheoutfile2
Definition: mod_Parameters.F90:1144
modparameters::lambda2_z20
real(8), public lambda2_z20
Definition: mod_Parameters.F90:797
modparameters::br_w_hadr
real(8), parameter, public br_w_hadr
Definition: mod_Parameters.F90:297
modparameters::gh2w3_prime5
complex(8), public gh2w3_prime5
Definition: mod_Parameters.F90:831
modparameters::lambda_w11
real(8), public lambda_w11
Definition: mod_Parameters.F90:500
modparameters::ghzpzp3
complex(8), public ghzpzp3
Definition: mod_Parameters.F90:561
modparameters::daawpwm
complex(8), public daawpwm
Definition: mod_Parameters.F90:874
modparameters::ghz1
complex(8), public ghz1
Definition: mod_Parameters.F90:381
modparameters::bzzp3
complex(8), public bzzp3
Definition: mod_Parameters.F90:949
modparameters::daz_qdn_right
real(8), public daz_qdn_right
Definition: mod_Parameters.F90:1006
modparameters::kappa2_tilde_4gen_bot
complex(8), public kappa2_tilde_4gen_bot
Definition: mod_Parameters.F90:905
modparameters::ghzpzp4_prime3
complex(8), public ghzpzp4_prime3
Definition: mod_Parameters.F90:591
modparameters::cl_qdn
real(8), public cl_qdn
Definition: mod_Parameters.F90:1075
modparameters::m_w_ps
real(8), public m_w_ps
Definition: mod_Parameters.F90:78
modparameters::setdefaultckm
subroutine setdefaultckm()
Definition: mod_Parameters.F90:2774
modparameters::a4
complex(8), public a4
Definition: mod_Parameters.F90:926
modparameters::ga_zprime
real(8), public ga_zprime
Definition: mod_Parameters.F90:620
modparameters::isawdecay
logical function isawdecay(DKMode)
Definition: mod_Parameters.F90:2278
modparameters::gev
real(8), parameter, public gev
Definition: mod_Parameters.F90:93
modparameters::daz_qup_right
real(8), public daz_qup_right
Definition: mod_Parameters.F90:1004
modparameters::collider
integer, public collider
Definition: mod_Parameters.F90:17
modparameters::percent
real(8), parameter, public percent
Definition: mod_Parameters.F90:90
modparameters::ghwpwp3_prime2
complex(8), public ghwpwp3_prime2
Definition: mod_Parameters.F90:684
modparameters::nue_
integer, target, public nue_
Definition: mod_Parameters.F90:1097
modparameters::bzgs3
complex(8), public bzgs3
Definition: mod_Parameters.F90:971
modparameters::ar_lep
real(8), public ar_lep
Definition: mod_Parameters.F90:1058
modparameters::ghw1_prime4
complex(8), public ghw1_prime4
Definition: mod_Parameters.F90:463
modparameters::ghzpzp4_prime6
complex(8), public ghzpzp4_prime6
Definition: mod_Parameters.F90:594
modparameters::lambda_z42
real(8), public lambda_z42
Definition: mod_Parameters.F90:445
modparameters::accepcounter_part2
integer(8), dimension(1:nmaxchannels), public accepcounter_part2
Definition: mod_Parameters.F90:55
modparameters::mininputmhstar
real(8), public mininputmhstar
Definition: mod_Parameters.F90:135
modparameters::milisecond
real(8), parameter, public milisecond
Definition: mod_Parameters.F90:106
modparameters::b7
complex(8), public b7
Definition: mod_Parameters.F90:942
modparameters::lambda2_w30
real(8), public lambda2_w30
Definition: mod_Parameters.F90:861
modparameters::rjet
real(8), public rjet
Definition: mod_Parameters.F90:208
modparameters::ewp_up_right
complex(8), public ewp_up_right
Definition: mod_Parameters.F90:706
modparameters::sethiggsmass
subroutine sethiggsmass(jH, mass)
Definition: mod_Parameters.F90:2099
modparameters::setmass
subroutine setmass(mass, ipart)
Definition: mod_Parameters.F90:1962
modparameters::bzzp8
complex(8), public bzzp8
Definition: mod_Parameters.F90:954
modparameters::ezp_str_right
complex(8), public ezp_str_right
Definition: mod_Parameters.F90:613
modparameters::zzpmode
integer, parameter, public zzpmode
Definition: mod_Parameters.F90:15
modparameters::mjjcut
real(8), public mjjcut
Definition: mod_Parameters.F90:209
modparameters::br_w_ll
real(8), parameter, public br_w_ll
Definition: mod_Parameters.F90:296
modparameters::czero
complex(8), parameter, public czero
Definition: mod_Parameters.F90:86
modparameters::ghzzp2_prime3
complex(8), public ghzzp2_prime3
Definition: mod_Parameters.F90:530
modparameters::t3ul
real(8), parameter, public t3ul
Definition: mod_Parameters.F90:1028
modparameters::isuptypequark
logical function isuptypequark(PartType)
Definition: mod_Parameters.F90:2343
modparameters::gszmode
integer, parameter, public gszmode
Definition: mod_Parameters.F90:13
modparameters::gh2w4_prime6
complex(8), public gh2w4_prime6
Definition: mod_Parameters.F90:840
modparameters::pdfup_
integer, parameter, public pdfup_
Definition: mod_Parameters.F90:1126
modparameters::ghzpzp2_prime2
complex(8), public ghzpzp2_prime2
Definition: mod_Parameters.F90:574
modparameters::ghzpzp1_prime
complex(8), public ghzpzp1_prime
Definition: mod_Parameters.F90:565
modparameters::couplzffsq
real(dp), public couplzffsq
Definition: mod_Parameters.F90:1079
modparameters::ezp_up_right
complex(8), public ezp_up_right
Definition: mod_Parameters.F90:607
modparameters::lambda2_z31
real(8), public lambda2_z31
Definition: mod_Parameters.F90:790
modparameters::ghzzp3
complex(8), public ghzzp3
Definition: mod_Parameters.F90:516
modparameters::aveqg
real(dp), parameter, public aveqg
Definition: mod_Parameters.F90:1020
modparameters::cone
complex(8), parameter, public cone
Definition: mod_Parameters.F90:87
modparameters::gh2w2_prime
complex(8), public gh2w2_prime
Definition: mod_Parameters.F90:819
modparameters::ptlepcut
real(8), public ptlepcut
Definition: mod_Parameters.F90:216
modparameters::renscheme
integer, public renscheme
Definition: mod_Parameters.F90:20
modparameters::lambda_z30
real(8), public lambda_z30
Definition: mod_Parameters.F90:448
modparameters::brlept_z_nn
real(8), parameter, public brlept_z_nn
Definition: mod_Parameters.F90:314
modparameters::bzpzp4
complex(8), public bzpzp4
Definition: mod_Parameters.F90:961
modparameters::ghzgs1_prime2
complex(8), public ghzgs1_prime2
Definition: mod_Parameters.F90:419
modparameters::ghwpwp2_prime4
complex(8), public ghwpwp2_prime4
Definition: mod_Parameters.F90:678
modparameters::c3
complex(8), parameter, public c3
Definition: mod_Parameters.F90:990
modparameters::su2flip
integer function su2flip(Part)
Definition: mod_Parameters.F90:2603
modparameters::vpffcoupling_pdg
complex(8) function vpffcoupling_pdg(pdgid, hel, useWp)
Definition: mod_Parameters.F90:1498
modparameters::elp_
integer, target, public elp_
Definition: mod_Parameters.F90:1090
modparameters::ghzpzp4_prime2
complex(8), public ghzpzp4_prime2
Definition: mod_Parameters.F90:590
modparameters::kappa_tilde
complex(8), public kappa_tilde
Definition: mod_Parameters.F90:883
modparameters::vegasnc1_default
integer, public vegasnc1_default
Definition: mod_Parameters.F90:22
modparameters::pdfaup_
integer, parameter, public pdfaup_
Definition: mod_Parameters.F90:1132
modparameters::ghwwp1_prime3
complex(8), public ghwwp1_prime3
Definition: mod_Parameters.F90:631
modparameters::nmaxseeds
integer, parameter, public nmaxseeds
Definition: mod_Parameters.F90:159
modparameters::converttopartindex
integer function converttopartindex(Part)
Definition: mod_Parameters.F90:1886
modparameters::findcrosssectionwithweights
logical, public findcrosssectionwithweights
Definition: mod_Parameters.F90:145
modparameters::daz_qup_left
real(8), public daz_qup_left
Definition: mod_Parameters.F90:1003
modparameters::lhapdfstring
character(len=100) lhapdfstring
Definition: mod_Parameters.F90:66
modparameters::readlhefile
logical, public readlhefile
Definition: mod_Parameters.F90:28
modparameters::gh2z2_prime3
complex(8), public gh2z2_prime3
Definition: mod_Parameters.F90:754
modparameters::ghzzp3_prime5
complex(8), public ghzzp3_prime5
Definition: mod_Parameters.F90:540
modparameters::kappa2_tilde_top
complex(8), public kappa2_tilde_top
Definition: mod_Parameters.F90:898
modparameters::qdl
real(8), parameter, public qdl
Definition: mod_Parameters.F90:1038
modparameters::kappa_tilde_4gen_top
complex(8), public kappa_tilde_4gen_top
Definition: mod_Parameters.F90:892
modparameters::ctauunit
real(8), parameter, public ctauunit
Definition: mod_Parameters.F90:100
modparameters::datafile
character(len=500) datafile
Definition: mod_Parameters.F90:61
modparameters::ghw2_prime2
complex(8), public ghw2_prime2
Definition: mod_Parameters.F90:469
modparameters::br_w_ud_counter
integer, public br_w_ud_counter
Definition: mod_Parameters.F90:125
modparameters::ghzpgs2
complex(8), public ghzpgs2
Definition: mod_Parameters.F90:554
modparameters::chargeflip
integer function chargeflip(Part)
Definition: mod_Parameters.F90:2689
modparameters::ghz4_prime7
complex(8), public ghz4_prime7
Definition: mod_Parameters.F90:417
modparameters::ghwpwp4_prime5
complex(8), public ghwpwp4_prime5
Definition: mod_Parameters.F90:695
modparameters::lambda_w20
real(8), public lambda_w20
Definition: mod_Parameters.F90:509
modparameters::c2w_q12sq
integer, public c2w_q12sq
Definition: mod_Parameters.F90:850
modparameters::ga_z
real(8), public ga_z
Definition: mod_Parameters.F90:227
modparameters::m4l_minmax
real(8), dimension(1:2), public m4l_minmax
Definition: mod_Parameters.F90:210
modparameters::ghz3_prime7
complex(8), public ghz3_prime7
Definition: mod_Parameters.F90:409
modparameters::br_z_hadr
real(8), parameter, public br_z_hadr
Definition: mod_Parameters.F90:289
modparameters::zpzmode
integer, parameter, public zpzmode
Definition: mod_Parameters.F90:15
modparameters::gh2z4_prime7
complex(8), public gh2z4_prime7
Definition: mod_Parameters.F90:774
modparameters::bzzp7
complex(8), public bzzp7
Definition: mod_Parameters.F90:953
modparameters::m_w
real(8), public m_w
Definition: mod_Parameters.F90:228
modparameters::atop_
integer, target, public atop_
Definition: mod_Parameters.F90:1110
modparameters::crosssecmax2
real(8), dimension(1:nmaxchannels), public crosssecmax2
Definition: mod_Parameters.F90:57
modparameters::lambda2_w22
real(8), public lambda2_w22
Definition: mod_Parameters.F90:856
modparameters::useunformattedread
logical, public useunformattedread
Definition: mod_Parameters.F90:176
modparameters::bzpzp7
complex(8), public bzpzp7
Definition: mod_Parameters.F90:964
modparameters::br_z_tt
real(8), parameter, public br_z_tt
Definition: mod_Parameters.F90:303
modparameters::qnl
real(8), parameter, public qnl
Definition: mod_Parameters.F90:1034
modparameters::scale_alpha_z_tt
real(8), public scale_alpha_z_tt
Definition: mod_Parameters.F90:330
modparameters::zmass_pdf
real(8), public zmass_pdf
Definition: mod_Parameters.F90:266
modparameters::m_v
real(8), public m_v
Definition: mod_Parameters.F90:78
modparameters::ghw4_prime6
complex(8), public ghw4_prime6
Definition: mod_Parameters.F90:489
modparameters::ghwpwp2_prime6
complex(8), public ghwpwp2_prime6
Definition: mod_Parameters.F90:680
modparameters::ghzzp2_prime7
complex(8), public ghzzp2_prime7
Definition: mod_Parameters.F90:534
modparameters::requestossf
integer, dimension(1:2), public requestossf
Definition: mod_Parameters.F90:27
modparameters::readcsmax
logical, public readcsmax
Definition: mod_Parameters.F90:29
modparameters::vh_pc
character(len=2), public vh_pc
Definition: mod_Parameters.F90:184
modparameters::spinavg
real(8), parameter, public spinavg
Definition: mod_Parameters.F90:91
modparameters::t3nl
real(8), parameter, public t3nl
Definition: mod_Parameters.F90:1026
modparameters::anum_
integer, target, public anum_
Definition: mod_Parameters.F90:1117
modparameters::ewp_chm_left
complex(8), public ewp_chm_left
Definition: mod_Parameters.F90:707
modparameters::ghwpwp4_prime7
complex(8), public ghwpwp4_prime7
Definition: mod_Parameters.F90:697
modparameters::foundhiggswidth
logical, public foundhiggswidth
Definition: mod_Parameters.F90:29
modparameters::bzzp4
complex(8), public bzzp4
Definition: mod_Parameters.F90:950
modparameters::lambda_w3
real(8), parameter, public lambda_w3
Definition: mod_Parameters.F90:494
modparameters::gh2w1
complex(8), public gh2w1
Definition: mod_Parameters.F90:805
modparameters::gh2w1_prime5
complex(8), public gh2w1_prime5
Definition: mod_Parameters.F90:815
modparameters::al_qup
real(8), public al_qup
Definition: mod_Parameters.F90:1063
modparameters::ghzpzp3_prime7
complex(8), public ghzpzp3_prime7
Definition: mod_Parameters.F90:587
modparameters::alpha_qed
real(8), public alpha_qed
Definition: mod_Parameters.F90:248
modparameters::lambda_q
real(8), parameter, public lambda_q
Definition: mod_Parameters.F90:432
modparameters::gh2w3
complex(8), public gh2w3
Definition: mod_Parameters.F90:807
modparameters::lambda2_w40
real(8), public lambda2_w40
Definition: mod_Parameters.F90:862
modparameters::isaboson
logical function isaboson(PartType)
Definition: mod_Parameters.F90:2453
modparameters::ilc_energy
real(8), parameter, public ilc_energy
Definition: mod_Parameters.F90:258
modparameters::gh2w3_prime3
complex(8), public gh2w3_prime3
Definition: mod_Parameters.F90:829
modparameters::lambda_w41
real(8), public lambda_w41
Definition: mod_Parameters.F90:503
modparameters::ewp_el_right
complex(8), public ewp_el_right
Definition: mod_Parameters.F90:700
modparameters::qnr
real(8), parameter, public qnr
Definition: mod_Parameters.F90:1035
modparameters::computeckmelements
subroutine computeckmelements(inVCKM_ud, inVCKM_us, inVCKM_cd, inVCKM_cs, inVCKM_ts, inVCKM_tb, inVCKM_ub, inVCKM_cb, inVCKM_td)
Definition: mod_Parameters.F90:2708
modparameters::bzzp6
complex(8), public bzzp6
Definition: mod_Parameters.F90:952
modparameters::kappa2_bot
complex(8), public kappa2_bot
Definition: mod_Parameters.F90:899
modparameters::ghzpzp2_prime5
complex(8), public ghzpzp2_prime5
Definition: mod_Parameters.F90:577
modparameters::lambda_w40
real(8), public lambda_w40
Definition: mod_Parameters.F90:511
modparameters::ghz1_prime3
complex(8), public ghz1_prime3
Definition: mod_Parameters.F90:389
modparameters::ghwwp3
complex(8), public ghwwp3
Definition: mod_Parameters.F90:625
modparameters::ga_v
real(8), public ga_v
Definition: mod_Parameters.F90:78
modparameters::ghzzp4_prime7
complex(8), public ghzzp4_prime7
Definition: mod_Parameters.F90:550
modparameters::gh2w4_prime7
complex(8), public gh2w4_prime7
Definition: mod_Parameters.F90:841
modparameters::ewp_up_left
complex(8), public ewp_up_left
Definition: mod_Parameters.F90:705
modparameters::anut_
integer, target, public anut_
Definition: mod_Parameters.F90:1118
modparameters::meter
real(8), parameter, public meter
Definition: mod_Parameters.F90:98
modparameters::pol_dk2mom
complex(dp) function, dimension(4) pol_dk2mom(plepton, antilepton, i, outgoing)
Definition: mod_Parameters.F90:3068
modparameters::sethiggsdecaywidth
subroutine sethiggsdecaywidth(jH, width)
Definition: mod_Parameters.F90:2110
modparameters::anue_
integer, target, public anue_
Definition: mod_Parameters.F90:1116
modparameters::br_z_inv
real(8), parameter, public br_z_inv
Definition: mod_Parameters.F90:290
modparameters::gh2w1_prime3
complex(8), public gh2w1_prime3
Definition: mod_Parameters.F90:813
modparameters::pdfchm_
integer, parameter, public pdfchm_
Definition: mod_Parameters.F90:1128
modparameters::gh2w4_prime4
complex(8), public gh2w4_prime4
Definition: mod_Parameters.F90:838
modparameters::ezp_dn_left
complex(8), public ezp_dn_left
Definition: mod_Parameters.F90:610
modparameters::dp_a
complex(8), public dp_a
Definition: mod_Parameters.F90:869
modparameters::brhadr_z_dd
real(8), parameter, public brhadr_z_dd
Definition: mod_Parameters.F90:321
modparameters::ghz2_prime6
complex(8), public ghz2_prime6
Definition: mod_Parameters.F90:400
modparameters::kappa2_4gen_top
complex(8), public kappa2_4gen_top
Definition: mod_Parameters.F90:902
modparameters::ghz4_prime
complex(8), public ghz4_prime
Definition: mod_Parameters.F90:411
modparameters::pthcut
real(8), public pthcut
Definition: mod_Parameters.F90:218
modparameters::coupltolhewp
logical function coupltolhewp(Part)
Definition: mod_Parameters.F90:2640
modparameters::gh2w3_prime2
complex(8), public gh2w3_prime2
Definition: mod_Parameters.F90:828
modparameters::br_w_tn
real(8), parameter, public br_w_tn
Definition: mod_Parameters.F90:307
modparameters::getdecaywidth
real(8) function getdecaywidth(Part)
Definition: mod_Parameters.F90:2074
modparameters::io_tmpfile
integer, parameter io_tmpfile
Definition: mod_Parameters.F90:1146
modparameters::islheuptypequark
logical function islheuptypequark(PartType)
Definition: mod_Parameters.F90:2349
modparameters::gh2w2_prime3
complex(8), public gh2w2_prime3
Definition: mod_Parameters.F90:821
modparameters::lambda_z40
real(8), public lambda_z40
Definition: mod_Parameters.F90:449
modparameters::ghwwp2_prime3
complex(8), public ghwwp2_prime3
Definition: mod_Parameters.F90:639
modparameters::ghg4_4gen
complex(8), public ghg4_4gen
Definition: mod_Parameters.F90:378
modparameters::vegasit1_default
integer, public vegasit1_default
Definition: mod_Parameters.F90:22
modparameters::gh2zgs3
complex(8), public gh2zgs3
Definition: mod_Parameters.F90:737
modparameters::jpart_sel
integer, public jpart_sel
Definition: mod_Parameters.F90:58
modparameters
Definition: mod_Parameters.F90:1
modparameters::warmup
logical, public warmup
Definition: mod_Parameters.F90:60
modparameters::ezp_nue_right
complex(8), public ezp_nue_right
Definition: mod_Parameters.F90:605
modparameters::accepcounter
integer(8), public accepcounter
Definition: mod_Parameters.F90:50
modparameters::gh2z1_prime
complex(8), public gh2z1_prime
Definition: mod_Parameters.F90:744
modparameters::reweightinterference
logical, public reweightinterference
Definition: mod_Parameters.F90:77
modparameters::countleptons
integer function countleptons(MY_IDUP)
Definition: mod_Parameters.F90:2572
modparameters::vckm_ud
real(8), public vckm_ud
Definition: mod_Parameters.F90:276
modparameters::ghz3_prime
complex(8), public ghz3_prime
Definition: mod_Parameters.F90:403
modparameters::ghwwp3_prime3
complex(8), public ghwwp3_prime3
Definition: mod_Parameters.F90:647
modparameters::c2z_q2sq
integer, public c2z_q2sq
Definition: mod_Parameters.F90:786
modparameters::scale_alpha_z_nn
real(8), public scale_alpha_z_nn
Definition: mod_Parameters.F90:331
modparameters::gh2z3_prime5
complex(8), public gh2z3_prime5
Definition: mod_Parameters.F90:764
modparameters::m_charm
real(8), public m_charm
Definition: mod_Parameters.F90:239
modparameters::wpwmode
integer, parameter, public wpwmode
Definition: mod_Parameters.F90:14
modparameters::lambda2_q
real(8), parameter, public lambda2_q
Definition: mod_Parameters.F90:783
modparameters::ghwpwp4
complex(8), public ghwpwp4
Definition: mod_Parameters.F90:664
modparameters::dfour_a
complex(8), public dfour_a
Definition: mod_Parameters.F90:877
modparameters::ghwpwp4_prime3
complex(8), public ghwpwp4_prime3
Definition: mod_Parameters.F90:693
modparameters::ghwpwp4_prime6
complex(8), public ghwpwp4_prime6
Definition: mod_Parameters.F90:696
modparameters::t3lr
real(8), parameter, public t3lr
Definition: mod_Parameters.F90:1025
modparameters::maxinputmhstar
real(8), public maxinputmhstar
Definition: mod_Parameters.F90:135
modparameters::ghz2_prime4
complex(8), public ghz2_prime4
Definition: mod_Parameters.F90:398
modparameters::ezp_ta_left
complex(8), public ezp_ta_left
Definition: mod_Parameters.F90:602
modparameters::ghzpzp4_prime7
complex(8), public ghzpzp4_prime7
Definition: mod_Parameters.F90:595
modparameters::gh2z4_prime4
complex(8), public gh2z4_prime4
Definition: mod_Parameters.F90:771
modparameters::ghw2_prime
complex(8), public ghw2_prime
Definition: mod_Parameters.F90:468
modparameters::lambda_z1
real(8), parameter, public lambda_z1
Definition: mod_Parameters.F90:427
modparameters::wppr_
integer, target, public wppr_
Definition: mod_Parameters.F90:1103
modparameters::scalefactor
real(8) function scalefactor(id1in, id2in)
Definition: mod_Parameters.F90:1509
modparameters::pdfbot_
integer, parameter, public pdfbot_
Definition: mod_Parameters.F90:1129
modparameters::ghw1_prime7
complex(8), public ghw1_prime7
Definition: mod_Parameters.F90:466
modparameters::femtobarn
real(8), parameter, public femtobarn
Definition: mod_Parameters.F90:103
modparameters::ghzzp1_prime2
complex(8), public ghzzp1_prime2
Definition: mod_Parameters.F90:521
modparameters::ghw4_prime
complex(8), public ghw4_prime
Definition: mod_Parameters.F90:484
modparameters::bzpzp10
complex(8), public bzpzp10
Definition: mod_Parameters.F90:967
modparameters::ezp_mu_left
complex(8), public ezp_mu_left
Definition: mod_Parameters.F90:600
modparameters::lhc_energy
real(8), parameter, public lhc_energy
Definition: mod_Parameters.F90:256
modparameters::ghwwp3_prime4
complex(8), public ghwwp3_prime4
Definition: mod_Parameters.F90:648
modparameters::ghzgs3
complex(8), public ghzgs3
Definition: mod_Parameters.F90:421
modparameters::vckm_cd
real(8), public vckm_cd
Definition: mod_Parameters.F90:279
modparameters::ghwwp4
complex(8), public ghwwp4
Definition: mod_Parameters.F90:626
modparameters::breitwignercutoff
real(8), public breitwignercutoff
Definition: mod_Parameters.F90:272
modparameters::lambda_z22
real(8), public lambda_z22
Definition: mod_Parameters.F90:443
modparameters::ghwwp2_prime
complex(8), public ghwwp2_prime
Definition: mod_Parameters.F90:637
modparameters::convertlhe
integer function convertlhe(Part)
Definition: mod_Parameters.F90:1798
modparameters::ezp_dn_right
complex(8), public ezp_dn_right
Definition: mod_Parameters.F90:611
modparameters::graviton_qq_right
complex(8), public graviton_qq_right
Definition: mod_Parameters.F90:929
modparameters::brhadr_z_uu
real(8), parameter, public brhadr_z_uu
Definition: mod_Parameters.F90:319
modparameters::vegasnc0
integer, public vegasnc0
Definition: mod_Parameters.F90:18
modparameters::nqflavors_pdf
integer, parameter, public nqflavors_pdf
Definition: mod_Parameters.F90:264
modparameters::gh2z1_prime6
complex(8), public gh2z1_prime6
Definition: mod_Parameters.F90:749
modparameters::cr_qup
real(8), public cr_qup
Definition: mod_Parameters.F90:1072
modparameters::bzpzp2
complex(8), public bzpzp2
Definition: mod_Parameters.F90:959
modparameters::a1
complex(8), public a1
Definition: mod_Parameters.F90:923
modparameters::isalheneutrino
logical function isalheneutrino(PartType)
Definition: mod_Parameters.F90:2402
modparameters::dv_z
complex(8), public dv_z
Definition: mod_Parameters.F90:871
modparameters::lambda2_w31
real(8), public lambda2_w31
Definition: mod_Parameters.F90:853
modparameters::ghwpwp2_prime5
complex(8), public ghwpwp2_prime5
Definition: mod_Parameters.F90:679
modparameters::dfour_z
complex(8), public dfour_z
Definition: mod_Parameters.F90:878
modparameters::scale_alpha_z_dd
real(8), public scale_alpha_z_dd
Definition: mod_Parameters.F90:328
modparameters::symmfac
real(8), parameter, public symmfac
Definition: mod_Parameters.F90:91
modparameters::qur
real(8), parameter, public qur
Definition: mod_Parameters.F90:1037
modparameters::ghzpzp1_prime2
complex(8), public ghzpzp1_prime2
Definition: mod_Parameters.F90:566
modparameters::b2
complex(8), public b2
Definition: mod_Parameters.F90:937
modparameters::ezp_nue_left
complex(8), public ezp_nue_left
Definition: mod_Parameters.F90:604
modparameters::lambda_w42
real(8), public lambda_w42
Definition: mod_Parameters.F90:507
modparameters::ghzzp3_prime2
complex(8), public ghzzp3_prime2
Definition: mod_Parameters.F90:537
modparameters::vpffcoupling
complex(8) function vpffcoupling(jhuid, hel, useWp)
Definition: mod_Parameters.F90:1370
modparameters::krenfacscheme_maxptj
integer, parameter, public krenfacscheme_maxptj
Definition: mod_Parameters.F90:42
modparameters::br
real(8), public br
Definition: mod_Parameters.F90:1067
modparameters::mhstarforphasespace
real(8), public mhstarforphasespace
Definition: mod_Parameters.F90:135
modparameters::hbarc
real(8), parameter, public hbarc
Definition: mod_Parameters.F90:110
modparameters::krenfacscheme_mj_mj_mhstar
integer, parameter, public krenfacscheme_mj_mj_mhstar
Definition: mod_Parameters.F90:36
modparameters::pdfachm_
integer, parameter, public pdfachm_
Definition: mod_Parameters.F90:1134
modparameters::ghwpwp1_prime6
complex(8), public ghwpwp1_prime6
Definition: mod_Parameters.F90:672
modparameters::gh2w2_prime7
complex(8), public gh2w2_prime7
Definition: mod_Parameters.F90:825
modparameters::lambda2
real(8), parameter, public lambda2
Definition: mod_Parameters.F90:355
modparameters::foundhiggsmass
logical, public foundhiggsmass
Definition: mod_Parameters.F90:29
modparameters::m_z
real(8), public m_z
Definition: mod_Parameters.F90:226
modparameters::lambda_z20
real(8), public lambda_z20
Definition: mod_Parameters.F90:447
modparameters::generate_as
logical, parameter, public generate_as
Definition: mod_Parameters.F90:361
modparameters::nut_
integer, target, public nut_
Definition: mod_Parameters.F90:1099
modparameters::pol_mass
complex(8) function, dimension(4) pol_mass(p, i, outgoing)
Definition: mod_Parameters.F90:3102
modparameters::wwmode
integer, parameter, public wwmode
Definition: mod_Parameters.F90:13
modparameters::dp
integer, parameter, public dp
Definition: mod_Parameters.F90:11
modparameters::b6
complex(8), public b6
Definition: mod_Parameters.F90:941
modparameters::ghz2_prime7
complex(8), public ghz2_prime7
Definition: mod_Parameters.F90:401
modparameters::ghw3
complex(8), public ghw3
Definition: mod_Parameters.F90:456
modparameters::ghzpzp2_prime4
complex(8), public ghzpzp2_prime4
Definition: mod_Parameters.F90:576
modparameters::vckm_ub
real(8), public vckm_ub
Definition: mod_Parameters.F90:283
modparameters::rejecounter_part
integer(8), dimension(-6:6,-6:6), public rejecounter_part
Definition: mod_Parameters.F90:52
modparameters::ghw4_prime5
complex(8), public ghw4_prime5
Definition: mod_Parameters.F90:488
modparameters::ghzpzp1_prime3
complex(8), public ghzpzp1_prime3
Definition: mod_Parameters.F90:567
modparameters::ghwwp4_prime6
complex(8), public ghwwp4_prime6
Definition: mod_Parameters.F90:658
modparameters::ghgsgs2
complex(8), public ghgsgs2
Definition: mod_Parameters.F90:423
modparameters::lambda2_z32
real(8), public lambda2_z32
Definition: mod_Parameters.F90:794
modparameters::m_z_ps
real(8), public m_z_ps
Definition: mod_Parameters.F90:78
modparameters::mu_ren
real(8), public mu_ren
Definition: mod_Parameters.F90:268
modparameters::mum_
integer, target, public mum_
Definition: mod_Parameters.F90:1113
modparameters::ga_w_ps
real(8), public ga_w_ps
Definition: mod_Parameters.F90:78
modparameters::b3
complex(8), public b3
Definition: mod_Parameters.F90:938
modparameters::zpr_
integer, target, public zpr_
Definition: mod_Parameters.F90:1101
modparameters::qlr
real(8), parameter, public qlr
Definition: mod_Parameters.F90:1033
modparameters::io_lheinfile
integer, parameter io_lheinfile
Definition: mod_Parameters.F90:1141
modparameters::ghzpgs1_prime2
complex(8), public ghzpgs1_prime2
Definition: mod_Parameters.F90:553
modparameters::c1
complex(8), parameter, public c1
Definition: mod_Parameters.F90:988
modparameters::ghzpzp3_prime4
complex(8), public ghzpzp3_prime4
Definition: mod_Parameters.F90:584
modparameters::ghw3_prime2
complex(8), public ghw3_prime2
Definition: mod_Parameters.F90:477
modparameters::isalightquark
logical function isalightquark(PartType)
Definition: mod_Parameters.F90:2382
modparameters::wp_
integer, target, public wp_
Definition: mod_Parameters.F90:1096
modparameters::mom_not_a_particle
real(8), dimension(1:4), parameter, public mom_not_a_particle
Definition: mod_Parameters.F90:1122
modparameters::channels_ratio_fix
real(8), public channels_ratio_fix
Definition: mod_Parameters.F90:170
modparameters::t3dr
real(8), parameter, public t3dr
Definition: mod_Parameters.F90:1031
modparameters::time_start
real(8) time_start
Definition: mod_Parameters.F90:59
modparameters::isalhejet
logical function isalhejet(PartType)
Definition: mod_Parameters.F90:2445
modparameters::gh2w3_prime
complex(8), public gh2w3_prime
Definition: mod_Parameters.F90:827
modparameters::lambda2_w20
real(8), public lambda2_w20
Definition: mod_Parameters.F90:860
modparameters::hbarc2_fbgev2
real(8), parameter, public hbarc2_fbgev2
Definition: mod_Parameters.F90:112
modparameters::planckconstant
real(8), parameter, public planckconstant
Definition: mod_Parameters.F90:108
modparameters::br_w_cs
real(8), parameter, public br_w_cs
Definition: mod_Parameters.F90:309
modparameters::m_bot
real(8), public m_bot
Definition: mod_Parameters.F90:238
modparameters::lambda2_z3
real(8), parameter, public lambda2_z3
Definition: mod_Parameters.F90:780
modparameters::gabarh
real(8), public gabarh
Definition: mod_Parameters.F90:134
modparameters::m_reso
real(8), public m_reso
Definition: mod_Parameters.F90:230
modparameters::isalheboson
logical function isalheboson(PartType)
Definition: mod_Parameters.F90:2460
modparameters::lenlhapdfstring
integer, public lenlhapdfstring
Definition: mod_Parameters.F90:68
modparameters::requestos
integer, dimension(1:2), public requestos
Definition: mod_Parameters.F90:26
modparameters::ghzzp2_prime4
complex(8), public ghzzp2_prime4
Definition: mod_Parameters.F90:531
modparameters::gh2w1_prime4
complex(8), public gh2w1_prime4
Definition: mod_Parameters.F90:814
modparameters::ghzzp3_prime6
complex(8), public ghzzp3_prime6
Definition: mod_Parameters.F90:541
modparameters::br_z_dd_counter
integer, public br_z_dd_counter
Definition: mod_Parameters.F90:123
modparameters::evalcounter
integer(8), public evalcounter
Definition: mod_Parameters.F90:48
modparameters::alphas
real(dp), public alphas
Definition: mod_Parameters.F90:269
modparameters::ghw3_prime3
complex(8), public ghw3_prime3
Definition: mod_Parameters.F90:478
modparameters::str_
integer, target, public str_
Definition: mod_Parameters.F90:1087
modparameters::ghzpzp2_prime
complex(8), public ghzpzp2_prime
Definition: mod_Parameters.F90:573
modparameters::lambda2_z41
real(8), public lambda2_z41
Definition: mod_Parameters.F90:791
modparameters::isaquark
logical function isaquark(PartType)
Definition: mod_Parameters.F90:2369
modparameters::ghwwp1_prime6
complex(8), public ghwwp1_prime6
Definition: mod_Parameters.F90:634
modparameters::vckm_us
real(8), public vckm_us
Definition: mod_Parameters.F90:277
modparameters::gh2w3_prime7
complex(8), public gh2w3_prime7
Definition: mod_Parameters.F90:833
modparameters::ewp_top_left
complex(8), public ewp_top_left
Definition: mod_Parameters.F90:709
modparameters::randomizevvdecays
logical, public randomizevvdecays
Definition: mod_Parameters.F90:174
modparameters::br_z_inv_counter
integer, public br_z_inv_counter
Definition: mod_Parameters.F90:121
modparameters::ghwwp1_prime2
complex(8), public ghwwp1_prime2
Definition: mod_Parameters.F90:630
modparameters::qll
real(8), parameter, public qll
Definition: mod_Parameters.F90:1032
modparameters::ghg4
complex(8), public ghg4
Definition: mod_Parameters.F90:373
modparameters::bzpgs1
complex(8), public bzpgs1
Definition: mod_Parameters.F90:981
modparameters::mev
real(8), parameter, public mev
Definition: mod_Parameters.F90:95
modparameters::pol_a
real(8), public pol_a
Definition: mod_Parameters.F90:260
modparameters::ghwwp2_prime6
complex(8), public ghwwp2_prime6
Definition: mod_Parameters.F90:642
modparameters::io_stdout
integer, parameter io_stdout
Definition: mod_Parameters.F90:1138
modparameters::hbarc2
real(8), parameter, public hbarc2
Definition: mod_Parameters.F90:111
modparameters::scale_alpha_z_uu
real(8), public scale_alpha_z_uu
Definition: mod_Parameters.F90:327
modparameters::ewp_mu_right
complex(8), public ewp_mu_right
Definition: mod_Parameters.F90:702
modparameters::etajetcut
real(8), public etajetcut
Definition: mod_Parameters.F90:206
modparameters::ghwpwp3_prime4
complex(8), public ghwpwp3_prime4
Definition: mod_Parameters.F90:686
modparameters::gh2z1_prime3
complex(8), public gh2z1_prime3
Definition: mod_Parameters.F90:746
modparameters::ghwwp3_prime7
complex(8), public ghwwp3_prime7
Definition: mod_Parameters.F90:651
modparameters::lambda2_z30
real(8), public lambda2_z30
Definition: mod_Parameters.F90:798
modparameters::ghw1_prime
complex(8), public ghw1_prime
Definition: mod_Parameters.F90:460
modparameters::astr_
integer, target, public astr_
Definition: mod_Parameters.F90:1109
modparameters::t3ur
real(8), parameter, public t3ur
Definition: mod_Parameters.F90:1029
modparameters::speedoflight
real(8), parameter, public speedoflight
Definition: mod_Parameters.F90:107
modparameters::maxcs
real(8), public maxcs
Definition: mod_Parameters.F90:119
modparameters::ghw4
complex(8), public ghw4
Definition: mod_Parameters.F90:457
modparameters::ghz1_prime6
complex(8), public ghz1_prime6
Definition: mod_Parameters.F90:392
modparameters::zpgmode
integer, parameter, public zpgmode
Definition: mod_Parameters.F90:16
modparameters::lambda2_z22
real(8), public lambda2_z22
Definition: mod_Parameters.F90:793
modparameters::lambda_w12
real(8), public lambda_w12
Definition: mod_Parameters.F90:504
modparameters::gh2z3_prime
complex(8), public gh2z3_prime
Definition: mod_Parameters.F90:760
modparameters::lambda_z4
real(8), parameter, public lambda_z4
Definition: mod_Parameters.F90:430
modparameters::brhadr_z_cc
real(8), parameter, public brhadr_z_cc
Definition: mod_Parameters.F90:320
modparameters::not_a_particle_
integer, parameter, public not_a_particle_
Definition: mod_Parameters.F90:1121
modparameters::bgsgs2
complex(8), public bgsgs2
Definition: mod_Parameters.F90:976
modparameters::ga_top
real(8), public ga_top
Definition: mod_Parameters.F90:225
modparameters::b4
complex(8), public b4
Definition: mod_Parameters.F90:939
modparameters::leptinevent
integer, dimension(0:8), public leptinevent
Definition: mod_Parameters.F90:127
modparameters::ggmode
integer, parameter, public ggmode
Definition: mod_Parameters.F90:13
modparameters::ghzpzp3_prime6
complex(8), public ghzpzp3_prime6
Definition: mod_Parameters.F90:586
modparameters::lambda_w32
real(8), public lambda_w32
Definition: mod_Parameters.F90:506
modparameters::isdowntypequark
logical function isdowntypequark(PartType)
Definition: mod_Parameters.F90:2330
modparameters::isazdecay
logical function isazdecay(DKMode)
Definition: mod_Parameters.F90:2243
modparameters::gf
real(8), public gf
Definition: mod_Parameters.F90:247
modparameters::decaymode1
integer, public decaymode1
Definition: mod_Parameters.F90:17
modparameters::al_lep
real(8), public al_lep
Definition: mod_Parameters.F90:1059
modparameters::esq
real(dp), public esq
Definition: mod_Parameters.F90:251
modparameters::ghwwp4_prime
complex(8), public ghwwp4_prime
Definition: mod_Parameters.F90:653
modparameters::offshellv2
logical, public offshellv2
Definition: mod_Parameters.F90:28
modparameters::ezp_str_left
complex(8), public ezp_str_left
Definition: mod_Parameters.F90:612
modparameters::ghwwp1_prime4
complex(8), public ghwwp1_prime4
Definition: mod_Parameters.F90:632
modparameters::ghw1_prime5
complex(8), public ghw1_prime5
Definition: mod_Parameters.F90:464
modparameters::num_
integer, target, public num_
Definition: mod_Parameters.F90:1098
modparameters::ewp_ta_right
complex(8), public ewp_ta_right
Definition: mod_Parameters.F90:704
modparameters::wpwpmode
integer, parameter, public wpwpmode
Definition: mod_Parameters.F90:14
modparameters::ezp_el_right
complex(8), public ezp_el_right
Definition: mod_Parameters.F90:599
modparameters::cz_q1sq
integer, public cz_q1sq
Definition: mod_Parameters.F90:434
modparameters::use_dynamic_mg
logical, parameter, public use_dynamic_mg
Definition: mod_Parameters.F90:933
modparameters::customupperscalebound
real(8), public customupperscalebound
Definition: mod_Parameters.F90:21
modparameters::barn
real(8), parameter, public barn
Definition: mod_Parameters.F90:101
modparameters::gh2z1_prime4
complex(8), public gh2z1_prime4
Definition: mod_Parameters.F90:747
modparameters::kappa_top
complex(8), public kappa_top
Definition: mod_Parameters.F90:886
modparameters::ghwpwp1_prime2
complex(8), public ghwpwp1_prime2
Definition: mod_Parameters.F90:668
modparameters::cl_qup
real(8), public cl_qup
Definition: mod_Parameters.F90:1073
modparameters::zprime_zz_2
complex(8), public zprime_zz_2
Definition: mod_Parameters.F90:915
modparameters::ghz4_prime2
complex(8), public ghz4_prime2
Definition: mod_Parameters.F90:412
modparameters::ghzzp1_prime5
complex(8), public ghzzp1_prime5
Definition: mod_Parameters.F90:524
modparameters::daz_lep_left
real(8), public daz_lep_left
Definition: mod_Parameters.F90:999
modparameters::hig_
integer, target, public hig_
Definition: mod_Parameters.F90:1100
modparameters::c2w_q2sq
integer, public c2w_q2sq
Definition: mod_Parameters.F90:849
modparameters::gszpmode
integer, parameter, public gszpmode
Definition: mod_Parameters.F90:16
modparameters::tap_
integer, target, public tap_
Definition: mod_Parameters.F90:1092
modparameters::ghwpwp4_prime
complex(8), public ghwpwp4_prime
Definition: mod_Parameters.F90:691
modparameters::collider_energy
real(8), public collider_energy
Definition: mod_Parameters.F90:19
modparameters::gh2w1_prime6
complex(8), public gh2w1_prime6
Definition: mod_Parameters.F90:816
modparameters::ghzpzp3_prime3
complex(8), public ghzpzp3_prime3
Definition: mod_Parameters.F90:583
modparameters::bzzp9
complex(8), public bzzp9
Definition: mod_Parameters.F90:955
modparameters::tr
real(dp), parameter, public tr
Definition: mod_Parameters.F90:1018
modparameters::invalidmode
integer, parameter, public invalidmode
Definition: mod_Parameters.F90:13
modparameters::countjets
integer function countjets(MY_IDUP)
Definition: mod_Parameters.F90:2587
modparameters::evalalphas
subroutine evalalphas()
Definition: mod_Parameters.F90:2850
modparameters::pol_mless
complex(dp) function, dimension(4) pol_mless(p, i, outgoing)
Definition: mod_Parameters.F90:2997
modparameters::ghzpzp4
complex(8), public ghzpzp4
Definition: mod_Parameters.F90:562
modparameters::hbar
real(8), parameter, public hbar
Definition: mod_Parameters.F90:109
modparameters::daz_lep_right
real(8), public daz_lep_right
Definition: mod_Parameters.F90:1000
modparameters::ga_z_ps
real(8), public ga_z_ps
Definition: mod_Parameters.F90:78
modparameters::ghwpwp4_prime4
complex(8), public ghwpwp4_prime4
Definition: mod_Parameters.F90:694
modparameters::ghwpwp1_prime3
complex(8), public ghwpwp1_prime3
Definition: mod_Parameters.F90:669
modparameters::br_z_uu_counter
integer, public br_z_uu_counter
Definition: mod_Parameters.F90:122
modparameters::decaymode2
integer, public decaymode2
Definition: mod_Parameters.F90:17
modparameters::krenfacscheme_mjj_mhstar
integer, parameter, public krenfacscheme_mjj_mhstar
Definition: mod_Parameters.F90:35
modparameters::gh2z4_prime3
complex(8), public gh2z4_prime3
Definition: mod_Parameters.F90:770
modparameters::top_
integer, target, public top_
Definition: mod_Parameters.F90:1088
modparameters::ghzzp2_prime2
complex(8), public ghzzp2_prime2
Definition: mod_Parameters.F90:529
modparameters::gsgsmode
integer, parameter, public gsgsmode
Definition: mod_Parameters.F90:13
modparameters::coupltolhewm
logical function coupltolhewm(Part)
Definition: mod_Parameters.F90:2663
modparameters::quarkcolavg
real(8), parameter, public quarkcolavg
Definition: mod_Parameters.F90:91
modparameters::userseed
integer, public userseed
Definition: mod_Parameters.F90:130
modparameters::ghwwp2_prime4
complex(8), public ghwwp2_prime4
Definition: mod_Parameters.F90:640
modparameters::overallcouplvffsq
real(8), public overallcouplvffsq
Definition: mod_Parameters.F90:1057
modparameters::printpmzz
complex(8), public printpmzz
Definition: mod_Parameters.F90:142
modparameters::setdecaywidth
subroutine setdecaywidth(width, ipart)
Definition: mod_Parameters.F90:1996
modparameters::wmpr_
integer, target, public wmpr_
Definition: mod_Parameters.F90:1119
modparameters::ghz2_prime2
complex(8), public ghz2_prime2
Definition: mod_Parameters.F90:396
modparameters::gh2gsgs4
complex(8), public gh2gsgs4
Definition: mod_Parameters.F90:741
modparameters::lambda2_z42
real(8), public lambda2_z42
Definition: mod_Parameters.F90:795
modparameters::ga_vprime
real(8), public ga_vprime
Definition: mod_Parameters.F90:78
modparameters::isalhegluon
logical function isalhegluon(PartType)
Definition: mod_Parameters.F90:2430
modparameters::hbbdecays
logical, public hbbdecays
Definition: mod_Parameters.F90:178
modparameters::ichann_sel
integer, public ichann_sel
Definition: mod_Parameters.F90:58
modparameters::ghw1_prime2
complex(8), public ghw1_prime2
Definition: mod_Parameters.F90:461
modparameters::gh2z2_prime4
complex(8), public gh2z2_prime4
Definition: mod_Parameters.F90:755
modparameters::ghz1_prime5
complex(8), public ghz1_prime5
Definition: mod_Parameters.F90:391
modparameters::cw_q12sq
integer, public cw_q12sq
Definition: mod_Parameters.F90:499
modparameters::higgsdecaylengthmm
real(8), public higgsdecaylengthmm
Definition: mod_Parameters.F90:235
modparameters::ghz2_prime3
complex(8), public ghz2_prime3
Definition: mod_Parameters.F90:397
modparameters::avgcs
real(8), public avgcs
Definition: mod_Parameters.F90:119
modparameters::ghzzp3_prime4
complex(8), public ghzzp3_prime4
Definition: mod_Parameters.F90:539
modparameters::c2z_q1sq
integer, public c2z_q1sq
Definition: mod_Parameters.F90:785
modparameters::pho_
integer, target, public pho_
Definition: mod_Parameters.F90:1094
modparameters::ghgsgs3
complex(8), public ghgsgs3
Definition: mod_Parameters.F90:424
modparameters::lambda_w21
real(8), public lambda_w21
Definition: mod_Parameters.F90:501
modparameters::ghw3_prime
complex(8), public ghw3_prime
Definition: mod_Parameters.F90:476
modparameters::z0_
integer, target, public z0_
Definition: mod_Parameters.F90:1095
modparameters::ghzzp4_prime6
complex(8), public ghzzp4_prime6
Definition: mod_Parameters.F90:549
modparameters::qdr
real(8), parameter, public qdr
Definition: mod_Parameters.F90:1039
modparameters::generateevents
logical, public generateevents
Definition: mod_Parameters.F90:29
modparameters::ewp_el_left
complex(8), public ewp_el_left
Definition: mod_Parameters.F90:699
modparameters::br_z_cc
real(8), parameter, public br_z_cc
Definition: mod_Parameters.F90:292
modparameters::ghwpwp1_prime4
complex(8), public ghwpwp1_prime4
Definition: mod_Parameters.F90:670
modparameters::pdfadn_
integer, parameter, public pdfadn_
Definition: mod_Parameters.F90:1131
modparameters::ghwwp3_prime6
complex(8), public ghwwp3_prime6
Definition: mod_Parameters.F90:650
modparameters::ghwpwp3_prime5
complex(8), public ghwpwp3_prime5
Definition: mod_Parameters.F90:687
modparameters::ghw3_prime7
complex(8), public ghw3_prime7
Definition: mod_Parameters.F90:482
modparameters::gh2w4_prime3
complex(8), public gh2w4_prime3
Definition: mod_Parameters.F90:837
modparameters::process
integer, public process
Definition: mod_Parameters.F90:17
modparameters::kappa_tilde_bot
complex(8), public kappa_tilde_bot
Definition: mod_Parameters.F90:889
modparameters::ghw4_prime2
complex(8), public ghw4_prime2
Definition: mod_Parameters.F90:485
modparameters::ne
complex(8), parameter, public ne
Definition: mod_Parameters.F90:89
modparameters::vegasnc0_default
integer, public vegasnc0_default
Definition: mod_Parameters.F90:22
modparameters::pmzz_mreso
real(8), public pmzz_mreso
Definition: mod_Parameters.F90:138
modparameters::jhugen_version
character(len= *), parameter jhugen_version
Definition: mod_Parameters.F90:6
modparameters::ahg3
complex(8), parameter, public ahg3
Definition: mod_Parameters.F90:364
modparameters::lambda2_z4
real(8), parameter, public lambda2_z4
Definition: mod_Parameters.F90:781
modparameters::br_z_ll
real(8), parameter, public br_z_ll
Definition: mod_Parameters.F90:288
modparameters::gh2w1_prime
complex(8), public gh2w1_prime
Definition: mod_Parameters.F90:811
modparameters::bgsgs3
complex(8), public bgsgs3
Definition: mod_Parameters.F90:977
modparameters::vprimedecaylengthmasscutofffactor
real(8), public vprimedecaylengthmasscutofffactor
Definition: mod_Parameters.F90:146
modparameters::ghz4_prime3
complex(8), public ghz4_prime3
Definition: mod_Parameters.F90:413
modparameters::gluoncolavg
real(8), parameter, public gluoncolavg
Definition: mod_Parameters.F90:91
modparameters::taudecays
integer, public taudecays
Definition: mod_Parameters.F90:17
modparameters::lambda2_z21
real(8), public lambda2_z21
Definition: mod_Parameters.F90:789
modparameters::filesopened
logical, public filesopened
Definition: mod_Parameters.F90:31
modparameters::lambda2_w1
real(8), parameter, public lambda2_w1
Definition: mod_Parameters.F90:843
modparameters::widthscheme
integer, public widthscheme
Definition: mod_Parameters.F90:131
modparameters::accepcounter_part
integer(8), dimension(-6:6,-6:6), public accepcounter_part
Definition: mod_Parameters.F90:52
modparameters::pdfabot_
integer, parameter, public pdfabot_
Definition: mod_Parameters.F90:1135
modparameters::ghwpwp1_prime7
complex(8), public ghwpwp1_prime7
Definition: mod_Parameters.F90:673
modparameters::ghwwp1_prime7
complex(8), public ghwwp1_prime7
Definition: mod_Parameters.F90:635
modparameters::ezp_ta_right
complex(8), public ezp_ta_right
Definition: mod_Parameters.F90:603
modparameters::io_histofile
integer, parameter io_histofile
Definition: mod_Parameters.F90:1140
modparameters::br_z_dd
real(8), parameter, public br_z_dd
Definition: mod_Parameters.F90:293
modparameters::ghw4_prime7
complex(8), public ghw4_prime7
Definition: mod_Parameters.F90:490
modparameters::cw_q1sq
integer, public cw_q1sq
Definition: mod_Parameters.F90:497
modparameters::br_z_uu
real(8), parameter, public br_z_uu
Definition: mod_Parameters.F90:291
modparameters::ezp_bot_left
complex(8), public ezp_bot_left
Definition: mod_Parameters.F90:614
modparameters::readpmzz
logical, public readpmzz
Definition: mod_Parameters.F90:136
modparameters::br_z_bb
real(8), parameter, public br_z_bb
Definition: mod_Parameters.F90:295
modparameters::brhadr_w_ud
real(8), parameter, public brhadr_w_ud
Definition: mod_Parameters.F90:324
modparameters::c2
complex(8), parameter, public c2
Definition: mod_Parameters.F90:989
modparameters::ga_tau
real(8), public ga_tau
Definition: mod_Parameters.F90:243
modparameters::ghzzp4_prime
complex(8), public ghzzp4_prime
Definition: mod_Parameters.F90:544
modparameters::al_neu
real(8), public al_neu
Definition: mod_Parameters.F90:1061
modparameters::mup_
integer, target, public mup_
Definition: mod_Parameters.F90:1091
modparameters::ar_qup
real(8), public ar_qup
Definition: mod_Parameters.F90:1062
modparameters::vckm_td
real(8), public vckm_td
Definition: mod_Parameters.F90:284
modparameters::ghz4
complex(8), public ghz4
Definition: mod_Parameters.F90:384
modparameters::requestnjets
integer, dimension(1:2), public requestnjets
Definition: mod_Parameters.F90:25
modparameters::lambda2_w21
real(8), public lambda2_w21
Definition: mod_Parameters.F90:852
modparameters::numhistograms
integer, public numhistograms
Definition: mod_Parameters.F90:23
modparameters::gethiggsmass
real(8) function gethiggsmass(jH)
Definition: mod_Parameters.F90:2121
modparameters::gh2z4_prime2
complex(8), public gh2z4_prime2
Definition: mod_Parameters.F90:769
Generated by   doxygen 1.8.17