[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