Commit 9ae752a8 authored by Prabhakar Palni's avatar Prabhakar Palni
Browse files

removed block of comments

parent 8e8dff56
......@@ -78,176 +78,6 @@ ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
include 'sAA.f'
include 'mphi.f'
ccccccc
cccccccccccccccccccccccccccccccccccccccccccccccccccc
cccccccccccccccccccccccccccccccccccccccccccccccccccc
c
c forward=.false.
c
c call length(outtag,outl)
c
c open(45,file='evrecs/evrec'//outtag(1:outl)//'.dat')
c wmax=0d0
c evnum=0
c
c if(genunw)then
c else
c readwt=.false.
c endif
c if(readwt)wmax=wtmax
c
c if(erec.eq.'hepmc')then
c erech=.true.
c erec='lhe'
c endif
c
c iw=0
c
c gf=1.16639d-5
c v=dsqrt(1d0/dsqrt(2d0)/gf)
c mt=173d0
c mb=4.75d0
c mc=1.4d0
c mmu=0.10566d0
c mpsi=3.096916d0
c mpsip=3.686109d0
c mups=9.46030d0
c mchic0=3.41475d0
c mchib0=9.85944d0
c mp=0.938272046d0
c mn=0.939565413d0
c mw=80.318d0
c me=0.511d-3
c mtau=1.77682d0
c mpip=0.13957018d0
c mkp=0.493677d0
c mphi=1.019461d0
c alpha=7.2974d-3
c pi=dacos(-1d0)
c conv=389379d3
c zi=(0d0,1d0)
c mup=0.062d0
c md=0.083d0
c ms=0.215d0
c
c rmf1( 1) = 1d-10
c rmf1( 2) = me
c rmf1( 3) = 1d-10
c rmf1( 4) = mmu
c rmf1( 5) = 1d-10
c rmf1( 6) = mtau
c rmf1( 7) = 0.062d0
c rmf1( 8) = 0.083d0
c rmf1( 9) = mc
c rmf1(10) = 0.215d0
c rmf1(11) = mt
c rmf1(12) = mb
c
c rmf1( 1) = me
c rmf1( 2) = mmu
c rmf1( 3) = mtau
c rmf1( 4) = md
c rmf1( 5) = mup
c rmf1( 6) = ms
c rmf1( 7) = mc
c rmf1( 8) = mb
c rmf1( 9) = mt
c
c
c mq=0d0
c hel=1
c mes=.false.
c mfact='mx'
c forward=.false.
c decay2=.false.
c decay3=.false.
c decay4=.false.
c decay6=.false.
c
cccccccccccc
c
c do i=1,20
c jdahep(1,i)=0
c jdahep(2,i)=0
c enddo
c
cccccccccccccccccccccccccc
c call inpdf
c call supinit
c call header
cccccccccccccccccccccccccc
c if(mes)then
c call calcmes
c call wfinit
c endif
ccccccccccccccccccccccccc
c if(beam.eq.'el')then
c if(sfaci)then
c print*,'Error : must have sfaci = .false. for initial-state
c &electrons'
c stop
c endif
c endif
ccccccccccccccccccccccccccccccccccccccccccccccccccc
ccccccccccccccccccccccccccccccccccccccccccccccccccc
c s=rts**2
c if(beam.eq.'prot')then
c beta=dsqrt(1d0-4d0*mp**2/s)
c elseif(beam.eq.'el')then
c beta=dsqrt(1d0-4d0*me**2/s)
c elseif(beam.eq.'ion'.or.beam.eq.'ionp')then
c mion=mp*an
c rtsi=rts
c si=s
c endif
c
c q(1,1)=0d0
c q(2,1)=0d0
c q(3,1)=rts/2d0*beta
c q(4,1)=rts/2d0
c
c q(1,2)=0d0
c q(2,2)=0d0
c q(3,2)=-rts/2d0*beta
c q(4,2)=rts/2d0
c
c if(beam.eq.'ionp')call pAinit
c if(beam.eq.'ion')call AAinit
c if(beam.eq.'prot')then
c pdgid(1)=2212
c pdgid(2)=2212
c pdgid(3)=2212
c pdgid(4)=2212
c elseif(beam.eq.'el')then
c pdgid(1)=11
c pdgid(2)=-11
c pdgid(3)=11
c pdgid(4)=-11
c elseif(beam.eq.'ion')then
c pdgid(1)=1000000000
c pdgid(1)=pdgid(1)+nint(az)*10000
c pdgid(1)=pdgid(1)+nint(an)*10
c pdgid(2)=pdgid(1)
c pdgid(3)=pdgid(1)
c pdgid(4)=pdgid(1)
c elseif(beam.eq.'ionp')then
c pdgid(2)=1000000000
c pdgid(2)=pdgid(2)+nint(az)*10000
c pdgid(2)=pdgid(2)+nint(an)*10
c pdgid(1)=2212
c pdgid(3)=pdgid(1)
c pdgid(4)=pdgid(2)
c endif
ccccccccccccccccccccccccccccccccccccccccccccccc
cccc HEPEVT
ccccccccccccccccccccccccccccccccccccccccccccccc
......
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