-/**************************************************************************\r
- * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *\r
- * *\r
- * Author: The ALICE Off-line Project. *\r
- * Contributors are mentioned in the code where appropriate. *\r
- * *\r
- * Permission to use, copy, modify and distribute this software and its *\r
- * documentation strictly for non-commercial purposes is hereby granted *\r
- * without fee, provided that the above copyright notice appears in all *\r
- * copies and that both the copyright notice and this permission notice *\r
- * appear in the supporting documentation. The authors make no claims *\r
- * about the suitability of this software for any purpose. It is *\r
- * provided "as is" without express or implied warranty. *\r
- **************************************************************************/\r
-//------------------------------------------------------------------------------\r
-// AlidNdPtCutAnalysisPbPb class. \r
-//\r
-// a. functionality:\r
-// - fills generic cut histograms\r
-// - generates cuts (selection criteria)\r
-//\r
-// b. data members:\r
-// - generic cut histograms\r
-// - control histograms\r
-//\r
-// Author: J.Otwinowski 04/11/2008 \r
-//------------------------------------------------------------------------------\r
-#include "TH1.h"\r
-#include "TH2.h"\r
-#include "TH3.h"\r
-#include "TMatrixD.h"\r
-\r
-#include <TPDGCode.h>\r
-\r
-#include "AliHeader.h" \r
-#include "AliGenEventHeader.h" \r
-#include "AliInputEventHandler.h" \r
-#include "AliCentrality.h" \r
-#include "AliStack.h" \r
-#include "AliESDEvent.h" \r
-#include "AliMCEvent.h" \r
-#include "AliESDtrackCuts.h" \r
-#include "AliLog.h" \r
-#include "AliTracker.h" \r
-\r
-#include "AlidNdPtEventCuts.h"\r
-#include "AlidNdPtAcceptanceCuts.h"\r
-#include "AlidNdPtBackgroundCuts.h"\r
-#include "AlidNdPtAnalysis.h"\r
-#include "AliPhysicsSelection.h"\r
-\r
-#include "AliPWG0Helper.h"\r
-#include "AlidNdPtHelper.h"\r
-#include "AlidNdPtCutAnalysisPbPb.h"\r
-\r
-using namespace std;\r
-\r
-ClassImp(AlidNdPtCutAnalysisPbPb)\r
-\r
-//_____________________________________________________________________________\r
- AlidNdPtCutAnalysisPbPb::AlidNdPtCutAnalysisPbPb(): AlidNdPt(),\r
- fAnalysisFolder(0),\r
- fEventCount(0),\r
- fRecEventHist(0),\r
- fMCEventHist(0),\r
- fRecMCEventHist(0),\r
- fRecMCTrackHist(0),\r
- fCentralityEstimator(0),\r
- fFolderObj(0)\r
-{\r
- // default constructor\r
- Init();\r
-}\r
-\r
-//_____________________________________________________________________________\r
-AlidNdPtCutAnalysisPbPb::AlidNdPtCutAnalysisPbPb(Char_t* name, Char_t* title): AlidNdPt(name,title),\r
- fAnalysisFolder(0),\r
- fEventCount(0),\r
- fRecEventHist(0),\r
- fMCEventHist(0),\r
- fRecMCEventHist(0),\r
- fRecMCTrackHist(0),\r
- fCentralityEstimator(0),\r
- fFolderObj(0)\r
-{\r
- // constructor\r
- Init();\r
-}\r
-\r
-//_____________________________________________________________________________\r
-AlidNdPtCutAnalysisPbPb::~AlidNdPtCutAnalysisPbPb() {\r
- // \r
- if(fEventCount) delete fEventCount; fEventCount=0;\r
- if(fRecEventHist) delete fRecEventHist; fRecEventHist=0;\r
- if(fMCEventHist) delete fMCEventHist; fMCEventHist=0;\r
- if(fRecMCEventHist) delete fRecMCEventHist; fRecMCEventHist=0;\r
- if(fRecMCTrackHist) delete fRecMCTrackHist; fRecMCTrackHist=0;\r
-\r
- if(fAnalysisFolder) delete fAnalysisFolder; fAnalysisFolder=0;\r
- if(fFolderObj) delete fFolderObj; fFolderObj=0;\r
-}\r
-\r
-//_____________________________________________________________________________\r
-void AlidNdPtCutAnalysisPbPb::Init(){\r
- //\r
- // Init histograms\r
- //\r
- //const Int_t ptNbins = 58; \r
- //const Double_t ptMin = 0.; \r
- //const Double_t ptMax = 20.; \r
- /*\r
- Double_t binsPt[ptNbins+1] = {0.,0.05,0.1,0.15,0.2,0.25,0.3,0.35,0.4,0.45,0.5,0.55,0.6,0.65,0.7,0.75,0.8,0.85,0.9,0.95,1.0,1.1,1.2,1.3,1.4,1.5,1.6,1.7,1.8,1.9,2.0,2.2,2.4,2.6,2.8,3.0,3.2,3.4,3.6,3.8,4.0,4.5,5.0,5.5,6.0,6.5,7.0,7.5,8.0,9.0,10.0,11.0,12.0,13.0,14.0,15.0,16.0,18.0, 20.};\r
- */\r
-\r
- // set pt bins\r
- const Int_t ptNbins = 50;\r
- const Double_t ptMin = 1.e-2, ptMax = 100.;\r
- Double_t *binsPt = CreateLogAxis(ptNbins,ptMin,ptMax);\r
-\r
- // centrality bins\r
- const Int_t centrNbins = 3;\r
- const Double_t centrMin = 0, centrMax = 1;\r
- Double_t binsCentr[centrNbins+1] = {0.0, 20., 50., 100.};\r
-\r
- // \r
- Int_t binsEventCount[2]={2,2};\r
- Double_t minEventCount[2]={0,0}; \r
- Double_t maxEventCount[2]={2,2}; \r
- fEventCount = new THnSparseF("fEventCount","trig vs trig+vertex",2,binsEventCount,minEventCount,maxEventCount);\r
- fEventCount->GetAxis(0)->SetTitle("trig");\r
- fEventCount->GetAxis(1)->SetTitle("trig+vert");\r
- fEventCount->Sumw2();\r
-\r
- //Xv:Yv:Zv:ResZv:Mult\r
- Double_t kFact = 0.1;\r
-\r
- Int_t binsRecEventHist[5]={80,80,100,80,150};\r
- Double_t minRecEventHist[5]={-1.*kFact,-1.*kFact,-35.,0.,0.}; \r
- Double_t maxRecEventHist[5]={1.*kFact,1.*kFact,35.,10.,3000.}; \r
- fRecEventHist = new THnSparseF("fRecEventHist","Xv:Yv:Zv:ResZv:Mult",5,binsRecEventHist,minRecEventHist,maxRecEventHist);\r
- fRecEventHist->GetAxis(0)->SetTitle("Xv (cm)");\r
- fRecEventHist->GetAxis(1)->SetTitle("Yv (cm)");\r
- fRecEventHist->GetAxis(2)->SetTitle("Zv (cm)");\r
- fRecEventHist->GetAxis(3)->SetTitle("ResZv (cm)");\r
- fRecEventHist->GetAxis(4)->SetTitle("Mult");\r
- fRecEventHist->Sumw2();\r
-\r
- //Xv:Yv:Zv\r
- Int_t binsMCEventHist[3]={80,80,100};\r
- Double_t minMCEventHist[3]={-0.1,-0.1,-35.}; \r
- Double_t maxMCEventHist[3]={0.1,0.1,35.}; \r
- fMCEventHist = new THnSparseF("fMCEventHist","mcXv:mcYv:mcZv",3,binsMCEventHist,minMCEventHist,maxMCEventHist);\r
- fMCEventHist->GetAxis(0)->SetTitle("mcXv (cm)");\r
- fMCEventHist->GetAxis(1)->SetTitle("mcYv (cm)");\r
- fMCEventHist->GetAxis(2)->SetTitle("mcZv (cm)");\r
- fMCEventHist->Sumw2();\r
-\r
- //Xv-mcXv:Yv-mcYv:Zv-mcZv:Mult\r
- Int_t binsRecMCEventHist[4]={100,100,100,150};\r
- Double_t minRecMCEventHist[4]={-1.0*kFact,-1.0*kFact,-1.0*kFact,0.}; \r
- Double_t maxRecMCEventHist[4]={1.0*kFact,1.0*kFact,1.0*kFact,3000.}; \r
- fRecMCEventHist = new THnSparseF("fRecMCEventHist","Xv-mcXv:Yv-mcYv:Zv-mcZv:Mult",4,binsRecMCEventHist,minRecMCEventHist,maxRecMCEventHist);\r
- fRecMCEventHist->GetAxis(0)->SetTitle("Xv-mcXv (cm)");\r
- fRecMCEventHist->GetAxis(1)->SetTitle("Yv-mcYv (cm)");\r
- fRecMCEventHist->GetAxis(2)->SetTitle("Zv-mcZv (cm)");\r
- fRecMCEventHist->GetAxis(3)->SetTitle("Mult");\r
- fRecMCEventHist->Sumw2();\r
-\r
- //\r
- // THnSparse track histograms\r
- //\r
- //\r
- //\r
- //nCrossRows:chi2PerClust:chi2PerClustITS:nCrossRows/nFindableClust:fracSharedClust:DCAy:DCAz:eta:phi:pt:hasStrangeMother:isFromConversion:isFromMaterial:isPrim:charge:centr:chi2ToTPCc\r
- Int_t binsRecMCTrackHist[17]= {160, 10, 70, 30, 20, 50, 50, 20, 90, ptNbins, 2, 2, 2, 2, 3, centrNbins,100};\r
- Double_t minRecMCTrackHist[17]={0., 0., 0., 0., 0.,-0.5,-2.5,-1.0, 0., ptMin, 0., 0., 0., 0.,-1., centrMin, 0.};\r
- Double_t maxRecMCTrackHist[17]={160.,10.,70.,1.5, 1., 0.5, 2.5, 1.0, 2.*TMath::Pi(), ptMax, 2., 2., 2., 2., 2., centrMax, 100.};\r
-\r
- fRecMCTrackHist = new THnSparseF("fRecMCTrackHist","nCrossRows:chi2PerClust:chi2PerClustITS:nCrossRows/nFindableClust:fracSharedClust:DCAy:DCAz:eta:phi:pt:hasStrangeMother:isFromConversion:isFromMaterial:isPrim:charge:centr:chi2ToTPCc",17,binsRecMCTrackHist,minRecMCTrackHist,maxRecMCTrackHist);\r
- fRecMCTrackHist->SetBinEdges(9,binsPt);\r
- fRecMCTrackHist->SetBinEdges(15,binsCentr);\r
-\r
- fRecMCTrackHist->GetAxis(0)->SetTitle("nCrossRows");\r
- fRecMCTrackHist->GetAxis(1)->SetTitle("chi2PerClust");\r
- fRecMCTrackHist->GetAxis(2)->SetTitle("chi2PerClustITS");\r
- fRecMCTrackHist->GetAxis(3)->SetTitle("nCrossRows/nFindableClust");\r
- fRecMCTrackHist->GetAxis(4)->SetTitle("fracSharedClust");\r
- fRecMCTrackHist->GetAxis(5)->SetTitle("DCAy (cm)");\r
- fRecMCTrackHist->GetAxis(6)->SetTitle("DCAz (cm)");\r
- fRecMCTrackHist->GetAxis(7)->SetTitle("#eta");\r
- fRecMCTrackHist->GetAxis(8)->SetTitle("#phi (rad)");\r
- fRecMCTrackHist->GetAxis(9)->SetTitle("p_{T} (GeV/c)");\r
- fRecMCTrackHist->GetAxis(10)->SetTitle("hasStrangeMother");\r
- fRecMCTrackHist->GetAxis(11)->SetTitle("isFromConversion");\r
- fRecMCTrackHist->GetAxis(12)->SetTitle("isFromMaterial");\r
- fRecMCTrackHist->GetAxis(13)->SetTitle("isPrim");\r
- fRecMCTrackHist->GetAxis(14)->SetTitle("charge");\r
- fRecMCTrackHist->GetAxis(15)->SetTitle("centrality");\r
- fRecMCTrackHist->GetAxis(16)->SetTitle("chi2ToTPCc");\r
- fRecMCTrackHist->Sumw2();\r
-\r
- // init output folder\r
- fAnalysisFolder = CreateFolder("folderdNdPt","Analysis dNdPt Folder");\r
-\r
-}\r
-\r
-//_____________________________________________________________________________\r
-void AlidNdPtCutAnalysisPbPb::Process(AliESDEvent *const esdEvent, AliMCEvent * const mcEvent)\r
-{\r
- //\r
- // Process real and/or simulated events\r
- //\r
- if(!esdEvent) {\r
- AliDebug(AliLog::kError, "esdEvent not available");\r
- return;\r
- }\r
-\r
- // get selection cuts\r
- AlidNdPtEventCuts *evtCuts = GetEventCuts(); \r
- AlidNdPtAcceptanceCuts *accCuts = GetAcceptanceCuts(); \r
- AliESDtrackCuts *esdTrackCuts = GetTrackCuts(); \r
-\r
- if(!evtCuts || !accCuts || !esdTrackCuts) {\r
- AliDebug(AliLog::kError, "cuts not available");\r
- return;\r
- }\r
-\r
- // trigger selection\r
-\r
- Bool_t isEventTriggered = kTRUE;\r
- AliPhysicsSelection *physicsSelection = NULL;\r
- AliTriggerAnalysis* triggerAnalysis = NULL;\r
-\r
- // \r
- AliInputEventHandler* inputHandler = (AliInputEventHandler*) AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler();\r
- if (!inputHandler)\r
- {\r
- Printf("ERROR: Could not receive input handler");\r
- return;\r
- }\r
-\r
- if(evtCuts->IsTriggerRequired()) \r
- {\r
- // always MB\r
- isEventTriggered = inputHandler->IsEventSelected() & AliVEvent::kMB;\r
-\r
- physicsSelection = static_cast<AliPhysicsSelection*> (inputHandler->GetEventSelection());\r
- if(!physicsSelection) return;\r
- //SetPhysicsTriggerSelection(physicsSelection);\r
-\r
- if (isEventTriggered && (GetTrigger() == AliTriggerAnalysis::kV0AND)) {\r
- // set trigger (V0AND)\r
- triggerAnalysis = physicsSelection->GetTriggerAnalysis();\r
- if(!triggerAnalysis) return;\r
- isEventTriggered = triggerAnalysis->IsOfflineTriggerFired(esdEvent, GetTrigger());\r
- }\r
- }\r
-\r
- // centrality determination\r
- Float_t centralityF = -1;\r
- AliCentrality *esdCentrality = esdEvent->GetCentrality();\r
- centralityF = esdCentrality->GetCentralityPercentile(fCentralityEstimator.Data());\r
-\r
- // use MC information\r
- AliHeader* header = 0;\r
- AliGenEventHeader* genHeader = 0;\r
- AliStack* stack = 0;\r
- TArrayF vtxMC(3);\r
-\r
- Int_t multMCTrueTracks = 0;\r
- if(IsUseMCInfo())\r
- {\r
- //\r
- if(!mcEvent) {\r
- AliDebug(AliLog::kError, "mcEvent not available");\r
- return;\r
- }\r
- // get MC event header\r
- header = mcEvent->Header();\r
- if (!header) {\r
- AliDebug(AliLog::kError, "Header not available");\r
- return;\r
- }\r
- // MC particle stack\r
- stack = mcEvent->Stack();\r
- if (!stack) {\r
- AliDebug(AliLog::kError, "Stack not available");\r
- return;\r
- }\r
-\r
- // get MC vertex\r
- genHeader = header->GenEventHeader();\r
- if (!genHeader) {\r
- AliDebug(AliLog::kError, "Could not retrieve genHeader from Header");\r
- return;\r
- }\r
- genHeader->PrimaryVertex(vtxMC);\r
-\r
- Double_t vMCEventHist[3]={vtxMC[0],vtxMC[1],vtxMC[2]};\r
- fMCEventHist->Fill(vMCEventHist);\r
-\r
- // multipliticy of all MC primary tracks\r
- // in Zv, pt and eta ranges)\r
- multMCTrueTracks = AlidNdPtHelper::GetMCTrueTrackMult(mcEvent,evtCuts,accCuts);\r
-\r
- } // end bUseMC\r
-\r
- // get reconstructed vertex \r
- const AliESDVertex* vtxESD = 0; \r
- if(GetAnalysisMode() == AlidNdPtHelper::kTPC) {\r
- vtxESD = esdEvent->GetPrimaryVertexTPC();\r
- }\r
- else if(GetAnalysisMode() == AlidNdPtHelper::kTPCITS) {\r
- vtxESD = esdEvent->GetPrimaryVertexTracks();\r
- }\r
- else {\r
- return;\r
- }\r
-\r
- if(!vtxESD) return;\r
-\r
- Bool_t isEventOK = evtCuts->AcceptEvent(esdEvent,mcEvent,vtxESD); \r
- //printf("isEventOK %d, isEventTriggered %d \n",isEventOK, isEventTriggered);\r
- //printf("GetAnalysisMode() %d \n",GetAnalysisMode());\r
-\r
- TObjArray *allChargedTracks=0;\r
- Int_t multAll=0;\r
- \r
- //\r
- // event counter\r
- // \r
- //printf("isEventOK %d, isEventTriggered %d \n",isEventOK,isEventTriggered);\r
-\r
- Bool_t isTrigAndVertex = isEventTriggered && isEventOK;\r
- Double_t vEventCount[2] = { isEventTriggered, isTrigAndVertex};\r
- fEventCount->Fill(vEventCount);\r
-\r
- // check event cuts\r
- if(isEventOK && isEventTriggered)\r
- {\r
- // get all charged tracks\r
- allChargedTracks = AlidNdPtHelper::GetAllChargedTracks(esdEvent,GetAnalysisMode());\r
- if(!allChargedTracks) return;\r
-\r
- Int_t entries = allChargedTracks->GetEntries();\r
- for(Int_t i=0; i<entries;++i) \r
- {\r
- AliESDtrack *track = (AliESDtrack*)allChargedTracks->At(i);\r
- if(!track) continue;\r
-\r
- if(!esdTrackCuts->AcceptTrack(track)) continue;\r
- if(!accCuts->AcceptTrack(track)) continue;\r
-\r
- FillHistograms(esdEvent, track, stack, centralityF);\r
- multAll++;\r
-\r
- }\r
-\r
- Double_t vRecEventHist[5] = {vtxESD->GetXv(),vtxESD->GetYv(),vtxESD->GetZv(),vtxESD->GetZRes(),multAll};\r
- fRecEventHist->Fill(vRecEventHist);\r
-\r
- if(IsUseMCInfo()) {\r
- Double_t vRecMCEventHist[5] = {vtxESD->GetXv()-vtxMC[0],vtxESD->GetYv()-vtxMC[1],vtxESD->GetZv()-vtxMC[2],multAll};\r
- fRecMCEventHist->Fill(vRecMCEventHist);\r
- }\r
- }\r
-\r
- if(allChargedTracks) delete allChargedTracks; allChargedTracks = 0;\r
-\r
-}\r
-\r
-//_____________________________________________________________________________\r
-void AlidNdPtCutAnalysisPbPb::FillHistograms(AliESDEvent *const esdEvent, AliESDtrack *const esdTrack, AliStack *const stack, Float_t centralityF) const\r
-{\r
- //\r
- // Fill ESD track and MC histograms \r
- //\r
- if(!esdTrack) return;\r
- if(esdTrack->Charge() == 0.) return;\r
-\r
- Float_t pt = esdTrack->Pt();\r
- Float_t eta = esdTrack->Eta();\r
- Float_t phi = esdTrack->Phi();\r
-\r
- Int_t nClust = 0;\r
- if(GetAnalysisMode() == AlidNdPtHelper::kTPC) { \r
- nClust = esdTrack->GetTPCNclsIter1();\r
- } else {\r
- nClust = esdTrack->GetTPCclusters(0);\r
- }\r
-\r
- Float_t chi2PerCluster = 0.;\r
- if(GetAnalysisMode() == AlidNdPtHelper::kTPC) { \r
- if(nClust>0.) chi2PerCluster = esdTrack->GetTPCchi2Iter1()/Float_t(nClust);\r
- } else {\r
- chi2PerCluster = esdTrack->GetTPCchi2()/Float_t(nClust);\r
- }\r
-\r
- Int_t nClustersITS = esdTrack->GetITSclusters(0);\r
- Float_t chi2PerClusterITS = -1;\r
- if (nClustersITS!=0) {\r
- chi2PerClusterITS = esdTrack->GetITSchi2()/Float_t(nClustersITS);\r
- }\r
-\r
- Float_t clustPerFindClust = 0.;\r
- Int_t nFindableClust = esdTrack->GetTPCNclsF();\r
- if(nFindableClust>0.) clustPerFindClust = Float_t(nClust)/nFindableClust;\r
-\r
- Float_t b[2], bCov[3];\r
- esdTrack->GetImpactParameters(b,bCov);\r
-\r
- // \r
- Float_t nCrossedRowsTPC = esdTrack->GetTPCClusterInfo(2,1);\r
-\r
- Float_t ratioCrossedRowsOverFindableClustersTPC = 1.0;\r
- if (esdTrack->GetTPCNclsF()>0) {\r
- ratioCrossedRowsOverFindableClustersTPC = esdTrack->GetTPCClusterInfo(2,1)/esdTrack->GetTPCNclsF();\r
- }\r
-\r
- //\r
- Int_t nClustersTPCShared = esdTrack->GetTPCnclsS();\r
- Float_t fracClustersTPCShared = -1.;\r
- fracClustersTPCShared = Float_t(nClustersTPCShared)/Float_t(nClust);\r
-\r
- // \r
- // kink idx\r
- Int_t kinkIdx = 0;\r
- //if(esdTrack->GetKinkIndex(0) > 0.) isKink = kTRUE;\r
- if(esdTrack->GetKinkIndex(0) > 0) kinkIdx = 1; // kink daughter\r
- else if(esdTrack->GetKinkIndex(0) < 0) kinkIdx = -1; // kink mother\r
- else kinkIdx = 0; // not kink\r
-\r
- //printf("esdTrack->GetKinkIndex(0) %d \n", esdTrack->GetKinkIndex(0));\r
- //printf("esdTrack->GetKinkIndex(1) %d \n", esdTrack->GetKinkIndex(1));\r
- //printf("esdTrack->GetKinkIndex(2) %d \n", esdTrack->GetKinkIndex(2));\r
- //printf("kinkIdx %d \n", kinkIdx);\r
-\r
- Bool_t isOK = kTRUE;\r
- AliExternalTrackParam * tpcInner = (AliExternalTrackParam *)(esdTrack->GetTPCInnerParam());\r
- if (!tpcInner) return;\r
- tpcInner->Rotate(esdTrack->GetAlpha());\r
- tpcInner->PropagateTo(esdTrack->GetX(),esdEvent->GetMagneticField());\r
- \r
- // tpc constrained\r
- AliExternalTrackParam * tpcInnerC = (AliExternalTrackParam *)(esdTrack->GetTPCInnerParam()->Clone());\r
- if (!tpcInnerC) return;\r
- tpcInnerC->Rotate(esdTrack->GetAlpha());\r
- tpcInnerC->PropagateTo(esdTrack->GetX(),esdEvent->GetMagneticField());\r
- Double_t dz[2],cov[3];\r
- AliESDVertex *vtx= (AliESDVertex *)esdEvent->GetPrimaryVertex();\r
-\r
- if (!tpcInnerC->PropagateToDCA(vtx, esdEvent->GetMagneticField(), 3, dz, cov)) \r
- { \r
- if(tpcInnerC) delete tpcInnerC; \r
- return;\r
- }\r
- Double_t covar[6]; vtx->GetCovMatrix(covar);\r
- Double_t p[2]={tpcInnerC->GetParameter()[0]-dz[0],tpcInnerC->GetParameter()[1]-dz[1]};\r
- Double_t c[3]={covar[2],0.,covar[5]};\r
- Double_t chi2C=tpcInnerC->GetPredictedChi2(p,c);\r
- isOK = tpcInnerC->Update(p,c);\r
- isOK = tpcInnerC->Rotate(esdTrack->GetAlpha());\r
- isOK = tpcInnerC->Propagate(esdTrack->GetAlpha(),esdTrack->GetX(),esdEvent->GetMagneticField());\r
-\r
- if(!isOK || chi2C>kVeryBig) { \r
- if(tpcInnerC) delete tpcInnerC; \r
- return;\r
- }\r
-\r
- //\r
- // calculate chi2 between vi and vj vectors\r
- // with covi and covj covariance matrices\r
- // chi2ij = (vi-vj)^(T)*(covi+covj)^(-1)*(vi-vj)\r
- //\r
- TMatrixD deltaT(5,1);\r
- TMatrixD delta(1,5);\r
- TMatrixD covarM(5,5);\r
-\r
- for (Int_t ipar=0; ipar<5; ipar++) deltaT(ipar,0)=tpcInnerC->GetParameter()[ipar]-esdTrack->GetParameter()[ipar];\r
- for (Int_t ipar=0; ipar<5; ipar++) delta(0,ipar)=tpcInnerC->GetParameter()[ipar]-esdTrack->GetParameter()[ipar];\r
-\r
- for (Int_t ipar=0; ipar<5; ipar++) {\r
- for (Int_t jpar=0; jpar<5; jpar++) {\r
- Int_t index=esdTrack->GetIndex(ipar,jpar);\r
- covarM(ipar,jpar)=esdTrack->GetCovariance()[index]+tpcInnerC->GetCovariance()[index];\r
- }\r
- }\r
- TMatrixD covarMInv = covarM.Invert();\r
- TMatrixD mat2 = covarMInv*deltaT;\r
- TMatrixD chi2 = delta*mat2; \r
-\r
- //\r
- // Fill rec vs MC information\r
- //\r
- Bool_t isPrim = kTRUE;\r
- Bool_t hasStrangeMother = kFALSE;\r
- Bool_t isFromMaterial = kFALSE;\r
- Bool_t isFromConversion = kFALSE;\r
-\r
- if(IsUseMCInfo()) \r
- {\r
- if(!stack) return;\r
- Int_t label = TMath::Abs(esdTrack->GetLabel()); \r
- TParticle* particle = stack->Particle(label);\r
- if(!particle) return;\r
- if(particle->GetPDG() && particle->GetPDG()->Charge()==0.) return;\r
- isPrim = stack->IsPhysicalPrimary(label);\r
-\r
- // check whether has stange mother\r
- //\r
- Int_t motherPdg = -1; \r
- TParticle* mother = 0; \r
- \r
- Int_t motherLabel = particle->GetMother(0); \r
- if(motherLabel>0) mother = stack->Particle(motherLabel); \r
- if(mother) motherPdg = TMath::Abs(mother->GetPdgCode()); // take abs for visualisation only \r
- Int_t mech = particle->GetUniqueID(); // production mechanism \r
-\r
- // K+-, lambda, antilambda, K0s decays\r
- if(!isPrim && mech==4 && \r
- (TMath::Abs(motherPdg)==kKPlus || TMath::Abs(motherPdg)==kLambda0 || motherPdg==kK0Short))\r
- {\r
- hasStrangeMother = kTRUE;\r
- } \r
-\r
- if(!isPrim && mech==5 && motherPdg==kGamma) { \r
- isFromConversion=kTRUE; \r
- }\r
-\r
- if(!isPrim && mech==13) { \r
- isFromMaterial=kTRUE; \r
- }\r
- }\r
- \r
- // fill histo\r
- Int_t charge = esdTrack->Charge();\r
- Double_t vRecMCTrackHist[17] = { nCrossedRowsTPC, chi2PerCluster, chi2PerClusterITS,ratioCrossedRowsOverFindableClustersTPC, fracClustersTPCShared, b[0], b[1], eta, phi, pt, hasStrangeMother, isFromConversion, isFromMaterial, isPrim, charge, centralityF, chi2(0,0) }; \r
- fRecMCTrackHist->Fill(vRecMCTrackHist);\r
-\r
- if(tpcInnerC) delete tpcInnerC;\r
-}\r
-\r
-//_____________________________________________________________________________\r
-Long64_t AlidNdPtCutAnalysisPbPb::Merge(TCollection* const list) \r
-{\r
- // Merge list of objects (needed by PROOF)\r
-\r
- if (!list)\r
- return 0;\r
-\r
- if (list->IsEmpty())\r
- return 1;\r
-\r
- TIterator* iter = list->MakeIterator();\r
- TObject* obj = 0;\r
- TObjArray* objArrayList = 0;\r
- objArrayList = new TObjArray();\r
-\r
- //TList *collPhysSelection = new TList;\r
-\r
- // collection of generated histograms\r
- Int_t count=0;\r
- while((obj = iter->Next()) != 0) {\r
- AlidNdPtCutAnalysisPbPb* entry = dynamic_cast<AlidNdPtCutAnalysisPbPb*>(obj);\r
- if (entry == 0) continue; \r
- \r
- // event histo\r
- if(GetMergeTHnSparse()) \r
- {\r
- fEventCount->Add(entry->fEventCount);\r
- fRecEventHist->Add(entry->fRecEventHist);\r
- fRecMCEventHist->Add(entry->fRecMCEventHist);\r
- fMCEventHist->Add(entry->fMCEventHist);\r
-\r
- // track histo\r
- fRecMCTrackHist->Add(entry->fRecMCTrackHist);\r
- }\r
-\r
- if (entry->fFolderObj) { objArrayList->Add(entry->fFolderObj); }\r
- \r
- count++;\r
- }\r
- if (fFolderObj) { fFolderObj->Merge(objArrayList); } \r
-\r
- // to signal that track histos were not merged: reset\r
- if (!GetMergeTHnSparse()) \r
- { \r
- // reset content \r
- // save memory\r
- fEventCount->Reset();\r
- fRecEventHist->Reset();\r
- fRecMCEventHist->Reset();\r
- fMCEventHist->Reset();\r
-\r
- // track histo\r
- fRecMCTrackHist->Reset();\r
- }\r
-\r
- // delete\r
- if (objArrayList) delete objArrayList; objArrayList=0;\r
-\r
-return count;\r
-}\r
-\r
-//_____________________________________________________________________________\r
-void AlidNdPtCutAnalysisPbPb::Analyse() \r
-{\r
- //\r
- // Analyse histograms\r
- //\r
- TH1::AddDirectory(kFALSE);\r
- TObjArray *aFolderObj = new TObjArray;\r
- if(!aFolderObj) return;\r
-\r
- TH1D *h1D = 0; \r
- TH2D *h2D = 0; \r
- TH3D *h3D = 0; \r
- \r
- //\r
- // make event level projection\r
- //\r
- //THnSparseF *fEventCount; //-> trig, trig + vertex\r
- //THnSparseF *fRecEventHist; //-> Xv:Yv:Zv:ResZv:Mult\r
- //THnSparseF *fMCEventHist; //-> mcXv:mcYv:mcZv\r
- //THnSparseF *fRecMCEventHist; //-> Xv-mcXv:Yv-mcYv:Zv-mcZv:Mult\r
-\r
- //\r
- // Event counters\r
- //\r
- h2D = (TH2D*)fEventCount->Projection(0,1);\r
- h2D->SetName("trig_vs_trigANDvertex");\r
- h2D->SetTitle("trig_vs_trigANDvertex");\r
- aFolderObj->Add(h2D);\r
-\r
- fEventCount->GetAxis(0)->SetRange(2,2); // triggered\r
- h1D = (TH1D*)fEventCount->Projection(1);\r
- h1D->SetTitle("rec. vertex for triggered events");\r
- h1D->SetName("trigANDvertex");\r
- h1D->SetTitle("trigANDvertex");\r
- aFolderObj->Add(h1D);\r
-\r
- //\r
- // Create rec. event histograms\r
- //\r
- h1D = (TH1D *)fRecEventHist->Projection(0);\r
- h1D->SetName("rec_xv");\r
- h1D->SetTitle("rec_xv");\r
- aFolderObj->Add(h1D);\r
-\r
- h1D = (TH1D *)fRecEventHist->Projection(1);\r
- h1D->SetName("rec_yv");\r
- h1D->SetTitle("rec_yv");\r
- aFolderObj->Add(h1D);\r
-\r
- h1D = (TH1D *)fRecEventHist->Projection(2);\r
- h1D->SetName("rec_zv");\r
- h1D->SetTitle("rec_zv");\r
- aFolderObj->Add(h1D);\r
-\r
- h2D = (TH2D *)fRecEventHist->Projection(0,1);\r
- h2D->SetName("rec_xv_vs_yv");\r
- h2D->SetTitle("rec_xv_vs_yv");\r
- aFolderObj->Add(h2D);\r
-\r
- h2D = (TH2D *)fRecEventHist->Projection(1,2);\r
- h2D->SetName("rec_yv_vs_zv");\r
- h2D->SetTitle("rec_yv_vs_zv");\r
- aFolderObj->Add(h2D);\r
-\r
- h2D = (TH2D *)fRecEventHist->Projection(0,2);\r
- h2D->SetName("rec_xv_vs_zv");\r
- h2D->SetTitle("rec_xv_vs_zv");\r
- aFolderObj->Add(h2D);\r
-\r
- h2D = (TH2D *)fRecEventHist->Projection(3,4);\r
- h2D->SetName("rec_resZv_vs_Mult");\r
- h2D->SetTitle("rec_resZv_vs_Mult");\r
- aFolderObj->Add(h2D);\r
-\r
- //\r
- // MC available\r
- //\r
- if(IsUseMCInfo()) {\r
-\r
- //\r
- // Create mc event histograms\r
- //\r
-\r
- h1D = (TH1D *)fMCEventHist->Projection(0);\r
- h1D->SetName("mc_xv");\r
- h1D->SetTitle("mc_xv");\r
- aFolderObj->Add(h1D);\r
-\r
- h1D = (TH1D *)fMCEventHist->Projection(1);\r
- h1D->SetName("mc_yv");\r
- h1D->SetTitle("mc_yv");\r
- aFolderObj->Add(h1D);\r
-\r
- h1D = (TH1D *)fMCEventHist->Projection(2);\r
- h1D->SetName("mc_zv");\r
- h1D->SetTitle("mc_zv");\r
- aFolderObj->Add(h1D);\r
-\r
-\r
- h2D = (TH2D *)fMCEventHist->Projection(0,1);\r
- h2D->SetName("mc_xv_vs_yv");\r
- h2D->SetTitle("mc_xv_vs_yv");\r
- aFolderObj->Add(h2D);\r
-\r
- h2D = (TH2D *)fMCEventHist->Projection(1,2);\r
- h2D->SetName("mc_yv_vs_zv");\r
- h2D->SetTitle("mc_yv_vs_zv");\r
- aFolderObj->Add(h2D);\r
-\r
- h2D = (TH2D *)fMCEventHist->Projection(0,2);\r
- h2D->SetName("mc_xv_vs_zv");\r
- h2D->SetTitle("mc_xv_vs_zv");\r
- aFolderObj->Add(h2D);\r
-\r
- //\r
- // Create rec-mc event histograms\r
- //\r
- h2D = (TH2D *)fRecMCEventHist->Projection(0,3);\r
- h2D->SetName("rec_mc_deltaXv_vs_mult");\r
- h2D->SetTitle("rec_mc_deltaXv_vs_mult");\r
- aFolderObj->Add(h2D);\r
-\r
- h2D = (TH2D *)fRecMCEventHist->Projection(1,3);\r
- h2D->SetName("rec_mc_deltaYv_vs_mult");\r
- h2D->SetTitle("rec_mc_deltaYv_vs_mult");\r
- aFolderObj->Add(h2D);\r
-\r
- h2D = (TH2D *)fRecMCEventHist->Projection(2,3);\r
- h2D->SetName("rec_mc_deltaZv_vs_mult");\r
- h2D->SetTitle("rec_mc_deltaZv_vs_mult");\r
- aFolderObj->Add(h2D);\r
-\r
- } // end use MC info \r
-\r
- //\r
- // make track level projection\r
- //\r
- // THnSparse track histograms\r
- // nCrossRows:chi2PerClust:chi2PerClustITS:nCrossRows/nFindableClust:fracSharedClust:DCAy:DCAz:eta:phi:pt:hasStrangeMother:isFromConversion:isFromMaterial:isPrim:charge:centr:chi2ToTPCc\r
- //\r
-\r
- // only TH3 histograms\r
- // x : pT : centr\r
- Double_t minNCrossRows = 120.;\r
- Double_t maxChi2PerClust = 4.;\r
- Double_t maxChi2PerClustITS = 36.;\r
- Double_t minNCrossRowsOverFindable = 0.8;\r
- Double_t maxFracSharedClust = 0.4;\r
- Double_t minDCAr = -0.2, maxDCAr = 0.2;\r
- Double_t minDCAz = -2., maxDCAz = 2.;\r
- Double_t minEta = -0.8, maxEta = 0.8;\r
- Double_t maxChi2ToTPCc = 50;\r
-\r
- // only TH3 histograms\r
- // x : pT : centr\r
-\r
- // set all cuts\r
- fRecMCTrackHist->GetAxis(0)->SetRangeUser(minNCrossRows,160);\r
- fRecMCTrackHist->GetAxis(1)->SetRangeUser(0.,maxChi2PerClust-0.0001);\r
- fRecMCTrackHist->GetAxis(2)->SetRangeUser(0.,maxChi2PerClustITS-0.0001);\r
- fRecMCTrackHist->GetAxis(3)->SetRangeUser(minNCrossRowsOverFindable,1.4999);\r
- fRecMCTrackHist->GetAxis(4)->SetRangeUser(0.,maxFracSharedClust-0.0001);\r
- fRecMCTrackHist->GetAxis(5)->SetRangeUser(minDCAr,maxDCAr-0.0001);\r
- fRecMCTrackHist->GetAxis(6)->SetRangeUser(minDCAz,maxDCAz-0.0001);\r
- fRecMCTrackHist->GetAxis(7)->SetRangeUser(minEta,maxEta-0.0001);\r
- fRecMCTrackHist->GetAxis(16)->SetRangeUser(0.,maxChi2ToTPCc-0.0001);\r
-\r
- // nCrossRows \r
- fRecMCTrackHist->GetAxis(0)->SetRangeUser(0,160);\r
-\r
- h3D = (TH3D *)fRecMCTrackHist->Projection(0,9,15);\r
- h3D->SetName("rec_nCrossRows_pT_centr");\r
- h3D->SetTitle("rec_nCrossRows_pT_centr");\r
- aFolderObj->Add(h3D);\r
-\r
-\r
- // Chi2PerClust \r
- fRecMCTrackHist->GetAxis(0)->SetRangeUser(minNCrossRows,160);\r
- fRecMCTrackHist->GetAxis(1)->SetRangeUser(0.,10);\r
-\r
- h3D = (TH3D *)fRecMCTrackHist->Projection(1,9,15);\r
- h3D->SetName("rec_Chi2PerClust_pT_centr");\r
- h3D->SetTitle("rec_Chi2PerClust_pT_centr");\r
- aFolderObj->Add(h3D);\r
-\r
-\r
- // Chi2PerClustITS \r
- fRecMCTrackHist->GetAxis(1)->SetRangeUser(0.,maxChi2PerClust-0.0001);\r
- fRecMCTrackHist->GetAxis(2)->SetRangeUser(0.,100);\r
-\r
- h3D = (TH3D *)fRecMCTrackHist->Projection(2,9,15);\r
- h3D->SetName("rec_Chi2PerClustITS_pT_centr");\r
- h3D->SetTitle("rec_Chi2PerClustITS_pT_centr");\r
- aFolderObj->Add(h3D);\r
-\r
-\r
- // NCrossRowsOverFindable \r
- fRecMCTrackHist->GetAxis(2)->SetRangeUser(0.,maxChi2PerClustITS-0.0001);\r
- fRecMCTrackHist->GetAxis(3)->SetRangeUser(0,1.4999);\r
-\r
- h3D = (TH3D *)fRecMCTrackHist->Projection(3,9,15);\r
- h3D->SetName("rec_NCrossRowsOverFindable_pT_centr");\r
- h3D->SetTitle("rec_NCrossRowsOverFindable_pT_centr");\r
- aFolderObj->Add(h3D);\r
-\r
-\r
- // FracSharedClust\r
- fRecMCTrackHist->GetAxis(3)->SetRangeUser(minNCrossRowsOverFindable,1.4999);\r
- fRecMCTrackHist->GetAxis(4)->SetRangeUser(0.,0.9999);\r
-\r
- h3D = (TH3D *)fRecMCTrackHist->Projection(4,9,15);\r
- h3D->SetName("rec_FracSharedClust_pT_centr");\r
- h3D->SetTitle("rec_FracSharedClust_pT_centr");\r
- aFolderObj->Add(h3D);\r
-\r
-\r
- // DCAr\r
- fRecMCTrackHist->GetAxis(4)->SetRangeUser(0.,maxFracSharedClust-0.0001);\r
- fRecMCTrackHist->GetAxis(5)->SetRangeUser(-10.,9.9999);\r
-\r
- h3D = (TH3D *)fRecMCTrackHist->Projection(5,9,15);\r
- h3D->SetName("rec_DCAr_pT_centr");\r
- h3D->SetTitle("rec_DCAr_pT_centr");\r
- aFolderObj->Add(h3D);\r
-\r
- // DCAz\r
- fRecMCTrackHist->GetAxis(5)->SetRangeUser(minDCAr,maxDCAr-0.0001);\r
- fRecMCTrackHist->GetAxis(6)->SetRangeUser(-10.,9.9999);\r
-\r
- h3D = (TH3D *)fRecMCTrackHist->Projection(6,9,15);\r
- h3D->SetName("rec_DCAz_pT_centr");\r
- h3D->SetTitle("rec_DCAz_pT_centr");\r
- aFolderObj->Add(h3D);\r
-\r
- // chi2ToTPCc\r
- fRecMCTrackHist->GetAxis(6)->SetRangeUser(minDCAz,maxDCAz-0.0001);\r
- fRecMCTrackHist->GetAxis(16)->SetRangeUser(0,100);\r
-\r
- h3D = (TH3D *)fRecMCTrackHist->Projection(16,9,15);\r
- h3D->SetName("rec_chi2ToTPCc_pT_centr");\r
- h3D->SetTitle("rec_chi2ToTPCc_pT_centr");\r
- aFolderObj->Add(h3D);\r
-\r
-\r
- //\r
- // MC available\r
- //\r
- // nCrossRows:chi2PerClust:chi2PerClustITS:nCrossRows/nFindableClust:fracSharedClust:DCAy:DCAz:eta:phi:pt:hasStrangeMother:isFromConversion:isFromMaterial:isPrim:charge:centr:chi2ToTPCc\r
- //\r
-\r
- if(IsUseMCInfo()) {\r
- //\r
- // only TH3 histograms \r
- // x : pT : centr\r
- //\r
- // comes from week decays\r
- //\r
- fRecMCTrackHist->GetAxis(10)->SetRange(2,2);\r
-\r
- // set all cuts\r
- fRecMCTrackHist->GetAxis(0)->SetRangeUser(minNCrossRows,160);\r
- fRecMCTrackHist->GetAxis(1)->SetRangeUser(0.,maxChi2PerClust-0.0001);\r
- fRecMCTrackHist->GetAxis(2)->SetRangeUser(0.,maxChi2PerClustITS-0.0001);\r
- fRecMCTrackHist->GetAxis(3)->SetRangeUser(minNCrossRowsOverFindable,1.4999);\r
- fRecMCTrackHist->GetAxis(4)->SetRangeUser(0.,maxFracSharedClust-0.0001);\r
- fRecMCTrackHist->GetAxis(5)->SetRangeUser(minDCAr,maxDCAr-0.0001);\r
- fRecMCTrackHist->GetAxis(6)->SetRangeUser(minDCAz,maxDCAz-0.0001);\r
- fRecMCTrackHist->GetAxis(7)->SetRangeUser(minEta,maxEta-0.0001);\r
- fRecMCTrackHist->GetAxis(16)->SetRangeUser(0.,maxChi2ToTPCc-0.0001);\r
-\r
- // nCrossRows \r
- fRecMCTrackHist->GetAxis(0)->SetRangeUser(0,160);\r
-\r
- h3D = (TH3D *)fRecMCTrackHist->Projection(0,9,15);\r
- h3D->SetName("rec_sec_strange_nCrossRows_pT_centr");\r
- h3D->SetTitle("rec_sec_strange_nCrossRows_pT_centr");\r
- aFolderObj->Add(h3D);\r
-\r
- // Chi2PerClust \r
- fRecMCTrackHist->GetAxis(0)->SetRangeUser(minNCrossRows,160);\r
- fRecMCTrackHist->GetAxis(1)->SetRangeUser(0.,10);\r
-\r
- h3D = (TH3D *)fRecMCTrackHist->Projection(1,9,15);\r
- h3D->SetName("rec_sec_strange_Chi2PerClust_pT_centr");\r
- h3D->SetTitle("rec_sec_strange_Chi2PerClust_pT_centr");\r
- aFolderObj->Add(h3D);\r
-\r
- // Chi2PerClustITS \r
- fRecMCTrackHist->GetAxis(1)->SetRangeUser(0.,maxChi2PerClust-0.0001);\r
- fRecMCTrackHist->GetAxis(2)->SetRangeUser(0.,100);\r
-\r
- h3D = (TH3D *)fRecMCTrackHist->Projection(2,9,15);\r
- h3D->SetName("rec_sec_strange_Chi2PerClustITS_pT_centr");\r
- h3D->SetTitle("rec_sec_strange_Chi2PerClustITS_pT_centr");\r
- aFolderObj->Add(h3D);\r
-\r
-\r
- // NCrossRowsOverFindable \r
- fRecMCTrackHist->GetAxis(2)->SetRangeUser(0.,maxChi2PerClustITS-0.0001);\r
- fRecMCTrackHist->GetAxis(3)->SetRangeUser(0,1.4999);\r
-\r
- h3D = (TH3D *)fRecMCTrackHist->Projection(3,9,15);\r
- h3D->SetName("rec_sec_strange_NCrossRowsOverFindable_pT_centr");\r
- h3D->SetTitle("rec_sec_strange_NCrossRowsOverFindable_pT_centr");\r
- aFolderObj->Add(h3D);\r
-\r
-\r
- // FracSharedClust\r
- fRecMCTrackHist->GetAxis(3)->SetRangeUser(minNCrossRowsOverFindable,1.4999);\r
- fRecMCTrackHist->GetAxis(4)->SetRangeUser(0.,0.9999);\r
-\r
- h3D = (TH3D *)fRecMCTrackHist->Projection(4,9,15);\r
- h3D->SetName("rec_sec_strange_FracSharedClust_pT_centr");\r
- h3D->SetTitle("rec_sec_strange_FracSharedClust_pT_centr");\r
- aFolderObj->Add(h3D);\r
-\r
-\r
- // DCAr\r
- fRecMCTrackHist->GetAxis(4)->SetRangeUser(0.,maxFracSharedClust-0.0001);\r
- fRecMCTrackHist->GetAxis(5)->SetRangeUser(-10.,9.9999);\r
-\r
- h3D = (TH3D *)fRecMCTrackHist->Projection(5,9,15);\r
- h3D->SetName("rec_sec_strange_DCAr_pT_centr");\r
- h3D->SetTitle("rec_sec_strange_DCAr_pT_centr");\r
- aFolderObj->Add(h3D);\r
-\r
-\r
- // DCAz\r
- fRecMCTrackHist->GetAxis(5)->SetRangeUser(minDCAr,maxDCAr-0.0001);\r
- fRecMCTrackHist->GetAxis(6)->SetRangeUser(-10.,9.9999);\r
-\r
- h3D = (TH3D *)fRecMCTrackHist->Projection(6,9,15);\r
- h3D->SetName("rec_sec_strange_DCAz_pT_centr");\r
- h3D->SetTitle("rec_sec_strange_DCAz_pT_centr");\r
- aFolderObj->Add(h3D);\r
-\r
- // chi2ToTPCc\r
- fRecMCTrackHist->GetAxis(6)->SetRangeUser(minDCAz,maxDCAz-0.0001);\r
- fRecMCTrackHist->GetAxis(16)->SetRangeUser(0,100);\r
-\r
- h3D = (TH3D *)fRecMCTrackHist->Projection(16,9,15);\r
- h3D->SetName("rec_sec_strange_chi2ToTPCc_pT_centr");\r
- h3D->SetTitle("rec_sec_strange_chi2ToTPCc_pT_centr");\r
- aFolderObj->Add(h3D);\r
-\r
-\r
- // only TH3 histograms \r
- // x : pT : centr\r
- //\r
-\r
- // comes from conversion\r
- //\r
- fRecMCTrackHist->GetAxis(10)->SetRange(1,2);\r
- fRecMCTrackHist->GetAxis(11)->SetRange(2,2);\r
-\r
- // set all cuts\r
- fRecMCTrackHist->GetAxis(0)->SetRangeUser(minNCrossRows,160);\r
- fRecMCTrackHist->GetAxis(1)->SetRangeUser(0.,maxChi2PerClust-0.0001);\r
- fRecMCTrackHist->GetAxis(2)->SetRangeUser(0.,maxChi2PerClustITS-0.0001);\r
- fRecMCTrackHist->GetAxis(3)->SetRangeUser(minNCrossRowsOverFindable,1.4999);\r
- fRecMCTrackHist->GetAxis(4)->SetRangeUser(0.,maxFracSharedClust-0.0001);\r
- fRecMCTrackHist->GetAxis(5)->SetRangeUser(minDCAr,maxDCAr-0.0001);\r
- fRecMCTrackHist->GetAxis(6)->SetRangeUser(minDCAz,maxDCAz-0.0001);\r
- fRecMCTrackHist->GetAxis(7)->SetRangeUser(minEta,maxEta-0.0001);\r
- fRecMCTrackHist->GetAxis(16)->SetRangeUser(0.,maxChi2ToTPCc-0.0001);\r
-\r
- // nCrossRows \r
- fRecMCTrackHist->GetAxis(0)->SetRangeUser(0,160);\r
-\r
- h3D = (TH3D *)fRecMCTrackHist->Projection(0,9,15);\r
- h3D->SetName("rec_sec_conversion_nCrossRows_pT_centr");\r
- h3D->SetTitle("rec_sec_conversion_nCrossRows_pT_centr");\r
- aFolderObj->Add(h3D);\r
-\r
- // Chi2PerClust \r
- fRecMCTrackHist->GetAxis(0)->SetRangeUser(minNCrossRows,160);\r
- fRecMCTrackHist->GetAxis(1)->SetRangeUser(0.,10);\r
-\r
- h3D = (TH3D *)fRecMCTrackHist->Projection(1,9,15);\r
- h3D->SetName("rec_sec_conversion_Chi2PerClust_pT_centr");\r
- h3D->SetTitle("rec_sec_conversion_Chi2PerClust_pT_centr");\r
- aFolderObj->Add(h3D);\r
-\r
- // Chi2PerClustITS \r
- fRecMCTrackHist->GetAxis(1)->SetRangeUser(0.,maxChi2PerClust-0.0001);\r
- fRecMCTrackHist->GetAxis(2)->SetRangeUser(0.,100);\r
-\r
- h3D = (TH3D *)fRecMCTrackHist->Projection(2,9,15);\r
- h3D->SetName("rec_sec_conversion_Chi2PerClustITS_pT_centr");\r
- h3D->SetTitle("rec_sec_conversion_Chi2PerClustITS_pT_centr");\r
- aFolderObj->Add(h3D);\r
-\r
-\r
- // NCrossRowsOverFindable \r
- fRecMCTrackHist->GetAxis(2)->SetRangeUser(0.,maxChi2PerClustITS-0.0001);\r
- fRecMCTrackHist->GetAxis(3)->SetRangeUser(0,1.4999);\r
-\r
- h3D = (TH3D *)fRecMCTrackHist->Projection(3,9,15);\r
- h3D->SetName("rec_sec_conversion_NCrossRowsOverFindable_pT_centr");\r
- h3D->SetTitle("rec_sec_conversion_NCrossRowsOverFindable_pT_centr");\r
- aFolderObj->Add(h3D);\r
-\r
-\r
- // FracSharedClust\r
- fRecMCTrackHist->GetAxis(3)->SetRangeUser(minNCrossRowsOverFindable,1.4999);\r
- fRecMCTrackHist->GetAxis(4)->SetRangeUser(0.,0.9999);\r
-\r
- h3D = (TH3D *)fRecMCTrackHist->Projection(4,9,15);\r
- h3D->SetName("rec_sec_conversion_FracSharedClust_pT_centr");\r
- h3D->SetTitle("rec_sec_conversion_FracSharedClust_pT_centr");\r
- aFolderObj->Add(h3D);\r
-\r
-\r
- // DCAr\r
- fRecMCTrackHist->GetAxis(4)->SetRangeUser(0.,maxFracSharedClust-0.0001);\r
- fRecMCTrackHist->GetAxis(5)->SetRangeUser(-10.,9.9999);\r
-\r
- h3D = (TH3D *)fRecMCTrackHist->Projection(5,9,15);\r
- h3D->SetName("rec_sec_conversion_DCAr_pT_centr");\r
- h3D->SetTitle("rec_sec_conversion_DCAr_pT_centr");\r
- aFolderObj->Add(h3D);\r
-\r
-\r
- // DCAz\r
- fRecMCTrackHist->GetAxis(5)->SetRangeUser(minDCAr,maxDCAr-0.0001);\r
- fRecMCTrackHist->GetAxis(6)->SetRangeUser(-10.,9.9999);\r
-\r
- h3D = (TH3D *)fRecMCTrackHist->Projection(6,9,15);\r
- h3D->SetName("rec_sec_conversion_DCAz_pT_centr");\r
- h3D->SetTitle("rec_sec_conversion_DCAz_pT_centr");\r
- aFolderObj->Add(h3D);\r
-\r
- // chi2ToTPCc\r
- fRecMCTrackHist->GetAxis(6)->SetRangeUser(minDCAz,maxDCAz-0.0001);\r
- fRecMCTrackHist->GetAxis(16)->SetRangeUser(0,100);\r
-\r
- h3D = (TH3D *)fRecMCTrackHist->Projection(16,9,15);\r
- h3D->SetName("rec_sec_conversion_chi2ToTPCc_pT_centr");\r
- h3D->SetTitle("rec_sec_conversion_chi2ToTPCc_pT_centr");\r
- aFolderObj->Add(h3D);\r
-\r
-\r
-\r
-\r
-\r
- // comes from interaction with material\r
- //\r
- fRecMCTrackHist->GetAxis(11)->SetRange(1,2);\r
- fRecMCTrackHist->GetAxis(12)->SetRange(2,2);\r
-\r
- // set all cuts\r
- fRecMCTrackHist->GetAxis(0)->SetRangeUser(minNCrossRows,160);\r
- fRecMCTrackHist->GetAxis(1)->SetRangeUser(0.,maxChi2PerClust-0.0001);\r
- fRecMCTrackHist->GetAxis(2)->SetRangeUser(0.,maxChi2PerClustITS-0.0001);\r
- fRecMCTrackHist->GetAxis(3)->SetRangeUser(minNCrossRowsOverFindable,1.4999);\r
- fRecMCTrackHist->GetAxis(4)->SetRangeUser(0.,maxFracSharedClust-0.0001);\r
- fRecMCTrackHist->GetAxis(5)->SetRangeUser(minDCAr,maxDCAr-0.0001);\r
- fRecMCTrackHist->GetAxis(6)->SetRangeUser(minDCAz,maxDCAz-0.0001);\r
- fRecMCTrackHist->GetAxis(7)->SetRangeUser(minEta,maxEta-0.0001);\r
- fRecMCTrackHist->GetAxis(16)->SetRangeUser(0.,maxChi2ToTPCc-0.0001);\r
-\r
- // nCrossRows \r
- fRecMCTrackHist->GetAxis(0)->SetRangeUser(0,160);\r
-\r
- h3D = (TH3D *)fRecMCTrackHist->Projection(0,9,15);\r
- h3D->SetName("rec_sec_material_nCrossRows_pT_centr");\r
- h3D->SetTitle("rec_sec_material_nCrossRows_pT_centr");\r
- aFolderObj->Add(h3D);\r
-\r
- // Chi2PerClust \r
- fRecMCTrackHist->GetAxis(0)->SetRangeUser(minNCrossRows,160);\r
- fRecMCTrackHist->GetAxis(1)->SetRangeUser(0.,10);\r
-\r
- h3D = (TH3D *)fRecMCTrackHist->Projection(1,9,15);\r
- h3D->SetName("rec_sec_material_Chi2PerClust_pT_centr");\r
- h3D->SetTitle("rec_sec_material_Chi2PerClust_pT_centr");\r
- aFolderObj->Add(h3D);\r
-\r
- // Chi2PerClustITS \r
- fRecMCTrackHist->GetAxis(1)->SetRangeUser(0.,maxChi2PerClust-0.0001);\r
- fRecMCTrackHist->GetAxis(2)->SetRangeUser(0.,100);\r
-\r
- h3D = (TH3D *)fRecMCTrackHist->Projection(2,9,15);\r
- h3D->SetName("rec_sec_material_Chi2PerClustITS_pT_centr");\r
- h3D->SetTitle("rec_sec_material_Chi2PerClustITS_pT_centr");\r
- aFolderObj->Add(h3D);\r
-\r
-\r
- // NCrossRowsOverFindable \r
- fRecMCTrackHist->GetAxis(2)->SetRangeUser(0.,maxChi2PerClustITS-0.0001);\r
- fRecMCTrackHist->GetAxis(3)->SetRangeUser(0,1.4999);\r
-\r
- h3D = (TH3D *)fRecMCTrackHist->Projection(3,9,15);\r
- h3D->SetName("rec_sec_material_NCrossRowsOverFindable_pT_centr");\r
- h3D->SetTitle("rec_sec_material_NCrossRowsOverFindable_pT_centr");\r
- aFolderObj->Add(h3D);\r
-\r
-\r
- // FracSharedClust\r
- fRecMCTrackHist->GetAxis(3)->SetRangeUser(minNCrossRowsOverFindable,1.4999);\r
- fRecMCTrackHist->GetAxis(4)->SetRangeUser(0.,0.9999);\r
-\r
- h3D = (TH3D *)fRecMCTrackHist->Projection(4,9,15);\r
- h3D->SetName("rec_sec_material_FracSharedClust_pT_centr");\r
- h3D->SetTitle("rec_sec_material_FracSharedClust_pT_centr");\r
- aFolderObj->Add(h3D);\r
-\r
-\r
- // DCAr\r
- fRecMCTrackHist->GetAxis(4)->SetRangeUser(0.,maxFracSharedClust-0.0001);\r
- fRecMCTrackHist->GetAxis(5)->SetRangeUser(-10.,9.9999);\r
-\r
- h3D = (TH3D *)fRecMCTrackHist->Projection(5,9,15);\r
- h3D->SetName("rec_sec_material_DCAr_pT_centr");\r
- h3D->SetTitle("rec_sec_material_DCAr_pT_centr");\r
- aFolderObj->Add(h3D);\r
-\r
-\r
- // DCAz\r
- fRecMCTrackHist->GetAxis(5)->SetRangeUser(minDCAr,maxDCAr-0.0001);\r
- fRecMCTrackHist->GetAxis(6)->SetRangeUser(-10.,9.9999);\r
-\r
- h3D = (TH3D *)fRecMCTrackHist->Projection(6,9,15);\r
- h3D->SetName("rec_sec_material_DCAz_pT_centr");\r
- h3D->SetTitle("rec_sec_material_DCAz_pT_centr");\r
- aFolderObj->Add(h3D);\r
-\r
- // chi2ToTPCc\r
- fRecMCTrackHist->GetAxis(6)->SetRangeUser(minDCAz,maxDCAz-0.0001);\r
- fRecMCTrackHist->GetAxis(16)->SetRangeUser(0,100);\r
-\r
- h3D = (TH3D *)fRecMCTrackHist->Projection(16,9,15);\r
- h3D->SetName("rec_sec_material_chi2ToTPCc_pT_centr");\r
- h3D->SetTitle("rec_sec_material_chi2ToTPCc_pT_centr");\r
- aFolderObj->Add(h3D);\r
-\r
-\r
-\r
- //\r
- // only TH3 histograms \r
- // x : pT : centr\r
- //\r
- // secondaries defined as AliStack::IsPhysicalPrimary()\r
- //\r
- fRecMCTrackHist->GetAxis(12)->SetRange(1,2);\r
- fRecMCTrackHist->GetAxis(13)->SetRange(1,1);\r
-\r
- // set all cuts\r
- fRecMCTrackHist->GetAxis(0)->SetRangeUser(minNCrossRows,160);\r
- fRecMCTrackHist->GetAxis(1)->SetRangeUser(0.,maxChi2PerClust-0.0001);\r
- fRecMCTrackHist->GetAxis(2)->SetRangeUser(0.,maxChi2PerClustITS-0.0001);\r
- fRecMCTrackHist->GetAxis(3)->SetRangeUser(minNCrossRowsOverFindable,1.4999);\r
- fRecMCTrackHist->GetAxis(4)->SetRangeUser(0.,maxFracSharedClust-0.0001);\r
- fRecMCTrackHist->GetAxis(5)->SetRangeUser(minDCAr,maxDCAr-0.0001);\r
- fRecMCTrackHist->GetAxis(6)->SetRangeUser(minDCAz,maxDCAz-0.0001);\r
- fRecMCTrackHist->GetAxis(7)->SetRangeUser(minEta,maxEta-0.0001);\r
- fRecMCTrackHist->GetAxis(16)->SetRangeUser(0.,maxChi2ToTPCc-0.0001);\r
-\r
- // nCrossRows \r
- fRecMCTrackHist->GetAxis(0)->SetRangeUser(0,160);\r
-\r
- h3D = (TH3D *)fRecMCTrackHist->Projection(0,9,15);\r
- h3D->SetName("rec_sec_nCrossRows_pT_centr");\r
- h3D->SetTitle("rec_sec_nCrossRows_pT_centr");\r
- aFolderObj->Add(h3D);\r
-\r
- // Chi2PerClust \r
- fRecMCTrackHist->GetAxis(0)->SetRangeUser(minNCrossRows,160);\r
- fRecMCTrackHist->GetAxis(1)->SetRangeUser(0.,10);\r
-\r
- h3D = (TH3D *)fRecMCTrackHist->Projection(1,9,15);\r
- h3D->SetName("rec_sec_Chi2PerClust_pT_centr");\r
- h3D->SetTitle("rec_sec_Chi2PerClust_pT_centr");\r
- aFolderObj->Add(h3D);\r
-\r
- // Chi2PerClustITS \r
- fRecMCTrackHist->GetAxis(1)->SetRangeUser(0.,maxChi2PerClust-0.0001);\r
- fRecMCTrackHist->GetAxis(2)->SetRangeUser(0.,100);\r
-\r
- h3D = (TH3D *)fRecMCTrackHist->Projection(2,9,15);\r
- h3D->SetName("rec_sec_Chi2PerClustITS_pT_centr");\r
- h3D->SetTitle("rec_sec_Chi2PerClustITS_pT_centr");\r
- aFolderObj->Add(h3D);\r
-\r
-\r
- // NCrossRowsOverFindable \r
- fRecMCTrackHist->GetAxis(2)->SetRangeUser(0.,maxChi2PerClustITS-0.0001);\r
- fRecMCTrackHist->GetAxis(3)->SetRangeUser(0,1.4999);\r
-\r
- h3D = (TH3D *)fRecMCTrackHist->Projection(3,9,15);\r
- h3D->SetName("rec_sec_NCrossRowsOverFindable_pT_centr");\r
- h3D->SetTitle("rec_sec_NCrossRowsOverFindable_pT_centr");\r
- aFolderObj->Add(h3D);\r
-\r
-\r
- // FracSharedClust\r
- fRecMCTrackHist->GetAxis(3)->SetRangeUser(minNCrossRowsOverFindable,1.4999);\r
- fRecMCTrackHist->GetAxis(4)->SetRangeUser(0.,0.9999);\r
-\r
- h3D = (TH3D *)fRecMCTrackHist->Projection(4,9,15);\r
- h3D->SetName("rec_sec_FracSharedClust_pT_centr");\r
- h3D->SetTitle("rec_sec_FracSharedClust_pT_centr");\r
- aFolderObj->Add(h3D);\r
-\r
-\r
- // DCAr\r
- fRecMCTrackHist->GetAxis(4)->SetRangeUser(0.,maxFracSharedClust-0.0001);\r
- fRecMCTrackHist->GetAxis(5)->SetRangeUser(-10.,9.9999);\r
-\r
- h3D = (TH3D *)fRecMCTrackHist->Projection(5,9,15);\r
- h3D->SetName("rec_sec_DCAr_pT_centr");\r
- h3D->SetTitle("rec_sec_DCAr_pT_centr");\r
- aFolderObj->Add(h3D);\r
-\r
-\r
- // DCAz\r
- fRecMCTrackHist->GetAxis(5)->SetRangeUser(minDCAr,maxDCAr-0.0001);\r
- fRecMCTrackHist->GetAxis(6)->SetRangeUser(-10.,9.9999);\r
-\r
- h3D = (TH3D *)fRecMCTrackHist->Projection(6,9,15);\r
- h3D->SetName("rec_sec_DCAz_pT_centr");\r
- h3D->SetTitle("rec_sec_DCAz_pT_centr");\r
- aFolderObj->Add(h3D);\r
-\r
- // chi2ToTPCc\r
- fRecMCTrackHist->GetAxis(6)->SetRangeUser(minDCAz,maxDCAz-0.0001);\r
- fRecMCTrackHist->GetAxis(16)->SetRangeUser(0,100);\r
-\r
- h3D = (TH3D *)fRecMCTrackHist->Projection(16,9,15);\r
- h3D->SetName("rec_sec_chi2ToTPCc_pT_centr");\r
- h3D->SetTitle("rec_sec_chi2ToTPCc_pT_centr");\r
- aFolderObj->Add(h3D);\r
-\r
- //\r
- // only TH3 histograms \r
- // x : pT : centr\r
- //\r
- // primaries defined as AliStack::IsPhysicalPrimary()\r
- //\r
- fRecMCTrackHist->GetAxis(12)->SetRange(1,2);\r
- fRecMCTrackHist->GetAxis(13)->SetRange(2,2);\r
-\r
- // set all cuts\r
- fRecMCTrackHist->GetAxis(0)->SetRangeUser(minNCrossRows,160);\r
- fRecMCTrackHist->GetAxis(1)->SetRangeUser(0.,maxChi2PerClust-0.0001);\r
- fRecMCTrackHist->GetAxis(2)->SetRangeUser(0.,maxChi2PerClustITS-0.0001);\r
- fRecMCTrackHist->GetAxis(3)->SetRangeUser(minNCrossRowsOverFindable,1.4999);\r
- fRecMCTrackHist->GetAxis(4)->SetRangeUser(0.,maxFracSharedClust-0.0001);\r
- fRecMCTrackHist->GetAxis(5)->SetRangeUser(minDCAr,maxDCAr-0.0001);\r
- fRecMCTrackHist->GetAxis(6)->SetRangeUser(minDCAz,maxDCAz-0.0001);\r
- fRecMCTrackHist->GetAxis(7)->SetRangeUser(minEta,maxEta-0.0001);\r
- fRecMCTrackHist->GetAxis(16)->SetRangeUser(0.,maxChi2ToTPCc-0.0001);\r
-\r
- // nCrossRows \r
- fRecMCTrackHist->GetAxis(0)->SetRangeUser(0,160);\r
-\r
- h3D = (TH3D *)fRecMCTrackHist->Projection(0,9,15);\r
- h3D->SetName("rec_prim_nCrossRows_pT_centr");\r
- h3D->SetTitle("rec_prim_nCrossRows_pT_centr");\r
- aFolderObj->Add(h3D);\r
-\r
- // Chi2PerClust \r
- fRecMCTrackHist->GetAxis(0)->SetRangeUser(minNCrossRows,160);\r
- fRecMCTrackHist->GetAxis(1)->SetRangeUser(0.,10);\r
-\r
- h3D = (TH3D *)fRecMCTrackHist->Projection(1,9,15);\r
- h3D->SetName("rec_prim_Chi2PerClust_pT_centr");\r
- h3D->SetTitle("rec_prim_Chi2PerClust_pT_centr");\r
- aFolderObj->Add(h3D);\r
-\r
- // Chi2PerClustITS \r
- fRecMCTrackHist->GetAxis(1)->SetRangeUser(0.,maxChi2PerClust-0.0001);\r
- fRecMCTrackHist->GetAxis(2)->SetRangeUser(0.,100);\r
-\r
- h3D = (TH3D *)fRecMCTrackHist->Projection(2,9,15);\r
- h3D->SetName("rec_prim_Chi2PerClustITS_pT_centr");\r
- h3D->SetTitle("rec_prim_Chi2PerClustITS_pT_centr");\r
- aFolderObj->Add(h3D);\r
-\r
-\r
- // NCrossRowsOverFindable \r
- fRecMCTrackHist->GetAxis(2)->SetRangeUser(0.,maxChi2PerClustITS-0.0001);\r
- fRecMCTrackHist->GetAxis(3)->SetRangeUser(0,1.4999);\r
-\r
- h3D = (TH3D *)fRecMCTrackHist->Projection(3,9,15);\r
- h3D->SetName("rec_prim_NCrossRowsOverFindable_pT_centr");\r
- h3D->SetTitle("rec_prim_NCrossRowsOverFindable_pT_centr");\r
- aFolderObj->Add(h3D);\r
-\r
-\r
- // FracSharedClust\r
- fRecMCTrackHist->GetAxis(3)->SetRangeUser(minNCrossRowsOverFindable,1.4999);\r
- fRecMCTrackHist->GetAxis(4)->SetRangeUser(0.,0.9999);\r
-\r
- h3D = (TH3D *)fRecMCTrackHist->Projection(4,9,15);\r
- h3D->SetName("rec_prim_FracSharedClust_pT_centr");\r
- h3D->SetTitle("rec_prim_FracSharedClust_pT_centr");\r
- aFolderObj->Add(h3D);\r
-\r
-\r
- // DCAr\r
- fRecMCTrackHist->GetAxis(4)->SetRangeUser(0.,maxFracSharedClust-0.0001);\r
- fRecMCTrackHist->GetAxis(5)->SetRangeUser(-10.,9.9999);\r
-\r
- h3D = (TH3D *)fRecMCTrackHist->Projection(5,9,15);\r
- h3D->SetName("rec_prim_DCAr_pT_centr");\r
- h3D->SetTitle("rec_prim_DCAr_pT_centr");\r
- aFolderObj->Add(h3D);\r
-\r
-\r
- // DCAz\r
- fRecMCTrackHist->GetAxis(5)->SetRangeUser(minDCAr,maxDCAr-0.0001);\r
- fRecMCTrackHist->GetAxis(6)->SetRangeUser(-10.,9.9999);\r
-\r
- h3D = (TH3D *)fRecMCTrackHist->Projection(6,9,15);\r
- h3D->SetName("rec_prim_DCAz_pT_centr");\r
- h3D->SetTitle("rec_prim_DCAz_pT_centr");\r
- aFolderObj->Add(h3D);\r
-\r
- // chi2ToTPCc\r
- fRecMCTrackHist->GetAxis(6)->SetRangeUser(minDCAz,maxDCAz-0.0001);\r
- fRecMCTrackHist->GetAxis(16)->SetRangeUser(0,100);\r
-\r
- h3D = (TH3D *)fRecMCTrackHist->Projection(16,9,15);\r
- h3D->SetName("rec_prim_chi2ToTPCc_pT_centr");\r
- h3D->SetTitle("rec_prim_chi2ToTPCc_pT_centr");\r
- aFolderObj->Add(h3D);\r
-\r
- }\r
-\r
- // export objects to analysis folder\r
- fAnalysisFolder = ExportToFolder(aFolderObj);\r
- if(!fAnalysisFolder) {\r
- if(aFolderObj) delete aFolderObj;\r
- return;\r
- }\r
-\r
- // Reset setings\r
- \r
- for(Int_t ndim = 0; ndim < fEventCount->GetNdimensions(); ndim++) {\r
- fEventCount->GetAxis(ndim)->SetRange(1,fEventCount->GetAxis(ndim)->GetNbins());\r
- }\r
-\r
- for(Int_t ndim = 0; ndim < fMCEventHist->GetNdimensions(); ndim++) {\r
- fMCEventHist->GetAxis(ndim)->SetRange(1,fMCEventHist->GetAxis(ndim)->GetNbins());\r
- }\r
-\r
- for(Int_t ndim = 0; ndim < fRecMCEventHist->GetNdimensions(); ndim++) {\r
- fRecMCEventHist->GetAxis(ndim)->SetRange(1,fRecMCEventHist->GetAxis(ndim)->GetNbins());\r
- }\r
-\r
- for(Int_t ndim = 0; ndim < fRecMCTrackHist->GetNdimensions(); ndim++) \r
- {\r
- fRecMCTrackHist->GetAxis(ndim)->SetRange(1,fRecMCTrackHist->GetAxis(ndim)->GetNbins());\r
- }\r
-\r
- // delete only TObjArray\r
- //if(aFolderObj) delete aFolderObj;\r
-\r
- if (fFolderObj) delete fFolderObj;\r
- fFolderObj = aFolderObj;\r
- aFolderObj=0;\r
-}\r
-\r
-//_____________________________________________________________________________\r
-TFolder* AlidNdPtCutAnalysisPbPb::ExportToFolder(TObjArray * const array) \r
-{\r
- // recreate folder avery time and export objects to new one\r
- //\r
- AlidNdPtCutAnalysisPbPb * comp=this;\r
- TFolder *folder = comp->GetAnalysisFolder();\r
-\r
- TString name, title;\r
- TFolder *newFolder = 0;\r
- Int_t i = 0;\r
- Int_t size = array->GetSize();\r
-\r
- if(folder) { \r
- // get name and title from old folder\r
- name = folder->GetName(); \r
- title = folder->GetTitle(); \r
-\r
- // delete old one\r
- delete folder;\r
-\r
- // create new one\r
- newFolder = CreateFolder(name.Data(),title.Data());\r
- newFolder->SetOwner();\r
-\r
- // add objects to folder\r
- while(i < size) {\r
- newFolder->Add(array->At(i));\r
- i++;\r
- }\r
- }\r
-\r
-return newFolder;\r
-}\r
-\r
-//_____________________________________________________________________________\r
-TFolder* AlidNdPtCutAnalysisPbPb::CreateFolder(TString name,TString title) { \r
-// create folder for analysed histograms\r
-//\r
-TFolder *folder = 0;\r
- folder = new TFolder(name.Data(),title.Data());\r
-\r
- return folder;\r
-}\r
+/**************************************************************************
+ * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
+ * *
+ * Author: The ALICE Off-line Project. *
+ * Contributors are mentioned in the code where appropriate. *
+ * *
+ * Permission to use, copy, modify and distribute this software and its *
+ * documentation strictly for non-commercial purposes is hereby granted *
+ * without fee, provided that the above copyright notice appears in all *
+ * copies and that both the copyright notice and this permission notice *
+ * appear in the supporting documentation. The authors make no claims *
+ * about the suitability of this software for any purpose. It is *
+ * provided "as is" without express or implied warranty. *
+ **************************************************************************/
+//------------------------------------------------------------------------------
+// AlidNdPtCutAnalysisPbPb class.
+//
+// a. functionality:
+// - fills generic cut histograms
+// - generates cuts (selection criteria)
+//
+// b. data members:
+// - generic cut histograms
+// - control histograms
+//
+// Author: J.Otwinowski 04/11/2008
+//------------------------------------------------------------------------------
+#include "TH1.h"
+#include "TH2.h"
+#include "TH3.h"
+#include "TMatrixD.h"
+
+#include <TPDGCode.h>
+
+#include "AliHeader.h"
+#include "AliGenEventHeader.h"
+#include "AliInputEventHandler.h"
+#include "AliCentrality.h"
+#include "AliStack.h"
+#include "AliESDEvent.h"
+#include "AliMCEvent.h"
+#include "AliESDtrackCuts.h"
+#include "AliLog.h"
+#include "AliTracker.h"
+
+#include "AlidNdPtEventCuts.h"
+#include "AlidNdPtAcceptanceCuts.h"
+#include "AlidNdPtBackgroundCuts.h"
+#include "AlidNdPtAnalysis.h"
+#include "AliPhysicsSelection.h"
+
+#include "AliPWG0Helper.h"
+#include "AlidNdPtHelper.h"
+#include "AlidNdPtCutAnalysisPbPb.h"
+
+using namespace std;
+
+ClassImp(AlidNdPtCutAnalysisPbPb)
+
+//_____________________________________________________________________________
+ AlidNdPtCutAnalysisPbPb::AlidNdPtCutAnalysisPbPb(): AlidNdPt(),
+ fAnalysisFolder(0),
+ fEventCount(0),
+ fRecEventHist(0),
+ fMCEventHist(0),
+ fRecMCEventHist(0),
+ fRecMCTrackHist(0),
+ fCentralityEstimator(0),
+ fFolderObj(0)
+{
+ // default constructor
+ Init();
+}
+
+//_____________________________________________________________________________
+AlidNdPtCutAnalysisPbPb::AlidNdPtCutAnalysisPbPb(Char_t* name, Char_t* title): AlidNdPt(name,title),
+ fAnalysisFolder(0),
+ fEventCount(0),
+ fRecEventHist(0),
+ fMCEventHist(0),
+ fRecMCEventHist(0),
+ fRecMCTrackHist(0),
+ fCentralityEstimator(0),
+ fFolderObj(0)
+{
+ // constructor
+ Init();
+}
+
+//_____________________________________________________________________________
+AlidNdPtCutAnalysisPbPb::~AlidNdPtCutAnalysisPbPb() {
+ //
+ if(fEventCount) delete fEventCount; fEventCount=0;
+ if(fRecEventHist) delete fRecEventHist; fRecEventHist=0;
+ if(fMCEventHist) delete fMCEventHist; fMCEventHist=0;
+ if(fRecMCEventHist) delete fRecMCEventHist; fRecMCEventHist=0;
+ if(fRecMCTrackHist) delete fRecMCTrackHist; fRecMCTrackHist=0;
+
+ if(fAnalysisFolder) delete fAnalysisFolder; fAnalysisFolder=0;
+ if(fFolderObj) delete fFolderObj; fFolderObj=0;
+}
+
+//_____________________________________________________________________________
+void AlidNdPtCutAnalysisPbPb::Init(){
+ //
+ // Init histograms
+ //
+ //const Int_t ptNbins = 58;
+ //const Double_t ptMin = 0.;
+ //const Double_t ptMax = 20.;
+ /*
+ Double_t binsPt[ptNbins+1] = {0.,0.05,0.1,0.15,0.2,0.25,0.3,0.35,0.4,0.45,0.5,0.55,0.6,0.65,0.7,0.75,0.8,0.85,0.9,0.95,1.0,1.1,1.2,1.3,1.4,1.5,1.6,1.7,1.8,1.9,2.0,2.2,2.4,2.6,2.8,3.0,3.2,3.4,3.6,3.8,4.0,4.5,5.0,5.5,6.0,6.5,7.0,7.5,8.0,9.0,10.0,11.0,12.0,13.0,14.0,15.0,16.0,18.0, 20.};
+ */
+
+ // set pt bins
+ const Int_t ptNbins = 50;
+ const Double_t ptMin = 1.e-2, ptMax = 100.;
+ Double_t *binsPt = CreateLogAxis(ptNbins,ptMin,ptMax);
+
+ // centrality bins
+ const Int_t centrNbins = 3;
+ const Double_t centrMin = 0, centrMax = 1;
+ Double_t binsCentr[centrNbins+1] = {0.0, 20., 50., 100.};
+
+ //
+ Int_t binsEventCount[2]={2,2};
+ Double_t minEventCount[2]={0,0};
+ Double_t maxEventCount[2]={2,2};
+ fEventCount = new THnSparseF("fEventCount","trig vs trig+vertex",2,binsEventCount,minEventCount,maxEventCount);
+ fEventCount->GetAxis(0)->SetTitle("trig");
+ fEventCount->GetAxis(1)->SetTitle("trig+vert");
+ fEventCount->Sumw2();
+
+ //Xv:Yv:Zv:ResZv:Mult
+ Double_t kFact = 0.1;
+
+ Int_t binsRecEventHist[5]={80,80,100,80,150};
+ Double_t minRecEventHist[5]={-1.*kFact,-1.*kFact,-35.,0.,0.};
+ Double_t maxRecEventHist[5]={1.*kFact,1.*kFact,35.,10.,3000.};
+ fRecEventHist = new THnSparseF("fRecEventHist","Xv:Yv:Zv:ResZv:Mult",5,binsRecEventHist,minRecEventHist,maxRecEventHist);
+ fRecEventHist->GetAxis(0)->SetTitle("Xv (cm)");
+ fRecEventHist->GetAxis(1)->SetTitle("Yv (cm)");
+ fRecEventHist->GetAxis(2)->SetTitle("Zv (cm)");
+ fRecEventHist->GetAxis(3)->SetTitle("ResZv (cm)");
+ fRecEventHist->GetAxis(4)->SetTitle("Mult");
+ fRecEventHist->Sumw2();
+
+ //Xv:Yv:Zv
+ Int_t binsMCEventHist[3]={80,80,100};
+ Double_t minMCEventHist[3]={-0.1,-0.1,-35.};
+ Double_t maxMCEventHist[3]={0.1,0.1,35.};
+ fMCEventHist = new THnSparseF("fMCEventHist","mcXv:mcYv:mcZv",3,binsMCEventHist,minMCEventHist,maxMCEventHist);
+ fMCEventHist->GetAxis(0)->SetTitle("mcXv (cm)");
+ fMCEventHist->GetAxis(1)->SetTitle("mcYv (cm)");
+ fMCEventHist->GetAxis(2)->SetTitle("mcZv (cm)");
+ fMCEventHist->Sumw2();
+
+ //Xv-mcXv:Yv-mcYv:Zv-mcZv:Mult
+ Int_t binsRecMCEventHist[4]={100,100,100,150};
+ Double_t minRecMCEventHist[4]={-1.0*kFact,-1.0*kFact,-1.0*kFact,0.};
+ Double_t maxRecMCEventHist[4]={1.0*kFact,1.0*kFact,1.0*kFact,3000.};
+ fRecMCEventHist = new THnSparseF("fRecMCEventHist","Xv-mcXv:Yv-mcYv:Zv-mcZv:Mult",4,binsRecMCEventHist,minRecMCEventHist,maxRecMCEventHist);
+ fRecMCEventHist->GetAxis(0)->SetTitle("Xv-mcXv (cm)");
+ fRecMCEventHist->GetAxis(1)->SetTitle("Yv-mcYv (cm)");
+ fRecMCEventHist->GetAxis(2)->SetTitle("Zv-mcZv (cm)");
+ fRecMCEventHist->GetAxis(3)->SetTitle("Mult");
+ fRecMCEventHist->Sumw2();
+
+ //
+ // THnSparse track histograms
+ //
+ //
+ //
+ //nCrossRows:chi2PerClust:chi2PerClustITS:nCrossRows/nFindableClust:fracSharedClust:DCAy:DCAz:eta:phi:pt:hasStrangeMother:isFromConversion:isFromMaterial:isPrim:charge:centr:chi2ToTPCc
+ Int_t binsRecMCTrackHist[17]= {160, 10, 70, 30, 20, 50, 50, 20, 90, ptNbins, 2, 2, 2, 2, 3, centrNbins,100};
+ Double_t minRecMCTrackHist[17]={0., 0., 0., 0., 0.,-0.5,-2.5,-1.0, 0., ptMin, 0., 0., 0., 0.,-1., centrMin, 0.};
+ Double_t maxRecMCTrackHist[17]={160.,10.,70.,1.5, 1., 0.5, 2.5, 1.0, 2.*TMath::Pi(), ptMax, 2., 2., 2., 2., 2., centrMax, 100.};
+
+ fRecMCTrackHist = new THnSparseF("fRecMCTrackHist","nCrossRows:chi2PerClust:chi2PerClustITS:nCrossRows/nFindableClust:fracSharedClust:DCAy:DCAz:eta:phi:pt:hasStrangeMother:isFromConversion:isFromMaterial:isPrim:charge:centr:chi2ToTPCc",17,binsRecMCTrackHist,minRecMCTrackHist,maxRecMCTrackHist);
+ fRecMCTrackHist->SetBinEdges(9,binsPt);
+ fRecMCTrackHist->SetBinEdges(15,binsCentr);
+
+ fRecMCTrackHist->GetAxis(0)->SetTitle("nCrossRows");
+ fRecMCTrackHist->GetAxis(1)->SetTitle("chi2PerClust");
+ fRecMCTrackHist->GetAxis(2)->SetTitle("chi2PerClustITS");
+ fRecMCTrackHist->GetAxis(3)->SetTitle("nCrossRows/nFindableClust");
+ fRecMCTrackHist->GetAxis(4)->SetTitle("fracSharedClust");
+ fRecMCTrackHist->GetAxis(5)->SetTitle("DCAy (cm)");
+ fRecMCTrackHist->GetAxis(6)->SetTitle("DCAz (cm)");
+ fRecMCTrackHist->GetAxis(7)->SetTitle("#eta");
+ fRecMCTrackHist->GetAxis(8)->SetTitle("#phi (rad)");
+ fRecMCTrackHist->GetAxis(9)->SetTitle("p_{T} (GeV/c)");
+ fRecMCTrackHist->GetAxis(10)->SetTitle("hasStrangeMother");
+ fRecMCTrackHist->GetAxis(11)->SetTitle("isFromConversion");
+ fRecMCTrackHist->GetAxis(12)->SetTitle("isFromMaterial");
+ fRecMCTrackHist->GetAxis(13)->SetTitle("isPrim");
+ fRecMCTrackHist->GetAxis(14)->SetTitle("charge");
+ fRecMCTrackHist->GetAxis(15)->SetTitle("centrality");
+ fRecMCTrackHist->GetAxis(16)->SetTitle("chi2ToTPCc");
+ fRecMCTrackHist->Sumw2();
+
+ // init output folder
+ fAnalysisFolder = CreateFolder("folderdNdPt","Analysis dNdPt Folder");
+
+}
+
+//_____________________________________________________________________________
+void AlidNdPtCutAnalysisPbPb::Process(AliESDEvent *const esdEvent, AliMCEvent * const mcEvent)
+{
+ //
+ // Process real and/or simulated events
+ //
+ if(!esdEvent) {
+ AliDebug(AliLog::kError, "esdEvent not available");
+ return;
+ }
+
+ // get selection cuts
+ AlidNdPtEventCuts *evtCuts = GetEventCuts();
+ AlidNdPtAcceptanceCuts *accCuts = GetAcceptanceCuts();
+ AliESDtrackCuts *esdTrackCuts = GetTrackCuts();
+
+ if(!evtCuts || !accCuts || !esdTrackCuts) {
+ AliDebug(AliLog::kError, "cuts not available");
+ return;
+ }
+
+ // trigger selection
+
+ Bool_t isEventTriggered = kTRUE;
+ AliPhysicsSelection *physicsSelection = NULL;
+ AliTriggerAnalysis* triggerAnalysis = NULL;
+
+ //
+ AliInputEventHandler* inputHandler = (AliInputEventHandler*) AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler();
+ if (!inputHandler)
+ {
+ Printf("ERROR: Could not receive input handler");
+ return;
+ }
+
+ if(evtCuts->IsTriggerRequired())
+ {
+ // always MB
+ isEventTriggered = inputHandler->IsEventSelected() & AliVEvent::kMB;
+
+ physicsSelection = static_cast<AliPhysicsSelection*> (inputHandler->GetEventSelection());
+ if(!physicsSelection) return;
+ //SetPhysicsTriggerSelection(physicsSelection);
+
+ if (isEventTriggered && (GetTrigger() == AliTriggerAnalysis::kV0AND)) {
+ // set trigger (V0AND)
+ triggerAnalysis = physicsSelection->GetTriggerAnalysis();
+ if(!triggerAnalysis) return;
+ isEventTriggered = triggerAnalysis->IsOfflineTriggerFired(esdEvent, GetTrigger());
+ }
+ }
+
+ // centrality determination
+ Float_t centralityF = -1;
+ AliCentrality *esdCentrality = esdEvent->GetCentrality();
+ centralityF = esdCentrality->GetCentralityPercentile(fCentralityEstimator.Data());
+
+ // use MC information
+ AliHeader* header = 0;
+ AliGenEventHeader* genHeader = 0;
+ AliStack* stack = 0;
+ TArrayF vtxMC(3);
+
+ Int_t multMCTrueTracks = 0;
+ if(IsUseMCInfo())
+ {
+ //
+ if(!mcEvent) {
+ AliDebug(AliLog::kError, "mcEvent not available");
+ return;
+ }
+ // get MC event header
+ header = mcEvent->Header();
+ if (!header) {
+ AliDebug(AliLog::kError, "Header not available");
+ return;
+ }
+ // MC particle stack
+ stack = mcEvent->Stack();
+ if (!stack) {
+ AliDebug(AliLog::kError, "Stack not available");
+ return;
+ }
+
+ // get MC vertex
+ genHeader = header->GenEventHeader();
+ if (!genHeader) {
+ AliDebug(AliLog::kError, "Could not retrieve genHeader from Header");
+ return;
+ }
+ genHeader->PrimaryVertex(vtxMC);
+
+ Double_t vMCEventHist[3]={vtxMC[0],vtxMC[1],vtxMC[2]};
+ fMCEventHist->Fill(vMCEventHist);
+
+ // multipliticy of all MC primary tracks
+ // in Zv, pt and eta ranges)
+ multMCTrueTracks = AlidNdPtHelper::GetMCTrueTrackMult(mcEvent,evtCuts,accCuts);
+
+ } // end bUseMC
+
+ // get reconstructed vertex
+ const AliESDVertex* vtxESD = 0;
+ if(GetAnalysisMode() == AlidNdPtHelper::kTPC) {
+ vtxESD = esdEvent->GetPrimaryVertexTPC();
+ }
+ else if(GetAnalysisMode() == AlidNdPtHelper::kTPCITS) {
+ vtxESD = esdEvent->GetPrimaryVertexTracks();
+ }
+ else {
+ return;
+ }
+
+ if(!vtxESD) return;
+
+ Bool_t isEventOK = evtCuts->AcceptEvent(esdEvent,mcEvent,vtxESD);
+ //printf("isEventOK %d, isEventTriggered %d \n",isEventOK, isEventTriggered);
+ //printf("GetAnalysisMode() %d \n",GetAnalysisMode());
+
+ TObjArray *allChargedTracks=0;
+ Int_t multAll=0;
+
+ //
+ // event counter
+ //
+ //printf("isEventOK %d, isEventTriggered %d \n",isEventOK,isEventTriggered);
+
+ Bool_t isTrigAndVertex = isEventTriggered && isEventOK;
+ Double_t vEventCount[2] = { isEventTriggered, isTrigAndVertex};
+ fEventCount->Fill(vEventCount);
+
+ // check event cuts
+ if(isEventOK && isEventTriggered)
+ {
+ // get all charged tracks
+ allChargedTracks = AlidNdPtHelper::GetAllChargedTracks(esdEvent,GetAnalysisMode());
+ if(!allChargedTracks) return;
+
+ Int_t entries = allChargedTracks->GetEntries();
+ for(Int_t i=0; i<entries;++i)
+ {
+ AliESDtrack *track = (AliESDtrack*)allChargedTracks->At(i);
+ if(!track) continue;
+
+ if(!esdTrackCuts->AcceptTrack(track)) continue;
+ if(!accCuts->AcceptTrack(track)) continue;
+
+ FillHistograms(esdEvent, track, stack, centralityF);
+ multAll++;
+
+ }
+
+ Double_t vRecEventHist[5] = {vtxESD->GetXv(),vtxESD->GetYv(),vtxESD->GetZv(),vtxESD->GetZRes(),multAll};
+ fRecEventHist->Fill(vRecEventHist);
+
+ if(IsUseMCInfo()) {
+ Double_t vRecMCEventHist[5] = {vtxESD->GetXv()-vtxMC[0],vtxESD->GetYv()-vtxMC[1],vtxESD->GetZv()-vtxMC[2],multAll};
+ fRecMCEventHist->Fill(vRecMCEventHist);
+ }
+ }
+
+ if(allChargedTracks) delete allChargedTracks; allChargedTracks = 0;
+
+}
+
+//_____________________________________________________________________________
+void AlidNdPtCutAnalysisPbPb::FillHistograms(AliESDEvent *const esdEvent, AliESDtrack *const esdTrack, AliStack *const stack, Float_t centralityF) const
+{
+ //
+ // Fill ESD track and MC histograms
+ //
+ if(!esdTrack) return;
+ if(esdTrack->Charge() == 0.) return;
+
+ Float_t pt = esdTrack->Pt();
+ Float_t eta = esdTrack->Eta();
+ Float_t phi = esdTrack->Phi();
+
+ Int_t nClust = 0;
+ if(GetAnalysisMode() == AlidNdPtHelper::kTPC) {
+ nClust = esdTrack->GetTPCNclsIter1();
+ } else {
+ nClust = esdTrack->GetTPCclusters(0);
+ }
+
+ Float_t chi2PerCluster = 0.;
+ if(GetAnalysisMode() == AlidNdPtHelper::kTPC) {
+ if(nClust>0.) chi2PerCluster = esdTrack->GetTPCchi2Iter1()/Float_t(nClust);
+ } else {
+ chi2PerCluster = esdTrack->GetTPCchi2()/Float_t(nClust);
+ }
+
+ Int_t nClustersITS = esdTrack->GetITSclusters(0);
+ Float_t chi2PerClusterITS = -1;
+ if (nClustersITS!=0) {
+ chi2PerClusterITS = esdTrack->GetITSchi2()/Float_t(nClustersITS);
+ }
+
+ Float_t clustPerFindClust = 0.;
+ Int_t nFindableClust = esdTrack->GetTPCNclsF();
+ if(nFindableClust>0.) clustPerFindClust = Float_t(nClust)/nFindableClust;
+
+ Float_t b[2], bCov[3];
+ esdTrack->GetImpactParameters(b,bCov);
+
+ //
+ Float_t nCrossedRowsTPC = esdTrack->GetTPCClusterInfo(2,1);
+
+ Float_t ratioCrossedRowsOverFindableClustersTPC = 1.0;
+ if (esdTrack->GetTPCNclsF()>0) {
+ ratioCrossedRowsOverFindableClustersTPC = esdTrack->GetTPCClusterInfo(2,1)/esdTrack->GetTPCNclsF();
+ }
+
+ //
+ Int_t nClustersTPCShared = esdTrack->GetTPCnclsS();
+ Float_t fracClustersTPCShared = -1.;
+ fracClustersTPCShared = Float_t(nClustersTPCShared)/Float_t(nClust);
+
+ //
+ // kink idx
+ Int_t kinkIdx = 0;
+ //if(esdTrack->GetKinkIndex(0) > 0.) isKink = kTRUE;
+ if(esdTrack->GetKinkIndex(0) > 0) kinkIdx = 1; // kink daughter
+ else if(esdTrack->GetKinkIndex(0) < 0) kinkIdx = -1; // kink mother
+ else kinkIdx = 0; // not kink
+
+ //printf("esdTrack->GetKinkIndex(0) %d \n", esdTrack->GetKinkIndex(0));
+ //printf("esdTrack->GetKinkIndex(1) %d \n", esdTrack->GetKinkIndex(1));
+ //printf("esdTrack->GetKinkIndex(2) %d \n", esdTrack->GetKinkIndex(2));
+ //printf("kinkIdx %d \n", kinkIdx);
+
+ Bool_t isOK = kTRUE;
+ AliExternalTrackParam * tpcInner = (AliExternalTrackParam *)(esdTrack->GetTPCInnerParam());
+ if (!tpcInner) return;
+ tpcInner->Rotate(esdTrack->GetAlpha());
+ tpcInner->PropagateTo(esdTrack->GetX(),esdEvent->GetMagneticField());
+
+ // tpc constrained
+ AliExternalTrackParam * tpcInnerC = (AliExternalTrackParam *)(esdTrack->GetTPCInnerParam()->Clone());
+ if (!tpcInnerC) return;
+ tpcInnerC->Rotate(esdTrack->GetAlpha());
+ tpcInnerC->PropagateTo(esdTrack->GetX(),esdEvent->GetMagneticField());
+ Double_t dz[2],cov[3];
+ AliESDVertex *vtx= (AliESDVertex *)esdEvent->GetPrimaryVertex();
+
+ if (!tpcInnerC->PropagateToDCA(vtx, esdEvent->GetMagneticField(), 3, dz, cov))
+ {
+ if(tpcInnerC) delete tpcInnerC;
+ return;
+ }
+ Double_t covar[6]; vtx->GetCovMatrix(covar);
+ Double_t p[2]={tpcInnerC->GetParameter()[0]-dz[0],tpcInnerC->GetParameter()[1]-dz[1]};
+ Double_t c[3]={covar[2],0.,covar[5]};
+ Double_t chi2C=tpcInnerC->GetPredictedChi2(p,c);
+ isOK = tpcInnerC->Update(p,c);
+ isOK = tpcInnerC->Rotate(esdTrack->GetAlpha());
+ isOK = tpcInnerC->Propagate(esdTrack->GetAlpha(),esdTrack->GetX(),esdEvent->GetMagneticField());
+
+ if(!isOK || chi2C>kVeryBig) {
+ if(tpcInnerC) delete tpcInnerC;
+ return;
+ }
+
+ //
+ // calculate chi2 between vi and vj vectors
+ // with covi and covj covariance matrices
+ // chi2ij = (vi-vj)^(T)*(covi+covj)^(-1)*(vi-vj)
+ //
+ TMatrixD deltaT(5,1);
+ TMatrixD delta(1,5);
+ TMatrixD covarM(5,5);
+
+ for (Int_t ipar=0; ipar<5; ipar++) deltaT(ipar,0)=tpcInnerC->GetParameter()[ipar]-esdTrack->GetParameter()[ipar];
+ for (Int_t ipar=0; ipar<5; ipar++) delta(0,ipar)=tpcInnerC->GetParameter()[ipar]-esdTrack->GetParameter()[ipar];
+
+ for (Int_t ipar=0; ipar<5; ipar++) {
+ for (Int_t jpar=0; jpar<5; jpar++) {
+ Int_t index=esdTrack->GetIndex(ipar,jpar);
+ covarM(ipar,jpar)=esdTrack->GetCovariance()[index]+tpcInnerC->GetCovariance()[index];
+ }
+ }
+ TMatrixD covarMInv = covarM.Invert();
+ TMatrixD mat2 = covarMInv*deltaT;
+ TMatrixD chi2 = delta*mat2;
+
+ //
+ // Fill rec vs MC information
+ //
+ Bool_t isPrim = kTRUE;
+ Bool_t hasStrangeMother = kFALSE;
+ Bool_t isFromMaterial = kFALSE;
+ Bool_t isFromConversion = kFALSE;
+
+ if(IsUseMCInfo())
+ {
+ if(!stack) return;
+ Int_t label = TMath::Abs(esdTrack->GetLabel());
+ TParticle* particle = stack->Particle(label);
+ if(!particle) return;
+ if(particle->GetPDG() && particle->GetPDG()->Charge()==0.) return;
+ isPrim = stack->IsPhysicalPrimary(label);
+
+ // check whether has stange mother
+ //
+ Int_t motherPdg = -1;
+ TParticle* mother = 0;
+
+ Int_t motherLabel = particle->GetMother(0);
+ if(motherLabel>0) mother = stack->Particle(motherLabel);
+ if(mother) motherPdg = TMath::Abs(mother->GetPdgCode()); // take abs for visualisation only
+ Int_t mech = particle->GetUniqueID(); // production mechanism
+
+ // K+-, lambda, antilambda, K0s decays
+ if(!isPrim && mech==4 &&
+ (TMath::Abs(motherPdg)==kKPlus || TMath::Abs(motherPdg)==kLambda0 || motherPdg==kK0Short))
+ {
+ hasStrangeMother = kTRUE;
+ }
+
+ if(!isPrim && mech==5 && motherPdg==kGamma) {
+ isFromConversion=kTRUE;
+ }
+
+ if(!isPrim && mech==13) {
+ isFromMaterial=kTRUE;
+ }
+ }
+
+ // fill histo
+ Int_t charge = esdTrack->Charge();
+ Double_t vRecMCTrackHist[17] = { nCrossedRowsTPC, chi2PerCluster, chi2PerClusterITS,ratioCrossedRowsOverFindableClustersTPC, fracClustersTPCShared, b[0], b[1], eta, phi, pt, hasStrangeMother, isFromConversion, isFromMaterial, isPrim, charge, centralityF, chi2(0,0) };
+ fRecMCTrackHist->Fill(vRecMCTrackHist);
+
+ if(tpcInnerC) delete tpcInnerC;
+}
+
+//_____________________________________________________________________________
+Long64_t AlidNdPtCutAnalysisPbPb::Merge(TCollection* const list)
+{
+ // Merge list of objects (needed by PROOF)
+
+ if (!list)
+ return 0;
+
+ if (list->IsEmpty())
+ return 1;
+
+ TIterator* iter = list->MakeIterator();
+ TObject* obj = 0;
+ TObjArray* objArrayList = 0;
+ objArrayList = new TObjArray();
+
+ //TList *collPhysSelection = new TList;
+
+ // collection of generated histograms
+ Int_t count=0;
+ while((obj = iter->Next()) != 0) {
+ AlidNdPtCutAnalysisPbPb* entry = dynamic_cast<AlidNdPtCutAnalysisPbPb*>(obj);
+ if (entry == 0) continue;
+
+ // event histo
+ if(GetMergeTHnSparse())
+ {
+ fEventCount->Add(entry->fEventCount);
+ fRecEventHist->Add(entry->fRecEventHist);
+ fRecMCEventHist->Add(entry->fRecMCEventHist);
+ fMCEventHist->Add(entry->fMCEventHist);
+
+ // track histo
+ fRecMCTrackHist->Add(entry->fRecMCTrackHist);
+ }
+
+ if (entry->fFolderObj) { objArrayList->Add(entry->fFolderObj); }
+
+ count++;
+ }
+ if (fFolderObj) { fFolderObj->Merge(objArrayList); }
+
+ // to signal that track histos were not merged: reset
+ if (!GetMergeTHnSparse())
+ {
+ // reset content
+ // save memory
+ fEventCount->Reset();
+ fRecEventHist->Reset();
+ fRecMCEventHist->Reset();
+ fMCEventHist->Reset();
+
+ // track histo
+ fRecMCTrackHist->Reset();
+ }
+
+ // delete
+ if (objArrayList) delete objArrayList; objArrayList=0;
+
+return count;
+}
+
+//_____________________________________________________________________________
+void AlidNdPtCutAnalysisPbPb::Analyse()
+{
+ //
+ // Analyse histograms
+ //
+ TH1::AddDirectory(kFALSE);
+ TObjArray *aFolderObj = new TObjArray;
+ if(!aFolderObj) return;
+
+ TH1D *h1D = 0;
+ TH2D *h2D = 0;
+ TH3D *h3D = 0;
+
+ //
+ // make event level projection
+ //
+ //THnSparseF *fEventCount; //-> trig, trig + vertex
+ //THnSparseF *fRecEventHist; //-> Xv:Yv:Zv:ResZv:Mult
+ //THnSparseF *fMCEventHist; //-> mcXv:mcYv:mcZv
+ //THnSparseF *fRecMCEventHist; //-> Xv-mcXv:Yv-mcYv:Zv-mcZv:Mult
+
+ //
+ // Event counters
+ //
+ h2D = (TH2D*)fEventCount->Projection(0,1);
+ h2D->SetName("trig_vs_trigANDvertex");
+ h2D->SetTitle("trig_vs_trigANDvertex");
+ aFolderObj->Add(h2D);
+
+ fEventCount->GetAxis(0)->SetRange(2,2); // triggered
+ h1D = (TH1D*)fEventCount->Projection(1);
+ h1D->SetTitle("rec. vertex for triggered events");
+ h1D->SetName("trigANDvertex");
+ h1D->SetTitle("trigANDvertex");
+ aFolderObj->Add(h1D);
+
+ //
+ // Create rec. event histograms
+ //
+ h1D = (TH1D *)fRecEventHist->Projection(0);
+ h1D->SetName("rec_xv");
+ h1D->SetTitle("rec_xv");
+ aFolderObj->Add(h1D);
+
+ h1D = (TH1D *)fRecEventHist->Projection(1);
+ h1D->SetName("rec_yv");
+ h1D->SetTitle("rec_yv");
+ aFolderObj->Add(h1D);
+
+ h1D = (TH1D *)fRecEventHist->Projection(2);
+ h1D->SetName("rec_zv");
+ h1D->SetTitle("rec_zv");
+ aFolderObj->Add(h1D);
+
+ h2D = (TH2D *)fRecEventHist->Projection(0,1);
+ h2D->SetName("rec_xv_vs_yv");
+ h2D->SetTitle("rec_xv_vs_yv");
+ aFolderObj->Add(h2D);
+
+ h2D = (TH2D *)fRecEventHist->Projection(1,2);
+ h2D->SetName("rec_yv_vs_zv");
+ h2D->SetTitle("rec_yv_vs_zv");
+ aFolderObj->Add(h2D);
+
+ h2D = (TH2D *)fRecEventHist->Projection(0,2);
+ h2D->SetName("rec_xv_vs_zv");
+ h2D->SetTitle("rec_xv_vs_zv");
+ aFolderObj->Add(h2D);
+
+ h2D = (TH2D *)fRecEventHist->Projection(3,4);
+ h2D->SetName("rec_resZv_vs_Mult");
+ h2D->SetTitle("rec_resZv_vs_Mult");
+ aFolderObj->Add(h2D);
+
+ //
+ // MC available
+ //
+ if(IsUseMCInfo()) {
+
+ //
+ // Create mc event histograms
+ //
+
+ h1D = (TH1D *)fMCEventHist->Projection(0);
+ h1D->SetName("mc_xv");
+ h1D->SetTitle("mc_xv");
+ aFolderObj->Add(h1D);
+
+ h1D = (TH1D *)fMCEventHist->Projection(1);
+ h1D->SetName("mc_yv");
+ h1D->SetTitle("mc_yv");
+ aFolderObj->Add(h1D);
+
+ h1D = (TH1D *)fMCEventHist->Projection(2);
+ h1D->SetName("mc_zv");
+ h1D->SetTitle("mc_zv");
+ aFolderObj->Add(h1D);
+
+
+ h2D = (TH2D *)fMCEventHist->Projection(0,1);
+ h2D->SetName("mc_xv_vs_yv");
+ h2D->SetTitle("mc_xv_vs_yv");
+ aFolderObj->Add(h2D);
+
+ h2D = (TH2D *)fMCEventHist->Projection(1,2);
+ h2D->SetName("mc_yv_vs_zv");
+ h2D->SetTitle("mc_yv_vs_zv");
+ aFolderObj->Add(h2D);
+
+ h2D = (TH2D *)fMCEventHist->Projection(0,2);
+ h2D->SetName("mc_xv_vs_zv");
+ h2D->SetTitle("mc_xv_vs_zv");
+ aFolderObj->Add(h2D);
+
+ //
+ // Create rec-mc event histograms
+ //
+ h2D = (TH2D *)fRecMCEventHist->Projection(0,3);
+ h2D->SetName("rec_mc_deltaXv_vs_mult");
+ h2D->SetTitle("rec_mc_deltaXv_vs_mult");
+ aFolderObj->Add(h2D);
+
+ h2D = (TH2D *)fRecMCEventHist->Projection(1,3);
+ h2D->SetName("rec_mc_deltaYv_vs_mult");
+ h2D->SetTitle("rec_mc_deltaYv_vs_mult");
+ aFolderObj->Add(h2D);
+
+ h2D = (TH2D *)fRecMCEventHist->Projection(2,3);
+ h2D->SetName("rec_mc_deltaZv_vs_mult");
+ h2D->SetTitle("rec_mc_deltaZv_vs_mult");
+ aFolderObj->Add(h2D);
+
+ } // end use MC info
+
+ //
+ // make track level projection
+ //
+ // THnSparse track histograms
+ // nCrossRows:chi2PerClust:chi2PerClustITS:nCrossRows/nFindableClust:fracSharedClust:DCAy:DCAz:eta:phi:pt:hasStrangeMother:isFromConversion:isFromMaterial:isPrim:charge:centr:chi2ToTPCc
+ //
+
+ // only TH3 histograms
+ // x : pT : centr
+ Double_t minNCrossRows = 120.;
+ Double_t maxChi2PerClust = 4.;
+ Double_t maxChi2PerClustITS = 36.;
+ Double_t minNCrossRowsOverFindable = 0.8;
+ Double_t maxFracSharedClust = 0.4;
+ Double_t minDCAr = -0.2, maxDCAr = 0.2;
+ Double_t minDCAz = -2., maxDCAz = 2.;
+ Double_t minEta = -0.8, maxEta = 0.8;
+ Double_t maxChi2ToTPCc = 50;
+
+ // only TH3 histograms
+ // x : pT : centr
+
+ // set all cuts
+ fRecMCTrackHist->GetAxis(0)->SetRangeUser(minNCrossRows,160);
+ fRecMCTrackHist->GetAxis(1)->SetRangeUser(0.,maxChi2PerClust-0.0001);
+ fRecMCTrackHist->GetAxis(2)->SetRangeUser(0.,maxChi2PerClustITS-0.0001);
+ fRecMCTrackHist->GetAxis(3)->SetRangeUser(minNCrossRowsOverFindable,1.4999);
+ fRecMCTrackHist->GetAxis(4)->SetRangeUser(0.,maxFracSharedClust-0.0001);
+ fRecMCTrackHist->GetAxis(5)->SetRangeUser(minDCAr,maxDCAr-0.0001);
+ fRecMCTrackHist->GetAxis(6)->SetRangeUser(minDCAz,maxDCAz-0.0001);
+ fRecMCTrackHist->GetAxis(7)->SetRangeUser(minEta,maxEta-0.0001);
+ fRecMCTrackHist->GetAxis(16)->SetRangeUser(0.,maxChi2ToTPCc-0.0001);
+
+ // nCrossRows
+ fRecMCTrackHist->GetAxis(0)->SetRangeUser(0,160);
+
+ h3D = (TH3D *)fRecMCTrackHist->Projection(0,9,15);
+ h3D->SetName("rec_nCrossRows_pT_centr");
+ h3D->SetTitle("rec_nCrossRows_pT_centr");
+ aFolderObj->Add(h3D);
+
+
+ // Chi2PerClust
+ fRecMCTrackHist->GetAxis(0)->SetRangeUser(minNCrossRows,160);
+ fRecMCTrackHist->GetAxis(1)->SetRangeUser(0.,10);
+
+ h3D = (TH3D *)fRecMCTrackHist->Projection(1,9,15);
+ h3D->SetName("rec_Chi2PerClust_pT_centr");
+ h3D->SetTitle("rec_Chi2PerClust_pT_centr");
+ aFolderObj->Add(h3D);
+
+
+ // Chi2PerClustITS
+ fRecMCTrackHist->GetAxis(1)->SetRangeUser(0.,maxChi2PerClust-0.0001);
+ fRecMCTrackHist->GetAxis(2)->SetRangeUser(0.,100);
+
+ h3D = (TH3D *)fRecMCTrackHist->Projection(2,9,15);
+ h3D->SetName("rec_Chi2PerClustITS_pT_centr");
+ h3D->SetTitle("rec_Chi2PerClustITS_pT_centr");
+ aFolderObj->Add(h3D);
+
+
+ // NCrossRowsOverFindable
+ fRecMCTrackHist->GetAxis(2)->SetRangeUser(0.,maxChi2PerClustITS-0.0001);
+ fRecMCTrackHist->GetAxis(3)->SetRangeUser(0,1.4999);
+
+ h3D = (TH3D *)fRecMCTrackHist->Projection(3,9,15);
+ h3D->SetName("rec_NCrossRowsOverFindable_pT_centr");
+ h3D->SetTitle("rec_NCrossRowsOverFindable_pT_centr");
+ aFolderObj->Add(h3D);
+
+
+ // FracSharedClust
+ fRecMCTrackHist->GetAxis(3)->SetRangeUser(minNCrossRowsOverFindable,1.4999);
+ fRecMCTrackHist->GetAxis(4)->SetRangeUser(0.,0.9999);
+
+ h3D = (TH3D *)fRecMCTrackHist->Projection(4,9,15);
+ h3D->SetName("rec_FracSharedClust_pT_centr");
+ h3D->SetTitle("rec_FracSharedClust_pT_centr");
+ aFolderObj->Add(h3D);
+
+
+ // DCAr
+ fRecMCTrackHist->GetAxis(4)->SetRangeUser(0.,maxFracSharedClust-0.0001);
+ fRecMCTrackHist->GetAxis(5)->SetRangeUser(-10.,9.9999);
+
+ h3D = (TH3D *)fRecMCTrackHist->Projection(5,9,15);
+ h3D->SetName("rec_DCAr_pT_centr");
+ h3D->SetTitle("rec_DCAr_pT_centr");
+ aFolderObj->Add(h3D);
+
+ // DCAz
+ fRecMCTrackHist->GetAxis(5)->SetRangeUser(minDCAr,maxDCAr-0.0001);
+ fRecMCTrackHist->GetAxis(6)->SetRangeUser(-10.,9.9999);
+
+ h3D = (TH3D *)fRecMCTrackHist->Projection(6,9,15);
+ h3D->SetName("rec_DCAz_pT_centr");
+ h3D->SetTitle("rec_DCAz_pT_centr");
+ aFolderObj->Add(h3D);
+
+ // chi2ToTPCc
+ fRecMCTrackHist->GetAxis(6)->SetRangeUser(minDCAz,maxDCAz-0.0001);
+ fRecMCTrackHist->GetAxis(16)->SetRangeUser(0,100);
+
+ h3D = (TH3D *)fRecMCTrackHist->Projection(16,9,15);
+ h3D->SetName("rec_chi2ToTPCc_pT_centr");
+ h3D->SetTitle("rec_chi2ToTPCc_pT_centr");
+ aFolderObj->Add(h3D);
+
+
+ //
+ // MC available
+ //
+ // nCrossRows:chi2PerClust:chi2PerClustITS:nCrossRows/nFindableClust:fracSharedClust:DCAy:DCAz:eta:phi:pt:hasStrangeMother:isFromConversion:isFromMaterial:isPrim:charge:centr:chi2ToTPCc
+ //
+
+ if(IsUseMCInfo()) {
+ //
+ // only TH3 histograms
+ // x : pT : centr
+ //
+ // comes from week decays
+ //
+ fRecMCTrackHist->GetAxis(10)->SetRange(2,2);
+
+ // set all cuts
+ fRecMCTrackHist->GetAxis(0)->SetRangeUser(minNCrossRows,160);
+ fRecMCTrackHist->GetAxis(1)->SetRangeUser(0.,maxChi2PerClust-0.0001);
+ fRecMCTrackHist->GetAxis(2)->SetRangeUser(0.,maxChi2PerClustITS-0.0001);
+ fRecMCTrackHist->GetAxis(3)->SetRangeUser(minNCrossRowsOverFindable,1.4999);
+ fRecMCTrackHist->GetAxis(4)->SetRangeUser(0.,maxFracSharedClust-0.0001);
+ fRecMCTrackHist->GetAxis(5)->SetRangeUser(minDCAr,maxDCAr-0.0001);
+ fRecMCTrackHist->GetAxis(6)->SetRangeUser(minDCAz,maxDCAz-0.0001);
+ fRecMCTrackHist->GetAxis(7)->SetRangeUser(minEta,maxEta-0.0001);
+ fRecMCTrackHist->GetAxis(16)->SetRangeUser(0.,maxChi2ToTPCc-0.0001);
+
+ // nCrossRows
+ fRecMCTrackHist->GetAxis(0)->SetRangeUser(0,160);
+
+ h3D = (TH3D *)fRecMCTrackHist->Projection(0,9,15);
+ h3D->SetName("rec_sec_strange_nCrossRows_pT_centr");
+ h3D->SetTitle("rec_sec_strange_nCrossRows_pT_centr");
+ aFolderObj->Add(h3D);
+
+ // Chi2PerClust
+ fRecMCTrackHist->GetAxis(0)->SetRangeUser(minNCrossRows,160);
+ fRecMCTrackHist->GetAxis(1)->SetRangeUser(0.,10);
+
+ h3D = (TH3D *)fRecMCTrackHist->Projection(1,9,15);
+ h3D->SetName("rec_sec_strange_Chi2PerClust_pT_centr");
+ h3D->SetTitle("rec_sec_strange_Chi2PerClust_pT_centr");
+ aFolderObj->Add(h3D);
+
+ // Chi2PerClustITS
+ fRecMCTrackHist->GetAxis(1)->SetRangeUser(0.,maxChi2PerClust-0.0001);
+ fRecMCTrackHist->GetAxis(2)->SetRangeUser(0.,100);
+
+ h3D = (TH3D *)fRecMCTrackHist->Projection(2,9,15);
+ h3D->SetName("rec_sec_strange_Chi2PerClustITS_pT_centr");
+ h3D->SetTitle("rec_sec_strange_Chi2PerClustITS_pT_centr");
+ aFolderObj->Add(h3D);
+
+
+ // NCrossRowsOverFindable
+ fRecMCTrackHist->GetAxis(2)->SetRangeUser(0.,maxChi2PerClustITS-0.0001);
+ fRecMCTrackHist->GetAxis(3)->SetRangeUser(0,1.4999);
+
+ h3D = (TH3D *)fRecMCTrackHist->Projection(3,9,15);
+ h3D->SetName("rec_sec_strange_NCrossRowsOverFindable_pT_centr");
+ h3D->SetTitle("rec_sec_strange_NCrossRowsOverFindable_pT_centr");
+ aFolderObj->Add(h3D);
+
+
+ // FracSharedClust
+ fRecMCTrackHist->GetAxis(3)->SetRangeUser(minNCrossRowsOverFindable,1.4999);
+ fRecMCTrackHist->GetAxis(4)->SetRangeUser(0.,0.9999);
+
+ h3D = (TH3D *)fRecMCTrackHist->Projection(4,9,15);
+ h3D->SetName("rec_sec_strange_FracSharedClust_pT_centr");
+ h3D->SetTitle("rec_sec_strange_FracSharedClust_pT_centr");
+ aFolderObj->Add(h3D);
+
+
+ // DCAr
+ fRecMCTrackHist->GetAxis(4)->SetRangeUser(0.,maxFracSharedClust-0.0001);
+ fRecMCTrackHist->GetAxis(5)->SetRangeUser(-10.,9.9999);
+
+ h3D = (TH3D *)fRecMCTrackHist->Projection(5,9,15);
+ h3D->SetName("rec_sec_strange_DCAr_pT_centr");
+ h3D->SetTitle("rec_sec_strange_DCAr_pT_centr");
+ aFolderObj->Add(h3D);
+
+
+ // DCAz
+ fRecMCTrackHist->GetAxis(5)->SetRangeUser(minDCAr,maxDCAr-0.0001);
+ fRecMCTrackHist->GetAxis(6)->SetRangeUser(-10.,9.9999);
+
+ h3D = (TH3D *)fRecMCTrackHist->Projection(6,9,15);
+ h3D->SetName("rec_sec_strange_DCAz_pT_centr");
+ h3D->SetTitle("rec_sec_strange_DCAz_pT_centr");
+ aFolderObj->Add(h3D);
+
+ // chi2ToTPCc
+ fRecMCTrackHist->GetAxis(6)->SetRangeUser(minDCAz,maxDCAz-0.0001);
+ fRecMCTrackHist->GetAxis(16)->SetRangeUser(0,100);
+
+ h3D = (TH3D *)fRecMCTrackHist->Projection(16,9,15);
+ h3D->SetName("rec_sec_strange_chi2ToTPCc_pT_centr");
+ h3D->SetTitle("rec_sec_strange_chi2ToTPCc_pT_centr");
+ aFolderObj->Add(h3D);
+
+
+ // only TH3 histograms
+ // x : pT : centr
+ //
+
+ // comes from conversion
+ //
+ fRecMCTrackHist->GetAxis(10)->SetRange(1,2);
+ fRecMCTrackHist->GetAxis(11)->SetRange(2,2);
+
+ // set all cuts
+ fRecMCTrackHist->GetAxis(0)->SetRangeUser(minNCrossRows,160);
+ fRecMCTrackHist->GetAxis(1)->SetRangeUser(0.,maxChi2PerClust-0.0001);
+ fRecMCTrackHist->GetAxis(2)->SetRangeUser(0.,maxChi2PerClustITS-0.0001);
+ fRecMCTrackHist->GetAxis(3)->SetRangeUser(minNCrossRowsOverFindable,1.4999);
+ fRecMCTrackHist->GetAxis(4)->SetRangeUser(0.,maxFracSharedClust-0.0001);
+ fRecMCTrackHist->GetAxis(5)->SetRangeUser(minDCAr,maxDCAr-0.0001);
+ fRecMCTrackHist->GetAxis(6)->SetRangeUser(minDCAz,maxDCAz-0.0001);
+ fRecMCTrackHist->GetAxis(7)->SetRangeUser(minEta,maxEta-0.0001);
+ fRecMCTrackHist->GetAxis(16)->SetRangeUser(0.,maxChi2ToTPCc-0.0001);
+
+ // nCrossRows
+ fRecMCTrackHist->GetAxis(0)->SetRangeUser(0,160);
+
+ h3D = (TH3D *)fRecMCTrackHist->Projection(0,9,15);
+ h3D->SetName("rec_sec_conversion_nCrossRows_pT_centr");
+ h3D->SetTitle("rec_sec_conversion_nCrossRows_pT_centr");
+ aFolderObj->Add(h3D);
+
+ // Chi2PerClust
+ fRecMCTrackHist->GetAxis(0)->SetRangeUser(minNCrossRows,160);
+ fRecMCTrackHist->GetAxis(1)->SetRangeUser(0.,10);
+
+ h3D = (TH3D *)fRecMCTrackHist->Projection(1,9,15);
+ h3D->SetName("rec_sec_conversion_Chi2PerClust_pT_centr");
+ h3D->SetTitle("rec_sec_conversion_Chi2PerClust_pT_centr");
+ aFolderObj->Add(h3D);
+
+ // Chi2PerClustITS
+ fRecMCTrackHist->GetAxis(1)->SetRangeUser(0.,maxChi2PerClust-0.0001);
+ fRecMCTrackHist->GetAxis(2)->SetRangeUser(0.,100);
+
+ h3D = (TH3D *)fRecMCTrackHist->Projection(2,9,15);
+ h3D->SetName("rec_sec_conversion_Chi2PerClustITS_pT_centr");
+ h3D->SetTitle("rec_sec_conversion_Chi2PerClustITS_pT_centr");
+ aFolderObj->Add(h3D);
+
+
+ // NCrossRowsOverFindable
+ fRecMCTrackHist->GetAxis(2)->SetRangeUser(0.,maxChi2PerClustITS-0.0001);
+ fRecMCTrackHist->GetAxis(3)->SetRangeUser(0,1.4999);
+
+ h3D = (TH3D *)fRecMCTrackHist->Projection(3,9,15);
+ h3D->SetName("rec_sec_conversion_NCrossRowsOverFindable_pT_centr");
+ h3D->SetTitle("rec_sec_conversion_NCrossRowsOverFindable_pT_centr");
+ aFolderObj->Add(h3D);
+
+
+ // FracSharedClust
+ fRecMCTrackHist->GetAxis(3)->SetRangeUser(minNCrossRowsOverFindable,1.4999);
+ fRecMCTrackHist->GetAxis(4)->SetRangeUser(0.,0.9999);
+
+ h3D = (TH3D *)fRecMCTrackHist->Projection(4,9,15);
+ h3D->SetName("rec_sec_conversion_FracSharedClust_pT_centr");
+ h3D->SetTitle("rec_sec_conversion_FracSharedClust_pT_centr");
+ aFolderObj->Add(h3D);
+
+
+ // DCAr
+ fRecMCTrackHist->GetAxis(4)->SetRangeUser(0.,maxFracSharedClust-0.0001);
+ fRecMCTrackHist->GetAxis(5)->SetRangeUser(-10.,9.9999);
+
+ h3D = (TH3D *)fRecMCTrackHist->Projection(5,9,15);
+ h3D->SetName("rec_sec_conversion_DCAr_pT_centr");
+ h3D->SetTitle("rec_sec_conversion_DCAr_pT_centr");
+ aFolderObj->Add(h3D);
+
+
+ // DCAz
+ fRecMCTrackHist->GetAxis(5)->SetRangeUser(minDCAr,maxDCAr-0.0001);
+ fRecMCTrackHist->GetAxis(6)->SetRangeUser(-10.,9.9999);
+
+ h3D = (TH3D *)fRecMCTrackHist->Projection(6,9,15);
+ h3D->SetName("rec_sec_conversion_DCAz_pT_centr");
+ h3D->SetTitle("rec_sec_conversion_DCAz_pT_centr");
+ aFolderObj->Add(h3D);
+
+ // chi2ToTPCc
+ fRecMCTrackHist->GetAxis(6)->SetRangeUser(minDCAz,maxDCAz-0.0001);
+ fRecMCTrackHist->GetAxis(16)->SetRangeUser(0,100);
+
+ h3D = (TH3D *)fRecMCTrackHist->Projection(16,9,15);
+ h3D->SetName("rec_sec_conversion_chi2ToTPCc_pT_centr");
+ h3D->SetTitle("rec_sec_conversion_chi2ToTPCc_pT_centr");
+ aFolderObj->Add(h3D);
+
+
+
+
+
+ // comes from interaction with material
+ //
+ fRecMCTrackHist->GetAxis(11)->SetRange(1,2);
+ fRecMCTrackHist->GetAxis(12)->SetRange(2,2);
+
+ // set all cuts
+ fRecMCTrackHist->GetAxis(0)->SetRangeUser(minNCrossRows,160);
+ fRecMCTrackHist->GetAxis(1)->SetRangeUser(0.,maxChi2PerClust-0.0001);
+ fRecMCTrackHist->GetAxis(2)->SetRangeUser(0.,maxChi2PerClustITS-0.0001);
+ fRecMCTrackHist->GetAxis(3)->SetRangeUser(minNCrossRowsOverFindable,1.4999);
+ fRecMCTrackHist->GetAxis(4)->SetRangeUser(0.,maxFracSharedClust-0.0001);
+ fRecMCTrackHist->GetAxis(5)->SetRangeUser(minDCAr,maxDCAr-0.0001);
+ fRecMCTrackHist->GetAxis(6)->SetRangeUser(minDCAz,maxDCAz-0.0001);
+ fRecMCTrackHist->GetAxis(7)->SetRangeUser(minEta,maxEta-0.0001);
+ fRecMCTrackHist->GetAxis(16)->SetRangeUser(0.,maxChi2ToTPCc-0.0001);
+
+ // nCrossRows
+ fRecMCTrackHist->GetAxis(0)->SetRangeUser(0,160);
+
+ h3D = (TH3D *)fRecMCTrackHist->Projection(0,9,15);
+ h3D->SetName("rec_sec_material_nCrossRows_pT_centr");
+ h3D->SetTitle("rec_sec_material_nCrossRows_pT_centr");
+ aFolderObj->Add(h3D);
+
+ // Chi2PerClust
+ fRecMCTrackHist->GetAxis(0)->SetRangeUser(minNCrossRows,160);
+ fRecMCTrackHist->GetAxis(1)->SetRangeUser(0.,10);
+
+ h3D = (TH3D *)fRecMCTrackHist->Projection(1,9,15);
+ h3D->SetName("rec_sec_material_Chi2PerClust_pT_centr");
+ h3D->SetTitle("rec_sec_material_Chi2PerClust_pT_centr");
+ aFolderObj->Add(h3D);
+
+ // Chi2PerClustITS
+ fRecMCTrackHist->GetAxis(1)->SetRangeUser(0.,maxChi2PerClust-0.0001);
+ fRecMCTrackHist->GetAxis(2)->SetRangeUser(0.,100);
+
+ h3D = (TH3D *)fRecMCTrackHist->Projection(2,9,15);
+ h3D->SetName("rec_sec_material_Chi2PerClustITS_pT_centr");
+ h3D->SetTitle("rec_sec_material_Chi2PerClustITS_pT_centr");
+ aFolderObj->Add(h3D);
+
+
+ // NCrossRowsOverFindable
+ fRecMCTrackHist->GetAxis(2)->SetRangeUser(0.,maxChi2PerClustITS-0.0001);
+ fRecMCTrackHist->GetAxis(3)->SetRangeUser(0,1.4999);
+
+ h3D = (TH3D *)fRecMCTrackHist->Projection(3,9,15);
+ h3D->SetName("rec_sec_material_NCrossRowsOverFindable_pT_centr");
+ h3D->SetTitle("rec_sec_material_NCrossRowsOverFindable_pT_centr");
+ aFolderObj->Add(h3D);
+
+
+ // FracSharedClust
+ fRecMCTrackHist->GetAxis(3)->SetRangeUser(minNCrossRowsOverFindable,1.4999);
+ fRecMCTrackHist->GetAxis(4)->SetRangeUser(0.,0.9999);
+
+ h3D = (TH3D *)fRecMCTrackHist->Projection(4,9,15);
+ h3D->SetName("rec_sec_material_FracSharedClust_pT_centr");
+ h3D->SetTitle("rec_sec_material_FracSharedClust_pT_centr");
+ aFolderObj->Add(h3D);
+
+
+ // DCAr
+ fRecMCTrackHist->GetAxis(4)->SetRangeUser(0.,maxFracSharedClust-0.0001);
+ fRecMCTrackHist->GetAxis(5)->SetRangeUser(-10.,9.9999);
+
+ h3D = (TH3D *)fRecMCTrackHist->Projection(5,9,15);
+ h3D->SetName("rec_sec_material_DCAr_pT_centr");
+ h3D->SetTitle("rec_sec_material_DCAr_pT_centr");
+ aFolderObj->Add(h3D);
+
+
+ // DCAz
+ fRecMCTrackHist->GetAxis(5)->SetRangeUser(minDCAr,maxDCAr-0.0001);
+ fRecMCTrackHist->GetAxis(6)->SetRangeUser(-10.,9.9999);
+
+ h3D = (TH3D *)fRecMCTrackHist->Projection(6,9,15);
+ h3D->SetName("rec_sec_material_DCAz_pT_centr");
+ h3D->SetTitle("rec_sec_material_DCAz_pT_centr");
+ aFolderObj->Add(h3D);
+
+ // chi2ToTPCc
+ fRecMCTrackHist->GetAxis(6)->SetRangeUser(minDCAz,maxDCAz-0.0001);
+ fRecMCTrackHist->GetAxis(16)->SetRangeUser(0,100);
+
+ h3D = (TH3D *)fRecMCTrackHist->Projection(16,9,15);
+ h3D->SetName("rec_sec_material_chi2ToTPCc_pT_centr");
+ h3D->SetTitle("rec_sec_material_chi2ToTPCc_pT_centr");
+ aFolderObj->Add(h3D);
+
+
+
+ //
+ // only TH3 histograms
+ // x : pT : centr
+ //
+ // secondaries defined as AliStack::IsPhysicalPrimary()
+ //
+ fRecMCTrackHist->GetAxis(12)->SetRange(1,2);
+ fRecMCTrackHist->GetAxis(13)->SetRange(1,1);
+
+ // set all cuts
+ fRecMCTrackHist->GetAxis(0)->SetRangeUser(minNCrossRows,160);
+ fRecMCTrackHist->GetAxis(1)->SetRangeUser(0.,maxChi2PerClust-0.0001);
+ fRecMCTrackHist->GetAxis(2)->SetRangeUser(0.,maxChi2PerClustITS-0.0001);
+ fRecMCTrackHist->GetAxis(3)->SetRangeUser(minNCrossRowsOverFindable,1.4999);
+ fRecMCTrackHist->GetAxis(4)->SetRangeUser(0.,maxFracSharedClust-0.0001);
+ fRecMCTrackHist->GetAxis(5)->SetRangeUser(minDCAr,maxDCAr-0.0001);
+ fRecMCTrackHist->GetAxis(6)->SetRangeUser(minDCAz,maxDCAz-0.0001);
+ fRecMCTrackHist->GetAxis(7)->SetRangeUser(minEta,maxEta-0.0001);
+ fRecMCTrackHist->GetAxis(16)->SetRangeUser(0.,maxChi2ToTPCc-0.0001);
+
+ // nCrossRows
+ fRecMCTrackHist->GetAxis(0)->SetRangeUser(0,160);
+
+ h3D = (TH3D *)fRecMCTrackHist->Projection(0,9,15);
+ h3D->SetName("rec_sec_nCrossRows_pT_centr");
+ h3D->SetTitle("rec_sec_nCrossRows_pT_centr");
+ aFolderObj->Add(h3D);
+
+ // Chi2PerClust
+ fRecMCTrackHist->GetAxis(0)->SetRangeUser(minNCrossRows,160);
+ fRecMCTrackHist->GetAxis(1)->SetRangeUser(0.,10);
+
+ h3D = (TH3D *)fRecMCTrackHist->Projection(1,9,15);
+ h3D->SetName("rec_sec_Chi2PerClust_pT_centr");
+ h3D->SetTitle("rec_sec_Chi2PerClust_pT_centr");
+ aFolderObj->Add(h3D);
+
+ // Chi2PerClustITS
+ fRecMCTrackHist->GetAxis(1)->SetRangeUser(0.,maxChi2PerClust-0.0001);
+ fRecMCTrackHist->GetAxis(2)->SetRangeUser(0.,100);
+
+ h3D = (TH3D *)fRecMCTrackHist->Projection(2,9,15);
+ h3D->SetName("rec_sec_Chi2PerClustITS_pT_centr");
+ h3D->SetTitle("rec_sec_Chi2PerClustITS_pT_centr");
+ aFolderObj->Add(h3D);
+
+
+ // NCrossRowsOverFindable
+ fRecMCTrackHist->GetAxis(2)->SetRangeUser(0.,maxChi2PerClustITS-0.0001);
+ fRecMCTrackHist->GetAxis(3)->SetRangeUser(0,1.4999);
+
+ h3D = (TH3D *)fRecMCTrackHist->Projection(3,9,15);
+ h3D->SetName("rec_sec_NCrossRowsOverFindable_pT_centr");
+ h3D->SetTitle("rec_sec_NCrossRowsOverFindable_pT_centr");
+ aFolderObj->Add(h3D);
+
+
+ // FracSharedClust
+ fRecMCTrackHist->GetAxis(3)->SetRangeUser(minNCrossRowsOverFindable,1.4999);
+ fRecMCTrackHist->GetAxis(4)->SetRangeUser(0.,0.9999);
+
+ h3D = (TH3D *)fRecMCTrackHist->Projection(4,9,15);
+ h3D->SetName("rec_sec_FracSharedClust_pT_centr");
+ h3D->SetTitle("rec_sec_FracSharedClust_pT_centr");
+ aFolderObj->Add(h3D);
+
+
+ // DCAr
+ fRecMCTrackHist->GetAxis(4)->SetRangeUser(0.,maxFracSharedClust-0.0001);
+ fRecMCTrackHist->GetAxis(5)->SetRangeUser(-10.,9.9999);
+
+ h3D = (TH3D *)fRecMCTrackHist->Projection(5,9,15);
+ h3D->SetName("rec_sec_DCAr_pT_centr");
+ h3D->SetTitle("rec_sec_DCAr_pT_centr");
+ aFolderObj->Add(h3D);
+
+
+ // DCAz
+ fRecMCTrackHist->GetAxis(5)->SetRangeUser(minDCAr,maxDCAr-0.0001);
+ fRecMCTrackHist->GetAxis(6)->SetRangeUser(-10.,9.9999);
+
+ h3D = (TH3D *)fRecMCTrackHist->Projection(6,9,15);
+ h3D->SetName("rec_sec_DCAz_pT_centr");
+ h3D->SetTitle("rec_sec_DCAz_pT_centr");
+ aFolderObj->Add(h3D);
+
+ // chi2ToTPCc
+ fRecMCTrackHist->GetAxis(6)->SetRangeUser(minDCAz,maxDCAz-0.0001);
+ fRecMCTrackHist->GetAxis(16)->SetRangeUser(0,100);
+
+ h3D = (TH3D *)fRecMCTrackHist->Projection(16,9,15);
+ h3D->SetName("rec_sec_chi2ToTPCc_pT_centr");
+ h3D->SetTitle("rec_sec_chi2ToTPCc_pT_centr");
+ aFolderObj->Add(h3D);
+
+ //
+ // only TH3 histograms
+ // x : pT : centr
+ //
+ // primaries defined as AliStack::IsPhysicalPrimary()
+ //
+ fRecMCTrackHist->GetAxis(12)->SetRange(1,2);
+ fRecMCTrackHist->GetAxis(13)->SetRange(2,2);
+
+ // set all cuts
+ fRecMCTrackHist->GetAxis(0)->SetRangeUser(minNCrossRows,160);
+ fRecMCTrackHist->GetAxis(1)->SetRangeUser(0.,maxChi2PerClust-0.0001);
+ fRecMCTrackHist->GetAxis(2)->SetRangeUser(0.,maxChi2PerClustITS-0.0001);
+ fRecMCTrackHist->GetAxis(3)->SetRangeUser(minNCrossRowsOverFindable,1.4999);
+ fRecMCTrackHist->GetAxis(4)->SetRangeUser(0.,maxFracSharedClust-0.0001);
+ fRecMCTrackHist->GetAxis(5)->SetRangeUser(minDCAr,maxDCAr-0.0001);
+ fRecMCTrackHist->GetAxis(6)->SetRangeUser(minDCAz,maxDCAz-0.0001);
+ fRecMCTrackHist->GetAxis(7)->SetRangeUser(minEta,maxEta-0.0001);
+ fRecMCTrackHist->GetAxis(16)->SetRangeUser(0.,maxChi2ToTPCc-0.0001);
+
+ // nCrossRows
+ fRecMCTrackHist->GetAxis(0)->SetRangeUser(0,160);
+
+ h3D = (TH3D *)fRecMCTrackHist->Projection(0,9,15);
+ h3D->SetName("rec_prim_nCrossRows_pT_centr");
+ h3D->SetTitle("rec_prim_nCrossRows_pT_centr");
+ aFolderObj->Add(h3D);
+
+ // Chi2PerClust
+ fRecMCTrackHist->GetAxis(0)->SetRangeUser(minNCrossRows,160);
+ fRecMCTrackHist->GetAxis(1)->SetRangeUser(0.,10);
+
+ h3D = (TH3D *)fRecMCTrackHist->Projection(1,9,15);
+ h3D->SetName("rec_prim_Chi2PerClust_pT_centr");
+ h3D->SetTitle("rec_prim_Chi2PerClust_pT_centr");
+ aFolderObj->Add(h3D);
+
+ // Chi2PerClustITS
+ fRecMCTrackHist->GetAxis(1)->SetRangeUser(0.,maxChi2PerClust-0.0001);
+ fRecMCTrackHist->GetAxis(2)->SetRangeUser(0.,100);
+
+ h3D = (TH3D *)fRecMCTrackHist->Projection(2,9,15);
+ h3D->SetName("rec_prim_Chi2PerClustITS_pT_centr");
+ h3D->SetTitle("rec_prim_Chi2PerClustITS_pT_centr");
+ aFolderObj->Add(h3D);
+
+
+ // NCrossRowsOverFindable
+ fRecMCTrackHist->GetAxis(2)->SetRangeUser(0.,maxChi2PerClustITS-0.0001);
+ fRecMCTrackHist->GetAxis(3)->SetRangeUser(0,1.4999);
+
+ h3D = (TH3D *)fRecMCTrackHist->Projection(3,9,15);
+ h3D->SetName("rec_prim_NCrossRowsOverFindable_pT_centr");
+ h3D->SetTitle("rec_prim_NCrossRowsOverFindable_pT_centr");
+ aFolderObj->Add(h3D);
+
+
+ // FracSharedClust
+ fRecMCTrackHist->GetAxis(3)->SetRangeUser(minNCrossRowsOverFindable,1.4999);
+ fRecMCTrackHist->GetAxis(4)->SetRangeUser(0.,0.9999);
+
+ h3D = (TH3D *)fRecMCTrackHist->Projection(4,9,15);
+ h3D->SetName("rec_prim_FracSharedClust_pT_centr");
+ h3D->SetTitle("rec_prim_FracSharedClust_pT_centr");
+ aFolderObj->Add(h3D);
+
+
+ // DCAr
+ fRecMCTrackHist->GetAxis(4)->SetRangeUser(0.,maxFracSharedClust-0.0001);
+ fRecMCTrackHist->GetAxis(5)->SetRangeUser(-10.,9.9999);
+
+ h3D = (TH3D *)fRecMCTrackHist->Projection(5,9,15);
+ h3D->SetName("rec_prim_DCAr_pT_centr");
+ h3D->SetTitle("rec_prim_DCAr_pT_centr");
+ aFolderObj->Add(h3D);
+
+
+ // DCAz
+ fRecMCTrackHist->GetAxis(5)->SetRangeUser(minDCAr,maxDCAr-0.0001);
+ fRecMCTrackHist->GetAxis(6)->SetRangeUser(-10.,9.9999);
+
+ h3D = (TH3D *)fRecMCTrackHist->Projection(6,9,15);
+ h3D->SetName("rec_prim_DCAz_pT_centr");
+ h3D->SetTitle("rec_prim_DCAz_pT_centr");
+ aFolderObj->Add(h3D);
+
+ // chi2ToTPCc
+ fRecMCTrackHist->GetAxis(6)->SetRangeUser(minDCAz,maxDCAz-0.0001);
+ fRecMCTrackHist->GetAxis(16)->SetRangeUser(0,100);
+
+ h3D = (TH3D *)fRecMCTrackHist->Projection(16,9,15);
+ h3D->SetName("rec_prim_chi2ToTPCc_pT_centr");
+ h3D->SetTitle("rec_prim_chi2ToTPCc_pT_centr");
+ aFolderObj->Add(h3D);
+
+ }
+
+ // export objects to analysis folder
+ fAnalysisFolder = ExportToFolder(aFolderObj);
+ if(!fAnalysisFolder) {
+ if(aFolderObj) delete aFolderObj;
+ return;
+ }
+
+ // Reset setings
+
+ for(Int_t ndim = 0; ndim < fEventCount->GetNdimensions(); ndim++) {
+ fEventCount->GetAxis(ndim)->SetRange(1,fEventCount->GetAxis(ndim)->GetNbins());
+ }
+
+ for(Int_t ndim = 0; ndim < fMCEventHist->GetNdimensions(); ndim++) {
+ fMCEventHist->GetAxis(ndim)->SetRange(1,fMCEventHist->GetAxis(ndim)->GetNbins());
+ }
+
+ for(Int_t ndim = 0; ndim < fRecMCEventHist->GetNdimensions(); ndim++) {
+ fRecMCEventHist->GetAxis(ndim)->SetRange(1,fRecMCEventHist->GetAxis(ndim)->GetNbins());
+ }
+
+ for(Int_t ndim = 0; ndim < fRecMCTrackHist->GetNdimensions(); ndim++)
+ {
+ fRecMCTrackHist->GetAxis(ndim)->SetRange(1,fRecMCTrackHist->GetAxis(ndim)->GetNbins());
+ }
+
+ // delete only TObjArray
+ //if(aFolderObj) delete aFolderObj;
+
+ if (fFolderObj) delete fFolderObj;
+ fFolderObj = aFolderObj;
+ aFolderObj=0;
+}
+
+//_____________________________________________________________________________
+TFolder* AlidNdPtCutAnalysisPbPb::ExportToFolder(TObjArray * const array)
+{
+ // recreate folder avery time and export objects to new one
+ //
+ AlidNdPtCutAnalysisPbPb * comp=this;
+ TFolder *folder = comp->GetAnalysisFolder();
+
+ TString name, title;
+ TFolder *newFolder = 0;
+ Int_t i = 0;
+ Int_t size = array->GetSize();
+
+ if(folder) {
+ // get name and title from old folder
+ name = folder->GetName();
+ title = folder->GetTitle();
+
+ // delete old one
+ delete folder;
+
+ // create new one
+ newFolder = CreateFolder(name.Data(),title.Data());
+ newFolder->SetOwner();
+
+ // add objects to folder
+ while(i < size) {
+ newFolder->Add(array->At(i));
+ i++;
+ }
+ }
+
+return newFolder;
+}
+
+//_____________________________________________________________________________
+TFolder* AlidNdPtCutAnalysisPbPb::CreateFolder(TString name,TString title) {
+// create folder for analysed histograms
+//
+TFolder *folder = 0;
+ folder = new TFolder(name.Data(),title.Data());
+
+ return folder;
+}