//
//
//
- //delete fHistNTracks; // histogram showing number of ESD tracks per event
- //delete fHistGainSector; // histogram showing the number of clusters per track
-
+ delete fHistNTracks; // histogram showing number of ESD tracks per event
+ delete fHistClusterShape; // histogram to check the cluster shape
+ delete fHistQA; // dE/dx histogram showing the final spectrum
+ //
+ delete fHistGainSector; // histogram which shows MIP peak for each of the 3x36 sectors (pad region)
+ delete fHistPadEqual; // histogram for the equalization of the gain in the different pad regions -> pass0
+ delete fHistGainMult; // histogram which shows decrease of MIP signal as a function
+
}
fCutMaxDz(25), // maximal distance in rfi ditection
fCutTheta(0.03), // maximal distan theta
fCutMinDir(-0.99), // direction vector products
- fCutTracks(100),
+ fCutTracks(2500),
fArrayLaserA(0), //laser fit parameters C
fArrayLaserC(0), //laser fit parameters A
fArrayDz(0), //NEW! Tmap of V drifts for different triggers
fCutMaxDz(40), // maximal distance in rfi ditection
fCutTheta(5*0.004644),// maximal distan theta
fCutMinDir(-0.99), // direction vector products
- fCutTracks(100),
+ fCutTracks(2500),
fArrayLaserA(new TObjArray(1000)), //laser fit parameters C
fArrayLaserC(new TObjArray(1000)), //laser fit parameters A
fArrayDz(0), //Tmap of V drifts for different triggers
TObjArray * GetLaserArrayA() const { return fArrayLaserA;}
TObjArray * GetLaserArrayC() const { return fArrayLaserC;}
+
+ void SetCutTracks(Int_t maxTracks) { fCutTracks = maxTracks; } // set maximal number of tracks
+ Int_t GetCutTracks() const { return fCutTracks; } // retun maximal number of tracks
+
+
+
protected:
void ResetCurrent(); // reset current values
Int_t fMemoryMode; // 0 -do not fill THnSparse with residuals 1- fill only important QA THn 2 - Fill all THnsparse for calibration
if (seed) {
Int_t particleCase = 0;
- if (meanP > 0.5 || meanP < 0.4) particleCase = 2; // MIP pions
- if (meanP > 0.56 || meanP < 0.57) particleCase = 3; // protons 1
- if (meanP > 0.64 || meanP < 0.66) particleCase = 4; // protons 2
+ if (meanP < 0.5 && meanP > 0.4) particleCase = 2; // MIP pions
+ if (meanP < 0.57 && meanP > 0.56) particleCase = 3; // protons 1
+ if (meanP < 0.66 && meanP > 0.65) particleCase = 4; // protons 2
//
if (fLowMemoryConsumption && particleCase == 0) continue;
//