Commit 35b97fc4 authored by Prabhakar Palni's avatar Prabhakar Palni
Browse files

improvement in event generation and ALP process fix

parent 013a0a89
......@@ -78,6 +78,91 @@ ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
include 'sAA.f'
include 'mphi.f'
cccccccccccccccccccccccccccccccccccccccccccccccccccc
cccccccccccccccccccccccccccccccccccccccccccccccccccc
forward=.false.
call length(outtag,outl)
open(45,file='evrecs/evrec'//outtag(1:outl)//'.dat')
wmax=0d0
evnum=0
if(genunw)then
else
readwt=.false.
endif
if(readwt)wmax=wtmax
if(erec.eq.'hepmc')then
erech=.true.
erec='lhe'
endif
iw=0
gf=1.16639d-5
v=dsqrt(1d0/dsqrt(2d0)/gf)
mt=173d0
mb=4.75d0
mc=1.4d0
mmu=0.10566d0
mpsi=3.096916d0
mpsip=3.686109d0
mups=9.46030d0
mchic0=3.41475d0
mchib0=9.85944d0
mp=0.938272046d0
mn=0.939565413d0
mw=80.318d0
me=0.511d-3
mtau=1.77682d0
mpip=0.13957018d0
mkp=0.493677d0
mphi=1.019461d0
alpha=7.2974d-3
pi=dacos(-1d0)
conv=389379d3
zi=(0d0,1d0)
mup=0.062d0
md=0.083d0
ms=0.215d0
rmf1( 1) = 1d-10
rmf1( 2) = me
rmf1( 3) = 1d-10
rmf1( 4) = mmu
rmf1( 5) = 1d-10
rmf1( 6) = mtau
rmf1( 7) = 0.062d0
rmf1( 8) = 0.083d0
rmf1( 9) = mc
rmf1(10) = 0.215d0
rmf1(11) = mt
rmf1(12) = mb
rmf1( 1) = me
rmf1( 2) = mmu
rmf1( 3) = mtau
rmf1( 4) = md
rmf1( 5) = mup
rmf1( 6) = ms
rmf1( 7) = mc
rmf1( 8) = mb
rmf1( 9) = mt
mq=0d0
hel=1
mes=.false.
mfact='mx'
forward=.false.
decay2=.false.
decay3=.false.
decay4=.false.
decay6=.false.
ccccccccccccccccccccccccccccccccccccccccccccccc
cccc HEPEVT
ccccccccccccccccccccccccccccccccccccccccccccccc
......@@ -184,7 +269,7 @@ ccc NEW LHE init
enddo
enddo
ccccccccc
cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
if(proc.eq.18.or.proc.eq.19.or.proc.eq.20)then
......@@ -221,41 +306,12 @@ ccccccccc
nhep=nup
nup=nup-2
if(proc.eq.68)then
nhep=7
nup=7
endif
ccccccccc
c
c surv=1d0
c if(beam.eq.'prot'.or.ionqcd.eq.'coh'.or.ionqcd.eq.'incoh')then
c call initparsr(isurv)
c call readscreen
c if(beam.eq.'prot'.or.ionqcd.eq.'incoh')surv=1d0/norm**2
c endif
c
c if(qcd)then
c call calcsud
c call calchg
c endif
c
c
c if(beam.eq.'ion'.or.beam.eq.'ionp')then
c call ioninit
c endif
c
cccccccccccc
c if(beam.eq.'ion'.or.beam.eq.'ionp')call ioninit
c if(beam.eq.'ionp')then
c rts=rtspa
c s=spa
c elseif(beam.eq.'ion')then
c rts=rtsaa
c beta=dsqrt(1d0-4d0*mp**2/s)
c saa=an**2*s
c rts=dsqrt(saa)
c rts=rts
c s=saa
c endif
cccccccccccccccccccccccccccccccccccccccccccccc
nhist=0
nhistmax=20
......@@ -289,7 +345,7 @@ c endif
print*,'*****************************************************'
print*,'Vegas:initialisation run (outputs *bare* crosssection)'
print*,'******************************************************'
ncall=ncall1+10000
CALL VEGAS(cs,AVGI,SD,CHI2A)
......@@ -321,13 +377,12 @@ c endif
endif
ITMX=ITMX1
NCALL=NCALL1
NCALL=NCALL1+1000
avgi1=avgi
sd1=sd
ncall=ncall*iinc
c inccall=inccall*iinc
inccall=inccall
inccall=inccall*iinc
779 bin=.true.
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment