TPC/AliTPCcalibGainMult.h

   1 #ifndef ALITPCCALIBGAINMULT_H
   2 #define ALITPCCALIBGAINMULT_H
   3
   4 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
   5  * See cxx source for full Copyright notice                               */
   6
   7 #include "AliTPCcalibBase.h"
   8 #include "AliTPCCalPad.h"
   9 #include "TH3F.h"
  10 #include "TF1.h"
  11 #include "THnSparse.h"
  12 #include "TMatrixD.h"
  13 #include "TVectorD.h"
  14 class TH1F;
  15 class TList;
  16 class AliESDEvent;
  17 class AliESDtrack;
  18 class AliTPCseed;
  19
  20 #include "TTreeStream.h"
  21
  22
  23 class AliTPCcalibGainMult:public AliTPCcalibBase {
  24 public:
  25   AliTPCcalibGainMult();
  26   AliTPCcalibGainMult(const Text_t *name, const Text_t *title);
  27   virtual ~AliTPCcalibGainMult();
  28   void SetBBParam(TVectorD * param) {fBBParam=param;}
  29   //  virtual void Terminate();
  30   //
  31   virtual void           Process(AliESDEvent *event);
  32   virtual void           ProcessV0s(AliESDEvent *event);
  33   virtual void           ProcessCosmic(const AliESDEvent *event);
  34   virtual void           ProcessKinks(const AliESDEvent *event);
  35   virtual void           ProcessTOF(const AliESDEvent *event);
  36   virtual void           DumpHPT(const AliESDEvent *event);
  37   virtual Long64_t       Merge(TCollection *li);
  38   virtual void           Analyze();
  39   void                   DumpTrack(AliESDtrack * track, AliESDfriendTrack *ftrack, AliTPCseed * seed, Int_t index);
  40   //
  41   TH1F   *          GetHistNTracks() const {return fHistNTracks;};
  42   TH1F   *          GetHistClusterShape() const {return fHistClusterShape;};
  43   TH3F   *          GetHistQA() const {return fHistQA;};
  44   //
  45   THnSparseF *      GetHistGainSector() const {return fHistGainSector;};
  46   THnSparseF *      GetHistPadEqual() const {return fHistPadEqual;};
  47   THnSparseF *      GetHistGainMult() const {return fHistGainMult;};
  48   THnSparseF * GetHistdEdxMap() const { return fHistdEdxMap;}      // 4D dedx histogram
  49   THnSparseF * GetHistdEdxMax() const { return fHistdEdxMax;}      // 4D dedx histogram
  50   THnSparseF * GetHistdEdxTot() const { return fHistdEdxTot;}      // 4D dedx histogram
  51   TTree *      GetdEdxTree() const {return fdEdxTree;}         // tree for the later minimization
  52
  53   TGraphErrors* GetGainPerChamber(Int_t padRegion=1, Bool_t plotQA=kFALSE);
  54   //
  55   void SetMIPvalue(Float_t mip){fMIP = mip;};
  56   void SetLowerTrunc(Float_t lowerTrunc){fLowerTrunc = lowerTrunc;};
  57   void SetUpperTrunc(Float_t upperTrunc){fUpperTrunc = upperTrunc;};
  58   void SetUseMax(Bool_t useMax){fUseMax = useMax;};
  59   //
  60   //
  61   void     Process(AliESDtrack *track, Int_t runNo=-1){AliTPCcalibBase::Process(track,runNo);};
  62   void     Process(AliTPCseed *track){return AliTPCcalibBase::Process(track);}
  63   //
  64   void     MakeLookup(THnSparse * hist, Char_t * outputFile);
  65   //
  66   void     UpdateGainMap();
  67   void     UpdateClusterParam();
  68
  69
  70 private:
  71   //
  72   // parameter specifications
  73   //
  74   Float_t fMIP;                  // MIP position to be in fMIP
  75   Float_t fLowerTrunc;           // lower truncation for dEdx
  76   Float_t fUpperTrunc;           // upper truncation for dEdx
  77   //
  78   Bool_t fUseMax;                 // flag if Qmax or Qtot should be used
  79   //
  80   // histograms
  81   //
  82   TH1F  *fHistNTracks;            //  histogram showing number of ESD tracks per event
  83   TH1F  *fHistClusterShape;       //  histogram to check the cluster shape
  84   TH3F  *fHistQA;                 //  dE/dx histogram showing the final spectrum
  85   //
  86   THnSparseF * fHistGainSector;   //  histogram which shows MIP peak for each of the 3x36 sectors (pad region)
  87   THnSparseF * fHistPadEqual;     //  histogram for the equalization of the gain in the different pad regions -> pass0
  88   THnSparseF * fHistGainMult;     //  histogram which shows decrease of MIP signal as a function
  89   TMatrixD *fPIDMatrix;           //! custom PID matrix
  90   //
  91   THnSparseF * fHistdEdxMap;      // 4D dedx histogram - per sector/phi
  92   THnSparseF * fHistdEdxMax;      // 5D dedx histogram - per 1/dedx, tan(theta), tan(phi), mult, pad-type
  93   THnSparseF * fHistdEdxTot;      // 5D dedx histogram - per 1/dedx, tan(theta), tan(phi), mult, pad-type
  94   TTree *      fdEdxTree;         // tree for the later minimization
  95   TVectorD    *fBBParam;          // BetheBloch parameterization used for the dedx expected calculation
  96   //
  97   AliTPCcalibGainMult(const AliTPCcalibGainMult&);
  98   AliTPCcalibGainMult& operator=(const AliTPCcalibGainMult&);
  99
 100   ClassDef(AliTPCcalibGainMult, 2);
 101 };
 102
 103 #endif
 104
 105