1 #ifndef ALIANALYSISTASKTOFQAID_h
2 #define ALIANALYSISTASKTOFQAID_h
6 class AliAnalysisFilter;
13 #include "AliAnalysisTaskSE.h"
15 class AliAnalysisTaskTOFqaID : public AliAnalysisTaskSE {
18 enum ETrackCutSetTOFqa_t { kRun1Cuts = 0,
25 AliAnalysisTaskTOFqaID();
26 AliAnalysisTaskTOFqaID(const char *name);
27 AliAnalysisTaskTOFqaID(const AliAnalysisTaskTOFqaID& copy);
28 AliAnalysisTaskTOFqaID& operator= (const AliAnalysisTaskTOFqaID& copy);
29 virtual ~AliAnalysisTaskTOFqaID();
31 virtual void UserCreateOutputObjects();
32 virtual void UserExec(Option_t *option);
33 virtual void Terminate(Option_t *);
35 Int_t GetStripIndex(const Int_t * in);
36 void SetTrackFilter(AliAnalysisFilter *filter) {fTrackFilter = filter;};
37 void EnableAdvancedCheck(Bool_t enable){fEnableAdvancedCheck=enable;};
38 void SetExpTimeHistoRange(Float_t min, Float_t max){fExpTimeRangeMin=min; fExpTimeRangeMax=max;return;};
39 void SetExpTimeHistoSmallRange(Float_t min, Float_t max){fExpTimeSmallRangeMin=min; fExpTimeSmallRangeMax=max;return;};
40 void SetExpTimeBinWidth(Float_t width){fExpTimeBinWidth=width;return;};
41 Bool_t SetSelectMCspecies(Bool_t enableMC, Int_t absPdgCode){fIsMC=enableMC; fSelectedPdg=absPdgCode; return kTRUE;};
42 TString GetSpeciesName(Int_t absPdgCode);
43 void HistogramMakeUp(TH1* hist, Color_t color, Int_t markerStyle, TString drawOpt, TString newName, TString newTitle, TString xTitle, TString yTitle);
44 Double_t GetPhiAtTPCouterRadius(AliESDtrack * track);
47 void AddTofBaseHisto(TList *list, Int_t charge, TString suffix);
48 void AddMatchingEffHisto(TList *list, Int_t charge, TString suffix);
49 void AddPidHisto(TList *list, Int_t charge, TString suffix);
50 void AddStartTimeHisto(TList *list, TString suffix);
52 void AddTofTrgHisto(TString suffix);
54 void FillTofBaseHisto(AliESDtrack * track, Int_t charge, TString suffix);
55 void FillMatchedTrkHisto(Int_t charge, TString suffix);
56 void FillPrimaryTrkHisto(Int_t charge, TString suffix);
57 void FillPidHisto(AliESDtrack * track,Int_t charge, TString suffix);
58 void FillStartTimeHisto(TString suffix);
59 void FillTrdHisto(AliESDtrack * track, Int_t charge);
60 void FillStartTimeMaskHisto(TString suffix);
61 void FillTofTrgHisto(TString suffix);
63 Bool_t ComputeTimeZeroByTOF1GeV();
64 Bool_t SelectMCspecies(AliMCEvent * ev, AliESDtrack * track);
65 Bool_t ComputeMatchingEfficiency(TList* list, TString variable);
66 Bool_t IsTPCTOFMatched(AliESDtrack * track);
67 Bool_t IsInTRD(AliESDtrack * track);
68 Bool_t IsEventSelected(AliESDEvent * event);
71 Int_t fRunNumber; //run number
72 AliESDEvent * fESD; //ESD object
73 AliMCEvent * fMCevent; //MC event
74 AliAnalysisFilter * fTrackFilter; //track filter object
75 AliESDVertex * fVertex; //pointer to the vertex object
76 AliESDpid * fESDpid; //pointer to the PID object
77 //AliTOFT0v1 * fTOFT0v1; // TOF-T0 v1
78 AliTOFHeader * fTOFHeader; // TOF header needed for trigger info
79 Int_t fNTOFtracks[3]; //number of tracks matching with TOF
80 //Int_t fNPrimaryTracks; //number of primary tracks
81 Float_t fT0[3]; //event time
82 Float_t fSigmaSpecie[5]; //number of TOF PID sigmas, ie.fSigmaPion, fSigmaKaon, fSigmaProton;
83 Double_t fTrkExpTimes[5]; //expected times from tracking for 5 mass hypothesis
84 Double_t fThExpTimes[5]; //theoretical expected times for 5 mass hypothesis
85 Bool_t fEnableAdvancedCheck; //flag to enable advanced checks
86 Float_t fExpTimeBinWidth;//bin width for t-texp histos
87 Float_t fExpTimeRangeMin, fExpTimeRangeMax; //range of t-texp histogram
88 Float_t fExpTimeSmallRangeMin, fExpTimeSmallRangeMax; //reduced range of t-texp histogram
90 Int_t fnExpTimeSmallBins;
92 Double_t fMyTimeZeroTOF, fMyTimeZeroTOFsigma; //timeZero by TOF recomputed
93 Int_t fMyTimeZeroTOFtracks; // number of tracks used to recompute TOF_T0
94 Bool_t fIsMC; //flag for MC
95 Int_t fSelectedPdg; //pdg code of the selected specie (for MC only)
96 Double_t fP; //momentum
97 Double_t fPt; //transverse momentum
98 Double_t fEta; //psedorapidity
99 Double_t fPhi; //phi at vertex
100 Double_t fTPCOuterPhi; //phi at outer tpc radius
101 Double_t fL; //integrated track lenght
102 Double_t fMatchingMomCut;//pT cut for eta, phi matching eff
105 TList * fHlist; //list of general histos
106 TList * fHlistTimeZero; //list of timeZero related histos
107 TList * fHlistPID; //list of PID-related histos
108 TList * fHlistTRD; //list of general histos for positive tracks
109 TList * fHlistTrigger; //list of general histos for TOF trg infos
111 ClassDef(AliAnalysisTaskTOFqaID, 1); // example of analysis