from Alex Kalweit:
[u/mrichter/AliRoot.git] / TPC / Calib / AliTPCcalibGainMult.h
index 1d4fa167c7aa0b4cb65f9df4612db2eb357d8dcf..4e0eefc73b61278a02189b054e23c2edbc992a01 100644 (file)
@@ -57,6 +57,16 @@ public:
   void SetUpperTrunc(Float_t upperTrunc){fUpperTrunc = upperTrunc;};
   void SetUseMax(Bool_t useMax){fUseMax = useMax;};
   //
   void SetUpperTrunc(Float_t upperTrunc){fUpperTrunc = upperTrunc;};
   void SetUseMax(Bool_t useMax){fUseMax = useMax;};
   //
+  void SetCutMinCrossRows(Int_t crossRows){fCutCrossRows = crossRows;};
+  void SetCutMaxEta(Float_t maxEta){fCutEtaWindow = maxEta;};
+  void SetCutRequireITSrefit(Bool_t requireItsRefit = kFALSE){fCutRequireITSrefit = requireItsRefit;};
+  void SetCutMaxDcaXY(Float_t maxXY){fCutMaxDcaXY = maxXY;}; 
+  void SetCutMaxDcaZ(Float_t maxZ){fCutMaxDcaZ = maxZ;}; 
+  //
+  void SetMinMomentumMIP(Float_t minMom = 0.4){fMinMomentumMIP = minMom;};
+  void SetMaxMomentumMIP(Float_t maxMom = 0.6){fMaxMomentumMIP = maxMom;};
+  void SetAlephParameters(Float_t * parameters){for(Int_t j=0;j<5;j++) fAlephParameters[j] = parameters[j];};
+  //
   //
   void     Process(AliESDtrack *track, Int_t runNo=-1){AliTPCcalibBase::Process(track,runNo);};
   void     Process(AliTPCseed *track){return AliTPCcalibBase::Process(track);}
   //
   void     Process(AliESDtrack *track, Int_t runNo=-1){AliTPCcalibBase::Process(track,runNo);};
   void     Process(AliTPCseed *track){return AliTPCcalibBase::Process(track);}
@@ -82,12 +92,27 @@ private:
   //
   Bool_t fUseMax;                 // flag if Qmax or Qtot should be used
   //
   //
   Bool_t fUseMax;                 // flag if Qmax or Qtot should be used
   //
+  // track cuts
+  //
+  Int_t   fCutCrossRows;                // minimum number of crossed rows 
+  Float_t fCutEtaWindow;                // maximum eta of tracks
+  Bool_t  fCutRequireITSrefit;          // if ITSrefit should be required (dangerous in cpass0)
+  Float_t fCutMaxDcaXY;                 // max dca_xy (only TPConly resolution is guaranteed!)
+  Float_t fCutMaxDcaZ;                  // max dca_z  (dangerous if vDrift is not calibrated)
+  //
+  // definition of MIP window
+  //
+  Float_t fMinMomentumMIP;              // minimum momentum of MIP region, e.g. 400 MeV
+  Float_t fMaxMomentumMIP;              // maximum momentum of MIP region, e.g. 600 MeV
+  Float_t fAlephParameters[5];          // parameters for equalization in MIP window, parameter set should be =1 at MIP
+  //
   // histograms
   //
   TH1F  *fHistNTracks;            //  histogram showing number of ESD tracks per event
   TH1F  *fHistClusterShape;       //  histogram to check the cluster shape
   TH3F  *fHistQA;                 //  dE/dx histogram showing the final spectrum
   //
   // histograms
   //
   TH1F  *fHistNTracks;            //  histogram showing number of ESD tracks per event
   TH1F  *fHistClusterShape;       //  histogram to check the cluster shape
   TH3F  *fHistQA;                 //  dE/dx histogram showing the final spectrum
   //
+  //
   THnSparseF * fHistGainSector;   //  histogram which shows MIP peak for each of the 3x36 sectors (pad region)
   THnSparseF * fHistPadEqual;     //  histogram for the equalization of the gain in the different pad regions -> pass0
   THnSparseF * fHistGainMult;     //  histogram which shows decrease of MIP signal as a function
   THnSparseF * fHistGainSector;   //  histogram which shows MIP peak for each of the 3x36 sectors (pad region)
   THnSparseF * fHistPadEqual;     //  histogram for the equalization of the gain in the different pad regions -> pass0
   THnSparseF * fHistGainMult;     //  histogram which shows decrease of MIP signal as a function
@@ -102,7 +127,7 @@ private:
   AliTPCcalibGainMult(const AliTPCcalibGainMult&); 
   AliTPCcalibGainMult& operator=(const AliTPCcalibGainMult&); 
 
   AliTPCcalibGainMult(const AliTPCcalibGainMult&); 
   AliTPCcalibGainMult& operator=(const AliTPCcalibGainMult&); 
 
-  ClassDef(AliTPCcalibGainMult, 2); 
+  ClassDef(AliTPCcalibGainMult, 3); 
 };
 
 #endif
 };
 
 #endif