Cerenkov photon update.

[u/mrichter/AliRoot.git] / TFluka / TFluka.cxx

diff --git a/TFluka/TFluka.cxx b/TFluka/TFluka.cxx
index 329e2b643b4c877f03c49655320237cf5564e12d..f9808f00a2acdae803ce2e5daa98c2e7929aff6e 100644 (file)
--- a/TFluka/TFluka.cxx
+++ b/TFluka/TFluka.cxx
@@ -206,8 +206,8 @@ TFluka::~TFluka()
         cout << "<== TFluka::~TFluka() destructor called." << endl;
     if (fMaterials) delete [] fMaterials;
     
-    delete fGeom;
-    delete fMCGeo;
+//    delete fGeom;
+//    delete fMCGeo;
     
     if (fUserConfig) {
         fUserConfig->Delete();
@@ -2095,8 +2095,11 @@ Int_t TFluka::StepProcesses(TArrayI &proc) const
     
     proc.Set(1);
     TMCProcess iproc;
-    if (caller == kBXEntering || caller == kBXExiting || caller == kEEDRAW) {
+    if (caller == kBXEntering || caller == kBXExiting || caller == kEEDRAW || caller == kSODRAW) {
        iproc = kPTransportation;
+    }
+    else if (caller == kUSTCKV) {
+       iproc = kPCerenkov;
     } else {
        switch (icode) {
        case kEMFSCO:
@@ -2113,6 +2116,16 @@ Int_t TFluka::StepProcesses(TArrayI &proc) const
                iproc = kPTransportation;
            }
            break;
+       case kKASHEA:
+       case kKASNEU:
+       case kKASOPH:
+       case kKASKADescape:
+       case kEMFSCOescape:
+       case kKASNEUescape:
+       case kKASHEAescape:
+       case kKASOPHescape:
+           iproc = kPTransportation;
+           break;
        case kKASKADtimekill:
        case kEMFSCOtimekill:
        case kKASNEUtimekill:
@@ -2121,14 +2134,9 @@ Int_t TFluka::StepProcesses(TArrayI &proc) const
            iproc =  kPTOFlimit;
            break;
        case kKASKADstopping:
-       case kKASKADescape:
        case kEMFSCOstopping1:
        case kEMFSCOstopping2:
-       case kEMFSCOescape:
        case kKASNEUstopping:
-       case kKASNEUescape:
-       case kKASHEAescape:
-       case kKASOPHescape:
            iproc = kPStop;
            break;
        case kKASOPHabsorption:
@@ -2136,6 +2144,7 @@ Int_t TFluka::StepProcesses(TArrayI &proc) const
            break;
        case kKASOPHrefraction:
            iproc = kPLightRefraction;
+           break;
        case kEMFSCOlocaldep : 
            iproc = kPPhotoelectric;
            break;