Ticket #793: E500.P4f_ww_h.eL.pR.sin

model = SM_CKM

alias uq = u:U
alias dq = d:D
alias cq = c:C
alias sq = s:S
alias bq = b:B

alias up_type_q = u:c:U:C
alias down_type_q = d:s:b:D:S:B
alias quarks = u:d:s:c:b:U:D:S:C:B

alias lepton = e1:e2:e3:E1:E2:E3
alias neutrino = n1:n2:n3:N1:N2:N3
alias fermions = quarks:lepton

alias neutral = pi0:dd1:a00:a20:rho316900:a420400:rho319900:rho322500:rho0:b10:pi216700:a217000:pi221000:rho14500:a10:rho17000:a116400:a0980:pi114000:pi116000:rho19000:rho2150:pi1300:pi1300:eta:f0600:f0980:f01500:f02020:f02100:f02200:eta2225:f02330:ETA1295:f1370:ETA1760:f20:omega:h10:f10:omega1420:f11510:h11595:omega1650:f2:f21430:f21565:f21640:ETA21645:f21810:f21910:f22150:f22300:f22340:f21950:omega3:f4:f42300:etaprime:f01710:eta1440:phi:h10p:f10p:phi1680:F2p:phi3:eta21870:f22010:fj2220:k0:k0star0:k14600:k18300:k0star19500:kL0:kS0:kstar0:k10:k1star0:kstar16800:k114100:k116500:k2star0:k217700:k218200:k2star19800:k215800:k222500:k3star0:k323200:k4star0:k425000:d0:dbar0:d0star0:dbar0star0:d2star0:dbar2star0:dstar0:d10:d1star0:b0:bbar0:b0star0:bbar0star0:bstar0:bstarbar0:b1l0:b1l0bar:b1h0:b1h0bar:b2star0:b2starbar0:bs0:bsbar0:bs0:bsbar0:bstars0:bstarsbar0:bs1l0:bs1lbar0:bs1h0:bs1hbar0:bstars20:bstarsbar0:etac:chic0:etac2s:jpsi:h10c:chic1:psiprime:psi3770:psi4040:psi4160:psi4415:chi2c:psi3836:etab:chib01p:etab2s:chib02p:etab3s:chib03p:upsilon:hb1p:chib11p:upsilon11d:upsilonprime:hb2p:chib12p:upsilon12d:upsilon3s:hb3p:chib13p:upsilon4s:upsilon10860:upsilon11020:chib21p:etab21d:upsilon21d:chib22p:etab22d:upsilon22d:chib23p:upsilon31d:upsilon32d:n:nbar:Delta0:Deltabar0:lambda:lambdabar:sigma:sigmabar:sigmastar0:sigmastarbar0:xi0:xibar0:xistar0:xistarbar0:sigmac:sigmacbar:sigmastarc:sigmastarcbar:xic0:xicbar0:xiprimec0:xiprimecbar0:xistarc0:xistarcbar0:omegac:omegacbar:omegastarc:omegastarcbar:lambdab:lambdabbar:sigmab:sigmabbar:sigmastarb:sigmastarbbar:sigmabp:xib0:xibbar0:xiprimeb0:xiprimebbar0:xistarb0:xistarbbar0:xibc0:xibcbar0:xibcprime0:xibcprimebar0:xibcstar0:xibcstarbar0:omegabc0:omegabcbar0:omegabcprime0:omegabcprimebar0:omegabcstar0:omegabcstarbar0:sigmabb0:sigmabbbar0:sigmabbstar0:sigmabbstarbar0:omegabb0:omegabbbar0:omegabbc0:omegabbcbar0:omegabbcstar0:omegabbcstarbar0


scale = mW

show (scale)

#   mu = 0 GeV
#   md = 0 GeV
  ms = 0 GeV
  mc = 0 GeV
  mb = 0 GeV 
  alphas = 0
#
#      e+e- --> mu+mu- and tau+tau-
#
polarized e1,E1,e2,E2,e3,E3,n1,N1,n2,N2,n3,N3,u,U,d,D,s,S,c,C,b,B,t,T, "W+", "W-", Z

process P4f_ww_h = e1, E1 => ( cq,down_type_q,uq,down_type_q ) + (uq,bq,uq,dq) + (uq,sq,uq,dq) + ( cq,bq,cq,sq ) + ( cq,dq,cq,sq ) { process_num_id = 250006 }

compile

#
logical ?is_higgs_process = false
include ("pythia6-parameters.sin")
 mH = 2000.
 wH = 10.
 alphas = 1.e-6
 isr_alpha = 0.0072993
 isr_mass = 0.000511
 epa_alpha = 0.0072993

#   mH = 2000 GeV
#   wH = 10 GeV
#   ?ps_fsr_active = true
#   ?hadronization_active = true
#   $shower_method = "PYTHIA6"

sqrts = 500 GeV

beams = e1, E1 => isr
# $beam_events_file = "/group/ilc/soft/gcc481/GP/beamdata/lumiee500-hs/lumiee500-run1to20.dat"
beams_pol_density = @(-1), @(+1)
?isr_recoil = true
?keep_remnants = true
?keep_beams = true

real default_M_cut = 4 GeV
real default_jet_cut = 10 GeV
real default_E_cut = 10 GeV
real default_Q_cut = 4 GeV
alias uncoloured=lepton:neutrino:photon
cuts =
  all M > default_jet_cut [quarks,quarks]

#   all M > default_M_cut [uncoloured,uncoloured]
#  and
#  all M > default_jet_cut [quarks,quarks]
#  and
#  all E > default_E_cut [photon:gluon]
#  and
#  all M < - default_Q_cut [incoming particle, fermions]
#

# real default_M_cut = 4 GeV
# cuts = 
#   all M > default_M_cut [ quarks, quarks ]


integrate ( P4f_ww_h ) { ?polarized_events = true }

printf "#Total cross section : P4f_ww_h %13e fb" ( integral( P4f_ww_h ) )

show (model)

# sample_format = lcio
sample_format = stdhep_ev4
$sample = "E500.P4f_ww_h.Gwhizard2.eL.pR"


histogram n_photon (0, 150, 2)  { $title = "Distribution of $\gamma$" }
histogram n_ch (0, 150, 2)      { $title = "Charged particles" }
histogram n_charged (0, 150, 2) { $title = "Charged hadrons, no protons" }
histogram n_pi0 (0, 150, 2)     { $title = "$\pi^0$ distribution" }
histogram n_neutral (-1, 150, 2) { $title = "Neutral particle distribution" }

model = SM_hadrons


analysis = record n_photon (count [A]);
           record n_ch (count [charged]);
           record n_charged (count if abs (PDG) <> 2212 [charged]);
           record n_pi0 (count [pi0]);
           record n_neutral (count [neutral])

n_events = 1
simulate ( P4f_ww_h ) { ?polarized_events = true } 

compile_analysis ()