#include "TH3F.h"
#include "TF1.h"
#include "THnSparse.h"
+#include "TMatrixD.h"
+#include "TVectorD.h"
class TH1F;
class TList;
class AliESDEvent;
AliTPCcalibGainMult();
AliTPCcalibGainMult(const Text_t *name, const Text_t *title);
virtual ~AliTPCcalibGainMult();
+ void SetBBParam(TVectorD * param) {fBBParam=param;}
+ // virtual void Terminate();
//
virtual void Process(AliESDEvent *event);
+ virtual void ProcessV0s(AliESDEvent *event);
+ virtual void ProcessCosmic(const AliESDEvent *event);
+ virtual void ProcessKinks(const AliESDEvent *event);
+ virtual void ProcessTOF(const AliESDEvent *event);
+ virtual void DumpHPT(const AliESDEvent *event);
virtual Long64_t Merge(TCollection *li);
virtual void Analyze();
+ void DumpTrack(AliESDtrack * track, AliESDfriendTrack *ftrack, AliTPCseed * seed, Int_t index);
//
TH1F * GetHistNTracks() const {return fHistNTracks;};
TH1F * GetHistClusterShape() const {return fHistClusterShape;};
//
THnSparseF * GetHistGainSector() const {return fHistGainSector;};
THnSparseF * GetHistPadEqual() const {return fHistPadEqual;};
- THnSparseF * GetHistGainMult() const {return fHistGainMult;};
+ THnSparseF * GetHistGainMult() const {return fHistGainMult;};
+ THnSparseF * GetHistdEdxMap() const { return fHistdEdxMap;} // 4D dedx histogram
+ THnSparseF * GetHistdEdxMax() const { return fHistdEdxMax;} // 4D dedx histogram
+ THnSparseF * GetHistdEdxTot() const { return fHistdEdxTot;} // 4D dedx histogram
+ TTree * GetdEdxTree() const {return fdEdxTree;} // tree for the later minimization
+
//
void SetMIPvalue(Float_t mip){fMIP = mip;};
void SetLowerTrunc(Float_t lowerTrunc){fLowerTrunc = lowerTrunc;};
void SetUseMax(Bool_t useMax){fUseMax = useMax;};
//
//
- static void BinLogX(const TH1 *h); // method for correct histogram binning
void Process(AliESDtrack *track, Int_t runNo=-1){AliTPCcalibBase::Process(track,runNo);};
void Process(AliTPCseed *track){return AliTPCcalibBase::Process(track);}
//
+ void MakeLookup(THnSparse * hist, Char_t * outputFile);
+ //
void UpdateGainMap();
void UpdateClusterParam();
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
+ TMatrixD *fPIDMatrix; //! custom PID matrix
+ //
+ THnSparseF * fHistdEdxMap; // 4D dedx histogram - per sector/phi
+ THnSparseF * fHistdEdxMax; // 5D dedx histogram - per 1/dedx, tan(theta), tan(phi), mult, pad-type
+ THnSparseF * fHistdEdxTot; // 5D dedx histogram - per 1/dedx, tan(theta), tan(phi), mult, pad-type
+ TTree * fdEdxTree; // tree for the later minimization
+ TVectorD *fBBParam; // BetheBloch parameterization used for the dedx expected calculation
//
AliTPCcalibGainMult(const AliTPCcalibGainMult&);
AliTPCcalibGainMult& operator=(const AliTPCcalibGainMult&);
- ClassDef(AliTPCcalibGainMult, 1);
+ ClassDef(AliTPCcalibGainMult, 2);
};
#endif