#include "TH1.h"
#include "TH2.h"
#include "TCanvas.h"
#include "THnSparse.h"
#include "AliHeader.h"
#include "AliInputEventHandler.h"
#include "AliAnalysisManager.h"
#include "AliGenEventHeader.h"
#include "AliStack.h"
#include "AliESDEvent.h"
#include "AliMCEvent.h"
#include "AliESDtrackCuts.h"
#include "AliLog.h"
#include "AliMultiplicity.h"
#include "AliTracker.h"
#include "AlidNdPtEventCuts.h"
#include "AlidNdPtAcceptanceCuts.h"
#include "AliPhysicsSelection.h"
#include "AliTriggerAnalysis.h"
#include "AliPWG0Helper.h"
#include "AlidNdPtHelper.h"
#include "AlidNdPtAnalysis.h"
using namespace std;
ClassImp(AlidNdPtAnalysis)
AlidNdPtAnalysis::AlidNdPtAnalysis(): AlidNdPt(),
fAnalysisFolder(0),
fHistogramsOn(kFALSE),
fEventMultCorrelationMatrix(0),
fTrackPtCorrelationMatrix(0),
fGenEventMatrix(0),
fGenSDEventMatrix(0),
fGenDDEventMatrix(0),
fGenNDEventMatrix(0),
fGenNSDEventMatrix(0),
fTriggerEventMatrix(0),
fTriggerSDEventMatrix(0),
fTriggerDDEventMatrix(0),
fTriggerNDEventMatrix(0),
fTriggerNSDEventMatrix(0),
fRecEventMatrix(0),
fRecSDEventMatrix(0),
fRecDDEventMatrix(0),
fRecNDEventMatrix(0),
fRecNSDEventMatrix(0),
fGenTrackEventMatrix(0),
fGenTrackSDEventMatrix(0),
fGenTrackDDEventMatrix(0),
fGenTrackNDEventMatrix(0),
fGenTrackNSDEventMatrix(0),
fTriggerTrackEventMatrix(0),
fTriggerTrackSDEventMatrix(0),
fTriggerTrackDDEventMatrix(0),
fTriggerTrackNDEventMatrix(0),
fTriggerTrackNSDEventMatrix(0),
fRecTrackEventMatrix(0),
fRecTrackSDEventMatrix(0),
fRecTrackDDEventMatrix(0),
fRecTrackNDEventMatrix(0),
fRecTrackNSDEventMatrix(0),
fGenTrackMatrix(0),
fGenPrimTrackMatrix(0),
fRecPrimTrackMatrix(0),
fRecTrackMatrix(0),
fRecSecTrackMatrix(0),
fRecMultTrackMatrix(0),
fMCEventHist1(0),
fRecEventHist1(0),
fRecEventHist2(0),
fRecMCEventHist1(0),
fRecMCEventHist2(0),
fRecMCEventHist3(0),
fRecMCTrackHist1(0),
fMCMultRecTrackHist1(0),
fRecTrackHist2(0),
fRecEventHist(0),
fRecTrackHist(0),
fEventCount(0),
fMCPrimTrackHist(0),
fRecCandleEventMatrix(0),
fMultNbins(0),
fPtNbins(0),
fPtCorrNbins(0),
fEtaNbins(0),
fZvNbins(0),
fBinsMult(0),
fBinsPt(0),
fBinsPtCorr(0),
fBinsEta(0),
fBinsZv(0),
fIsInit(kFALSE)
{
for(Int_t i=0; i<AlidNdPtHelper::kCutSteps; i++) {
fMCTrackHist1[i]=0;
fMCPrimTrackHist1[i]=0;
fMCPrimTrackHist2[i]=0;
fMCSecTrackHist1[i]=0;
fRecTrackHist1[i]=0;
fRecTrackMultHist1[i]=0;
}
}
AlidNdPtAnalysis::AlidNdPtAnalysis(Char_t* name, Char_t* title): AlidNdPt(name,title),
fAnalysisFolder(0),
fHistogramsOn(kFALSE),
fEventMultCorrelationMatrix(0),
fTrackPtCorrelationMatrix(0),
fGenEventMatrix(0),
fGenSDEventMatrix(0),
fGenDDEventMatrix(0),
fGenNDEventMatrix(0),
fGenNSDEventMatrix(0),
fTriggerEventMatrix(0),
fTriggerSDEventMatrix(0),
fTriggerDDEventMatrix(0),
fTriggerNDEventMatrix(0),
fTriggerNSDEventMatrix(0),
fRecEventMatrix(0),
fRecSDEventMatrix(0),
fRecDDEventMatrix(0),
fRecNDEventMatrix(0),
fRecNSDEventMatrix(0),
fGenTrackEventMatrix(0),
fGenTrackSDEventMatrix(0),
fGenTrackDDEventMatrix(0),
fGenTrackNDEventMatrix(0),
fGenTrackNSDEventMatrix(0),
fTriggerTrackEventMatrix(0),
fTriggerTrackSDEventMatrix(0),
fTriggerTrackDDEventMatrix(0),
fTriggerTrackNDEventMatrix(0),
fTriggerTrackNSDEventMatrix(0),
fRecTrackEventMatrix(0),
fRecTrackSDEventMatrix(0),
fRecTrackDDEventMatrix(0),
fRecTrackNDEventMatrix(0),
fRecTrackNSDEventMatrix(0),
fGenTrackMatrix(0),
fGenPrimTrackMatrix(0),
fRecPrimTrackMatrix(0),
fRecTrackMatrix(0),
fRecSecTrackMatrix(0),
fRecMultTrackMatrix(0),
fMCEventHist1(0),
fRecEventHist1(0),
fRecEventHist2(0),
fRecMCEventHist1(0),
fRecMCEventHist2(0),
fRecMCEventHist3(0),
fRecMCTrackHist1(0),
fMCMultRecTrackHist1(0),
fRecTrackHist2(0),
fRecEventHist(0),
fRecTrackHist(0),
fEventCount(0),
fMCPrimTrackHist(0),
fRecCandleEventMatrix(0),
fMultNbins(0),
fPtNbins(0),
fPtCorrNbins(0),
fEtaNbins(0),
fZvNbins(0),
fBinsMult(0),
fBinsPt(0),
fBinsPtCorr(0),
fBinsEta(0),
fBinsZv(0),
fIsInit(kFALSE)
{
for(Int_t i=0; i<AlidNdPtHelper::kCutSteps; i++) {
fMCTrackHist1[i]=0;
fMCPrimTrackHist1[i]=0;
fMCPrimTrackHist2[i]=0;
fMCSecTrackHist1[i]=0;
fRecTrackHist1[i]=0;
fRecTrackMultHist1[i]=0;
}
}
AlidNdPtAnalysis::~AlidNdPtAnalysis() {
if(fEventMultCorrelationMatrix) delete fEventMultCorrelationMatrix; fEventMultCorrelationMatrix=0;
if(fTrackPtCorrelationMatrix) delete fTrackPtCorrelationMatrix; fTrackPtCorrelationMatrix=0;
if(fGenEventMatrix) delete fGenEventMatrix; fGenEventMatrix=0;
if(fGenSDEventMatrix) delete fGenSDEventMatrix; fGenSDEventMatrix=0;
if(fGenDDEventMatrix) delete fGenDDEventMatrix; fGenDDEventMatrix=0;
if(fGenNDEventMatrix) delete fGenNDEventMatrix; fGenNDEventMatrix=0;
if(fGenNSDEventMatrix) delete fGenNSDEventMatrix; fGenNSDEventMatrix=0;
if(fTriggerEventMatrix) delete fTriggerEventMatrix; fTriggerEventMatrix=0;
if(fTriggerSDEventMatrix) delete fTriggerSDEventMatrix; fTriggerSDEventMatrix=0;
if(fTriggerDDEventMatrix) delete fTriggerDDEventMatrix; fTriggerDDEventMatrix=0;
if(fTriggerNDEventMatrix) delete fTriggerNDEventMatrix; fTriggerNDEventMatrix=0;
if(fTriggerNSDEventMatrix) delete fTriggerNSDEventMatrix; fTriggerNSDEventMatrix=0;
if(fRecEventMatrix) delete fRecEventMatrix; fRecEventMatrix=0;
if(fRecSDEventMatrix) delete fRecSDEventMatrix; fRecSDEventMatrix=0;
if(fRecDDEventMatrix) delete fRecDDEventMatrix; fRecDDEventMatrix=0;
if(fRecNDEventMatrix) delete fRecNDEventMatrix; fRecNDEventMatrix=0;
if(fRecNSDEventMatrix) delete fRecNSDEventMatrix; fRecNSDEventMatrix=0;
if(fRecCandleEventMatrix) delete fRecCandleEventMatrix; fRecCandleEventMatrix=0;
if(fGenTrackEventMatrix) delete fGenTrackEventMatrix; fGenTrackEventMatrix=0;
if(fGenTrackSDEventMatrix) delete fGenTrackSDEventMatrix; fGenTrackSDEventMatrix=0;
if(fGenTrackDDEventMatrix) delete fGenTrackDDEventMatrix; fGenTrackDDEventMatrix=0;
if(fGenTrackNDEventMatrix) delete fGenTrackNDEventMatrix; fGenTrackNDEventMatrix=0;
if(fGenTrackNSDEventMatrix) delete fGenTrackNSDEventMatrix; fGenTrackNSDEventMatrix=0;
if(fTriggerTrackEventMatrix) delete fTriggerTrackEventMatrix; fTriggerTrackEventMatrix=0;
if(fTriggerTrackSDEventMatrix) delete fTriggerTrackSDEventMatrix; fTriggerTrackSDEventMatrix=0;
if(fTriggerTrackDDEventMatrix) delete fTriggerTrackDDEventMatrix; fTriggerTrackDDEventMatrix=0;
if(fTriggerTrackNDEventMatrix) delete fTriggerTrackNDEventMatrix; fTriggerTrackNDEventMatrix=0;
if(fTriggerTrackNSDEventMatrix) delete fTriggerTrackNSDEventMatrix; fTriggerTrackNSDEventMatrix=0;
if(fRecTrackEventMatrix) delete fRecTrackEventMatrix; fRecTrackEventMatrix=0;
if(fRecTrackSDEventMatrix) delete fRecTrackSDEventMatrix; fRecTrackSDEventMatrix=0;
if(fRecTrackDDEventMatrix) delete fRecTrackDDEventMatrix; fRecTrackDDEventMatrix=0;
if(fRecTrackNDEventMatrix) delete fRecTrackNDEventMatrix; fRecTrackNDEventMatrix=0;
if(fRecTrackNSDEventMatrix) delete fRecTrackNSDEventMatrix; fRecTrackNSDEventMatrix=0;
if(fGenTrackMatrix) delete fGenTrackMatrix; fGenTrackMatrix=0;
if(fGenPrimTrackMatrix) delete fGenPrimTrackMatrix; fGenPrimTrackMatrix=0;
if(fRecPrimTrackMatrix) delete fRecPrimTrackMatrix; fRecPrimTrackMatrix=0;
if(fRecTrackMatrix) delete fRecTrackMatrix; fRecTrackMatrix=0;
if(fRecSecTrackMatrix) delete fRecSecTrackMatrix; fRecSecTrackMatrix=0;
if(fRecMultTrackMatrix) delete fRecMultTrackMatrix; fRecMultTrackMatrix=0;
if(fMCEventHist1) delete fMCEventHist1; fMCEventHist1=0;
if(fRecEventHist1) delete fRecEventHist1; fRecEventHist1=0;
if(fRecEventHist2) delete fRecEventHist2; fRecEventHist2=0;
if(fRecMCEventHist1) delete fRecMCEventHist1; fRecMCEventHist1=0;
if(fRecMCEventHist2) delete fRecMCEventHist2; fRecMCEventHist2=0;
if(fRecMCEventHist3) delete fRecMCEventHist3; fRecMCEventHist3=0;
for(Int_t i=0; i<AlidNdPtHelper::kCutSteps; i++) {
if(fMCTrackHist1[i]) delete fMCTrackHist1[i]; fMCTrackHist1[i]=0;
if(fMCPrimTrackHist1[i]) delete fMCPrimTrackHist1[i]; fMCPrimTrackHist1[i]=0;
if(fMCPrimTrackHist2[i]) delete fMCPrimTrackHist2[i]; fMCPrimTrackHist2[i]=0;
if(fMCSecTrackHist1[i]) delete fMCSecTrackHist1[i]; fMCSecTrackHist1[i]=0;
if(fRecTrackHist1[i]) delete fRecTrackHist1[i]; fRecTrackHist1[i]=0;
if(fRecTrackMultHist1[i]) delete fRecTrackMultHist1[i]; fRecTrackMultHist1[i]=0;
}
if(fRecMCTrackHist1) delete fRecMCTrackHist1; fRecMCTrackHist1=0;
if(fMCMultRecTrackHist1) delete fMCMultRecTrackHist1; fMCMultRecTrackHist1=0;
if(fRecTrackHist2) delete fRecTrackHist2; fRecTrackHist2=0;
if(fRecEventHist) delete fRecEventHist; fRecEventHist=0;
if(fRecTrackHist) delete fRecTrackHist; fRecTrackHist=0;
if(fEventCount) delete fEventCount; fEventCount=0;
if(fMCPrimTrackHist) delete fMCPrimTrackHist; fMCPrimTrackHist=0;
if(fAnalysisFolder) delete fAnalysisFolder; fAnalysisFolder=0;
if (fBinsMult) delete[] fBinsMult; fBinsMult=0;
if (fBinsPt) delete[] fBinsPt; fBinsPt=0;
if (fBinsPtCorr) delete[] fBinsPtCorr; fBinsPtCorr=0;
if (fBinsEta) delete[] fBinsEta; fBinsEta=0;
if (fBinsZv) delete[] fBinsMult; fBinsZv=0;
}
void AlidNdPtAnalysis::Init()
{
Double_t binsMultDefault[15] = {-0.5, 0.5 , 1.5 , 2.5 , 3.5 , 4.5 , 5.5 , 6.5 , 7.5 , 8.5, 9.5, 10.5, 20.5, 50.5, 150.5};
Double_t binsPtDefault[69] = {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,8.0,9.0,10.0,11.0,12.0,13.0,14.0,15.0,16.0,18.0,20.0,22.0,24.0,26.0,28.0,30.0,32.0,34.0,36.0,40.0,45.0,50.0};
Double_t binsPtCorrDefault[37] = {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,3.0,4.0,50.0};
Double_t binsEtaDefault[31] = {-1.5,-1.4,-1.3,-1.2,-1.1,-1.0,-0.9,-0.8,-0.7,-0.6,-0.5,-0.4,-0.3,-0.2,-0.1,0.,0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9,1.0,1.1,1.2,1.3,1.4,1.5};
Double_t binsZvDefault[13] = {-30.,-25.,-20.,-15.,-10.,-5.,0.,5.,10.,15.,20.,25.,30.};
if (!fBinsMult) { SetBinsMult(14,binsMultDefault); }
if (!fBinsPt) { SetBinsPt(68,binsPtDefault); }
if (!fBinsPtCorr) { SetBinsPtCorr(36,binsPtCorrDefault); }
if (!fBinsEta) { SetBinsEta(30,binsEtaDefault); }
if (!fBinsZv) { SetBinsZv(12,binsZvDefault); }
Int_t binsTrackEventCorrMatrix[4]={fZvNbins,fPtCorrNbins,fEtaNbins,fMultNbins};
Int_t binsEventHist[2]={fZvNbins,fMultNbins};
fRecEventHist = new THnSparseF("fRecEventHist","Zv:multMB",2,binsEventHist);
fRecEventHist->SetBinEdges(0,fBinsZv);
fRecEventHist->SetBinEdges(1,fBinsMult);
fRecEventHist->GetAxis(0)->SetTitle("Zv (cm)");
fRecEventHist->GetAxis(1)->SetTitle("multiplicity MB");
fRecEventHist->Sumw2();
Int_t binsTrackHist[4]={fZvNbins,fPtNbins,fEtaNbins,fMultNbins};
fRecTrackHist = new THnSparseF("fRecTrackHist","Zv:pT:eta:multRecMult",4,binsTrackHist);
fRecTrackHist->SetBinEdges(0,fBinsZv);
fRecTrackHist->SetBinEdges(1,fBinsPt);
fRecTrackHist->SetBinEdges(2,fBinsEta);
fRecTrackHist->SetBinEdges(3,fBinsMult);
fRecTrackHist->GetAxis(0)->SetTitle("Zv (cm)");
fRecTrackHist->GetAxis(1)->SetTitle("p_{T} (GeV/c)");
fRecTrackHist->GetAxis(2)->SetTitle("#eta");
fRecTrackHist->GetAxis(3)->SetTitle("multiplicity (multCuts)");
fRecTrackHist->Sumw2();
Int_t binsMCPrimTrackHist[4]={fZvNbins,fPtNbins,fEtaNbins,fMultNbins};
fMCPrimTrackHist = new THnSparseF("fMCPrimTrackHist","Zv:mcpT:mceta:multTrueMC",4,binsMCPrimTrackHist);
fMCPrimTrackHist->SetBinEdges(0,fBinsZv);
fMCPrimTrackHist->SetBinEdges(1,fBinsPt);
fMCPrimTrackHist->SetBinEdges(2,fBinsEta);
fMCPrimTrackHist->SetBinEdges(3,fBinsMult);
fMCPrimTrackHist->GetAxis(0)->SetTitle("Zv (cm)");
fMCPrimTrackHist->GetAxis(1)->SetTitle("MC p_{T} (GeV/c)");
fMCPrimTrackHist->GetAxis(2)->SetTitle("#eta (MC)");
fMCPrimTrackHist->GetAxis(3)->SetTitle("true multiplicity (MC)");
fMCPrimTrackHist->Sumw2();
Int_t binsMultTrueEventMatrix[3]={fMultNbins,fMultNbins,fMultNbins};
fEventMultCorrelationMatrix = new THnSparseF("fEventMultCorrelationMatrix","multRecMult:multTrueMC:multMB",3,binsMultTrueEventMatrix);
fEventMultCorrelationMatrix->SetBinEdges(0,fBinsMult);
fEventMultCorrelationMatrix->SetBinEdges(1,fBinsMult);
fEventMultCorrelationMatrix->SetBinEdges(2,fBinsMult);
fEventMultCorrelationMatrix->GetAxis(0)->SetTitle("multiplicity (multCuts)");
fEventMultCorrelationMatrix->GetAxis(1)->SetTitle("true multiplicity (MC)");
fEventMultCorrelationMatrix->GetAxis(2)->SetTitle("MB multiplicity");
fEventMultCorrelationMatrix->Sumw2();
Int_t binsTrackPtCorrelationMatrix[3]={fPtCorrNbins,fPtCorrNbins,fEtaNbins};
fTrackPtCorrelationMatrix = new THnSparseF("fTrackPtCorrelationMatrix","Pt:mcPt:mcEta",3,binsTrackPtCorrelationMatrix);
fTrackPtCorrelationMatrix->SetBinEdges(0,fBinsPtCorr);
fTrackPtCorrelationMatrix->SetBinEdges(1,fBinsPtCorr);
fTrackPtCorrelationMatrix->SetBinEdges(2,fBinsEta);
fTrackPtCorrelationMatrix->GetAxis(0)->SetTitle("Pt (GeV/c)");
fTrackPtCorrelationMatrix->GetAxis(1)->SetTitle("mcPt (GeV/c)");
fTrackPtCorrelationMatrix->GetAxis(2)->SetTitle("mcEta");
fTrackPtCorrelationMatrix->Sumw2();
Int_t binsEventMatrix[2]={fZvNbins,fMultNbins};
fGenEventMatrix = new THnSparseF("fGenEventMatrix","mcZv:multTrueMC",2,binsEventMatrix);
fGenEventMatrix->SetBinEdges(0,fBinsZv);
fGenEventMatrix->SetBinEdges(1,fBinsMult);
fGenEventMatrix->GetAxis(0)->SetTitle("mcZv (cm)");
fGenEventMatrix->GetAxis(1)->SetTitle("true multiplicity (MC)");
fGenEventMatrix->Sumw2();
fGenSDEventMatrix = new THnSparseF("fGenSDEventMatrix","mcZv:multTrueMC",2,binsEventMatrix);
fGenSDEventMatrix->SetBinEdges(0,fBinsZv);
fGenSDEventMatrix->SetBinEdges(1,fBinsMult);
fGenSDEventMatrix->GetAxis(0)->SetTitle("mcZv (cm)");
fGenSDEventMatrix->GetAxis(1)->SetTitle("true multiplicity (MC)");
fGenSDEventMatrix->Sumw2();
fGenDDEventMatrix = new THnSparseF("fGenDDEventMatrix","mcZv:multTrueMC",2,binsEventMatrix);
fGenDDEventMatrix->SetBinEdges(0,fBinsZv);
fGenDDEventMatrix->SetBinEdges(1,fBinsMult);
fGenDDEventMatrix->GetAxis(0)->SetTitle("mcZv (cm)");
fGenDDEventMatrix->GetAxis(1)->SetTitle("true multiplicity (MC)");
fGenDDEventMatrix->Sumw2();
fGenNDEventMatrix = new THnSparseF("fGenNDEventMatrix","mcZv:multTrueMC",2,binsEventMatrix);
fGenNDEventMatrix->SetBinEdges(0,fBinsZv);
fGenNDEventMatrix->SetBinEdges(1,fBinsMult);
fGenNDEventMatrix->GetAxis(0)->SetTitle("mcZv (cm)");
fGenNDEventMatrix->GetAxis(1)->SetTitle("true multiplicity (MC)");
fGenNDEventMatrix->Sumw2();
fGenNSDEventMatrix = new THnSparseF("fGenNSDEventMatrix","mcZv:multTrueMC",2,binsEventMatrix);
fGenNSDEventMatrix->SetBinEdges(0,fBinsZv);
fGenNSDEventMatrix->SetBinEdges(1,fBinsMult);
fGenNSDEventMatrix->GetAxis(0)->SetTitle("mcZv (cm)");
fGenNSDEventMatrix->GetAxis(1)->SetTitle("true multiplicity (MC)");
fGenNSDEventMatrix->Sumw2();
fTriggerEventMatrix = new THnSparseF("fTriggerEventMatrix","mcZv:multTrueMC",2,binsEventMatrix);
fTriggerEventMatrix->SetBinEdges(0,fBinsZv);
fTriggerEventMatrix->SetBinEdges(1,fBinsMult);
fTriggerEventMatrix->GetAxis(0)->SetTitle("mcZv (cm)");
fTriggerEventMatrix->GetAxis(1)->SetTitle("true multiplicity (MC)");
fTriggerEventMatrix->Sumw2();
fTriggerSDEventMatrix = new THnSparseF("fTriggerSDEventMatrix","mcZv:multTrueMC",2,binsEventMatrix);
fTriggerSDEventMatrix->SetBinEdges(0,fBinsZv);
fTriggerSDEventMatrix->SetBinEdges(1,fBinsMult);
fTriggerSDEventMatrix->GetAxis(0)->SetTitle("mcZv (cm)");
fTriggerSDEventMatrix->GetAxis(1)->SetTitle("true multiplicity (MC)");
fTriggerSDEventMatrix->Sumw2();
fTriggerDDEventMatrix = new THnSparseF("fTriggerDDEventMatrix","mcZv:multTrueMC",2,binsEventMatrix);
fTriggerDDEventMatrix->SetBinEdges(0,fBinsZv);
fTriggerDDEventMatrix->SetBinEdges(1,fBinsMult);
fTriggerDDEventMatrix->GetAxis(0)->SetTitle("mcZv (cm)");
fTriggerDDEventMatrix->GetAxis(1)->SetTitle("true multiplicity (MC)");
fTriggerDDEventMatrix->Sumw2();
fTriggerNDEventMatrix = new THnSparseF("fTriggerNDEventMatrix","mcZv:multTrueMC",2,binsEventMatrix);
fTriggerNDEventMatrix->SetBinEdges(0,fBinsZv);
fTriggerNDEventMatrix->SetBinEdges(1,fBinsMult);
fTriggerNDEventMatrix->GetAxis(0)->SetTitle("mcZv (cm)");
fTriggerNDEventMatrix->GetAxis(1)->SetTitle("true multiplicity (MC)");
fTriggerNDEventMatrix->Sumw2();
fTriggerNSDEventMatrix = new THnSparseF("fTriggerNSDEventMatrix","mcZv:multTrueMC",2,binsEventMatrix);
fTriggerNSDEventMatrix->SetBinEdges(0,fBinsZv);
fTriggerNSDEventMatrix->SetBinEdges(1,fBinsMult);
fTriggerNSDEventMatrix->GetAxis(0)->SetTitle("mcZv (cm)");
fTriggerNSDEventMatrix->GetAxis(1)->SetTitle("true multiplicity (MC)");
fTriggerNSDEventMatrix->Sumw2();
fRecEventMatrix = new THnSparseF("fRecEventMatrix","mcZv:multTrueMC",2,binsEventMatrix);
fRecEventMatrix->SetBinEdges(0,fBinsZv);
fRecEventMatrix->SetBinEdges(1,fBinsMult);
fRecEventMatrix->GetAxis(0)->SetTitle("mcZv (cm)");
fRecEventMatrix->GetAxis(1)->SetTitle("true multiplicity (MC)");
fRecEventMatrix->Sumw2();
fRecSDEventMatrix = new THnSparseF("fRecSDEventMatrix","mcZv:multTrueMC",2,binsEventMatrix);
fRecSDEventMatrix->SetBinEdges(0,fBinsZv);
fRecSDEventMatrix->SetBinEdges(1,fBinsMult);
fRecSDEventMatrix->GetAxis(0)->SetTitle("mcZv (cm)");
fRecSDEventMatrix->GetAxis(1)->SetTitle("true multiplicity (MC)");
fRecSDEventMatrix->Sumw2();
fRecDDEventMatrix = new THnSparseF("fRecDDEventMatrix","mcZv:multTrueMC",2,binsEventMatrix);
fRecDDEventMatrix->SetBinEdges(0,fBinsZv);
fRecDDEventMatrix->SetBinEdges(1,fBinsMult);
fRecDDEventMatrix->GetAxis(0)->SetTitle("mcZv (cm)");
fRecDDEventMatrix->GetAxis(1)->SetTitle("true multiplicity (MC)");
fRecDDEventMatrix->Sumw2();
fRecNDEventMatrix = new THnSparseF("fRecNDEventMatrix","mcZv:multTrueMC",2,binsEventMatrix);
fRecNDEventMatrix->SetBinEdges(0,fBinsZv);
fRecNDEventMatrix->SetBinEdges(1,fBinsMult);
fRecNDEventMatrix->GetAxis(0)->SetTitle("mcZv (cm)");
fRecNDEventMatrix->GetAxis(1)->SetTitle("true multiplicity (MC)");
fRecNDEventMatrix->Sumw2();
fRecNSDEventMatrix = new THnSparseF("fRecNSDEventMatrix","mcZv:multTrueMC",2,binsEventMatrix);
fRecNSDEventMatrix->SetBinEdges(0,fBinsZv);
fRecNSDEventMatrix->SetBinEdges(1,fBinsMult);
fRecNSDEventMatrix->GetAxis(0)->SetTitle("mcZv (cm)");
fRecNSDEventMatrix->GetAxis(1)->SetTitle("true multiplicity (MC)");
fRecNSDEventMatrix->Sumw2();
fRecCandleEventMatrix = new THnSparseF("fRecCandleEventMatrix","mcZv:multTrueMC",2,binsEventMatrix);
fRecCandleEventMatrix->SetBinEdges(0,fBinsZv);
fRecCandleEventMatrix->SetBinEdges(1,fBinsMult);
fRecCandleEventMatrix->GetAxis(0)->SetTitle("mcZv (cm)");
fRecCandleEventMatrix->GetAxis(1)->SetTitle("true multiplicity (MC)");
fRecCandleEventMatrix->Sumw2();
fGenTrackEventMatrix = new THnSparseF("fGenTrackEventMatrix","mcZv:mcPt:mcEta:multTrueMC",4,binsTrackEventCorrMatrix);
fGenTrackEventMatrix->SetBinEdges(0,fBinsZv);
fGenTrackEventMatrix->SetBinEdges(1,fBinsPtCorr);
fGenTrackEventMatrix->SetBinEdges(2,fBinsEta);
fGenTrackEventMatrix->SetBinEdges(3,fBinsMult);
fGenTrackEventMatrix->GetAxis(0)->SetTitle("mcZv (cm)");
fGenTrackEventMatrix->GetAxis(1)->SetTitle("mcPt (GeV/c)");
fGenTrackEventMatrix->GetAxis(2)->SetTitle("mcEta");
fGenTrackEventMatrix->GetAxis(3)->SetTitle("true multiplicity (MC)");
fGenTrackEventMatrix->Sumw2();
fGenTrackSDEventMatrix = new THnSparseF("fGenTrackSDEventMatrix","mcZv:mcPt:mcEta:multTrueMC",4,binsTrackEventCorrMatrix);
fGenTrackSDEventMatrix->SetBinEdges(0,fBinsZv);
fGenTrackSDEventMatrix->SetBinEdges(1,fBinsPtCorr);
fGenTrackSDEventMatrix->SetBinEdges(2,fBinsEta);
fGenTrackSDEventMatrix->SetBinEdges(3,fBinsMult);
fGenTrackSDEventMatrix->GetAxis(0)->SetTitle("mcZv (cm)");
fGenTrackSDEventMatrix->GetAxis(1)->SetTitle("mcPt (GeV/c)");
fGenTrackSDEventMatrix->GetAxis(2)->SetTitle("mcEta");
fGenTrackSDEventMatrix->GetAxis(3)->SetTitle("true multiplicity (MC)");
fGenTrackSDEventMatrix->Sumw2();
fGenTrackDDEventMatrix = new THnSparseF("fGenTrackDDEventMatrix","mcZv:mcPt:mcEta:multTrueMC",4,binsTrackEventCorrMatrix);
fGenTrackDDEventMatrix->SetBinEdges(0,fBinsZv);
fGenTrackDDEventMatrix->SetBinEdges(1,fBinsPtCorr);
fGenTrackDDEventMatrix->SetBinEdges(2,fBinsEta);
fGenTrackDDEventMatrix->SetBinEdges(3,fBinsMult);
fGenTrackDDEventMatrix->GetAxis(0)->SetTitle("mcZv (cm)");
fGenTrackDDEventMatrix->GetAxis(1)->SetTitle("mcPt (GeV/c)");
fGenTrackDDEventMatrix->GetAxis(2)->SetTitle("mcEta");
fGenTrackDDEventMatrix->GetAxis(3)->SetTitle("true multiplicity (MC)");
fGenTrackDDEventMatrix->Sumw2();
fGenTrackNDEventMatrix = new THnSparseF("fGenTrackNDEventMatrix","mcZv:mcPt:mcEta:multTrueMC",4,binsTrackEventCorrMatrix);
fGenTrackNDEventMatrix->SetBinEdges(0,fBinsZv);
fGenTrackNDEventMatrix->SetBinEdges(1,fBinsPtCorr);
fGenTrackNDEventMatrix->SetBinEdges(2,fBinsEta);
fGenTrackNDEventMatrix->SetBinEdges(3,fBinsMult);
fGenTrackNDEventMatrix->GetAxis(0)->SetTitle("mcZv (cm)");
fGenTrackNDEventMatrix->GetAxis(1)->SetTitle("mcPt (GeV/c)");
fGenTrackNDEventMatrix->GetAxis(2)->SetTitle("mcEta");
fGenTrackNDEventMatrix->GetAxis(3)->SetTitle("true multiplicity (MC)");
fGenTrackNDEventMatrix->Sumw2();
fGenTrackNSDEventMatrix = new THnSparseF("fGenTrackNSDEventMatrix","mcZv:mcPt:mcEta:multTrueMC",4,binsTrackEventCorrMatrix);
fGenTrackNSDEventMatrix->SetBinEdges(0,fBinsZv);
fGenTrackNSDEventMatrix->SetBinEdges(1,fBinsPtCorr);
fGenTrackNSDEventMatrix->SetBinEdges(2,fBinsEta);
fGenTrackNSDEventMatrix->SetBinEdges(3,fBinsMult);
fGenTrackNSDEventMatrix->GetAxis(0)->SetTitle("mcZv (cm)");
fGenTrackNSDEventMatrix->GetAxis(1)->SetTitle("mcPt (GeV/c)");
fGenTrackNSDEventMatrix->GetAxis(2)->SetTitle("mcEta");
fGenTrackNSDEventMatrix->GetAxis(3)->SetTitle("true multiplicity (MC)");
fGenTrackNSDEventMatrix->Sumw2();
fTriggerTrackEventMatrix = new THnSparseF("fTriggerTrackEventMatrix","mcZv:mcPt:mcEta:multTrueMC",4,binsTrackEventCorrMatrix);
fTriggerTrackEventMatrix->SetBinEdges(0,fBinsZv);
fTriggerTrackEventMatrix->SetBinEdges(1,fBinsPtCorr);
fTriggerTrackEventMatrix->SetBinEdges(2,fBinsEta);
fTriggerTrackEventMatrix->SetBinEdges(3,fBinsMult);
fTriggerTrackEventMatrix->GetAxis(0)->SetTitle("mcZv (cm)");
fTriggerTrackEventMatrix->GetAxis(1)->SetTitle("mcPt (GeV/c)");
fTriggerTrackEventMatrix->GetAxis(2)->SetTitle("mcEta");
fTriggerTrackEventMatrix->GetAxis(3)->SetTitle("true multiplicity (MC)");
fTriggerTrackEventMatrix->Sumw2();
fTriggerTrackSDEventMatrix = new THnSparseF("fTriggerTrackSDEventMatrix","mcZv:mcPt:mcEta:multTrueMC",4,binsTrackEventCorrMatrix);
fTriggerTrackSDEventMatrix->SetBinEdges(0,fBinsZv);
fTriggerTrackSDEventMatrix->SetBinEdges(1,fBinsPtCorr);
fTriggerTrackSDEventMatrix->SetBinEdges(2,fBinsEta);
fTriggerTrackSDEventMatrix->SetBinEdges(3,fBinsMult);
fTriggerTrackSDEventMatrix->GetAxis(0)->SetTitle("mcZv (cm)");
fTriggerTrackSDEventMatrix->GetAxis(1)->SetTitle("mcPt (GeV/c)");
fTriggerTrackSDEventMatrix->GetAxis(2)->SetTitle("mcEta");
fTriggerTrackSDEventMatrix->GetAxis(3)->SetTitle("true multiplicity (MC)");
fTriggerTrackSDEventMatrix->Sumw2();
fTriggerTrackDDEventMatrix = new THnSparseF("fTriggerTrackDDEventMatrix","mcZv:mcPt:mcEta:multTrueMC",4,binsTrackEventCorrMatrix);
fTriggerTrackDDEventMatrix->SetBinEdges(0,fBinsZv);
fTriggerTrackDDEventMatrix->SetBinEdges(1,fBinsPtCorr);
fTriggerTrackDDEventMatrix->SetBinEdges(2,fBinsEta);
fTriggerTrackDDEventMatrix->SetBinEdges(3,fBinsMult);
fTriggerTrackDDEventMatrix->GetAxis(0)->SetTitle("mcZv (cm)");
fTriggerTrackDDEventMatrix->GetAxis(1)->SetTitle("mcPt (GeV/c)");
fTriggerTrackDDEventMatrix->GetAxis(2)->SetTitle("mcEta");
fTriggerTrackDDEventMatrix->GetAxis(3)->SetTitle("true multiplicity (MC)");
fTriggerTrackDDEventMatrix->Sumw2();
fTriggerTrackNDEventMatrix = new THnSparseF("fTriggerTrackNDEventMatrix","mcZv:mcPt:mcEta:multTrueMC",4,binsTrackEventCorrMatrix);
fTriggerTrackNDEventMatrix->SetBinEdges(0,fBinsZv);
fTriggerTrackNDEventMatrix->SetBinEdges(1,fBinsPtCorr);
fTriggerTrackNDEventMatrix->SetBinEdges(2,fBinsEta);
fTriggerTrackNDEventMatrix->SetBinEdges(3,fBinsMult);
fTriggerTrackNDEventMatrix->GetAxis(0)->SetTitle("mcZv (cm)");
fTriggerTrackNDEventMatrix->GetAxis(1)->SetTitle("mcPt (GeV/c)");
fTriggerTrackNDEventMatrix->GetAxis(2)->SetTitle("mcEta");
fTriggerTrackNDEventMatrix->GetAxis(3)->SetTitle("true multiplicity (MC)");
fTriggerTrackNDEventMatrix->Sumw2();
fTriggerTrackNSDEventMatrix = new THnSparseF("fTriggerTrackNSDEventMatrix","mcZv:mcPt:mcEta:multTrueMC",4,binsTrackEventCorrMatrix);
fTriggerTrackNSDEventMatrix->SetBinEdges(0,fBinsZv);
fTriggerTrackNSDEventMatrix->SetBinEdges(1,fBinsPtCorr);
fTriggerTrackNSDEventMatrix->SetBinEdges(2,fBinsEta);
fTriggerTrackNSDEventMatrix->SetBinEdges(3,fBinsMult);
fTriggerTrackNSDEventMatrix->GetAxis(0)->SetTitle("mcZv (cm)");
fTriggerTrackNSDEventMatrix->GetAxis(1)->SetTitle("mcPt (GeV/c)");
fTriggerTrackNSDEventMatrix->GetAxis(2)->SetTitle("mcEta");
fTriggerTrackNSDEventMatrix->GetAxis(3)->SetTitle("true multiplicity (MC)");
fTriggerTrackNSDEventMatrix->Sumw2();
fRecTrackEventMatrix = new THnSparseF("fRecTrackEventMatrix","mcZv:mcPt:mcEta:multTrueMC",4,binsTrackEventCorrMatrix);
fRecTrackEventMatrix->SetBinEdges(0,fBinsZv);
fRecTrackEventMatrix->SetBinEdges(1,fBinsPtCorr);
fRecTrackEventMatrix->SetBinEdges(2,fBinsEta);
fRecTrackEventMatrix->SetBinEdges(3,fBinsMult);
fRecTrackEventMatrix->GetAxis(0)->SetTitle("mcZv (cm)");
fRecTrackEventMatrix->GetAxis(1)->SetTitle("mcPt (GeV/c)");
fRecTrackEventMatrix->GetAxis(2)->SetTitle("mcEta");
fRecTrackEventMatrix->GetAxis(3)->SetTitle("true multiplicity (MC)");
fRecTrackEventMatrix->Sumw2();
fRecTrackSDEventMatrix = new THnSparseF("fRecTrackSDEventMatrix","mcZv:mcPt:mcEta:multTrueMC",4,binsTrackEventCorrMatrix);
fRecTrackSDEventMatrix->SetBinEdges(0,fBinsZv);
fRecTrackSDEventMatrix->SetBinEdges(1,fBinsPtCorr);
fRecTrackSDEventMatrix->SetBinEdges(2,fBinsEta);
fRecTrackSDEventMatrix->SetBinEdges(3,fBinsMult);
fRecTrackSDEventMatrix->GetAxis(0)->SetTitle("mcZv (cm)");
fRecTrackSDEventMatrix->GetAxis(1)->SetTitle("mcPt (GeV/c)");
fRecTrackSDEventMatrix->GetAxis(2)->SetTitle("mcEta");
fRecTrackSDEventMatrix->GetAxis(3)->SetTitle("true multiplicity (MC)");
fRecTrackSDEventMatrix->Sumw2();
fRecTrackDDEventMatrix = new THnSparseF("fRecTrackDDEventMatrix","mcZv:mcPt:mcEta:multTrueMC",4,binsTrackEventCorrMatrix);
fRecTrackDDEventMatrix->SetBinEdges(0,fBinsZv);
fRecTrackDDEventMatrix->SetBinEdges(1,fBinsPtCorr);
fRecTrackDDEventMatrix->SetBinEdges(2,fBinsEta);
fRecTrackDDEventMatrix->SetBinEdges(3,fBinsMult);
fRecTrackDDEventMatrix->GetAxis(0)->SetTitle("mcZv (cm)");
fRecTrackDDEventMatrix->GetAxis(1)->SetTitle("mcPt (GeV/c)");
fRecTrackDDEventMatrix->GetAxis(2)->SetTitle("mcEta");
fRecTrackDDEventMatrix->GetAxis(3)->SetTitle("true multiplicity (MC)");
fRecTrackDDEventMatrix->Sumw2();
fRecTrackNDEventMatrix = new THnSparseF("fRecTrackNDEventMatrix","mcZv:mcPt:mcEta:multTrueMC",4,binsTrackEventCorrMatrix);
fRecTrackNDEventMatrix->SetBinEdges(0,fBinsZv);
fRecTrackNDEventMatrix->SetBinEdges(1,fBinsPtCorr);
fRecTrackNDEventMatrix->SetBinEdges(2,fBinsEta);
fRecTrackNDEventMatrix->SetBinEdges(3,fBinsMult);
fRecTrackNDEventMatrix->GetAxis(0)->SetTitle("mcZv (cm)");
fRecTrackNDEventMatrix->GetAxis(1)->SetTitle("mcPt (GeV/c)");
fRecTrackNDEventMatrix->GetAxis(2)->SetTitle("mcEta");
fRecTrackNDEventMatrix->GetAxis(3)->SetTitle("true multiplicity (MC)");
fRecTrackNDEventMatrix->Sumw2();
fRecTrackNSDEventMatrix = new THnSparseF("fRecTrackNSDEventMatrix","mcZv:mcPt:mcEta:multTrueMC",4,binsTrackEventCorrMatrix);
fRecTrackNSDEventMatrix->SetBinEdges(0,fBinsZv);
fRecTrackNSDEventMatrix->SetBinEdges(1,fBinsPtCorr);
fRecTrackNSDEventMatrix->SetBinEdges(2,fBinsEta);
fRecTrackNSDEventMatrix->SetBinEdges(3,fBinsMult);
fRecTrackNSDEventMatrix->GetAxis(0)->SetTitle("mcZv (cm)");
fRecTrackNSDEventMatrix->GetAxis(1)->SetTitle("mcPt (GeV/c)");
fRecTrackNSDEventMatrix->GetAxis(2)->SetTitle("mcEta");
fRecTrackNSDEventMatrix->GetAxis(3)->SetTitle("true multiplicity (MC)");
fRecTrackNSDEventMatrix->Sumw2();
fGenTrackMatrix = new THnSparseF("fGenTrackMatrix","mcZv:mcPt:mcEta",3,binsTrackEventCorrMatrix);
fGenTrackMatrix->SetBinEdges(0,fBinsZv);
fGenTrackMatrix->SetBinEdges(1,fBinsPtCorr);
fGenTrackMatrix->SetBinEdges(2,fBinsEta);
fGenTrackMatrix->GetAxis(0)->SetTitle("mcZv (cm)");
fGenTrackMatrix->GetAxis(1)->SetTitle("mcPt (GeV/c)");
fGenTrackMatrix->GetAxis(2)->SetTitle("mcEta");
fGenTrackMatrix->Sumw2();
fGenPrimTrackMatrix = new THnSparseF("fGenPrimTrackMatrix","mcZv:mcPt:mcEta",3,binsTrackEventCorrMatrix);
fGenPrimTrackMatrix->SetBinEdges(0,fBinsZv);
fGenPrimTrackMatrix->SetBinEdges(1,fBinsPtCorr);
fGenPrimTrackMatrix->SetBinEdges(2,fBinsEta);
fGenPrimTrackMatrix->GetAxis(0)->SetTitle("mcZv (cm)");
fGenPrimTrackMatrix->GetAxis(1)->SetTitle("mcPt (GeV/c)");
fGenPrimTrackMatrix->GetAxis(2)->SetTitle("mcEta");
fGenPrimTrackMatrix->Sumw2();
fRecPrimTrackMatrix = new THnSparseF("fRecPrimTrackMatrix","mcZv:mcPt:mcEta",3,binsTrackEventCorrMatrix);
fRecPrimTrackMatrix->SetBinEdges(0,fBinsZv);
fRecPrimTrackMatrix->SetBinEdges(1,fBinsPtCorr);
fRecPrimTrackMatrix->SetBinEdges(2,fBinsEta);
fRecPrimTrackMatrix->GetAxis(0)->SetTitle("mcZv (cm)");
fRecPrimTrackMatrix->GetAxis(1)->SetTitle("mcPt (GeV/c)");
fRecPrimTrackMatrix->GetAxis(2)->SetTitle("mcEta");
fRecPrimTrackMatrix->Sumw2();
fRecTrackMatrix = new THnSparseF("fRecTrackMatrix","mcZv:mcPt:mcEta",3,binsTrackEventCorrMatrix);
fRecTrackMatrix->SetBinEdges(0,fBinsZv);
fRecTrackMatrix->SetBinEdges(1,fBinsPtCorr);
fRecTrackMatrix->SetBinEdges(2,fBinsEta);
fRecTrackMatrix->GetAxis(0)->SetTitle("mcZv (cm)");
fRecTrackMatrix->GetAxis(1)->SetTitle("Pt (GeV/c)");
fRecTrackMatrix->GetAxis(2)->SetTitle("Eta");
fRecTrackMatrix->Sumw2();
fRecSecTrackMatrix = new THnSparseF("fRecSecTrackMatrix","mcZv:mcPt:mcEta",3,binsTrackEventCorrMatrix);
fRecSecTrackMatrix->SetBinEdges(0,fBinsZv);
fRecSecTrackMatrix->SetBinEdges(1,fBinsPtCorr);
fRecSecTrackMatrix->SetBinEdges(2,fBinsEta);
fRecSecTrackMatrix->GetAxis(0)->SetTitle("mcZv (cm)");
fRecSecTrackMatrix->GetAxis(1)->SetTitle("Pt (GeV/c)");
fRecSecTrackMatrix->GetAxis(2)->SetTitle("Eta");
fRecSecTrackMatrix->Sumw2();
fRecMultTrackMatrix = new THnSparseF("fRecMultTrackMatrix","mcZv:mcPt:mcEta",3,binsTrackEventCorrMatrix);
fRecMultTrackMatrix->SetBinEdges(0,fBinsZv);
fRecMultTrackMatrix->SetBinEdges(1,fBinsPtCorr);
fRecMultTrackMatrix->SetBinEdges(2,fBinsEta);
fRecMultTrackMatrix->GetAxis(0)->SetTitle("mcZv (cm)");
fRecMultTrackMatrix->GetAxis(1)->SetTitle("mcPt (GeV/c)");
fRecMultTrackMatrix->GetAxis(2)->SetTitle("mcEta");
fRecMultTrackMatrix->Sumw2();
Int_t binsMCEventHist1[3]={100,100,140};
Double_t minMCEventHist1[3]={-0.1,-0.1,-35.};
Double_t maxMCEventHist1[3]={0.1,0.1,35.};
fMCEventHist1 = new THnSparseF("fMCEventHist1","mcXv:mcYv:mcZv",3,binsMCEventHist1,minMCEventHist1,maxMCEventHist1);
fMCEventHist1->GetAxis(0)->SetTitle("mcXv (cm)");
fMCEventHist1->GetAxis(1)->SetTitle("mcYv (cm)");
fMCEventHist1->GetAxis(2)->SetTitle("mcZv (cm)");
fMCEventHist1->Sumw2();
Int_t binsRecEventHist1[3]={100,100,140};
Double_t minRecEventHist1[3]={-3.,-3.,-35.};
Double_t maxRecEventHist1[3]={3.,3.,35.};
fRecEventHist1 = new THnSparseF("fRecEventHist1","Xv:Yv:Zv",3,binsRecEventHist1,minRecEventHist1,maxRecEventHist1);
fRecEventHist1->GetAxis(0)->SetTitle("Xv (cm)");
fRecEventHist1->GetAxis(1)->SetTitle("Yv (cm)");
fRecEventHist1->GetAxis(2)->SetTitle("Zv (cm)");
fRecEventHist1->Sumw2();
Int_t binsRecEventHist2[3]={fZvNbins,fMultNbins,fMultNbins};
fRecEventHist2 = new THnSparseF("fRecEventHist2","Zv:multMB:multRecMult",3,binsRecEventHist2);
fRecEventHist2->SetBinEdges(0,fBinsZv);
fRecEventHist2->SetBinEdges(1,fBinsMult);
fRecEventHist2->SetBinEdges(2,fBinsMult);
fRecEventHist2->GetAxis(0)->SetTitle("Zv (cm)");
fRecEventHist2->GetAxis(1)->SetTitle("multiplicity MB");
fRecEventHist2->GetAxis(2)->SetTitle("multiplicity (multCuts)");
fRecEventHist2->Sumw2();
Double_t kFact = 0.1;
Int_t binsRecMCEventHist1[3]={100,100,100};
Double_t minRecMCEventHist1[3]={-10.0*kFact,-10.0*kFact,-10.0*kFact};
Double_t maxRecMCEventHist1[3]={10.0*kFact,10.0*kFact,10.0*kFact};
fRecMCEventHist1 = new THnSparseF("fRecMCEventHist1","Xv-mcXv:Yv-mcYv:Zv-mcZv",3,binsRecMCEventHist1,minRecMCEventHist1,maxRecMCEventHist1);
fRecMCEventHist1->GetAxis(0)->SetTitle("Xv-mcXv (cm)");
fRecMCEventHist1->GetAxis(1)->SetTitle("Yv-mcYv (cm)");
fRecMCEventHist1->GetAxis(2)->SetTitle("Zv-mcZv (cm)");
fRecMCEventHist1->Sumw2();
Int_t binsRecMCEventHist2[3]={100,100,fMultNbins};
Double_t minRecMCEventHist2[3]={-10.0*kFact,-10.0*kFact,0.0};
Double_t maxRecMCEventHist2[3]={10.0*kFact,10.0*kFact,149.50};
fRecMCEventHist2 = new THnSparseF("fRecMCEventHist2","Xv-mcXv:Zv-mcZv:multMB",3,binsRecMCEventHist2,minRecMCEventHist2,maxRecMCEventHist2);
fRecMCEventHist2->SetBinEdges(2,fBinsMult);
fRecMCEventHist2->GetAxis(0)->SetTitle("Xv-mcXv (cm)");
fRecMCEventHist2->GetAxis(1)->SetTitle("Zv-mcZv (cm)");
fRecMCEventHist2->GetAxis(2)->SetTitle("multiplicity MB");
fRecMCEventHist2->Sumw2();
Int_t binsRecMCEventHist3[2]={fMultNbins,5};
Double_t minRecMCEventHist3[2]={-0.5,0.0};
Double_t maxRecMCEventHist3[2]={149.50,5.0};
fRecMCEventHist3 = new THnSparseF("fRecMCEventHist3","multRecMult:EventType (ND, DD, SD)",2,binsRecMCEventHist3,minRecMCEventHist3,maxRecMCEventHist3);
fRecMCEventHist3->SetBinEdges(0,fBinsMult);
fRecMCEventHist3->GetAxis(0)->SetTitle("multiplicity (multCuts)");
fRecMCEventHist3->GetAxis(1)->SetTitle("EventType");
fRecMCEventHist3->Sumw2();
char name[256];
char title[256];
for(Int_t i=0; i<AlidNdPtHelper::kCutSteps; i++)
{
Int_t binsMCTrackHist1[3]= {fPtCorrNbins, fEtaNbins, 90};
Double_t minMCTrackHist1[3]={0.,-1.,0.};
Double_t maxMCTrackHist1[3]={10.,1.,2.*TMath::Pi()};
snprintf(name,256,"fMCTrackHist1_%d",i);
snprintf(title,256,"mcPt:mcEta:mcPhi");
fMCTrackHist1[i] = new THnSparseF(name,title,3,binsMCTrackHist1,minMCTrackHist1,maxMCTrackHist1);
fMCTrackHist1[i]->SetBinEdges(0,fBinsPtCorr);
fMCTrackHist1[i]->SetBinEdges(1,fBinsEta);
fMCTrackHist1[i]->GetAxis(0)->SetTitle("mcPt (GeV/c)");
fMCTrackHist1[i]->GetAxis(1)->SetTitle("mcEta");
fMCTrackHist1[i]->GetAxis(2)->SetTitle("mcPhi (rad)");
fMCTrackHist1[i]->Sumw2();
Int_t binsMCPrimTrackHist1[5]= {fPtCorrNbins,fEtaNbins,6,20,4000};
Double_t minMCPrimTrackHist1[5]={0.,-1.,0.,0.,0.};
Double_t maxMCPrimTrackHist1[5]={10.,1.,6.,20.,4000.};
snprintf(name,256,"fMCPrimTrackHist1_%d",i);
snprintf(title,256,"mcPt:mcEta:pid:mech:mother");
fMCPrimTrackHist1[i] = new THnSparseF(name,title,5,binsMCPrimTrackHist1,minMCPrimTrackHist1,maxMCPrimTrackHist1);
fMCPrimTrackHist1[i]->SetBinEdges(0,fBinsPtCorr);
fMCPrimTrackHist1[i]->SetBinEdges(1,fBinsEta);
fMCPrimTrackHist1[i]->GetAxis(0)->SetTitle("mcPt (GeV/c)");
fMCPrimTrackHist1[i]->GetAxis(1)->SetTitle("mcEta");
fMCPrimTrackHist1[i]->GetAxis(2)->SetTitle("pid");
fMCPrimTrackHist1[i]->GetAxis(3)->SetTitle("mech");
fMCPrimTrackHist1[i]->GetAxis(4)->SetTitle("mother");
fMCPrimTrackHist1[i]->Sumw2();
Int_t binsMCPrimTrackHist2[3]= {4000,20,4000};
Double_t minMCPrimTrackHist2[3]={0.,0.,0.};
Double_t maxMCPrimTrackHist2[3]={4000.,20.,4000.};
snprintf(name,256,"fMCPrimTrackHist2_%d",i);
snprintf(title,256,"pdg:mech:mother");
fMCPrimTrackHist2[i] = new THnSparseF(name,title,3,binsMCPrimTrackHist2,minMCPrimTrackHist2,maxMCPrimTrackHist2);
fMCPrimTrackHist2[i]->GetAxis(0)->SetTitle("pdg");
fMCPrimTrackHist2[i]->GetAxis(1)->SetTitle("mech");
fMCPrimTrackHist2[i]->GetAxis(2)->SetTitle("mother");
fMCPrimTrackHist2[i]->Sumw2();
Int_t binsMCSecTrackHist1[5]= {fPtCorrNbins,fEtaNbins,6,20,4000};
Double_t minMCSecTrackHist1[5]={0.,-1.,0.,0.,0.};
Double_t maxMCSecTrackHist1[5]={10.,1.,6.,20.,4000.};
snprintf(name,256,"fMCSecTrackHist1_%d",i);
snprintf(title,256,"mcPt:mcEta:pid:mech:mother");
fMCSecTrackHist1[i] = new THnSparseF(name,title,5,binsMCSecTrackHist1,minMCSecTrackHist1,maxMCSecTrackHist1);
fMCSecTrackHist1[i]->SetBinEdges(0,fBinsPtCorr);
fMCSecTrackHist1[i]->SetBinEdges(1,fBinsEta);
fMCSecTrackHist1[i]->GetAxis(0)->SetTitle("mcPt (GeV/c)");
fMCSecTrackHist1[i]->GetAxis(1)->SetTitle("mcEta");
fMCSecTrackHist1[i]->GetAxis(2)->SetTitle("pid");
fMCSecTrackHist1[i]->GetAxis(3)->SetTitle("mech");
fMCSecTrackHist1[i]->GetAxis(4)->SetTitle("mother");
fMCSecTrackHist1[i]->Sumw2();
Int_t binsRecTrackHist1[3]={fPtNbins,fEtaNbins,90};
Double_t minRecTrackHist1[3]={0.,-1.,0.};
Double_t maxRecTrackHist1[3]={10.,1.,2.*TMath::Pi()};
snprintf(name,256,"fRecTrackHist1_%d",i);
snprintf(title,256,"Pt:Eta:Phi");
fRecTrackHist1[i] = new THnSparseF(name,title,3,binsRecTrackHist1,minRecTrackHist1,maxRecTrackHist1);
fRecTrackHist1[i]->SetBinEdges(0,fBinsPt);
fRecTrackHist1[i]->SetBinEdges(1,fBinsEta);
fRecTrackHist1[i]->GetAxis(0)->SetTitle("p_{T} (GeV/c)");
fRecTrackHist1[i]->GetAxis(1)->SetTitle("#eta");
fRecTrackHist1[i]->GetAxis(2)->SetTitle("#phi (rad)");
fRecTrackHist1[i]->Sumw2();
Int_t binsRecTrackMultHist1[2]={fPtNbins,fMultNbins};
snprintf(name,256,"fRecTrackMultHist_%d",i);
snprintf(title,256,"Pt:Mult");
fRecTrackMultHist1[i] = new THnSparseF(name,title,2,binsRecTrackMultHist1);
fRecTrackMultHist1[i]->SetBinEdges(0,fBinsPt);
fRecTrackMultHist1[i]->SetBinEdges(1,fBinsMult);
fRecTrackMultHist1[i]->GetAxis(0)->SetTitle("Pt (GeV/c)");
fRecTrackMultHist1[i]->GetAxis(1)->SetTitle("multiplicity");
fRecTrackMultHist1[i]->Sumw2();
}
Int_t binsRecMCTrackHist1[4] = {fPtCorrNbins,fEtaNbins,100,100};
Double_t minRecMCTrackHist1[4]={0.,-1.,-0.5,-0.5};
Double_t maxRecMCTrackHist1[4]={20.,1.,0.5,0.5};
snprintf(name,256,"fRecMCTrackHist1");
snprintf(title,256,"mcPt:mcEta:(Pt-mcPt)/mcPt:(Eta-mcEta)");
fRecMCTrackHist1 = new THnSparseF(name,title,4,binsRecMCTrackHist1,minRecMCTrackHist1,maxRecMCTrackHist1);
fRecMCTrackHist1->SetBinEdges(0,fBinsPtCorr);
fRecMCTrackHist1->SetBinEdges(1,fBinsEta);
fRecMCTrackHist1->GetAxis(0)->SetTitle("mcPt (GeV/c)");
fRecMCTrackHist1->GetAxis(1)->SetTitle("mcEta");
fRecMCTrackHist1->GetAxis(2)->SetTitle("(Pt-mcPt)/mcPt");
fRecMCTrackHist1->GetAxis(3)->SetTitle("Eta-mcEta");
Int_t binsMCMultRecTrackHist1[3] = {fPtCorrNbins,fEtaNbins,6};
Double_t minMCMultRecTrackHist1[3]={0.,-1.,0.};
Double_t maxMCMultRecTrackHist1[3]={20.,1.,6.};
snprintf(name,256,"fMCMultRecTrackHist1");
snprintf(title,256,"mcPt:mcEta:pid");
fMCMultRecTrackHist1 = new THnSparseF(name,title,3,binsMCMultRecTrackHist1,minMCMultRecTrackHist1,maxMCMultRecTrackHist1);
fMCMultRecTrackHist1->SetBinEdges(0,fBinsPtCorr);
fMCMultRecTrackHist1->SetBinEdges(1,fBinsEta);
fMCMultRecTrackHist1->GetAxis(0)->SetTitle("mcPt (GeV/c)");
fMCMultRecTrackHist1->GetAxis(1)->SetTitle("mcEta");
fMCMultRecTrackHist1->GetAxis(2)->SetTitle("pid");
Int_t binsRecTrackHist2[5]={160,100,fPtNbins,fEtaNbins,90};
Double_t minRecTrackHist2[5]={0., 0., 0, -1.5, 0.};
Double_t maxRecRecTrackHist2[5]={160.,10., 16, 1.5, 2.*TMath::Pi()};
fRecTrackHist2 = new THnSparseF("fRecTrackHist2","nClust:chi2PerClust:pt:eta:phi",5,binsRecTrackHist2,minRecTrackHist2,maxRecRecTrackHist2);
fRecTrackHist2->SetBinEdges(2,fBinsPt);
fRecTrackHist2->SetBinEdges(3,fBinsEta);
fRecTrackHist2->GetAxis(0)->SetTitle("nClust");
fRecTrackHist2->GetAxis(1)->SetTitle("chi2PerClust");
fRecTrackHist2->GetAxis(2)->SetTitle("p_{T} (GeV/c)");
fRecTrackHist2->GetAxis(3)->SetTitle("#eta");
fRecTrackHist2->GetAxis(4)->SetTitle("#phi (rad)");
fRecTrackHist2->Sumw2();
Int_t binsEventCount[3]={2,2,2};
Double_t minEventCount[3]={0,0,0};
Double_t maxEventCount[3]={2,2,2};
fEventCount = new THnSparseF("fEventCount","trig vs trig+vertex",3,binsEventCount,minEventCount,maxEventCount);
fEventCount->GetAxis(0)->SetTitle("trig");
fEventCount->GetAxis(1)->SetTitle("trig+vert");
fEventCount->GetAxis(2)->SetTitle("selected");
fEventCount->Sumw2();
fAnalysisFolder = CreateFolder("folderdNdPt","Analysis dNdPt Folder");
fIsInit = kTRUE;
}
void AlidNdPtAnalysis::Process(AliESDEvent *const esdEvent, AliMCEvent *const mcEvent)
{
if (!fIsInit) { Init(); }
if(!esdEvent) {
AliDebug(AliLog::kError, "esdEvent not available");
return;
}
AlidNdPtEventCuts *evtCuts = GetEventCuts();
AlidNdPtAcceptanceCuts *accCuts = GetAcceptanceCuts();
AliESDtrackCuts *esdTrackCuts = GetTrackCuts();
AlidNdPtAcceptanceCuts *multAccCuts = GetMultAcceptanceCuts();
if(!evtCuts || !accCuts || !esdTrackCuts) {
AliDebug(AliLog::kError, "cuts not available");
return;
}
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())
{
isEventTriggered = inputHandler->IsEventSelected() & GetTriggerMask();
physicsSelection = static_cast<AliPhysicsSelection*> (inputHandler->GetEventSelection());
if(!physicsSelection) return;
if (isEventTriggered && (GetTrigger() == AliTriggerAnalysis::kV0AND)) {
triggerAnalysis = physicsSelection->GetTriggerAnalysis();
if(!triggerAnalysis) return;
isEventTriggered = triggerAnalysis->IsOfflineTriggerFired(esdEvent, GetTrigger());
}
if(!IsUseMCInfo())
{
}
}
AliHeader* header = 0;
AliGenEventHeader* genHeader = 0;
AliStack* stack = 0;
TArrayF vtxMC(3);
AliPWG0Helper::MCProcessType evtType = AliPWG0Helper::kInvalidProcess;
Int_t multMCTrueTracks = 0;
if(IsUseMCInfo())
{
if(!mcEvent) {
AliDebug(AliLog::kError, "mcEvent not available");
return;
}
header = mcEvent->Header();
if (!header) {
AliDebug(AliLog::kError, "Header not available");
return;
}
stack = mcEvent->Stack();
if (!stack) {
AliDebug(AliLog::kError, "Stack not available");
return;
}
evtType = AliPWG0Helper::GetEventProcessType(header);
AliDebug(AliLog::kDebug+1, Form("Found process type %d", evtType));
genHeader = header->GenEventHeader();
if (!genHeader) {
AliDebug(AliLog::kError, "Could not retrieve genHeader from Header");
return;
}
genHeader->PrimaryVertex(vtxMC);
Double_t vMCEventHist1[3]={vtxMC[0],vtxMC[1],vtxMC[2]};
fMCEventHist1->Fill(vMCEventHist1);
multMCTrueTracks = AlidNdPtHelper::GetMCTrueTrackMult(mcEvent,evtCuts,multAccCuts);
}
const AliESDVertex* vtxESD = 0;
Bool_t isRecVertex = kFALSE;
if(evtCuts->IsRecVertexRequired())
{
Bool_t bRedoTPCVertex = evtCuts->IsRedoTPCVertex();
Bool_t bUseConstraints = evtCuts->IsUseBeamSpotConstraint();
vtxESD = AlidNdPtHelper::GetVertex(esdEvent,evtCuts,accCuts,esdTrackCuts,GetAnalysisMode(),kFALSE,bRedoTPCVertex,bUseConstraints);
if(!vtxESD) return;
isRecVertex = AlidNdPtHelper::TestRecVertex(vtxESD, esdEvent->GetPrimaryVertexSPD(), GetAnalysisMode(), kFALSE);
}
if( IsUseMCInfo() && !evtCuts->IsRecVertexRequired() ) {
vtxESD = new AliESDVertex(vtxMC[2],10.,genHeader->NProduced(),"smearMC");
if(!vtxESD) return;
isRecVertex = kTRUE;
}
if(!vtxESD) return;
Bool_t isEventOK = evtCuts->AcceptEvent(esdEvent,mcEvent,vtxESD) && isRecVertex;
Bool_t isTrigAndVertex = isEventTriggered && isEventOK;
Double_t vEventCount[3] = { static_cast<Double_t>((isEventTriggered && kTRUE)) , static_cast<Double_t>(isTrigAndVertex), static_cast<Double_t>(isTrigAndVertex && (TMath::Abs(vtxESD->GetZ()) < 10.))};
fEventCount->Fill(vEventCount);
Int_t multMBTracks = 0;
if(GetAnalysisMode() == AlidNdPtHelper::kTPC)
{
if(vtxESD && vtxESD->GetStatus() && isRecVertex)
multMBTracks = AlidNdPtHelper::GetTPCMBTrackMult(esdEvent,evtCuts,accCuts,esdTrackCuts);
}
else if( GetAnalysisMode() == AlidNdPtHelper::kTPCSPDvtx ||
GetAnalysisMode()==AlidNdPtHelper::kTPCSPDvtxUpdate )
{
const AliMultiplicity* mult = esdEvent->GetMultiplicity();
if(mult && vtxESD && vtxESD->GetStatus() && isRecVertex)
multMBTracks = mult->GetNumberOfTracklets();
}
else if( GetAnalysisMode() == AlidNdPtHelper::kTPCITS ||
GetAnalysisMode()==AlidNdPtHelper::kTPCITSHybrid ||
GetAnalysisMode() == AlidNdPtHelper::kTPCTrackSPDvtx ||
GetAnalysisMode() == AlidNdPtHelper::kTPCTrackSPDvtxUpdate ||
GetAnalysisMode() == AlidNdPtHelper::kTPCITSHybridTrackSPDvtx ||
GetAnalysisMode() == AlidNdPtHelper::kTPCITSHybridTrackSPDvtxDCArPt ||
GetAnalysisMode() == AlidNdPtHelper::kITSStandAloneTrackSPDvtx ||
GetAnalysisMode() == AlidNdPtHelper::kITSStandAloneTPCTrackSPDvtx
)
{
if(vtxESD && vtxESD->GetStatus() && isRecVertex)
multMBTracks = vtxESD->GetNContributors();
}
else {
AliDebug(AliLog::kError, Form("Found analysis type %d", GetAnalysisMode()));
return;
}
Bool_t isEventSelected = kTRUE;
if(evtCuts->IsEventSelectedRequired())
{
isEventSelected = AlidNdPtHelper::SelectEvent(esdEvent,esdTrackCuts);
}
TObjArray *allChargedTracks=0;
Int_t multAll=0, multRec=0;
Int_t *labelsAll=0, *labelsAcc=0, *labelsRec=0;
Int_t multRecMult=0;
if(isEventOK && isEventTriggered && isEventSelected)
{
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(track->Charge()==0) continue;
if(GetParticleMode() == AlidNdPtHelper::kPlus && track->Charge() < 0)
continue;
if(GetParticleMode() == AlidNdPtHelper::kMinus && track->Charge() > 0)
continue;
if(GetMultTrackCuts()->AcceptTrack(track))
{
if(GetMultAcceptanceCuts()->AcceptTrack(track)) multRecMult++;
}
}
labelsAll = new Int_t[entries];
labelsAcc = new Int_t[entries];
labelsRec = new Int_t[entries];
for(Int_t i=0; i<entries;++i)
{
AliESDtrack *track = (AliESDtrack*)allChargedTracks->At(i);
if(!track) continue;
if(track->Charge()==0) continue;
if(GetParticleMode() == AlidNdPtHelper::kPlus && track->Charge() < 0)
continue;
if(GetParticleMode() == AlidNdPtHelper::kMinus && track->Charge() > 0)
continue;
FillHistograms(track,stack,vtxESD->GetZ(),AlidNdPtHelper::kAllTracks, multRecMult);
labelsAll[multAll] = TMath::Abs(track->GetLabel());
multAll++;
if(GetAnalysisMode() == AlidNdPtHelper::kITSStandAloneTrackSPDvtx)
{
if(!(track->GetStatus() & AliESDtrack::kITSpureSA)) continue;
if(!(track->GetStatus() & AliESDtrack::kITSrefit)) continue;
if(track->GetNcls(0)<4) continue;
if(!track->HasPointOnITSLayer(0) && !track->HasPointOnITSLayer(1)) continue;
}
else if(GetAnalysisMode() == AlidNdPtHelper::kITSStandAloneTPCTrackSPDvtx)
{
if(!(track->GetStatus() & AliESDtrack::kITSpureSA)) continue;
if(!(track->GetStatus() & AliESDtrack::kITSrefit)) continue;
if(track->GetNcls(0)<4) continue;
if(!track->HasPointOnITSLayer(0) && !track->HasPointOnITSLayer(1)) continue;
Bool_t hasMatch = kFALSE;
for(Int_t j=0; j<entries;++j)
{
if(i==j) continue;
AliESDtrack *track2 = (AliESDtrack*)allChargedTracks->At(j);
if(!track2) continue;
if(track2->Charge()==0) continue;
if (!track2->GetTPCInnerParam()) continue;
if(!esdTrackCuts->AcceptTrack(track2)) continue;
AliExternalTrackParam *innerTPC = new AliExternalTrackParam(*(track2->GetTPCInnerParam()));
if(!innerTPC) continue;
Double_t x2[3]; track2->GetXYZ(x2);
Double_t b2[3]; AliTracker::GetBxByBz(x2,b2);
Double_t dz[2],cov[3];
Bool_t isPropOK = innerTPC->PropagateToDCABxByBz(vtxESD,b2,kVeryBig,dz,cov);
if(!isPropOK && innerTPC) {delete innerTPC; continue;}
if (TMath::Abs(track->GetY() - innerTPC->GetY()) > 3) { delete innerTPC; continue; }
if (TMath::Abs(track->GetSnp() - innerTPC->GetSnp()) > 0.2) { delete innerTPC; continue; }
if (TMath::Abs(track->GetTgl() - innerTPC->GetTgl()) > 0.2) { delete innerTPC; continue; }
hasMatch = kTRUE;
if(innerTPC) delete innerTPC;
}
if(!hasMatch) continue;
}
else {
if(!esdTrackCuts->AcceptTrack(track))
continue;
}
Bool_t isOK = kFALSE;
Double_t x[3]; track->GetXYZ(x);
Double_t b[3]; AliTracker::GetBxByBz(x,b);
if( GetAnalysisMode() == AlidNdPtHelper::kTPCITSHybrid ||
GetAnalysisMode() == AlidNdPtHelper::kTPCITSHybridTrackSPDvtx ||
GetAnalysisMode() == AlidNdPtHelper::kTPCITSHybridTrackSPDvtxDCArPt)
{
isOK = track->RelateToVertexTPCBxByBz(vtxESD, b, kVeryBig);
if(!isOK) continue;
AliExternalTrackParam *tpcTrack = new AliExternalTrackParam(*(track->GetTPCInnerParam()));
if (tpcTrack) {
track->Set(tpcTrack->GetX(),tpcTrack->GetAlpha(),tpcTrack->GetParameter(),tpcTrack->GetCovariance());
}
if(tpcTrack) delete tpcTrack;
}
if( GetAnalysisMode() == AlidNdPtHelper::kTPCSPDvtxUpdate ||
GetAnalysisMode() == AlidNdPtHelper::kTPCTrackSPDvtxUpdate )
{
AliExternalTrackParam cParam;
isOK = track->RelateToVertexTPCBxByBz(vtxESD, b, kVeryBig, &cParam);
if(!isOK) continue;
track->Set(cParam.GetX(),cParam.GetAlpha(),cParam.GetParameter(),cParam.GetCovariance());
if(accCuts->AcceptTrack(track)) {
FillHistograms(track,stack,vtxESD->GetZ(),AlidNdPtHelper::kRecTracks,multRecMult);
labelsRec[multRec] = TMath::Abs(track->GetLabel());
multRec++;
}
}
else {
if(accCuts->AcceptTrack(track))
{
FillHistograms(track,stack,vtxESD->GetZ(),AlidNdPtHelper::kRecTracks,multRecMult);
labelsRec[multRec] = TMath::Abs(track->GetLabel());
multRec++;
}
}
}
Double_t vRecEventHist1[3] = {vtxESD->GetX(),vtxESD->GetY(),vtxESD->GetZ()};
fRecEventHist1->Fill(vRecEventHist1);
Double_t vRecEventHist2[3] = {vtxESD->GetZ(),static_cast<Double_t>(multMBTracks),static_cast<Double_t>(multRecMult)};
fRecEventHist2->Fill(vRecEventHist2);
Double_t vRecEventHist[2] = {vtxESD->GetZ(),static_cast<Double_t>(multMBTracks)};
fRecEventHist->Fill(vRecEventHist);
}
if(IsUseMCInfo())
{
if(mcEvent) {
if(evtCuts->IsEventSelectedRequired())
{
Bool_t isMCEventSelected = AlidNdPtHelper::SelectMCEvent(mcEvent);
if(!isMCEventSelected) {
if(allChargedTracks) delete allChargedTracks; allChargedTracks = 0;
if(labelsAll) delete [] labelsAll; labelsAll = 0;
if(labelsAcc) delete [] labelsAcc; labelsAcc = 0;
if(labelsRec) delete [] labelsRec; labelsRec = 0;
return;
}
}
Double_t vMultTrueEventMatrix[3] = { static_cast<Double_t>(multRecMult), static_cast<Double_t>(multMCTrueTracks), static_cast<Double_t>(multMBTracks)};
if(isEventOK && isEventTriggered) {
if(TMath::Abs(vtxMC[2]) < 10.0)
fEventMultCorrelationMatrix->Fill(vMultTrueEventMatrix);
}
Double_t vEventMatrix[2] = {vtxMC[2],static_cast<Double_t>(multMCTrueTracks)};
fGenEventMatrix->Fill(vEventMatrix);
if(isEventTriggered) fTriggerEventMatrix->Fill(vEventMatrix);
if(isEventOK && isEventTriggered) fRecEventMatrix->Fill(vEventMatrix);
if(evtType == AliPWG0Helper::kSD) {
fGenSDEventMatrix->Fill(vEventMatrix);
if(isEventTriggered) fTriggerSDEventMatrix->Fill(vEventMatrix);
if(isEventOK && isEventTriggered) fRecSDEventMatrix->Fill(vEventMatrix);
}
if(evtType == AliPWG0Helper::kDD) {
fGenDDEventMatrix->Fill(vEventMatrix);
if(isEventTriggered) fTriggerDDEventMatrix->Fill(vEventMatrix);
if(isEventOK && isEventTriggered) fRecDDEventMatrix->Fill(vEventMatrix);
}
if(evtType == AliPWG0Helper::kND) {
fGenNDEventMatrix->Fill(vEventMatrix);
if(isEventTriggered) fTriggerNDEventMatrix->Fill(vEventMatrix);
if(isEventOK && isEventTriggered) fRecNDEventMatrix->Fill(vEventMatrix);
}
if(evtType != AliPWG0Helper::kSD) {
fGenNSDEventMatrix->Fill(vEventMatrix);
if(isEventTriggered) fTriggerNSDEventMatrix->Fill(vEventMatrix);
if(isEventOK && isEventTriggered) fRecNSDEventMatrix->Fill(vEventMatrix);
}
for (Int_t iMc = 0; iMc < stack->GetNtrack(); ++iMc)
{
TParticle* particle = stack->Particle(iMc);
if (!particle)
continue;
if(!particle->GetPDG()) continue;
Double_t charge = particle->GetPDG()->Charge()/3.;
if ( TMath::Abs(charge) < 0.001 )
continue;
if(GetParticleMode() == AlidNdPtHelper::kPlus && charge < 0.)
continue;
if(GetParticleMode() == AlidNdPtHelper::kMinus && charge > 0.)
continue;
Bool_t prim = stack->IsPhysicalPrimary(iMc);
if(!prim) continue;
if(accCuts->AcceptTrack(particle))
{
Double_t vTrackEventMatrix[4] = {vtxMC[2], particle->Pt(), particle->Eta(),static_cast<Double_t>(multMCTrueTracks)};
fGenTrackEventMatrix->Fill(vTrackEventMatrix);
if(evtType == AliPWG0Helper::kSD) {
fGenTrackSDEventMatrix->Fill(vTrackEventMatrix);
}
if(evtType == AliPWG0Helper::kDD) {
fGenTrackDDEventMatrix->Fill(vTrackEventMatrix);
}
if(evtType == AliPWG0Helper::kND) {
fGenTrackNDEventMatrix->Fill(vTrackEventMatrix);
}
if(evtType != AliPWG0Helper::kSD) {
fGenTrackNSDEventMatrix->Fill(vTrackEventMatrix);
}
if(!isEventTriggered) continue;
fTriggerTrackEventMatrix->Fill(vTrackEventMatrix);
if(evtType == AliPWG0Helper::kSD) {
fTriggerTrackSDEventMatrix->Fill(vTrackEventMatrix);
}
if(evtType == AliPWG0Helper::kDD) {
fTriggerTrackDDEventMatrix->Fill(vTrackEventMatrix);
}
if(evtType == AliPWG0Helper::kND) {
fTriggerTrackNDEventMatrix->Fill(vTrackEventMatrix);
}
if(evtType != AliPWG0Helper::kSD) {
fTriggerTrackNSDEventMatrix->Fill(vTrackEventMatrix);
}
if(!isEventOK) continue;
fRecTrackEventMatrix->Fill(vTrackEventMatrix);
if(evtType == AliPWG0Helper::kSD) {
fRecTrackSDEventMatrix->Fill(vTrackEventMatrix);
}
if(evtType == AliPWG0Helper::kDD) {
fRecTrackDDEventMatrix->Fill(vTrackEventMatrix);
}
if(evtType == AliPWG0Helper::kND) {
fRecTrackNDEventMatrix->Fill(vTrackEventMatrix);
}
if(evtType != AliPWG0Helper::kSD) {
fRecTrackNSDEventMatrix->Fill(vTrackEventMatrix);
}
}
}
if(isEventOK && isEventTriggered)
{
if(fHistogramsOn) {
Double_t vRecMCEventHist1[3] = {vtxESD->GetX()-vtxMC[0],vtxESD->GetY()-vtxMC[1],vtxESD->GetZ()-vtxMC[2]};
fRecMCEventHist1->Fill(vRecMCEventHist1);
Double_t vRecMCEventHist2[3] = {vtxESD->GetX()-vtxMC[0],vtxESD->GetZ()-vtxMC[2],static_cast<Double_t>(multMBTracks)};
fRecMCEventHist2->Fill(vRecMCEventHist2);
Double_t vRecMCEventHist3[2] = {static_cast<Double_t>(multRecMult),evtType};
fRecMCEventHist3->Fill(vRecMCEventHist3);
}
Int_t countRecCandle = 0;
for (Int_t iMc = 0; iMc < stack->GetNtrack(); ++iMc)
{
TParticle* particle = stack->Particle(iMc);
if (!particle)
continue;
Double_t vTrackMatrix[3] = {vtxMC[2],particle->Pt(),particle->Eta()};
if(!particle->GetPDG()) continue;
Double_t charge = particle->GetPDG()->Charge()/3.;
if (TMath::Abs(charge) < 0.001)
continue;
if(GetParticleMode() == AlidNdPtHelper::kPlus && charge < 0.)
continue;
if(GetParticleMode() == AlidNdPtHelper::kMinus && charge > 0.)
continue;
Bool_t prim = stack->IsPhysicalPrimary(iMc);
if(accCuts->AcceptTrack(particle))
{
if( AlidNdPtHelper::IsPrimaryParticle(stack, iMc, GetParticleMode()) ) {
fGenPrimTrackMatrix->Fill(vTrackMatrix);
Double_t vMCPrimTrackHist[4] = {vtxMC[2],particle->Pt(),particle->Eta(),static_cast<Double_t>(multMCTrueTracks)};
fMCPrimTrackHist->Fill(vMCPrimTrackHist);
}
if(fHistogramsOn)
FillHistograms(stack,iMc,AlidNdPtHelper::kAccTracks);
Int_t multCount = 0;
for(Int_t iRec=0; iRec<multRec; ++iRec)
{
if(iMc == labelsRec[iRec])
{
multCount++;
if(multCount>1)
{
fRecMultTrackMatrix->Fill(vTrackMatrix);
if(fHistogramsOn) {
Int_t pid = AlidNdPtHelper::ConvertPdgToPid(particle);
Double_t vMCMultRecTrackHist1[3] = {particle->Pt(), particle->Eta(), static_cast<Double_t>(pid)};
fMCMultRecTrackHist1->Fill(vMCMultRecTrackHist1);
}
}
}
}
for(Int_t iRec=0; iRec<multRec; ++iRec)
{
if(iMc == labelsRec[iRec])
{
fRecTrackMatrix->Fill(vTrackMatrix);
if( AlidNdPtHelper::IsPrimaryParticle(stack, iMc, GetParticleMode()) ) {
fRecPrimTrackMatrix->Fill(vTrackMatrix);
countRecCandle++;
}
if(!prim) fRecSecTrackMatrix->Fill(vTrackMatrix);
if(fHistogramsOn)
FillHistograms(stack,iMc,AlidNdPtHelper::kRecTracks);
break;
}
}
}
}
if(countRecCandle>0) fRecCandleEventMatrix->Fill(vEventMatrix);
}
}
}
if(allChargedTracks) delete allChargedTracks; allChargedTracks = 0;
if(labelsAll) delete [] labelsAll; labelsAll = 0;
if(labelsAcc) delete [] labelsAcc; labelsAcc = 0;
if(labelsRec) delete [] labelsRec; labelsRec = 0;
if(!evtCuts->IsRecVertexRequired() && vtxESD != NULL) delete vtxESD;
}
void AlidNdPtAnalysis::FillHistograms(TObjArray *const allChargedTracks,Int_t *const labelsAll,Int_t multAll,Int_t *const labelsAcc,Int_t multAcc,Int_t *const labelsRec,Int_t multRec) {
if(!allChargedTracks) return;
if(!labelsAll) return;
if(!labelsAcc) return;
if(!labelsRec) return;
Int_t entries = allChargedTracks->GetEntries();
for(Int_t i=0; i<entries; ++i)
{
AliESDtrack *track = (AliESDtrack*)allChargedTracks->At(i);
if(!track) continue;
if(track->Charge() == 0) continue;
Int_t label = TMath::Abs(track->GetLabel());
for(Int_t iAll=0; iAll<multAll; ++iAll) {
if(label == labelsAll[iAll]) {
Double_t v1[2] = {track->Pt(), static_cast<Double_t>(multAll)};
fRecTrackMultHist1[AlidNdPtHelper::kAllTracks]->Fill(v1);
}
}
for(Int_t iAcc=0; iAcc<multAcc; ++iAcc) {
if(label == labelsAcc[iAcc]) {
Double_t v2[2] = {track->Pt(), static_cast<Double_t>(multAcc)};
fRecTrackMultHist1[AlidNdPtHelper::kAccTracks]->Fill(v2);
}
}
for(Int_t iRec=0; iRec<multRec; ++iRec) {
if(label == labelsRec[iRec]) {
Double_t v3[2] = {track->Pt(), static_cast<Double_t>(multRec)};
fRecTrackMultHist1[AlidNdPtHelper::kRecTracks]->Fill(v3);
}
}
}
}
void AlidNdPtAnalysis::FillHistograms(AliESDtrack *const esdTrack, AliStack *const stack, const Double_t zv, AlidNdPtHelper::TrackObject trackObj, Int_t multRecMult)
{
if(!esdTrack) return;
Float_t q = esdTrack->Charge();
if(TMath::Abs(q) < 0.001) return;
Float_t pt = esdTrack->Pt();
Float_t eta = esdTrack->Eta();
Float_t phi = esdTrack->Phi();
Float_t dca[2], bCov[3];
esdTrack->GetImpactParameters(dca,bCov);
Int_t nClust = esdTrack->GetTPCclusters(0);
Float_t chi2PerCluster = 0.;
if(nClust>0.) chi2PerCluster = esdTrack->GetTPCchi2()/Float_t(nClust);
Double_t values1[3] = {pt,eta,phi};
fRecTrackHist1[trackObj]->Fill(values1);
Double_t values[4] = {zv, pt,eta, static_cast<Double_t>(multRecMult)};
if(trackObj == AlidNdPtHelper::kRecTracks) {
fRecTrackHist->Fill(values);
}
if(!stack) return;
Int_t label = TMath::Abs(esdTrack->GetLabel());
if(label > stack->GetNtrack()) return;
TParticle* particle = stack->Particle(label);
if(!particle) return;
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());
if(!particle->GetPDG()) return;
Double_t gq = particle->GetPDG()->Charge()/3.0;
if(TMath::Abs(gq)<0.001) return;
Float_t gpt = particle->Pt();
Float_t geta = particle->Eta();
Double_t dpt=0;
if(gpt) dpt = (pt-gpt)/gpt;
Double_t deta = (eta-geta);
if(trackObj == AlidNdPtHelper::kRecTracks)
{
Double_t vTrackPtCorrelationMatrix[3]={pt,gpt,geta};
fTrackPtCorrelationMatrix->Fill(vTrackPtCorrelationMatrix);
Double_t vRecMCTrackHist1[4]={gpt,geta,dpt,deta};
fRecMCTrackHist1->Fill(vRecMCTrackHist1);
}
}
void AlidNdPtAnalysis::FillHistograms(AliStack *const stack, Int_t label, AlidNdPtHelper::TrackObject trackObj)
{
if(!stack) return;
if(label > stack->GetNtrack()) return;
TParticle* particle = stack->Particle(label);
if(!particle) return;
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());
Int_t mech = particle->GetUniqueID();
if(!particle->GetPDG()) return;
Double_t gq = particle->GetPDG()->Charge()/3.0;
if(TMath::Abs(gq) < 0.001) return;
Float_t gpt = particle->Pt();
Float_t geta = particle->Eta();
Float_t gphi = particle->Phi();
Bool_t prim = stack->IsPhysicalPrimary(label);
Int_t pid = AlidNdPtHelper::ConvertPdgToPid(particle);
Double_t vMCTrackHist1[3] = {gpt,geta,gphi};
fMCTrackHist1[trackObj]->Fill(vMCTrackHist1);
Double_t vMCPrimTrackHist1[5] = {gpt,geta,static_cast<Double_t>(pid),static_cast<Double_t>(mech),static_cast<Double_t>(motherPdg)};
Double_t vMCPrimTrackHist2[5] = {static_cast<Double_t>(TMath::Abs(particle->GetPdgCode())),static_cast<Double_t>(mech),static_cast<Double_t>(motherPdg)};
if(prim) {
fMCPrimTrackHist1[trackObj]->Fill(vMCPrimTrackHist1);
if(pid == 5) fMCPrimTrackHist2[trackObj]->Fill(vMCPrimTrackHist2);
}
else {
fMCSecTrackHist1[trackObj]->Fill(vMCPrimTrackHist1);
}
}
Long64_t AlidNdPtAnalysis::Merge(TCollection* const list)
{
if (!fIsInit) { Init(); }
if (!list)
return 0;
if (list->IsEmpty())
return 1;
TIterator* iter = list->MakeIterator();
TObject* obj = 0;
Int_t count=0;
while((obj = iter->Next()) != 0) {
AlidNdPtAnalysis* entry = dynamic_cast<AlidNdPtAnalysis*>(obj);
if (entry == 0) continue;
fRecEventHist->Add(entry->fRecEventHist);
fRecTrackHist->Add(entry->fRecTrackHist);
fEventMultCorrelationMatrix->Add(entry->fEventMultCorrelationMatrix);
fTrackPtCorrelationMatrix->Add(entry->fTrackPtCorrelationMatrix);
fGenEventMatrix->Add(entry->fGenEventMatrix);
fGenSDEventMatrix->Add(entry->fGenSDEventMatrix);
fGenDDEventMatrix->Add(entry->fGenDDEventMatrix);
fGenNDEventMatrix->Add(entry->fGenNDEventMatrix);
fGenNSDEventMatrix->Add(entry->fGenNSDEventMatrix);
fTriggerEventMatrix->Add(entry->fTriggerEventMatrix);
fTriggerSDEventMatrix->Add(entry->fTriggerSDEventMatrix);
fTriggerDDEventMatrix->Add(entry->fTriggerDDEventMatrix);
fTriggerNDEventMatrix->Add(entry->fTriggerNDEventMatrix);
fTriggerNSDEventMatrix->Add(entry->fTriggerNSDEventMatrix);
fRecEventMatrix->Add(entry->fRecEventMatrix);
fRecSDEventMatrix->Add(entry->fRecSDEventMatrix);
fRecDDEventMatrix->Add(entry->fRecDDEventMatrix);
fRecNDEventMatrix->Add(entry->fRecNDEventMatrix);
fRecNSDEventMatrix->Add(entry->fRecNSDEventMatrix);
fRecCandleEventMatrix->Add(entry->fRecCandleEventMatrix);
fGenTrackEventMatrix->Add(entry->fGenTrackEventMatrix);
fGenTrackSDEventMatrix->Add(entry->fGenTrackSDEventMatrix);
fGenTrackDDEventMatrix->Add(entry->fGenTrackDDEventMatrix);
fGenTrackNDEventMatrix->Add(entry->fGenTrackNDEventMatrix);
fGenTrackNSDEventMatrix->Add(entry->fGenTrackNSDEventMatrix);
fTriggerTrackEventMatrix->Add(entry->fTriggerTrackEventMatrix);
fTriggerTrackSDEventMatrix->Add(entry->fTriggerTrackSDEventMatrix);
fTriggerTrackDDEventMatrix->Add(entry->fTriggerTrackDDEventMatrix);
fTriggerTrackNDEventMatrix->Add(entry->fTriggerTrackNDEventMatrix);
fTriggerTrackNSDEventMatrix->Add(entry->fTriggerTrackNSDEventMatrix);
fRecTrackEventMatrix->Add(entry->fRecTrackEventMatrix);
fRecTrackSDEventMatrix->Add(entry->fRecTrackSDEventMatrix);
fRecTrackDDEventMatrix->Add(entry->fRecTrackDDEventMatrix);
fRecTrackNDEventMatrix->Add(entry->fRecTrackNDEventMatrix);
fRecTrackNSDEventMatrix->Add(entry->fRecTrackNSDEventMatrix);
fGenTrackMatrix->Add(entry->fGenTrackMatrix);
fGenPrimTrackMatrix->Add(entry->fGenPrimTrackMatrix);
fRecPrimTrackMatrix->Add(entry->fRecPrimTrackMatrix);
fRecTrackMatrix->Add(entry->fRecTrackMatrix);
fRecSecTrackMatrix->Add(entry->fRecSecTrackMatrix);
fRecMultTrackMatrix->Add(entry->fRecMultTrackMatrix);
fMCEventHist1->Add(entry->fMCEventHist1);
fRecEventHist1->Add(entry->fRecEventHist1);
fRecEventHist2->Add(entry->fRecEventHist2);
fRecMCEventHist1->Add(entry->fRecMCEventHist1);
fRecMCEventHist2->Add(entry->fRecMCEventHist2);
fRecMCEventHist3->Add(entry->fRecMCEventHist3);
for(Int_t i=0; i<AlidNdPtHelper::kCutSteps; i++) {
fMCTrackHist1[i]->Add(entry->fMCTrackHist1[i]);
fMCPrimTrackHist1[i]->Add(entry->fMCPrimTrackHist1[i]);
fMCPrimTrackHist2[i]->Add(entry->fMCPrimTrackHist2[i]);
fMCSecTrackHist1[i]->Add(entry->fMCSecTrackHist1[i]);
fRecTrackHist1[i]->Add(entry->fRecTrackHist1[i]);
fRecTrackMultHist1[i]->Add(entry->fRecTrackMultHist1[i]);
}
fRecMCTrackHist1->Add(entry->fRecMCTrackHist1);
fMCMultRecTrackHist1->Add(entry->fMCMultRecTrackHist1);
fRecTrackHist2->Add(entry->fRecTrackHist2);
fEventCount->Add(entry->fEventCount);
fMCPrimTrackHist->Add(entry->fMCPrimTrackHist);
count++;
}
return count;
}
Bool_t AlidNdPtAnalysis::IsBinZeroTrackSPDvtx(const AliESDEvent* esdEvent) {
if(!esdEvent) return kFALSE;
const AliESDVertex *vertex = esdEvent->GetPrimaryVertexTracks();
if(!vertex) return kTRUE;
if(vertex->GetNContributors() < 1) {
vertex = esdEvent->GetPrimaryVertexSPD();
if(!vertex) return kTRUE;
}
if(TMath::Abs(vertex->GetZ()) > 15.0) return kTRUE;
if( !AlidNdPtHelper::TestRecVertex(vertex, esdEvent->GetPrimaryVertexSPD(), AlidNdPtHelper::kTPCITSHybridTrackSPDvtx) )
return kTRUE;
return kFALSE;
}
Bool_t AlidNdPtAnalysis::IsBinZeroSPDvtx(const AliESDEvent* esdEvent) {
if(!esdEvent) return kFALSE;
const AliESDVertex* vertex = esdEvent->GetPrimaryVertexSPD();
if(!vertex) return kTRUE;
if(TMath::Abs(vertex->GetZ()) > 15.0) return kTRUE;
if( !AlidNdPtHelper::TestRecVertex(vertex, esdEvent->GetPrimaryVertexSPD(), AlidNdPtHelper::kTPCSPDvtx) ) return kTRUE;
return kFALSE;
}
void AlidNdPtAnalysis::Analyse()
{
if (!fIsInit) { Init(); }
TH1::AddDirectory(kFALSE);
TH1 *h=0, *h1=0, *h2=0, *h2c = 0;
THnSparse *hs=0;
TH2 *h2D=0;
char name[256];
TObjArray *aFolderObj = new TObjArray;
if(!aFolderObj) return;
h = fRecEventHist1->Projection(0);
h->SetName("Xv");
aFolderObj->Add(h);
h = fRecEventHist1->Projection(1);
h->SetName("Yv");
aFolderObj->Add(h);
h = fRecEventHist1->Projection(2);
h->SetName("Zv");
aFolderObj->Add(h);
h = fRecEventHist2->Projection(1);
h->SetName("multMB");
aFolderObj->Add(h);
h = fRecEventHist2->Projection(2);
h->SetName("multiplicity");
aFolderObj->Add(h);
h2D = fRecEventHist2->Projection(0,1);
h2D->SetName("Zv_vs_multiplicity_MB");
aFolderObj->Add(h2D);
h = fRecTrackHist1[0]->Projection(0);
h->Scale(1.,"width");
h->SetName("pt_all_ch");
aFolderObj->Add(h);
h = fRecTrackHist1[1]->Projection(0);
h->Scale(1.,"width");
h->SetName("pt_acc");
aFolderObj->Add(h);
h = fRecTrackHist1[2]->Projection(0);
h->Scale(1.,"width");
h->SetName("pt_rec");
aFolderObj->Add(h);
h = fRecTrackHist1[0]->Projection(1);
h->SetName("eta_all_ch");
aFolderObj->Add(h);
h = fRecTrackHist1[1]->Projection(1);
h->SetName("eta_acc");
aFolderObj->Add(h);
h = fRecTrackHist1[2]->Projection(1);
h->SetName("eta_rec");
aFolderObj->Add(h);
h = fRecTrackHist1[0]->Projection(2);
h->SetName("phi_all_ch");
aFolderObj->Add(h);
h = fRecTrackHist1[1]->Projection(2);
h->SetName("phi_acc");
aFolderObj->Add(h);
h = fRecTrackHist1[2]->Projection(2);
h->SetName("phi_rec");
aFolderObj->Add(h);
h2D = fRecTrackHist1[0]->Projection(1,0);
h2D->SetName("pt_eta_all_ch");
aFolderObj->Add(h2D);
h2D = fRecTrackHist1[1]->Projection(1,0);
h2D->SetName("pt_eta_acc");
aFolderObj->Add(h2D);
h2D = fRecTrackHist1[2]->Projection(1,0);
h2D->SetName("pt_eta_rec");
aFolderObj->Add(h2D);
h2D = fRecTrackHist1[0]->Projection(2,0);
h2D->SetName("pt_phi_all_ch");
aFolderObj->Add(h2D);
h2D = fRecTrackHist1[1]->Projection(2,0);
h2D->SetName("pt_phi_acc");
aFolderObj->Add(h2D);
h2D = fRecTrackHist1[2]->Projection(2,0);
h2D->SetName("pt_phi_rec");
aFolderObj->Add(h2D);
h2D = fRecTrackHist1[0]->Projection(2,1);
h2D->SetName("eta_phi_all_ch");
aFolderObj->Add(h2D);
h2D = fRecTrackHist1[1]->Projection(2,1);
h2D->SetName("eta_phi_acc");
aFolderObj->Add(h2D);
h2D = fRecTrackHist1[2]->Projection(2,1);
h2D->SetName("eta_phi_rec");
aFolderObj->Add(h2D);
if(fHistogramsOn) {
h2D = fRecTrackHist2->Projection(0,1);
h2D->SetName("nClust_chi2_rec");
aFolderObj->Add(h2D);
h2D = fRecTrackHist2->Projection(0,2);
h2D->SetName("nClust_pt_rec");
aFolderObj->Add(h2D);
h2D = fRecTrackHist2->Projection(0,3);
h2D->SetName("nClust_eta_rec");
aFolderObj->Add(h2D);
h2D = fRecTrackHist2->Projection(0,4);
h2D->SetName("nClust_phi_rec");
aFolderObj->Add(h2D);
h2D = fRecTrackHist2->Projection(1,2);
h2D->SetName("chi2_pt_rec");
aFolderObj->Add(h2D);
h2D = fRecTrackHist2->Projection(1,3);
h2D->SetName("chi2_eta_rec");
aFolderObj->Add(h2D);
h2D = fRecTrackHist2->Projection(1,4);
h2D->SetName("chi2_phi_rec");
aFolderObj->Add(h2D);
}
h = fRecEventHist2->Projection(0);
if( h->Integral() ) h->Scale(1./h->Integral());
h->SetName("zv_distribution_norm");
aFolderObj->Add(h);
if(IsUseMCInfo()) {
h2D = fRecMCEventHist2->Projection(0,2);
h2D->SetName("DeltaXv_vs_mult");
aFolderObj->Add(h2D);
h2D = fRecMCEventHist2->Projection(1,2);
h2D->SetName("DeltaZv_vs_mult");
aFolderObj->Add(h2D);
fTriggerEventMatrix->GetAxis(1)->SetRangeUser(0.,0.);
h = fTriggerEventMatrix->Projection(0);
h2D = fTriggerEventMatrix->Projection(0,1);
if(h2D->Integral()) h->Scale(1./h2D->Integral());
h1 = fRecEventMatrix->Projection(0);
h2D = fRecEventMatrix->Projection(0,1);
if(h2D->Integral()) h1->Scale(1./h2D->Integral());
h->Divide(h1);
h->SetName("zv_empty_events_norm");
aFolderObj->Add(h);
fTriggerEventMatrix->GetAxis(1)->SetRange(1,fTriggerEventMatrix->GetAxis(1)->GetNbins());
hs = (THnSparse*)fEventMultCorrelationMatrix->Clone("event_mult_correlation_matrix");
aFolderObj->Add(hs);
hs = (THnSparse*)fTrackPtCorrelationMatrix->Clone("track_pt_correlation_matrix");
aFolderObj->Add(hs);
h = AlidNdPtHelper::GenerateCorrMatrix(fTriggerEventMatrix->Projection(0),fGenEventMatrix->Projection(0),"zv_trig_INEL_eff_matrix");
aFolderObj->Add(h);
h = AlidNdPtHelper::GenerateCorrMatrix(fTriggerNSDEventMatrix->Projection(0),fGenNSDEventMatrix->Projection(0),"zv_trig_NSD_eff_matrix");
aFolderObj->Add(h);
hs = AlidNdPtHelper::GenerateCorrMatrix(fGenNDEventMatrix,fTriggerEventMatrix,"zv_mult_trig_MBtoND_corr_matrix");
aFolderObj->Add(hs);
h = AlidNdPtHelper::GenerateCorrMatrix(fGenNDEventMatrix->Projection(0),fTriggerEventMatrix->Projection(0),"zv_trig_MBtoND_corr_matrix");
aFolderObj->Add(h);
h = AlidNdPtHelper::GenerateCorrMatrix(fGenNDEventMatrix->Projection(1),fTriggerEventMatrix->Projection(1),"mult_trig_MBtoND_corr_matrix");
aFolderObj->Add(h);
hs = AlidNdPtHelper::GenerateCorrMatrix(fGenNSDEventMatrix,fTriggerEventMatrix,"zv_mult_trig_MBtoNSD_corr_matrix");
aFolderObj->Add(hs);
h2D = AlidNdPtHelper::GenerateCorrMatrix(fGenNSDEventMatrix->Projection(0,1),fTriggerEventMatrix->Projection(0,1),"zv_mult_trig_MBtoNSD_corr_matrix_2D");
aFolderObj->Add(h2D);
h = AlidNdPtHelper::GenerateCorrMatrix(fGenNSDEventMatrix->Projection(0),fTriggerEventMatrix->Projection(0),"zv_trig_MBtoNSD_corr_matrix");
aFolderObj->Add(h);
h = AlidNdPtHelper::GenerateCorrMatrix(fGenNSDEventMatrix->Projection(1),fTriggerEventMatrix->Projection(1),"mult_trig_MBtoNSD_corr_matrix");
aFolderObj->Add(h);
hs = AlidNdPtHelper::GenerateCorrMatrix(fGenEventMatrix,fTriggerEventMatrix,"zv_mult_trig_MBtoInel_corr_matrix");
aFolderObj->Add(hs);
h2D = AlidNdPtHelper::GenerateCorrMatrix(fGenEventMatrix->Projection(0,1),fTriggerEventMatrix->Projection(0,1),"zv_mult_trig_MBtoInel_corr_matrix_2D");
aFolderObj->Add(h2D);
h = AlidNdPtHelper::GenerateCorrMatrix(fGenEventMatrix->Projection(0),fTriggerEventMatrix->Projection(0),"zv_trig_MBtoInel_corr_matrix");
aFolderObj->Add(h);
h = AlidNdPtHelper::GenerateCorrMatrix(fGenEventMatrix->Projection(1),fTriggerEventMatrix->Projection(1),"mult_trig_MBtoInel_corr_matrix");
aFolderObj->Add(h);
hs = AlidNdPtHelper::GenerateCorrMatrix(fTriggerEventMatrix,fRecEventMatrix,"zv_mult_event_corr_matrix");
aFolderObj->Add(hs);
h2D = AlidNdPtHelper::GenerateCorrMatrix(fTriggerEventMatrix->Projection(0,1),fRecEventMatrix->Projection(0,1),"zv_mult_event_corr_matrix_2D");
aFolderObj->Add(h2D);
h = AlidNdPtHelper::GenerateCorrMatrix(fTriggerEventMatrix->Projection(1),fRecEventMatrix->Projection(1),"mult_event_corr_matrix");
aFolderObj->Add(h);
h = AlidNdPtHelper::GenerateCorrMatrix(fTriggerEventMatrix->Projection(0),fRecEventMatrix->Projection(0),"zv_event_corr_matrix");
aFolderObj->Add(h);
hs = AlidNdPtHelper::GenerateCorrMatrix(fGenTrackNDEventMatrix,fTriggerTrackEventMatrix,"zv_pt_eta_track_trig_MBtoND_corr_matrix");
aFolderObj->Add(hs);
h2D = AlidNdPtHelper::GenerateCorrMatrix(fGenTrackNDEventMatrix->Projection(1,2),fTriggerTrackEventMatrix->Projection(1,2),"eta_pt_track_trig_MBtoND_corr_matrix");
aFolderObj->Add(h2D);
h2D = AlidNdPtHelper::GenerateCorrMatrix(fGenTrackNDEventMatrix->Projection(1,0),fTriggerTrackEventMatrix->Projection(1,0),"pt_zv_track_trig_MBtoND_corr_matrix");
aFolderObj->Add(h2D);
h2D = AlidNdPtHelper::GenerateCorrMatrix(fGenTrackNDEventMatrix->Projection(2,0),fTriggerTrackEventMatrix->Projection(2,0),"zv_eta_track_trig_MBtoND_corr_matrix");
aFolderObj->Add(h2D);
hs = AlidNdPtHelper::GenerateCorrMatrix(fGenTrackNSDEventMatrix,fTriggerTrackEventMatrix,"zv_pt_eta_track_trig_MBtoNSD_corr_matrix");
aFolderObj->Add(hs);
h2D = AlidNdPtHelper::GenerateCorrMatrix(fGenTrackNSDEventMatrix->Projection(1,2),fTriggerTrackEventMatrix->Projection(1,2),"eta_pt_track_trig_MBtoNSD_corr_matrix");
aFolderObj->Add(h2D);
h2D = AlidNdPtHelper::GenerateCorrMatrix(fGenTrackNSDEventMatrix->Projection(1,0),fTriggerTrackEventMatrix->Projection(1,0),"pt_zv_track_trig_MBtoNSD_corr_matrix");
aFolderObj->Add(h2D);
h2D = AlidNdPtHelper::GenerateCorrMatrix(fGenTrackNSDEventMatrix->Projection(2,0),fTriggerTrackEventMatrix->Projection(2,0),"zv_eta_track_trig_MBtoNSD_corr_matrix");
aFolderObj->Add(h2D);
hs = AlidNdPtHelper::GenerateCorrMatrix(fGenTrackEventMatrix,fTriggerTrackEventMatrix,"zv_pt_eta_track_trig_MBtoInel_corr_matrix");
aFolderObj->Add(hs);
h2D = AlidNdPtHelper::GenerateCorrMatrix(fGenTrackEventMatrix->Projection(1,2),fTriggerTrackEventMatrix->Projection(1,2),"eta_pt_track_trig_MBtoInel_corr_matrix");
aFolderObj->Add(h2D);
h2D = AlidNdPtHelper::GenerateCorrMatrix(fGenTrackEventMatrix->Projection(1,0),fTriggerTrackEventMatrix->Projection(1,0),"pt_zv_track_trig_MBtoInel_corr_matrix");
aFolderObj->Add(h2D);
h2D = AlidNdPtHelper::GenerateCorrMatrix(fGenTrackEventMatrix->Projection(2,0),fTriggerTrackEventMatrix->Projection(2,0),"zv_eta_track_trig_MBtoInel_corr_matrix");
aFolderObj->Add(h2D);
h = AlidNdPtHelper::GenerateCorrMatrix(fTriggerTrackEventMatrix->Projection(1),fGenTrackEventMatrix->Projection(1),"pt_track_trig_MBtoInel_eff_matrix");
aFolderObj->Add(h);
hs = AlidNdPtHelper::GenerateCorrMatrix(fTriggerTrackEventMatrix,fRecTrackEventMatrix,"zv_pt_eta_track_event_corr_matrix");
aFolderObj->Add(hs);
h2D = AlidNdPtHelper::GenerateCorrMatrix(fTriggerTrackEventMatrix->Projection(1,2),fRecTrackEventMatrix->Projection(1,2),"eta_pt_track_event_corr_matrix");
aFolderObj->Add(h2D);
h2D = AlidNdPtHelper::GenerateCorrMatrix(fTriggerTrackEventMatrix->Projection(1,0),fRecTrackEventMatrix->Projection(1,0),"pt_zv_track_event_corr_matrix");
aFolderObj->Add(h2D);
h2D = AlidNdPtHelper::GenerateCorrMatrix(fTriggerTrackEventMatrix->Projection(2,0),fRecTrackEventMatrix->Projection(2,0),"zv_eta_track_event_corr_matrix");
aFolderObj->Add(h2D);
h = AlidNdPtHelper::GenerateCorrMatrix(fRecTrackEventMatrix->Projection(1),fTriggerTrackEventMatrix->Projection(1),"pt_track_event_eff_matrix");
aFolderObj->Add(h);
hs = AlidNdPtHelper::GenerateCorrMatrix(fGenPrimTrackMatrix,fRecPrimTrackMatrix,"zv_pt_eta_track_corr_matrix");
aFolderObj->Add(hs);
h2D = AlidNdPtHelper::GenerateCorrMatrix(fGenPrimTrackMatrix->Projection(1,2),fRecPrimTrackMatrix->Projection(1,2),"eta_pt_track_corr_matrix");
aFolderObj->Add(h2D);
h2D = AlidNdPtHelper::GenerateCorrMatrix(fGenPrimTrackMatrix->Projection(1,0),fRecPrimTrackMatrix->Projection(1,0),"pt_zv_track_corr_matrix");
aFolderObj->Add(h2D);
h2D = AlidNdPtHelper::GenerateCorrMatrix(fGenPrimTrackMatrix->Projection(2,0),fRecPrimTrackMatrix->Projection(2,0),"zv_eta_track_corr_matrix");
aFolderObj->Add(h2D);
h = AlidNdPtHelper::GenerateCorrMatrix(fGenPrimTrackMatrix->Projection(0),fRecPrimTrackMatrix->Projection(0),"zv_track_corr_matrix");
aFolderObj->Add(h);
h = AlidNdPtHelper::GenerateCorrMatrix(fGenPrimTrackMatrix->Projection(1),fRecPrimTrackMatrix->Projection(1),"pt_track_corr_matrix");
aFolderObj->Add(h);
h = AlidNdPtHelper::GenerateCorrMatrix(fRecPrimTrackMatrix->Projection(1), fGenPrimTrackMatrix->Projection(1),"pt_track_eff_matrix");
aFolderObj->Add(h);
h = AlidNdPtHelper::GenerateCorrMatrix(fGenPrimTrackMatrix->Projection(2),fRecPrimTrackMatrix->Projection(2),"eta_track_corr_matrix");
aFolderObj->Add(h);
hs = AlidNdPtHelper::GenerateCorrMatrix(fRecSecTrackMatrix,fRecTrackMatrix,"zv_pt_eta_track_cont_matrix");
aFolderObj->Add(hs);
h2D = AlidNdPtHelper::GenerateCorrMatrix(fRecSecTrackMatrix->Projection(1,2),fRecTrackMatrix->Projection(1,2),"eta_pt_track_cont_matrix");
aFolderObj->Add(h2D);
h2D = AlidNdPtHelper::GenerateCorrMatrix(fRecSecTrackMatrix->Projection(1,0),fRecTrackMatrix->Projection(1,0),"pt_zv_track_cont_matrix");
aFolderObj->Add(h2D);
h2D = AlidNdPtHelper::GenerateCorrMatrix(fRecSecTrackMatrix->Projection(2,0),fRecTrackMatrix->Projection(2,0),"zv_eta_track_cont_matrix");
aFolderObj->Add(h2D);
h = AlidNdPtHelper::GenerateCorrMatrix(fRecSecTrackMatrix->Projection(0),fRecTrackMatrix->Projection(0),"zv_track_cont_matrix");
aFolderObj->Add(h);
h = AlidNdPtHelper::GenerateCorrMatrix(fRecSecTrackMatrix->Projection(1),fRecTrackMatrix->Projection(1),"pt_track_cont_matrix");
aFolderObj->Add(h);
h = AlidNdPtHelper::GenerateCorrMatrix(fRecSecTrackMatrix->Projection(2),fRecTrackMatrix->Projection(2),"eta_track_cont_matrix");
aFolderObj->Add(h);
hs = AlidNdPtHelper::GenerateCorrMatrix(fRecMultTrackMatrix,fRecTrackMatrix,"zv_pt_eta_mult_track_cont_matrix");
aFolderObj->Add(hs);
h2D = AlidNdPtHelper::GenerateCorrMatrix(fRecMultTrackMatrix->Projection(1,2),fRecTrackMatrix->Projection(1,2),"eta_pt_mult_track_cont_matrix");
aFolderObj->Add(h2D);
h2D = AlidNdPtHelper::GenerateCorrMatrix(fRecMultTrackMatrix->Projection(1,0),fRecTrackMatrix->Projection(1,0),"pt_zv_mult_track_cont_matrix");
aFolderObj->Add(h2D);
h2D = AlidNdPtHelper::GenerateCorrMatrix(fRecMultTrackMatrix->Projection(2,0),fRecTrackMatrix->Projection(2,0),"zv_eta_mult_track_cont_matrix");
aFolderObj->Add(h2D);
h = AlidNdPtHelper::GenerateCorrMatrix(fRecMultTrackMatrix->Projection(0),fRecTrackMatrix->Projection(0),"zv_mult_track_cont_matrix");
aFolderObj->Add(h);
h = AlidNdPtHelper::GenerateCorrMatrix(fRecMultTrackMatrix->Projection(1),fRecTrackMatrix->Projection(1),"pt_mult_track_cont_matrix");
aFolderObj->Add(h);
h = AlidNdPtHelper::GenerateCorrMatrix(fRecMultTrackMatrix->Projection(2),fRecTrackMatrix->Projection(2),"eta_mult_track_cont_matrix");
aFolderObj->Add(h);
if(fHistogramsOn) {
fMCPrimTrackHist1[1]->GetAxis(2)->SetRange(1,1);
fMCPrimTrackHist1[2]->GetAxis(2)->SetRange(1,1);
h1 = fMCPrimTrackHist1[1]->Projection(0);
h2 = fMCPrimTrackHist1[2]->Projection(0);
h2c = (TH1D *)h2->Clone();
h2c->Divide(h1);
h2c->SetName("eff_pt_electrons");
aFolderObj->Add(h2c);
fMCPrimTrackHist1[1]->GetAxis(2)->SetRange(2,2);
fMCPrimTrackHist1[2]->GetAxis(2)->SetRange(2,2);
h1 = fMCPrimTrackHist1[1]->Projection(0);
h2 = fMCPrimTrackHist1[2]->Projection(0);
h2c = (TH1D *)h2->Clone();
h2c->Divide(h1);
h2c->SetName("eff_pt_muons");
aFolderObj->Add(h2c);
fMCPrimTrackHist1[1]->GetAxis(2)->SetRange(3,3);
fMCPrimTrackHist1[2]->GetAxis(2)->SetRange(3,3);
h1 = fMCPrimTrackHist1[1]->Projection(0);
h2 = fMCPrimTrackHist1[2]->Projection(0);
h2c = (TH1D *)h2->Clone();
h2c->Divide(h1);
h2c->SetName("eff_pt_pions");
aFolderObj->Add(h2c);
fMCPrimTrackHist1[1]->GetAxis(2)->SetRange(4,4);
fMCPrimTrackHist1[2]->GetAxis(2)->SetRange(4,4);
h1 = fMCPrimTrackHist1[1]->Projection(0);
h2 = fMCPrimTrackHist1[2]->Projection(0);
h2c = (TH1D *)h2->Clone();
h2c->Divide(h1);
h2c->SetName("eff_pt_kaons");
aFolderObj->Add(h2c);
fMCPrimTrackHist1[1]->GetAxis(2)->SetRange(5,5);
fMCPrimTrackHist1[2]->GetAxis(2)->SetRange(5,5);
h1 = fMCPrimTrackHist1[1]->Projection(0);
h2 = fMCPrimTrackHist1[2]->Projection(0);
h2c = (TH1D *)h2->Clone();
h2c->Divide(h1);
h2c->SetName("eff_pt_protons");
aFolderObj->Add(h2c);
fMCPrimTrackHist1[1]->GetAxis(2)->SetRange(1,5);
fMCPrimTrackHist1[2]->GetAxis(2)->SetRange(1,5);
h1 = fMCPrimTrackHist1[1]->Projection(0);
h2 = fMCPrimTrackHist1[2]->Projection(0);
h2c = (TH1D *)h2->Clone();
h2c->Divide(h1);
h2c->SetName("eff_pt_selected");
aFolderObj->Add(h2c);
fMCPrimTrackHist1[1]->GetAxis(2)->SetRange(1,6);
fMCPrimTrackHist1[2]->GetAxis(2)->SetRange(1,6);
h1 = fMCPrimTrackHist1[1]->Projection(0);
h2 = fMCPrimTrackHist1[2]->Projection(0);
h2c = (TH1D *)h2->Clone();
h2c->Divide(h1);
h2c->SetName("eff_pt_all");
aFolderObj->Add(h2c);
fMCPrimTrackHist1[1]->GetAxis(1)->SetRange(1,fMCPrimTrackHist1[1]->GetAxis(1)->GetNbins());
fMCPrimTrackHist1[2]->GetAxis(1)->SetRange(1,fMCPrimTrackHist1[2]->GetAxis(1)->GetNbins());
TH1D *mcPtAccall = fMCTrackHist1[1]->Projection(0);
mcPtAccall->SetName("mc_pt_acc_all");
aFolderObj->Add(mcPtAccall);
TH1D *mcPtAccprim = fMCPrimTrackHist1[1]->Projection(0);
mcPtAccprim->SetName("mc_pt_acc_prim");
aFolderObj->Add(mcPtAccprim);
TH1D *mcPtRecall = fMCTrackHist1[2]->Projection(0);
mcPtRecall->SetName("mc_pt_rec_all");
aFolderObj->Add(mcPtRecall);
TH1D *mcPtRecprim = fMCPrimTrackHist1[2]->Projection(0);
mcPtRecprim->SetName("mc_pt_rec_prim");
aFolderObj->Add(mcPtRecprim);
TH1D *mcPtRecsec = fMCSecTrackHist1[2]->Projection(0);
mcPtRecsec->SetName("mc_pt_rec_sec");
aFolderObj->Add(mcPtRecsec);
for(Int_t i = 0; i<6; i++)
{
snprintf(name,256,"mc_pt_rec_prim_pid_%d",i);
fMCPrimTrackHist1[2]->GetAxis(2)->SetRange(i+1,i+1);
h = fMCPrimTrackHist1[2]->Projection(0);
h->SetName(name);
aFolderObj->Add(h);
snprintf(name,256,"mc_pt_rec_sec_pid_%d",i);
fMCSecTrackHist1[2]->GetAxis(2)->SetRange(i+1,i+1);
h = fMCSecTrackHist1[2]->Projection(0);
h->SetName(name);
aFolderObj->Add(h);
fMCSecTrackHist1[2]->GetAxis(2)->SetRange(i+1,i+1);
for(Int_t j=0; j<20; j++) {
if(j == 4) {
snprintf(name,256,"mc_pt_rec_sec_pid_%d_decay",i);
fMCSecTrackHist1[2]->GetAxis(3)->SetRange(j+1,j+1);
h = fMCSecTrackHist1[2]->Projection(0);
h->SetName(name);
aFolderObj->Add(h);
snprintf(name,256,"mc_eta_rec_sec_pid_%d_decay",i);
fMCSecTrackHist1[2]->GetAxis(3)->SetRange(j+1,j+1);
h = fMCSecTrackHist1[2]->Projection(1);
h->SetName(name);
aFolderObj->Add(h);
snprintf(name,256,"mc_mother_rec_sec_pid_%d_decay",i);
fMCSecTrackHist1[2]->GetAxis(3)->SetRange(j+1,j+1);
h = fMCSecTrackHist1[2]->Projection(4);
h->SetName(name);
aFolderObj->Add(h);
} else if (j == 5) {
snprintf(name,256,"mc_pt_rec_sec_pid_%d_conv",i);
fMCSecTrackHist1[2]->GetAxis(3)->SetRange(j+1,j+1);
h = fMCSecTrackHist1[2]->Projection(0);
h->SetName(name);
aFolderObj->Add(h);
snprintf(name,256,"mc_eta_rec_sec_pid_%d_conv",i);
fMCSecTrackHist1[2]->GetAxis(3)->SetRange(j+1,j+1);
h = fMCSecTrackHist1[2]->Projection(1);
h->SetName(name);
aFolderObj->Add(h);
snprintf(name,256,"mc_mother_rec_sec_pid_%d_conv",i);
fMCSecTrackHist1[2]->GetAxis(3)->SetRange(j+1,j+1);
h = fMCSecTrackHist1[2]->Projection(4);
h->SetName(name);
aFolderObj->Add(h);
} else if (j == 13) {
snprintf(name,256,"mc_pt_rec_sec_pid_%d_mat",i);
fMCSecTrackHist1[2]->GetAxis(3)->SetRange(j+1,j+1);
h = fMCSecTrackHist1[2]->Projection(0);
h->SetName(name);
aFolderObj->Add(h);
snprintf(name,256,"mc_eta_rec_sec_pid_%d_mat",i);
fMCSecTrackHist1[2]->GetAxis(3)->SetRange(j+1,j+1);
h = fMCSecTrackHist1[2]->Projection(1);
h->SetName(name);
aFolderObj->Add(h);
snprintf(name,256,"mc_eta_mother_rec_sec_pid_%d_mat",i);
fMCSecTrackHist1[2]->GetAxis(3)->SetRange(j+1,j+1);
h = fMCSecTrackHist1[2]->Projection(4,1);
h->SetName(name);
aFolderObj->Add(h);
snprintf(name,256,"mc_mother_rec_sec_pid_%d_mat",i);
fMCSecTrackHist1[2]->GetAxis(3)->SetRange(j+1,j+1);
h = fMCSecTrackHist1[2]->Projection(4);
h->SetName(name);
aFolderObj->Add(h);
snprintf(name,256,"mc_pt_mother_rec_sec_pid_%d_mat",i);
fMCSecTrackHist1[2]->GetAxis(3)->SetRange(j+1,j+1);
h = fMCSecTrackHist1[2]->Projection(4,0);
h->SetName(name);
aFolderObj->Add(h);
} else {
continue;
}
}
}
}
TH2F *h2F=0;
TCanvas * c = new TCanvas("resol","resol");
c->cd();
fRecMCTrackHist1->GetAxis(1)->SetRangeUser(-0.8,0.79);
h2F = (TH2F*)fRecMCTrackHist1->Projection(2,0);
h = AlidNdPtHelper::MakeResol(h2F,1,0,kTRUE,10);
h->SetXTitle("p_{tmc} (GeV/c)");
h->SetYTitle("(p_{t}-p_{tmc})/p_{tmc} resolution");
h->Draw();
h->SetName("pt_resolution_vs_mcpt");
aFolderObj->Add(h);
h2F = (TH2F*)fRecMCTrackHist1->Projection(2,0);
h = AlidNdPtHelper::MakeResol(h2F,1,1,kTRUE,10);
h->SetXTitle("p_{tmc} (GeV/c)");
h->SetYTitle("(p_{t}-p_{tmc})/p_{tmc} mean");
h->Draw();
h->SetName("dpt_mean_vs_mcpt");
aFolderObj->Add(h);
h2F = (TH2F*)fRecMCTrackHist1->Projection(3,0);
h = AlidNdPtHelper::MakeResol(h2F,1,0,kTRUE,10);
h->SetXTitle("p_{tmc} (GeV/c)");
h->SetYTitle("(#eta-#eta_{mc}) resolution");
h->Draw();
h->SetName("eta_resolution_vs_mcpt");
aFolderObj->Add(h);
h2F = (TH2F*)fRecMCTrackHist1->Projection(3,0);
h = AlidNdPtHelper::MakeResol(h2F,1,1,kTRUE,10);
h->SetXTitle("p_{tmc} (GeV/c)");
h->SetYTitle("(#eta-mc#eta) mean");
h->Draw();
h->SetName("deta_mean_vs_mcpt");
aFolderObj->Add(h);
fRecMCTrackHist1->GetAxis(1)->SetRange(1,fRecMCTrackHist1->GetAxis(1)->GetNbins());
h2F = (TH2F*)fRecMCTrackHist1->Projection(2,1);
h = AlidNdPtHelper::MakeResol(h2F,1,0,kTRUE,10);
h->SetXTitle("#eta_{mc}");
h->SetYTitle("(p_{t}-p_{tmc})/p_{tmc} resolution");
h->Draw();
h->SetName("pt_resolution_vs_mceta");
aFolderObj->Add(h);
h2F = (TH2F*)fRecMCTrackHist1->Projection(2,1);
h = AlidNdPtHelper::MakeResol(h2F,1,1,kTRUE,10);
h->SetXTitle("#eta_{mc}");
h->SetYTitle("(p_{t}-p_{tmc})/p_{tmc} mean");
h->Draw();
h->SetName("dpt_mean_vs_mceta");
aFolderObj->Add(h);
h2F = (TH2F*)fRecMCTrackHist1->Projection(3,1);
h = AlidNdPtHelper::MakeResol(h2F,1,0,kTRUE,10);
h->SetXTitle("#eta_{mc}");
h->SetYTitle("(#eta-#eta_{mc}) resolution");
h->Draw();
h->SetName("eta_resolution_vs_mceta");
aFolderObj->Add(h);
h2F = (TH2F*)fRecMCTrackHist1->Projection(3,1);
h = AlidNdPtHelper::MakeResol(h2F,1,1,kTRUE,10);
h->SetXTitle("#eta_{mc}");
h->SetYTitle("(#eta-mc#eta) mean");
h->Draw();
h->SetName("deta_mean_vs_mceta");
aFolderObj->Add(h);
fRecMCTrackHist1->GetAxis(0)->SetRange(1,fRecMCTrackHist1->GetAxis(0)->GetNbins());
}
fAnalysisFolder = ExportToFolder(aFolderObj);
if(aFolderObj) delete aFolderObj;
}
TFolder* AlidNdPtAnalysis::ExportToFolder(TObjArray * const array)
{
AlidNdPtAnalysis * comp=this;
TFolder *folder = comp->GetAnalysisFolder();
TString name, title;
TFolder *newFolder = 0;
Int_t i = 0;
Int_t size = array->GetSize();
if(folder) {
name = folder->GetName();
title = folder->GetTitle();
delete folder;
newFolder = CreateFolder(name.Data(),title.Data());
newFolder->SetOwner();
while(i < size) {
newFolder->Add(array->At(i));
i++;
}
}
return newFolder;
}
TFolder* AlidNdPtAnalysis::CreateFolder(TString name,TString title) {
TFolder *folder = 0;
folder = new TFolder(name.Data(),title.Data());
return folder;
}
AlidNdPtAnalysis.cxx:1000 AlidNdPtAnalysis.cxx:1001 AlidNdPtAnalysis.cxx:1002 AlidNdPtAnalysis.cxx:1003 AlidNdPtAnalysis.cxx:1004 AlidNdPtAnalysis.cxx:1005 AlidNdPtAnalysis.cxx:1006 AlidNdPtAnalysis.cxx:1007 AlidNdPtAnalysis.cxx:1008 AlidNdPtAnalysis.cxx:1009 AlidNdPtAnalysis.cxx:1010 AlidNdPtAnalysis.cxx:1011 AlidNdPtAnalysis.cxx:1012 AlidNdPtAnalysis.cxx:1013 AlidNdPtAnalysis.cxx:1014 AlidNdPtAnalysis.cxx:1015 AlidNdPtAnalysis.cxx:1016 AlidNdPtAnalysis.cxx:1017 AlidNdPtAnalysis.cxx:1018 AlidNdPtAnalysis.cxx:1019 AlidNdPtAnalysis.cxx:1020 AlidNdPtAnalysis.cxx:1021 AlidNdPtAnalysis.cxx:1022 AlidNdPtAnalysis.cxx:1023 AlidNdPtAnalysis.cxx:1024 AlidNdPtAnalysis.cxx:1025 AlidNdPtAnalysis.cxx:1026 AlidNdPtAnalysis.cxx:1027 AlidNdPtAnalysis.cxx:1028 AlidNdPtAnalysis.cxx:1029 AlidNdPtAnalysis.cxx:1030 AlidNdPtAnalysis.cxx:1031 AlidNdPtAnalysis.cxx:1032 AlidNdPtAnalysis.cxx:1033 AlidNdPtAnalysis.cxx:1034 AlidNdPtAnalysis.cxx:1035 AlidNdPtAnalysis.cxx:1036 AlidNdPtAnalysis.cxx:1037 AlidNdPtAnalysis.cxx:1038 AlidNdPtAnalysis.cxx:1039 AlidNdPtAnalysis.cxx:1040 AlidNdPtAnalysis.cxx:1041 AlidNdPtAnalysis.cxx:1042 AlidNdPtAnalysis.cxx:1043 AlidNdPtAnalysis.cxx:1044 AlidNdPtAnalysis.cxx:1045 AlidNdPtAnalysis.cxx:1046 AlidNdPtAnalysis.cxx:1047 AlidNdPtAnalysis.cxx:1048 AlidNdPtAnalysis.cxx:1049 AlidNdPtAnalysis.cxx:1050 AlidNdPtAnalysis.cxx:1051 AlidNdPtAnalysis.cxx:1052 AlidNdPtAnalysis.cxx:1053 AlidNdPtAnalysis.cxx:1054 AlidNdPtAnalysis.cxx:1055 AlidNdPtAnalysis.cxx:1056 AlidNdPtAnalysis.cxx:1057 AlidNdPtAnalysis.cxx:1058 AlidNdPtAnalysis.cxx:1059 AlidNdPtAnalysis.cxx:1060 AlidNdPtAnalysis.cxx:1061 AlidNdPtAnalysis.cxx:1062 AlidNdPtAnalysis.cxx:1063 AlidNdPtAnalysis.cxx:1064 AlidNdPtAnalysis.cxx:1065 AlidNdPtAnalysis.cxx:1066 AlidNdPtAnalysis.cxx:1067 AlidNdPtAnalysis.cxx:1068 AlidNdPtAnalysis.cxx:1069 AlidNdPtAnalysis.cxx:1070 AlidNdPtAnalysis.cxx:1071 AlidNdPtAnalysis.cxx:1072 AlidNdPtAnalysis.cxx:1073 AlidNdPtAnalysis.cxx:1074 AlidNdPtAnalysis.cxx:1075 AlidNdPtAnalysis.cxx:1076 AlidNdPtAnalysis.cxx:1077 AlidNdPtAnalysis.cxx:1078 AlidNdPtAnalysis.cxx:1079 AlidNdPtAnalysis.cxx:1080 AlidNdPtAnalysis.cxx:1081 AlidNdPtAnalysis.cxx:1082 AlidNdPtAnalysis.cxx:1083 AlidNdPtAnalysis.cxx:1084 AlidNdPtAnalysis.cxx:1085 AlidNdPtAnalysis.cxx:1086 AlidNdPtAnalysis.cxx:1087 AlidNdPtAnalysis.cxx:1088 AlidNdPtAnalysis.cxx:1089 AlidNdPtAnalysis.cxx:1090 AlidNdPtAnalysis.cxx:1091 AlidNdPtAnalysis.cxx:1092 AlidNdPtAnalysis.cxx:1093 AlidNdPtAnalysis.cxx:1094 AlidNdPtAnalysis.cxx:1095 AlidNdPtAnalysis.cxx:1096 AlidNdPtAnalysis.cxx:1097 AlidNdPtAnalysis.cxx:1098 AlidNdPtAnalysis.cxx:1099 AlidNdPtAnalysis.cxx:1100 AlidNdPtAnalysis.cxx:1101 AlidNdPtAnalysis.cxx:1102 AlidNdPtAnalysis.cxx:1103 AlidNdPtAnalysis.cxx:1104 AlidNdPtAnalysis.cxx:1105 AlidNdPtAnalysis.cxx:1106 AlidNdPtAnalysis.cxx:1107 AlidNdPtAnalysis.cxx:1108 AlidNdPtAnalysis.cxx:1109 AlidNdPtAnalysis.cxx:1110 AlidNdPtAnalysis.cxx:1111 AlidNdPtAnalysis.cxx:1112 AlidNdPtAnalysis.cxx:1113 AlidNdPtAnalysis.cxx:1114 AlidNdPtAnalysis.cxx:1115 AlidNdPtAnalysis.cxx:1116 AlidNdPtAnalysis.cxx:1117 AlidNdPtAnalysis.cxx:1118 AlidNdPtAnalysis.cxx:1119 AlidNdPtAnalysis.cxx:1120 AlidNdPtAnalysis.cxx:1121 AlidNdPtAnalysis.cxx:1122 AlidNdPtAnalysis.cxx:1123 AlidNdPtAnalysis.cxx:1124 AlidNdPtAnalysis.cxx:1125 AlidNdPtAnalysis.cxx:1126 AlidNdPtAnalysis.cxx:1127 AlidNdPtAnalysis.cxx:1128 AlidNdPtAnalysis.cxx:1129 AlidNdPtAnalysis.cxx:1130 AlidNdPtAnalysis.cxx:1131 AlidNdPtAnalysis.cxx:1132 AlidNdPtAnalysis.cxx:1133 AlidNdPtAnalysis.cxx:1134 AlidNdPtAnalysis.cxx:1135 AlidNdPtAnalysis.cxx:1136 AlidNdPtAnalysis.cxx:1137 AlidNdPtAnalysis.cxx:1138 AlidNdPtAnalysis.cxx:1139 AlidNdPtAnalysis.cxx:1140 AlidNdPtAnalysis.cxx:1141 AlidNdPtAnalysis.cxx:1142 AlidNdPtAnalysis.cxx:1143 AlidNdPtAnalysis.cxx:1144 AlidNdPtAnalysis.cxx:1145 AlidNdPtAnalysis.cxx:1146 AlidNdPtAnalysis.cxx:1147 AlidNdPtAnalysis.cxx:1148 AlidNdPtAnalysis.cxx:1149 AlidNdPtAnalysis.cxx:1150 AlidNdPtAnalysis.cxx:1151 AlidNdPtAnalysis.cxx:1152 AlidNdPtAnalysis.cxx:1153 AlidNdPtAnalysis.cxx:1154 AlidNdPtAnalysis.cxx:1155 AlidNdPtAnalysis.cxx:1156 AlidNdPtAnalysis.cxx:1157 AlidNdPtAnalysis.cxx:1158 AlidNdPtAnalysis.cxx:1159 AlidNdPtAnalysis.cxx:1160 AlidNdPtAnalysis.cxx:1161 AlidNdPtAnalysis.cxx:1162 AlidNdPtAnalysis.cxx:1163 AlidNdPtAnalysis.cxx:1164 AlidNdPtAnalysis.cxx:1165 AlidNdPtAnalysis.cxx:1166 AlidNdPtAnalysis.cxx:1167 AlidNdPtAnalysis.cxx:1168 AlidNdPtAnalysis.cxx:1169 AlidNdPtAnalysis.cxx:1170 AlidNdPtAnalysis.cxx:1171 AlidNdPtAnalysis.cxx:1172 AlidNdPtAnalysis.cxx:1173 AlidNdPtAnalysis.cxx:1174 AlidNdPtAnalysis.cxx:1175 AlidNdPtAnalysis.cxx:1176 AlidNdPtAnalysis.cxx:1177 AlidNdPtAnalysis.cxx:1178 AlidNdPtAnalysis.cxx:1179 AlidNdPtAnalysis.cxx:1180 AlidNdPtAnalysis.cxx:1181 AlidNdPtAnalysis.cxx:1182 AlidNdPtAnalysis.cxx:1183 AlidNdPtAnalysis.cxx:1184 AlidNdPtAnalysis.cxx:1185 AlidNdPtAnalysis.cxx:1186 AlidNdPtAnalysis.cxx:1187 AlidNdPtAnalysis.cxx:1188 AlidNdPtAnalysis.cxx:1189 AlidNdPtAnalysis.cxx:1190 AlidNdPtAnalysis.cxx:1191 AlidNdPtAnalysis.cxx:1192 AlidNdPtAnalysis.cxx:1193 AlidNdPtAnalysis.cxx:1194 AlidNdPtAnalysis.cxx:1195 AlidNdPtAnalysis.cxx:1196 AlidNdPtAnalysis.cxx:1197 AlidNdPtAnalysis.cxx:1198 AlidNdPtAnalysis.cxx:1199 AlidNdPtAnalysis.cxx:1200 AlidNdPtAnalysis.cxx:1201 AlidNdPtAnalysis.cxx:1202 AlidNdPtAnalysis.cxx:1203 AlidNdPtAnalysis.cxx:1204 AlidNdPtAnalysis.cxx:1205 AlidNdPtAnalysis.cxx:1206 AlidNdPtAnalysis.cxx:1207 AlidNdPtAnalysis.cxx:1208 AlidNdPtAnalysis.cxx:1209 AlidNdPtAnalysis.cxx:1210 AlidNdPtAnalysis.cxx:1211 AlidNdPtAnalysis.cxx:1212 AlidNdPtAnalysis.cxx:1213 AlidNdPtAnalysis.cxx:1214 AlidNdPtAnalysis.cxx:1215 AlidNdPtAnalysis.cxx:1216 AlidNdPtAnalysis.cxx:1217 AlidNdPtAnalysis.cxx:1218 AlidNdPtAnalysis.cxx:1219 AlidNdPtAnalysis.cxx:1220 AlidNdPtAnalysis.cxx:1221 AlidNdPtAnalysis.cxx:1222 AlidNdPtAnalysis.cxx:1223 AlidNdPtAnalysis.cxx:1224 AlidNdPtAnalysis.cxx:1225 AlidNdPtAnalysis.cxx:1226 AlidNdPtAnalysis.cxx:1227 AlidNdPtAnalysis.cxx:1228 AlidNdPtAnalysis.cxx:1229 AlidNdPtAnalysis.cxx:1230 AlidNdPtAnalysis.cxx:1231 AlidNdPtAnalysis.cxx:1232 AlidNdPtAnalysis.cxx:1233 AlidNdPtAnalysis.cxx:1234 AlidNdPtAnalysis.cxx:1235 AlidNdPtAnalysis.cxx:1236 AlidNdPtAnalysis.cxx:1237 AlidNdPtAnalysis.cxx:1238 AlidNdPtAnalysis.cxx:1239 AlidNdPtAnalysis.cxx:1240 AlidNdPtAnalysis.cxx:1241 AlidNdPtAnalysis.cxx:1242 AlidNdPtAnalysis.cxx:1243 AlidNdPtAnalysis.cxx:1244 AlidNdPtAnalysis.cxx:1245 AlidNdPtAnalysis.cxx:1246 AlidNdPtAnalysis.cxx:1247 AlidNdPtAnalysis.cxx:1248 AlidNdPtAnalysis.cxx:1249 AlidNdPtAnalysis.cxx:1250 AlidNdPtAnalysis.cxx:1251 AlidNdPtAnalysis.cxx:1252 AlidNdPtAnalysis.cxx:1253 AlidNdPtAnalysis.cxx:1254 AlidNdPtAnalysis.cxx:1255 AlidNdPtAnalysis.cxx:1256 AlidNdPtAnalysis.cxx:1257 AlidNdPtAnalysis.cxx:1258 AlidNdPtAnalysis.cxx:1259 AlidNdPtAnalysis.cxx:1260 AlidNdPtAnalysis.cxx:1261 AlidNdPtAnalysis.cxx:1262 AlidNdPtAnalysis.cxx:1263 AlidNdPtAnalysis.cxx:1264 AlidNdPtAnalysis.cxx:1265 AlidNdPtAnalysis.cxx:1266 AlidNdPtAnalysis.cxx:1267 AlidNdPtAnalysis.cxx:1268 AlidNdPtAnalysis.cxx:1269 AlidNdPtAnalysis.cxx:1270 AlidNdPtAnalysis.cxx:1271 AlidNdPtAnalysis.cxx:1272 AlidNdPtAnalysis.cxx:1273 AlidNdPtAnalysis.cxx:1274 AlidNdPtAnalysis.cxx:1275 AlidNdPtAnalysis.cxx:1276 AlidNdPtAnalysis.cxx:1277 AlidNdPtAnalysis.cxx:1278 AlidNdPtAnalysis.cxx:1279 AlidNdPtAnalysis.cxx:1280 AlidNdPtAnalysis.cxx:1281 AlidNdPtAnalysis.cxx:1282 AlidNdPtAnalysis.cxx:1283 AlidNdPtAnalysis.cxx:1284 AlidNdPtAnalysis.cxx:1285 AlidNdPtAnalysis.cxx:1286 AlidNdPtAnalysis.cxx:1287 AlidNdPtAnalysis.cxx:1288 AlidNdPtAnalysis.cxx:1289 AlidNdPtAnalysis.cxx:1290 AlidNdPtAnalysis.cxx:1291 AlidNdPtAnalysis.cxx:1292 AlidNdPtAnalysis.cxx:1293 AlidNdPtAnalysis.cxx:1294 AlidNdPtAnalysis.cxx:1295 AlidNdPtAnalysis.cxx:1296 AlidNdPtAnalysis.cxx:1297 AlidNdPtAnalysis.cxx:1298 AlidNdPtAnalysis.cxx:1299 AlidNdPtAnalysis.cxx:1300 AlidNdPtAnalysis.cxx:1301 AlidNdPtAnalysis.cxx:1302 AlidNdPtAnalysis.cxx:1303 AlidNdPtAnalysis.cxx:1304 AlidNdPtAnalysis.cxx:1305 AlidNdPtAnalysis.cxx:1306 AlidNdPtAnalysis.cxx:1307 AlidNdPtAnalysis.cxx:1308 AlidNdPtAnalysis.cxx:1309 AlidNdPtAnalysis.cxx:1310 AlidNdPtAnalysis.cxx:1311 AlidNdPtAnalysis.cxx:1312 AlidNdPtAnalysis.cxx:1313 AlidNdPtAnalysis.cxx:1314 AlidNdPtAnalysis.cxx:1315 AlidNdPtAnalysis.cxx:1316 AlidNdPtAnalysis.cxx:1317 AlidNdPtAnalysis.cxx:1318 AlidNdPtAnalysis.cxx:1319 AlidNdPtAnalysis.cxx:1320 AlidNdPtAnalysis.cxx:1321 AlidNdPtAnalysis.cxx:1322 AlidNdPtAnalysis.cxx:1323 AlidNdPtAnalysis.cxx:1324 AlidNdPtAnalysis.cxx:1325 AlidNdPtAnalysis.cxx:1326 AlidNdPtAnalysis.cxx:1327 AlidNdPtAnalysis.cxx:1328 AlidNdPtAnalysis.cxx:1329 AlidNdPtAnalysis.cxx:1330 AlidNdPtAnalysis.cxx:1331 AlidNdPtAnalysis.cxx:1332 AlidNdPtAnalysis.cxx:1333 AlidNdPtAnalysis.cxx:1334 AlidNdPtAnalysis.cxx:1335 AlidNdPtAnalysis.cxx:1336 AlidNdPtAnalysis.cxx:1337 AlidNdPtAnalysis.cxx:1338 AlidNdPtAnalysis.cxx:1339 AlidNdPtAnalysis.cxx:1340 AlidNdPtAnalysis.cxx:1341 AlidNdPtAnalysis.cxx:1342 AlidNdPtAnalysis.cxx:1343 AlidNdPtAnalysis.cxx:1344 AlidNdPtAnalysis.cxx:1345 AlidNdPtAnalysis.cxx:1346 AlidNdPtAnalysis.cxx:1347 AlidNdPtAnalysis.cxx:1348 AlidNdPtAnalysis.cxx:1349 AlidNdPtAnalysis.cxx:1350 AlidNdPtAnalysis.cxx:1351 AlidNdPtAnalysis.cxx:1352 AlidNdPtAnalysis.cxx:1353 AlidNdPtAnalysis.cxx:1354 AlidNdPtAnalysis.cxx:1355 AlidNdPtAnalysis.cxx:1356 AlidNdPtAnalysis.cxx:1357 AlidNdPtAnalysis.cxx:1358 AlidNdPtAnalysis.cxx:1359 AlidNdPtAnalysis.cxx:1360 AlidNdPtAnalysis.cxx:1361 AlidNdPtAnalysis.cxx:1362 AlidNdPtAnalysis.cxx:1363 AlidNdPtAnalysis.cxx:1364 AlidNdPtAnalysis.cxx:1365 AlidNdPtAnalysis.cxx:1366 AlidNdPtAnalysis.cxx:1367 AlidNdPtAnalysis.cxx:1368 AlidNdPtAnalysis.cxx:1369 AlidNdPtAnalysis.cxx:1370 AlidNdPtAnalysis.cxx:1371 AlidNdPtAnalysis.cxx:1372 AlidNdPtAnalysis.cxx:1373 AlidNdPtAnalysis.cxx:1374 AlidNdPtAnalysis.cxx:1375 AlidNdPtAnalysis.cxx:1376 AlidNdPtAnalysis.cxx:1377 AlidNdPtAnalysis.cxx:1378 AlidNdPtAnalysis.cxx:1379 AlidNdPtAnalysis.cxx:1380 AlidNdPtAnalysis.cxx:1381 AlidNdPtAnalysis.cxx:1382 AlidNdPtAnalysis.cxx:1383 AlidNdPtAnalysis.cxx:1384 AlidNdPtAnalysis.cxx:1385 AlidNdPtAnalysis.cxx:1386 AlidNdPtAnalysis.cxx:1387 AlidNdPtAnalysis.cxx:1388 AlidNdPtAnalysis.cxx:1389 AlidNdPtAnalysis.cxx:1390 AlidNdPtAnalysis.cxx:1391 AlidNdPtAnalysis.cxx:1392 AlidNdPtAnalysis.cxx:1393 AlidNdPtAnalysis.cxx:1394 AlidNdPtAnalysis.cxx:1395 AlidNdPtAnalysis.cxx:1396 AlidNdPtAnalysis.cxx:1397 AlidNdPtAnalysis.cxx:1398 AlidNdPtAnalysis.cxx:1399 AlidNdPtAnalysis.cxx:1400 AlidNdPtAnalysis.cxx:1401 AlidNdPtAnalysis.cxx:1402 AlidNdPtAnalysis.cxx:1403 AlidNdPtAnalysis.cxx:1404 AlidNdPtAnalysis.cxx:1405 AlidNdPtAnalysis.cxx:1406 AlidNdPtAnalysis.cxx:1407 AlidNdPtAnalysis.cxx:1408 AlidNdPtAnalysis.cxx:1409 AlidNdPtAnalysis.cxx:1410 AlidNdPtAnalysis.cxx:1411 AlidNdPtAnalysis.cxx:1412 AlidNdPtAnalysis.cxx:1413 AlidNdPtAnalysis.cxx:1414 AlidNdPtAnalysis.cxx:1415 AlidNdPtAnalysis.cxx:1416 AlidNdPtAnalysis.cxx:1417 AlidNdPtAnalysis.cxx:1418 AlidNdPtAnalysis.cxx:1419 AlidNdPtAnalysis.cxx:1420 AlidNdPtAnalysis.cxx:1421 AlidNdPtAnalysis.cxx:1422 AlidNdPtAnalysis.cxx:1423 AlidNdPtAnalysis.cxx:1424 AlidNdPtAnalysis.cxx:1425 AlidNdPtAnalysis.cxx:1426 AlidNdPtAnalysis.cxx:1427 AlidNdPtAnalysis.cxx:1428 AlidNdPtAnalysis.cxx:1429 AlidNdPtAnalysis.cxx:1430 AlidNdPtAnalysis.cxx:1431 AlidNdPtAnalysis.cxx:1432 AlidNdPtAnalysis.cxx:1433 AlidNdPtAnalysis.cxx:1434 AlidNdPtAnalysis.cxx:1435 AlidNdPtAnalysis.cxx:1436 AlidNdPtAnalysis.cxx:1437 AlidNdPtAnalysis.cxx:1438 AlidNdPtAnalysis.cxx:1439 AlidNdPtAnalysis.cxx:1440 AlidNdPtAnalysis.cxx:1441 AlidNdPtAnalysis.cxx:1442 AlidNdPtAnalysis.cxx:1443 AlidNdPtAnalysis.cxx:1444 AlidNdPtAnalysis.cxx:1445 AlidNdPtAnalysis.cxx:1446 AlidNdPtAnalysis.cxx:1447 AlidNdPtAnalysis.cxx:1448 AlidNdPtAnalysis.cxx:1449 AlidNdPtAnalysis.cxx:1450 AlidNdPtAnalysis.cxx:1451 AlidNdPtAnalysis.cxx:1452 AlidNdPtAnalysis.cxx:1453 AlidNdPtAnalysis.cxx:1454 AlidNdPtAnalysis.cxx:1455 AlidNdPtAnalysis.cxx:1456 AlidNdPtAnalysis.cxx:1457 AlidNdPtAnalysis.cxx:1458 AlidNdPtAnalysis.cxx:1459 AlidNdPtAnalysis.cxx:1460 AlidNdPtAnalysis.cxx:1461 AlidNdPtAnalysis.cxx:1462 AlidNdPtAnalysis.cxx:1463 AlidNdPtAnalysis.cxx:1464 AlidNdPtAnalysis.cxx:1465 AlidNdPtAnalysis.cxx:1466 AlidNdPtAnalysis.cxx:1467 AlidNdPtAnalysis.cxx:1468 AlidNdPtAnalysis.cxx:1469 AlidNdPtAnalysis.cxx:1470 AlidNdPtAnalysis.cxx:1471 AlidNdPtAnalysis.cxx:1472 AlidNdPtAnalysis.cxx:1473 AlidNdPtAnalysis.cxx:1474 AlidNdPtAnalysis.cxx:1475 AlidNdPtAnalysis.cxx:1476 AlidNdPtAnalysis.cxx:1477 AlidNdPtAnalysis.cxx:1478 AlidNdPtAnalysis.cxx:1479 AlidNdPtAnalysis.cxx:1480 AlidNdPtAnalysis.cxx:1481 AlidNdPtAnalysis.cxx:1482 AlidNdPtAnalysis.cxx:1483 AlidNdPtAnalysis.cxx:1484 AlidNdPtAnalysis.cxx:1485 AlidNdPtAnalysis.cxx:1486 AlidNdPtAnalysis.cxx:1487 AlidNdPtAnalysis.cxx:1488 AlidNdPtAnalysis.cxx:1489 AlidNdPtAnalysis.cxx:1490 AlidNdPtAnalysis.cxx:1491 AlidNdPtAnalysis.cxx:1492 AlidNdPtAnalysis.cxx:1493 AlidNdPtAnalysis.cxx:1494 AlidNdPtAnalysis.cxx:1495 AlidNdPtAnalysis.cxx:1496 AlidNdPtAnalysis.cxx:1497 AlidNdPtAnalysis.cxx:1498 AlidNdPtAnalysis.cxx:1499 AlidNdPtAnalysis.cxx:1500 AlidNdPtAnalysis.cxx:1501 AlidNdPtAnalysis.cxx:1502 AlidNdPtAnalysis.cxx:1503 AlidNdPtAnalysis.cxx:1504 AlidNdPtAnalysis.cxx:1505 AlidNdPtAnalysis.cxx:1506 AlidNdPtAnalysis.cxx:1507 AlidNdPtAnalysis.cxx:1508 AlidNdPtAnalysis.cxx:1509 AlidNdPtAnalysis.cxx:1510 AlidNdPtAnalysis.cxx:1511 AlidNdPtAnalysis.cxx:1512 AlidNdPtAnalysis.cxx:1513 AlidNdPtAnalysis.cxx:1514 AlidNdPtAnalysis.cxx:1515 AlidNdPtAnalysis.cxx:1516 AlidNdPtAnalysis.cxx:1517 AlidNdPtAnalysis.cxx:1518 AlidNdPtAnalysis.cxx:1519 AlidNdPtAnalysis.cxx:1520 AlidNdPtAnalysis.cxx:1521 AlidNdPtAnalysis.cxx:1522 AlidNdPtAnalysis.cxx:1523 AlidNdPtAnalysis.cxx:1524 AlidNdPtAnalysis.cxx:1525 AlidNdPtAnalysis.cxx:1526 AlidNdPtAnalysis.cxx:1527 AlidNdPtAnalysis.cxx:1528 AlidNdPtAnalysis.cxx:1529 AlidNdPtAnalysis.cxx:1530 AlidNdPtAnalysis.cxx:1531 AlidNdPtAnalysis.cxx:1532 AlidNdPtAnalysis.cxx:1533 AlidNdPtAnalysis.cxx:1534 AlidNdPtAnalysis.cxx:1535 AlidNdPtAnalysis.cxx:1536 AlidNdPtAnalysis.cxx:1537 AlidNdPtAnalysis.cxx:1538 AlidNdPtAnalysis.cxx:1539 AlidNdPtAnalysis.cxx:1540 AlidNdPtAnalysis.cxx:1541 AlidNdPtAnalysis.cxx:1542 AlidNdPtAnalysis.cxx:1543 AlidNdPtAnalysis.cxx:1544 AlidNdPtAnalysis.cxx:1545 AlidNdPtAnalysis.cxx:1546 AlidNdPtAnalysis.cxx:1547 AlidNdPtAnalysis.cxx:1548 AlidNdPtAnalysis.cxx:1549 AlidNdPtAnalysis.cxx:1550 AlidNdPtAnalysis.cxx:1551 AlidNdPtAnalysis.cxx:1552 AlidNdPtAnalysis.cxx:1553 AlidNdPtAnalysis.cxx:1554 AlidNdPtAnalysis.cxx:1555 AlidNdPtAnalysis.cxx:1556 AlidNdPtAnalysis.cxx:1557 AlidNdPtAnalysis.cxx:1558 AlidNdPtAnalysis.cxx:1559 AlidNdPtAnalysis.cxx:1560 AlidNdPtAnalysis.cxx:1561 AlidNdPtAnalysis.cxx:1562 AlidNdPtAnalysis.cxx:1563 AlidNdPtAnalysis.cxx:1564 AlidNdPtAnalysis.cxx:1565 AlidNdPtAnalysis.cxx:1566 AlidNdPtAnalysis.cxx:1567 AlidNdPtAnalysis.cxx:1568 AlidNdPtAnalysis.cxx:1569 AlidNdPtAnalysis.cxx:1570 AlidNdPtAnalysis.cxx:1571 AlidNdPtAnalysis.cxx:1572 AlidNdPtAnalysis.cxx:1573 AlidNdPtAnalysis.cxx:1574 AlidNdPtAnalysis.cxx:1575 AlidNdPtAnalysis.cxx:1576 AlidNdPtAnalysis.cxx:1577 AlidNdPtAnalysis.cxx:1578 AlidNdPtAnalysis.cxx:1579 AlidNdPtAnalysis.cxx:1580 AlidNdPtAnalysis.cxx:1581 AlidNdPtAnalysis.cxx:1582 AlidNdPtAnalysis.cxx:1583 AlidNdPtAnalysis.cxx:1584 AlidNdPtAnalysis.cxx:1585 AlidNdPtAnalysis.cxx:1586 AlidNdPtAnalysis.cxx:1587 AlidNdPtAnalysis.cxx:1588 AlidNdPtAnalysis.cxx:1589 AlidNdPtAnalysis.cxx:1590 AlidNdPtAnalysis.cxx:1591 AlidNdPtAnalysis.cxx:1592 AlidNdPtAnalysis.cxx:1593 AlidNdPtAnalysis.cxx:1594 AlidNdPtAnalysis.cxx:1595 AlidNdPtAnalysis.cxx:1596 AlidNdPtAnalysis.cxx:1597 AlidNdPtAnalysis.cxx:1598 AlidNdPtAnalysis.cxx:1599 AlidNdPtAnalysis.cxx:1600 AlidNdPtAnalysis.cxx:1601 AlidNdPtAnalysis.cxx:1602 AlidNdPtAnalysis.cxx:1603 AlidNdPtAnalysis.cxx:1604 AlidNdPtAnalysis.cxx:1605 AlidNdPtAnalysis.cxx:1606 AlidNdPtAnalysis.cxx:1607 AlidNdPtAnalysis.cxx:1608 AlidNdPtAnalysis.cxx:1609 AlidNdPtAnalysis.cxx:1610 AlidNdPtAnalysis.cxx:1611 AlidNdPtAnalysis.cxx:1612 AlidNdPtAnalysis.cxx:1613 AlidNdPtAnalysis.cxx:1614 AlidNdPtAnalysis.cxx:1615 AlidNdPtAnalysis.cxx:1616 AlidNdPtAnalysis.cxx:1617 AlidNdPtAnalysis.cxx:1618 AlidNdPtAnalysis.cxx:1619 AlidNdPtAnalysis.cxx:1620 AlidNdPtAnalysis.cxx:1621 AlidNdPtAnalysis.cxx:1622 AlidNdPtAnalysis.cxx:1623 AlidNdPtAnalysis.cxx:1624 AlidNdPtAnalysis.cxx:1625 AlidNdPtAnalysis.cxx:1626 AlidNdPtAnalysis.cxx:1627 AlidNdPtAnalysis.cxx:1628 AlidNdPtAnalysis.cxx:1629 AlidNdPtAnalysis.cxx:1630 AlidNdPtAnalysis.cxx:1631 AlidNdPtAnalysis.cxx:1632 AlidNdPtAnalysis.cxx:1633 AlidNdPtAnalysis.cxx:1634 AlidNdPtAnalysis.cxx:1635 AlidNdPtAnalysis.cxx:1636 AlidNdPtAnalysis.cxx:1637 AlidNdPtAnalysis.cxx:1638 AlidNdPtAnalysis.cxx:1639 AlidNdPtAnalysis.cxx:1640 AlidNdPtAnalysis.cxx:1641 AlidNdPtAnalysis.cxx:1642 AlidNdPtAnalysis.cxx:1643 AlidNdPtAnalysis.cxx:1644 AlidNdPtAnalysis.cxx:1645 AlidNdPtAnalysis.cxx:1646 AlidNdPtAnalysis.cxx:1647 AlidNdPtAnalysis.cxx:1648 AlidNdPtAnalysis.cxx:1649 AlidNdPtAnalysis.cxx:1650 AlidNdPtAnalysis.cxx:1651 AlidNdPtAnalysis.cxx:1652 AlidNdPtAnalysis.cxx:1653 AlidNdPtAnalysis.cxx:1654 AlidNdPtAnalysis.cxx:1655 AlidNdPtAnalysis.cxx:1656 AlidNdPtAnalysis.cxx:1657 AlidNdPtAnalysis.cxx:1658 AlidNdPtAnalysis.cxx:1659 AlidNdPtAnalysis.cxx:1660 AlidNdPtAnalysis.cxx:1661 AlidNdPtAnalysis.cxx:1662 AlidNdPtAnalysis.cxx:1663 AlidNdPtAnalysis.cxx:1664 AlidNdPtAnalysis.cxx:1665 AlidNdPtAnalysis.cxx:1666 AlidNdPtAnalysis.cxx:1667 AlidNdPtAnalysis.cxx:1668 AlidNdPtAnalysis.cxx:1669 AlidNdPtAnalysis.cxx:1670 AlidNdPtAnalysis.cxx:1671 AlidNdPtAnalysis.cxx:1672 AlidNdPtAnalysis.cxx:1673 AlidNdPtAnalysis.cxx:1674 AlidNdPtAnalysis.cxx:1675 AlidNdPtAnalysis.cxx:1676 AlidNdPtAnalysis.cxx:1677 AlidNdPtAnalysis.cxx:1678 AlidNdPtAnalysis.cxx:1679 AlidNdPtAnalysis.cxx:1680 AlidNdPtAnalysis.cxx:1681 AlidNdPtAnalysis.cxx:1682 AlidNdPtAnalysis.cxx:1683 AlidNdPtAnalysis.cxx:1684 AlidNdPtAnalysis.cxx:1685 AlidNdPtAnalysis.cxx:1686 AlidNdPtAnalysis.cxx:1687 AlidNdPtAnalysis.cxx:1688 AlidNdPtAnalysis.cxx:1689 AlidNdPtAnalysis.cxx:1690 AlidNdPtAnalysis.cxx:1691 AlidNdPtAnalysis.cxx:1692 AlidNdPtAnalysis.cxx:1693 AlidNdPtAnalysis.cxx:1694 AlidNdPtAnalysis.cxx:1695 AlidNdPtAnalysis.cxx:1696 AlidNdPtAnalysis.cxx:1697 AlidNdPtAnalysis.cxx:1698 AlidNdPtAnalysis.cxx:1699 AlidNdPtAnalysis.cxx:1700 AlidNdPtAnalysis.cxx:1701 AlidNdPtAnalysis.cxx:1702 AlidNdPtAnalysis.cxx:1703 AlidNdPtAnalysis.cxx:1704 AlidNdPtAnalysis.cxx:1705 AlidNdPtAnalysis.cxx:1706 AlidNdPtAnalysis.cxx:1707 AlidNdPtAnalysis.cxx:1708 AlidNdPtAnalysis.cxx:1709 AlidNdPtAnalysis.cxx:1710 AlidNdPtAnalysis.cxx:1711 AlidNdPtAnalysis.cxx:1712 AlidNdPtAnalysis.cxx:1713 AlidNdPtAnalysis.cxx:1714 AlidNdPtAnalysis.cxx:1715 AlidNdPtAnalysis.cxx:1716 AlidNdPtAnalysis.cxx:1717 AlidNdPtAnalysis.cxx:1718 AlidNdPtAnalysis.cxx:1719 AlidNdPtAnalysis.cxx:1720 AlidNdPtAnalysis.cxx:1721 AlidNdPtAnalysis.cxx:1722 AlidNdPtAnalysis.cxx:1723 AlidNdPtAnalysis.cxx:1724 AlidNdPtAnalysis.cxx:1725 AlidNdPtAnalysis.cxx:1726 AlidNdPtAnalysis.cxx:1727 AlidNdPtAnalysis.cxx:1728 AlidNdPtAnalysis.cxx:1729 AlidNdPtAnalysis.cxx:1730 AlidNdPtAnalysis.cxx:1731 AlidNdPtAnalysis.cxx:1732 AlidNdPtAnalysis.cxx:1733 AlidNdPtAnalysis.cxx:1734 AlidNdPtAnalysis.cxx:1735 AlidNdPtAnalysis.cxx:1736 AlidNdPtAnalysis.cxx:1737 AlidNdPtAnalysis.cxx:1738 AlidNdPtAnalysis.cxx:1739 AlidNdPtAnalysis.cxx:1740 AlidNdPtAnalysis.cxx:1741 AlidNdPtAnalysis.cxx:1742 AlidNdPtAnalysis.cxx:1743 AlidNdPtAnalysis.cxx:1744 AlidNdPtAnalysis.cxx:1745 AlidNdPtAnalysis.cxx:1746 AlidNdPtAnalysis.cxx:1747 AlidNdPtAnalysis.cxx:1748 AlidNdPtAnalysis.cxx:1749 AlidNdPtAnalysis.cxx:1750 AlidNdPtAnalysis.cxx:1751 AlidNdPtAnalysis.cxx:1752 AlidNdPtAnalysis.cxx:1753 AlidNdPtAnalysis.cxx:1754 AlidNdPtAnalysis.cxx:1755 AlidNdPtAnalysis.cxx:1756 AlidNdPtAnalysis.cxx:1757 AlidNdPtAnalysis.cxx:1758 AlidNdPtAnalysis.cxx:1759 AlidNdPtAnalysis.cxx:1760 AlidNdPtAnalysis.cxx:1761 AlidNdPtAnalysis.cxx:1762 AlidNdPtAnalysis.cxx:1763 AlidNdPtAnalysis.cxx:1764 AlidNdPtAnalysis.cxx:1765 AlidNdPtAnalysis.cxx:1766 AlidNdPtAnalysis.cxx:1767 AlidNdPtAnalysis.cxx:1768 AlidNdPtAnalysis.cxx:1769 AlidNdPtAnalysis.cxx:1770 AlidNdPtAnalysis.cxx:1771 AlidNdPtAnalysis.cxx:1772 AlidNdPtAnalysis.cxx:1773 AlidNdPtAnalysis.cxx:1774 AlidNdPtAnalysis.cxx:1775 AlidNdPtAnalysis.cxx:1776 AlidNdPtAnalysis.cxx:1777 AlidNdPtAnalysis.cxx:1778 AlidNdPtAnalysis.cxx:1779 AlidNdPtAnalysis.cxx:1780 AlidNdPtAnalysis.cxx:1781 AlidNdPtAnalysis.cxx:1782 AlidNdPtAnalysis.cxx:1783 AlidNdPtAnalysis.cxx:1784 AlidNdPtAnalysis.cxx:1785 AlidNdPtAnalysis.cxx:1786 AlidNdPtAnalysis.cxx:1787 AlidNdPtAnalysis.cxx:1788 AlidNdPtAnalysis.cxx:1789 AlidNdPtAnalysis.cxx:1790 AlidNdPtAnalysis.cxx:1791 AlidNdPtAnalysis.cxx:1792 AlidNdPtAnalysis.cxx:1793 AlidNdPtAnalysis.cxx:1794 AlidNdPtAnalysis.cxx:1795 AlidNdPtAnalysis.cxx:1796 AlidNdPtAnalysis.cxx:1797 AlidNdPtAnalysis.cxx:1798 AlidNdPtAnalysis.cxx:1799 AlidNdPtAnalysis.cxx:1800 AlidNdPtAnalysis.cxx:1801 AlidNdPtAnalysis.cxx:1802 AlidNdPtAnalysis.cxx:1803 AlidNdPtAnalysis.cxx:1804 AlidNdPtAnalysis.cxx:1805 AlidNdPtAnalysis.cxx:1806 AlidNdPtAnalysis.cxx:1807 AlidNdPtAnalysis.cxx:1808 AlidNdPtAnalysis.cxx:1809 AlidNdPtAnalysis.cxx:1810 AlidNdPtAnalysis.cxx:1811 AlidNdPtAnalysis.cxx:1812 AlidNdPtAnalysis.cxx:1813 AlidNdPtAnalysis.cxx:1814 AlidNdPtAnalysis.cxx:1815 AlidNdPtAnalysis.cxx:1816 AlidNdPtAnalysis.cxx:1817 AlidNdPtAnalysis.cxx:1818 AlidNdPtAnalysis.cxx:1819 AlidNdPtAnalysis.cxx:1820 AlidNdPtAnalysis.cxx:1821 AlidNdPtAnalysis.cxx:1822 AlidNdPtAnalysis.cxx:1823 AlidNdPtAnalysis.cxx:1824 AlidNdPtAnalysis.cxx:1825 AlidNdPtAnalysis.cxx:1826 AlidNdPtAnalysis.cxx:1827 AlidNdPtAnalysis.cxx:1828 AlidNdPtAnalysis.cxx:1829 AlidNdPtAnalysis.cxx:1830 AlidNdPtAnalysis.cxx:1831 AlidNdPtAnalysis.cxx:1832 AlidNdPtAnalysis.cxx:1833 AlidNdPtAnalysis.cxx:1834 AlidNdPtAnalysis.cxx:1835 AlidNdPtAnalysis.cxx:1836 AlidNdPtAnalysis.cxx:1837 AlidNdPtAnalysis.cxx:1838 AlidNdPtAnalysis.cxx:1839 AlidNdPtAnalysis.cxx:1840 AlidNdPtAnalysis.cxx:1841 AlidNdPtAnalysis.cxx:1842 AlidNdPtAnalysis.cxx:1843 AlidNdPtAnalysis.cxx:1844 AlidNdPtAnalysis.cxx:1845 AlidNdPtAnalysis.cxx:1846 AlidNdPtAnalysis.cxx:1847 AlidNdPtAnalysis.cxx:1848 AlidNdPtAnalysis.cxx:1849 AlidNdPtAnalysis.cxx:1850 AlidNdPtAnalysis.cxx:1851 AlidNdPtAnalysis.cxx:1852 AlidNdPtAnalysis.cxx:1853 AlidNdPtAnalysis.cxx:1854 AlidNdPtAnalysis.cxx:1855 AlidNdPtAnalysis.cxx:1856 AlidNdPtAnalysis.cxx:1857 AlidNdPtAnalysis.cxx:1858 AlidNdPtAnalysis.cxx:1859 AlidNdPtAnalysis.cxx:1860 AlidNdPtAnalysis.cxx:1861 AlidNdPtAnalysis.cxx:1862 AlidNdPtAnalysis.cxx:1863 AlidNdPtAnalysis.cxx:1864 AlidNdPtAnalysis.cxx:1865 AlidNdPtAnalysis.cxx:1866 AlidNdPtAnalysis.cxx:1867 AlidNdPtAnalysis.cxx:1868 AlidNdPtAnalysis.cxx:1869 AlidNdPtAnalysis.cxx:1870 AlidNdPtAnalysis.cxx:1871 AlidNdPtAnalysis.cxx:1872 AlidNdPtAnalysis.cxx:1873 AlidNdPtAnalysis.cxx:1874 AlidNdPtAnalysis.cxx:1875 AlidNdPtAnalysis.cxx:1876 AlidNdPtAnalysis.cxx:1877 AlidNdPtAnalysis.cxx:1878 AlidNdPtAnalysis.cxx:1879 AlidNdPtAnalysis.cxx:1880 AlidNdPtAnalysis.cxx:1881 AlidNdPtAnalysis.cxx:1882 AlidNdPtAnalysis.cxx:1883 AlidNdPtAnalysis.cxx:1884 AlidNdPtAnalysis.cxx:1885 AlidNdPtAnalysis.cxx:1886 AlidNdPtAnalysis.cxx:1887 AlidNdPtAnalysis.cxx:1888 AlidNdPtAnalysis.cxx:1889 AlidNdPtAnalysis.cxx:1890 AlidNdPtAnalysis.cxx:1891 AlidNdPtAnalysis.cxx:1892 AlidNdPtAnalysis.cxx:1893 AlidNdPtAnalysis.cxx:1894 AlidNdPtAnalysis.cxx:1895 AlidNdPtAnalysis.cxx:1896 AlidNdPtAnalysis.cxx:1897 AlidNdPtAnalysis.cxx:1898 AlidNdPtAnalysis.cxx:1899 AlidNdPtAnalysis.cxx:1900 AlidNdPtAnalysis.cxx:1901 AlidNdPtAnalysis.cxx:1902 AlidNdPtAnalysis.cxx:1903 AlidNdPtAnalysis.cxx:1904 AlidNdPtAnalysis.cxx:1905 AlidNdPtAnalysis.cxx:1906 AlidNdPtAnalysis.cxx:1907 AlidNdPtAnalysis.cxx:1908 AlidNdPtAnalysis.cxx:1909 AlidNdPtAnalysis.cxx:1910 AlidNdPtAnalysis.cxx:1911 AlidNdPtAnalysis.cxx:1912 AlidNdPtAnalysis.cxx:1913 AlidNdPtAnalysis.cxx:1914 AlidNdPtAnalysis.cxx:1915 AlidNdPtAnalysis.cxx:1916 AlidNdPtAnalysis.cxx:1917 AlidNdPtAnalysis.cxx:1918 AlidNdPtAnalysis.cxx:1919 AlidNdPtAnalysis.cxx:1920 AlidNdPtAnalysis.cxx:1921 AlidNdPtAnalysis.cxx:1922 AlidNdPtAnalysis.cxx:1923 AlidNdPtAnalysis.cxx:1924 AlidNdPtAnalysis.cxx:1925 AlidNdPtAnalysis.cxx:1926 AlidNdPtAnalysis.cxx:1927 AlidNdPtAnalysis.cxx:1928 AlidNdPtAnalysis.cxx:1929 AlidNdPtAnalysis.cxx:1930 AlidNdPtAnalysis.cxx:1931 AlidNdPtAnalysis.cxx:1932 AlidNdPtAnalysis.cxx:1933 AlidNdPtAnalysis.cxx:1934 AlidNdPtAnalysis.cxx:1935 AlidNdPtAnalysis.cxx:1936 AlidNdPtAnalysis.cxx:1937 AlidNdPtAnalysis.cxx:1938 AlidNdPtAnalysis.cxx:1939 AlidNdPtAnalysis.cxx:1940 AlidNdPtAnalysis.cxx:1941 AlidNdPtAnalysis.cxx:1942 AlidNdPtAnalysis.cxx:1943 AlidNdPtAnalysis.cxx:1944 AlidNdPtAnalysis.cxx:1945 AlidNdPtAnalysis.cxx:1946 AlidNdPtAnalysis.cxx:1947 AlidNdPtAnalysis.cxx:1948 AlidNdPtAnalysis.cxx:1949 AlidNdPtAnalysis.cxx:1950 AlidNdPtAnalysis.cxx:1951 AlidNdPtAnalysis.cxx:1952 AlidNdPtAnalysis.cxx:1953 AlidNdPtAnalysis.cxx:1954 AlidNdPtAnalysis.cxx:1955 AlidNdPtAnalysis.cxx:1956 AlidNdPtAnalysis.cxx:1957 AlidNdPtAnalysis.cxx:1958 AlidNdPtAnalysis.cxx:1959 AlidNdPtAnalysis.cxx:1960 AlidNdPtAnalysis.cxx:1961 AlidNdPtAnalysis.cxx:1962 AlidNdPtAnalysis.cxx:1963 AlidNdPtAnalysis.cxx:1964 AlidNdPtAnalysis.cxx:1965 AlidNdPtAnalysis.cxx:1966 AlidNdPtAnalysis.cxx:1967 AlidNdPtAnalysis.cxx:1968 AlidNdPtAnalysis.cxx:1969 AlidNdPtAnalysis.cxx:1970 AlidNdPtAnalysis.cxx:1971 AlidNdPtAnalysis.cxx:1972 AlidNdPtAnalysis.cxx:1973 AlidNdPtAnalysis.cxx:1974 AlidNdPtAnalysis.cxx:1975 AlidNdPtAnalysis.cxx:1976 AlidNdPtAnalysis.cxx:1977 AlidNdPtAnalysis.cxx:1978 AlidNdPtAnalysis.cxx:1979 AlidNdPtAnalysis.cxx:1980 AlidNdPtAnalysis.cxx:1981 AlidNdPtAnalysis.cxx:1982 AlidNdPtAnalysis.cxx:1983 AlidNdPtAnalysis.cxx:1984 AlidNdPtAnalysis.cxx:1985 AlidNdPtAnalysis.cxx:1986 AlidNdPtAnalysis.cxx:1987 AlidNdPtAnalysis.cxx:1988 AlidNdPtAnalysis.cxx:1989 AlidNdPtAnalysis.cxx:1990 AlidNdPtAnalysis.cxx:1991 AlidNdPtAnalysis.cxx:1992 AlidNdPtAnalysis.cxx:1993 AlidNdPtAnalysis.cxx:1994 AlidNdPtAnalysis.cxx:1995 AlidNdPtAnalysis.cxx:1996 AlidNdPtAnalysis.cxx:1997 AlidNdPtAnalysis.cxx:1998 AlidNdPtAnalysis.cxx:1999 AlidNdPtAnalysis.cxx:2000 AlidNdPtAnalysis.cxx:2001 AlidNdPtAnalysis.cxx:2002 AlidNdPtAnalysis.cxx:2003 AlidNdPtAnalysis.cxx:2004 AlidNdPtAnalysis.cxx:2005 AlidNdPtAnalysis.cxx:2006 AlidNdPtAnalysis.cxx:2007 AlidNdPtAnalysis.cxx:2008 AlidNdPtAnalysis.cxx:2009 AlidNdPtAnalysis.cxx:2010 AlidNdPtAnalysis.cxx:2011 AlidNdPtAnalysis.cxx:2012 AlidNdPtAnalysis.cxx:2013 AlidNdPtAnalysis.cxx:2014 AlidNdPtAnalysis.cxx:2015 AlidNdPtAnalysis.cxx:2016 AlidNdPtAnalysis.cxx:2017 AlidNdPtAnalysis.cxx:2018 AlidNdPtAnalysis.cxx:2019 AlidNdPtAnalysis.cxx:2020 AlidNdPtAnalysis.cxx:2021 AlidNdPtAnalysis.cxx:2022 AlidNdPtAnalysis.cxx:2023 AlidNdPtAnalysis.cxx:2024 AlidNdPtAnalysis.cxx:2025 AlidNdPtAnalysis.cxx:2026 AlidNdPtAnalysis.cxx:2027 AlidNdPtAnalysis.cxx:2028 AlidNdPtAnalysis.cxx:2029 AlidNdPtAnalysis.cxx:2030 AlidNdPtAnalysis.cxx:2031 AlidNdPtAnalysis.cxx:2032 AlidNdPtAnalysis.cxx:2033 AlidNdPtAnalysis.cxx:2034 AlidNdPtAnalysis.cxx:2035 AlidNdPtAnalysis.cxx:2036 AlidNdPtAnalysis.cxx:2037 AlidNdPtAnalysis.cxx:2038 AlidNdPtAnalysis.cxx:2039 AlidNdPtAnalysis.cxx:2040 AlidNdPtAnalysis.cxx:2041 AlidNdPtAnalysis.cxx:2042 AlidNdPtAnalysis.cxx:2043 AlidNdPtAnalysis.cxx:2044 AlidNdPtAnalysis.cxx:2045 AlidNdPtAnalysis.cxx:2046 AlidNdPtAnalysis.cxx:2047 AlidNdPtAnalysis.cxx:2048 AlidNdPtAnalysis.cxx:2049 AlidNdPtAnalysis.cxx:2050 AlidNdPtAnalysis.cxx:2051 AlidNdPtAnalysis.cxx:2052 AlidNdPtAnalysis.cxx:2053 AlidNdPtAnalysis.cxx:2054 AlidNdPtAnalysis.cxx:2055 AlidNdPtAnalysis.cxx:2056 AlidNdPtAnalysis.cxx:2057 AlidNdPtAnalysis.cxx:2058 AlidNdPtAnalysis.cxx:2059 AlidNdPtAnalysis.cxx:2060 AlidNdPtAnalysis.cxx:2061 AlidNdPtAnalysis.cxx:2062 AlidNdPtAnalysis.cxx:2063 AlidNdPtAnalysis.cxx:2064 AlidNdPtAnalysis.cxx:2065 AlidNdPtAnalysis.cxx:2066 AlidNdPtAnalysis.cxx:2067 AlidNdPtAnalysis.cxx:2068 AlidNdPtAnalysis.cxx:2069 AlidNdPtAnalysis.cxx:2070 AlidNdPtAnalysis.cxx:2071 AlidNdPtAnalysis.cxx:2072 AlidNdPtAnalysis.cxx:2073 AlidNdPtAnalysis.cxx:2074 AlidNdPtAnalysis.cxx:2075 AlidNdPtAnalysis.cxx:2076 AlidNdPtAnalysis.cxx:2077 AlidNdPtAnalysis.cxx:2078 AlidNdPtAnalysis.cxx:2079 AlidNdPtAnalysis.cxx:2080 AlidNdPtAnalysis.cxx:2081 AlidNdPtAnalysis.cxx:2082 AlidNdPtAnalysis.cxx:2083 AlidNdPtAnalysis.cxx:2084 AlidNdPtAnalysis.cxx:2085 AlidNdPtAnalysis.cxx:2086 AlidNdPtAnalysis.cxx:2087 AlidNdPtAnalysis.cxx:2088 AlidNdPtAnalysis.cxx:2089 AlidNdPtAnalysis.cxx:2090 AlidNdPtAnalysis.cxx:2091 AlidNdPtAnalysis.cxx:2092 AlidNdPtAnalysis.cxx:2093 AlidNdPtAnalysis.cxx:2094 AlidNdPtAnalysis.cxx:2095 AlidNdPtAnalysis.cxx:2096 AlidNdPtAnalysis.cxx:2097 AlidNdPtAnalysis.cxx:2098 AlidNdPtAnalysis.cxx:2099 AlidNdPtAnalysis.cxx:2100 AlidNdPtAnalysis.cxx:2101 AlidNdPtAnalysis.cxx:2102 AlidNdPtAnalysis.cxx:2103 AlidNdPtAnalysis.cxx:2104 AlidNdPtAnalysis.cxx:2105 AlidNdPtAnalysis.cxx:2106 AlidNdPtAnalysis.cxx:2107 AlidNdPtAnalysis.cxx:2108 AlidNdPtAnalysis.cxx:2109 AlidNdPtAnalysis.cxx:2110 AlidNdPtAnalysis.cxx:2111 AlidNdPtAnalysis.cxx:2112 AlidNdPtAnalysis.cxx:2113 AlidNdPtAnalysis.cxx:2114 AlidNdPtAnalysis.cxx:2115 AlidNdPtAnalysis.cxx:2116 AlidNdPtAnalysis.cxx:2117 AlidNdPtAnalysis.cxx:2118 AlidNdPtAnalysis.cxx:2119 AlidNdPtAnalysis.cxx:2120 AlidNdPtAnalysis.cxx:2121 AlidNdPtAnalysis.cxx:2122 AlidNdPtAnalysis.cxx:2123 AlidNdPtAnalysis.cxx:2124 AlidNdPtAnalysis.cxx:2125 AlidNdPtAnalysis.cxx:2126 AlidNdPtAnalysis.cxx:2127 AlidNdPtAnalysis.cxx:2128 AlidNdPtAnalysis.cxx:2129 AlidNdPtAnalysis.cxx:2130 AlidNdPtAnalysis.cxx:2131 AlidNdPtAnalysis.cxx:2132 AlidNdPtAnalysis.cxx:2133 AlidNdPtAnalysis.cxx:2134 AlidNdPtAnalysis.cxx:2135 AlidNdPtAnalysis.cxx:2136 AlidNdPtAnalysis.cxx:2137 AlidNdPtAnalysis.cxx:2138 AlidNdPtAnalysis.cxx:2139 AlidNdPtAnalysis.cxx:2140 AlidNdPtAnalysis.cxx:2141 AlidNdPtAnalysis.cxx:2142 AlidNdPtAnalysis.cxx:2143 AlidNdPtAnalysis.cxx:2144 AlidNdPtAnalysis.cxx:2145 AlidNdPtAnalysis.cxx:2146 AlidNdPtAnalysis.cxx:2147 AlidNdPtAnalysis.cxx:2148 AlidNdPtAnalysis.cxx:2149 AlidNdPtAnalysis.cxx:2150 AlidNdPtAnalysis.cxx:2151 AlidNdPtAnalysis.cxx:2152 AlidNdPtAnalysis.cxx:2153 AlidNdPtAnalysis.cxx:2154 AlidNdPtAnalysis.cxx:2155 AlidNdPtAnalysis.cxx:2156 AlidNdPtAnalysis.cxx:2157 AlidNdPtAnalysis.cxx:2158 AlidNdPtAnalysis.cxx:2159 AlidNdPtAnalysis.cxx:2160 AlidNdPtAnalysis.cxx:2161 AlidNdPtAnalysis.cxx:2162 AlidNdPtAnalysis.cxx:2163 AlidNdPtAnalysis.cxx:2164 AlidNdPtAnalysis.cxx:2165 AlidNdPtAnalysis.cxx:2166 AlidNdPtAnalysis.cxx:2167 AlidNdPtAnalysis.cxx:2168 AlidNdPtAnalysis.cxx:2169 AlidNdPtAnalysis.cxx:2170 AlidNdPtAnalysis.cxx:2171 AlidNdPtAnalysis.cxx:2172 AlidNdPtAnalysis.cxx:2173 AlidNdPtAnalysis.cxx:2174 AlidNdPtAnalysis.cxx:2175 AlidNdPtAnalysis.cxx:2176 AlidNdPtAnalysis.cxx:2177 AlidNdPtAnalysis.cxx:2178 AlidNdPtAnalysis.cxx:2179 AlidNdPtAnalysis.cxx:2180 AlidNdPtAnalysis.cxx:2181 AlidNdPtAnalysis.cxx:2182 AlidNdPtAnalysis.cxx:2183 AlidNdPtAnalysis.cxx:2184 AlidNdPtAnalysis.cxx:2185 AlidNdPtAnalysis.cxx:2186 AlidNdPtAnalysis.cxx:2187 AlidNdPtAnalysis.cxx:2188 AlidNdPtAnalysis.cxx:2189 AlidNdPtAnalysis.cxx:2190 AlidNdPtAnalysis.cxx:2191 AlidNdPtAnalysis.cxx:2192 AlidNdPtAnalysis.cxx:2193 AlidNdPtAnalysis.cxx:2194 AlidNdPtAnalysis.cxx:2195 AlidNdPtAnalysis.cxx:2196 AlidNdPtAnalysis.cxx:2197 AlidNdPtAnalysis.cxx:2198 AlidNdPtAnalysis.cxx:2199 AlidNdPtAnalysis.cxx:2200 AlidNdPtAnalysis.cxx:2201 AlidNdPtAnalysis.cxx:2202 AlidNdPtAnalysis.cxx:2203 AlidNdPtAnalysis.cxx:2204 AlidNdPtAnalysis.cxx:2205 AlidNdPtAnalysis.cxx:2206 AlidNdPtAnalysis.cxx:2207 AlidNdPtAnalysis.cxx:2208 AlidNdPtAnalysis.cxx:2209 AlidNdPtAnalysis.cxx:2210 AlidNdPtAnalysis.cxx:2211 AlidNdPtAnalysis.cxx:2212 AlidNdPtAnalysis.cxx:2213 AlidNdPtAnalysis.cxx:2214 AlidNdPtAnalysis.cxx:2215 AlidNdPtAnalysis.cxx:2216 AlidNdPtAnalysis.cxx:2217 AlidNdPtAnalysis.cxx:2218 AlidNdPtAnalysis.cxx:2219 AlidNdPtAnalysis.cxx:2220 AlidNdPtAnalysis.cxx:2221 AlidNdPtAnalysis.cxx:2222 AlidNdPtAnalysis.cxx:2223 AlidNdPtAnalysis.cxx:2224 AlidNdPtAnalysis.cxx:2225 AlidNdPtAnalysis.cxx:2226 AlidNdPtAnalysis.cxx:2227 AlidNdPtAnalysis.cxx:2228 AlidNdPtAnalysis.cxx:2229 AlidNdPtAnalysis.cxx:2230 AlidNdPtAnalysis.cxx:2231 AlidNdPtAnalysis.cxx:2232 AlidNdPtAnalysis.cxx:2233 AlidNdPtAnalysis.cxx:2234 AlidNdPtAnalysis.cxx:2235 AlidNdPtAnalysis.cxx:2236 AlidNdPtAnalysis.cxx:2237 AlidNdPtAnalysis.cxx:2238 AlidNdPtAnalysis.cxx:2239 AlidNdPtAnalysis.cxx:2240 AlidNdPtAnalysis.cxx:2241 AlidNdPtAnalysis.cxx:2242 AlidNdPtAnalysis.cxx:2243 AlidNdPtAnalysis.cxx:2244 AlidNdPtAnalysis.cxx:2245 AlidNdPtAnalysis.cxx:2246 AlidNdPtAnalysis.cxx:2247 AlidNdPtAnalysis.cxx:2248 AlidNdPtAnalysis.cxx:2249 AlidNdPtAnalysis.cxx:2250 AlidNdPtAnalysis.cxx:2251 AlidNdPtAnalysis.cxx:2252 AlidNdPtAnalysis.cxx:2253 AlidNdPtAnalysis.cxx:2254 AlidNdPtAnalysis.cxx:2255 AlidNdPtAnalysis.cxx:2256 AlidNdPtAnalysis.cxx:2257 AlidNdPtAnalysis.cxx:2258 AlidNdPtAnalysis.cxx:2259 AlidNdPtAnalysis.cxx:2260 AlidNdPtAnalysis.cxx:2261 AlidNdPtAnalysis.cxx:2262 AlidNdPtAnalysis.cxx:2263 AlidNdPtAnalysis.cxx:2264 AlidNdPtAnalysis.cxx:2265 AlidNdPtAnalysis.cxx:2266 AlidNdPtAnalysis.cxx:2267 AlidNdPtAnalysis.cxx:2268 AlidNdPtAnalysis.cxx:2269 AlidNdPtAnalysis.cxx:2270 AlidNdPtAnalysis.cxx:2271 AlidNdPtAnalysis.cxx:2272 AlidNdPtAnalysis.cxx:2273 AlidNdPtAnalysis.cxx:2274 AlidNdPtAnalysis.cxx:2275 AlidNdPtAnalysis.cxx:2276 AlidNdPtAnalysis.cxx:2277 AlidNdPtAnalysis.cxx:2278 AlidNdPtAnalysis.cxx:2279 AlidNdPtAnalysis.cxx:2280 AlidNdPtAnalysis.cxx:2281 AlidNdPtAnalysis.cxx:2282 AlidNdPtAnalysis.cxx:2283 AlidNdPtAnalysis.cxx:2284 AlidNdPtAnalysis.cxx:2285 AlidNdPtAnalysis.cxx:2286 AlidNdPtAnalysis.cxx:2287 AlidNdPtAnalysis.cxx:2288 AlidNdPtAnalysis.cxx:2289 AlidNdPtAnalysis.cxx:2290 AlidNdPtAnalysis.cxx:2291 AlidNdPtAnalysis.cxx:2292 AlidNdPtAnalysis.cxx:2293 AlidNdPtAnalysis.cxx:2294 AlidNdPtAnalysis.cxx:2295 AlidNdPtAnalysis.cxx:2296 AlidNdPtAnalysis.cxx:2297 AlidNdPtAnalysis.cxx:2298 AlidNdPtAnalysis.cxx:2299 AlidNdPtAnalysis.cxx:2300 AlidNdPtAnalysis.cxx:2301 AlidNdPtAnalysis.cxx:2302 AlidNdPtAnalysis.cxx:2303 AlidNdPtAnalysis.cxx:2304 AlidNdPtAnalysis.cxx:2305 AlidNdPtAnalysis.cxx:2306 AlidNdPtAnalysis.cxx:2307 AlidNdPtAnalysis.cxx:2308 AlidNdPtAnalysis.cxx:2309 AlidNdPtAnalysis.cxx:2310 AlidNdPtAnalysis.cxx:2311 AlidNdPtAnalysis.cxx:2312 AlidNdPtAnalysis.cxx:2313 AlidNdPtAnalysis.cxx:2314 AlidNdPtAnalysis.cxx:2315 AlidNdPtAnalysis.cxx:2316 AlidNdPtAnalysis.cxx:2317 AlidNdPtAnalysis.cxx:2318 AlidNdPtAnalysis.cxx:2319 AlidNdPtAnalysis.cxx:2320 AlidNdPtAnalysis.cxx:2321 AlidNdPtAnalysis.cxx:2322 AlidNdPtAnalysis.cxx:2323 AlidNdPtAnalysis.cxx:2324 AlidNdPtAnalysis.cxx:2325 AlidNdPtAnalysis.cxx:2326 AlidNdPtAnalysis.cxx:2327 AlidNdPtAnalysis.cxx:2328 AlidNdPtAnalysis.cxx:2329 AlidNdPtAnalysis.cxx:2330 AlidNdPtAnalysis.cxx:2331 AlidNdPtAnalysis.cxx:2332 AlidNdPtAnalysis.cxx:2333 AlidNdPtAnalysis.cxx:2334 AlidNdPtAnalysis.cxx:2335 AlidNdPtAnalysis.cxx:2336 AlidNdPtAnalysis.cxx:2337 AlidNdPtAnalysis.cxx:2338 AlidNdPtAnalysis.cxx:2339 AlidNdPtAnalysis.cxx:2340 AlidNdPtAnalysis.cxx:2341 AlidNdPtAnalysis.cxx:2342 AlidNdPtAnalysis.cxx:2343 AlidNdPtAnalysis.cxx:2344 AlidNdPtAnalysis.cxx:2345 AlidNdPtAnalysis.cxx:2346 AlidNdPtAnalysis.cxx:2347 AlidNdPtAnalysis.cxx:2348 AlidNdPtAnalysis.cxx:2349 AlidNdPtAnalysis.cxx:2350 AlidNdPtAnalysis.cxx:2351 AlidNdPtAnalysis.cxx:2352 AlidNdPtAnalysis.cxx:2353 AlidNdPtAnalysis.cxx:2354 AlidNdPtAnalysis.cxx:2355 AlidNdPtAnalysis.cxx:2356 AlidNdPtAnalysis.cxx:2357 AlidNdPtAnalysis.cxx:2358 AlidNdPtAnalysis.cxx:2359 AlidNdPtAnalysis.cxx:2360 AlidNdPtAnalysis.cxx:2361 AlidNdPtAnalysis.cxx:2362 AlidNdPtAnalysis.cxx:2363 AlidNdPtAnalysis.cxx:2364 AlidNdPtAnalysis.cxx:2365 AlidNdPtAnalysis.cxx:2366 AlidNdPtAnalysis.cxx:2367 AlidNdPtAnalysis.cxx:2368 AlidNdPtAnalysis.cxx:2369 AlidNdPtAnalysis.cxx:2370 AlidNdPtAnalysis.cxx:2371 AlidNdPtAnalysis.cxx:2372 AlidNdPtAnalysis.cxx:2373 AlidNdPtAnalysis.cxx:2374 AlidNdPtAnalysis.cxx:2375 AlidNdPtAnalysis.cxx:2376 AlidNdPtAnalysis.cxx:2377 AlidNdPtAnalysis.cxx:2378 AlidNdPtAnalysis.cxx:2379 AlidNdPtAnalysis.cxx:2380 AlidNdPtAnalysis.cxx:2381 AlidNdPtAnalysis.cxx:2382 AlidNdPtAnalysis.cxx:2383 AlidNdPtAnalysis.cxx:2384 AlidNdPtAnalysis.cxx:2385 AlidNdPtAnalysis.cxx:2386 AlidNdPtAnalysis.cxx:2387 AlidNdPtAnalysis.cxx:2388 AlidNdPtAnalysis.cxx:2389 AlidNdPtAnalysis.cxx:2390 AlidNdPtAnalysis.cxx:2391 AlidNdPtAnalysis.cxx:2392 AlidNdPtAnalysis.cxx:2393 AlidNdPtAnalysis.cxx:2394 AlidNdPtAnalysis.cxx:2395 AlidNdPtAnalysis.cxx:2396 AlidNdPtAnalysis.cxx:2397 AlidNdPtAnalysis.cxx:2398 AlidNdPtAnalysis.cxx:2399 AlidNdPtAnalysis.cxx:2400 AlidNdPtAnalysis.cxx:2401 AlidNdPtAnalysis.cxx:2402 AlidNdPtAnalysis.cxx:2403 AlidNdPtAnalysis.cxx:2404 AlidNdPtAnalysis.cxx:2405 AlidNdPtAnalysis.cxx:2406 AlidNdPtAnalysis.cxx:2407 AlidNdPtAnalysis.cxx:2408 AlidNdPtAnalysis.cxx:2409 AlidNdPtAnalysis.cxx:2410 AlidNdPtAnalysis.cxx:2411 AlidNdPtAnalysis.cxx:2412 AlidNdPtAnalysis.cxx:2413 AlidNdPtAnalysis.cxx:2414 AlidNdPtAnalysis.cxx:2415 AlidNdPtAnalysis.cxx:2416 AlidNdPtAnalysis.cxx:2417 AlidNdPtAnalysis.cxx:2418 AlidNdPtAnalysis.cxx:2419 AlidNdPtAnalysis.cxx:2420 AlidNdPtAnalysis.cxx:2421 AlidNdPtAnalysis.cxx:2422 AlidNdPtAnalysis.cxx:2423 AlidNdPtAnalysis.cxx:2424 AlidNdPtAnalysis.cxx:2425 AlidNdPtAnalysis.cxx:2426 AlidNdPtAnalysis.cxx:2427 AlidNdPtAnalysis.cxx:2428 AlidNdPtAnalysis.cxx:2429 AlidNdPtAnalysis.cxx:2430 AlidNdPtAnalysis.cxx:2431 AlidNdPtAnalysis.cxx:2432 AlidNdPtAnalysis.cxx:2433 AlidNdPtAnalysis.cxx:2434 AlidNdPtAnalysis.cxx:2435 AlidNdPtAnalysis.cxx:2436 AlidNdPtAnalysis.cxx:2437 AlidNdPtAnalysis.cxx:2438 AlidNdPtAnalysis.cxx:2439 AlidNdPtAnalysis.cxx:2440 AlidNdPtAnalysis.cxx:2441 AlidNdPtAnalysis.cxx:2442 AlidNdPtAnalysis.cxx:2443 AlidNdPtAnalysis.cxx:2444 AlidNdPtAnalysis.cxx:2445 AlidNdPtAnalysis.cxx:2446 AlidNdPtAnalysis.cxx:2447 AlidNdPtAnalysis.cxx:2448 AlidNdPtAnalysis.cxx:2449 AlidNdPtAnalysis.cxx:2450 AlidNdPtAnalysis.cxx:2451 AlidNdPtAnalysis.cxx:2452 AlidNdPtAnalysis.cxx:2453 AlidNdPtAnalysis.cxx:2454 AlidNdPtAnalysis.cxx:2455 AlidNdPtAnalysis.cxx:2456 AlidNdPtAnalysis.cxx:2457 AlidNdPtAnalysis.cxx:2458 AlidNdPtAnalysis.cxx:2459 AlidNdPtAnalysis.cxx:2460 AlidNdPtAnalysis.cxx:2461 AlidNdPtAnalysis.cxx:2462 AlidNdPtAnalysis.cxx:2463 AlidNdPtAnalysis.cxx:2464 AlidNdPtAnalysis.cxx:2465 AlidNdPtAnalysis.cxx:2466 AlidNdPtAnalysis.cxx:2467 AlidNdPtAnalysis.cxx:2468 AlidNdPtAnalysis.cxx:2469 AlidNdPtAnalysis.cxx:2470 AlidNdPtAnalysis.cxx:2471 AlidNdPtAnalysis.cxx:2472 AlidNdPtAnalysis.cxx:2473 AlidNdPtAnalysis.cxx:2474 AlidNdPtAnalysis.cxx:2475 AlidNdPtAnalysis.cxx:2476 AlidNdPtAnalysis.cxx:2477 AlidNdPtAnalysis.cxx:2478 AlidNdPtAnalysis.cxx:2479 AlidNdPtAnalysis.cxx:2480 AlidNdPtAnalysis.cxx:2481 AlidNdPtAnalysis.cxx:2482 AlidNdPtAnalysis.cxx:2483 AlidNdPtAnalysis.cxx:2484 AlidNdPtAnalysis.cxx:2485 AlidNdPtAnalysis.cxx:2486 AlidNdPtAnalysis.cxx:2487 AlidNdPtAnalysis.cxx:2488 AlidNdPtAnalysis.cxx:2489 AlidNdPtAnalysis.cxx:2490 AlidNdPtAnalysis.cxx:2491 AlidNdPtAnalysis.cxx:2492 AlidNdPtAnalysis.cxx:2493 AlidNdPtAnalysis.cxx:2494 AlidNdPtAnalysis.cxx:2495 AlidNdPtAnalysis.cxx:2496 AlidNdPtAnalysis.cxx:2497 AlidNdPtAnalysis.cxx:2498 AlidNdPtAnalysis.cxx:2499 AlidNdPtAnalysis.cxx:2500 AlidNdPtAnalysis.cxx:2501 AlidNdPtAnalysis.cxx:2502 AlidNdPtAnalysis.cxx:2503 AlidNdPtAnalysis.cxx:2504 AlidNdPtAnalysis.cxx:2505 AlidNdPtAnalysis.cxx:2506 AlidNdPtAnalysis.cxx:2507 AlidNdPtAnalysis.cxx:2508 AlidNdPtAnalysis.cxx:2509 AlidNdPtAnalysis.cxx:2510 AlidNdPtAnalysis.cxx:2511 AlidNdPtAnalysis.cxx:2512 AlidNdPtAnalysis.cxx:2513 AlidNdPtAnalysis.cxx:2514 AlidNdPtAnalysis.cxx:2515 AlidNdPtAnalysis.cxx:2516 AlidNdPtAnalysis.cxx:2517 AlidNdPtAnalysis.cxx:2518 AlidNdPtAnalysis.cxx:2519 AlidNdPtAnalysis.cxx:2520 AlidNdPtAnalysis.cxx:2521 AlidNdPtAnalysis.cxx:2522 AlidNdPtAnalysis.cxx:2523 AlidNdPtAnalysis.cxx:2524 AlidNdPtAnalysis.cxx:2525 AlidNdPtAnalysis.cxx:2526 AlidNdPtAnalysis.cxx:2527 AlidNdPtAnalysis.cxx:2528 AlidNdPtAnalysis.cxx:2529 AlidNdPtAnalysis.cxx:2530 AlidNdPtAnalysis.cxx:2531 AlidNdPtAnalysis.cxx:2532 AlidNdPtAnalysis.cxx:2533 AlidNdPtAnalysis.cxx:2534 AlidNdPtAnalysis.cxx:2535 AlidNdPtAnalysis.cxx:2536 AlidNdPtAnalysis.cxx:2537 AlidNdPtAnalysis.cxx:2538 AlidNdPtAnalysis.cxx:2539 AlidNdPtAnalysis.cxx:2540 AlidNdPtAnalysis.cxx:2541 AlidNdPtAnalysis.cxx:2542 AlidNdPtAnalysis.cxx:2543 AlidNdPtAnalysis.cxx:2544 AlidNdPtAnalysis.cxx:2545 AlidNdPtAnalysis.cxx:2546 AlidNdPtAnalysis.cxx:2547 AlidNdPtAnalysis.cxx:2548 AlidNdPtAnalysis.cxx:2549 AlidNdPtAnalysis.cxx:2550 AlidNdPtAnalysis.cxx:2551 AlidNdPtAnalysis.cxx:2552 AlidNdPtAnalysis.cxx:2553 AlidNdPtAnalysis.cxx:2554 AlidNdPtAnalysis.cxx:2555 AlidNdPtAnalysis.cxx:2556 AlidNdPtAnalysis.cxx:2557 AlidNdPtAnalysis.cxx:2558 AlidNdPtAnalysis.cxx:2559 AlidNdPtAnalysis.cxx:2560 AlidNdPtAnalysis.cxx:2561 AlidNdPtAnalysis.cxx:2562 AlidNdPtAnalysis.cxx:2563 AlidNdPtAnalysis.cxx:2564 AlidNdPtAnalysis.cxx:2565 AlidNdPtAnalysis.cxx:2566 AlidNdPtAnalysis.cxx:2567 AlidNdPtAnalysis.cxx:2568 AlidNdPtAnalysis.cxx:2569 AlidNdPtAnalysis.cxx:2570 AlidNdPtAnalysis.cxx:2571 AlidNdPtAnalysis.cxx:2572 AlidNdPtAnalysis.cxx:2573 AlidNdPtAnalysis.cxx:2574 AlidNdPtAnalysis.cxx:2575 AlidNdPtAnalysis.cxx:2576 AlidNdPtAnalysis.cxx:2577 AlidNdPtAnalysis.cxx:2578 AlidNdPtAnalysis.cxx:2579 AlidNdPtAnalysis.cxx:2580 AlidNdPtAnalysis.cxx:2581 AlidNdPtAnalysis.cxx:2582 AlidNdPtAnalysis.cxx:2583 AlidNdPtAnalysis.cxx:2584 AlidNdPtAnalysis.cxx:2585 AlidNdPtAnalysis.cxx:2586 AlidNdPtAnalysis.cxx:2587 AlidNdPtAnalysis.cxx:2588 AlidNdPtAnalysis.cxx:2589 AlidNdPtAnalysis.cxx:2590 AlidNdPtAnalysis.cxx:2591 AlidNdPtAnalysis.cxx:2592 AlidNdPtAnalysis.cxx:2593 AlidNdPtAnalysis.cxx:2594 AlidNdPtAnalysis.cxx:2595 AlidNdPtAnalysis.cxx:2596 AlidNdPtAnalysis.cxx:2597 AlidNdPtAnalysis.cxx:2598 AlidNdPtAnalysis.cxx:2599 AlidNdPtAnalysis.cxx:2600 AlidNdPtAnalysis.cxx:2601 AlidNdPtAnalysis.cxx:2602 AlidNdPtAnalysis.cxx:2603 AlidNdPtAnalysis.cxx:2604 AlidNdPtAnalysis.cxx:2605 AlidNdPtAnalysis.cxx:2606 AlidNdPtAnalysis.cxx:2607 AlidNdPtAnalysis.cxx:2608 AlidNdPtAnalysis.cxx:2609 AlidNdPtAnalysis.cxx:2610 AlidNdPtAnalysis.cxx:2611 AlidNdPtAnalysis.cxx:2612 AlidNdPtAnalysis.cxx:2613 AlidNdPtAnalysis.cxx:2614 AlidNdPtAnalysis.cxx:2615 AlidNdPtAnalysis.cxx:2616 AlidNdPtAnalysis.cxx:2617 AlidNdPtAnalysis.cxx:2618 AlidNdPtAnalysis.cxx:2619 AlidNdPtAnalysis.cxx:2620 AlidNdPtAnalysis.cxx:2621 AlidNdPtAnalysis.cxx:2622 AlidNdPtAnalysis.cxx:2623 AlidNdPtAnalysis.cxx:2624 AlidNdPtAnalysis.cxx:2625 AlidNdPtAnalysis.cxx:2626 AlidNdPtAnalysis.cxx:2627 AlidNdPtAnalysis.cxx:2628 AlidNdPtAnalysis.cxx:2629 AlidNdPtAnalysis.cxx:2630 AlidNdPtAnalysis.cxx:2631 AlidNdPtAnalysis.cxx:2632 AlidNdPtAnalysis.cxx:2633 AlidNdPtAnalysis.cxx:2634 AlidNdPtAnalysis.cxx:2635 AlidNdPtAnalysis.cxx:2636 AlidNdPtAnalysis.cxx:2637 AlidNdPtAnalysis.cxx:2638 AlidNdPtAnalysis.cxx:2639 AlidNdPtAnalysis.cxx:2640 AlidNdPtAnalysis.cxx:2641 AlidNdPtAnalysis.cxx:2642 AlidNdPtAnalysis.cxx:2643 AlidNdPtAnalysis.cxx:2644 AlidNdPtAnalysis.cxx:2645 AlidNdPtAnalysis.cxx:2646 AlidNdPtAnalysis.cxx:2647 AlidNdPtAnalysis.cxx:2648 AlidNdPtAnalysis.cxx:2649 AlidNdPtAnalysis.cxx:2650 AlidNdPtAnalysis.cxx:2651 AlidNdPtAnalysis.cxx:2652 AlidNdPtAnalysis.cxx:2653 AlidNdPtAnalysis.cxx:2654 AlidNdPtAnalysis.cxx:2655 AlidNdPtAnalysis.cxx:2656 AlidNdPtAnalysis.cxx:2657 AlidNdPtAnalysis.cxx:2658 AlidNdPtAnalysis.cxx:2659 AlidNdPtAnalysis.cxx:2660 AlidNdPtAnalysis.cxx:2661 AlidNdPtAnalysis.cxx:2662 AlidNdPtAnalysis.cxx:2663 AlidNdPtAnalysis.cxx:2664 AlidNdPtAnalysis.cxx:2665 AlidNdPtAnalysis.cxx:2666 AlidNdPtAnalysis.cxx:2667 AlidNdPtAnalysis.cxx:2668 AlidNdPtAnalysis.cxx:2669 AlidNdPtAnalysis.cxx:2670 AlidNdPtAnalysis.cxx:2671 AlidNdPtAnalysis.cxx:2672 AlidNdPtAnalysis.cxx:2673 AlidNdPtAnalysis.cxx:2674 AlidNdPtAnalysis.cxx:2675 AlidNdPtAnalysis.cxx:2676 AlidNdPtAnalysis.cxx:2677 AlidNdPtAnalysis.cxx:2678 AlidNdPtAnalysis.cxx:2679 AlidNdPtAnalysis.cxx:2680 AlidNdPtAnalysis.cxx:2681 AlidNdPtAnalysis.cxx:2682 AlidNdPtAnalysis.cxx:2683 AlidNdPtAnalysis.cxx:2684 AlidNdPtAnalysis.cxx:2685 AlidNdPtAnalysis.cxx:2686 AlidNdPtAnalysis.cxx:2687 AlidNdPtAnalysis.cxx:2688 AlidNdPtAnalysis.cxx:2689 AlidNdPtAnalysis.cxx:2690 AlidNdPtAnalysis.cxx:2691 AlidNdPtAnalysis.cxx:2692 AlidNdPtAnalysis.cxx:2693 AlidNdPtAnalysis.cxx:2694 AlidNdPtAnalysis.cxx:2695 AlidNdPtAnalysis.cxx:2696 AlidNdPtAnalysis.cxx:2697 AlidNdPtAnalysis.cxx:2698 AlidNdPtAnalysis.cxx:2699 AlidNdPtAnalysis.cxx:2700 AlidNdPtAnalysis.cxx:2701 AlidNdPtAnalysis.cxx:2702 AlidNdPtAnalysis.cxx:2703 AlidNdPtAnalysis.cxx:2704 AlidNdPtAnalysis.cxx:2705 AlidNdPtAnalysis.cxx:2706 AlidNdPtAnalysis.cxx:2707 AlidNdPtAnalysis.cxx:2708 AlidNdPtAnalysis.cxx:2709 AlidNdPtAnalysis.cxx:2710 AlidNdPtAnalysis.cxx:2711 AlidNdPtAnalysis.cxx:2712 AlidNdPtAnalysis.cxx:2713 AlidNdPtAnalysis.cxx:2714 AlidNdPtAnalysis.cxx:2715 AlidNdPtAnalysis.cxx:2716 AlidNdPtAnalysis.cxx:2717 AlidNdPtAnalysis.cxx:2718 AlidNdPtAnalysis.cxx:2719 AlidNdPtAnalysis.cxx:2720 AlidNdPtAnalysis.cxx:2721 AlidNdPtAnalysis.cxx:2722 AlidNdPtAnalysis.cxx:2723 AlidNdPtAnalysis.cxx:2724 AlidNdPtAnalysis.cxx:2725 AlidNdPtAnalysis.cxx:2726 AlidNdPtAnalysis.cxx:2727 AlidNdPtAnalysis.cxx:2728 AlidNdPtAnalysis.cxx:2729 AlidNdPtAnalysis.cxx:2730 AlidNdPtAnalysis.cxx:2731 AlidNdPtAnalysis.cxx:2732 AlidNdPtAnalysis.cxx:2733 AlidNdPtAnalysis.cxx:2734 AlidNdPtAnalysis.cxx:2735 AlidNdPtAnalysis.cxx:2736 AlidNdPtAnalysis.cxx:2737 AlidNdPtAnalysis.cxx:2738 AlidNdPtAnalysis.cxx:2739 AlidNdPtAnalysis.cxx:2740 AlidNdPtAnalysis.cxx:2741 AlidNdPtAnalysis.cxx:2742 AlidNdPtAnalysis.cxx:2743 AlidNdPtAnalysis.cxx:2744 AlidNdPtAnalysis.cxx:2745 AlidNdPtAnalysis.cxx:2746 AlidNdPtAnalysis.cxx:2747 AlidNdPtAnalysis.cxx:2748 AlidNdPtAnalysis.cxx:2749 AlidNdPtAnalysis.cxx:2750 AlidNdPtAnalysis.cxx:2751 AlidNdPtAnalysis.cxx:2752 AlidNdPtAnalysis.cxx:2753 AlidNdPtAnalysis.cxx:2754 AlidNdPtAnalysis.cxx:2755 AlidNdPtAnalysis.cxx:2756 AlidNdPtAnalysis.cxx:2757 AlidNdPtAnalysis.cxx:2758 AlidNdPtAnalysis.cxx:2759 AlidNdPtAnalysis.cxx:2760 AlidNdPtAnalysis.cxx:2761 AlidNdPtAnalysis.cxx:2762 AlidNdPtAnalysis.cxx:2763 AlidNdPtAnalysis.cxx:2764 AlidNdPtAnalysis.cxx:2765 AlidNdPtAnalysis.cxx:2766 AlidNdPtAnalysis.cxx:2767 AlidNdPtAnalysis.cxx:2768 AlidNdPtAnalysis.cxx:2769 AlidNdPtAnalysis.cxx:2770 AlidNdPtAnalysis.cxx:2771 AlidNdPtAnalysis.cxx:2772 AlidNdPtAnalysis.cxx:2773 AlidNdPtAnalysis.cxx:2774 AlidNdPtAnalysis.cxx:2775 AlidNdPtAnalysis.cxx:2776 AlidNdPtAnalysis.cxx:2777 AlidNdPtAnalysis.cxx:2778 AlidNdPtAnalysis.cxx:2779 AlidNdPtAnalysis.cxx:2780 AlidNdPtAnalysis.cxx:2781 AlidNdPtAnalysis.cxx:2782 AlidNdPtAnalysis.cxx:2783 AlidNdPtAnalysis.cxx:2784 AlidNdPtAnalysis.cxx:2785 AlidNdPtAnalysis.cxx:2786 AlidNdPtAnalysis.cxx:2787 AlidNdPtAnalysis.cxx:2788 AlidNdPtAnalysis.cxx:2789 AlidNdPtAnalysis.cxx:2790 AlidNdPtAnalysis.cxx:2791 AlidNdPtAnalysis.cxx:2792 AlidNdPtAnalysis.cxx:2793 AlidNdPtAnalysis.cxx:2794 AlidNdPtAnalysis.cxx:2795 AlidNdPtAnalysis.cxx:2796 AlidNdPtAnalysis.cxx:2797 AlidNdPtAnalysis.cxx:2798 AlidNdPtAnalysis.cxx:2799 AlidNdPtAnalysis.cxx:2800 AlidNdPtAnalysis.cxx:2801 AlidNdPtAnalysis.cxx:2802 AlidNdPtAnalysis.cxx:2803 AlidNdPtAnalysis.cxx:2804 AlidNdPtAnalysis.cxx:2805 AlidNdPtAnalysis.cxx:2806 AlidNdPtAnalysis.cxx:2807 AlidNdPtAnalysis.cxx:2808 AlidNdPtAnalysis.cxx:2809 AlidNdPtAnalysis.cxx:2810 AlidNdPtAnalysis.cxx:2811 AlidNdPtAnalysis.cxx:2812 AlidNdPtAnalysis.cxx:2813 AlidNdPtAnalysis.cxx:2814 AlidNdPtAnalysis.cxx:2815 AlidNdPtAnalysis.cxx:2816 AlidNdPtAnalysis.cxx:2817 AlidNdPtAnalysis.cxx:2818 AlidNdPtAnalysis.cxx:2819 AlidNdPtAnalysis.cxx:2820 AlidNdPtAnalysis.cxx:2821 AlidNdPtAnalysis.cxx:2822 AlidNdPtAnalysis.cxx:2823 AlidNdPtAnalysis.cxx:2824 AlidNdPtAnalysis.cxx:2825 AlidNdPtAnalysis.cxx:2826 AlidNdPtAnalysis.cxx:2827 AlidNdPtAnalysis.cxx:2828 AlidNdPtAnalysis.cxx:2829 AlidNdPtAnalysis.cxx:2830 AlidNdPtAnalysis.cxx:2831 AlidNdPtAnalysis.cxx:2832 AlidNdPtAnalysis.cxx:2833 AlidNdPtAnalysis.cxx:2834 AlidNdPtAnalysis.cxx:2835 AlidNdPtAnalysis.cxx:2836 AlidNdPtAnalysis.cxx:2837 AlidNdPtAnalysis.cxx:2838 AlidNdPtAnalysis.cxx:2839 AlidNdPtAnalysis.cxx:2840 AlidNdPtAnalysis.cxx:2841 AlidNdPtAnalysis.cxx:2842 AlidNdPtAnalysis.cxx:2843 AlidNdPtAnalysis.cxx:2844 AlidNdPtAnalysis.cxx:2845 AlidNdPtAnalysis.cxx:2846 AlidNdPtAnalysis.cxx:2847 AlidNdPtAnalysis.cxx:2848 AlidNdPtAnalysis.cxx:2849 AlidNdPtAnalysis.cxx:2850 AlidNdPtAnalysis.cxx:2851 AlidNdPtAnalysis.cxx:2852 AlidNdPtAnalysis.cxx:2853 AlidNdPtAnalysis.cxx:2854 AlidNdPtAnalysis.cxx:2855 AlidNdPtAnalysis.cxx:2856 AlidNdPtAnalysis.cxx:2857 AlidNdPtAnalysis.cxx:2858 AlidNdPtAnalysis.cxx:2859 AlidNdPtAnalysis.cxx:2860 AlidNdPtAnalysis.cxx:2861 AlidNdPtAnalysis.cxx:2862 AlidNdPtAnalysis.cxx:2863 AlidNdPtAnalysis.cxx:2864 AlidNdPtAnalysis.cxx:2865 AlidNdPtAnalysis.cxx:2866 AlidNdPtAnalysis.cxx:2867 AlidNdPtAnalysis.cxx:2868 AlidNdPtAnalysis.cxx:2869 AlidNdPtAnalysis.cxx:2870 AlidNdPtAnalysis.cxx:2871 AlidNdPtAnalysis.cxx:2872 AlidNdPtAnalysis.cxx:2873 AlidNdPtAnalysis.cxx:2874 AlidNdPtAnalysis.cxx:2875 AlidNdPtAnalysis.cxx:2876 AlidNdPtAnalysis.cxx:2877 AlidNdPtAnalysis.cxx:2878 AlidNdPtAnalysis.cxx:2879 AlidNdPtAnalysis.cxx:2880 AlidNdPtAnalysis.cxx:2881 AlidNdPtAnalysis.cxx:2882 AlidNdPtAnalysis.cxx:2883 AlidNdPtAnalysis.cxx:2884 AlidNdPtAnalysis.cxx:2885 AlidNdPtAnalysis.cxx:2886 AlidNdPtAnalysis.cxx:2887 AlidNdPtAnalysis.cxx:2888 AlidNdPtAnalysis.cxx:2889 AlidNdPtAnalysis.cxx:2890 AlidNdPtAnalysis.cxx:2891 AlidNdPtAnalysis.cxx:2892 AlidNdPtAnalysis.cxx:2893 AlidNdPtAnalysis.cxx:2894 AlidNdPtAnalysis.cxx:2895 AlidNdPtAnalysis.cxx:2896 AlidNdPtAnalysis.cxx:2897 AlidNdPtAnalysis.cxx:2898 AlidNdPtAnalysis.cxx:2899 AlidNdPtAnalysis.cxx:2900 AlidNdPtAnalysis.cxx:2901 AlidNdPtAnalysis.cxx:2902