[Prev][Next][Index][Thread]
PH corrections
ciao Adriano e Monica,
le correzioni PH di Jon Link sul momento si implementano in
usrgeom
==================================================
Es:
SUBROUTINE USRGEOM
c for post hoc momentum correction
call phcmom_geom
return
end
=====================================================
e naturalmente vanno applicate sul vettore impulsi
===================================================
c momentum cor for D daughters
do ii=1,nezcdau(icnd)
if(iezmlst(ii,imndx).eq.4)then !4=K
I = iezcdau(ii,icnd)
cf = phcmom(I) ******************************************
PKX = cf * TRAK(13,I)
PKY = cf * TRAK(14,I)
PKZ = cf * TRAK(15,I)
PKT = PKX*PKX + PKY*PKY + PKZ*PKZ
EKAON = sqrt( AKMAS*AKMAS + PKT)
dmc_px(ii) = PKX
dmc_py(ii) = PKY
dmc_pz(ii) = PKZ
dmc_en(ii) = EKAON
endif
if(iezmlst(ii,imndx).eq.3)then !3=pi
I = iezcdau(ii,icnd)
cf = phcmom(I)
PKX = cf * TRAK(13,I)
PKY = cf * TRAK(14,I)
PKZ = cf * TRAK(15,I)
PKT = PKX*PKX + PKY*PKY + PKZ*PKZ
EKAON = sqrt( PIMAS*PIMAS + PKT)
dmc_px(ii) = PKX
dmc_py(ii) = PKY
dmc_pz(ii) = PKZ
dmc_en(ii) = EKAON
endif
enddo
ciao,
s
=========================================================
Stefano Bianco
Alte Energie - Laboratori Nazionali di Frascati dell'INFN
v.E.Fermi, 40 - 00044 Frascati (Roma) ITALY
ph. +39-06-94032793 fax +39-06-94032427