#include "AliTRDhit.h"
#include "AliTRDrawData.h"
#include "AliTRDSimParam.h"
-#include "TTreeStream.h"
ClassImp(AliTRD)
,fGasDensity(0)
,fFoilDensity(0)
,fGasNobleFraction(0)
- ,fDebugStream(0x0)
{
//
// Default constructor
,fGasDensity(0)
,fFoilDensity(0)
,fGasNobleFraction(0)
- ,fDebugStream(0x0)
{
//
// Standard constructor for the TRD
// Allocate the hit array
fHits = new TClonesArray("AliTRDhit",405);
gAlice->GetMCApp()->AddHitList(fHits);
+
}
//_____________________________________________________________________________
//
// TRD destructor
//
- delete fDebugStream;
+
if (fGeometry) {
delete fGeometry;
fGeometry = 0;
// Initialization
digitizer.InitDetector();
-
- TDirectory *savedir = gDirectory;
- fDebugStream = new TTreeSRedirector("TRD.SimDebug.root");
- savedir->cd();
- digitizer.SetDebugStream(fDebugStream);
-
+
if (!fLoader->TreeH()) {
fLoader->LoadHits("read");
}
// Initialization
digitizer.InitDetector();
-
- TDirectory *savedir = gDirectory;
- fDebugStream = new TTreeSRedirector("TRD.SimDebug.root");
- savedir->cd();
- digitizer.SetDebugStream(fDebugStream);
-
+
if (!fLoader->TreeH()) {
fLoader->LoadHits("read");
}
gMC->Gstpar((* fIdtmed)[9],"DRAY" , 1.0);
gMC->Gstpar((* fIdtmed)[9],"STRA" , 1.0);
gMC->Gstpar((* fIdtmed)[9],"LOSS" ,13.0); // Specific energy loss
- // 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
+ gMC->Gstpar((* fIdtmed)[9],"PRIMIO_E",23.53); // 1st ionisation potential
+ gMC->Gstpar((* fIdtmed)[9],"PRIMIO_N",19.344431); // Number of primaries
}