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;}
44 TH1F *fHistPtESD; //Pt spectrum of all ESD inside eta, Pt cuts
45 TH1F *fHistPt; //Pt spectrum of all ESD tracks
46 TH1F *fHistQtAll; //Qt spectrum of all kinks
47 TH1F *fHistQt1; //Qt spectrum of Kaon selected sample
48 TH1F *fHistQt2; //Qt spectrum in Qt region of kaons
49 TH1F *fHistPtKaon; //Pt Kaon spectrum of clean sample
50 TH1F *fHistPtKPDG; //Pt Kaon spectrum , confirmed by PDG,inside kaon Qt region
51 TH1F *fHistEta; //Eta spectrum of all kinks
52 TH1F *fHistEtaK; //Eta spectrum of kaons selected by kink topology
53 TH1F *fptKMC; //Pt Kaon spectrum MC, inside eta and pt cuts
54 TH1F *fMultiplMC; //charge multipl MC
55 TH1F *fESDMult; //ESD charged mult
56 TH1F *fgenpt; //Pt Kaon-Kink->mu spectrum , MC, inside eta, Pt, radius cuts
57 TH1F *frad; //radius of kinks, MC , inside the eta nad Pt cuts
58 TH1F *fKinkKaon; //Pt of PDG Kaons inside the selcted ones by the KInk topology
59 TH1F *fKinKRbn; //Pt of PDG Kaons inside the selcted ones by the KInk topology
60 TH1F *fKinkKaonBg; //Pt of the BG inside the kink-Kaon identified spectrum
61 TH1F *fM1kaon; //inv mass of kink-tracks taken as kaons decaying to mu + neutrino
62 TH1F *fPtKink; //Pt spectrum of all kinks from track bank
63 TH1F *fptKink; //Pt spectrum of all kinks from kink bank
64 TH2F *fAngMomK; // Decay angle vrs Mother Mom for pdg kaons
65 TH2F *fAngMomPi; // Decay angle vrs Mother Mom for pdg pions
66 TH2F *fAngMomKC; //Decay angle vrs Mother Mom for pdg kaons, inside the selected sample
67 TH1F *fMultESDK; //ESD charged mult
68 TH1F *fMultMCK; //MC K charged mult
69 TH2F *fSignPtNcl;//signPt vrs number of clusters in TPC for kaons from kink sele sample
70 TH2F *fSignPtEta;//signPt vrs Eta in TPC for kaons from kink sele sample
71 TH2F *fEtaNcl;//Eta vrs Nclu in TPC for kaons from kink sele sample
72 TH1F *fSignPt;//signPt in TPC for kaons from kink sele sample
73 TH2F *fChi2NclTPC;//chi2 vrs TPC Nclusters for kaons from kink sele sample
74 TH1F *fRatChi2Ncl;// Ratio chi2/ Ncl TPC for kaons from kink sele sample
75 TH2F *fRadiusNcl;//kink Radius Ncl TPC for kaons from kink sele sample
76 TH2F *fTPCSgnlP;//kink Radius Ncl TPC for kaons from kink sele sample
77 TH2F *fTPCSgnlPa;//kink Radius Ncl TPC for kink sele sample
78 TH1D *fRpr; // Radius of VTX at Y , X plane
79 TH1D *fZpr; //Z distrio of main vertex
80 TH1D *fdcatoVxXY; //dca to Vertex XY distrio
81 TH1D *fnSigmToVx; //nSigma to Vertex distrio of main vertex
82 TH2F *fKinkMothDau; //Mother vrs Daughter
83 TH2F *fZvXv; //two dime of Z vrs X of vtx main
84 TH2F *fZvYv; // two dime of Z vrs Y of vtx main
85 TH2F *fXvYv; // two dime of X vrs Y of main tracks vtx main
86 TH1F *fHistPtKaoP; //Pt Kaon spectrum of clean sample pos
87 TH1F *fHistPtKaoN; //Pt Kaon spectrum of clean sample neg
88 TH1F *frapiKESD;// rapidi K
89 TH1F *flifetime;//radius of kinks, MC , inside the eta nad Pt cuts
90 TH1F *fradLK;//Length of kinks, MC , inside the eta nad Pt cuts
91 TH3F *fradPtRpDt;//radius of kinks, MC , inside the eta nad Pt cuts
92 TH1F *fInvMuNuAll;//radius of kinks, MC , inside the eta nad Pt cuts
94 TH1F *fDCAkink;//!MC dcs kink
95 TH2F *fPosiKink;//!MC position kink
96 TH2F *fPosiKinkK;//!MC position kink
97 TH2F *fPosiKinKXZ;//!MC position kink
98 TH2F *fPosiKinKYZ;//!MC position kink
99 TH2F *fPosiKinKBg;//!MC position kink
100 TH2F *fQtMothP;//!qt vrs p mother
101 TH2F *fTPCSgnlPtpc;//Kink mother moment vrs TPC signal
102 TH2F *fTPCMomNSgnl;//kink mother TPC momentum vrs nsigmas of dEdx
103 TH2F *fMothKinkMomSgnl;//kink mother TPC momentum vrs nsigmas of dEdx
104 TH1F *fNSigmTPC;//kink mother TPC momentum vrs nsigmas of dEdx
105 TH2F *fTPCSgnlKinkDau;//Kink mother moment vrs TPC signal
106 TH1F *fPtKinkPos; //Pos K Pt spectrum of all kinks from track bank, K0 bins
107 TH1F *fPtKinkNeg; //Neg K Pt spectrum of all kinks from track bank, K0 bins
108 TH2F *fRadNclCln;//kink Radius Ncl TPC for kaons from kink clean sample
109 TH1F *fRatioCrossedRows; //ratio crossed rows
110 TH1F *fRatioCrossedRowsKink; //ratio crossed rows for kinks
111 TH2F *fRadiusPt;//kinks, Radius vs Pt
112 TH2F *fRadiusPtcln;//kinks, Radius vs Pt for clean kaons
113 TH2F *fInvMassMuNuPt;//kinks,Invariant Mass MuNu vs Pt
114 TH1F *fPtCut1; //K Pt spectrum of all kinks from track bank, K0 bins
115 TH1F *fPtCut2; //K Pt spectrum of all kinks from track bank, K0 bins
116 TH1F *fPtCut3; //K Pt spectrum of all kinks from track bank, K0 bins
117 TH2F *fAngMomKKinks;//kinks, Angle vs Momentum for K-kinks
123 TList *fListOfHistos; // list of histos
124 // Marek Multiplicity
129 AliESDtrackCuts* fCutsMul;
130 AliESDtrackCuts* fMaxDCAtoVtxCut;
131 AliPIDResponse *fPIDResponse; //! PID response object
133 AliAnalysisKinkESDat(const AliAnalysisKinkESDat&); // not implemented
134 AliAnalysisKinkESDat& operator=(const AliAnalysisKinkESDat&); // not implemented
136 ClassDef(AliAnalysisKinkESDat, 1); // example of analysis