1 /**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * Authors: Svein Lindal, Daniel Lohner *
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 ////////////////////////////////////////////////
17 //---------------------------------------------
18 // QA Task for V0 Reader V1
19 //---------------------------------------------
20 ////////////////////////////////////////////////
22 #include "AliAnalysisTaskConversionQA.h"
24 #include "AliAnalysisManager.h"
25 #include "TParticle.h"
27 #include "AliPIDResponse.h"
29 #include "AliESDtrackCuts.h"
35 ClassImp(AliAnalysisTaskConversionQA)
37 //________________________________________________________________________
38 AliAnalysisTaskConversionQA::AliAnalysisTaskConversionQA(const char *name) : AliAnalysisTaskSE(name),
40 fConversionGammas(NULL),
41 fConversionCuts(NULL),
43 fNumberOfESDTracks(0),
49 ffillHistograms(kFALSE),
53 hNGoodESDTracks(NULL),
55 hNContributorsVertex(NULL),
60 hGammaChi2perNDF(NULL),
63 hGammaCosinePointingAngle(NULL),
69 hElecNfindableClsTPC(NULL),
70 hPosiNfindableClsTPC(NULL),
72 hTrueResoulutionR(NULL),
73 hTrueResoulutionZ(NULL),
74 hTrueResoulutionPhi(NULL),
79 hTrueGammaChi2perNDF(NULL),
80 hTrueGammaPsiPair(NULL),
82 hTrueGammaCosinePointingAngle(NULL),
88 hTrueElecNfindableClsTPC(NULL),
89 hTruePosiNfindableClsTPC(NULL)
91 // Default constructor
93 DefineInput(0, TChain::Class());
94 DefineOutput(1, TList::Class());
97 //________________________________________________________________________
98 AliAnalysisTaskConversionQA::~AliAnalysisTaskConversionQA()
100 // default deconstructor
106 //________________________________________________________________________
107 void AliAnalysisTaskConversionQA::UserCreateOutputObjects()
109 // Create User Output Objects
111 if(fOutputList != NULL){
115 if(fOutputList == NULL){
116 fOutputList = new TList();
117 fOutputList->SetOwner(kTRUE);
122 TH1::SetDefaultSumw2(kTRUE);
124 fESDList = new TList();
125 fESDList->SetOwner(kTRUE);
126 fESDList->SetName("ESD QA");
127 fOutputList->Add(fESDList);
129 hVertexZ = new TH1F("Vertex_Z","Vertex_Z",300,-15,15);
130 fESDList->Add(hVertexZ);
131 hNContributorsVertex = new TH1I("ContrVertex_Z","ContrVertex_Z",3000,0,3000);
132 fESDList->Add(hNContributorsVertex);
133 if(fIsHeavyIon) hNGoodESDTracks = new TH1I("GoodESDTracks","GoodESDTracks",3000,0,3000);
134 else hNGoodESDTracks = new TH1I("GoodESDTracks","GoodESDTracks",200,0,200);
135 fESDList->Add(hNGoodESDTracks);
136 if(fIsHeavyIon) hNV0Tracks = new TH1I("V0 Multiplicity","V0 Multiplicity",30000,0,30000);
137 else hNV0Tracks = new TH1I("V0 Multiplicity","V0 Multiplicity",2000,0,2000);
138 fESDList->Add(hNV0Tracks);
140 hITSClusterPhi = new TH2F("ITSClusterPhi","hITSClusterPhi",72,0,2*TMath::Pi(),7,0,7);
141 fESDList->Add(hITSClusterPhi);
142 hGammaPt = new TH1F("Gamma_Pt","Gamma_Pt",250,0,25);
143 fESDList->Add(hGammaPt);
144 hGammaPhi = new TH1F("Gamma_Phi","Gamma_Phi",360,0,2*TMath::Pi());
145 fESDList->Add(hGammaPhi);
146 hGammaEta = new TH1F("Gamma_Eta","Gamma_Eta",400,-1.2,1.2);
147 fESDList->Add(hGammaEta);
148 hGammaChi2perNDF = new TH1F("Gamma_Chi2perNDF","Gamma_Chi2perNDF",500,0,100);
149 fESDList->Add(hGammaChi2perNDF);
150 hGammaPsiPair = new TH1F("Gamma_PsiPair","Gamma_PsiPair",500,0,2);
151 fESDList->Add(hGammaPsiPair);
152 hGammaQt = new TH1F("Gamma_Qt","Gamma_Qt",400,0,0.1);
153 fESDList->Add(hGammaQt);
154 hGammaCosinePointingAngle = new TH1F("Gamma_CosinePointingAngle","Gamma_CosinePointingAngle",900,0.7,1.);
155 fESDList->Add(hGammaCosinePointingAngle);
156 hGammaXY = new TH2F("Gamma_ConversionPoint_XY","Gamma_ConversionPoint_XY",960,-120,120,960,-120,120);
157 fESDList->Add(hGammaXY);
158 hGammaZR= new TH2F("Gamma_ConversionPoint_ZR","Gamma_ConversionPoint_ZR",1200,-150,150,480,0,120);
159 fESDList->Add(hGammaZR);
161 hElecPt = new TH2F("Electron_Positron_Pt","Electron_Positron_Pt",250,0,25,250,0,25);
162 fESDList->Add(hElecPt);
163 hElecEta = new TH2F("Electron_Positron_Eta","Electron_Positron_Eta",400,-1.2,1.2,400,-1.2,1.2);
164 fESDList->Add(hElecEta);
165 hElecPhi = new TH2F("Electron_Positron_Phi","Electron_Positron_Phi",360,0,2*TMath::Pi(),360,0,2*TMath::Pi());
166 fESDList->Add(hElecPhi);
167 hElecNfindableClsTPC = new TH1F("Electron_findableClusterTPC","Electron_findableClusterTPC",100,0,1);
168 fESDList->Add(hElecNfindableClsTPC);
169 hPosiNfindableClsTPC = new TH1F("Positron_findableClusterTPC","Positron_findableClusterTPC",100,0,1);
170 fESDList->Add(hPosiNfindableClsTPC);
173 fTrueList = new TList();
174 fTrueList->SetOwner(kTRUE);
175 fTrueList->SetName("True QA");
176 fOutputList->Add(fTrueList);
178 hTrueResoulutionR = new TH2F("True_ConversionPointResolution_R","True_ConversionPointResolution_R",240,0,120,200,-20,20);
179 fTrueList->Add(hTrueResoulutionR);
180 hTrueResoulutionZ = new TH2F("True_ConversionPointResolution_Z","True_ConversionPointResolution_Z",480,-120,120,200,-20,20);
181 fTrueList->Add(hTrueResoulutionZ);
182 hTrueResoulutionPhi = new TH2F("True_ConversionPointResolution_Phi","True_ConversionPointResolution_Phi",360,0,2*TMath::Pi(),200,-20,20);
183 fTrueList->Add(hTrueResoulutionPhi);
185 hTrueGammaPt = new TH1F("True_Gamma_Pt","True_Gamma_Pt",250,0,25);
186 fTrueList->Add(hTrueGammaPt);
187 hTrueGammaPhi = new TH1F("True_Gamma_Phi","True_Gamma_Phi",360,0,2*TMath::Pi());
188 fTrueList->Add(hTrueGammaPhi);
189 hTrueGammaEta = new TH1F("True_Gamma_Eta","True_Gamma_Eta",400,-1.2,1.2);
190 fTrueList->Add(hTrueGammaEta);
191 hTrueGammaMass = new TH1F("True_Gamma_Mass","True_Gamma_Mass",1000,0,0.3);
192 fTrueList->Add(hTrueGammaMass);
193 hTrueGammaChi2perNDF = new TH1F("True_Gamma_Chi2perNDF","True_Gamma_Chi2perNDF",500,0,100);
194 fTrueList->Add(hTrueGammaChi2perNDF);
195 hTrueGammaPsiPair = new TH1F("True_Gamma_PsiPair","True_Gamma_PsiPair",500,0,2);
196 fTrueList->Add(hTrueGammaPsiPair);
197 hTrueGammaQt = new TH1F("True_Gamma_Qt","True_Gamma_Qt",400,0,0.1);
198 fTrueList->Add(hTrueGammaQt);
199 hTrueGammaCosinePointingAngle = new TH1F("True_Gamma_CosinePointingAngle","True_Gamma_CosinePointingAngle",900,0.7,1.);
200 fTrueList->Add(hTrueGammaCosinePointingAngle);
201 hTrueGammaXY = new TH2F("True_Gamma_ConversionPoint_XY","True_Gamma_ConversionPoint_XY",960,-120,120,960,-120,120);
202 fTrueList->Add(hTrueGammaXY);
203 hTrueGammaZR= new TH2F("TrueGamma_ConversionPoint_ZR","TrueGamma_ConversionPoint_ZR",1200,-150,150,480,0,120);
204 fTrueList->Add(hTrueGammaZR);
206 hTrueElecPt = new TH2F("True_Electron_Positron_Pt","True_Electron_Positron_Pt",250,0,25,250,0,25);
207 fTrueList->Add(hTrueElecPt);
208 hTrueElecEta = new TH2F("True_Electron_Positron_Eta","True_Electron_Positron_Eta",400,-1.2,1.2,400,-1.2,1.2);
209 fTrueList->Add(hTrueElecEta);
210 hTrueElecPhi = new TH2F("True_Electron_Positron_Phi","True_Electron_Positron_Phi",360,0,2*TMath::Pi(),360,0,2*TMath::Pi());
211 fTrueList->Add(hTrueElecPhi);
212 hTrueElecNfindableClsTPC = new TH1F("True_Electron_findableClusterTPC","True_Electron_findableClusterTPC",100,0,1);
213 fTrueList->Add(hTrueElecNfindableClsTPC);
214 hTruePosiNfindableClsTPC = new TH1F("True_Positron_findableClusterTPC","True_Positron_findableClusterTPC",100,0,1);
215 fTrueList->Add(hTruePosiNfindableClsTPC);
217 if(fConversionCuts->GetCutHistograms()){
218 fOutputList->Add(fConversionCuts->GetCutHistograms());
220 TH1::SetDefaultSumw2(kFALSE);
224 TString cutnumber = fConversionCuts->GetCutNumber();
225 fStreamQA = new TTreeSRedirector(Form("GammaConvV1_QATree_%s.root",cutnumber.Data()),"recreate");
228 PostData(1, fOutputList);
231 //________________________________________________________________________
232 void AliAnalysisTaskConversionQA::UserExec(Option_t *){
234 fV0Reader=(AliV0ReaderV1*)AliAnalysisManager::GetAnalysisManager()->GetTask("V0ReaderV1");
236 Int_t eventQuality = ((AliConversionCuts*)fV0Reader->GetConversionCuts())->GetEventQuality();
237 if(eventQuality != 0){// Event Not Accepted
240 fInputEvent = InputEvent();
241 fMCEvent = MCEvent();
242 if(fMCEvent) fMCStack = fMCEvent->Stack();
244 Int_t eventNotAccepted =
245 fConversionCuts->IsEventAcceptedByConversionCut(fV0Reader->GetConversionCuts(),fInputEvent,fMCEvent,fIsHeavyIon);
246 if(eventNotAccepted) return; // Check Centrality, PileUp, SDD and V0AND --> Not Accepted => eventQuality = 1
248 fConversionGammas=fV0Reader->GetReconstructedGammas();
251 if(fConversionCuts->GetSignalRejection() != 0){
252 fConversionCuts->GetNotRejectedParticles(fConversionCuts->GetSignalRejection(),
253 fConversionCuts->GetAcceptedHeader(),
260 hVertexZ->Fill(fInputEvent->GetPrimaryVertex()->GetZ());
261 hNContributorsVertex->Fill(fConversionCuts->GetNumberOfContributorsVtx(fInputEvent));
262 hNGoodESDTracks->Fill(fNumberOfESDTracks);
263 hNV0Tracks->Fill(fInputEvent->GetVZEROData()->GetMTotV0A()+fInputEvent->GetVZEROData()->GetMTotV0C());
266 for(Int_t firstGammaIndex=0;firstGammaIndex<fConversionGammas->GetEntriesFast();firstGammaIndex++){
267 AliAODConversionPhoton *gamma=dynamic_cast<AliAODConversionPhoton*>(fConversionGammas->At(firstGammaIndex));
268 if(fMCEvent && fConversionCuts->GetSignalRejection() != 0){
269 if(!fConversionCuts->IsParticleFromBGEvent(gamma->GetMCLabelPositive(), fMCStack))
271 if(!fConversionCuts->IsParticleFromBGEvent(gamma->GetMCLabelNegative(), fMCStack))
274 if(!fConversionCuts->PhotonIsSelected(gamma,fInputEvent)){
278 if(ffillTree) ProcessQATree(gamma);
279 if(ffillHistograms) ProcessQA(gamma);
282 PostData(1, fOutputList);
286 ///________________________________________________________________________
287 void AliAnalysisTaskConversionQA::ProcessQATree(AliAODConversionPhoton *gamma){
289 // Fill Histograms for QA and MC
290 AliESDEvent* event = (AliESDEvent*) InputEvent();
291 AliPIDResponse* pidResonse = ((AliConversionCuts*)fV0Reader->GetConversionCuts())->GetPIDResponse();
293 Float_t gammaPt = gamma->GetPhotonPt();
294 Float_t gammaPhi = gamma->GetPhotonPhi();
295 Float_t gammaTheta = gamma->Theta();
296 Float_t gammaChi2NDF = gamma->GetChi2perNDF();
297 Float_t gammaQt = gamma->GetArmenterosQt();
298 Float_t gammaAlpha = gamma->GetArmenterosAlpha();
299 Float_t gammaPsiPair = gamma->GetPsiPair();
300 Float_t gammaCosPointing = fConversionCuts->GetCosineOfPointingAngle(gamma,event);
301 TVectorF conversionPoint(3);
302 conversionPoint(0) = gamma->GetConversionX();
303 conversionPoint(1) = gamma->GetConversionY();
304 conversionPoint(2) = gamma->GetConversionZ();
305 TVectorF daughterProp(18);
306 AliESDtrack * negTrack = fConversionCuts->GetESDTrack(event, gamma->GetTrackLabelNegative());
307 AliESDtrack * posTrack = fConversionCuts->GetESDTrack(event, gamma->GetTrackLabelPositive());
309 if(!negTrack||!posTrack)return;
311 Bool_t isTruePhoton = kFALSE;
313 if(IsTruePhoton(gamma)) isTruePhoton = kTRUE;
316 daughterProp(0) = posTrack->Pt();
317 daughterProp(7) = negTrack->Pt();
318 daughterProp(1) = posTrack->Theta();
319 daughterProp(8) = negTrack->Theta();
320 Double32_t signalPos[4] = {0,0,0,0};
323 if (posTrack->GetTPCdEdxInfo()) {
324 posTrack->GetTPCdEdxInfo()->GetTPCSignalRegionInfo(signalPos,nclPos,nrowsPos);
325 daughterProp(2) = signalPos[0];
326 daughterProp(14) = signalPos[1];
327 daughterProp(16) = signalPos[2];
329 daughterProp(2) = posTrack->GetTPCsignal();
330 daughterProp(14) = 0;
331 daughterProp(16) = 0;
333 Double32_t signalNeg[4] = {0,0,0,0};
336 if (negTrack->GetTPCdEdxInfo()) {
337 negTrack->GetTPCdEdxInfo()->GetTPCSignalRegionInfo(signalNeg,nclNeg,nrowsNeg);
338 daughterProp(9) = signalNeg[0];
339 daughterProp(15) = signalNeg[1];
340 daughterProp(17) = signalNeg[2];
342 daughterProp(9) = negTrack->GetTPCsignal();
343 daughterProp(15) = 0;
344 daughterProp(17) = 0;
346 daughterProp(3) = pidResonse->NumberOfSigmasTPC(posTrack,AliPID::kElectron);
347 daughterProp(10) = pidResonse->NumberOfSigmasTPC(negTrack,AliPID::kElectron);
348 if((posTrack->GetStatus() & AliESDtrack::kTOFpid) && !(posTrack->GetStatus() & AliESDtrack::kTOFmismatch)){
349 Double_t t0pos = pidResonse->GetTOFResponse().GetStartTime(posTrack->P());
350 Double_t timesPos[5];
351 posTrack->GetIntegratedTimes(timesPos);
352 Double_t TOFsignalPos = posTrack->GetTOFsignal();
353 Double_t dTpos = TOFsignalPos - t0pos - timesPos[0];
354 daughterProp(4) = dTpos;
355 daughterProp(5) = pidResonse->NumberOfSigmasTOF(posTrack, AliPID::kElectron);
357 daughterProp(4) = 20000;
358 daughterProp(5) = -20;
360 if((negTrack->GetStatus() & AliESDtrack::kTOFpid) && !(negTrack->GetStatus() & AliESDtrack::kTOFmismatch)){
361 Double_t t0neg = pidResonse->GetTOFResponse().GetStartTime(negTrack->P());
362 Double_t timesNeg[5];
363 negTrack->GetIntegratedTimes(timesNeg);
364 Double_t TOFsignalNeg = negTrack->GetTOFsignal();
365 Double_t dTneg = TOFsignalNeg - t0neg - timesNeg[0];
366 daughterProp(11) = dTneg;
367 daughterProp(12) = pidResonse->NumberOfSigmasTOF(negTrack, AliPID::kElectron);
369 daughterProp(11) = 20000;
370 daughterProp(12) = -20;
373 daughterProp(6) = (Float_t)posTrack->GetTPCClusterInfo(2,0,fConversionCuts->GetFirstTPCRow(gamma->GetConversionRadius()));
374 //posTrack->GetNcls(1)/(Float_t)posTrack->GetTPCNclsF();
375 daughterProp(13) = (Float_t)negTrack->GetTPCClusterInfo(2,0,fConversionCuts->GetFirstTPCRow(gamma->GetConversionRadius()));
376 //negTrack->GetNcls(1)/(Float_t)negTrack->GetTPCNclsF();
380 (*fStreamQA)<<"PhotonQA"
382 << "phi=" << gammaPhi
383 << "theta=" << gammaTheta
384 << "chi2ndf=" << gammaChi2NDF
386 << "alpha=" << gammaAlpha
387 << "psipair=" << gammaPsiPair
388 << "cosPoint=" << gammaCosPointing
389 << "TruePhoton=" << isTruePhoton
390 << "conversionPoint=" << &conversionPoint
391 << "daugtherProp.=" << &daughterProp
395 (*fStreamQA)<<"PhotonQA"
397 << "phi=" << gammaPhi
398 << "theta=" << gammaTheta
399 << "chi2ndf=" << gammaChi2NDF
401 << "alpha=" << gammaAlpha
402 << "psipair=" << gammaPsiPair
403 << "cosPoint=" << gammaCosPointing
404 << "conversionPoint=" << &conversionPoint
405 << "daugtherProp.=" << &daughterProp
411 //_____________________________________________________________________________________________________
412 void AliAnalysisTaskConversionQA::ProcessQA(AliAODConversionPhoton *gamma){
415 // Fill Histograms for QA and MC
416 hGammaPt->Fill(gamma->GetPhotonPt());
417 hGammaPhi->Fill(gamma->GetPhotonPhi());
418 hGammaEta->Fill(gamma->Eta());
419 hGammaChi2perNDF->Fill(gamma->GetChi2perNDF());
420 hGammaPsiPair->Fill(gamma->GetPsiPair());
421 hGammaQt->Fill(gamma->GetArmenterosQt());
422 hGammaCosinePointingAngle->Fill(fConversionCuts->GetCosineOfPointingAngle(gamma,fInputEvent));
423 hGammaXY->Fill(gamma->GetConversionX(),gamma->GetConversionY());
424 hGammaZR->Fill(gamma->GetConversionZ(),gamma->GetConversionRadius());
426 AliESDtrack * negTrack = fConversionCuts->GetESDTrack((AliESDEvent*)fInputEvent, gamma->GetTrackLabelNegative());
427 AliESDtrack * posTrack = fConversionCuts->GetESDTrack((AliESDEvent*)fInputEvent, gamma->GetTrackLabelPositive());
428 if(!negTrack||!posTrack)return;
430 Double_t negtrackPhi = negTrack->Phi();
431 Double_t postrackPhi = posTrack->Phi();
432 hITSClusterPhi->Fill(negtrackPhi,6);
433 hITSClusterPhi->Fill(postrackPhi,6);
434 for(Int_t itsLayer = 0; itsLayer<6;itsLayer++){
435 if(TESTBIT(negTrack->GetITSClusterMap(),itsLayer)){
436 hITSClusterPhi->Fill(negtrackPhi,itsLayer);
438 if(TESTBIT(posTrack->GetITSClusterMap(),itsLayer)){
439 hITSClusterPhi->Fill(postrackPhi,itsLayer);
443 hElecPt->Fill(negTrack->Pt(),posTrack->Pt());
444 hElecEta->Fill(negTrack->Eta(),posTrack->Eta());
445 hElecPhi->Fill(negTrack->Phi(),posTrack->Phi());
447 hElecNfindableClsTPC->Fill((Float_t)posTrack->GetTPCClusterInfo(2,0,fConversionCuts->GetFirstTPCRow(gamma->GetConversionRadius())));
448 hPosiNfindableClsTPC->Fill((Float_t)negTrack->GetTPCClusterInfo(2,0,fConversionCuts->GetFirstTPCRow(gamma->GetConversionRadius())));
450 // hElecNfindableClsTPC->Fill((Float_t)negTrack->GetNcls(1)/(Float_t)negTrack->GetTPCNclsF());
451 // hPosiNfindableClsTPC->Fill((Float_t)posTrack->GetNcls(1)/(Float_t)posTrack->GetTPCNclsF());
453 if(fMCEvent) ProcessTrueQA(gamma,negTrack,posTrack);
456 //________________________________________________________________________
457 void AliAnalysisTaskConversionQA::ProcessTrueQA(AliAODConversionPhoton *TruePhotonCandidate, AliESDtrack *elec, AliESDtrack *posi)
460 if(!IsTruePhoton(TruePhotonCandidate)) return;
462 TParticle *negDaughter = TruePhotonCandidate->GetNegativeMCDaughter(fMCStack);
463 TParticle *mcPhoton = TruePhotonCandidate->GetMCParticle(fMCStack);
466 hTrueResoulutionR->Fill(TruePhotonCandidate->GetConversionRadius(),
467 TruePhotonCandidate->GetConversionRadius()-negDaughter->R());
468 hTrueResoulutionZ->Fill(TruePhotonCandidate->GetConversionZ(),
469 TruePhotonCandidate->GetConversionZ()-negDaughter->Vz());
470 hTrueResoulutionPhi->Fill(TruePhotonCandidate->Phi(),
471 TruePhotonCandidate->Phi()-mcPhoton->Phi());
472 hTrueGammaPt->Fill(TruePhotonCandidate->Pt());
473 hTrueGammaPhi->Fill(TruePhotonCandidate->Phi());
474 hTrueGammaEta->Fill(TruePhotonCandidate->Eta());
475 hTrueGammaMass->Fill(TruePhotonCandidate->GetMass());
476 hTrueGammaChi2perNDF->Fill(TruePhotonCandidate->GetChi2perNDF());
477 hTrueGammaPsiPair->Fill(TruePhotonCandidate->GetPsiPair());
478 hTrueGammaQt->Fill(TruePhotonCandidate->GetArmenterosQt());
479 hTrueGammaCosinePointingAngle->Fill(fConversionCuts->GetCosineOfPointingAngle(TruePhotonCandidate,fInputEvent));
480 hTrueGammaXY->Fill(TruePhotonCandidate->GetConversionX(),TruePhotonCandidate->GetConversionY());
481 hTrueGammaZR->Fill(TruePhotonCandidate->GetConversionZ(),TruePhotonCandidate->GetConversionRadius());
483 hTrueElecPt->Fill(elec->Pt(),posi->Pt());
484 hTrueElecEta->Fill(elec->Eta(),posi->Eta());
485 hTrueElecPhi->Fill(elec->Phi(),posi->Phi());
486 hTrueElecNfindableClsTPC
487 ->Fill((Float_t)elec->GetTPCClusterInfo(2,0,fConversionCuts->GetFirstTPCRow(TruePhotonCandidate->GetConversionRadius())));
488 hTruePosiNfindableClsTPC
489 ->Fill((Float_t)posi->GetTPCClusterInfo(2,0,fConversionCuts->GetFirstTPCRow(TruePhotonCandidate->GetConversionRadius())));
490 // hTrueElecNfindableClsTPC->Fill((Float_t)elec->GetNcls(1)/(Float_t)elec->GetTPCNclsF());
491 // hTruePosiNfindableClsTPC->Fill((Float_t)posi->GetNcls(1)/(Float_t)posi->GetTPCNclsF());
494 //________________________________________________________________________
495 void AliAnalysisTaskConversionQA::CountESDTracks(){
497 AliESDtrackCuts *EsdTrackCuts = new AliESDtrackCuts("AliESDtrackCuts");
498 // Using standard function for setting Cuts
499 Bool_t selectPrimaries=kTRUE;
500 EsdTrackCuts = AliESDtrackCuts::GetStandardITSTPCTrackCuts2010(selectPrimaries);
501 EsdTrackCuts->SetMaxDCAToVertexZ(2);
502 EsdTrackCuts->SetEtaRange(-0.8, 0.8);
503 EsdTrackCuts->SetPtRange(0.15);
505 fNumberOfESDTracks = 0;
506 for(Int_t iTracks = 0; iTracks < fInputEvent->GetNumberOfTracks(); iTracks++){
507 AliESDtrack* curTrack = (AliESDtrack*)fInputEvent->GetTrack(iTracks);
508 if(!curTrack) continue;
509 // if(fMCEvent && fConversionCuts->GetSignalRejection() != 0){
510 // if(!fConversionCuts->IsParticleFromBGEvent(abs(curTrack->GetLabel()), fMCStack)) continue;
512 if(EsdTrackCuts->AcceptTrack(curTrack) ) fNumberOfESDTracks++;
519 //________________________________________________________________________
520 Bool_t AliAnalysisTaskConversionQA::IsTruePhoton(AliAODConversionPhoton *TruePhotonCandidate)
522 TParticle *posDaughter = TruePhotonCandidate->GetPositiveMCDaughter(fMCStack);
523 TParticle *negDaughter = TruePhotonCandidate->GetNegativeMCDaughter(fMCStack);
525 if(posDaughter == NULL || negDaughter == NULL) return kFALSE; // One particle does not exist
526 Int_t pdgCode[2] = {abs(posDaughter->GetPdgCode()),abs(negDaughter->GetPdgCode())};
527 if(posDaughter->GetMother(0) != negDaughter->GetMother(0)) return kFALSE;
528 else if(posDaughter->GetMother(0) == -1) return kFALSE;
530 if(pdgCode[0]!=11 || pdgCode[1]!=11) return kFALSE; //One Particle is not electron
531 if(posDaughter->GetPdgCode()==negDaughter->GetPdgCode()) return kFALSE; // Same Charge
532 if(posDaughter->GetUniqueID() != 5 || negDaughter->GetUniqueID() !=5) return kFALSE;// check if the daughters come from a conversion
534 TParticle *Photon = TruePhotonCandidate->GetMCParticle(fMCStack);
535 if(Photon->GetPdgCode() != 22) return kFALSE; // Mother is no Photon
540 //________________________________________________________________________
541 void AliAnalysisTaskConversionQA::Terminate(Option_t *)