remove Fluka related Gstpar parameters. Waiting for new tests.

diff --git a/TRD/AliTRD.cxx b/TRD/AliTRD.cxx
index 9b2adbe88bf8b18ed5450817ab7d599ee0f74941..3fd1b38f6b9dc816adbefd9e4a0fda92743047c7 100644 (file)
--- a/TRD/AliTRD.cxx
+++ b/TRD/AliTRD.cxx
@@ -37,6 +37,7 @@
 #include "AliTRDhit.h"
 #include "AliTRDrawData.h"
 #include "AliTRDSimParam.h"
+#include "TTreeStream.h"
 
 ClassImp(AliTRD)
  
@@ -47,6 +48,7 @@ AliTRD::AliTRD()
   ,fGasDensity(0)
   ,fFoilDensity(0)
   ,fGasNobleFraction(0)
+  ,fDebugStream(0x0)
 {
   //
   // Default constructor
@@ -61,6 +63,7 @@ AliTRD::AliTRD(const char *name, const char *title)
   ,fGasDensity(0)
   ,fFoilDensity(0)
   ,fGasNobleFraction(0)
+  ,fDebugStream(0x0)
 {
   //
   // Standard constructor for the TRD
@@ -86,7 +89,6 @@ AliTRD::AliTRD(const char *name, const char *title)
   // Allocate the hit array
   fHits = new TClonesArray("AliTRDhit",405);
   gAlice->GetMCApp()->AddHitList(fHits);
-
 }
 
 //_____________________________________________________________________________
@@ -95,7 +97,7 @@ AliTRD::~AliTRD()
   //
   // TRD destructor
   //
-
+  delete fDebugStream;
   if (fGeometry) {
     delete fGeometry;
     fGeometry = 0;
@@ -120,7 +122,12 @@ void AliTRD::Hits2Digits()
 
   // Initialization
   digitizer.InitDetector();
-    
+
+  TDirectory *savedir = gDirectory;
+  fDebugStream = new TTreeSRedirector("TRD.SimDebug.root");
+  savedir->cd();
+  digitizer.SetDebugStream(fDebugStream);
+  
   if (!fLoader->TreeH()) {
     fLoader->LoadHits("read");
   }
@@ -154,7 +161,12 @@ void AliTRD::Hits2SDigits()
 
   // Initialization
   digitizer.InitDetector();
-    
+
+  TDirectory *savedir = gDirectory;
+  fDebugStream = new TTreeSRedirector("TRD.SimDebug.root");
+  savedir->cd();
+  digitizer.SetDebugStream(fDebugStream);
+  
   if (!fLoader->TreeH()) {
     fLoader->LoadHits("read");
   }
@@ -530,8 +542,9 @@ void AliTRD::Init()
   gMC->Gstpar((* fIdtmed)[9],"DRAY"    , 1.0);
   gMC->Gstpar((* fIdtmed)[9],"STRA"    , 1.0); 
   gMC->Gstpar((* fIdtmed)[9],"LOSS"    ,13.0);      // Specific energy loss
-  gMC->Gstpar((* fIdtmed)[9],"PRIMIO_E",23.53);     // 1st ionisation potential
-  gMC->Gstpar((* fIdtmed)[9],"PRIMIO_N",19.344431); // Number of primaries
+  // Parameters related to Fluka MC machine
+  //gMC->Gstpar((* fIdtmed)[9],"PRIMIO_E",23.53);     // 1st ionisation potential
+  //gMC->Gstpar((* fIdtmed)[9],"PRIMIO_N",19.344431); // Number of primaries
 
 }