]> git.uio.no Git - u/mrichter/AliRoot.git/blame - PWG3/vertexingHF/AliAnalysisTaskSESignificance.cxx
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
Changes related to PID mode and parameters configuration (Rossella, Francesco)
[u/mrichter/AliRoot.git] / PWG3 / vertexingHF / AliAnalysisTaskSESignificance.cxx
CommitLineData
cbedddce 1/**************************************************************************
2 * Copyright(c) 1998-2008, ALICE Experiment at CERN, All rights reserved. *
3 * *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
6 * *
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 **************************************************************************/
15
16/////////////////////////////////////////////////////////////
17//
18// AliAnalysisTaskSESignificane to calculate effects on
19// significance of D mesons cut
20// Authors: G. Ortona, ortona@to.infn.it
21// F. Prino, prino@to.infn.it
22// Renu Bala, bala@to.infn.it
23// Chiara Bianchin, cbianchi@pd.infn.it
24/////////////////////////////////////////////////////////////
25
26#include <Riostream.h>
27#include <TClonesArray.h>
28#include <TCanvas.h>
29#include <TList.h>
30#include <TFile.h>
31#include <TString.h>
32#include <TH1F.h>
33#include <TDatabasePDG.h>
34
83e7d8d8 35#include <AliLog.h>
cbedddce 36#include "AliAnalysisManager.h"
37#include "AliAODHandler.h"
38#include "AliAODEvent.h"
39#include "AliAODVertex.h"
40#include "AliAODTrack.h"
41#include "AliAODMCHeader.h"
42#include "AliAODMCParticle.h"
43#include "AliAODRecoDecayHF3Prong.h"
44#include "AliAODRecoDecayHF.h"
45#include "AliAODRecoDecayHF2Prong.h"
46#include "AliAODRecoDecayHF4Prong.h"
47#include "AliAODRecoCascadeHF.h"
48
49#include "AliAnalysisTaskSE.h"
50#include "AliRDHFCutsDplustoKpipi.h"
83e7d8d8 51#include "AliRDHFCutsD0toKpipipi.h"
52#include "AliRDHFCutsDstoKKpi.h"
53#include "AliRDHFCutsDStartoKpipi.h"
cbedddce 54#include "AliRDHFCutsD0toKpi.h"
83e7d8d8 55#include "AliRDHFCutsLctopKpi.h"
cbedddce 56#include "AliMultiDimVector.h"
57
58#include "AliAnalysisTaskSESignificance.h"
59
60ClassImp(AliAnalysisTaskSESignificance)
61
62
63//________________________________________________________________________
64AliAnalysisTaskSESignificance::AliAnalysisTaskSESignificance():
65 AliAnalysisTaskSE(),
66 fOutput(0),
67 fCutList(0),
68 fHistNEvents(0),
69 fUpmasslimit(1.965),
70 fLowmasslimit(1.765),
71 fRDCuts(0),
72 fNPtBins(0),
73 fReadMC(kFALSE),
74 fDecChannel(0),
75 fSelectionlevel(0)
76{
77 // Default constructor
78}
79
80//________________________________________________________________________
81AliAnalysisTaskSESignificance::AliAnalysisTaskSESignificance(const char *name, TList* listMDV,AliRDHFCuts *rdCuts,Int_t decaychannel,Int_t selectionlevel):
82 AliAnalysisTaskSE(name),
83 fOutput(0),
84 fCutList(listMDV),
85 fHistNEvents(0),
86 fUpmasslimit(0),
87 fLowmasslimit(0),
88 fRDCuts(rdCuts),
89 fNPtBins(0),
90 fReadMC(kFALSE),
91 fDecChannel(decaychannel),
92 fSelectionlevel(selectionlevel)
93{
94
95 Int_t pdg=421;
96 switch(fDecChannel){
97 case 0:
98 pdg=411;
99 break;
100 case 1:
101 pdg=421;
102 break;
103 case 2:
104 pdg=413;
105 break;
106 case 3:
107 pdg=431;
108 break;
109 case 4:
110 pdg=421;
111 break;
112 case 5:
113 pdg=4122;
114 break;
115 }
116
117 SetMassLimits(0.1,pdg); //check range
118 fNPtBins=fRDCuts->GetNPtBins();
83e7d8d8 119
cbedddce 120 if(fDebug>1)fRDCuts->PrintAll();
121 // Output slot #1 writes into a TList container
122 DefineOutput(1,TList::Class()); //My private output
123 DefineOutput(2,TList::Class());
3aa6ce53 124 CheckConsistency();
cbedddce 125}
126
127 //________________________________________________________________________
128AliAnalysisTaskSESignificance::~AliAnalysisTaskSESignificance()
129{
130 // Destructor
131 if (fOutput) {
132 delete fOutput;
133 fOutput = 0;
134 }
135 if (fCutList) {
136 delete fCutList;
137 fCutList = 0;
138 }
139 if(fHistNEvents){
140 delete fHistNEvents;
141 fHistNEvents=0;
142 }
143
144/*
145 if(fRDCuts) {
146 delete fRDCuts;
147 fRDCuts= 0;
148 }
149*/
150
3aa6ce53 151}
152//_________________________________________________________________
153Bool_t AliAnalysisTaskSESignificance::CheckConsistency(){
154
155 Bool_t result = kTRUE;
156
157 const Int_t nvars=fRDCuts->GetNVars();//ForOpt();
158 //Float_t *vars = new Float_t[nvars];
159 Bool_t *varsforopt = fRDCuts->GetVarsForOpt();
160 Bool_t *uppervars = fRDCuts->GetIsUpperCut();
161 TString *names = fRDCuts->GetVarNames();
162
163 for(Int_t i=0;i<fNPtBins;i++){
164 TString mdvname=Form("multiDimVectorPtBin%d",i);
165 Int_t ic=0;
166
167 for(Int_t ivar=0;ivar<nvars;ivar++){
168 if(varsforopt[ivar]){
169 Float_t min = ((AliMultiDimVector*)fCutList->FindObject(mdvname.Data()))->GetMinLimit(ic);
170 Float_t max = ((AliMultiDimVector*)fCutList->FindObject(mdvname.Data()))->GetMaxLimit(ic);
171 if(min==max){
172 printf("AliAnalysisTaskSESignificance::CheckConsistency: ERROR! \n tight and loose cut for optimization variable number %d are the same in ptbin %d\n",ic,i);
173 result = kFALSE;
174 }
175 Bool_t lowermdv = ((AliMultiDimVector*)fCutList->FindObject(mdvname.Data()))->GetGreaterThan(ic);
176 if(uppervars[ivar]&&lowermdv){
177 printf("AliAnalysisTaskSESignificance::CheckConsistency: WARNING! %s is declared as uppercut, but as been given tighter cut larger then loose cut in ptbin %d \n ---Task will use swapped Tight/Loose cuts \n ",names[ivar].Data(),i);
178 ((AliMultiDimVector*)fCutList->FindObject(mdvname.Data()))->SwapLimits(ic);
179 result = kTRUE;
180 }
181 if(!uppervars[ivar]&&!lowermdv){
182 printf("AliAnalysisTaskSESignificance::CheckConsistency: WARNING! %s is declared as lower cut, but as been given tighter cut smaller then loose cut in ptbin %d \n ---Task will use swapped Tight/Loose cuts \n",names[ivar].Data(),i);
183 ((AliMultiDimVector*)fCutList->FindObject(mdvname.Data()))->SwapLimits(ic);
184 result = kTRUE;
185 }
186 ic++;
187 }
188 }
189 }
190 return result;
191}
cbedddce 192//_________________________________________________________________
193void AliAnalysisTaskSESignificance::SetMassLimits(Float_t range, Int_t pdg){
194 Float_t mass=0;
195 Int_t abspdg=TMath::Abs(pdg);
5ecc6e9b 196 mass=TDatabasePDG::Instance()->GetParticle(abspdg)->Mass();
cbedddce 197 fUpmasslimit = mass+range;
198 fLowmasslimit = mass-range;
199}
200//_________________________________________________________________
201void AliAnalysisTaskSESignificance::SetMassLimits(Float_t lowlimit, Float_t uplimit){
202 if(uplimit>lowlimit)
203 {
204 fUpmasslimit = lowlimit;
205 fLowmasslimit = uplimit;
206 }
207}
208
209
210
211//________________________________________________________________________
212void AliAnalysisTaskSESignificance::LocalInit()
213{
214 // Initialization
215
216 if(fDebug > 1) printf("AnalysisTaskSESignificance::Init() \n");
217
218 TList *mdvList = new TList();
219 mdvList->SetOwner();
220 mdvList = fCutList;
221
222 PostData(2,mdvList);
223
224 return;
225}
226//________________________________________________________________________
227void AliAnalysisTaskSESignificance::UserCreateOutputObjects()
228{
229 // Create the output container
230
231 if(fDebug > 1) printf("AnalysisTaskSESignificance::UserCreateOutputObjects() \n");
232
233 // Several histograms are more conveniently managed in a TList
234 fOutput = new TList();
235 fOutput->SetOwner();
236 fOutput->SetName("OutputHistos");
237
238 //same number of steps in each multiDimVectorPtBin%d !
239 Int_t nHist=((AliMultiDimVector*)fCutList->FindObject("multiDimVectorPtBin0"))->GetNTotCells();
83e7d8d8 240 cout<<"ncells = "<<nHist<<" n ptbins = "<<fNPtBins<<endl;
cbedddce 241 nHist=nHist*fNPtBins;
83e7d8d8 242 cout<<"Total = "<<nHist<<endl;
cbedddce 243 for(Int_t i=0;i<nHist;i++){
244
245 TString hisname;
246 TString signame;
247 TString bkgname;
248 TString rflname;
249 TString title;
250
251 hisname.Form("hMass_%d",i);
252 signame.Form("hSig_%d",i);
253 bkgname.Form("hBkg_%d",i);
254 rflname.Form("hRfl_%d",i);
255
256 title.Form("Invariant mass;M[GeV/c^{2}];Entries");
257
258 fMassHist[i]=new TH1F(hisname.Data(),title.Data(),100,fLowmasslimit,fUpmasslimit);
259 fMassHist[i]->Sumw2();
260 fOutput->Add(fMassHist[i]);
261
262 if(fReadMC){
263 fSigHist[i]=new TH1F(signame.Data(),title.Data(),100,fLowmasslimit,fUpmasslimit);
264 fSigHist[i]->Sumw2();
265 fOutput->Add(fSigHist[i]);
266
267 fBkgHist[i]=new TH1F(bkgname.Data(),title.Data(),100,fLowmasslimit,fUpmasslimit);
268 fBkgHist[i]->Sumw2();
269 fOutput->Add(fBkgHist[i]);
270
271 fRflHist[i]=new TH1F(rflname.Data(),title.Data(),100,fLowmasslimit,fUpmasslimit);
272 fRflHist[i]->Sumw2();
273 fOutput->Add(fRflHist[i]);
274
275 }
276 }
277
5ecc6e9b 278 fHistNEvents=new TH1F("fHistNEvents","Number of AODs scanned",4,0,4.);
cbedddce 279 fHistNEvents->GetXaxis()->SetBinLabel(1,"nEventsAnal");
83e7d8d8 280 fHistNEvents->GetXaxis()->SetBinLabel(2,"nEvNotSelected");
5ecc6e9b 281 fHistNEvents->GetXaxis()->SetBinLabel(3,"nCandidatesSelected");
282 fHistNEvents->GetXaxis()->SetBinLabel(4,"nTotEntries Mass hists");
cbedddce 283 fHistNEvents->GetXaxis()->SetNdivisions(1,kFALSE);
284 fOutput->Add(fHistNEvents);
285
286
287 return;
288}
289
290//________________________________________________________________________
291void AliAnalysisTaskSESignificance::UserExec(Option_t */*option*/)
292{
293 // Execute analysis for current event:
294 // heavy flavor candidates association to MC truth
295
296 AliAODEvent *aod = dynamic_cast<AliAODEvent*> (InputEvent());
83e7d8d8 297 if(fDebug>2) printf("Analysing decay %d\n",fDecChannel);
cbedddce 298 // Post the data already here
299 PostData(1,fOutput);
300 TClonesArray *arrayProng =0;
301
302 if(!aod && AODEvent() && IsStandardAOD()) {
303 // In case there is an AOD handler writing a standard AOD, use the AOD
304 // event in memory rather than the input (ESD) event.
305 aod = dynamic_cast<AliAODEvent*> (AODEvent());
306 // in this case the braches in the deltaAOD (AliAOD.VertexingHF.root)
307 // have to taken from the AOD event hold by the AliAODExtension
308 AliAODHandler* aodHandler = (AliAODHandler*)
309 ((AliAnalysisManager::GetAnalysisManager())->GetOutputEventHandler());
310 if(aodHandler->GetExtensions()) {
311
312 AliAODExtension *ext = (AliAODExtension*)aodHandler->GetExtensions()->FindObject("AliAOD.VertexingHF.root");
313 AliAODEvent *aodFromExt = ext->GetAOD();
314
315
316 switch(fDecChannel){
317 case 0:
318 arrayProng=(TClonesArray*)aodFromExt->GetList()->FindObject("Charm3Prong");
319 break;
320 case 1:
321 arrayProng=(TClonesArray*)aodFromExt->GetList()->FindObject("D0toKpi");
322 break;
323 case 2:
324 arrayProng=(TClonesArray*)aodFromExt->GetList()->FindObject("Dstar");
325 break;
326 case 3:
327 arrayProng=(TClonesArray*)aodFromExt->GetList()->FindObject("Charm3Prong");
328 break;
329 case 4:
330 arrayProng=(TClonesArray*)aodFromExt->GetList()->FindObject("Charm4Prong");
331 break;
332 case 5:
333 arrayProng=(TClonesArray*)aodFromExt->GetList()->FindObject("Charm3Prong");
334 break;
335 }
336 }
337 } else {
338 switch(fDecChannel){
339 case 0:
340 arrayProng=(TClonesArray*)aod->GetList()->FindObject("Charm3Prong");
341 break;
342 case 1:
343 arrayProng=(TClonesArray*)aod->GetList()->FindObject("D0toKpi");
344 break;
345 case 2:
346 arrayProng=(TClonesArray*)aod->GetList()->FindObject("Dstar");
347 break;
348 case 3:
349 arrayProng=(TClonesArray*)aod->GetList()->FindObject("Charm3Prong");
350 break;
351 case 4:
352 arrayProng=(TClonesArray*)aod->GetList()->FindObject("Charm4Prong");
353 break;
354 case 5:
355 arrayProng=(TClonesArray*)aod->GetList()->FindObject("Charm3Prong");
356 break;
357 }
358 }
359 if(!arrayProng) {
360 printf("AliAnalysisTaskSESignificance::UserExec: branch not found!\n");
361 return;
362 }
363
364 TClonesArray *arrayMC=0;
365 AliAODMCHeader *mcHeader=0;
366
367 // load MC particles
368 if(fReadMC){
369
370 arrayMC = (TClonesArray*)aod->GetList()->FindObject(AliAODMCParticle::StdBranchName());
371 if(!arrayMC) {
372 printf("AliAnalysisTaskSESignificance::UserExec: MC particles branch not found!\n");
373 // return;
374 }
375
376 // load MC header
377 mcHeader = (AliAODMCHeader*)aod->GetList()->FindObject(AliAODMCHeader::StdBranchName());
378 if(!mcHeader) {
379 printf("AliAnalysisTaskSESignificance::UserExec: MC header branch not found!\n");
380 return;
381 }
382 }
383
384 fHistNEvents->Fill(0); // count event
385
386 AliAODRecoDecayHF *d=0;
387 Int_t nprongs=1;
388 Int_t *pdgdaughters=0x0;
389 Int_t prongpdg=211;
390
391
392 switch(fDecChannel){
393 case 0:
394 //D+
395 pdgdaughters =new Int_t[3];
396 pdgdaughters[0]=211;//pi
397 pdgdaughters[1]=321;//K
398 pdgdaughters[2]=211;//pi
399 nprongs=3;
400 prongpdg=411;
401 break;
402 case 1:
403 //D0
404 pdgdaughters =new Int_t[2];
405 pdgdaughters[0]=211;//pi
406 pdgdaughters[1]=321;//K
407 nprongs=2;
408 prongpdg=421;
409 break;
410 case 2:
411 //D*
412 pdgdaughters =new Int_t[3];
413 pdgdaughters[1]=211;//pi
414 pdgdaughters[0]=321;//K
415 pdgdaughters[2]=211;//pi (soft?)
416 nprongs=3;
417 prongpdg=413;
418 break;
419 case 3:
420 //Ds
421 pdgdaughters =new Int_t[3];
422 pdgdaughters[0]=211;//pi
423 pdgdaughters[1]=321;//K
424 pdgdaughters[2]=321;//K
425 nprongs=3;
426 prongpdg=431;
427 break;
428 case 4:
429 //D0
430 pdgdaughters =new Int_t[4];
431 pdgdaughters[0]=321;
432 pdgdaughters[1]=211;
433 pdgdaughters[2]=211;
434 pdgdaughters[3]=211;
435 nprongs=4;
436 prongpdg=421;
437 break;
438 case 5:
439 //Lambda_c
440 pdgdaughters =new Int_t[3];
441 //check order
442 pdgdaughters[0]=2212;//p
443 pdgdaughters[1]=321;//K
444 pdgdaughters[2]=211;//pi
445 nprongs=3;
446 prongpdg=4122;
447 break;
448 }
449
450 Int_t nHistpermv=((AliMultiDimVector*)fCutList->FindObject("multiDimVectorPtBin0"))->GetNTotCells();
451 Int_t nProng = arrayProng->GetEntriesFast();
452 if(fDebug>1) printf("Number of D2H: %d\n",nProng);
453
5ecc6e9b 454 if(!fRDCuts->IsEventSelected(aod)){
455 fHistNEvents->Fill(1);
456 return;
457 }
458
cbedddce 459 for (Int_t iProng = 0; iProng < nProng; iProng++) {
460 d=(AliAODRecoDecayHF*)arrayProng->UncheckedAt(iProng);
461 Double_t invMass=d->InvMass(nprongs,(UInt_t*)pdgdaughters);
462 Double_t invMassC=0;
463
464 Int_t isSelected=fRDCuts->IsSelected(d,fSelectionlevel);
465
466 if(isSelected) {
5ecc6e9b 467 fHistNEvents->Fill(2); // count selected with loosest cuts
83e7d8d8 468 if(fDebug>1) printf("+++++++Is Selected\n");
cbedddce 469
470 const Int_t nvars=fRDCuts->GetNVarsForOpt();
471 Float_t *vars = new Float_t[nvars];
472 Int_t nVals=0;
473
474 fRDCuts->GetCutVarsForOpt(d,vars,nvars,pdgdaughters);
475 Int_t ptbin=fRDCuts->PtBin(d->Pt());
476 if(ptbin==-1) continue;
477 TString mdvname=Form("multiDimVectorPtBin%d",ptbin);
478 ULong64_t *addresses = ((AliMultiDimVector*)fCutList->FindObject(mdvname.Data()))->GetGlobalAddressesAboveCuts(vars,(Float_t)d->Pt(),nVals);
479 if(fDebug>1)printf("nvals = %d\n",nVals);
cbedddce 480 for(Int_t ivals=0;ivals<nVals;ivals++){
481 if(addresses[ivals]>=((AliMultiDimVector*)fCutList->FindObject(mdvname.Data()))->GetNTotCells()){
482 if (fDebug>1) printf("Overflow!!\n");
483 return;
484 }
485 fMassHist[ptbin*nHistpermv+addresses[ivals]]->Fill(invMass);
5ecc6e9b 486 fHistNEvents->Fill(3);
cbedddce 487 if(fReadMC){
488 Int_t lab=-1;
489 lab = d->MatchToMC(prongpdg,arrayMC,nprongs,pdgdaughters);
490 if(lab>=0){ //signal
491 AliAODMCParticle *dMC = (AliAODMCParticle*)arrayMC->At(lab);
492 Int_t pdgMC = dMC->GetPdgCode();
493
494 if(pdgMC==+prongpdg) fSigHist[ptbin*nHistpermv+addresses[ivals]]->Fill(invMass);
495 else fRflHist[ptbin*nHistpermv+addresses[ivals]]->Fill(invMass);
496 }
497 else{ //background
498 fBkgHist[ptbin*nHistpermv+addresses[ivals]]->Fill(invMass);
499 }
500
501 }
502 }
503
504
83e7d8d8 505 if(fDecChannel==AliAnalysisTaskSESignificance::kD0toKpi || fDecChannel==AliAnalysisTaskSESignificance::kLambdactopKpi){ //repeat the cycle to checks cuts passed as D0bar
506 if(fDecChannel==AliAnalysisTaskSESignificance::kD0toKpi){
507 pdgdaughters[0]=321;
508 pdgdaughters[1]=211;
509 }
510 if(fDecChannel==AliAnalysisTaskSESignificance::kLambdactopKpi){
511 pdgdaughters[0]=211;
512 pdgdaughters[1]=321;
513 pdgdaughters[2]=2212;
514 }
515
cbedddce 516 invMass=d->InvMass(nprongs,(UInt_t*)pdgdaughters);
517 fRDCuts->GetCutVarsForOpt(d,vars,nvars,pdgdaughters);
518 addresses = ((AliMultiDimVector*)fCutList->FindObject(mdvname.Data()))->GetGlobalAddressesAboveCuts(vars,(Float_t)d->Pt(),nVals);
519 for(Int_t ivals=0;ivals<nVals;ivals++){
520 fMassHist[ptbin*nHistpermv+addresses[ivals]]->Fill(invMassC);
521 if(fReadMC){
522 Int_t lab=-1;
523 lab = d->MatchToMC(-prongpdg,arrayMC,nprongs,pdgdaughters);
524 if(lab>=0){ //signal
525 AliAODMCParticle *dMC = (AliAODMCParticle*)arrayMC->At(lab);
526 Int_t pdgMC = dMC->GetPdgCode();
527 if(pdgMC==-prongpdg) fSigHist[ptbin*nHistpermv+addresses[ivals]]->Fill(invMassC);
528 else fRflHist[ptbin*nHistpermv+addresses[ivals]]->Fill(invMassC);
529 }
530 else{ //background
531 fBkgHist[ptbin*nHistpermv+addresses[ivals]]->Fill(invMassC);
532 }
533
534 }
535 }
536 }
ab3e5f67 537//inserisci qui
538 if(fDecChannel==AliAnalysisTaskSESignificance::kLambdactopKpi){ //repeat the cycle to checks cuts passed as Lambdacbar
539 pdgdaughters[0]=211;
540 pdgdaughters[1]=321;
541 pdgdaughters[2]=2212;
542 invMass=d->InvMass(nprongs,(UInt_t*)pdgdaughters);
543 fRDCuts->GetCutVarsForOpt(d,vars,nvars,pdgdaughters);
544 addresses = ((AliMultiDimVector*)fCutList->FindObject(mdvname.Data()))->GetGlobalAddressesAboveCuts(vars,(Float_t)d->Pt(),nVals);
545 for(Int_t ivals=0;ivals<nVals;ivals++){
546 fMassHist[ptbin*nHistpermv+addresses[ivals]]->Fill(invMassC);
547 if(fReadMC){
548 Int_t lab=-1;
549 lab = d->MatchToMC(-prongpdg,arrayMC,nprongs,pdgdaughters);
550 //cambia il match
551 if(lab>=0){ //signal
552 AliAODMCParticle *dMC = (AliAODMCParticle*)arrayMC->At(lab);
553 Int_t pdgMC = dMC->GetPdgCode();
554 if(pdgMC==-prongpdg) fSigHist[ptbin*nHistpermv+addresses[ivals]]->Fill(invMassC);
555 else fRflHist[ptbin*nHistpermv+addresses[ivals]]->Fill(invMassC);
556 }
557 else{ //background
558 fBkgHist[ptbin*nHistpermv+addresses[ivals]]->Fill(invMassC);
559 }
560 }
561 }
562 }
cbedddce 563
564 }// end if selected
565
566
567 }
568
569 PostData(1,fOutput);
570 return;
571}
572
573
574//________________________________________________________________________
575void AliAnalysisTaskSESignificance::Terminate(Option_t */*option*/)
576{
577 // Terminate analysis
578 //
579 if(fDebug > 1) printf("AnalysisTaskSESignificance: Terminate() \n");
580
581
582 fOutput = dynamic_cast<TList*> (GetOutputData(1));
583 if (!fOutput) {
584 printf("ERROR: fOutput not available\n");
585 return;
586 }
587
588 fCutList = dynamic_cast<TList*> (GetOutputData(2));
589 if (!fCutList) {
590 printf("ERROR: fCutList not available\n");
591 return;
592 }
593
594 AliMultiDimVector* mdvtmp=(AliMultiDimVector*)fCutList->FindObject("multiDimVectorPtBin0");
595 if (!mdvtmp){
596 cout<<"multidimvec not found in TList"<<endl;
597 fCutList->ls();
598 return;
599 }
600 Int_t nHist=mdvtmp->GetNTotCells();
601 TCanvas *c1=new TCanvas("c1","Invariant mass distribution - loose cuts",500,500);
602 Bool_t drawn=kFALSE;
603 for(Int_t i=0;i<nHist;i++){
604
605 TString hisname;
606 TString signame;
607 TString bkgname;
608 TString rflname;
609
610 hisname.Form("hMass_%d",i);
611 signame.Form("hSig_%d",i);
612 bkgname.Form("hBkg_%d",i);
613 rflname.Form("hRfl_%d",i);
614
615 fMassHist[i]=dynamic_cast<TH1F*>(fOutput->FindObject(hisname.Data()));
616
617 if (!drawn && fMassHist[i]->GetEntries() > 0){
618 c1->cd();
619 fMassHist[i]->Draw();
620 drawn=kTRUE;
621 }
622
623 if(fReadMC){
624 fSigHist[i]=dynamic_cast<TH1F*>(fOutput->FindObject(signame.Data()));
625 fBkgHist[i]=dynamic_cast<TH1F*>(fOutput->FindObject(bkgname.Data()));
626 fRflHist[i]=dynamic_cast<TH1F*>(fOutput->FindObject(bkgname.Data()));
627 }
628
629 }
630
631
632
633
634 return;
635}
ab3e5f67 636//-------------------------------------------
cbedddce 637