#include "TChain.h"
#include "TList.h"
#include "TFile.h"
#include "TTree.h"
#include "TH1D.h"
#include "TH2F.h"
#include "THnSparse.h"
#include "TCanvas.h"
#include "TParticle.h"
#include <TDatabasePDG.h>
#include <TParticlePDG.h>
#include "AliAnalysisTaskSE.h"
#include "AliAnalysisManager.h"
#include "AliESDVertex.h"
#include "AliESDEvent.h"
#include "AliESDInputHandler.h"
#include "AliAODEvent.h"
#include "AliAODTrack.h"
#include "AliAODInputHandler.h"
#include "AliESD.h"
#include "AliESDEvent.h"
#include "AliAODEvent.h"
#include "AliStack.h"
#include "AliESDtrackCuts.h"
#include "AliAODMCHeader.h"
#include "AliAODMCParticle.h"
#include "AliMCEventHandler.h"
#include "AliMCEvent.h"
#include "AliMCParticle.h"
#include "AliStack.h"
#include "AliGenHijingEventHeader.h"
#include "AliGenEventHeader.h"
#include "AliPID.h"
#include "AliAODPid.h"
#include "AliPIDResponse.h"
#include "AliAODpidUtil.h"
#include "AliPIDCombined.h"
#include "AliHigherMomentsToyModel.h"
using std::cout;
using std::endl;
ClassImp(AliHigherMomentsToyModel)
AliHigherMomentsToyModel::AliHigherMomentsToyModel( const char *name )
: AliAnalysisTaskSE( name ),
fListOfHistosQA(0),
fListOfHistos(0),
fAOD(0),
fArrayMC(0),
fPIDResponse(0),
fParticleSpecies(AliPID::kProton),
fAnalysisType("AOD"),
fCentralityEstimator("V0M"),
fCentrality(0),
fVxMax(3.),
fVyMax(3.),
fVzMax(10.),
fDCAxy(2.4),
fDCAz(3.),
fPtLowerLimit(0.3),
fPtHigherLimit(1.5),
fEtaLowerLimit(-0.8),
fEtaHigherLimit(0.8),
fRapidityCut(0.5),
fNSigmaCut(3.),
fTPCNClus(80),
fChi2perNDF(4.),
fAODtrackCutBit(1024),
fUsePid(kFALSE),
fEventCounter(0),
fHistDCA(0),
fTPCSig(0),
fTPCSigA(0),
fTHnCentNplusNminusCh(0),
fTHnCentNplusNminusChTruth(0),
fTHnCentNplusNminus(0)
{
for ( Int_t i = 0; i < 13; i++) {
fHistQA[i] = NULL;
}
for ( Int_t i = 0; i < 5; i++ ){
fTHnCentNplusNminusPid[i] = NULL;
fTHnCentNplusNminusPidTruth[i] = NULL;
}
DefineOutput(1, TList::Class());
DefineOutput(2, TList::Class());
}
AliHigherMomentsToyModel::~AliHigherMomentsToyModel() {
if(fListOfHistosQA) delete fListOfHistosQA;
if(fListOfHistos) delete fListOfHistos;
}
void AliHigherMomentsToyModel::UserCreateOutputObjects() {
fListOfHistosQA = new TList();
fListOfHistosQA->SetOwner(kTRUE);
fListOfHistos = new TList();
fListOfHistos->SetOwner(kTRUE);
fEventCounter = new TH1D("fEventCounter","EventCounter", 10, 0.5,10.5);
fEventCounter->GetXaxis()->SetBinLabel(1,"Event Accesed");
fEventCounter->GetXaxis()->SetBinLabel(2,"Within 0-90% centrality");
fEventCounter->GetXaxis()->SetBinLabel(5,"Have a vertex");
fEventCounter->GetXaxis()->SetBinLabel(6,"After vertex Cut");
fEventCounter->GetXaxis()->SetBinLabel(7,"Event Analyzed");
fEventCounter->GetXaxis()->SetBinLabel(8,"Event Analysis finished");
fListOfHistosQA->Add(fEventCounter);
fHistQA[0] = new TH1D("fHistQAvx", "Histo Vx After Cut", 400, -4., 4.);
fHistQA[1] = new TH1D("fHistQAvy", "Histo Vy After Cut", 400, -4., 4.);
fHistQA[2] = new TH1D("fHistQAvz", "Histo Vz After Cut", 500, -25., 25.);
fHistQA[3] = new TH1D("fHistQAvxA", "Histo Vx All", 500, -5., 5.);
fHistQA[4] = new TH1D("fHistQAvyA", "Histo Vy All", 500, -5., 5.);
fHistQA[5] = new TH1D("fHistQAvzA", "Histo Vz All", 500, -25., 25.);
fHistQA[6] = new TH1D("fHistQADcaXyC", "Histo DCAxy after cut", 500, -5., 5.);
fHistQA[7] = new TH1D("fHistQADcaZC", "Histo DCAz after cut", 500, -5., 5.);
fHistQA[8] = new TH1D("fHistQAPt","p_{T} distribution",600,0,6);
fHistQA[9] = new TH1D("fHistQAEta","#eta distribution",240,-1.2,1.2);
fHistQA[10] = new TH1D("fHistQAPhi","#phi distribution",340,0,6.8);
fHistQA[11] = new TH1D("fHistQANCls","Number of TPC cluster",200,0,200);
fHistQA[12] = new TH1D("fHistQAChi2","Chi2 per NDF",100,0,10);
for(Int_t i = 0; i < 13; i++)
{
fListOfHistosQA->Add(fHistQA[i]);
}
fHistDCA = new TH2D("fHistDCA","DCAxy Vs DCAz", 500, -5., 5., 500, -5., 5.);
fTPCSig = new TH2D("fTPCSig","TPC signal",200, 0.0, 10. ,1000,0.,1000);
fTPCSig->SetMarkerColor(kRed);
fTPCSigA = new TH2D("fTPCSigA","TPC signal all ",200, 0.0, 10. ,1000,0.,1000);
fListOfHistosQA->Add(fHistDCA);
fListOfHistosQA->Add(fTPCSig);
fListOfHistosQA->Add(fTPCSigA);
const Int_t nDim = 3;
const Int_t nPid = 5;
TString Species[nPid] = {"Electron","Muon","Pion","Kaon","Proton"};
Int_t fBins[nPid][nDim];
Double_t fMin[nPid][nDim];
Double_t fMax[nPid][nDim];
for( Int_t i = 0; i < nPid; i++ ){
for( Int_t j = 0; j < nDim; j++ ){
fBins[i][j] = 0;
fMin[i][j] = 0.;
fMax[i][j] = 0.;
}
}
Int_t fBinsCh[nDim] = {100, 1500, 1500};
Double_t fMinCh[nDim] = { -0.5, -0.5, -0.5 };
Double_t fMaxCh[nDim] = { 99.5, 1499.5, 1499.5};
fTHnCentNplusNminusCh = new THnSparseD("fTHnCentNplusNminusCh","Cent-NplusChrg-NminusChrg", nDim, fBinsCh, fMinCh, fMaxCh);
fTHnCentNplusNminusCh->GetAxis(0)->SetTitle("Centrality");
fTHnCentNplusNminusCh->GetAxis(1)->SetTitle("Nplus");
fTHnCentNplusNminusCh->GetAxis(2)->SetTitle("Nminus");
fListOfHistos->Add(fTHnCentNplusNminusCh);
if( fAnalysisType == "MCAOD" ){
fTHnCentNplusNminusChTruth = new THnSparseD("fTHnCentNplusNminusChTruth","Cent-NplusChrg-NminusChrg", nDim, fBinsCh, fMinCh, fMaxCh);
fTHnCentNplusNminusChTruth->GetAxis(0)->SetTitle("Centrality");
fTHnCentNplusNminusChTruth->GetAxis(1)->SetTitle("Nplus");
fTHnCentNplusNminusChTruth->GetAxis(2)->SetTitle("Nminus");
fListOfHistos->Add(fTHnCentNplusNminusChTruth);
}
TString hname1, hname11;
TString htitle1, axisTitle1, axisTitle2;
if( fUsePid ){
Int_t gPid = (Int_t)fParticleSpecies;
if( gPid > 1 && gPid < 5 ){
fBins[2][0] = 100; fBins[2][1] = 1000; fBins[2][2] = 1000;
fMin[2][0] = -0.5; fMin[2][1] = -0.5; fMin[2][2] = -0.5;
fMax[2][0] = 99.5; fMax[2][1] = 999.5; fMax[2][2] = 999.5;
fBins[3][0] = 100; fBins[3][1] = 600; fBins[3][2] = 600;
fMin[3][0] = -0.5; fMin[3][1] = -0.5; fMin[3][2] = -0.5;
fMax[3][0] = 99.5; fMax[3][1] = 599.5; fMax[3][2] = 599.5;
fBins[4][0] = 100; fBins[4][1] = 400; fBins[4][2] = 400;
fMin[4][0] = -0.5; fMin[4][1] = -0.5; fMin[4][2] = -0.5;
fMax[4][0] = 99.5; fMax[4][1] = 399.5; fMax[4][2] = 399.5;
hname1 = "fCentNplusNminusPid"; hname1 +=gPid;
htitle1 = Species[gPid]; htitle1 +=" And Neg-"; htitle1 +=Species[gPid];
axisTitle1 = Species[gPid];
axisTitle2 ="Neg-"; axisTitle2 += Species[gPid];
fTHnCentNplusNminusPid[gPid] = new THnSparseD(hname1.Data(),htitle1.Data(), nDim, fBins[gPid], fMin[gPid],fMax[gPid]);
fTHnCentNplusNminusPid[gPid]->GetAxis(0)->SetTitle("Centrality");
fTHnCentNplusNminusPid[gPid]->GetAxis(1)->SetTitle(axisTitle1.Data());
fTHnCentNplusNminusPid[gPid]->GetAxis(1)->SetTitle(axisTitle2.Data());
fListOfHistos->Add(fTHnCentNplusNminusPid[gPid]);
if( fAnalysisType == "MCAOD" ){
hname11 = "fCentNplusNminusPidTruth"; hname11 +=gPid;
fTHnCentNplusNminusPidTruth[gPid] = new THnSparseD(hname11.Data(),htitle1.Data(), nDim, fBins[gPid], fMin[gPid],fMax[gPid]);
fTHnCentNplusNminusPidTruth[gPid]->GetAxis(0)->SetTitle("Centrality");
fTHnCentNplusNminusPidTruth[gPid]->GetAxis(1)->SetTitle(axisTitle1.Data());
fTHnCentNplusNminusPidTruth[gPid]->GetAxis(1)->SetTitle(axisTitle2.Data());
fListOfHistos->Add(fTHnCentNplusNminusPidTruth[gPid]);
}
}
else{
Int_t fBinsX[nDim] = {100, 1500, 1500};
Double_t fMinX[nDim] = { -0.5, -0.5, -0.5 };
Double_t fMaxX[nDim] = { 99.5, 1499.5, 1499.5};
fTHnCentNplusNminus = new THnSparseD("fTHnCentNplusNminus","Cent-NplusChrg-NminusChrg", nDim, fBinsX, fMinX, fMaxX);
fTHnCentNplusNminus->GetAxis(0)->SetTitle("Centrality");
fTHnCentNplusNminus->GetAxis(1)->SetTitle("UnwantedPlus");
fTHnCentNplusNminus->GetAxis(2)->SetTitle("UnwantedMinus");
fListOfHistos->Add(fTHnCentNplusNminus);
}
}
PostData(1, fListOfHistosQA);
PostData(2, fListOfHistos);
}
void AliHigherMomentsToyModel::UserExec( Option_t * ){
fEventCounter->Fill(1);
if(fAnalysisType == "AOD") {
doAODEvent();
}
else if(fAnalysisType == "MCAOD") {
doMCAODEvent();
}
else return;
fEventCounter->Fill(8);
PostData(1, fListOfHistosQA);
PostData(2, fListOfHistos);
}
void AliHigherMomentsToyModel::doAODEvent(){
Double_t nPlusCharge = 0.;
Double_t nMinusCharge = 0.;
Double_t nPartile = 0.;
Double_t nAntiParticle = 0.;
Int_t gPid = 0.;
AliAnalysisManager* manager = AliAnalysisManager::GetAnalysisManager();
if (!manager) {
cout<<"ERROR: Analysis manager not found."<<endl;
return;
}
AliAODInputHandler* inputHandler = dynamic_cast<AliAODInputHandler*> (manager->GetInputEventHandler());
if (!inputHandler) {
cout<<"ERROR: Input handler not found."<<endl;
return;
}
fAOD = dynamic_cast<AliAODEvent*>(InputEvent());
if (!fAOD)
{
cout<< "ERROR 01: AOD not found " <<endl;
return;
}
fPIDResponse =inputHandler->GetPIDResponse();
if (!fPIDResponse){
AliFatal("This Task needs the PID response attached to the inputHandler");
return;
}
AliAODHeader *aodHeader = dynamic_cast<AliAODHeader*>(fAOD->GetHeader());
if(!aodHeader) AliFatal("Not a standard AOD");
fCentrality = (Int_t)aodHeader->GetCentralityP()->GetCentralityPercentile(fCentralityEstimator.Data());
if(fCentrality < 0 || fCentrality >= 91) return;
fEventCounter->Fill(2);
if(!ProperVertex()) return;
Int_t nTracks = fAOD->GetNumberOfTracks();
TExMap *trackMap = new TExMap();
for(Int_t i = 0; i < nTracks; i++) {
AliAODTrack* aodTrack = dynamic_cast<AliAODTrack *>(fAOD->GetTrack(i));
if(!aodTrack) {
AliError(Form("ERROR: Could not retrieve AODtrack %d",i));
continue;
}
Double_t tpcSignalAll = aodTrack->GetTPCsignal();
fTPCSigA->Fill(aodTrack->GetTPCmomentum(),tpcSignalAll);
Int_t gID = aodTrack->GetID();
if (aodTrack->TestFilterBit(1)) trackMap->Add(gID, i);
}
AliAODTrack* newAodTrack;
for( Int_t j = 0; j < nTracks; j++ )
{
AliAODTrack* aodTrack1 = dynamic_cast<AliAODTrack *>(fAOD->GetTrack(j));
if(!aodTrack1) {
AliError(Form("ERROR: Could not retrieve AODtrack %d",j));
continue;
}
if(!aodTrack1->TestFilterBit(fAODtrackCutBit)) continue;
Int_t gID = aodTrack1->GetID();
newAodTrack = gID >= 0 ? aodTrack1 : dynamic_cast<AliAODTrack*>(fAOD->GetTrack(trackMap->GetValue(-1-gID)));
if(!newAodTrack) AliFatal("Not a standard AOD");
Float_t dxy = 0., dz = 0.;
dxy = aodTrack1->DCA();
dz = aodTrack1->ZAtDCA();
Double_t pt = aodTrack1->Pt();
Double_t eta = aodTrack1->Eta();
Double_t nclus = aodTrack1->GetTPCClusterInfo(2,1);
Double_t chi2ndf = aodTrack1->Chi2perNDF();
if( pt < fPtLowerLimit || pt > fPtHigherLimit ) continue;
if( eta < fEtaLowerLimit || eta > fEtaHigherLimit ) continue;
if( nclus < fTPCNClus ) continue;
if( chi2ndf > fChi2perNDF ) continue;
fHistQA[6]->Fill(dxy);
fHistQA[7]->Fill(dz);
fHistQA[8]->Fill(pt);
fHistQA[9]->Fill(eta);
fHistQA[10]->Fill(aodTrack1->Phi());
fHistQA[11]->Fill(nclus);
fHistQA[12]->Fill(chi2ndf);
fHistDCA->Fill(dxy,dz);
Short_t gCharge = aodTrack1->Charge();
if(gCharge > 0) nPlusCharge += 1.;
if(gCharge < 0) nMinusCharge += 1.;
if( fUsePid ) {
gPid = (Int_t)fParticleSpecies;
Double_t rap = newAodTrack->Y(AliPID::ParticleMass(fParticleSpecies));
Double_t tpcSignal = newAodTrack->GetTPCsignal();
if( fabs(rap ) > fRapidityCut ) continue;
fTPCSig->Fill(newAodTrack->GetTPCmomentum(),tpcSignal);
Float_t nsigmaTPCPID = -999.;
nsigmaTPCPID = TMath::Abs(fPIDResponse->NumberOfSigmasTPC(newAodTrack,fParticleSpecies));
if ( nsigmaTPCPID < fNSigmaCut ){
if (gCharge > 0) nPartile +=1.;
if( gCharge < 0 ) nAntiParticle +=1.;
}
}
}
Double_t fContainerCh[3] = { static_cast<Double_t>(fCentrality), nPlusCharge, nMinusCharge};
Double_t fContainerPid[3] = { static_cast<Double_t>(fCentrality), nPartile, nAntiParticle};
fTHnCentNplusNminusCh->Fill(fContainerCh);
if( fUsePid ){
gPid = (Int_t)fParticleSpecies;
fTHnCentNplusNminusPid[gPid]->Fill(fContainerPid);
}
fEventCounter->Fill(7);
return;
}
void AliHigherMomentsToyModel::doMCAODEvent(){
Double_t nPlusCharge = 0.;
Double_t nMinusCharge = 0.;
Double_t nPlusChargeTruth = 0.;
Double_t nMinusChargeTruth = 0.;
Double_t nPartile = 0.;
Double_t nAntiParticle = 0.;
Double_t nPartileTruth = 0.;
Double_t nAntiParticleTruth = 0.;
Int_t gPid = 0;
Int_t gPdgCode = AliPID::ParticleCode(fParticleSpecies);
AliAnalysisManager* manager = AliAnalysisManager::GetAnalysisManager();
if (!manager) {
cout<<"ERROR: Analysis manager not found."<<endl;
return;
}
AliAODInputHandler* inputHandler = dynamic_cast<AliAODInputHandler*> (manager->GetInputEventHandler());
if (!inputHandler) {
cout<<"ERROR: Input handler not found."<<endl;
return;
}
fAOD = dynamic_cast<AliAODEvent*>(InputEvent());
if (!fAOD)
{
cout<< "ERROR 01: AOD not found " <<endl;
return;
}
fPIDResponse =inputHandler->GetPIDResponse();
if (!fPIDResponse){
AliFatal("This Task needs the PID response attached to the inputHandler");
return;
}
fArrayMC = dynamic_cast<TClonesArray*>(fAOD->FindListObject(AliAODMCParticle::StdBranchName()));
if (!fArrayMC) {
AliFatal("Error: MC particles branch not found!\n");
return;
}
AliAODMCHeader *mcHdr=NULL;
mcHdr=(AliAODMCHeader*)fAOD->GetList()->FindObject(AliAODMCHeader::StdBranchName());
if(!mcHdr) {
printf("MC header branch not found!\n");
return;
}
AliAODHeader *aodHeader = dynamic_cast<AliAODHeader*>(fAOD->GetHeader());
if(!aodHeader) AliFatal("Not a standard AOD");
fCentrality = (Int_t)aodHeader->GetCentralityP()->GetCentralityPercentile(fCentralityEstimator.Data());
if( fCentrality < 0 || fCentrality >= 91) return;
fEventCounter->Fill(2);
if(!ProperVertex()) return;
Int_t nTracks = fAOD->GetNumberOfTracks();
TExMap *trackMap = new TExMap();
for(Int_t i = 0; i < nTracks; i++) {
AliAODTrack* aodTrack = dynamic_cast<AliAODTrack *>(fAOD->GetTrack(i));
if(!aodTrack) {
AliError(Form("ERROR: Could not retrieve AODtrack %d",i));
continue;
}
Double_t tpcSignalAll = aodTrack->GetTPCsignal();
fTPCSigA->Fill(aodTrack->GetTPCmomentum(),tpcSignalAll);
Int_t gID = aodTrack->GetID();
if (aodTrack->TestFilterBit(1)) trackMap->Add(gID, i);
}
AliAODTrack* newAodTrack;
for( Int_t j = 0; j < nTracks; j++ ){
AliAODTrack* aodTrack1 = dynamic_cast<AliAODTrack *>(fAOD->GetTrack(j));
if(!aodTrack1) {
AliError(Form("ERROR: Could not retrieve AODtrack %d",j));
continue;
}
if(!aodTrack1->TestFilterBit(fAODtrackCutBit)) continue;
Int_t gID = aodTrack1->GetID();
newAodTrack = gID >= 0 ? aodTrack1 : dynamic_cast<AliAODTrack*>(fAOD->GetTrack(trackMap->GetValue(-1-gID)));
if(!newAodTrack) AliFatal("Not a standard AOD");
Double_t pt = aodTrack1->Pt();
Double_t eta = aodTrack1->Eta();
Double_t nclus = aodTrack1->GetTPCClusterInfo(2,1);
Double_t chi2ndf = aodTrack1->Chi2perNDF();
Float_t dxy = 0., dz = 0.;
if( fAODtrackCutBit == 128 ){
dxy = aodTrack1->DCA();
dz = aodTrack1->ZAtDCA();
if( fabs(dxy) > fDCAxy ) continue;
if( fabs(dz) > fDCAz ) continue;
if( fDCAxy !=0 && fDCAz !=0 ){
if( TMath::Sqrt((dxy*dxy)/(fDCAxy*fDCAxy)+(dz*dz)/(fDCAz*fDCAz)) > 1. ) continue;
}
fHistQA[6]->Fill(dxy);
fHistQA[7]->Fill(dz);
fHistDCA->Fill(dxy,dz);
}
if( pt < fPtLowerLimit || pt > fPtHigherLimit ) continue;
if( eta < fEtaLowerLimit || eta > fEtaHigherLimit ) continue;
if( nclus < fTPCNClus ) continue;
if( chi2ndf > fChi2perNDF ) continue;
fHistQA[8]->Fill(pt);
fHistQA[9]->Fill(eta);
fHistQA[10]->Fill(aodTrack1->Phi());
fHistQA[11]->Fill(nclus);
fHistQA[12]->Fill(chi2ndf);
Short_t gCharge = aodTrack1->Charge();
if( gCharge == 0 ) continue;
if(gCharge > 0) nPlusCharge += 1.;
if(gCharge < 0) nMinusCharge += 1.;
if( fUsePid ) {
gPid = (Int_t)fParticleSpecies;;
Double_t rap = newAodTrack->Y(AliPID::ParticleMass(fParticleSpecies));
Double_t tpcSignal = newAodTrack->GetTPCsignal();
if( fabs(rap ) > fRapidityCut ) continue;
fTPCSig->Fill(newAodTrack->GetTPCmomentum(),tpcSignal);
Float_t nsigmaTPCPID = -999.;
nsigmaTPCPID = TMath::Abs(fPIDResponse->NumberOfSigmasTPC(newAodTrack,fParticleSpecies));
if ( nsigmaTPCPID < fNSigmaCut ){
if (gCharge > 0) nPartile +=1.;
if( gCharge < 0 ) nAntiParticle +=1.;
}
}
}
Double_t fContainerCh[3] = { static_cast<Double_t>(fCentrality), nPlusCharge, nMinusCharge};
Double_t fContainerPid[3] = { static_cast<Double_t>(fCentrality), nPartile, nAntiParticle};
fTHnCentNplusNminusCh->Fill(fContainerCh);
if( fUsePid ){
gPid = (Int_t)fParticleSpecies;
fTHnCentNplusNminusPid[gPid]->Fill(fContainerPid);
}
Int_t nMCTrack = fArrayMC->GetEntriesFast();
for (Int_t iMC = 0; iMC < nMCTrack; iMC++) {
AliAODMCParticle *partMC = (AliAODMCParticle*) fArrayMC->At(iMC);
if(!partMC){
AliError(Form("ERROR: Could not retrieve AODMCtrack %d",iMC));
continue;
}
if(partMC->Charge() == 0) continue;
if(!partMC->IsPhysicalPrimary()) continue;
if (partMC->Eta() < fEtaLowerLimit || partMC->Eta() > fEtaHigherLimit) continue;
if (partMC->Pt() < fPtLowerLimit || partMC->Pt() > fPtHigherLimit) continue;
Short_t gCharge = partMC->Charge();
if(gCharge > 0) nPlusChargeTruth += 1.;
if(gCharge < 0) nMinusChargeTruth += 1.;
if(fUsePid){
Double_t rap = partMC->Y();
if( fabs(rap ) > fRapidityCut ) continue;
if(TMath::Abs(partMC->GetPdgCode()) != gPdgCode) continue;
if(gCharge > 0) nPartileTruth += 1.;
if(gCharge < 0) nAntiParticleTruth += 1.;
}
}
Double_t fContainerChTruth[3] = { static_cast<Double_t>(fCentrality), nPlusChargeTruth, nMinusChargeTruth };
Double_t fContainerPidTruth[3] = { static_cast<Double_t>(fCentrality), nPartileTruth, nAntiParticleTruth };
fTHnCentNplusNminusChTruth->Fill(fContainerChTruth);
if( fUsePid ){
gPid = (Int_t)fParticleSpecies;
fTHnCentNplusNminusPidTruth[gPid]->Fill(fContainerPidTruth);
}
fEventCounter->Fill(7);
return;
}
Bool_t AliHigherMomentsToyModel::ProperVertex(){
Bool_t IsVtx = kFALSE;
const AliAODVertex *vertex = fAOD->GetPrimaryVertex();
if(vertex) {
Double32_t fCov[6];
vertex->GetCovarianceMatrix(fCov);
if(vertex->GetNContributors() > 0) {
if(fCov[5] != 0) {
Double_t lvx = vertex->GetX();
Double_t lvy = vertex->GetY();
Double_t lvz = vertex->GetZ();
fEventCounter->Fill(5);
fHistQA[3]->Fill(lvx);fHistQA[4]->Fill(lvy);fHistQA[5]->Fill(lvz);
if(TMath::Abs(lvx) < fVxMax) {
if(TMath::Abs(lvy) < fVyMax) {
if(TMath::Abs(lvz) < fVzMax) {
fEventCounter->Fill(6);
fHistQA[0]->Fill(lvx);fHistQA[1]->Fill(lvy);fHistQA[2]->Fill(lvz);
IsVtx = kTRUE;
}
}
}
}
}
}
return IsVtx;
}
void AliHigherMomentsToyModel::Terminate( Option_t * ){
Info("AliHigherMomentsToyModel"," Task Successfully finished");
}
AliHigherMomentsToyModel.cxx:1 AliHigherMomentsToyModel.cxx:2 AliHigherMomentsToyModel.cxx:3 AliHigherMomentsToyModel.cxx:4 AliHigherMomentsToyModel.cxx:5 AliHigherMomentsToyModel.cxx:6 AliHigherMomentsToyModel.cxx:7 AliHigherMomentsToyModel.cxx:8 AliHigherMomentsToyModel.cxx:9 AliHigherMomentsToyModel.cxx:10 AliHigherMomentsToyModel.cxx:11 AliHigherMomentsToyModel.cxx:12 AliHigherMomentsToyModel.cxx:13 AliHigherMomentsToyModel.cxx:14 AliHigherMomentsToyModel.cxx:15 AliHigherMomentsToyModel.cxx:16 AliHigherMomentsToyModel.cxx:17 AliHigherMomentsToyModel.cxx:18 AliHigherMomentsToyModel.cxx:19 AliHigherMomentsToyModel.cxx:20 AliHigherMomentsToyModel.cxx:21 AliHigherMomentsToyModel.cxx:22 AliHigherMomentsToyModel.cxx:23 AliHigherMomentsToyModel.cxx:24 AliHigherMomentsToyModel.cxx:25 AliHigherMomentsToyModel.cxx:26 AliHigherMomentsToyModel.cxx:27 AliHigherMomentsToyModel.cxx:28 AliHigherMomentsToyModel.cxx:29 AliHigherMomentsToyModel.cxx:30 AliHigherMomentsToyModel.cxx:31 AliHigherMomentsToyModel.cxx:32 AliHigherMomentsToyModel.cxx:33 AliHigherMomentsToyModel.cxx:34 AliHigherMomentsToyModel.cxx:35 AliHigherMomentsToyModel.cxx:36 AliHigherMomentsToyModel.cxx:37 AliHigherMomentsToyModel.cxx:38 AliHigherMomentsToyModel.cxx:39 AliHigherMomentsToyModel.cxx:40 AliHigherMomentsToyModel.cxx:41 AliHigherMomentsToyModel.cxx:42 AliHigherMomentsToyModel.cxx:43 AliHigherMomentsToyModel.cxx:44 AliHigherMomentsToyModel.cxx:45 AliHigherMomentsToyModel.cxx:46 AliHigherMomentsToyModel.cxx:47 AliHigherMomentsToyModel.cxx:48 AliHigherMomentsToyModel.cxx:49 AliHigherMomentsToyModel.cxx:50 AliHigherMomentsToyModel.cxx:51 AliHigherMomentsToyModel.cxx:52 AliHigherMomentsToyModel.cxx:53 AliHigherMomentsToyModel.cxx:54 AliHigherMomentsToyModel.cxx:55 AliHigherMomentsToyModel.cxx:56 AliHigherMomentsToyModel.cxx:57 AliHigherMomentsToyModel.cxx:58 AliHigherMomentsToyModel.cxx:59 AliHigherMomentsToyModel.cxx:60 AliHigherMomentsToyModel.cxx:61 AliHigherMomentsToyModel.cxx:62 AliHigherMomentsToyModel.cxx:63 AliHigherMomentsToyModel.cxx:64 AliHigherMomentsToyModel.cxx:65 AliHigherMomentsToyModel.cxx:66 AliHigherMomentsToyModel.cxx:67 AliHigherMomentsToyModel.cxx:68 AliHigherMomentsToyModel.cxx:69 AliHigherMomentsToyModel.cxx:70 AliHigherMomentsToyModel.cxx:71 AliHigherMomentsToyModel.cxx:72 AliHigherMomentsToyModel.cxx:73 AliHigherMomentsToyModel.cxx:74 AliHigherMomentsToyModel.cxx:75 AliHigherMomentsToyModel.cxx:76 AliHigherMomentsToyModel.cxx:77 AliHigherMomentsToyModel.cxx:78 AliHigherMomentsToyModel.cxx:79 AliHigherMomentsToyModel.cxx:80 AliHigherMomentsToyModel.cxx:81 AliHigherMomentsToyModel.cxx:82 AliHigherMomentsToyModel.cxx:83 AliHigherMomentsToyModel.cxx:84 AliHigherMomentsToyModel.cxx:85 AliHigherMomentsToyModel.cxx:86 AliHigherMomentsToyModel.cxx:87 AliHigherMomentsToyModel.cxx:88 AliHigherMomentsToyModel.cxx:89 AliHigherMomentsToyModel.cxx:90 AliHigherMomentsToyModel.cxx:91 AliHigherMomentsToyModel.cxx:92 AliHigherMomentsToyModel.cxx:93 AliHigherMomentsToyModel.cxx:94 AliHigherMomentsToyModel.cxx:95 AliHigherMomentsToyModel.cxx:96 AliHigherMomentsToyModel.cxx:97 AliHigherMomentsToyModel.cxx:98 AliHigherMomentsToyModel.cxx:99 AliHigherMomentsToyModel.cxx:100 AliHigherMomentsToyModel.cxx:101 AliHigherMomentsToyModel.cxx:102 AliHigherMomentsToyModel.cxx:103 AliHigherMomentsToyModel.cxx:104 AliHigherMomentsToyModel.cxx:105 AliHigherMomentsToyModel.cxx:106 AliHigherMomentsToyModel.cxx:107 AliHigherMomentsToyModel.cxx:108 AliHigherMomentsToyModel.cxx:109 AliHigherMomentsToyModel.cxx:110 AliHigherMomentsToyModel.cxx:111 AliHigherMomentsToyModel.cxx:112 AliHigherMomentsToyModel.cxx:113 AliHigherMomentsToyModel.cxx:114 AliHigherMomentsToyModel.cxx:115 AliHigherMomentsToyModel.cxx:116 AliHigherMomentsToyModel.cxx:117 AliHigherMomentsToyModel.cxx:118 AliHigherMomentsToyModel.cxx:119 AliHigherMomentsToyModel.cxx:120 AliHigherMomentsToyModel.cxx:121 AliHigherMomentsToyModel.cxx:122 AliHigherMomentsToyModel.cxx:123 AliHigherMomentsToyModel.cxx:124 AliHigherMomentsToyModel.cxx:125 AliHigherMomentsToyModel.cxx:126 AliHigherMomentsToyModel.cxx:127 AliHigherMomentsToyModel.cxx:128 AliHigherMomentsToyModel.cxx:129 AliHigherMomentsToyModel.cxx:130 AliHigherMomentsToyModel.cxx:131 AliHigherMomentsToyModel.cxx:132 AliHigherMomentsToyModel.cxx:133 AliHigherMomentsToyModel.cxx:134 AliHigherMomentsToyModel.cxx:135 AliHigherMomentsToyModel.cxx:136 AliHigherMomentsToyModel.cxx:137 AliHigherMomentsToyModel.cxx:138 AliHigherMomentsToyModel.cxx:139 AliHigherMomentsToyModel.cxx:140 AliHigherMomentsToyModel.cxx:141 AliHigherMomentsToyModel.cxx:142 AliHigherMomentsToyModel.cxx:143 AliHigherMomentsToyModel.cxx:144 AliHigherMomentsToyModel.cxx:145 AliHigherMomentsToyModel.cxx:146 AliHigherMomentsToyModel.cxx:147 AliHigherMomentsToyModel.cxx:148 AliHigherMomentsToyModel.cxx:149 AliHigherMomentsToyModel.cxx:150 AliHigherMomentsToyModel.cxx:151 AliHigherMomentsToyModel.cxx:152 AliHigherMomentsToyModel.cxx:153 AliHigherMomentsToyModel.cxx:154 AliHigherMomentsToyModel.cxx:155 AliHigherMomentsToyModel.cxx:156 AliHigherMomentsToyModel.cxx:157 AliHigherMomentsToyModel.cxx:158 AliHigherMomentsToyModel.cxx:159 AliHigherMomentsToyModel.cxx:160 AliHigherMomentsToyModel.cxx:161 AliHigherMomentsToyModel.cxx:162 AliHigherMomentsToyModel.cxx:163 AliHigherMomentsToyModel.cxx:164 AliHigherMomentsToyModel.cxx:165 AliHigherMomentsToyModel.cxx:166 AliHigherMomentsToyModel.cxx:167 AliHigherMomentsToyModel.cxx:168 AliHigherMomentsToyModel.cxx:169 AliHigherMomentsToyModel.cxx:170 AliHigherMomentsToyModel.cxx:171 AliHigherMomentsToyModel.cxx:172 AliHigherMomentsToyModel.cxx:173 AliHigherMomentsToyModel.cxx:174 AliHigherMomentsToyModel.cxx:175 AliHigherMomentsToyModel.cxx:176 AliHigherMomentsToyModel.cxx:177 AliHigherMomentsToyModel.cxx:178 AliHigherMomentsToyModel.cxx:179 AliHigherMomentsToyModel.cxx:180 AliHigherMomentsToyModel.cxx:181 AliHigherMomentsToyModel.cxx:182 AliHigherMomentsToyModel.cxx:183 AliHigherMomentsToyModel.cxx:184 AliHigherMomentsToyModel.cxx:185 AliHigherMomentsToyModel.cxx:186 AliHigherMomentsToyModel.cxx:187 AliHigherMomentsToyModel.cxx:188 AliHigherMomentsToyModel.cxx:189 AliHigherMomentsToyModel.cxx:190 AliHigherMomentsToyModel.cxx:191 AliHigherMomentsToyModel.cxx:192 AliHigherMomentsToyModel.cxx:193 AliHigherMomentsToyModel.cxx:194 AliHigherMomentsToyModel.cxx:195 AliHigherMomentsToyModel.cxx:196 AliHigherMomentsToyModel.cxx:197 AliHigherMomentsToyModel.cxx:198 AliHigherMomentsToyModel.cxx:199 AliHigherMomentsToyModel.cxx:200 AliHigherMomentsToyModel.cxx:201 AliHigherMomentsToyModel.cxx:202 AliHigherMomentsToyModel.cxx:203 AliHigherMomentsToyModel.cxx:204 AliHigherMomentsToyModel.cxx:205 AliHigherMomentsToyModel.cxx:206 AliHigherMomentsToyModel.cxx:207 AliHigherMomentsToyModel.cxx:208 AliHigherMomentsToyModel.cxx:209 AliHigherMomentsToyModel.cxx:210 AliHigherMomentsToyModel.cxx:211 AliHigherMomentsToyModel.cxx:212 AliHigherMomentsToyModel.cxx:213 AliHigherMomentsToyModel.cxx:214 AliHigherMomentsToyModel.cxx:215 AliHigherMomentsToyModel.cxx:216 AliHigherMomentsToyModel.cxx:217 AliHigherMomentsToyModel.cxx:218 AliHigherMomentsToyModel.cxx:219 AliHigherMomentsToyModel.cxx:220 AliHigherMomentsToyModel.cxx:221 AliHigherMomentsToyModel.cxx:222 AliHigherMomentsToyModel.cxx:223 AliHigherMomentsToyModel.cxx:224 AliHigherMomentsToyModel.cxx:225 AliHigherMomentsToyModel.cxx:226 AliHigherMomentsToyModel.cxx:227 AliHigherMomentsToyModel.cxx:228 AliHigherMomentsToyModel.cxx:229 AliHigherMomentsToyModel.cxx:230 AliHigherMomentsToyModel.cxx:231 AliHigherMomentsToyModel.cxx:232 AliHigherMomentsToyModel.cxx:233 AliHigherMomentsToyModel.cxx:234 AliHigherMomentsToyModel.cxx:235 AliHigherMomentsToyModel.cxx:236 AliHigherMomentsToyModel.cxx:237 AliHigherMomentsToyModel.cxx:238 AliHigherMomentsToyModel.cxx:239 AliHigherMomentsToyModel.cxx:240 AliHigherMomentsToyModel.cxx:241 AliHigherMomentsToyModel.cxx:242 AliHigherMomentsToyModel.cxx:243 AliHigherMomentsToyModel.cxx:244 AliHigherMomentsToyModel.cxx:245 AliHigherMomentsToyModel.cxx:246 AliHigherMomentsToyModel.cxx:247 AliHigherMomentsToyModel.cxx:248 AliHigherMomentsToyModel.cxx:249 AliHigherMomentsToyModel.cxx:250 AliHigherMomentsToyModel.cxx:251 AliHigherMomentsToyModel.cxx:252 AliHigherMomentsToyModel.cxx:253 AliHigherMomentsToyModel.cxx:254 AliHigherMomentsToyModel.cxx:255 AliHigherMomentsToyModel.cxx:256 AliHigherMomentsToyModel.cxx:257 AliHigherMomentsToyModel.cxx:258 AliHigherMomentsToyModel.cxx:259 AliHigherMomentsToyModel.cxx:260 AliHigherMomentsToyModel.cxx:261 AliHigherMomentsToyModel.cxx:262 AliHigherMomentsToyModel.cxx:263 AliHigherMomentsToyModel.cxx:264 AliHigherMomentsToyModel.cxx:265 AliHigherMomentsToyModel.cxx:266 AliHigherMomentsToyModel.cxx:267 AliHigherMomentsToyModel.cxx:268 AliHigherMomentsToyModel.cxx:269 AliHigherMomentsToyModel.cxx:270 AliHigherMomentsToyModel.cxx:271 AliHigherMomentsToyModel.cxx:272 AliHigherMomentsToyModel.cxx:273 AliHigherMomentsToyModel.cxx:274 AliHigherMomentsToyModel.cxx:275 AliHigherMomentsToyModel.cxx:276 AliHigherMomentsToyModel.cxx:277 AliHigherMomentsToyModel.cxx:278 AliHigherMomentsToyModel.cxx:279 AliHigherMomentsToyModel.cxx:280 AliHigherMomentsToyModel.cxx:281 AliHigherMomentsToyModel.cxx:282 AliHigherMomentsToyModel.cxx:283 AliHigherMomentsToyModel.cxx:284 AliHigherMomentsToyModel.cxx:285 AliHigherMomentsToyModel.cxx:286 AliHigherMomentsToyModel.cxx:287 AliHigherMomentsToyModel.cxx:288 AliHigherMomentsToyModel.cxx:289 AliHigherMomentsToyModel.cxx:290 AliHigherMomentsToyModel.cxx:291 AliHigherMomentsToyModel.cxx:292 AliHigherMomentsToyModel.cxx:293 AliHigherMomentsToyModel.cxx:294 AliHigherMomentsToyModel.cxx:295 AliHigherMomentsToyModel.cxx:296 AliHigherMomentsToyModel.cxx:297 AliHigherMomentsToyModel.cxx:298 AliHigherMomentsToyModel.cxx:299 AliHigherMomentsToyModel.cxx:300 AliHigherMomentsToyModel.cxx:301 AliHigherMomentsToyModel.cxx:302 AliHigherMomentsToyModel.cxx:303 AliHigherMomentsToyModel.cxx:304 AliHigherMomentsToyModel.cxx:305 AliHigherMomentsToyModel.cxx:306 AliHigherMomentsToyModel.cxx:307 AliHigherMomentsToyModel.cxx:308 AliHigherMomentsToyModel.cxx:309 AliHigherMomentsToyModel.cxx:310 AliHigherMomentsToyModel.cxx:311 AliHigherMomentsToyModel.cxx:312 AliHigherMomentsToyModel.cxx:313 AliHigherMomentsToyModel.cxx:314 AliHigherMomentsToyModel.cxx:315 AliHigherMomentsToyModel.cxx:316 AliHigherMomentsToyModel.cxx:317 AliHigherMomentsToyModel.cxx:318 AliHigherMomentsToyModel.cxx:319 AliHigherMomentsToyModel.cxx:320 AliHigherMomentsToyModel.cxx:321 AliHigherMomentsToyModel.cxx:322 AliHigherMomentsToyModel.cxx:323 AliHigherMomentsToyModel.cxx:324 AliHigherMomentsToyModel.cxx:325 AliHigherMomentsToyModel.cxx:326 AliHigherMomentsToyModel.cxx:327 AliHigherMomentsToyModel.cxx:328 AliHigherMomentsToyModel.cxx:329 AliHigherMomentsToyModel.cxx:330 AliHigherMomentsToyModel.cxx:331 AliHigherMomentsToyModel.cxx:332 AliHigherMomentsToyModel.cxx:333 AliHigherMomentsToyModel.cxx:334 AliHigherMomentsToyModel.cxx:335 AliHigherMomentsToyModel.cxx:336 AliHigherMomentsToyModel.cxx:337 AliHigherMomentsToyModel.cxx:338 AliHigherMomentsToyModel.cxx:339 AliHigherMomentsToyModel.cxx:340 AliHigherMomentsToyModel.cxx:341 AliHigherMomentsToyModel.cxx:342 AliHigherMomentsToyModel.cxx:343 AliHigherMomentsToyModel.cxx:344 AliHigherMomentsToyModel.cxx:345 AliHigherMomentsToyModel.cxx:346 AliHigherMomentsToyModel.cxx:347 AliHigherMomentsToyModel.cxx:348 AliHigherMomentsToyModel.cxx:349 AliHigherMomentsToyModel.cxx:350 AliHigherMomentsToyModel.cxx:351 AliHigherMomentsToyModel.cxx:352 AliHigherMomentsToyModel.cxx:353 AliHigherMomentsToyModel.cxx:354 AliHigherMomentsToyModel.cxx:355 AliHigherMomentsToyModel.cxx:356 AliHigherMomentsToyModel.cxx:357 AliHigherMomentsToyModel.cxx:358 AliHigherMomentsToyModel.cxx:359 AliHigherMomentsToyModel.cxx:360 AliHigherMomentsToyModel.cxx:361 AliHigherMomentsToyModel.cxx:362 AliHigherMomentsToyModel.cxx:363 AliHigherMomentsToyModel.cxx:364 AliHigherMomentsToyModel.cxx:365 AliHigherMomentsToyModel.cxx:366 AliHigherMomentsToyModel.cxx:367 AliHigherMomentsToyModel.cxx:368 AliHigherMomentsToyModel.cxx:369 AliHigherMomentsToyModel.cxx:370 AliHigherMomentsToyModel.cxx:371 AliHigherMomentsToyModel.cxx:372 AliHigherMomentsToyModel.cxx:373 AliHigherMomentsToyModel.cxx:374 AliHigherMomentsToyModel.cxx:375 AliHigherMomentsToyModel.cxx:376 AliHigherMomentsToyModel.cxx:377 AliHigherMomentsToyModel.cxx:378 AliHigherMomentsToyModel.cxx:379 AliHigherMomentsToyModel.cxx:380 AliHigherMomentsToyModel.cxx:381 AliHigherMomentsToyModel.cxx:382 AliHigherMomentsToyModel.cxx:383 AliHigherMomentsToyModel.cxx:384 AliHigherMomentsToyModel.cxx:385 AliHigherMomentsToyModel.cxx:386 AliHigherMomentsToyModel.cxx:387 AliHigherMomentsToyModel.cxx:388 AliHigherMomentsToyModel.cxx:389 AliHigherMomentsToyModel.cxx:390 AliHigherMomentsToyModel.cxx:391 AliHigherMomentsToyModel.cxx:392 AliHigherMomentsToyModel.cxx:393 AliHigherMomentsToyModel.cxx:394 AliHigherMomentsToyModel.cxx:395 AliHigherMomentsToyModel.cxx:396 AliHigherMomentsToyModel.cxx:397 AliHigherMomentsToyModel.cxx:398 AliHigherMomentsToyModel.cxx:399 AliHigherMomentsToyModel.cxx:400 AliHigherMomentsToyModel.cxx:401 AliHigherMomentsToyModel.cxx:402 AliHigherMomentsToyModel.cxx:403 AliHigherMomentsToyModel.cxx:404 AliHigherMomentsToyModel.cxx:405 AliHigherMomentsToyModel.cxx:406 AliHigherMomentsToyModel.cxx:407 AliHigherMomentsToyModel.cxx:408 AliHigherMomentsToyModel.cxx:409 AliHigherMomentsToyModel.cxx:410 AliHigherMomentsToyModel.cxx:411 AliHigherMomentsToyModel.cxx:412 AliHigherMomentsToyModel.cxx:413 AliHigherMomentsToyModel.cxx:414 AliHigherMomentsToyModel.cxx:415 AliHigherMomentsToyModel.cxx:416 AliHigherMomentsToyModel.cxx:417 AliHigherMomentsToyModel.cxx:418 AliHigherMomentsToyModel.cxx:419 AliHigherMomentsToyModel.cxx:420 AliHigherMomentsToyModel.cxx:421 AliHigherMomentsToyModel.cxx:422 AliHigherMomentsToyModel.cxx:423 AliHigherMomentsToyModel.cxx:424 AliHigherMomentsToyModel.cxx:425 AliHigherMomentsToyModel.cxx:426 AliHigherMomentsToyModel.cxx:427 AliHigherMomentsToyModel.cxx:428 AliHigherMomentsToyModel.cxx:429 AliHigherMomentsToyModel.cxx:430 AliHigherMomentsToyModel.cxx:431 AliHigherMomentsToyModel.cxx:432 AliHigherMomentsToyModel.cxx:433 AliHigherMomentsToyModel.cxx:434 AliHigherMomentsToyModel.cxx:435 AliHigherMomentsToyModel.cxx:436 AliHigherMomentsToyModel.cxx:437 AliHigherMomentsToyModel.cxx:438 AliHigherMomentsToyModel.cxx:439 AliHigherMomentsToyModel.cxx:440 AliHigherMomentsToyModel.cxx:441 AliHigherMomentsToyModel.cxx:442 AliHigherMomentsToyModel.cxx:443 AliHigherMomentsToyModel.cxx:444 AliHigherMomentsToyModel.cxx:445 AliHigherMomentsToyModel.cxx:446 AliHigherMomentsToyModel.cxx:447 AliHigherMomentsToyModel.cxx:448 AliHigherMomentsToyModel.cxx:449 AliHigherMomentsToyModel.cxx:450 AliHigherMomentsToyModel.cxx:451 AliHigherMomentsToyModel.cxx:452 AliHigherMomentsToyModel.cxx:453 AliHigherMomentsToyModel.cxx:454 AliHigherMomentsToyModel.cxx:455 AliHigherMomentsToyModel.cxx:456 AliHigherMomentsToyModel.cxx:457 AliHigherMomentsToyModel.cxx:458 AliHigherMomentsToyModel.cxx:459 AliHigherMomentsToyModel.cxx:460 AliHigherMomentsToyModel.cxx:461 AliHigherMomentsToyModel.cxx:462 AliHigherMomentsToyModel.cxx:463 AliHigherMomentsToyModel.cxx:464 AliHigherMomentsToyModel.cxx:465 AliHigherMomentsToyModel.cxx:466 AliHigherMomentsToyModel.cxx:467 AliHigherMomentsToyModel.cxx:468 AliHigherMomentsToyModel.cxx:469 AliHigherMomentsToyModel.cxx:470 AliHigherMomentsToyModel.cxx:471 AliHigherMomentsToyModel.cxx:472 AliHigherMomentsToyModel.cxx:473 AliHigherMomentsToyModel.cxx:474 AliHigherMomentsToyModel.cxx:475 AliHigherMomentsToyModel.cxx:476 AliHigherMomentsToyModel.cxx:477 AliHigherMomentsToyModel.cxx:478 AliHigherMomentsToyModel.cxx:479 AliHigherMomentsToyModel.cxx:480 AliHigherMomentsToyModel.cxx:481 AliHigherMomentsToyModel.cxx:482 AliHigherMomentsToyModel.cxx:483 AliHigherMomentsToyModel.cxx:484 AliHigherMomentsToyModel.cxx:485 AliHigherMomentsToyModel.cxx:486 AliHigherMomentsToyModel.cxx:487 AliHigherMomentsToyModel.cxx:488 AliHigherMomentsToyModel.cxx:489 AliHigherMomentsToyModel.cxx:490 AliHigherMomentsToyModel.cxx:491 AliHigherMomentsToyModel.cxx:492 AliHigherMomentsToyModel.cxx:493 AliHigherMomentsToyModel.cxx:494 AliHigherMomentsToyModel.cxx:495 AliHigherMomentsToyModel.cxx:496 AliHigherMomentsToyModel.cxx:497 AliHigherMomentsToyModel.cxx:498 AliHigherMomentsToyModel.cxx:499 AliHigherMomentsToyModel.cxx:500 AliHigherMomentsToyModel.cxx:501 AliHigherMomentsToyModel.cxx:502 AliHigherMomentsToyModel.cxx:503 AliHigherMomentsToyModel.cxx:504 AliHigherMomentsToyModel.cxx:505 AliHigherMomentsToyModel.cxx:506 AliHigherMomentsToyModel.cxx:507 AliHigherMomentsToyModel.cxx:508 AliHigherMomentsToyModel.cxx:509 AliHigherMomentsToyModel.cxx:510 AliHigherMomentsToyModel.cxx:511 AliHigherMomentsToyModel.cxx:512 AliHigherMomentsToyModel.cxx:513 AliHigherMomentsToyModel.cxx:514 AliHigherMomentsToyModel.cxx:515 AliHigherMomentsToyModel.cxx:516 AliHigherMomentsToyModel.cxx:517 AliHigherMomentsToyModel.cxx:518 AliHigherMomentsToyModel.cxx:519 AliHigherMomentsToyModel.cxx:520 AliHigherMomentsToyModel.cxx:521 AliHigherMomentsToyModel.cxx:522 AliHigherMomentsToyModel.cxx:523 AliHigherMomentsToyModel.cxx:524 AliHigherMomentsToyModel.cxx:525 AliHigherMomentsToyModel.cxx:526 AliHigherMomentsToyModel.cxx:527 AliHigherMomentsToyModel.cxx:528 AliHigherMomentsToyModel.cxx:529 AliHigherMomentsToyModel.cxx:530 AliHigherMomentsToyModel.cxx:531 AliHigherMomentsToyModel.cxx:532 AliHigherMomentsToyModel.cxx:533 AliHigherMomentsToyModel.cxx:534 AliHigherMomentsToyModel.cxx:535 AliHigherMomentsToyModel.cxx:536 AliHigherMomentsToyModel.cxx:537 AliHigherMomentsToyModel.cxx:538 AliHigherMomentsToyModel.cxx:539 AliHigherMomentsToyModel.cxx:540 AliHigherMomentsToyModel.cxx:541 AliHigherMomentsToyModel.cxx:542 AliHigherMomentsToyModel.cxx:543 AliHigherMomentsToyModel.cxx:544 AliHigherMomentsToyModel.cxx:545 AliHigherMomentsToyModel.cxx:546 AliHigherMomentsToyModel.cxx:547 AliHigherMomentsToyModel.cxx:548 AliHigherMomentsToyModel.cxx:549 AliHigherMomentsToyModel.cxx:550 AliHigherMomentsToyModel.cxx:551 AliHigherMomentsToyModel.cxx:552 AliHigherMomentsToyModel.cxx:553 AliHigherMomentsToyModel.cxx:554 AliHigherMomentsToyModel.cxx:555 AliHigherMomentsToyModel.cxx:556 AliHigherMomentsToyModel.cxx:557 AliHigherMomentsToyModel.cxx:558 AliHigherMomentsToyModel.cxx:559 AliHigherMomentsToyModel.cxx:560 AliHigherMomentsToyModel.cxx:561 AliHigherMomentsToyModel.cxx:562 AliHigherMomentsToyModel.cxx:563 AliHigherMomentsToyModel.cxx:564 AliHigherMomentsToyModel.cxx:565 AliHigherMomentsToyModel.cxx:566 AliHigherMomentsToyModel.cxx:567 AliHigherMomentsToyModel.cxx:568 AliHigherMomentsToyModel.cxx:569 AliHigherMomentsToyModel.cxx:570 AliHigherMomentsToyModel.cxx:571 AliHigherMomentsToyModel.cxx:572 AliHigherMomentsToyModel.cxx:573 AliHigherMomentsToyModel.cxx:574 AliHigherMomentsToyModel.cxx:575 AliHigherMomentsToyModel.cxx:576 AliHigherMomentsToyModel.cxx:577 AliHigherMomentsToyModel.cxx:578 AliHigherMomentsToyModel.cxx:579 AliHigherMomentsToyModel.cxx:580 AliHigherMomentsToyModel.cxx:581 AliHigherMomentsToyModel.cxx:582 AliHigherMomentsToyModel.cxx:583 AliHigherMomentsToyModel.cxx:584 AliHigherMomentsToyModel.cxx:585 AliHigherMomentsToyModel.cxx:586 AliHigherMomentsToyModel.cxx:587 AliHigherMomentsToyModel.cxx:588 AliHigherMomentsToyModel.cxx:589 AliHigherMomentsToyModel.cxx:590 AliHigherMomentsToyModel.cxx:591 AliHigherMomentsToyModel.cxx:592 AliHigherMomentsToyModel.cxx:593 AliHigherMomentsToyModel.cxx:594 AliHigherMomentsToyModel.cxx:595 AliHigherMomentsToyModel.cxx:596 AliHigherMomentsToyModel.cxx:597 AliHigherMomentsToyModel.cxx:598 AliHigherMomentsToyModel.cxx:599 AliHigherMomentsToyModel.cxx:600 AliHigherMomentsToyModel.cxx:601 AliHigherMomentsToyModel.cxx:602 AliHigherMomentsToyModel.cxx:603 AliHigherMomentsToyModel.cxx:604 AliHigherMomentsToyModel.cxx:605 AliHigherMomentsToyModel.cxx:606 AliHigherMomentsToyModel.cxx:607 AliHigherMomentsToyModel.cxx:608 AliHigherMomentsToyModel.cxx:609 AliHigherMomentsToyModel.cxx:610 AliHigherMomentsToyModel.cxx:611 AliHigherMomentsToyModel.cxx:612 AliHigherMomentsToyModel.cxx:613 AliHigherMomentsToyModel.cxx:614 AliHigherMomentsToyModel.cxx:615 AliHigherMomentsToyModel.cxx:616 AliHigherMomentsToyModel.cxx:617 AliHigherMomentsToyModel.cxx:618 AliHigherMomentsToyModel.cxx:619 AliHigherMomentsToyModel.cxx:620 AliHigherMomentsToyModel.cxx:621 AliHigherMomentsToyModel.cxx:622 AliHigherMomentsToyModel.cxx:623 AliHigherMomentsToyModel.cxx:624 AliHigherMomentsToyModel.cxx:625 AliHigherMomentsToyModel.cxx:626 AliHigherMomentsToyModel.cxx:627 AliHigherMomentsToyModel.cxx:628 AliHigherMomentsToyModel.cxx:629 AliHigherMomentsToyModel.cxx:630 AliHigherMomentsToyModel.cxx:631 AliHigherMomentsToyModel.cxx:632 AliHigherMomentsToyModel.cxx:633 AliHigherMomentsToyModel.cxx:634 AliHigherMomentsToyModel.cxx:635 AliHigherMomentsToyModel.cxx:636 AliHigherMomentsToyModel.cxx:637 AliHigherMomentsToyModel.cxx:638 AliHigherMomentsToyModel.cxx:639 AliHigherMomentsToyModel.cxx:640 AliHigherMomentsToyModel.cxx:641 AliHigherMomentsToyModel.cxx:642 AliHigherMomentsToyModel.cxx:643 AliHigherMomentsToyModel.cxx:644 AliHigherMomentsToyModel.cxx:645 AliHigherMomentsToyModel.cxx:646 AliHigherMomentsToyModel.cxx:647 AliHigherMomentsToyModel.cxx:648 AliHigherMomentsToyModel.cxx:649 AliHigherMomentsToyModel.cxx:650 AliHigherMomentsToyModel.cxx:651 AliHigherMomentsToyModel.cxx:652 AliHigherMomentsToyModel.cxx:653 AliHigherMomentsToyModel.cxx:654 AliHigherMomentsToyModel.cxx:655 AliHigherMomentsToyModel.cxx:656 AliHigherMomentsToyModel.cxx:657 AliHigherMomentsToyModel.cxx:658 AliHigherMomentsToyModel.cxx:659 AliHigherMomentsToyModel.cxx:660 AliHigherMomentsToyModel.cxx:661 AliHigherMomentsToyModel.cxx:662 AliHigherMomentsToyModel.cxx:663 AliHigherMomentsToyModel.cxx:664 AliHigherMomentsToyModel.cxx:665 AliHigherMomentsToyModel.cxx:666 AliHigherMomentsToyModel.cxx:667 AliHigherMomentsToyModel.cxx:668 AliHigherMomentsToyModel.cxx:669 AliHigherMomentsToyModel.cxx:670 AliHigherMomentsToyModel.cxx:671 AliHigherMomentsToyModel.cxx:672 AliHigherMomentsToyModel.cxx:673 AliHigherMomentsToyModel.cxx:674 AliHigherMomentsToyModel.cxx:675 AliHigherMomentsToyModel.cxx:676 AliHigherMomentsToyModel.cxx:677 AliHigherMomentsToyModel.cxx:678 AliHigherMomentsToyModel.cxx:679 AliHigherMomentsToyModel.cxx:680 AliHigherMomentsToyModel.cxx:681 AliHigherMomentsToyModel.cxx:682 AliHigherMomentsToyModel.cxx:683 AliHigherMomentsToyModel.cxx:684 AliHigherMomentsToyModel.cxx:685 AliHigherMomentsToyModel.cxx:686 AliHigherMomentsToyModel.cxx:687 AliHigherMomentsToyModel.cxx:688 AliHigherMomentsToyModel.cxx:689 AliHigherMomentsToyModel.cxx:690 AliHigherMomentsToyModel.cxx:691 AliHigherMomentsToyModel.cxx:692 AliHigherMomentsToyModel.cxx:693 AliHigherMomentsToyModel.cxx:694 AliHigherMomentsToyModel.cxx:695 AliHigherMomentsToyModel.cxx:696 AliHigherMomentsToyModel.cxx:697 AliHigherMomentsToyModel.cxx:698 AliHigherMomentsToyModel.cxx:699 AliHigherMomentsToyModel.cxx:700 AliHigherMomentsToyModel.cxx:701 AliHigherMomentsToyModel.cxx:702 AliHigherMomentsToyModel.cxx:703 AliHigherMomentsToyModel.cxx:704 AliHigherMomentsToyModel.cxx:705 AliHigherMomentsToyModel.cxx:706 AliHigherMomentsToyModel.cxx:707 AliHigherMomentsToyModel.cxx:708 AliHigherMomentsToyModel.cxx:709 AliHigherMomentsToyModel.cxx:710 AliHigherMomentsToyModel.cxx:711 AliHigherMomentsToyModel.cxx:712 AliHigherMomentsToyModel.cxx:713 AliHigherMomentsToyModel.cxx:714 AliHigherMomentsToyModel.cxx:715 AliHigherMomentsToyModel.cxx:716 AliHigherMomentsToyModel.cxx:717 AliHigherMomentsToyModel.cxx:718 AliHigherMomentsToyModel.cxx:719 AliHigherMomentsToyModel.cxx:720 AliHigherMomentsToyModel.cxx:721 AliHigherMomentsToyModel.cxx:722 AliHigherMomentsToyModel.cxx:723 AliHigherMomentsToyModel.cxx:724 AliHigherMomentsToyModel.cxx:725 AliHigherMomentsToyModel.cxx:726 AliHigherMomentsToyModel.cxx:727 AliHigherMomentsToyModel.cxx:728 AliHigherMomentsToyModel.cxx:729 AliHigherMomentsToyModel.cxx:730 AliHigherMomentsToyModel.cxx:731 AliHigherMomentsToyModel.cxx:732 AliHigherMomentsToyModel.cxx:733 AliHigherMomentsToyModel.cxx:734 AliHigherMomentsToyModel.cxx:735 AliHigherMomentsToyModel.cxx:736 AliHigherMomentsToyModel.cxx:737 AliHigherMomentsToyModel.cxx:738 AliHigherMomentsToyModel.cxx:739 AliHigherMomentsToyModel.cxx:740 AliHigherMomentsToyModel.cxx:741 AliHigherMomentsToyModel.cxx:742 AliHigherMomentsToyModel.cxx:743 AliHigherMomentsToyModel.cxx:744 AliHigherMomentsToyModel.cxx:745 AliHigherMomentsToyModel.cxx:746 AliHigherMomentsToyModel.cxx:747 AliHigherMomentsToyModel.cxx:748 AliHigherMomentsToyModel.cxx:749 AliHigherMomentsToyModel.cxx:750 AliHigherMomentsToyModel.cxx:751 AliHigherMomentsToyModel.cxx:752 AliHigherMomentsToyModel.cxx:753 AliHigherMomentsToyModel.cxx:754 AliHigherMomentsToyModel.cxx:755 AliHigherMomentsToyModel.cxx:756 AliHigherMomentsToyModel.cxx:757 AliHigherMomentsToyModel.cxx:758 AliHigherMomentsToyModel.cxx:759 AliHigherMomentsToyModel.cxx:760 AliHigherMomentsToyModel.cxx:761 AliHigherMomentsToyModel.cxx:762 AliHigherMomentsToyModel.cxx:763 AliHigherMomentsToyModel.cxx:764 AliHigherMomentsToyModel.cxx:765 AliHigherMomentsToyModel.cxx:766 AliHigherMomentsToyModel.cxx:767 AliHigherMomentsToyModel.cxx:768 AliHigherMomentsToyModel.cxx:769 AliHigherMomentsToyModel.cxx:770 AliHigherMomentsToyModel.cxx:771 AliHigherMomentsToyModel.cxx:772 AliHigherMomentsToyModel.cxx:773 AliHigherMomentsToyModel.cxx:774 AliHigherMomentsToyModel.cxx:775 AliHigherMomentsToyModel.cxx:776 AliHigherMomentsToyModel.cxx:777 AliHigherMomentsToyModel.cxx:778 AliHigherMomentsToyModel.cxx:779 AliHigherMomentsToyModel.cxx:780 AliHigherMomentsToyModel.cxx:781 AliHigherMomentsToyModel.cxx:782 AliHigherMomentsToyModel.cxx:783 AliHigherMomentsToyModel.cxx:784 AliHigherMomentsToyModel.cxx:785 AliHigherMomentsToyModel.cxx:786 AliHigherMomentsToyModel.cxx:787 AliHigherMomentsToyModel.cxx:788 AliHigherMomentsToyModel.cxx:789 AliHigherMomentsToyModel.cxx:790 AliHigherMomentsToyModel.cxx:791 AliHigherMomentsToyModel.cxx:792 AliHigherMomentsToyModel.cxx:793 AliHigherMomentsToyModel.cxx:794 AliHigherMomentsToyModel.cxx:795 AliHigherMomentsToyModel.cxx:796 AliHigherMomentsToyModel.cxx:797 AliHigherMomentsToyModel.cxx:798 AliHigherMomentsToyModel.cxx:799 AliHigherMomentsToyModel.cxx:800 AliHigherMomentsToyModel.cxx:801 AliHigherMomentsToyModel.cxx:802 AliHigherMomentsToyModel.cxx:803 AliHigherMomentsToyModel.cxx:804 AliHigherMomentsToyModel.cxx:805 AliHigherMomentsToyModel.cxx:806 AliHigherMomentsToyModel.cxx:807 AliHigherMomentsToyModel.cxx:808 AliHigherMomentsToyModel.cxx:809 AliHigherMomentsToyModel.cxx:810 AliHigherMomentsToyModel.cxx:811 AliHigherMomentsToyModel.cxx:812 AliHigherMomentsToyModel.cxx:813 AliHigherMomentsToyModel.cxx:814 AliHigherMomentsToyModel.cxx:815 AliHigherMomentsToyModel.cxx:816 AliHigherMomentsToyModel.cxx:817 AliHigherMomentsToyModel.cxx:818 AliHigherMomentsToyModel.cxx:819 AliHigherMomentsToyModel.cxx:820 AliHigherMomentsToyModel.cxx:821 AliHigherMomentsToyModel.cxx:822 AliHigherMomentsToyModel.cxx:823 AliHigherMomentsToyModel.cxx:824 AliHigherMomentsToyModel.cxx:825 AliHigherMomentsToyModel.cxx:826 AliHigherMomentsToyModel.cxx:827 AliHigherMomentsToyModel.cxx:828 AliHigherMomentsToyModel.cxx:829 AliHigherMomentsToyModel.cxx:830 AliHigherMomentsToyModel.cxx:831 AliHigherMomentsToyModel.cxx:832 AliHigherMomentsToyModel.cxx:833 AliHigherMomentsToyModel.cxx:834 AliHigherMomentsToyModel.cxx:835 AliHigherMomentsToyModel.cxx:836 AliHigherMomentsToyModel.cxx:837 AliHigherMomentsToyModel.cxx:838 AliHigherMomentsToyModel.cxx:839 AliHigherMomentsToyModel.cxx:840 AliHigherMomentsToyModel.cxx:841 AliHigherMomentsToyModel.cxx:842 AliHigherMomentsToyModel.cxx:843 AliHigherMomentsToyModel.cxx:844 AliHigherMomentsToyModel.cxx:845 AliHigherMomentsToyModel.cxx:846 AliHigherMomentsToyModel.cxx:847 AliHigherMomentsToyModel.cxx:848 AliHigherMomentsToyModel.cxx:849 AliHigherMomentsToyModel.cxx:850 AliHigherMomentsToyModel.cxx:851 AliHigherMomentsToyModel.cxx:852 AliHigherMomentsToyModel.cxx:853 AliHigherMomentsToyModel.cxx:854 AliHigherMomentsToyModel.cxx:855 AliHigherMomentsToyModel.cxx:856