1 #ifndef ALITPCCALIBGAINMULT_H
2 #define ALITPCCALIBGAINMULT_H
4 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
5 * See cxx source for full Copyright notice */
7 #include "AliTPCcalibBase.h"
8 #include "AliTPCCalPad.h"
11 #include "THnSparse.h"
20 #include "TTreeStream.h"
23 class AliTPCcalibGainMult:public AliTPCcalibBase {
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();
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);
41 TH1F * GetHistNTracks() const {return fHistNTracks;};
42 TH1F * GetHistClusterShape() const {return fHistClusterShape;};
43 TH3F * GetHistQA() const {return fHistQA;};
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
54 void SetMIPvalue(Float_t mip){fMIP = mip;};
55 void SetLowerTrunc(Float_t lowerTrunc){fLowerTrunc = lowerTrunc;};
56 void SetUpperTrunc(Float_t upperTrunc){fUpperTrunc = upperTrunc;};
57 void SetUseMax(Bool_t useMax){fUseMax = useMax;};
60 void Process(AliESDtrack *track, Int_t runNo=-1){AliTPCcalibBase::Process(track,runNo);};
61 void Process(AliTPCseed *track){return AliTPCcalibBase::Process(track);}
63 void MakeLookup(THnSparse * hist, Char_t * outputFile);
66 void UpdateClusterParam();
71 // parameter specifications
73 Float_t fMIP; // MIP position to be in fMIP
74 Float_t fLowerTrunc; // lower truncation for dEdx
75 Float_t fUpperTrunc; // upper truncation for dEdx
77 Bool_t fUseMax; // flag if Qmax or Qtot should be used
81 TH1F *fHistNTracks; // histogram showing number of ESD tracks per event
82 TH1F *fHistClusterShape; // histogram to check the cluster shape
83 TH3F *fHistQA; // dE/dx histogram showing the final spectrum
85 THnSparseF * fHistGainSector; // histogram which shows MIP peak for each of the 3x36 sectors (pad region)
86 THnSparseF * fHistPadEqual; // histogram for the equalization of the gain in the different pad regions -> pass0
87 THnSparseF * fHistGainMult; // histogram which shows decrease of MIP signal as a function
88 TMatrixD *fPIDMatrix; //! custom PID matrix
90 THnSparseF * fHistdEdxMap; // 4D dedx histogram - per sector/phi
91 THnSparseF * fHistdEdxMax; // 5D dedx histogram - per 1/dedx, tan(theta), tan(phi), mult, pad-type
92 THnSparseF * fHistdEdxTot; // 5D dedx histogram - per 1/dedx, tan(theta), tan(phi), mult, pad-type
93 TTree * fdEdxTree; // tree for the later minimization
94 TVectorD *fBBParam; // BetheBloch parameterization used for the dedx expected calculation
96 AliTPCcalibGainMult(const AliTPCcalibGainMult&);
97 AliTPCcalibGainMult& operator=(const AliTPCcalibGainMult&);
99 ClassDef(AliTPCcalibGainMult, 2);