]>
Commit | Line | Data |
---|---|---|
828a6e0c | 1 | #ifndef ALIANALYSISKINKESDat_H |
2 | #define ALIANALYSISKINKESDat_H | |
3 | ||
4 | /* See cxx source for full Copyright notice */ | |
5 | ||
6 | //----------------------------------------------------------------- | |
7 | // AliAnalysisKinkESDat class | |
8 | // This task is an example of an analysis task | |
9 | // for kink topology Study | |
10 | // Authors: Martha Spyropoulou-Stassinaki | |
11 | // and members of the Greek group at the | |
12 | // Physics Department of Athens University | |
13 | // mspyrop@phys.uoa.gr | |
14 | //----------------------------------------------------------------- | |
fcd98e58 | 15 | class AliPIDResponse; |
828a6e0c | 16 | class AliESDVertex; |
17 | class AliESDtrack; | |
18 | class TF1; | |
19 | class TH1F; | |
20 | class TH2F; | |
21 | class TH3F; | |
22 | class TH1D; | |
23 | class TH2D; | |
24 | class TList; | |
25 | class AliESDtrackCuts; | |
26 | ||
27 | #include "AliAnalysisTaskSE.h" | |
28 | ||
29 | class AliAnalysisKinkESDat : public AliAnalysisTaskSE { | |
30 | public: | |
828a6e0c | 31 | AliAnalysisKinkESDat(const char *name = "AliAnalysisKinkESDat"); |
32 | virtual ~AliAnalysisKinkESDat() {} | |
33 | ||
34 | virtual void UserCreateOutputObjects(); | |
35 | virtual void UserExec(Option_t *option); | |
36 | virtual void Terminate(Option_t *); | |
37 | ||
38 | const AliESDVertex *GetEventVertex(const AliESDEvent* esd) const; | |
39 | // Marek multiplicity bins | |
40 | void SetMulCut(Int_t low, Int_t up){fLowMulcut=low;fUpMulcut=up;} | |
b2fb60b8 | 41 | void SetKinkRadius(Float_t lRadiusKLow, Float_t lRadiusKUp) { fKinkRadLow=lRadiusKLow; fKinkRadUp=lRadiusKUp;} |
828a6e0c | 42 | |
5220a87d | 43 | void SetNClusterCut(Int_t lowCluster){fLowCluster=lowCluster;} |
44 | void SetQtCut(Float_t lowQt){fLowQt=lowQt;} | |
828a6e0c | 45 | private: |
46 | TH1F *fHistPtESD; //Pt spectrum of all ESD inside eta, Pt cuts | |
47 | TH1F *fHistPt; //Pt spectrum of all ESD tracks | |
48 | TH1F *fHistQtAll; //Qt spectrum of all kinks | |
49 | TH1F *fHistQt1; //Qt spectrum of Kaon selected sample | |
50 | TH1F *fHistQt2; //Qt spectrum in Qt region of kaons | |
51 | TH1F *fHistPtKaon; //Pt Kaon spectrum of clean sample | |
52 | TH1F *fHistPtKPDG; //Pt Kaon spectrum , confirmed by PDG,inside kaon Qt region | |
53 | TH1F *fHistEta; //Eta spectrum of all kinks | |
54 | TH1F *fHistEtaK; //Eta spectrum of kaons selected by kink topology | |
55 | TH1F *fptKMC; //Pt Kaon spectrum MC, inside eta and pt cuts | |
56 | TH1F *fMultiplMC; //charge multipl MC | |
57 | TH1F *fESDMult; //ESD charged mult | |
58 | TH1F *fgenpt; //Pt Kaon-Kink->mu spectrum , MC, inside eta, Pt, radius cuts | |
59 | TH1F *frad; //radius of kinks, MC , inside the eta nad Pt cuts | |
60 | TH1F *fKinkKaon; //Pt of PDG Kaons inside the selcted ones by the KInk topology | |
61 | TH1F *fKinKRbn; //Pt of PDG Kaons inside the selcted ones by the KInk topology | |
62 | TH1F *fKinkKaonBg; //Pt of the BG inside the kink-Kaon identified spectrum | |
63 | TH1F *fM1kaon; //inv mass of kink-tracks taken as kaons decaying to mu + neutrino | |
828a6e0c | 64 | TH1F *fPtKink; //Pt spectrum of all kinks from track bank |
65 | TH1F *fptKink; //Pt spectrum of all kinks from kink bank | |
828a6e0c | 66 | TH2F *fAngMomK; // Decay angle vrs Mother Mom for pdg kaons |
67 | TH2F *fAngMomPi; // Decay angle vrs Mother Mom for pdg pions | |
68 | TH2F *fAngMomKC; //Decay angle vrs Mother Mom for pdg kaons, inside the selected sample | |
69 | TH1F *fMultESDK; //ESD charged mult | |
70 | TH1F *fMultMCK; //MC K charged mult | |
71 | TH2F *fSignPtNcl;//signPt vrs number of clusters in TPC for kaons from kink sele sample | |
72 | TH2F *fSignPtEta;//signPt vrs Eta in TPC for kaons from kink sele sample | |
73 | TH2F *fEtaNcl;//Eta vrs Nclu in TPC for kaons from kink sele sample | |
74 | TH1F *fSignPt;//signPt in TPC for kaons from kink sele sample | |
75 | TH2F *fChi2NclTPC;//chi2 vrs TPC Nclusters for kaons from kink sele sample | |
76 | TH1F *fRatChi2Ncl;// Ratio chi2/ Ncl TPC for kaons from kink sele sample | |
77 | TH2F *fRadiusNcl;//kink Radius Ncl TPC for kaons from kink sele sample | |
78 | TH2F *fTPCSgnlP;//kink Radius Ncl TPC for kaons from kink sele sample | |
79 | TH2F *fTPCSgnlPa;//kink Radius Ncl TPC for kink sele sample | |
80 | TH1D *fRpr; // Radius of VTX at Y , X plane | |
81 | TH1D *fZpr; //Z distrio of main vertex | |
82 | TH1D *fdcatoVxXY; //dca to Vertex XY distrio | |
83 | TH1D *fnSigmToVx; //nSigma to Vertex distrio of main vertex | |
84 | TH2F *fKinkMothDau; //Mother vrs Daughter | |
85 | TH2F *fZvXv; //two dime of Z vrs X of vtx main | |
86 | TH2F *fZvYv; // two dime of Z vrs Y of vtx main | |
87 | TH2F *fXvYv; // two dime of X vrs Y of main tracks vtx main | |
828a6e0c | 88 | TH1F *fHistPtKaoP; //Pt Kaon spectrum of clean sample pos |
89 | TH1F *fHistPtKaoN; //Pt Kaon spectrum of clean sample neg | |
90 | TH1F *frapiKESD;// rapidi K | |
91 | TH1F *flifetime;//radius of kinks, MC , inside the eta nad Pt cuts | |
92 | TH1F *fradLK;//Length of kinks, MC , inside the eta nad Pt cuts | |
93 | TH3F *fradPtRpDt;//radius of kinks, MC , inside the eta nad Pt cuts | |
94 | TH1F *fInvMuNuAll;//radius of kinks, MC , inside the eta nad Pt cuts | |
95 | TH2F *fQtInvM;// | |
96 | TH1F *fDCAkink;//!MC dcs kink | |
97 | TH2F *fPosiKink;//!MC position kink | |
98 | TH2F *fPosiKinkK;//!MC position kink | |
99 | TH2F *fPosiKinKXZ;//!MC position kink | |
100 | TH2F *fPosiKinKYZ;//!MC position kink | |
101 | TH2F *fPosiKinKBg;//!MC position kink | |
102 | TH2F *fQtMothP;//!qt vrs p mother | |
ad1036b2 | 103 | TH2F *fTPCSgnlPtpc;//Kink mother moment vrs TPC signal |
104 | TH2F *fTPCMomNSgnl;//kink mother TPC momentum vrs nsigmas of dEdx | |
105 | TH2F *fMothKinkMomSgnl;//kink mother TPC momentum vrs nsigmas of dEdx | |
106 | TH1F *fNSigmTPC;//kink mother TPC momentum vrs nsigmas of dEdx | |
107 | TH2F *fTPCSgnlKinkDau;//Kink mother moment vrs TPC signal | |
108 | TH1F *fPtKinkPos; //Pos K Pt spectrum of all kinks from track bank, K0 bins | |
109 | TH1F *fPtKinkNeg; //Neg K Pt spectrum of all kinks from track bank, K0 bins | |
110 | TH2F *fRadNclCln;//kink Radius Ncl TPC for kaons from kink clean sample | |
8e57b720 | 111 | TH1F *fRatioCrossedRows; //ratio crossed rows |
112 | TH1F *fRatioCrossedRowsKink; //ratio crossed rows for kinks | |
113 | TH2F *fRadiusPt;//kinks, Radius vs Pt | |
114 | TH2F *fRadiusPtcln;//kinks, Radius vs Pt for clean kaons | |
115 | TH2F *fInvMassMuNuPt;//kinks,Invariant Mass MuNu vs Pt | |
f9016b0c | 116 | TH2F *fInvMassMuNuPtAll;//kinks,Invariant Mass MuNu vs Pt |
8e57b720 | 117 | TH1F *fPtCut1; //K Pt spectrum of all kinks from track bank, K0 bins |
118 | TH1F *fPtCut2; //K Pt spectrum of all kinks from track bank, K0 bins | |
119 | TH1F *fPtCut3; //K Pt spectrum of all kinks from track bank, K0 bins | |
120 | TH2F *fAngMomKKinks;//kinks, Angle vs Momentum for K-kinks | |
121 | ||
122 | ||
828a6e0c | 123 | |
124 | TF1 *f1; | |
125 | TF1 *f2; | |
126 | TList *fListOfHistos; // list of histos | |
127 | // Marek Multiplicity | |
128 | Int_t fLowMulcut; // | |
129 | Int_t fUpMulcut; | |
b2fb60b8 | 130 | Int_t fKinkRadUp; |
131 | Int_t fKinkRadLow; | |
5220a87d | 132 | Int_t fLowCluster; |
133 | Float_t fLowQt; | |
828a6e0c | 134 | AliESDtrackCuts* fCutsMul; |
8e57b720 | 135 | AliESDtrackCuts* fMaxDCAtoVtxCut; |
fcd98e58 | 136 | AliPIDResponse *fPIDResponse; //! PID response object |
828a6e0c | 137 | |
138 | AliAnalysisKinkESDat(const AliAnalysisKinkESDat&); // not implemented | |
139 | AliAnalysisKinkESDat& operator=(const AliAnalysisKinkESDat&); // not implemented | |
140 | ||
141 | ClassDef(AliAnalysisKinkESDat, 1); // example of analysis | |
142 | }; | |
143 | ||
ad1036b2 | 144 | #endif |