From: sjena Date: Fri, 31 Oct 2014 14:43:36 +0000 (+0100) Subject: Update in Two Particle Corr: debojit X-Git-Url: http://git.uio.no/git/?p=u%2Fmrichter%2FAliRoot.git;a=commitdiff_plain;h=73792dfe5103b320df0c5e3996df96224f0b3802 Update in Two Particle Corr: debojit --- diff --git a/PWGCF/Correlations/DPhi/TriggerPID/AliTwoParticlePIDCorr.cxx b/PWGCF/Correlations/DPhi/TriggerPID/AliTwoParticlePIDCorr.cxx index 698c84c6b4c..e6d28de348d 100644 --- a/PWGCF/Correlations/DPhi/TriggerPID/AliTwoParticlePIDCorr.cxx +++ b/PWGCF/Correlations/DPhi/TriggerPID/AliTwoParticlePIDCorr.cxx @@ -30,6 +30,9 @@ #include "AliCentrality.h" #include "Riostream.h" +#include "AliAnalysisDataSlot.h" + #include "AliAnalysisDataContainer.h" + #include "AliTHn.h" #include "AliCFContainer.h" #include "THn.h" @@ -96,6 +99,7 @@ AliTwoParticlePIDCorr::AliTwoParticlePIDCorr() // All data members should be ini fPPVsMultUtils(kFALSE), fSampleType("pPb"), fRequestEventPlane(kFALSE), + fRequestEventPlanemixing(kFALSE), fnTracksVertex(1), // QA tracks pointing to principal vertex (= 3 default) trkVtx(0), zvtx(0), @@ -372,6 +376,7 @@ AliTwoParticlePIDCorr::AliTwoParticlePIDCorr(const char *name) // All data membe fPPVsMultUtils(kFALSE), fSampleType("pPb"), fRequestEventPlane(kFALSE), + fRequestEventPlanemixing(kFALSE), fnTracksVertex(1), // QA tracks pointing to principal vertex (= 3 default) trkVtx(0), zvtx(0), @@ -948,13 +953,18 @@ for(Int_t i = 0; i < 16; i++) "delta_phi: -1.570796, -1.483530, -1.396263, -1.308997, -1.221730, -1.134464, -1.047198, -0.959931, -0.872665, -0.785398, -0.698132, -0.610865, -0.523599, -0.436332, -0.349066, -0.261799, -0.174533, -0.087266, 0.0, 0.087266, 0.174533, 0.261799, 0.349066, 0.436332, 0.523599, 0.610865, 0.698132, 0.785398, 0.872665, 0.959931, 1.047198, 1.134464, 1.221730, 1.308997, 1.396263, 1.483530, 1.570796, 1.658063, 1.745329, 1.832596, 1.919862, 2.007129, 2.094395, 2.181662, 2.268928, 2.356194, 2.443461, 2.530727, 2.617994, 2.705260, 2.792527, 2.879793, 2.967060, 3.054326, 3.141593, 3.228859, 3.316126, 3.403392, 3.490659, 3.577925, 3.665191, 3.752458, 3.839724, 3.926991, 4.014257, 4.101524, 4.188790, 4.276057, 4.363323, 4.450590, 4.537856, 4.625123, 4.712389\n" // this binning starts at -pi/2 and is modulo 3 "delta_eta: -2.4, -2.3, -2.2, -2.1, -2.0, -1.9, -1.8, -1.7, -1.6, -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, 1.6, 1.7, 1.8, 1.9, 2.0,2.1, 2.2, 2.3, 2.4\n" "multiplicity: 0, 1, 2, 3, 4, 5, 10, 20, 30, 40, 50, 60, 70, 80, 90, 100.1\n" - "multiplicity_mixing: 0., 1., 2., 3., 4., 5., 10., 20., 30., 40., 50., 60., 70., 80., 90., 100.1\n" - "InvariantMass:0.200,0.300,0.398,0.399,0.4,0.401,0.402,0.403,0.404,0.405,0.406,0.407,0.408,0.409,0.41,0.411,0.412,0.413,0.414,0.415,0.416,0.417,0.418,0.419,0.42,0.421,0.422,0.423,0.424,0.425,0.426,0.427,0.428,0.429,0.43,0.431,0.432,0.433,0.434,0.435,0.436,0.437,0.438,0.439,0.44,0.441,0.442,0.443,0.444,0.445,0.446,0.447,0.448,0.449,0.45,0.451,0.452,0.453,0.454,0.455,0.456,0.457,0.458,0.459,0.46,0.461,0.462,0.463,0.464,0.465,0.466,0.467,0.468,0.469,0.47,0.471,0.472,0.473,0.474,0.475,0.476,0.477,0.478,0.479,0.48,0.481,0.482,0.483,0.484,0.485,0.486,0.487,0.488,0.489,0.49,0.491,0.492,0.493,0.494,0.495,0.496,0.497,0.498,0.499,0.5,0.501,0.502,0.503,0.504,0.505,0.506,0.507,0.508,0.509,0.51,0.511,0.512,0.513,0.514,0.515,0.516,0.517,0.518,0.519,0.52,0.521,0.522,0.523,0.524,0.525,0.526,0.527,0.528,0.529,0.53,0.531,0.532,0.533,0.534,0.535,0.536,0.537,0.538,0.539,0.54,0.541,0.542,0.543,0.544,0.545,0.546,0.547,0.548,0.549,0.55,0.551,0.552,0.553,0.554,0.555,0.556,0.557,0.558,0.559,0.56,0.561,0.562,0.563,0.564,0.565,0.566,0.567,0.568,0.569,0.57,0.571,0.572,0.573,0.574,0.575,0.576,0.577,0.578,0.579,0.58,0.581,0.582,0.583,0.584,0.585,0.586,0.587,0.588,0.589,0.59,0.591,0.592,0.593,0.594,0.595,0.596,0.597,0.598,0.599,0.600,0.700,0.800,0.900,1.000,1.065,1.066,1.067,1.068,1.069,1.07,1.071,1.072,1.073,1.074,1.075,1.076,1.077,1.078,1.079,1.08,1.081,1.082,1.083,1.084,1.085,1.086,1.087,1.088,1.089,1.09,1.091,1.092,1.093,1.094,1.095,1.096,1.097,1.098,1.099,1.1,1.101,1.102,1.103,1.104,1.105,1.106,1.107,1.108,1.109,1.11,1.111,1.112,1.113,1.114,1.115,1.116,1.117,1.118,1.119,1.12,1.121,1.122,1.123,1.124,1.125,1.126,1.127,1.128,1.129,1.13,1.131,1.132,1.133,1.134,1.135,1.136,1.137,1.138,1.139,1.14,1.141,1.142,1.143,1.144,1.145,1.146,1.147,1.148,1.149,1.15,1.151,1.152,1.153,1.154,1.155,1.156,1.157,1.158,1.159,1.16,1.161,1.162,1.163,1.164,1.165\n"; + "multiplicity_mixing: 0., 1., 2., 3., 4., 5., 10., 20., 30., 40., 50., 60., 70., 80., 90., 100.1\n"; + + if(fV0TrigCorr){ +defaultBinningStr += "InvariantMass:0.200,0.300,0.398,0.399,0.4,0.401,0.402,0.403,0.404,0.405,0.406,0.407,0.408,0.409,0.41,0.411,0.412,0.413,0.414,0.415,0.416,0.417,0.418,0.419,0.42,0.421,0.422,0.423,0.424,0.425,0.426,0.427,0.428,0.429,0.43,0.431,0.432,0.433,0.434,0.435,0.436,0.437,0.438,0.439,0.44,0.441,0.442,0.443,0.444,0.445,0.446,0.447,0.448,0.449,0.45,0.451,0.452,0.453,0.454,0.455,0.456,0.457,0.458,0.459,0.46,0.461,0.462,0.463,0.464,0.465,0.466,0.467,0.468,0.469,0.47,0.471,0.472,0.473,0.474,0.475,0.476,0.477,0.478,0.479,0.48,0.481,0.482,0.483,0.484,0.485,0.486,0.487,0.488,0.489,0.49,0.491,0.492,0.493,0.494,0.495,0.496,0.497,0.498,0.499,0.5,0.501,0.502,0.503,0.504,0.505,0.506,0.507,0.508,0.509,0.51,0.511,0.512,0.513,0.514,0.515,0.516,0.517,0.518,0.519,0.52,0.521,0.522,0.523,0.524,0.525,0.526,0.527,0.528,0.529,0.53,0.531,0.532,0.533,0.534,0.535,0.536,0.537,0.538,0.539,0.54,0.541,0.542,0.543,0.544,0.545,0.546,0.547,0.548,0.549,0.55,0.551,0.552,0.553,0.554,0.555,0.556,0.557,0.558,0.559,0.56,0.561,0.562,0.563,0.564,0.565,0.566,0.567,0.568,0.569,0.57,0.571,0.572,0.573,0.574,0.575,0.576,0.577,0.578,0.579,0.58,0.581,0.582,0.583,0.584,0.585,0.586,0.587,0.588,0.589,0.59,0.591,0.592,0.593,0.594,0.595,0.596,0.597,0.598,0.599,0.600,0.700,0.800,0.900,1.000,1.065,1.066,1.067,1.068,1.069,1.07,1.071,1.072,1.073,1.074,1.075,1.076,1.077,1.078,1.079,1.08,1.081,1.082,1.083,1.084,1.085,1.086,1.087,1.088,1.089,1.09,1.091,1.092,1.093,1.094,1.095,1.096,1.097,1.098,1.099,1.1,1.101,1.102,1.103,1.104,1.105,1.106,1.107,1.108,1.109,1.11,1.111,1.112,1.113,1.114,1.115,1.116,1.117,1.118,1.119,1.12,1.121,1.122,1.123,1.124,1.125,1.126,1.127,1.128,1.129,1.13,1.131,1.132,1.133,1.134,1.135,1.136,1.137,1.138,1.139,1.14,1.141,1.142,1.143,1.144,1.145,1.146,1.147,1.148,1.149,1.15,1.151,1.152,1.153,1.154,1.155,1.156,1.157,1.158,1.159,1.16,1.161,1.162,1.163,1.164,1.165\n"; + } if(fRequestEventPlane){ defaultBinningStr += "eventPlane: -0.5,0.5,1.5,2.5,3.5\n"; // Event Plane Bins (Psi: -0.5->0.5 (in plane), 0.5->1.5 (intermediate), 1.5->2.5 (out of plane), 2.5->3.5 (rest)) } - + if(fRequestEventPlanemixing){ +defaultBinningStr += "eventPlanemixing: 0.0*TMath::DegToRad(), 30.0*TMath::DegToRad(), 60.0*TMath::DegToRad(), 90.0*TMath::DegToRad(), 120.0*TMath::DegToRad(),150.0*TMath::DegToRad(),180.1*TMath::DegToRad()\n"; + } if(fcontainPIDtrig){ defaultBinningStr += "PIDTrig: -0.5,0.5,1.5,2.5,3.5,4.5,5.5,6.5,7.5,8.5,9.5,10.5\n"; // course } @@ -1104,17 +1114,16 @@ Double_t ZvrtxBins[NofVrtxBins+1]={ -10, -8, -6, -4, -2, 0, 2, 4, 6 if(fSampleType=="pp_2_76" || fCentralityMethod.EndsWith("_MANUAL") || (fSampleType=="pp_7" && fPPVsMultUtils==kFALSE))//mainly Tracks manual method { -if(fRequestEventPlane){ - // Event plane angle (Psi) bins - /* - Double_t* psibins = NULL; - Int_t nPsiBins; - psibins = GetBinning(fBinningString, "eventPlane", nPsiBins); - */ - const Int_t nPsiBins=6; - Double_t psibins[nPsiBins+1]={0.0*TMath::DegToRad(), 30.0*TMath::DegToRad(), 60.0*TMath::DegToRad(), 90.0*TMath::DegToRad(), 120.0*TMath::DegToRad(),150.0*TMath::DegToRad(),180.1*TMath::DegToRad()}; +if(fRequestEventPlanemixing){ + // Event plane angle (Psi) bins for event mixing + + Int_t nPsiBins=-1;; + Double_t* psibins = GetBinning(fBinningString, "eventPlanemixing", nPsiBins); + + //const Int_t nPsiBins=6; + //Double_t psibins[nPsiBins+1]={0.0*TMath::DegToRad(), 30.0*TMath::DegToRad(), 60.0*TMath::DegToRad(), 90.0*TMath::DegToRad(), 120.0*TMath::DegToRad(),150.0*TMath::DegToRad(),180.1*TMath::DegToRad()}; fPoolMgr = new AliEventPoolManager(MaxNofEvents,fMaxNofMixingTracks,multmixbin,multmix,NofVrtxBins,ZvrtxBins, nPsiBins, psibins); -// if(psibins) delete [] psibins; + if(psibins) delete [] psibins; } else{ @@ -1122,7 +1131,6 @@ fPoolMgr = new AliEventPoolManager(MaxNofEvents,fMaxNofMixingTracks,multmixbin,m Double_t psibinsd[nPsiBinsd+1]={0.0, 2000.0}; fPoolMgr = new AliEventPoolManager(MaxNofEvents,fMaxNofMixingTracks,multmixbin,multmix,NofVrtxBins,ZvrtxBins, nPsiBinsd, psibinsd); -// fPoolMgr = new AliEventPoolManager(MaxNofEvents,fMaxNofMixingTracks,NofCentBins,CentralityBins,NofVrtxBins,ZvrtxBins); } fPoolMgr->SetTargetValues(fMaxNofMixingTracks, 0.1, 5); @@ -1130,17 +1138,15 @@ fPoolMgr->SetTargetValues(fMaxNofMixingTracks, 0.1, 5); else//mainle centrality or quantile or Impactparameter method { - if(fRequestEventPlane){ - // Event plane angle (Psi) bins - /* - Double_t* psibins = NULL; - Int_t nPsiBins; - psibins = GetBinning(fBinningString, "eventPlane", nPsiBins); - */ - const Int_t nPsiBins=6; - Double_t psibins[nPsiBins+1]={0.0*TMath::DegToRad(), 30.0*TMath::DegToRad(), 60.0*TMath::DegToRad(), 90.0*TMath::DegToRad(), 120.0*TMath::DegToRad(),150.0*TMath::DegToRad(),180.1*TMath::DegToRad()}; + if(fRequestEventPlanemixing){ + // Event plane angle (Psi) bins for mixing + Int_t nPsiBins=-1;; + Double_t* psibins = GetBinning(fBinningString, "eventPlanemixing", nPsiBins); + + //const Int_t nPsiBins=6; + //Double_t psibins[nPsiBins+1]={0.0*TMath::DegToRad(), 30.0*TMath::DegToRad(), 60.0*TMath::DegToRad(), 90.0*TMath::DegToRad(), 120.0*TMath::DegToRad(),150.0*TMath::DegToRad(),180.1*TMath::DegToRad()}; fPoolMgr = new AliEventPoolManager(MaxNofEvents,fMaxNofMixingTracks,multmixbin,multmix,NofVrtxBins,ZvrtxBins, nPsiBins, psibins); -// if(psibins) delete [] psibins; + if(psibins) delete [] psibins; } else{ const Int_t nPsiBinsd=1; @@ -1334,7 +1340,7 @@ axisTitleTrig[dim_val_trig+1]=axisTitlePair[dim_val+2]; fOutput->Add(fTHnTrigcount); } - if((fAnalysisType =="MCAOD") && ffilltrigIDassoIDMCTRUTH) { + if((fAnalysisType =="MCAOD" || fAnalysisType =="MC") && ffilltrigIDassoIDMCTRUTH) { //AliTHns for trigger counting(truth MC) fTHnTrigcountMCTruthPrim = new AliTHn("fTHnTrigcountMCTruthPrim", "fTHnTrigcountMCTruthPrim", 2, dims, fBinst); //2 steps;;;;0->same event;;;;;1->mixed event for(Int_t i=0; iAdd(fTHnTrigcountMCTruthPrim); } -if(fAnalysisType=="MCAOD"){ +if(fAnalysisType=="MCAOD" || fAnalysisType=="MC"){ if(ffillhistQATruth) { MCtruthpt=new TH1F ("MCtruthpt","ptdistributiontruthprim",100,0.,10.); @@ -1842,13 +1848,13 @@ skipParticlesAbove = eventHeader->NProduced(); if (fSampleType=="pp_2_76" || fCentralityMethod.EndsWith("_MANUAL") || (fSampleType=="pp_7" && fPPVsMultUtils==kFALSE)) { //make the event selection with reco vertex cut and centrality cut and return the value of the centrality - Double_t refmultTruth = GetAcceptedEventMultiplicity(aod,kTRUE); //incase of ref multiplicity it will return the truth MC ref mullt value; need to determine the ref mult value separately for reco Mc case; in case of centrality this is final and fine - refmultReco = GetAcceptedEventMultiplicity(aod,kFALSE); + Double_t refmultTruth = GetAcceptedEventMultiplicity((AliVEvent*)aod,kTRUE); //incase of ref multiplicity it will return the truth MC ref mullt value; need to determine the ref mult value separately for reco Mc case; in case of centrality this is final and fine + refmultReco = GetAcceptedEventMultiplicity((AliVEvent*)aod,kFALSE); if(refmultTruth<=0 || refmultReco<=0) return; cent_v0=refmultTruth; } else { - cent_v0=GetAcceptedEventMultiplicity(aod,kFALSE); //centrality value; 2nd argument has no meaning + cent_v0=GetAcceptedEventMultiplicity((AliVEvent*)aod,kFALSE); //centrality value; 2nd argument has no meaning } if(cent_v0<0.) return; @@ -1856,7 +1862,7 @@ skipParticlesAbove = eventHeader->NProduced(); //get the event plane in case of PbPb if(fRequestEventPlane){ - gReactionPlane=GetEventPlane(aod,kTRUE,cent_v0);//get the truth event plane + gReactionPlane=GetEventPlane((AliVEvent*)aod,kTRUE,cent_v0);//get the truth event plane if(gReactionPlane==999.) return; } @@ -1937,6 +1943,7 @@ Int_t particletypeTruth=-999; particletypeTruth=SpPion; if(ffillhistQATruth) { + MCtruthpionpt->Fill(partMC->Pt()); MCtruthpioneta->Fill(partMC->Eta()); MCtruthpionphi->Fill(partMC->Phi()); @@ -2060,7 +2067,7 @@ if(tracksMCtruth) delete tracksMCtruth; effcent=cent_v0;// This will be required for efficiency THn filling(specially in case of pp) if(fRequestEventPlane){ - gReactionPlane = GetEventPlane(aod,kFALSE,cent_v0);//get the reconstructed event plane + gReactionPlane = GetEventPlane((AliVEvent*)aod,kFALSE,cent_v0);//get the reconstructed event plane if(gReactionPlane==999.) return; } @@ -2526,15 +2533,16 @@ if(fV0TrigCorr) { TArrayF gVertexArray; header->PrimaryVertex(gVertexArray); Float_t zVtxmc =gVertexArray.At(2); + //cout<<"*****************************************************************************************************hi I am here"<GetNumberOfPrimaries(); iTracks++) { if(partMC->Charge() == 0) continue; - //give only kinematic cuts at the generator level if (partMC->Eta() < fmineta || partMC->Eta() > fmaxeta) continue; if (partMC->Pt() < fminPt || partMC->Pt() > fmaxPt) continue; @@ -2768,13 +2775,13 @@ if (!fPID) return;//this should be available with each event even if we don't do // check event cuts and fill event histograms and return the centrality or reference multiplicity value - if((cent_v0 = GetAcceptedEventMultiplicity(aod,kFALSE)) < 0){ + if((cent_v0 = GetAcceptedEventMultiplicity((AliVEvent*)aod,kFALSE)) < 0){ return; } effcent=cent_v0;//required for efficiency correction case********Extremely Important //get the event plane in case of PbPb if(fRequestEventPlane){ - gReactionPlane = GetEventPlane(aod,kFALSE,cent_v0); + gReactionPlane = GetEventPlane((AliVEvent*)aod,kFALSE,cent_v0); if(gReactionPlane==999.) return; } @@ -4821,7 +4828,7 @@ if(fSampleType=="pp_7" && fPPVsMultUtils) else if(fCentralityMethod=="V0M_MANUAL" || fCentralityMethod=="V0A_MANUAL" || fCentralityMethod=="V0C_MANUAL" || fCentralityMethod=="TRACKS_MANUAL" || shift_to_TRACKS_MANUAL)//data or RecoMc and also for TRUTH { if(!truth){//for data or RecoMC - cent_v0 = GetReferenceMultiplicityVZEROFromAOD(event); + cent_v0 = GetReferenceMultiplicityVZEROFromAOD((AliVEvent*)event); }//for data or RecoMC if(truth && (fAnalysisType == "MCAOD")){//condition for TRUTH case @@ -5030,7 +5037,7 @@ fHistQA[0]->Fill((gVertexArray.At(0)));fHistQA[1]->Fill((gVertexArray.At(1)));fH fHistQA[3]->Fill((gVertexArray.At(0)));fHistQA[4]->Fill((gVertexArray.At(1)));fHistQA[5]->Fill((gVertexArray.At(2)));//after vertex cut,for trkVtx only // get the reference multiplicty or centrality - gRefMultiplicity = GetRefMultiOrCentrality(event,kFALSE);//2nd argument has no meaning + gRefMultiplicity = GetRefMultiOrCentrality((AliVEvent*)mcevent,kFALSE);//2nd argument has no meaning if(gRefMultiplicity<0) return -1; @@ -5138,9 +5145,9 @@ fHistQA[0]->Fill((trkVtx->GetX()));fHistQA[1]->Fill((trkVtx->GetY()));fHistQA[2] fHistQA[3]->Fill((trkVtx->GetX()));fHistQA[4]->Fill((trkVtx->GetY()));fHistQA[5]->Fill((trkVtx->GetZ()));//after vertex cut,for trkVtx only //get the centrality or multiplicity - if(truth) {gRefMultiplicity = GetRefMultiOrCentrality(aod,kTRUE);}//kTRUE-->for Truth case(only meaningful in case of ref multiplicity,in case of centrality it has no meaning) + if(truth) {gRefMultiplicity = GetRefMultiOrCentrality((AliVEvent*)aod,kTRUE);}//kTRUE-->for Truth case(only meaningful in case of ref multiplicity,in case of centrality it has no meaning) - else {gRefMultiplicity = GetRefMultiOrCentrality(aod,kFALSE);}//kFALSE-->for data and RecoMc case(only meaningful in case of ref multiplicity,in case of centrality it has no meaning) + else {gRefMultiplicity = GetRefMultiOrCentrality((AliVEvent*)aod,kFALSE);}//kFALSE-->for data and RecoMc case(only meaningful in case of ref multiplicity,in case of centrality it has no meaning) if(gRefMultiplicity<0) return -1; diff --git a/PWGCF/Correlations/DPhi/TriggerPID/AliTwoParticlePIDCorr.h b/PWGCF/Correlations/DPhi/TriggerPID/AliTwoParticlePIDCorr.h index 5061c82b85c..78b6cbc1b17 100644 --- a/PWGCF/Correlations/DPhi/TriggerPID/AliTwoParticlePIDCorr.h +++ b/PWGCF/Correlations/DPhi/TriggerPID/AliTwoParticlePIDCorr.h @@ -135,8 +135,9 @@ fVzMax_MC=VzMax_MC; void SetCentralityEstimator(TString CentralityMethod) { fCentralityMethod = CentralityMethod;} void SetPPVsMultUtils(Bool_t val) {fPPVsMultUtils = val;} void SetSampleType(TString SampleType) {fSampleType=SampleType;} - void SetRequestEventPlane(Bool_t RequestEventPlane,Bool_t V2,Bool_t V3,TString EPdetector,Bool_t IsAfter2011){ + void SetRequestEventPlane(Bool_t RequestEventPlane,Bool_t RequestEventPlanemixing,Bool_t V2,Bool_t V3,TString EPdetector,Bool_t IsAfter2011){ fRequestEventPlane=RequestEventPlane; +fRequestEventPlanemixing=RequestEventPlanemixing; fV2=V2; fV3=V3; fEPdet=EPdetector; @@ -296,6 +297,7 @@ fFracTPCcls=FracSharedTPCcls;//0.4 Bool_t fPPVsMultUtils;//switch to ON quantile information for pp 7 TeV case TString fSampleType; // pp,p-Pb,Pb-Pb Bool_t fRequestEventPlane; //only for PbPb + Bool_t fRequestEventPlanemixing; //only for PbPb Int_t fnTracksVertex; // QA tracks pointing to principal vertex AliAODVertex* trkVtx;//! Float_t zvtx; @@ -704,3 +706,4 @@ fCentralityCorrelation //ParticlePID_InvMass //particlepidtrig +//fRequestEventPlanemixing