search

impy-project / chromo

Hadronic Interaction Model interface in PYthon
Other
30 stars 7 forks source link

Dpmjet does not boost intermediate particles #177

Open jncots opened 1 year ago

jncots commented 1 year ago

This issue is related to #170. In DpmJet only final particles are boosted to lab frame. Thus, on graph illustrating history the energies of intermediate (decayed) particles are wrong. I found the place where it could be changed (in Dpmjet code itself):

C in SUBROUTINE DT_LT2LAB

C This

      DO i = NPOint(4) , NHKk
         IF ( (ABS(ISThkk(i)).EQ.1) .OR. (ISThkk(i).EQ.1000) .OR. 
     &        (ISThkk(i).EQ.1001) ) THEN
            CALL DT_LTNUC(PHKk(3,i),PHKk(4,i),pz,pe,-3)
            PHKk(3,i) = pz
            PHKk(4,i) = pe
         END IF
      END DO

C should be changed to

      DO i = 1 , NHKk
            CALL DT_LTNUC(PHKk(3,i),PHKk(4,i),pz,pe,-3)
            PHKk(3,i) = pz
            PHKk(4,i) = pe
      END DO