1 #ifndef ALITPCCALIBTIMEGAIN_H
2 #define ALITPCCALIBTIMEGAIN_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"
10 #include "THnSparse.h"
12 #include "TObjArray.h"
13 #include "AliSplineFit.h"
22 class AliTPCcalibLaser;
25 #include "TTreeStream.h"
28 class AliTPCcalibTimeGain:public AliTPCcalibBase {
30 AliTPCcalibTimeGain();
31 AliTPCcalibTimeGain(const Text_t *name, const Text_t *title, UInt_t StartTime, UInt_t EndTime, Int_t deltaIntegrationTimeGain);
32 virtual ~AliTPCcalibTimeGain();
34 virtual void Process(AliESDEvent *event);
35 virtual Long64_t Merge(TCollection *li);
36 virtual void AnalyzeRun(Int_t minEntries);
38 void ProcessCosmicEvent(AliESDEvent *event);
39 void ProcessBeamEvent(AliESDEvent *event);
41 void CalculateBetheAlephParams(TH2F *hist, Double_t * ini);
42 static void BinLogX(THnSparse *h, Int_t axisDim);
43 static void BinLogX(TH1 *h);
45 THnSparse * GetHistGainTime() const {return (THnSparse*) fHistGainTime;};
46 TH2F * GetHistDeDxTotal() const {return (TH2F*) fHistDeDxTotal;};
48 TGraphErrors * GetGraphGainVsTime(Int_t runNumber = 0, Int_t minEntries = 2000);
49 static AliSplineFit * MakeSplineFit(TGraphErrors * graph);
50 TGraphErrors * GetGraphAttachment(Int_t minEntries, Int_t nmaxBin, Float_t fracLow=0.1, Float_t fracUp=0.9);
52 void SetMIP(Float_t MIP){fMIP = MIP;};
53 void SetUseMax(Bool_t UseMax){fUseMax = UseMax;};
54 void SetLowerTrunc(Float_t LowerTrunc){fLowerTrunc = LowerTrunc;};
55 void SetUpperTrunc(Float_t UpperTrunc){fUpperTrunc = UpperTrunc;};
56 void SetUseShapeNorm(Bool_t UseShapeNorm){fUseShapeNorm = UseShapeNorm;};
57 void SetUsePosNorm(Bool_t UsePosNorm){fUsePosNorm = UsePosNorm;};
58 void SetUsePadNorm(Int_t UsePadNorm){fUsePadNorm = UsePadNorm;};
59 void SetIsCosmic(Bool_t IsCosmic){fIsCosmic = IsCosmic;};
60 void SetLowMemoryConsumption(Bool_t LowMemoryConsumption){fLowMemoryConsumption = LowMemoryConsumption;};
61 void SetUseCookAnalytical(Bool_t UseCookAnalytical){fUseCookAnalytical = UseCookAnalytical;};
63 static void SetMergeEntriesCut(Double_t entriesCut){fgMergeEntriesCut = entriesCut;}
65 Double_t GetEntries() const {return fHistGainTime->GetEntries();}
68 static Double_t fgMergeEntriesCut; //maximal number of entries for merging -can be modified via setter
70 Float_t GetTPCdEdx(AliTPCseed * seed); // wrapper for CookdEdxNorm or analytical
72 THnSparse * fHistGainTime; // dEdx vs. time, type, Driftlength, momentum P
73 TGraphErrors * fGainVsTime; // multiplication factor vs. time
74 TH2F * fHistDeDxTotal; // dEdx vs. momentum for quality assurance
76 Float_t fIntegrationTimeDeDx; // required statistics for each dEdx time bin
78 Float_t fMIP; // rough MIP position in order to have scaleable histograms
80 Bool_t fUseMax; // true: use max charge for dE/dx calculation, false: use total charge for dE/dx calculation
81 Float_t fLowerTrunc; // lower truncation of dE/dx ; at most 5%
82 Float_t fUpperTrunc; // upper truncation of dE/dx ; ca. 70%
83 Bool_t fUseShapeNorm; // use empirical correction of dependencies
84 Bool_t fUsePosNorm; // charge correction (analytical?)
85 Int_t fUsePadNorm; // normalization of pad geometries
86 Bool_t fUseCookAnalytical; // true if CookdEdxAnalytical should be used
88 Bool_t fIsCosmic; // kTRUE if the analyzed runs contain cosmic events
89 Bool_t fLowMemoryConsumption; // set this option kTRUE if the momenta information should not be stored in order to save memory
91 AliTPCcalibTimeGain(const AliTPCcalibTimeGain&);
92 AliTPCcalibTimeGain& operator=(const AliTPCcalibTimeGain&);
93 void Process(AliESDtrack *track, Int_t runNo=-1){AliTPCcalibBase::Process(track,runNo);};
94 void Process(AliTPCseed *track){return AliTPCcalibBase::Process(track);}
96 ClassDef(AliTPCcalibTimeGain, 1);