fCorrelEventplaneMCDATA(0),
//global
fIsMonteCarlo(0),
+fEPselector("Q"),
// event cut variables
fCutMaxZVertex(10.),
// track kinematic cut variables
dMCEventZv = mcHdr->GetVtxZ();
dMCTrackZvPtEtaCent[0] = dMCEventZv;
- dMCEventplaneAngle = MoveMCEventplane(genHijingHeader->ReactionPlaneAngle());
+ dMCEventplaneAngle = MoveEventplane(genHijingHeader->ReactionPlaneAngle());
fEventStatistics->Fill("MC all events",1);
fMCEventplaneDist->Fill(dMCEventplaneAngle);
}
// get event plane Angle from AODHeader, default is Q
ep = const_cast<AliAODEvent*>(eventAOD)->GetEventplane();
if(ep) {
- dEventplaneAngle = ep->GetEventplane("V0",eventAOD);
+ dEventplaneAngle = MoveEventplane(ep->GetEventplane(GetEventplaneSelector().Data(),eventAOD));
}
// cout << dEventplaneAngle << endl;
if(bIsEventSelectedSemi) fEventStatisticsCentralityTrigger->Fill(dCentrality, 1);
if(bIsEventSelectedCentral) fEventStatisticsCentralityTrigger->Fill(dCentrality, 2);
- Double_t dEventZvMultCent[3] = {dEventZv, iAcceptedMultiplicity, dCentrality};
+ Double_t dEventZvMultCent[3] = {dEventZv, static_cast<Double_t>(iAcceptedMultiplicity), dCentrality};
fZvMultCent->Fill(dEventZvMultCent);
// store correlation between data and MC eventplane
}
-Double_t AlidNdPtAnalysisPbPbAOD::MoveMCEventplane(Double_t dMCEP)
+Double_t AlidNdPtAnalysisPbPbAOD::MoveEventplane(Double_t dMCEP)
{
Double_t retval = 0;
retval = dMCEP;
fCutSettings->Fill("fCutMaxChi2TPCConstrainedGlobal", fCutMaxChi2TPCConstrainedGlobal);
if(fCutLengthInTPCPtDependent) fCutSettings->Fill("fCutLengthInTPCPtDependent", 1);
fCutSettings->Fill("fPrefactorLengthInTPCPtDependent", fPrefactorLengthInTPCPtDependent);
+ fCutSettings->Fill(Form("EP selector %s", fEPselector.Data()), 1);
}
Bool_t AlidNdPtAnalysisPbPbAOD::GetDCA(const AliAODTrack *track, AliAODEvent *evt, Double_t d0z0[2])