pythiaeRHIC/radgen/radgen_event.f

0001
0002       subroutine radgen_event
0003
0004       implicit none
0005
0006        include "mc_set.inc"
0007        include "mconsp.inc"
0008        include "phiout.inc"
0009        include "tailcom.inc"
0010        include "cmpcom.inc"
0011        include "radgen.inc"
0012        include "mcRadCor.inc"
0013
0014       COMMON/PYPARS/MSTP(200),PARP(200),MSTI(200),PARI(200)
0015       INTEGER MSTP,MSTI
0016       REAL PARP,PARI
0017       SAVE /PYPARS/
0018
0019       real PhRAD(4),q2true,nutrue,radweight
0020       real nu,q2,phi,yys,xxs
0021
0022 ! calculate radiative corrections
0023       nu=sngl(gennu)
0024       q2=sngl(genq2)
0025       phi=sngl(genphi)
0026       yys=sngl(geny)
0027       xxs=sngl(genx)
0028       call RADGEN(mcSet_EneBeam,q2,nu,yys,xxs,phi,PhRAD,q2true,nutrue,
0029      +     radweight)
0030
0031 ! fill mcRadCor WCB with ...
0032
0033       mcRadCor_ID=1
0034
0035 ! ... true kinematics
0036
0037 * by definition we calculate xbj using the proton mass
0038 * ---> xbj for elastic events is A
0039 *     mcRadCor_XTrue=q2true/(2.0d0*amp*nutrue)
0040       mcRadCor_NuTrue=nutrue
0041       mcRadCor_Q2True=q2true
0042       mcRadCor_YTrue=nutrue/mcSet_EneBeam
0043       mcRadCor_XTrue=q2true/mcRadCor_YTrue/(4.*ebeam*pbeam)
0044       mcRadCor_W2True=amp2 + (q2true*(1./mcRadCor_XTrue-1.))
0045 !      mcRadCor_XTrue=q2true/(2.0d0*0.938272d0*nutrue)
0046 !      mcRadCor_W2True=amp2 - q2true + 2.*amp*nutrue
0047
0048 ! ... kinematics of real photon
0049       mcRadCor_EBrems=phrad(4)
0050       mcRadCor_ThetaBrems=0.
0051       if(phrad(4).gt.0.)
0052      +    mcRadCor_ThetaBrems = acos(phrad(3)/phrad(4))
0053       mcRadCor_PhiBrems=0.
0054       if (.not.(phrad(1).eq.0..and.phrad(2).eq.0.)) then
0055         mcRadCor_PhiBrems = atan2(phrad(2),phrad(1))
0056         if (mcRadCor_PhiBrems.lt.0.)
0057      +       mcRadCor_PhiBrems = mcRadCor_PhiBrems + twopi
0058       endif
0059
0060 c...if we would like to have the TSAI system angles for the real gamma than
0061 c     mcRadCor_ThetaBrems = dthg
0062 c     mcRadCor_PhiBrems = dphig
0063
0064 ! ... radiative contributions
0065
0066       mcRadCor_Sigrad=sigrad
0067       mcRadCor_Sigcor=sigcor
0068       mcRadCor_Sigcorerr=0.
0069       mcRadCor_TailIne=tine
0070       mcRadCor_TailEla=tpro
0071       mcRadCor_TailCoh=tnuc
0072       mcRadCor_Vacuum=vac
0073       mcRadCor_Vertex=vertex
0074       mcRadCor_Small=small
0075       mcRadCor_Redfac=redfac
0076
0077 ! ... radiative correction type
0078
0079       if (ita.eq.2) then
0080         mcRadCor_cType='elas'
0081       else if (ita.eq.3) then
0082         mcRadCor_cType='qela'
0083       else
0084         mcRadCor_cType='inel'
0085       endif
0086
0087       end