1 #ifndef ALIANALYSISKINKESDat_H
2 #define ALIANALYSISKINKESDat_H
4 /* See cxx source for full Copyright notice */
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 //-----------------------------------------------------------------
25 class AliESDtrackCuts;
27 #include "AliAnalysisTaskSE.h"
29 class AliAnalysisKinkESDat : public AliAnalysisTaskSE {
31 AliAnalysisKinkESDat(const char *name = "AliAnalysisKinkESDat");
32 virtual ~AliAnalysisKinkESDat() {}
34 virtual void UserCreateOutputObjects();
35 virtual void UserExec(Option_t *option);
36 virtual void Terminate(Option_t *);
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;}
41 void SetKinkRadius(Float_t lRadiusKLow, Float_t lRadiusKUp) { fKinkRadLow=lRadiusKLow; fKinkRadUp=lRadiusKUp;}
43 void SetNClusterCut(Int_t lowCluster){fLowCluster=lowCluster;}
44 void SetQtCut(Float_t lowQt){fLowQt=lowQt;}
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
64 TH1F *fPtKink; //Pt spectrum of all kinks from track bank
65 TH1F *fptKink; //Pt spectrum of all kinks from kink bank
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
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
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
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
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
116 TH2F *fInvMassMuNuPtAll;//kinks,Invariant Mass MuNu vs Pt
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
126 TList *fListOfHistos; // list of histos
127 // Marek Multiplicity
134 AliESDtrackCuts* fCutsMul;
135 AliESDtrackCuts* fMaxDCAtoVtxCut;
136 AliPIDResponse *fPIDResponse; //! PID response object
138 AliAnalysisKinkESDat(const AliAnalysisKinkESDat&); // not implemented
139 AliAnalysisKinkESDat& operator=(const AliAnalysisKinkESDat&); // not implemented
141 ClassDef(AliAnalysisKinkESDat, 1); // example of analysis