1 /**************************************************************************
2 * This file is property of and copyright by the ALICE HLT Project *
3 * ALICE Experiment at CERN, All rights reserved. *
5 * Primary Author: Svein Lindal <slindal@fys.uio.no> *
7 * Permission to use, copy, modify and distribute this software and its *
8 * documentation strictly for non-commercial purposes is hereby granted *
9 * without fee, provided that the above copyright notice appears in all *
10 * copies and that both the copyright notice and this permission notice *
11 * appear in the supporting documentation. The authors make no claims *
12 * about the suitability of this software for any purpose. It is *
13 * provided "as is" without express or implied warranty. *
14 **************************************************************************/
16 /// @file AliAnalysisTaskGammaJet.cxx
17 /// @author Svein Lindal
18 /// @brief Class used to run isolation studies of conversion gamma/pions
21 #include "AliAnaConvIsolation.h"
22 #include "AliAODTrack.h"
25 #include "TClonesArray.h"
28 #include "AliAODConversionPhoton.h"
33 // Gamma - jet correlation analysis task
34 // Author: Svein Lindal
38 ClassImp(AliAnaConvIsolation)
40 //________________________________________________________________________
41 AliAnaConvIsolation::AliAnaConvIsolation () : TObject(),
43 fCurveFunction("0.1*x"),
51 // fHistogramMaxPt(50)
54 // for(Int_t i = 0; i < 2; i++){
55 // fhMaxPtInCone[i] = NULL;
56 // fhSumPtInCone[i] = NULL;
57 // fhSumPtVsMaxPt[i] = NULL;
58 // fhPtCandidates[i] = NULL;
59 // fhTrackMult[i] = NULL;
65 //________________________________________________________________________
66 AliAnaConvIsolation::AliAnaConvIsolation(Float_t coneSize, Float_t maxPtThreshold, Float_t sumPtThreshold, Float_t maxPtFraction, Float_t sumPtFraction) :
69 fCurveFunction("0.1*x"),
72 fMaxPtThreshold(maxPtThreshold),
73 fSumPtThreshold(sumPtThreshold),
74 fMaxPtFraction(maxPtFraction),
75 fSumPtFraction(sumPtFraction)
77 // fHistogramMaxPt(50)
80 // for(Int_t i = 0; i < 2; i++){
81 // fhMaxPtInCone[i] = NULL;
82 // fhSumPtInCone[i] = NULL;
83 // fhSumPtVsMaxPt[i] = NULL;
84 // fhPtCandidates[i] = NULL;
85 // fhTrackMult[i] = NULL;
88 fIsoCurve = new TF1("Isolation_curve", fCurveFunction.Data(), 0, 100);
93 //________________________________________________________________________________
94 AliAnaConvIsolation::~AliAnaConvIsolation() {
100 //________________________________________________________________________
101 // void AliAnaConvIsolation::CreateHistograms()
104 // if(!fHistograms) fHistograms = new TList();
105 // fHistograms->SetName(Form("Isolation_histo_cone_%f_maxPt_%f_sumPt_%f", fConeSize, fSumPtThreshold, fMaxPtThreshold));
107 //fHistograms->Add(fIsoCurve);
109 //cout << "Creatin isolation histograms conesize :" << fConeSize << endl;
111 //Create histograms add, to outputlis
112 // for(Int_t i = 0; i < 2; i++){
113 // fhMaxPtInCone[i] = new TH2F(Form("fhMaxPtInCone_%s_%f", (i==0)? "nonIso" : "isolated", fConeSize),
114 // Form("Max pt nonIso particle in cone %f vs candidate pt", fConeSize),
115 // 200, 0, fHistogramMaxPt, 200, 0, fHistogramMaxPt);
116 // fHistograms->Add(fhMaxPtInCone[i]);
118 // fhSumPtInCone[i] = new TH2F(Form("fhSumPtInCone_%s_%f", (i==0)? "nonIso" : "isolated", fConeSize),
119 // Form("Sum pt in cone %f vs candidate pt %s", fConeSize, (i==0)? "nonIsoay" : "isolated"),
120 // 200, 0, fHistogramMaxPt, 200, 0, fHistogramMaxPt);
121 // fHistograms->Add(fhSumPtInCone[i]);
123 // fhSumPtVsMaxPt[i] = new TH2F(Form("fhSumPtVsMaxPt_%s_%f", (i==0)? "nonIso" : "isolated", fConeSize),
124 // Form("fhSumPtVsMaxPt_%s_%f", (i==0)? "nonIso" : "isolated", fConeSize),
125 // 200, 0, fHistogramMaxPt, 200, 0, fHistogramMaxPt);
126 // fHistograms->Add(fhSumPtVsMaxPt[i]);
129 // for(Int_t iIso = 0; iIso < 2; iIso++){
130 // fhPtCandidates[iIso] = new TH1F(Form("fhPtCandidates_%s_%f", (iIso==0)? "nonIso" : "isolated", fConeSize),
131 // Form("Pt of %s candidates, cone size %f", (iIso==0)? "nonIsolated" : "isolated", fConeSize),
132 // 200, 0 , fHistogramMaxPt);
133 // fhPtCandidates[iIso]->GetXaxis()->SetTitle("P_{T} (GeV/c)");
134 // fhPtCandidates[iIso]->GetYaxis()->SetTitle("dN/dP_{T} (c/GeV)");
135 // fhPtCandidates[iIso]->SetMarkerStyle(kFullCircle);
136 // fHistograms->Add(fhPtCandidates[iIso]);
139 // for(int iDec = 0; iDec < 2; iDec++) {
140 // for(int iIso = 0; iIso < 2; iIso++) {
142 // fHistSumPt[iIso][iDec] = new TH1F(Form("fHistSumPt_%f_%s_%s", fConeSize, (iIso==0)?"nonIso":"isolated" , (iDec==0)?"noDec":"decay" ),
143 // Form("P_{T} distribution cone %f %s %s", fConeSize, (iIso==0)?"nonIso":"isolated" , (iDec==0)?"noDec":"decay" ),
145 // fHistSumPt[iIso][iDec]->GetXaxis()->SetTitle("P_{T} (GeV/c)");
146 // fHistSumPt[iIso][iDec]->GetYaxis()->SetTitle("dN/dP_{T} (c/GeV)");
147 // fHistSumPt[iIso][iDec]->SetMarkerStyle(kFullCircle);
148 // fHistograms->Add(fHistSumPt[iIso][iDec]);
151 // fHistMaxPt[iIso][iDec] = new TH1F(Form("fHistMaxPt_%f_%s_%s", fConeSize, (iIso==0)?"nonIso":"isolated" , (iDec==0)?"noDec":"decay" ),
152 // Form("P_{T} distribution cone %f %s %s", fConeSize, (iIso==0)?"nonIso":"isolated" , (iDec==0)?"noDec":"decay" ),
154 // fHistMaxPt[iIso][iDec]->GetXaxis()->SetTitle("P_{T} (GeV/c)");
155 // fHistMaxPt[iIso][iDec]->GetYaxis()->SetTitle("dN/dP_{T} (c/GeV)");
156 // fHistMaxPt[iIso][iDec]->SetMarkerStyle(kFullCircle);
157 // fHistograms->Add(fHistMaxPt[iIso][iDec]);
164 // for(Int_t iIso = 0; iIso < 2; iIso++){
165 // fhTrackMult[iIso] = new TH1F(Form("fhTrackMult_%s_%f", (iIso==0)? "nonIso" : "isolated", fConeSize),
166 // Form("Pt of %s candidates, cone size %f", (iIso==0)? "nonIsolated" : "isolated", fConeSize),
168 // fhTrackMult[iIso]->GetXaxis()->SetTitle("n tracks in event");
169 // fhTrackMult[iIso]->GetYaxis()->SetTitle("dN/dNTracks");
170 // fhTrackMult[iIso]->SetMarkerStyle(kFullCircle);
171 // fHistograms->Add(fhTrackMult[iIso]);
176 //________________________________________________________________________
177 Bool_t AliAnaConvIsolation::IsLeading(AliAODConversionParticle * particle, const TObjArray * tracks, const TObjArray * aodParticles) {
178 //Check if particle is highest pt particle in cone
179 for(Int_t ip = 0; ip < aodParticles->GetEntriesFast(); ip++) {
180 AliAODConversionParticle * oparticle = static_cast<AliAODConversionParticle*>(aodParticles->UncheckedAt(ip));
181 if(oparticle == particle) continue;
182 if(oparticle->Pt() > particle->Pt() &&
183 IsInCone(particle->Eta() - oparticle->Eta(),
184 particle->Phi() - oparticle->Phi(),
190 for(Int_t it = 0; it < tracks->GetEntriesFast(); it++) {
191 AliAODTrack * track = static_cast<AliAODTrack*>(tracks->UncheckedAt(it));
192 if(track->Pt() > particle->Pt() &&
193 IsInCone(particle->Eta() - track->Eta(), particle->Phi() - track->Phi(), fConeSize) ) return kFALSE;
199 ///________________________________________________________________________
200 Int_t AliAnaConvIsolation::IsLeading(const AliAODConversionParticle * particle, const TObjArray * tracks, const Int_t tIDs[4]) const {
202 Bool_t leadingEvent = kTRUE;
204 //Is there a higher pt particle within cone ?
205 Float_t tPt = particle->Pt();
206 for(int it = 0; it < tracks->GetEntriesFast(); it++) {
207 AliVTrack * track = static_cast<AliVTrack*>(tracks->UncheckedAt(it));
208 Int_t tid = track->GetID();
209 if(tid == tIDs[0] || tid == tIDs[1] || tid == tIDs[2] || tid == tIDs[3]) continue;
210 if(track->Pt() > tPt) {
211 if(IsInCone(particle->Eta() - track->Eta(), particle->Phi() - track->Phi(), 0.8)) return 0;
212 leadingEvent = kFALSE;
225 //_________________________________________________________________________
226 Bool_t AliAnaConvIsolation::IsIsolated(AliAODConversionPhoton * particle, const TClonesArray * const tracks, const Int_t nSpawn, const Int_t * const spawn, Bool_t &leading) {
227 //See header file for documentation
234 for(int it = 0; it < tracks->GetEntriesFast(); it++) {
235 AliAODTrack * track = dynamic_cast<AliAODTrack*>(tracks->At(it));
239 if(track->Pt() < fMinPt) continue;
241 if(nSpawn && spawn) { ;}
243 ///Ignore tracks that are grandchildren of pion
244 // if ( particle->IsMySpawn(track->GetID(), nSpawn, spawn))
248 if ( IsInCone(particle->Eta() - track->Eta(), particle->Phi() - track->Phi(), fConeSize) ) {
249 ptSum += track->Pt();
250 if(track->Pt() > ptMax) ptMax = track->Pt();
251 if(track->Pt() > particle->Pt()) leading = kFALSE;
254 cout << "Bad track"<<endl;
259 Bool_t isolated = EvaluateIsolationCriteria( ptSum, particle->Pt());
261 //FillHistograms(particle->Pt(), ptMax, ptSum, isolated, tracks->GetEntriesFast());
268 //_________________________________________________________________________
269 Bool_t AliAnaConvIsolation::IsIsolated(const AliAODConversionPhoton * const particle, const TClonesArray * const tracks, Bool_t &leading ) {
270 //See header file for documentation
277 for(int it = 0; it < tracks->GetEntriesFast(); it++) {
279 AliAODTrack * track = dynamic_cast<AliAODTrack*>(tracks->At(it));
283 if(track->Pt() < fMinPt) continue;
285 if ( (track->GetID() == particle->GetTrackLabel(0)) || track->GetID() == particle->GetTrackLabel(1) )
288 if ( IsInCone(particle->Eta() - track->Eta(), particle->Phi() - track->Phi(), fConeSize) ) {
289 ptSum += track->Pt();
290 if(track->Pt() > ptMax) ptMax = track->Pt();
291 if(track->Pt() > particle->Pt()) leading = kFALSE;
294 cout << "Bad track"<<endl;
299 Bool_t isolated = EvaluateIsolationCriteria( ptSum, particle->Pt());
301 //FillHistograms(particle->Pt(), ptMax, ptSum, isolated, tracks->GetEntriesFast());
308 // ///___________________________________________________________________________
309 // void AliAnaConvIsolation::FillHistograms(Float_t pt, Float_t ptMax, Float_t ptSum, Bool_t isolated, Int_t nTracks) {
311 // fhMaxPtInCone[isolated]->Fill(pt, ptMax);
312 // fhSumPtInCone[isolated]->Fill(pt, ptSum);
313 // fhSumPtVsMaxPt[isolated]->Fill(ptMax, ptSum);
314 // fhPtCandidates[isolated]->Fill(pt);
315 // fhTrackMult[isolated]->Fill(nTracks);
319 ///_____________________________________________________________________________
320 Bool_t AliAnaConvIsolation::EvaluateIsolationCriteria(Float_t ptSum, Float_t pt) const {
321 //Evaluate isolation criteria
322 return (ptSum < fIsoCurve->Eval(pt));