Add Cerenkov photons to TVirtualMCStack.

[u/mrichter/AliRoot.git] / TFluka / crnkvp.f

diff --git a/TFluka/crnkvp.f b/TFluka/crnkvp.f
index 987769fac2783114b9781a15ce329ffae02b4b79..9f596766f84b3d6cd2628992edf8cc4bd329e5c7 100644 (file)
--- a/TFluka/crnkvp.f
+++ b/TFluka/crnkvp.f
@@ -56,8 +56,8 @@
 *      IF ( LEMAGN )
 *     &   CALL FLABRT ( 'CRNKVP', ' STOP:LEMAGN-NOT-YET-IMPLEMENTED' )
 *  No change of lattice check implemented for the moment:
-      IF ( LT1TRK .NE. LT2TRK ) CALL FLABRT ( 'CRNKVP',
-     &    ' STOP:LT1TRK.NE.LT2TRK-NOT-YET-IMPLEMENTED' )
+*      IF ( LT1TRK .NE. LT2TRK ) CALL FLABRT ( 'CRNKVP',
+*     &    ' STOP:LT1TRK.NE.LT2TRK-NOT-YET-IMPLEMENTED' )
       DEDXCK = ZERZER
       MMAT   = MEDIUM (MREG)
       IF ( JTRACK .GE. -6 ) THEN
@@ -410,8 +410,23 @@ D    &              SCADOT
          AGOPPH (LSTOPP) = ATRKCR
          CMPOPP (LSTOPP) = ZERZER
          LOOPPH (LSTOPP) = LTRACK + 1
+*
+*
+*        Hook to TFluka
+*
+         PXCR =  EPHSMP * TXOPPH (LSTOPP)
+         PYCR =  EPHSMP * TYOPPH (LSTOPP)
+         PZCR =  EPHSMP * TZOPPH (LSTOPP)
+         POX  = TXPOPP(LSTOPP)
+         POY  = TYPOPP(LSTOPP)
+         POZ  = TZPOPP(LSTOPP)
+         CALL PushCerenkovPhoton(PXCR, PYCR, PZCR, EPHSMP, XTRKCR, 
+     &        YTRKCR , ZTRKCR, ATRKCR, POX, POY, POZ, WTRACK, ITFL)
+*
+*
+*
 *  |  !!!!!! Here Stuprf should be used !!!!!!
-         LOUOPP (LSTOPP) = LLOUSE
+         LOUOPP (LSTOPP) = ITFL
          DO 2100 ISPR = 1, MKBMX1
             SPAROK (ISPR,LSTOPP) = SPAUSR (ISPR)
  2100    CONTINUE