1 //------------------------------------------------------------------------------
2 // Implementation of AliPerformanceMatch class. It checks the track matching
3 // quality (residuals, pulls) bewteen tracking detectors TPC-ITS and TPC-TRD
4 // with TPC as the reference detector. In addition, the ITS and TRD track
5 // reconstruction efficiency with respect to TPC is calculated.
7 // The matchinig quality parameters are stored in the THnSparse histograms.
8 // The analysis of these histograms to extract reference information is done in
9 // the Analyse() function (post-processing).
11 // The histograms with reference information can be exported to the ROOT folder.
13 // Author: J.Otwinowski 17/10/2009
14 // Changes by M.Knichel 22/10/2010
15 //------------------------------------------------------------------------------
19 // after running comparison task, read the file, and get component
20 gROOT->LoadMacro("$ALICE_ROOT/PWGPP/Macros/LoadMyLibs.C");
23 TFile f("Output.root");
24 AliPerformanceMatch * compObj = (AliPerformanceMatch*)coutput->FindObject("AliPerformanceMatchTPCTRD");
26 // analyse comparison data
29 // the output histograms/graphs will be stored in the folder "folderMatch"
30 compObj->GetAnalysisFolder()->ls("*");
32 // user can save whole comparison object (or only folder with anlysed histograms)
33 // in the seperate output file (e.g.)
34 TFile fout("Analysed_Match.root","recreate");
35 compObj->Write(); // compObj->GetAnalysisFolder()->Write();
46 #include "AliPerformanceMatch.h"
47 #include "AliESDEvent.h"
48 #include "AliESDVertex.h"
49 #include "AliESDtrack.h"
50 #include "AliESDfriendTrack.h"
51 #include "AliESDfriend.h"
53 #include "AliMCEvent.h"
54 #include "AliMCParticle.h"
55 #include "AliHeader.h"
56 #include "AliGenEventHeader.h"
58 #include "AliMCInfoCuts.h"
59 #include "AliRecInfoCuts.h"
60 #include "AliTracker.h"
61 #include "AliTRDtrackV1.h"
62 #include "AliTreeDraw.h"
66 ClassImp(AliPerformanceMatch)
68 Bool_t AliPerformanceMatch::fgMergeTHnSparse = kFALSE;
69 Bool_t AliPerformanceMatch::fgUseMergeTHnSparse = kFALSE;
71 //_____________________________________________________________________________
72 AliPerformanceMatch::AliPerformanceMatch():
73 AliPerformanceObject("AliPerformanceMatch"),
91 //_____________________________________________________________________________
92 AliPerformanceMatch::AliPerformanceMatch(Char_t* name="AliPerformanceMatch", Char_t* title="AliPerformanceMatch",Int_t analysisMode=0,Bool_t hptGenerator=kFALSE):
93 AliPerformanceObject(name,title),
109 SetAnalysisMode(analysisMode);
110 SetHptGenerator(hptGenerator);
115 //_____________________________________________________________________________
116 AliPerformanceMatch::~AliPerformanceMatch()
120 if(fResolHisto) delete fResolHisto; fResolHisto=0;
121 if(fPullHisto) delete fPullHisto; fPullHisto=0;
122 if(fTrackingEffHisto) delete fTrackingEffHisto; fTrackingEffHisto = 0x0;
123 if(fTPCConstrain) delete fTPCConstrain; fTPCConstrain = 0x0;
125 if(fAnalysisFolder) delete fAnalysisFolder; fAnalysisFolder=0;
126 if(fFolderObj) delete fFolderObj; fFolderObj=0;
129 //_____________________________________________________________________________
130 void AliPerformanceMatch::Init(){
133 // Make performance histogrms
138 Double_t ptMin = 1.e-2, ptMax = 20.;
140 Double_t *binsPt = 0;
142 if (IsHptGenerator()) {
145 binsPt = CreateLogAxis(nPtBins,ptMin,ptMax);
148 // Double_t yMin = -0.02, yMax = 0.02;
149 // Double_t zMin = -12.0, zMax = 12.0;
150 // if(GetAnalysisMode() == 0 || GetAnalysisMode() == 1) {
151 // yMin = -100.; yMax = 100.;
152 // zMin = -100.; zMax = 100.;
156 //init ITS TPC Mactching
158 if(GetAnalysisMode()==1){
159 // res_y:res_z:res_phi,res_lambda:res_pt:y:z:eta:phi:pt:isRec
160 Int_t binsResolHisto[9]={100,100,100,100,100,90,30,nPtBins,2};
161 Double_t minResolHisto[9]={-1.,-1.,-0.03,-0.03,-0.2, 0., -1.5, ptMin,0};
162 Double_t maxResolHisto[9]={ 1., 1., 0.03, 0.03, 0.2, 2.*TMath::Pi(), 1.5, ptMax,2};
164 fResolHisto = new THnSparseF("fResolHisto","res_y:res_z:res_phi:res_lambda:res_pt:phi:eta:pt:isRec",9,binsResolHisto,minResolHisto,maxResolHisto);
165 fResolHisto->SetBinEdges(7,binsPt);
167 fResolHisto->GetAxis(0)->SetTitle("y-y_{ref} (cm)");
168 fResolHisto->GetAxis(1)->SetTitle("z-z_{ref} (cm)");
169 fResolHisto->GetAxis(2)->SetTitle("#phi-#phi_{ref} (rad)");
170 fResolHisto->GetAxis(3)->SetTitle("#lambda-#lambda_{ref} (rad)");
171 fResolHisto->GetAxis(4)->SetTitle("(p_{T}/p_{Tref}-1)");
172 fResolHisto->GetAxis(5)->SetTitle("#phi_{ref} (rad)");
173 fResolHisto->GetAxis(6)->SetTitle("#eta_{ref}");
174 fResolHisto->GetAxis(7)->SetTitle("p_{Tref} (GeV/c)");
175 fResolHisto->GetAxis(8)->SetTitle("isReconstructed");
176 fResolHisto->Sumw2();
178 //pull_y:pull_z:pull_snp:pull_tgl:pull_1pt:y:z:snp:tgl:1pt:isRec
179 Int_t binsPullHisto[9]={100,100,100,100,100,90,30,nPtBins,2};
180 Double_t minPullHisto[9]={-5.,-5.,-5.,-5.,-5., 0.,-1.5, ptMin,0};
181 Double_t maxPullHisto[9]={ 5., 5., 5., 5., 5., 2.*TMath::Pi(), 1.5, ptMax,2};
182 fPullHisto = new THnSparseF("fPullHisto","pull_y:pull_z:pull_snp:pull_tgl:pull_1pt:phi:eta:1pt:isRec",9,binsPullHisto,minPullHisto,maxPullHisto);
184 fPullHisto->GetAxis(0)->SetTitle("(y-y_{ref})/#sigma");
185 fPullHisto->GetAxis(1)->SetTitle("(z-z_{ref})/#sigma");
186 fPullHisto->GetAxis(2)->SetTitle("(sin#phi-sin#phi_{ref})/#sigma");
187 fPullHisto->GetAxis(3)->SetTitle("(tan#lambda-tan#lambda_{ref})/#sigma");
188 fPullHisto->GetAxis(4)->SetTitle("(p_{Tref}/p_{T}-1)/#sigma");
189 fPullHisto->GetAxis(5)->SetTitle("#phi_{ref} (rad)");
190 fPullHisto->GetAxis(6)->SetTitle("#eta_{ref}");
191 fPullHisto->GetAxis(7)->SetTitle("1/p_{Tref} (GeV/c)^{-1}");
192 fPullHisto->GetAxis(8)->SetTitle("isReconstructed");
199 if(GetAnalysisMode()==0){
200 // -> has match:y:z:snp:tgl:phi:pt:ITSclusters
201 Int_t binsTrackingEffHisto[5] = { 2, 90, nPtBins, 30, 7 };
202 Double_t minTrackingEffHisto[5] = {-0.5, 0., ptMin, -1.5, -0.5 };
203 Double_t maxTrackingEffHisto[5] = { 1.5, 2*TMath::Pi(), ptMax, 1.5, 6.5 };
205 fTrackingEffHisto = new THnSparseF("fTrackingEffHisto","has match:phi:pt:eta:ITSclusters",5,binsTrackingEffHisto,minTrackingEffHisto,maxTrackingEffHisto);
206 fTrackingEffHisto->SetBinEdges(2,binsPt);
207 fTrackingEffHisto->GetAxis(0)->SetTitle("IsMatching");
208 fTrackingEffHisto->GetAxis(1)->SetTitle("phi (rad)");
209 fTrackingEffHisto->GetAxis(2)->SetTitle("p_{T}");
210 fTrackingEffHisto->GetAxis(3)->SetTitle("eta");
211 fTrackingEffHisto->GetAxis(4)->SetTitle("number of ITS clusters");
212 fTrackingEffHisto->Sumw2();
216 //TPC constraining to vertex
218 if(GetAnalysisMode()==2){
219 //initilization of fTPCConstrain
220 Int_t binsTPCConstrain[4] = {100, 90, nPtBins, 30};
221 Double_t minTPCConstrain[4] = {-5, 0, ptMin, -1.5};
222 Double_t maxTPCConstrain[4] = {5, 2*TMath::Pi(), ptMax, 1.5};
224 fTPCConstrain = new THnSparseF("fTPCConstrain","pull_phi:phi:pt:eta",4, binsTPCConstrain,minTPCConstrain,maxTPCConstrain);
225 fTPCConstrain->SetBinEdges(2,binsPt);
226 fTPCConstrain->GetAxis(0)->SetTitle("(#phi-#phi_{ref})/#sigma");
227 fTPCConstrain->GetAxis(1)->SetTitle("phi (rad)");
228 fTPCConstrain->GetAxis(2)->SetTitle("p_{T}");
229 fTPCConstrain->GetAxis(3)->SetTitle("eta");
230 fTPCConstrain->Sumw2();
235 AliDebug(AliLog::kError, "ERROR: Cannot find AliMCInfoCuts object");
237 AliDebug(AliLog::kError, "ERROR: Cannot find AliRecInfoCuts object");
240 fAnalysisFolder = CreateFolder("folderMatch","Analysis Matching Folder");
242 // save merge status in object
243 fMergeTHnSparseObj = fgMergeTHnSparse;
249 //_____________________________________________________________________________
250 void AliPerformanceMatch::ProcessITSTPC(Int_t iTrack, AliESDEvent *const esdEvent, AliStack* /*const stack*/, AliESDtrack *const esdTrack)
253 // addition to standard analysis - check if ITS stand-alone tracks have a match in the TPC
254 // Origin: A. Kalwait
255 // Modified: J. Otwinowski
256 if(!esdEvent) return;
257 if(!esdTrack) return;
258 // if(!esdFriendTrack) return;
260 // ITS stand alone tracks with SPD layers
261 if(!(esdTrack->GetStatus() & AliESDtrack::kITSpureSA)) return;
262 if(!(esdTrack->GetStatus() & AliESDtrack::kITSrefit)) return;
263 if(esdTrack->GetNcls(0)<4) return;
264 if(!esdTrack->HasPointOnITSLayer(0) && !esdTrack->HasPointOnITSLayer(1)) return;
266 const AliESDVertex* vtxESD = esdEvent->GetPrimaryVertexTracks();
267 // const AliESDVertex *vtxESD = esdEvent->GetPrimaryVertexTPC();
268 AliESDtrack* tpcTrack2 = NULL;
269 Bool_t hasMatch = kFALSE;
270 for (Int_t jTrack = 0; jTrack < esdEvent->GetNumberOfTracks(); jTrack++) {
271 // loop for all those tracks and check if the corresponding TPC track is found
272 if (jTrack==iTrack) continue;
273 AliESDtrack *trackTPC = esdEvent->GetTrack(jTrack);
274 if (!trackTPC) continue;
275 if (!trackTPC->GetTPCInnerParam()) continue;
276 if(!(trackTPC->GetStatus() & AliESDtrack::kTPCrefit)) continue;
278 // TPC nClust/track after first tracking pass
279 // if(trackTPC->GetTPCNclsIter1()<fCutsRC->GetMinNClustersTPC()) continue;
280 tpcTrack2 = AliESDtrackCuts::GetTPCOnlyTrack(esdEvent, jTrack);
281 if(!tpcTrack2) continue;
282 if(!tpcTrack2->RelateToVertex(vtxESD,esdEvent->GetMagneticField(),100.)) { delete tpcTrack2; continue; }
284 if(!fCutsRC->AcceptTrack(tpcTrack2)) { delete tpcTrack2; continue; }
286 if (TMath::Abs(esdTrack->GetY() - tpcTrack2->GetY()) > 3) { delete tpcTrack2; continue; }
287 if (TMath::Abs(esdTrack->GetSnp() - tpcTrack2->GetSnp()) > 0.2) { delete tpcTrack2; continue; }
288 if (TMath::Abs(esdTrack->GetTgl() - tpcTrack2->GetTgl()) > 0.2) { delete tpcTrack2; continue; }
294 FillHistograms(tpcTrack2,esdTrack,hasMatch);
301 //_____________________________________________________________________________
302 void AliPerformanceMatch::ProcessTPCITS(AliStack* /*const stack*/, AliESDtrack *const esdTrack)
305 // Match TPC and ITS min-bias tracks
306 // at radius between detectors
308 if(!esdTrack) return;
309 // if(!esdFriendTrack) return;
311 Bool_t isTPC = kFALSE;
312 Bool_t isMatch = kFALSE;
315 if(esdTrack->Charge()==0) return;
316 if(!esdTrack->GetTPCInnerParam()) return;
317 if(!(esdTrack->GetStatus()&AliESDtrack::kTPCrefit)) return;
319 if(!fCutsRC->AcceptTrack(esdTrack)) return;
323 if( (esdTrack->GetStatus()&AliESDtrack::kITSrefit))
327 Double_t vecTrackingEff[5] = { isMatch,esdTrack->Phi(), esdTrack->Pt(),esdTrack->Eta(),esdTrack->GetITSclusters(0) };
328 fTrackingEffHisto->Fill(vecTrackingEff);
332 //_____________________________________________________________________________
333 /*void AliPerformanceMatch::ProcessTPCTRD(AliStack* , AliESDtrack *const esdTrack, AliESDfriendTrack *const esdFriendTrack)
338 void AliPerformanceMatch::ProcessTPCConstrain(AliStack* /*const stack*/, AliESDEvent *const esdEvent, AliESDtrack *const esdTrack){
340 // Contrain TPC inner track to the vertex
341 // then compare to the global tracks
343 if(!esdTrack) return;
344 if(esdTrack->Charge()==0) return;
346 if(!esdTrack->GetTPCInnerParam()) return;
347 if(!(esdTrack->GetStatus()&AliESDtrack::kITSrefit)) return;
348 if(!(esdTrack->GetStatus()&AliESDtrack::kTPCrefit)) return;
349 if(!fCutsRC->AcceptTrack(esdTrack)) return;
351 Double_t x[3]; esdTrack->GetXYZ(x);
352 Double_t b[3]; AliTracker::GetBxByBz(x,b);
353 Bool_t isOK = kFALSE;
354 const AliESDVertex* vtxESD = esdEvent->GetPrimaryVertexTracks();
356 AliExternalTrackParam * TPCinner = (AliExternalTrackParam *)esdTrack->GetTPCInnerParam();
357 if(!TPCinner) return;
359 // isOK = TPCinner->Rotate(esdTrack->GetAlpha());
360 //isOK = TPCinner->PropagateTo(esdTrack->GetX(),esdEvent->GetMagneticField());
362 AliExternalTrackParam * TPCinnerC = new AliExternalTrackParam(*TPCinner);
363 isOK = TPCinnerC->ConstrainToVertex(vtxESD, b);
365 // transform to the track reference frame
366 isOK = TPCinnerC->Rotate(esdTrack->GetAlpha());
367 isOK = TPCinnerC->PropagateTo(esdTrack->GetX(),esdEvent->GetMagneticField());
370 Double_t sigmaPhi=0,deltaPhi=0,pullPhi=0;
371 deltaPhi = TPCinnerC->GetSnp() - esdTrack->GetSnp();
372 sigmaPhi = TMath::Sqrt(esdTrack->GetSigmaSnp2()+TPCinnerC->GetSigmaSnp2());
374 pullPhi = deltaPhi/sigmaPhi;
376 Double_t vTPCConstrain[4] = {pullPhi,esdTrack->Phi(),esdTrack->Pt(),esdTrack->Eta()};
377 fTPCConstrain->Fill(vTPCConstrain);
384 //_____________________________________________________________________________
385 void AliPerformanceMatch::FillHistograms(AliESDtrack *const refParam, AliESDtrack *const param, Bool_t isRec)
388 // fill performance histograms
389 // (TPC always as reference)
392 if(!refParam) return;
397 // Deltas (dy,dz,dphi,dtheta,dpt)
399 Float_t delta[5] = {0};
401 delta[0] = param->GetY()-refParam->GetY();
402 delta[1] = param->GetZ()-refParam->GetZ();
403 delta[2] = TMath::ATan2(param->Py(),param->Px())-TMath::ATan2(refParam->Py(),refParam->Px());
404 delta[3] = TMath::ATan2(param->Pz(),param->Pt())-TMath::ATan2(refParam->Pz(),refParam->Pt());
405 if(refParam->Pt()) delta[4] = (param->Pt()-refParam->Pt())/refParam->Pt();
408 // Pulls (y,z,snp,tanl,1/pt)
410 Float_t sigma[5] = {0};
411 Float_t pull[5] = {0};
413 sigma[0] = TMath::Sqrt(param->GetSigmaY2()+refParam->GetSigmaY2());
414 sigma[1] = TMath::Sqrt(param->GetSigmaZ2()+refParam->GetSigmaZ2());
415 sigma[2] = TMath::Sqrt(param->GetSigmaSnp2()+refParam->GetSigmaSnp2());
416 sigma[3] = TMath::Sqrt(param->GetSigmaTgl2()+refParam->GetSigmaTgl2());
417 sigma[4] = TMath::Sqrt(param->GetSigma1Pt2()+refParam->GetSigma1Pt2());
418 if(sigma[0]) pull[0] = delta[0] / sigma[0];
419 if(sigma[1]) pull[1] = delta[1] / sigma[1];
420 if(sigma[2]) pull[2] = (param->GetSnp()-refParam->GetSnp()) / sigma[2];
421 if(sigma[3]) pull[3] = (param->GetTgl()-refParam->GetTgl()) / sigma[3];
422 if(sigma[4]) pull[4] = (param->OneOverPt()-refParam->OneOverPt()) / sigma[4];
426 Double_t vResolHisto[9] = {delta[0],delta[1],delta[2],delta[3],delta[4],refParam->Phi(),refParam->Eta(),refParam->Pt(),isRec};
428 fResolHisto->Fill(vResolHisto);
430 Double_t vPullHisto[9] = {pull[0],pull[1],pull[2],pull[3],pull[4],refParam->Phi(),refParam->Eta(),refParam->OneOverPt(),isRec};
432 fPullHisto->Fill(vPullHisto);
435 //_____________________________________________________________________________
436 void AliPerformanceMatch::Exec(AliMCEvent* const mcEvent, AliESDEvent *const esdEvent, AliESDfriend */*const esdFriend*/, const Bool_t bUseMC, const Bool_t /*bUseESDfriend*/)
438 // Process comparison information
442 Error("Exec","esdEvent not available");
445 AliHeader* header = 0;
446 AliGenEventHeader* genHeader = 0;
453 Error("Exec","mcEvent not available");
456 // get MC event header
457 header = mcEvent->Header();
459 Error("Exec","Header not available");
463 stack = mcEvent->Stack();
465 Error("Exec","Stack not available");
469 genHeader = header->GenEventHeader();
471 Error("Exec","Could not retrieve genHeader from Header");
474 genHeader->PrimaryVertex(vtxMC);
478 /* if(bUseESDfriend) {
480 Error("Exec","esdFriend not available");
486 if(!bUseMC && GetTriggerClass()) {
487 Bool_t isEventTriggered = esdEvent->IsTriggerClassFired(GetTriggerClass());
488 if(!isEventTriggered) return;
491 // get TPC event vertex
492 const AliESDVertex *vtxESD = esdEvent->GetPrimaryVertexTPC();
493 if(vtxESD && (vtxESD->GetStatus()<=0)) return;
496 for (Int_t iTrack = 0; iTrack < esdEvent->GetNumberOfTracks(); iTrack++)
498 AliESDtrack *track = esdEvent->GetTrack(iTrack);
501 /*AliESDfriendTrack *friendTrack=0;
503 friendTrack=esdFriend->GetTrack(iTrack);
504 if(!friendTrack) continue;
507 if(GetAnalysisMode() == 0){
508 if(!IsUseTOFBunchCrossing()){
509 ProcessTPCITS(stack,track);
512 if( track->GetTOFBunchCrossing(esdEvent->GetMagneticField())==0) {
513 ProcessTPCITS(stack,track);
516 /* else if(GetAnalysisMode() == 2) ProcessTPCTRD(stack,track,friendTrack);*/
517 else if(GetAnalysisMode() == 1) {ProcessITSTPC(iTrack,esdEvent,stack,track);}
518 else if(GetAnalysisMode() == 2){
519 if(!IsUseTOFBunchCrossing()){
520 ProcessTPCConstrain(stack,esdEvent,track);
523 if( track->GetTOFBunchCrossing(esdEvent->GetMagneticField())==0) {
524 ProcessTPCConstrain(stack,esdEvent,track);
528 printf("ERROR: AnalysisMode %d \n",fAnalysisMode);
535 //_____________________________________________________________________________
536 TH1F* AliPerformanceMatch::MakeResol(TH2F * his, Int_t integ, Bool_t type, Int_t cut){
537 // Create resolution histograms
540 if (!gPad) new TCanvas;
541 hisr = AliTreeDraw::CreateResHistoII(his,&hism,integ,kTRUE,cut);
542 if (type) return hism;
547 //_____________________________________________________________________________
548 void AliPerformanceMatch::Analyse() {
549 // Analyse comparison information and store output histograms
550 // in the folder "folderMatch"
554 TH1::AddDirectory(kFALSE);
559 TObjArray *aFolderObj = new TObjArray;
561 // write results in the folder
562 // TCanvas * c = new TCanvas("Phi resol Tan","Phi resol Tan");
568 if(GetAnalysisMode()==1) {
570 fResolHisto->GetAxis(8)->SetRangeUser(1.0,2.0); // only reconstructed
571 fPullHisto->GetAxis(8)->SetRangeUser(1.0,2.0); // only reconstructed
572 for(Int_t i=0; i<5; i++)
574 for(Int_t j=5; j<8; j++)
576 //if(j!=8) fResolHisto->GetAxis(8)->SetRangeUser(-0.9,0.89); // eta window
577 if(j!=6) fResolHisto->GetAxis(6)->SetRangeUser(0.0,1.49); // eta window
578 else fResolHisto->GetAxis(6)->SetRangeUser(-1.5,1.49);
579 fResolHisto->GetAxis(7)->SetRangeUser(0.01,10.); // pt threshold
582 AddProjection(aFolderObj, "match", fResolHisto, i, j, &selString);
585 if(j!=6) fPullHisto->GetAxis(6)->SetRangeUser(0.0,1.49); // eta window
586 else fPullHisto->GetAxis(6)->SetRangeUser(-1.5,1.49); // eta window
587 fPullHisto->GetAxis(7)->SetRangeUser(0.01,10.); // pt threshold
589 AddProjection(aFolderObj, "match", fPullHisto, i, j, &selString);
595 // TPC efficiency wrt ITS
597 if(GetAnalysisMode()==0) {
598 selString = "trackingeff";
599 AddProjection(aFolderObj, "match", fTrackingEffHisto, 0, &selString);
601 for(Int_t i=1; i<5; i++)
603 // all ITS standalone tracks
604 fTrackingEffHisto->GetAxis(0)->SetRange(1,fTrackingEffHisto->GetAxis(0)->GetNbins());
605 selString = "trackingeff_all";
606 AddProjection(aFolderObj, "match", fTrackingEffHisto, i, 3,&selString);
608 // TPC tracks which has matching with TPC
609 fTrackingEffHisto->GetAxis(0)->SetRange(2,2);
610 selString = "trackingeff_tpc";
611 AddProjection(aFolderObj, "match", fTrackingEffHisto, i, 3,&selString);
616 //TPC constrained to vertex
618 if(GetAnalysisMode()==2) {
620 // for(Int_t i=1; i<4; i++)
621 AddProjection(aFolderObj, "constrain", fTPCConstrain, 0, 1, 2, &selString);
622 AddProjection(aFolderObj, "constrain", fTPCConstrain, 0, 1, 3, &selString);
623 AddProjection(aFolderObj, "constrain", fTPCConstrain, 0, 2, 3, &selString);
626 printf("exportToFolder\n");
627 fAnalysisFolder = ExportToFolder(aFolderObj);
629 // delete only TObjArray
630 if(fFolderObj) delete fFolderObj;
631 fFolderObj = aFolderObj;
637 //_____________________________________________________________________________
638 TFolder* AliPerformanceMatch::ExportToFolder(TObjArray * array)
640 // recreate folder avery time and export objects to new one
642 AliPerformanceMatch * comp=this;
643 TFolder *folder = comp->GetAnalysisFolder();
646 TFolder *newFolder = 0;
648 Int_t size = array->GetSize();
651 // get name and title from old folder
652 name = folder->GetName();
653 title = folder->GetTitle();
659 newFolder = CreateFolder(name.Data(),title.Data());
660 newFolder->SetOwner();
662 // add objects to folder
664 newFolder->Add(array->At(i));
672 //_____________________________________________________________________________
673 TFolder* AliPerformanceMatch::CreateFolder(TString name,TString title) {
674 // create folder for analysed histograms
677 folder = new TFolder(name.Data(),title.Data());
682 //_____________________________________________________________________________
683 Long64_t AliPerformanceMatch::Merge(TCollection* const list)
685 // Merge list of objects (needed by PROOF)
693 Bool_t merge = ((fgUseMergeTHnSparse && fgMergeTHnSparse) || (!fgUseMergeTHnSparse && fMergeTHnSparseObj));
695 TIterator* iter = list->MakeIterator();
697 TObjArray* objArrayList = 0;
698 objArrayList = new TObjArray();
700 // collection of generated histograms
702 while((obj = iter->Next()) != 0)
704 AliPerformanceMatch* entry = dynamic_cast<AliPerformanceMatch*>(obj);
705 if (entry == 0) continue;
707 if ((fResolHisto) && (entry->fResolHisto)) { fResolHisto->Add(entry->fResolHisto); }
708 if ((fPullHisto) && (entry->fPullHisto)) { fPullHisto->Add(entry->fPullHisto); }
709 if ((fTrackingEffHisto) && (entry->fTrackingEffHisto)) { fTrackingEffHisto->Add(entry->fTrackingEffHisto); }
711 if ((fTPCConstrain) && (entry->fTPCConstrain)) { fTPCConstrain->Add(entry->fTPCConstrain); }
713 // the analysisfolder is only merged if present
714 if (entry->fFolderObj) { objArrayList->Add(entry->fFolderObj); }
718 if (fFolderObj) { fFolderObj->Merge(objArrayList); }
719 // to signal that track histos were not merged: reset
720 if (!merge) { fResolHisto->Reset(); fPullHisto->Reset(); fTrackingEffHisto->Reset(); fTPCConstrain->Reset(); }
722 if (objArrayList) delete objArrayList; objArrayList=0;