]>
Commit | Line | Data |
---|---|---|
0e743ea7 | 1 | #ifndef ALIANALYSISTPCTOFPA_H |
2 | #define ALIANALYSISTPCTOFPA_H | |
3 | ||
4 | ////////////////////////////////////////////////////////////////////////////// | |
5 | // // | |
6 | // This analysis extracts pT-spectra of charged kaons, protons, and pions. // | |
7 | // It is based on particles identifation via the dE/dx signal of the TPC. // | |
8 | // // | |
9 | ////////////////////////////////////////////////////////////////////////////// | |
10 | ||
11 | class TH1; | |
12 | class TH1F; | |
13 | class TH2D; | |
14 | class TH3D; | |
15 | class TList; | |
16 | class TObjArray; | |
17 | class AliESDEvent; | |
18 | class AliESDtrack; | |
19 | class AliESDtrackCuts; | |
20 | class AliHeader; | |
21 | class AliESDpid; | |
22 | ||
23 | ||
24 | #include "AliAnalysisTaskSE.h" | |
25 | #include "THnSparse.h" | |
26 | ||
27 | class AliAnalysisTPCTOFpA : public AliAnalysisTaskSE { | |
28 | public: | |
29 | AliAnalysisTPCTOFpA(const char *name); | |
30 | AliAnalysisTPCTOFpA(); | |
31 | virtual ~AliAnalysisTPCTOFpA() {} | |
32 | // | |
33 | virtual void UserCreateOutputObjects(); | |
34 | virtual void UserExec(Option_t *option); | |
35 | virtual void Terminate(Option_t *); | |
36 | // | |
37 | Bool_t SelectOnImpPar(AliESDtrack* t); | |
38 | // | |
39 | void SetESDtrackCuts(AliESDtrackCuts * trackCuts){fESDtrackCuts = trackCuts;}; | |
40 | void SetAlephParameters(const Double_t * parameters){for(Int_t j=0;j<5;j++) fAlephParameters[j] = parameters[j]; Initialize();}; | |
41 | void SetIsMCtrue(Bool_t isMCdata = kTRUE){fMCtrue = isMCdata;}; | |
42 | void SetUseHBTmultiplicity(Bool_t useHBTmultiplicity = kTRUE){fUseHBTmultiplicity = useHBTmultiplicity;}; | |
43 | void SetUseTPConlyTracks(Bool_t useTPConlyTracks = kTRUE){fUseTPConlyTracks = useTPConlyTracks;}; | |
44 | void SetSaveMotherPDG(Bool_t saveMotherPDG =kTRUE){fSaveMotherPDG = saveMotherPDG;}; | |
45 | void SetSmallTHnSparse(Bool_t smallTHnSparse = kFALSE) {fSmallTHnSparse = smallTHnSparse;}; | |
46 | void SetTPCnSigmaCuts(Double_t nSigmaTPCLow = -3., Double_t nSigmaTPCHigh = 3.){fTPCnSigmaCutLow = nSigmaTPCLow; fTPCnSigmaCutHigh = nSigmaTPCHigh;}; | |
47 | void SetRapidityCuts(Double_t rapidityLow = -0.2, Double_t rapidityHigh = 0.2){fRapidityCutLow = rapidityLow; fRapidityCutHigh = rapidityHigh;}; | |
48 | void SetEvenDCAbinning(Bool_t EvenDCAbinning = kTRUE) {fEvenDCAbinning = EvenDCAbinning;}; | |
49 | void SetIspA(Bool_t ispA = kFALSE) {fIspA = ispA;}; | |
50 | void Initialize(); | |
51 | // | |
52 | ||
53 | private: | |
54 | // | |
55 | void BinLogAxis(const TH1 *h); | |
56 | Int_t GetPythiaEventProcessType(const AliHeader* aHeader, const Bool_t adebug = kFALSE) const; | |
57 | // | |
58 | AliESDEvent *fESD; //! ESD object | |
59 | TList *fListHist; //! list for histograms | |
60 | // | |
61 | AliESDtrackCuts * fESDtrackCuts; // basic cut variables | |
62 | AliESDtrackCuts * fESDTrackCutsMult; // cuts for the MULTIPLICITY DETERMINATION | |
63 | AliESDpid * fESDpid; // basic TPC object for n-sigma cuts | |
64 | Bool_t fMCtrue; // flag if real data or MC is processed | |
65 | Bool_t fOnlyQA; // flag if only QA histograms should be filled | |
66 | Bool_t fUseHBTmultiplicity; // flag if multiplicity determination should be done as in the HBT paper | |
67 | Bool_t fUseTPConlyTracks; // flag if TPConly-track should be used | |
68 | Bool_t fSaveMotherPDG; // flag if PDG of mother should be saved (weak decays) | |
69 | Bool_t fSmallTHnSparse; // flag if to do cuts on TPC n-sigma and rapidity in task or not | |
70 | Bool_t fIspA; // flag for pA analysis | |
71 | Double_t fTPCnSigmaCutLow; // low border for TPC n-sigma cut | |
72 | Double_t fTPCnSigmaCutHigh; // high border for TPC n-sigma cut | |
73 | Double_t fRapidityCutLow; // low border for rapidity cut | |
74 | Double_t fRapidityCutHigh; // high border for rapidity cut | |
75 | Double_t fEvenDCAbinning; // same bin width for all bins in DCA xy | |
76 | Double_t fAlephParameters[5]; // Aleph Parameters for Bethe-Bloch | |
77 | ||
78 | // | |
79 | // | |
80 | // | |
81 | THnSparseF * fHistRealTracks; //! histogram with all necessary information for real tracks | |
82 | THnSparseF * fHistMCparticles; //! histogram with all necessary information for MC particles | |
83 | // | |
84 | TH3D * fHistPidQA; //! histogram for the QA of the PID | |
85 | TH2D * fHistMult; //! control histogram for multiplicity | |
86 | TH1D * fHistCentrality; //! control histogram for centrality | |
87 | // | |
88 | AliAnalysisTPCTOFpA(const AliAnalysisTPCTOFpA&); | |
89 | AliAnalysisTPCTOFpA& operator=(const AliAnalysisTPCTOFpA&); | |
90 | ||
91 | ClassDef(AliAnalysisTPCTOFpA, 1); | |
92 | }; | |
93 | ||
94 | #endif |