From: mfloris Date: Wed, 3 Sep 2014 20:03:47 +0000 (+0200) Subject: AliAODEvent::GetHeader() returns AliVHeader X-Git-Url: http://git.uio.no/git/?p=u%2Fmrichter%2FAliRoot.git;a=commitdiff_plain;h=0a918d8d80e9900008464c3bb616f57b0c154abf AliAODEvent::GetHeader() returns AliVHeader Transition completed, added dynamic or c-style casts where needed --- diff --git a/ANALYSIS/AliAnalysisUtils.cxx b/ANALYSIS/AliAnalysisUtils.cxx index 8890165bc7b..15fef1d74d7 100644 --- a/ANALYSIS/AliAnalysisUtils.cxx +++ b/ANALYSIS/AliAnalysisUtils.cxx @@ -114,7 +114,8 @@ Bool_t AliAnalysisUtils::IsFirstEventInChunk(AliVEvent *event) if(fisAOD){ AliAODHeader *aodheader = 0x0; - aodheader = aod->GetHeader(); + aodheader = dynamic_cast(aod->GetHeader()); + if(!aodheader) AliFatal("Not a standard AOD"); if(!aodheader){ AliFatal("AOD header not there ?!"); return kFALSE; @@ -216,11 +217,11 @@ Bool_t AliAnalysisUtils::IsOutOfBunchPileUp(AliVEvent *event) AliFatal("Event is neither of AOD nor ESD type"); return kFALSE; } - Int_t bc2 = (aod)?aod->GetHeader()->GetIRInt2ClosestInteractionMap():esd->GetHeader()->GetIRInt2ClosestInteractionMap(); + Int_t bc2 = (aod)?((AliVAODHeader*)aod->GetHeader())->GetIRInt2ClosestInteractionMap():esd->GetHeader()->GetIRInt2ClosestInteractionMap(); if (bc2 != 0) return kTRUE; - Int_t bc1 = (aod)?aod->GetHeader()->GetIRInt1ClosestInteractionMap():esd->GetHeader()->GetIRInt1ClosestInteractionMap(); + Int_t bc1 = (aod)?((AliVAODHeader*)aod->GetHeader())->GetIRInt1ClosestInteractionMap():esd->GetHeader()->GetIRInt1ClosestInteractionMap(); if (bc1 != 0) return kTRUE; diff --git a/ANALYSIS/AliCentralitySelectionTask.cxx b/ANALYSIS/AliCentralitySelectionTask.cxx index bb06834bdef..75c0e8f5132 100644 --- a/ANALYSIS/AliCentralitySelectionTask.cxx +++ b/ANALYSIS/AliCentralitySelectionTask.cxx @@ -1444,7 +1444,8 @@ void AliCentralitySelectionTask::UserExec(Option_t */*option*/) //nTracks = event->GetNumberOfTracks(); nTracks = fTrackCuts ? (Short_t)fTrackCuts->GetReferenceMultiplicity(esd,kTRUE):-1; } else { - AliAODHeader *h = aod->GetHeader(); + AliAODHeader *h = dynamic_cast(aod->GetHeader()); + if(!h) AliFatal("Not a standard AOD"); nTracks = h!=0 ? (Short_t)h->GetTPConlyRefMultiplicity():-1; } @@ -1493,7 +1494,8 @@ void AliCentralitySelectionTask::UserExec(Option_t */*option*/) } else { AliAODTracklets *mult = aod->GetTracklets(); nTracklets = mult->GetNumberOfTracklets(); - AliAODHeader *h = aod->GetHeader(); + AliAODHeader *h = dynamic_cast(aod->GetHeader()); + if(!h) AliFatal("Not a standard AOD"); for(Int_t ilay=0; ilay<6; ilay++){ nClusters[ilay] = h->GetNumberOfITSClusters(ilay); } @@ -1580,7 +1582,8 @@ void AliCentralitySelectionTask::UserExec(Option_t */*option*/) if (zpcFired) zpcTower = ZPCtower[0]; } else { - AliAODHeader *h = aod->GetHeader(); + AliAODHeader *h = dynamic_cast(aod->GetHeader()); + if(!h) AliFatal("Not a standard AOD"); zncEnergy = (Float_t) (h->GetZDCN1Energy()); zpcEnergy = (Float_t) (h->GetZDCP1Energy()); znaEnergy = (Float_t) (h->GetZDCN2Energy()); diff --git a/ANALYSIS/AliEPSelectionTask.cxx b/ANALYSIS/AliEPSelectionTask.cxx index 601035b61be..cf7bd3e4a0c 100644 --- a/ANALYSIS/AliEPSelectionTask.cxx +++ b/ANALYSIS/AliEPSelectionTask.cxx @@ -365,7 +365,9 @@ void AliEPSelectionTask::UserExec(Option_t */*option*/) if (headerH) fRP = headerH->GetReactionPlaneAngle(); } - esdEP = aod->GetHeader()->GetEventplaneP(); + AliAODHeader * header = dynamic_cast(aod->GetHeader()); + if(!header) AliFatal("Not a standard AOD"); + esdEP = header->GetEventplaneP(); if(!esdEP) return; // protection against missing EP branch (nanoAODs) esdEP->Reset(); diff --git a/ANALYSIS/ESDfilter/AliAnalysisTaskESDfilter.cxx b/ANALYSIS/ESDfilter/AliAnalysisTaskESDfilter.cxx index 24c0dd6ea2d..5e1d52d88f1 100644 --- a/ANALYSIS/ESDfilter/AliAnalysisTaskESDfilter.cxx +++ b/ANALYSIS/ESDfilter/AliAnalysisTaskESDfilter.cxx @@ -367,7 +367,8 @@ AliAODHeader* AliAnalysisTaskESDfilter::ConvertHeader(const AliESDEvent& esd) AliCodeTimerAuto("",0); - AliAODHeader* header = AODEvent()->GetHeader(); + AliAODHeader* header = dynamic_cast(AODEvent()->GetHeader()); + if(!header) AliFatal("Not a standard AOD"); header->SetRunNumber(esd.GetRunNumber()); header->SetOfflineTrigger(fInputHandler->IsEventSelected()); // propagate the decision of the physics selection diff --git a/MUON/AliMUONAlignmentTask.cxx b/MUON/AliMUONAlignmentTask.cxx index 08e96aac320..dc0961283f9 100644 --- a/MUON/AliMUONAlignmentTask.cxx +++ b/MUON/AliMUONAlignmentTask.cxx @@ -364,7 +364,8 @@ void AliMUONAlignmentTask::UserExec(Option_t *) if( handler ) { AliAODEvent* aod = handler->GetAOD(); - AliAODHeader* header = aod->GetHeader(); + AliAODHeader* header = dynamic_cast(aod->GetHeader()); + if(!header) AliFatal("Not a standard AOD"); header->SetRunNumber(lESD->GetRunNumber()); AliAnalysisManager::GetAnalysisManager()->GetOutputEventHandler()->SetFillAOD(kTRUE); diff --git a/PWG/DevNanoAOD/AliAnalysisNanoAODCuts.cxx b/PWG/DevNanoAOD/AliAnalysisNanoAODCuts.cxx index 72dd3c1ddbc..ae5b4e86191 100644 --- a/PWG/DevNanoAOD/AliAnalysisNanoAODCuts.cxx +++ b/PWG/DevNanoAOD/AliAnalysisNanoAODCuts.cxx @@ -51,9 +51,12 @@ Bool_t AliAnalysisNanoAODEventCuts::IsSelected(TObject* obj) void AliNanoAODSimpleSetter::SetNanoAODHeader(const AliAODEvent * event , AliNanoAODHeader * head ) { + + AliAODHeader * header = dynamic_cast(event->GetHeader()); + if (!header) AliFatal("Not a standard AOD"); // Set custom nano aod vars - Double_t centr = event->GetHeader()->GetCentralityP()->GetCentralityPercentile("V0M"); - Double_t magfield = event->GetHeader()->GetMagneticField(); + Double_t centr = header->GetCentralityP()->GetCentralityPercentile("V0M"); + Double_t magfield = header->GetMagneticField(); head->SetVar(0, centr); head->SetVar(1, magfield); diff --git a/PWG/EMCAL/AliAnalysisTaskEMCALClusterizeFast.cxx b/PWG/EMCAL/AliAnalysisTaskEMCALClusterizeFast.cxx index 1e304629a47..6a78c619aef 100644 --- a/PWG/EMCAL/AliAnalysisTaskEMCALClusterizeFast.cxx +++ b/PWG/EMCAL/AliAnalysisTaskEMCALClusterizeFast.cxx @@ -223,7 +223,7 @@ void AliAnalysisTaskEMCALClusterizeFast::UserExec(Option_t *) AliAnalysisManager *am = AliAnalysisManager::GetAnalysisManager(); offtrigger = ((AliInputEventHandler*)(am->GetInputEventHandler()))->IsEventSelected(); } else { - offtrigger = fAod->GetHeader()->GetOfflineTrigger(); + offtrigger = ((AliVAODHeader*)fAod->GetHeader())->GetOfflineTrigger(); } if (!MCEvent()) { @@ -829,7 +829,8 @@ void AliAnalysisTaskEMCALClusterizeFast::Init() if (fEsd) { gm = fEsd->GetEMCALMatrix(mod); } else { - AliAODHeader *aodheader = fAod->GetHeader(); + AliAODHeader *aodheader = dynamic_cast(fAod->GetHeader()); + if(!aodheader) AliFatal("Not a standard AOD"); if (aodheader) { gm = aodheader->GetEMCALMatrix(mod); } diff --git a/PWG/EMCAL/AliAnalysisTaskEmcal.cxx b/PWG/EMCAL/AliAnalysisTaskEmcal.cxx index 06ed828d5ae..46a9cbc0813 100644 --- a/PWG/EMCAL/AliAnalysisTaskEmcal.cxx +++ b/PWG/EMCAL/AliAnalysisTaskEmcal.cxx @@ -782,7 +782,7 @@ Bool_t AliAnalysisTaskEmcal::IsEventSelected() } else { const AliAODEvent *aev = dynamic_cast(InputEvent()); if (aev) { - res = aev->GetHeader()->GetOfflineTrigger(); + res = ((AliVAODHeader*)aev->GetHeader())->GetOfflineTrigger(); } } if ((res & fOffTrigger) == 0) { diff --git a/PWG/EMCAL/AliEmcalPhysicsSelection.cxx b/PWG/EMCAL/AliEmcalPhysicsSelection.cxx index d194f07cf47..a7e68ce5782 100644 --- a/PWG/EMCAL/AliEmcalPhysicsSelection.cxx +++ b/PWG/EMCAL/AliEmcalPhysicsSelection.cxx @@ -56,14 +56,14 @@ UInt_t AliEmcalPhysicsSelection::GetSelectionMask(const TObject* obj) am->LoadBranch("AliESDRun."); TString title(eev->GetHeader()->GetTitle()); if (1&&(title.Length()>0)) { - res = eev->GetHeader()->GetUniqueID(); + res = ((AliVAODHeader*)eev->GetHeader())->GetUniqueID(); res &= 0x4FFFFFFF; } else { res = IsCollisionCandidate(eev); } } else { aev = dynamic_cast(obj); - res = aev->GetHeader()->GetOfflineTrigger(); + res = ((AliVAODHeader*)aev->GetHeader())->GetOfflineTrigger(); } // return 0, if 0 found diff --git a/PWG/FLOW/Tasks/AliAnalysisTaskFlowEPCascade.cxx b/PWG/FLOW/Tasks/AliAnalysisTaskFlowEPCascade.cxx index 81d4263c112..fed211d7972 100644 --- a/PWG/FLOW/Tasks/AliAnalysisTaskFlowEPCascade.cxx +++ b/PWG/FLOW/Tasks/AliAnalysisTaskFlowEPCascade.cxx @@ -768,7 +768,8 @@ void AliAnalysisTaskFlowEPCascade::UserExec(Option_t *) //At the momment the cutting class does not handle AOD event properly //so we are doing the cuts explicitly here - AliAODHeader *aodHeader = fAOD->GetHeader(); + AliAODHeader *aodHeader = dynamic_cast(fAOD->GetHeader()); + if(!aodHeader) AliFatal("Not a standard AOD"); if(!aodHeader) return; AliCentrality *centrality = aodHeader->GetCentralityP(); if(!centrality) return; @@ -1433,7 +1434,7 @@ void AliAnalysisTaskFlowEPCascade::ReadFromESDv0(AliESDEvent *fESD){ void AliAnalysisTaskFlowEPCascade::ReadFromAODv0(AliAODEvent *fAOD){ - AliEventplane * ep = (fAOD->GetHeader())->GetEventplaneP(); + AliEventplane * ep = ((AliVAODHeader*)fAOD->GetHeader())->GetEventplaneP(); Double_t psiTPC = ep->GetEventplane("Q", fAOD, 2); // in range of [0, pi] // if(psiTPC > TMath::PiOver2()) // psiTPC -= TMath::Pi(); diff --git a/PWG/FLOW/Tasks/AliAnalysisTaskFlowEventforRP.cxx b/PWG/FLOW/Tasks/AliAnalysisTaskFlowEventforRP.cxx index 2c6ffccff5a..f294786714e 100644 --- a/PWG/FLOW/Tasks/AliAnalysisTaskFlowEventforRP.cxx +++ b/PWG/FLOW/Tasks/AliAnalysisTaskFlowEventforRP.cxx @@ -192,7 +192,8 @@ void AliAnalysisTaskFlowEventforRP::UserExec(Option_t *) // Update the header - AliAODHeader* header = AODEvent()->GetHeader(); + AliAODHeader* header = dynamic_cast(AODEvent()->GetHeader()); + if(!header) AliFatal("Not a standard AOD"); header->SetRunNumber(esd->GetRunNumber()); header->SetQTheta(dRP,1); diff --git a/PWG/FLOW/Tasks/AliAnalysisTaskFlowStrange.cxx b/PWG/FLOW/Tasks/AliAnalysisTaskFlowStrange.cxx index bc2e95e2bd6..e502a60d288 100644 --- a/PWG/FLOW/Tasks/AliAnalysisTaskFlowStrange.cxx +++ b/PWG/FLOW/Tasks/AliAnalysisTaskFlowStrange.cxx @@ -876,7 +876,7 @@ Bool_t AliAnalysisTaskFlowStrange::AcceptAAEvent(AliESDEvent*) { //======================================================================= Bool_t AliAnalysisTaskFlowStrange::MinimumRequirementsAA(AliAODEvent *tAOD) { fRunNumber = tAOD->GetRunNumber(); - AliCentrality *cent = tAOD->GetHeader()->GetCentralityP(); + AliCentrality *cent = ((AliVAODHeader*)tAOD->GetHeader())->GetCentralityP(); fV0M = cent->GetCentralityPercentile("V0M"); fTRK = cent->GetCentralityPercentile("TRK"); TString mycent = fCentMethod; @@ -989,9 +989,9 @@ Bool_t AliAnalysisTaskFlowStrange::AcceptPAEvent(AliAODEvent*) { // PA reading discontinued: TO BE UPDATED /* //if(aod->GetHeader()->GetEventNumberESDFile() == 0) return; //rejecting first chunk NOT NEEDED ANYMORE - Int_t bc2 = tAOD->GetHeader()->GetIRInt2ClosestInteractionMap(); + Int_t bc2 = ((AliVAODHeader*)tAOD->GetHeader())->GetIRInt2ClosestInteractionMap(); if(bc2!=0) return kFALSE; - Int_t bc1 = tAOD->GetHeader()->GetIRInt1ClosestInteractionMap(); + Int_t bc1 = ((AliVAODHeader*)tAOD->GetHeader())->GetIRInt1ClosestInteractionMap(); if(bc1!=0) return kFALSE; Short_t isPileup = tAOD->IsPileupFromSPD(5); if(isPileup!=0) return kFALSE; diff --git a/PWG/FLOW/Tasks/AliFlowEventCuts.cxx b/PWG/FLOW/Tasks/AliFlowEventCuts.cxx index e6bdcf55452..87b605451fe 100644 --- a/PWG/FLOW/Tasks/AliFlowEventCuts.cxx +++ b/PWG/FLOW/Tasks/AliFlowEventCuts.cxx @@ -426,7 +426,7 @@ Bool_t AliFlowEventCuts::PassesCuts(AliVEvent *event, AliMCEvent *mcevent) Double_t tSPDVtxZ = aodevent->GetPrimaryVertexSPD()->GetZ(); if( TMath::Abs(tVtxZ-tSPDVtxZ) > 0.5 ) pass = kFALSE; } - AliCentrality* centr = aodevent->GetHeader()->GetCentralityP(); + AliCentrality* centr = ((AliVAODHeader*)aodevent->GetHeader())->GetCentralityP(); if(fCutTPCmultiplicityOutliers || fCutTPCmultiplicityOutliersAOD){ Double_t v0Centr = centr->GetCentralityPercentile("V0M"); Double_t trkCentr = centr->GetCentralityPercentile("TRK"); @@ -503,7 +503,7 @@ Float_t AliFlowEventCuts::GetCentrality(AliVEvent* event, AliMCEvent* /*mcEvent* if (esdEvent) centr = esdEvent->GetCentrality(); if (aodEvent) - centr = aodEvent->GetHeader()->GetCentralityP(); + centr = ((AliVAODHeader*)aodEvent->GetHeader())->GetCentralityP(); if (!centr) return -1.; diff --git a/PWG/muon/AliAODMuonReplicator.cxx b/PWG/muon/AliAODMuonReplicator.cxx index 8f7a1bf6c5b..548f098e770 100644 --- a/PWG/muon/AliAODMuonReplicator.cxx +++ b/PWG/muon/AliAODMuonReplicator.cxx @@ -390,7 +390,9 @@ void AliAODMuonReplicator::ReplicateAndFilter(const AliAODEvent& source) if (fReplicateHeader) { - *fHeader = *(source.GetHeader()); + AliAODHeader * header = dynamic_cast(source.GetHeader()); + if(!header) AliFatal("Not a standard AOD"); + *fHeader = *(header); } if (fReplicateTracklets) diff --git a/PWG/muon/AliAnalysisMuonUtility.cxx b/PWG/muon/AliAnalysisMuonUtility.cxx index 8a9eb4dc401..1b039ce3e40 100644 --- a/PWG/muon/AliAnalysisMuonUtility.cxx +++ b/PWG/muon/AliAnalysisMuonUtility.cxx @@ -51,6 +51,7 @@ // CORRFW includes #include "AliCFGridSparse.h" +#include "assert.h" /// \cond CLASSIMP ClassImp(AliAnalysisMuonUtility) // Class implementation in ROOT context @@ -443,7 +444,9 @@ TString AliAnalysisMuonUtility::GetPassName ( const AliInputEventHandler* eventH // try first to find the info in the AOD header // (which is a priori safer because it works even on local copies of AODs) // and if it does not work, directly check the path to the AOD - filePath = static_cast (event)->GetHeader()->GetESDFileName(); + AliAODHeader * header = dynamic_cast(static_cast (event)->GetHeader()); + assert(header && "Not a standard AOD"); + filePath = header->GetESDFileName(); TString passName = GetPassName(filePath.Data()); if ( passName.IsNull() ) AliWarningClass("Check again with the AOD path"); else return passName; diff --git a/PWG/muon/AliAnalysisTaskCreateMixedDimuons.cxx b/PWG/muon/AliAnalysisTaskCreateMixedDimuons.cxx index c7bd3fa500b..d61d8ec23d7 100644 --- a/PWG/muon/AliAnalysisTaskCreateMixedDimuons.cxx +++ b/PWG/muon/AliAnalysisTaskCreateMixedDimuons.cxx @@ -197,9 +197,12 @@ void AliAnalysisTaskCreateMixedDimuons::UserExec(Option_t *) { } } - fOutputUserAOD->GetHeader()->SetRefMultiplicity(nFWMUonsAdded); - fOutputUserAOD->GetHeader()->SetRefMultiplicityPos(nPosTracksAdded); - fOutputUserAOD->GetHeader()->SetRefMultiplicityNeg(nNegTracksAdded); + + AliAODHeader * header = dynamic_cast(fOutputUserAOD->GetHeader()); + if(!header) AliFatal("Not a standard AOD"); + header->SetRefMultiplicity(nFWMUonsAdded); + header->SetRefMultiplicityPos(nPosTracksAdded); + header->SetRefMultiplicityNeg(nNegTracksAdded); fOutputUserHandler -> FinishEvent(); diff --git a/PWG/muon/AliAnalysisTaskESDMuonFilter.cxx b/PWG/muon/AliAnalysisTaskESDMuonFilter.cxx index df66140f82a..d4ee6e9fcd0 100644 --- a/PWG/muon/AliAnalysisTaskESDMuonFilter.cxx +++ b/PWG/muon/AliAnalysisTaskESDMuonFilter.cxx @@ -238,7 +238,8 @@ void AliAnalysisTaskESDMuonFilter::ConvertESDtoAOD() // for (Int_t i = 0; i < 10; pid[i++] = 0.) {} // pid[AliAODTrack::kMuon]=1.; - AliAODHeader* header = AODEvent()->GetHeader(); + AliAODHeader* header = dynamic_cast(AODEvent()->GetHeader()); + if(!header) AliFatal("Not a standard AOD"); AliAODTrack *aodTrack = 0x0; AliESDMuonTrack *esdMuTrack = 0x0; diff --git a/PWG/muon/AliAnalysisTaskLinkToMC.cxx b/PWG/muon/AliAnalysisTaskLinkToMC.cxx index dead2524ca5..49f57ce9f45 100644 --- a/PWG/muon/AliAnalysisTaskLinkToMC.cxx +++ b/PWG/muon/AliAnalysisTaskLinkToMC.cxx @@ -589,7 +589,8 @@ void AliAnalysisTaskLinkToMC::CreateAODTracks(TMap& links) // for (Int_t i = 0; i < 10; pid[i++] = 0.) {} // pid[AliAODTrack::kMuon]=1.; - AliAODHeader* header = AODEvent()->GetHeader(); + AliAODHeader* header = dynamic_cast(AODEvent()->GetHeader()); + if(!header) AliFatal("Not a standard AOD"); AliAODTrack *aodTrack = 0x0; AliESDMuonTrack *esdMuTrack = 0x0; diff --git a/PWG/muon/AliAnalysisTaskMuonAODfromGeneral.cxx b/PWG/muon/AliAnalysisTaskMuonAODfromGeneral.cxx index f28982b8193..0663a67a0a5 100644 --- a/PWG/muon/AliAnalysisTaskMuonAODfromGeneral.cxx +++ b/PWG/muon/AliAnalysisTaskMuonAODfromGeneral.cxx @@ -112,7 +112,8 @@ void AliAnalysisTaskMuonAODfromGeneral::Exec(Option_t *) { Double_t covVtx[6]; // Access to the header - AliAODHeader *header = fNewAOD->GetHeader(); + AliAODHeader *header = dynamic_cast(fNewAOD->GetHeader()); + if(!header) AliFatal("Not a standard AOD"); // fill the header header->SetRunNumber (fOrgAOD->GetRunNumber() ); diff --git a/PWGCF/Correlations/Base/AliAnalyseLeadingTrackUE.cxx b/PWGCF/Correlations/Base/AliAnalyseLeadingTrackUE.cxx index 9e55828df91..d156f4e14aa 100644 --- a/PWGCF/Correlations/Base/AliAnalyseLeadingTrackUE.cxx +++ b/PWGCF/Correlations/Base/AliAnalyseLeadingTrackUE.cxx @@ -725,7 +725,7 @@ AliVParticle* AliAnalyseLeadingTrackUE::ParticleWithCuts(TObject* obj, Int_t ip Double_t pos[2]; Double_t covar[3]; AliAODTrack* clone = (AliAODTrack*) part->Clone(); - Bool_t success = clone->PropagateToDCA(vertex, aodEvent->GetHeader()->GetMagneticField(), 3, pos, covar); + Bool_t success = clone->PropagateToDCA(vertex, ((AliVAODHeader*)aodEvent->GetHeader())->GetMagneticField(), 3, pos, covar); delete clone; if (!success) return 0; diff --git a/PWGCF/Correlations/DPhi/AliAnalysisTask3PCorrelations.cxx b/PWGCF/Correlations/DPhi/AliAnalysisTask3PCorrelations.cxx index e5676b79dbd..54111c2ea3e 100644 --- a/PWGCF/Correlations/DPhi/AliAnalysisTask3PCorrelations.cxx +++ b/PWGCF/Correlations/DPhi/AliAnalysisTask3PCorrelations.cxx @@ -1054,7 +1054,7 @@ void AliAnalysisTask3PCorrelations::UserExec(Option_t */*option*/) _eventAccounting->Fill(1);// count all calls to this function with a valid pointer //Centrality - AliCentrality* centralityObject = fAODEvent->GetHeader()->GetCentralityP(); + AliCentrality* centralityObject = ((AliVAODHeader*)fAODEvent->GetHeader())->GetCentralityP(); if (centralityObject) { // cout << "AliAnalysisTask3PCorrelations::UserExec(Option_t *option) - 6" << endl; diff --git a/PWGCF/Correlations/DPhi/AliAnalysisTaskContMC.cxx b/PWGCF/Correlations/DPhi/AliAnalysisTaskContMC.cxx index b3174553534..bfe2f953f5c 100644 --- a/PWGCF/Correlations/DPhi/AliAnalysisTaskContMC.cxx +++ b/PWGCF/Correlations/DPhi/AliAnalysisTaskContMC.cxx @@ -113,7 +113,7 @@ void AliAnalysisTaskContMC::UserExec(Option_t *) } Double_t centrality = 0; - AliCentrality *centralityObj = fAOD->GetHeader()->GetCentralityP(); + AliCentrality *centralityObj = ((AliVAODHeader*)fAOD->GetHeader())->GetCentralityP(); if (centralityObj) { centrality = centralityObj->GetCentralityPercentileUnchecked("V0M"); diff --git a/PWGCF/Correlations/DPhi/AliAnalysisTaskDiHadron.cxx b/PWGCF/Correlations/DPhi/AliAnalysisTaskDiHadron.cxx index 4b9ba1b557e..ac95e4852a3 100644 --- a/PWGCF/Correlations/DPhi/AliAnalysisTaskDiHadron.cxx +++ b/PWGCF/Correlations/DPhi/AliAnalysisTaskDiHadron.cxx @@ -2130,7 +2130,8 @@ void AliAnalysisTaskDiHadron::Exec(Option_t *) else tMult=fESD->GetVZEROData()->GetMTotV0A()+fESD->GetVZEROData()->GetMTotV0C(); } else{ - AliAODHeader *tHeader=fAOD->GetHeader(); + AliAODHeader *tHeader=dynamic_cast(fAOD->GetHeader()); + if(!tHeader) AliFatal("Not a standard AOD"); tMult=tHeader->GetCentrality(); } diff --git a/PWGCF/Correlations/DPhi/AliAnalysisTaskDptDptCorrelations.cxx b/PWGCF/Correlations/DPhi/AliAnalysisTaskDptDptCorrelations.cxx index 05868835a4d..7a286ef1e66 100644 --- a/PWGCF/Correlations/DPhi/AliAnalysisTaskDptDptCorrelations.cxx +++ b/PWGCF/Correlations/DPhi/AliAnalysisTaskDptDptCorrelations.cxx @@ -1087,7 +1087,7 @@ void AliAnalysisTaskDptDptCorrelations::UserExec(Option_t */*option*/) if(fAODEvent) { //Centrality - AliCentrality* centralityObject = fAODEvent->GetHeader()->GetCentralityP(); + AliCentrality* centralityObject = ((AliVAODHeader*)fAODEvent->GetHeader())->GetCentralityP(); if (centralityObject) { //cout << "AliAnalysisTaskDptDptCorrelations::UserExec(Option_t *option) - 6" << endl; diff --git a/PWGCF/Correlations/DPhi/AliAnalysisTaskDptDptQA.cxx b/PWGCF/Correlations/DPhi/AliAnalysisTaskDptDptQA.cxx index 784e45356a1..914233cf83c 100644 --- a/PWGCF/Correlations/DPhi/AliAnalysisTaskDptDptQA.cxx +++ b/PWGCF/Correlations/DPhi/AliAnalysisTaskDptDptQA.cxx @@ -1053,7 +1053,7 @@ void AliAnalysisTaskDptDptQA::UserExec(Option_t */*option*/) //cout << "AliAnalysisTaskDptDptQA::UserExec(Option_t *option) - 5" << endl; //Centrality - AliCentrality* centralityObject = fAODEvent->GetHeader()->GetCentralityP(); + AliCentrality* centralityObject = ((AliVAODHeader*)fAODEvent->GetHeader())->GetCentralityP(); if (centralityObject) { //cout << "AliAnalysisTaskDptDptQA::UserExec(Option_t *option) - 6" << endl; diff --git a/PWGCF/Correlations/DPhi/AliAnalysisTaskLongRangeCorrelations.cxx b/PWGCF/Correlations/DPhi/AliAnalysisTaskLongRangeCorrelations.cxx index bf875b67189..4d26c93fc11 100644 --- a/PWGCF/Correlations/DPhi/AliAnalysisTaskLongRangeCorrelations.cxx +++ b/PWGCF/Correlations/DPhi/AliAnalysisTaskLongRangeCorrelations.cxx @@ -156,7 +156,8 @@ void AliAnalysisTaskLongRangeCorrelations::UserExec(Option_t* ) { AliAODEvent* pAOD(dynamic_cast(InputEvent())); if (NULL == pAOD) return; - AliAODHeader *pAODHeader = pAOD->GetHeader(); + AliAODHeader *pAODHeader = dynamic_cast(pAOD->GetHeader()); + if(!pAODHeader) AliFatal("Not a standard AOD"); if (NULL == pAODHeader) return; AliAODMCHeader *pAODMCHeader(dynamic_cast(pAOD->FindListObject(AliAODMCHeader::StdBranchName()))); diff --git a/PWGCF/Correlations/DPhi/AliAnalysisTaskMinijet.cxx b/PWGCF/Correlations/DPhi/AliAnalysisTaskMinijet.cxx index 16a969349cb..852ba5baab5 100644 --- a/PWGCF/Correlations/DPhi/AliAnalysisTaskMinijet.cxx +++ b/PWGCF/Correlations/DPhi/AliAnalysisTaskMinijet.cxx @@ -552,7 +552,7 @@ void AliAnalysisTaskMinijet::UserExec(Option_t *) { AliCentrality *centralityObj = 0; if (fAODEvent) - centralityObj = fAODEvent->GetHeader()->GetCentralityP(); + centralityObj = ((AliVAODHeader*)fAODEvent->GetHeader())->GetCentralityP(); else if (fESDEvent) centralityObj = fESDEvent->GetCentrality(); if (centralityObj) diff --git a/PWGCF/Correlations/DPhi/AliAnalysisTaskPhiCorrelations.cxx b/PWGCF/Correlations/DPhi/AliAnalysisTaskPhiCorrelations.cxx index 7761f70b6b4..c74e0e75083 100644 --- a/PWGCF/Correlations/DPhi/AliAnalysisTaskPhiCorrelations.cxx +++ b/PWGCF/Correlations/DPhi/AliAnalysisTaskPhiCorrelations.cxx @@ -518,7 +518,7 @@ void AliAnalysisTaskPhiCorrelations::AnalyseCorrectionMode() { AliCentrality *centralityObj = 0; if (fAOD) - centralityObj = fAOD->GetHeader()->GetCentralityP(); + centralityObj = ((AliVAODHeader*)fAOD->GetHeader())->GetCentralityP(); else if (fESD) centralityObj = fESD->GetCentrality(); @@ -1106,7 +1106,7 @@ void AliAnalysisTaskPhiCorrelations::AnalyseDataMode() else { if (fAOD) - centralityObj = fAOD->GetHeader()->GetCentralityP(); + centralityObj = ((AliVAODHeader*)fAOD->GetHeader())->GetCentralityP(); else if (fESD) centralityObj = fESD->GetCentrality(); @@ -1263,7 +1263,7 @@ void AliAnalysisTaskPhiCorrelations::AnalyseDataMode() referenceMultiplicity = AliESDtrackCuts::GetReferenceMultiplicity(fESD); else if (fAOD) referenceMultiplicity = tracks->GetEntriesFast(); // TODO to be replaced by the estimator once available in the AOD -// referenceMultiplicity = fAOD->GetHeader()->GetRefMultiplicityComb05(); +// referenceMultiplicity = ((AliVAODHeader*)fAOD->GetHeader())->GetRefMultiplicityComb05(); ((TH2F*) fListOfHistos->FindObject("referenceMultiplicity"))->Fill(centrality, referenceMultiplicity); diff --git a/PWGCF/Correlations/DPhi/AliAnalysisTaskTwoPlusOne.cxx b/PWGCF/Correlations/DPhi/AliAnalysisTaskTwoPlusOne.cxx index d4b469ca9d4..3537fc002fe 100644 --- a/PWGCF/Correlations/DPhi/AliAnalysisTaskTwoPlusOne.cxx +++ b/PWGCF/Correlations/DPhi/AliAnalysisTaskTwoPlusOne.cxx @@ -123,7 +123,7 @@ void AliAnalysisTaskTwoPlusOne::UserExec(Option_t *) Double_t centrality = 0; AliCentrality *centralityObj = 0; - centralityObj = fAOD->GetHeader()->GetCentralityP(); + centralityObj = ((AliVAODHeader*)fAOD->GetHeader())->GetCentralityP(); if (centralityObj) diff --git a/PWGCF/Correlations/DPhi/AliDptDptInMC.cxx b/PWGCF/Correlations/DPhi/AliDptDptInMC.cxx index 368121be211..4ca3f12802e 100644 --- a/PWGCF/Correlations/DPhi/AliDptDptInMC.cxx +++ b/PWGCF/Correlations/DPhi/AliDptDptInMC.cxx @@ -1219,8 +1219,8 @@ void AliDptDptInMC::UserExec(Option_t */*option*/) { //Centrality - //AliAODHeader* centralityObject = fAODEvent->GetHeader()->GetCentralityP(); - AliCentrality* centralityObject = fAODEvent->GetHeader()->GetCentralityP(); + //AliAODHeader* centralityObject = ((AliVAODHeader*)fAODEvent->GetHeader())->GetCentralityP(); + AliCentrality* centralityObject = ((AliVAODHeader*)fAODEvent->GetHeader())->GetCentralityP(); if (centralityObject) { v0Centr = centralityObject->GetCentralityPercentile("V0M"); diff --git a/PWGCF/Correlations/DPhi/AliLeadingV0Correlation.cxx b/PWGCF/Correlations/DPhi/AliLeadingV0Correlation.cxx index 1fdc1c75ee4..ecbf1620830 100644 --- a/PWGCF/Correlations/DPhi/AliLeadingV0Correlation.cxx +++ b/PWGCF/Correlations/DPhi/AliLeadingV0Correlation.cxx @@ -441,7 +441,9 @@ void AliLeadingV0Correlation::UserExec(Option_t *) // Before Physics Selection //------------------------------------------------ ltrackMultiplicity = (InputEvent())->GetNumberOfTracks(); - lrefMultiplicity =fAODEvent->GetHeader()->GetRefMultiplicity(); + AliAODHeader * header = dynamic_cast(fAODEvent->GetHeader()); + if(!header) AliFatal("Not a standard AOD"); + lrefMultiplicity = header->GetRefMultiplicity(); fHist_Mult_B4_Trg_Sel->Fill(ltrackMultiplicity,lrefMultiplicity); diff --git a/PWGCF/Correlations/DPhi/DiHadronPID/AliAODEventCutsDiHadronPID.cxx b/PWGCF/Correlations/DPhi/DiHadronPID/AliAODEventCutsDiHadronPID.cxx index 7cd03d912c1..abd9a5b0624 100644 --- a/PWGCF/Correlations/DPhi/DiHadronPID/AliAODEventCutsDiHadronPID.cxx +++ b/PWGCF/Correlations/DPhi/DiHadronPID/AliAODEventCutsDiHadronPID.cxx @@ -329,7 +329,8 @@ Bool_t AliAODEventCutsDiHadronPID::IsSelected(AliAODEvent* event) { } // Get the event header. - AliAODHeader* CurrentHeader = event->GetHeader(); + AliAODHeader* CurrentHeader = dynamic_cast(event->GetHeader()); + if(!CurrentHeader) AliFatal("Not a standard AOD"); // Test minimum reference multiplicity. Int_t CurrentRefMultiplicity = CurrentHeader->GetRefMultiplicity(); diff --git a/PWGCF/Correlations/DPhi/FourierDecomposition/AliDhcTask.cxx b/PWGCF/Correlations/DPhi/FourierDecomposition/AliDhcTask.cxx index 00ea2591f94..91d04170402 100644 --- a/PWGCF/Correlations/DPhi/FourierDecomposition/AliDhcTask.cxx +++ b/PWGCF/Correlations/DPhi/FourierDecomposition/AliDhcTask.cxx @@ -491,7 +491,10 @@ void AliDhcTask::UserExec(Option_t *) Info("Exec()", "Event REJECTED (AOD vertex not OK). z = %.1f", fZVertex); return; } - const AliCentrality *aodCent = fAOD->GetHeader()->GetCentralityP(); + AliAODHeader * header = dynamic_cast(fAOD); + if(!header) AliFatal("Not a standard AOD"); + + const AliCentrality *aodCent = header->GetCentralityP(); if (aodCent) { fCentrality = aodCent->GetCentralityPercentile(fCentMethod); } diff --git a/PWGCF/Correlations/DPhi/MuonHadron/AliAnalysisTaskDiMuonCorrelations.cxx b/PWGCF/Correlations/DPhi/MuonHadron/AliAnalysisTaskDiMuonCorrelations.cxx index be6457e314a..f0832b80dbf 100644 --- a/PWGCF/Correlations/DPhi/MuonHadron/AliAnalysisTaskDiMuonCorrelations.cxx +++ b/PWGCF/Correlations/DPhi/MuonHadron/AliAnalysisTaskDiMuonCorrelations.cxx @@ -511,7 +511,7 @@ Int_t AliAnalysisTaskDiMuonCorrelations::GetCentBin() { Double_t AliAnalysisTaskDiMuonCorrelations::GetITSMultiplicity() { - Double_t multiplicity = fAOD->GetHeader()->GetNumberOfITSClusters(1); + Double_t multiplicity = ((AliVAODHeader*)fAOD->GetHeader())->GetNumberOfITSClusters(1); return multiplicity; diff --git a/PWGCF/Correlations/DPhi/MuonHadron/AliAnalysisTaskMuonHadronCorrelations.cxx b/PWGCF/Correlations/DPhi/MuonHadron/AliAnalysisTaskMuonHadronCorrelations.cxx index 652d1c6e1bd..bee7237c970 100644 --- a/PWGCF/Correlations/DPhi/MuonHadron/AliAnalysisTaskMuonHadronCorrelations.cxx +++ b/PWGCF/Correlations/DPhi/MuonHadron/AliAnalysisTaskMuonHadronCorrelations.cxx @@ -576,7 +576,7 @@ Int_t AliAnalysisTaskMuonHadronCorrelations::GetCentBin() { Double_t AliAnalysisTaskMuonHadronCorrelations::GetITSMultiplicity() { - Double_t multiplicity = fAOD->GetHeader()->GetNumberOfITSClusters(1); + Double_t multiplicity = ((AliVAODHeader*)fAOD->GetHeader())->GetNumberOfITSClusters(1); return multiplicity; diff --git a/PWGCF/Correlations/DPhi/TriggerPID/AliTwoParticlePIDCorr.cxx b/PWGCF/Correlations/DPhi/TriggerPID/AliTwoParticlePIDCorr.cxx index 877f679dbfc..5d8735e3e1e 100644 --- a/PWGCF/Correlations/DPhi/TriggerPID/AliTwoParticlePIDCorr.cxx +++ b/PWGCF/Correlations/DPhi/TriggerPID/AliTwoParticlePIDCorr.cxx @@ -1927,7 +1927,7 @@ if(tracksMCtruth) delete tracksMCtruth; //now deal with reco tracks - Float_t bSign1=aod->GetHeader()->GetMagneticField() ;//used for reconstructed track dca cut + Float_t bSign1=((AliVAODHeader*)aod->GetHeader())->GetMagneticField() ;//used for reconstructed track dca cut //detrmine the ref mult in case of Reco(not required if we get centrality info from AliCentrality) if (fSampleType=="pp_2_76" || fCentralityMethod.EndsWith("_MANUAL") || (fSampleType=="pp_7" && fPPVsMultUtils==kFALSE)) cent_v0=refmultReco; @@ -2039,10 +2039,10 @@ isduplicate2=kTRUE; AliAODTrack *PIDtrack=track;//for PID purpose, mainly important for TPC only tracks if(fFilterBit==128){ -Int_t gid1 = track->GetID(); -//if(gid1>=0) PIDtrack = track; - PIDtrack = aod->GetTrack(trackMap->GetValue(-1-gid1)); -if(!PIDtrack) continue;//for safety; so that each of the TPC only tracks have corresponding global track along with it + Int_t gid1 = track->GetID(); + //if(gid1>=0) PIDtrack = track; + PIDtrack = dynamic_cast(aod->GetTrack(trackMap->GetValue(-1-gid1))); + if(!PIDtrack) continue;//for safety; so that each of the TPC only tracks have corresponding global track along with it } trackscount++; @@ -2393,7 +2393,7 @@ if (!fPID) return;//this should be available with each event even if we don't do bSign = (aod->GetMagneticField() > 0) ? 1 : -1;//for two track efficiency cut in correlation function calculation - Float_t bSign1=aod->GetHeader()->GetMagneticField() ;//for dca cut in ClassifyTrack(), i.e in track loop + Float_t bSign1=((AliVAODHeader*)aod->GetHeader())->GetMagneticField() ;//for dca cut in ClassifyTrack(), i.e in track loop // check event cuts and fill event histograms and return the centrality or reference multiplicity value @@ -2453,10 +2453,10 @@ TExMap *trackMap = new TExMap(); AliAODTrack *PIDtrack=track;//for PID purpose, mainly important for TPC only tracks if(fFilterBit==128){ -Int_t gid1 = track->GetID(); -//if(gid1>=0) PIDtrack = track; - PIDtrack = aod->GetTrack(trackMap->GetValue(-1-gid1)); -if(!PIDtrack) continue;//for safety; so that each of the TPC only tracks have corresponding global track along with it + Int_t gid1 = track->GetID(); + //if(gid1>=0) PIDtrack = track; + PIDtrack = dynamic_cast(aod->GetTrack(trackMap->GetValue(-1-gid1))); + if(!PIDtrack) continue;//for safety; so that each of the TPC only tracks have corresponding global track along with it } //check for eta , phi holes @@ -4782,7 +4782,7 @@ if (v0Centr < 80){ // analysis only for 0-80% centrality classes for(Int_t iT = 0; iT < nAODTracks; iT++) { - AliAODTrack* aodTrack = event->GetTrack(iT); + AliAODTrack* aodTrack = dynamic_cast(event->GetTrack(iT)); if (!aodTrack){ continue; diff --git a/PWGCF/EBYE/BalanceFunctions/AliAnalysisTaskBF.cxx b/PWGCF/EBYE/BalanceFunctions/AliAnalysisTaskBF.cxx index 3390ec7e28d..c6dda4dae65 100755 --- a/PWGCF/EBYE/BalanceFunctions/AliAnalysisTaskBF.cxx +++ b/PWGCF/EBYE/BalanceFunctions/AliAnalysisTaskBF.cxx @@ -670,7 +670,8 @@ void AliAnalysisTaskBF::UserExec(Option_t *) { // for HBT like cuts need magnetic field sign bSign = (gAOD->GetMagneticField() > 0) ? 1 : -1; - AliAODHeader *aodHeader = gAOD->GetHeader(); + AliAODHeader *aodHeader = dynamic_cast(gAOD->GetHeader()); + if(!aodHeader) AliFatal("Not a standard AOD"); // store offline trigger bits fHistTriggerStats->Fill(aodHeader->GetOfflineTrigger()); diff --git a/PWGCF/EBYE/BalanceFunctions/AliAnalysisTaskEventMixingBF.cxx b/PWGCF/EBYE/BalanceFunctions/AliAnalysisTaskEventMixingBF.cxx index ee6c1d29b86..62a2e795ddb 100755 --- a/PWGCF/EBYE/BalanceFunctions/AliAnalysisTaskEventMixingBF.cxx +++ b/PWGCF/EBYE/BalanceFunctions/AliAnalysisTaskEventMixingBF.cxx @@ -434,7 +434,8 @@ void AliAnalysisTaskEventMixingBF::UserExecMix(Option_t *) // for HBT like cuts need magnetic field sign bSign = (aodEventMain->GetMagneticField() > 0) ? 1 : -1; - AliAODHeader *aodHeaderMain = aodEventMain->GetHeader(); + AliAODHeader *aodHeaderMain = dynamic_cast(aodEventMain->GetHeader()); + if(!aodHeaderMain) AliFatal("Not a standard AOD"); // event selection done in AliAnalysisTaskSE::Exec() --> this is not used fHistEventStats->Fill(1); //all events diff --git a/PWGCF/EBYE/Fluctuations/AliEbyEHigherMomentsEffContTask.cxx b/PWGCF/EBYE/Fluctuations/AliEbyEHigherMomentsEffContTask.cxx index abfe26dc1c5..f44d3f93bf1 100644 --- a/PWGCF/EBYE/Fluctuations/AliEbyEHigherMomentsEffContTask.cxx +++ b/PWGCF/EBYE/Fluctuations/AliEbyEHigherMomentsEffContTask.cxx @@ -423,7 +423,8 @@ void AliEbyEHigherMomentsEffContTask::doAODEvent(){ } - AliAODHeader *aodHeader = fAOD->GetHeader(); + AliAODHeader *aodHeader = dynamic_cast(fAOD->GetHeader()); + if(!aodHeader) AliFatal("Not a standard AOD"); Int_t cent = -1; cent = aodHeader->GetCentralityP()->GetCentralityClass10(fCentralityEstimator.Data()); @@ -630,7 +631,8 @@ void AliEbyEHigherMomentsEffContTask::doMCAODEvent(){ return; } - AliAODHeader *aodHeader = fAOD->GetHeader(); + AliAODHeader *aodHeader = dynamic_cast(fAOD->GetHeader()); + if(!aodHeader) AliFatal("Not a standard AOD"); Int_t cent = -1; cent = aodHeader->GetCentralityP()->GetCentralityClass10(fCentralityEstimator.Data()); diff --git a/PWGCF/EBYE/Fluctuations/AliEbyEHigherMomentsTask.cxx b/PWGCF/EBYE/Fluctuations/AliEbyEHigherMomentsTask.cxx index 472e620d834..3dd76a58d35 100644 --- a/PWGCF/EBYE/Fluctuations/AliEbyEHigherMomentsTask.cxx +++ b/PWGCF/EBYE/Fluctuations/AliEbyEHigherMomentsTask.cxx @@ -240,7 +240,8 @@ void AliEbyEHigherMomentsTask::doAODEvent(){ if(!ProperVertex(fAOD)) return; - AliAODHeader *aodHeader = fAOD->GetHeader(); + AliAODHeader *aodHeader = dynamic_cast(fAOD->GetHeader()); + if(!aodHeader) AliFatal("Not a standard AOD"); fCentrality = (Int_t)aodHeader->GetCentralityP()->GetCentralityPercentile(fCentralityEstimator.Data()); @@ -361,7 +362,8 @@ void AliEbyEHigherMomentsTask::doMCAODEvent(){ if(!ProperVertex(fAOD)) return; - AliAODHeader *aodHeader = fAOD->GetHeader(); + AliAODHeader *aodHeader = dynamic_cast(fAOD->GetHeader()); + if(!aodHeader) AliFatal("Not a standard AOD"); fCentrality = (Int_t)aodHeader->GetCentralityP()->GetCentralityPercentile(fCentralityEstimator.Data()); diff --git a/PWGCF/EBYE/Fluctuations/AliEbyEHigherMomentsTaskPID.cxx b/PWGCF/EBYE/Fluctuations/AliEbyEHigherMomentsTaskPID.cxx index fc21036ac30..e36618ea3c2 100644 --- a/PWGCF/EBYE/Fluctuations/AliEbyEHigherMomentsTaskPID.cxx +++ b/PWGCF/EBYE/Fluctuations/AliEbyEHigherMomentsTaskPID.cxx @@ -316,7 +316,8 @@ void AliEbyEHigherMomentsTaskPID::doAODEvent(){ if(!ProperVertex(fAOD)) return; - AliAODHeader *aodHeader = fAOD->GetHeader(); + AliAODHeader *aodHeader = dynamic_cast(fAOD->GetHeader()); + if(!aodHeader) AliFatal("Not a standard AOD"); fCentrality = (Int_t)aodHeader->GetCentralityP()->GetCentralityPercentile(fCentralityEstimator.Data()); @@ -462,7 +463,8 @@ void AliEbyEHigherMomentsTaskPID::doMCAODEvent(){ if(!ProperVertex(fAOD)) return; - AliAODHeader *aodHeader = fAOD->GetHeader(); + AliAODHeader *aodHeader = dynamic_cast(fAOD->GetHeader()); + if(!aodHeader) AliFatal("Not a standard AOD"); fCentrality = (Int_t)aodHeader->GetCentralityP()->GetCentralityPercentile(fCentralityEstimator.Data()); diff --git a/PWGCF/EBYE/Fluctuations/AliEbyEMultFluctuationTask.cxx b/PWGCF/EBYE/Fluctuations/AliEbyEMultFluctuationTask.cxx index 17ed978e1ce..aaae6559c5a 100644 --- a/PWGCF/EBYE/Fluctuations/AliEbyEMultFluctuationTask.cxx +++ b/PWGCF/EBYE/Fluctuations/AliEbyEMultFluctuationTask.cxx @@ -317,7 +317,8 @@ if(!isSelected) return; - AliAODHeader *fHeader = fAOD->GetHeader(); + AliAODHeader *fHeader = dynamic_cast(fAOD->GetHeader()); + if(!fHeader) AliFatal("Not a standard AOD"); Int_t spdmult0 = fHeader->GetNumberOfITSClusters(0); Int_t spdmult1 = fHeader->GetNumberOfITSClusters(1); Int_t run = fHeader->GetRunNumber(); diff --git a/PWGCF/EBYE/Fluctuations/AliHigherMomentsToyModel.cxx b/PWGCF/EBYE/Fluctuations/AliHigherMomentsToyModel.cxx index ea81e8ba8a5..fc92bdd8e7c 100644 --- a/PWGCF/EBYE/Fluctuations/AliHigherMomentsToyModel.cxx +++ b/PWGCF/EBYE/Fluctuations/AliHigherMomentsToyModel.cxx @@ -356,7 +356,8 @@ void AliHigherMomentsToyModel::doAODEvent(){ } - AliAODHeader *aodHeader = fAOD->GetHeader(); + AliAODHeader *aodHeader = dynamic_cast(fAOD->GetHeader()); + if(!aodHeader) AliFatal("Not a standard AOD"); fCentrality = (Int_t)aodHeader->GetCentralityP()->GetCentralityPercentile(fCentralityEstimator.Data()); /* Int_t cent = -1; @@ -573,7 +574,8 @@ void AliHigherMomentsToyModel::doMCAODEvent(){ return; } - AliAODHeader *aodHeader = fAOD->GetHeader(); + AliAODHeader *aodHeader = dynamic_cast(fAOD->GetHeader()); + if(!aodHeader) AliFatal("Not a standard AOD"); fCentrality = (Int_t)aodHeader->GetCentralityP()->GetCentralityPercentile(fCentralityEstimator.Data()); diff --git a/PWGCF/EBYE/NetParticle/AliAnalysisNetParticleHelper.cxx b/PWGCF/EBYE/NetParticle/AliAnalysisNetParticleHelper.cxx index 14dfda9809f..24d6b294b4e 100644 --- a/PWGCF/EBYE/NetParticle/AliAnalysisNetParticleHelper.cxx +++ b/PWGCF/EBYE/NetParticle/AliAnalysisNetParticleHelper.cxx @@ -377,7 +377,7 @@ Int_t AliAnalysisNetParticleHelper::SetupEvent(AliESDInputHandler *esdHandler, A if(esdHandler) centrality = fESD->GetCentrality(); else if(aodHandler) - centrality = fAOD->GetHeader()->GetCentralityP(); + centrality = ((AliVAODHeader*)fAOD->GetHeader())->GetCentralityP(); if (!centrality) { AliError("Centrality not available"); diff --git a/PWGCF/EBYE/NetParticle/AliAnalysisTaskNetParticle.cxx b/PWGCF/EBYE/NetParticle/AliAnalysisTaskNetParticle.cxx index cfd82d5c06a..5f9466e070f 100644 --- a/PWGCF/EBYE/NetParticle/AliAnalysisTaskNetParticle.cxx +++ b/PWGCF/EBYE/NetParticle/AliAnalysisTaskNetParticle.cxx @@ -536,7 +536,7 @@ Int_t AliAnalysisTaskNetParticle::SetupAODEvent() { // -- Check Centrality // ------------------------------------------------------------------ - if (!fAOD->GetHeader()->GetCentralityP()) { + if (!((AliVAODHeader*)fAOD->GetHeader())->GetCentralityP()) { AliError("Could not get centrality"); return -1; } diff --git a/PWGCF/EBYE/PIDFluctuation/task/AliEbyEParticleRatioFluctuationTask.cxx b/PWGCF/EBYE/PIDFluctuation/task/AliEbyEParticleRatioFluctuationTask.cxx index 7d6eb9b11e0..6abe3158ce7 100644 --- a/PWGCF/EBYE/PIDFluctuation/task/AliEbyEParticleRatioFluctuationTask.cxx +++ b/PWGCF/EBYE/PIDFluctuation/task/AliEbyEParticleRatioFluctuationTask.cxx @@ -143,7 +143,8 @@ void AliEbyEParticleRatioFluctuationTask::UserExec( Option_t * ){ Int_t gCent = -1; Float_t gRefMul = -1; - AliAODHeader *aodHeader = event->GetHeader(); + AliAODHeader *aodHeader = dynamic_cast(event->GetHeader()); + if(!aodHeader) AliFatal("Not a standard AOD"); gCent = (Int_t)aodHeader->GetCentralityP()->GetCentralityPercentile(fCentralityEstimator.Data()); gRefMul = aodHeader->GetRefMultiplicity(); if (gCent < 0 || gCent > 100) return; diff --git a/PWGCF/EBYE/PIDFluctuation/task/AliEbyEPidTTask.cxx b/PWGCF/EBYE/PIDFluctuation/task/AliEbyEPidTTask.cxx index 9e7b330b1ab..99c62752644 100644 --- a/PWGCF/EBYE/PIDFluctuation/task/AliEbyEPidTTask.cxx +++ b/PWGCF/EBYE/PIDFluctuation/task/AliEbyEPidTTask.cxx @@ -152,7 +152,8 @@ void AliEbyEPidTTask::UserExec( Option_t * ){ Int_t gCent = -1; // Float_t gRefMul = -1; - AliAODHeader *aodHeader = event->GetHeader(); + AliAODHeader *aodHeader = dynamic_cast(event->GetHeader()); + if(!aodHeader) AliFatal("Not a standard AOD"); gCent = (Int_t)aodHeader->GetCentralityP()->GetCentralityPercentile(fCentralityEstimator.Data()); // gRefMul = aodHeader->GetRefMultiplicity(); if (gCent < 0 || gCent > 100) return; diff --git a/PWGCF/EBYE/PIDFluctuation/task/AliEbyEPidTTaskMC.cxx b/PWGCF/EBYE/PIDFluctuation/task/AliEbyEPidTTaskMC.cxx index 78461a8c77b..3bb5fe109e8 100644 --- a/PWGCF/EBYE/PIDFluctuation/task/AliEbyEPidTTaskMC.cxx +++ b/PWGCF/EBYE/PIDFluctuation/task/AliEbyEPidTTaskMC.cxx @@ -149,7 +149,8 @@ void AliEbyEPidTTaskMC::UserExec( Option_t * ){ Int_t gCent = -1; //Float_t gRefMul = -1; - AliAODHeader *aodHeader = event->GetHeader(); + AliAODHeader *aodHeader = dynamic_cast(event->GetHeader()); + if(!aodHeader) AliFatal("Not a standard AOD"); gCent = (Int_t)aodHeader->GetCentralityP()->GetCentralityPercentile(fCentralityEstimator.Data()); // gRefMul = aodHeader->GetRefMultiplicity(); if (gCent < 0 || gCent > 100) return; diff --git a/PWGCF/FEMTOSCOPY/AliFemto/AliAnalysisTaskFemto.cxx b/PWGCF/FEMTOSCOPY/AliFemto/AliAnalysisTaskFemto.cxx index 02b3317c6a0..c237b2cb237 100644 --- a/PWGCF/FEMTOSCOPY/AliFemto/AliAnalysisTaskFemto.cxx +++ b/PWGCF/FEMTOSCOPY/AliFemto/AliAnalysisTaskFemto.cxx @@ -265,7 +265,8 @@ void AliAnalysisTaskFemto::ConnectInputData(Option_t *) { cout<<"AliAnalysisTaskFemto::AodpidUtil:"<SetAODpidUtil(fAODpidUtil); - fAODheader = fAOD->GetHeader(); + fAODheader = dynamic_cast(fAOD->GetHeader()); + if(!fAODheader) AliFatal("Not a standard AOD"); femtoReaderAOD->SetAODheader(fAODheader); } diff --git a/PWGCF/FEMTOSCOPY/AliFemto/AliFemtoEventReaderAOD.cxx b/PWGCF/FEMTOSCOPY/AliFemto/AliFemtoEventReaderAOD.cxx index 88de7f9b745..59071268ed1 100644 --- a/PWGCF/FEMTOSCOPY/AliFemto/AliFemtoEventReaderAOD.cxx +++ b/PWGCF/FEMTOSCOPY/AliFemto/AliFemtoEventReaderAOD.cxx @@ -30,6 +30,7 @@ #include "AliAODpidUtil.h" #include "AliAnalysisUtils.h" +#include "assert.h" #include "AliGenHijingEventHeader.h" ClassImp(AliFemtoEventReaderAOD) @@ -313,7 +314,10 @@ AliFemtoEvent* AliFemtoEventReaderAOD::CopyAODtoFemtoEvent() } } - tEvent->SetReactionPlaneAngle(fEvent->GetHeader()->GetQTheta(0)/2.0); + AliAODHeader * header = dynamic_cast(fEvent->GetHeader()); + assert(header&&"Not a standard AOD"); + + tEvent->SetReactionPlaneAngle(header->GetQTheta(0)/2.0); // Int_t *motherids=0; // if (mcP) { // const int motherTabSize = ((AliAODMCParticle *) mcP->At(mcP->GetEntries()-1))->GetLabel(); @@ -402,7 +406,8 @@ AliFemtoEvent* AliFemtoEventReaderAOD::CopyAODtoFemtoEvent() // looking for global tracks and saving their numbers to copy from them PID information to TPC-only tracks in the main loop over tracks for (int i=0;iGetTrack(i); + const AliAODTrack *aodtrack=dynamic_cast(fEvent->GetTrack(i)); + assert(aodtrack&&"Not a standard AOD"); if (!aodtrack->TestFilterBit(fFilterBit)) { if(aodtrack->GetID() < 0) continue; labels[aodtrack->GetID()] = i; @@ -539,8 +544,9 @@ AliFemtoEvent* AliFemtoEventReaderAOD::CopyAODtoFemtoEvent() // No additional information exists // Read in the normal AliAODTracks - // const AliAODTrack *aodtrack=fEvent->GetTrack(i); // getting the AODtrack directly - AliAODTrack *aodtrack=fEvent->GetTrack(i); // getting the AODtrack directly + // const AliAODTrack *aodtrack=dynamic_cast(fEvent->GetTrack(i)); + AliAODTrack *aodtrack=dynamic_cast(fEvent->GetTrack(i)); + assert(aodtrack&&"Not a standard AOD"); // getting the AODtrack directly @@ -585,11 +591,12 @@ AliFemtoEvent* AliFemtoEventReaderAOD::CopyAODtoFemtoEvent() AliAODTrack *aodtrackpid; if((fFilterBit == (1 << (7))) || fFilterMask == 128) {//for TPC Only tracks we have to copy PID information from corresponding global tracks - aodtrackpid = fEvent->GetTrack(labels[-1-fEvent->GetTrack(i)->GetID()]); + aodtrackpid = dynamic_cast(fEvent->GetTrack(labels[-1-fEvent->GetTrack(i)->GetID()])); } else { - aodtrackpid = fEvent->GetTrack(i); + aodtrackpid = dynamic_cast(fEvent->GetTrack(i)); } + assert(aodtrackpid&&"Not a standard AOD"); CopyPIDtoFemtoTrack(aodtrackpid, trackCopy); diff --git a/PWGCF/FEMTOSCOPY/AliFemto/AliFemtoEventReaderStandard.cxx b/PWGCF/FEMTOSCOPY/AliFemto/AliFemtoEventReaderStandard.cxx index 5c2aa673040..cf8a4dbe0ed 100644 --- a/PWGCF/FEMTOSCOPY/AliFemto/AliFemtoEventReaderStandard.cxx +++ b/PWGCF/FEMTOSCOPY/AliFemto/AliFemtoEventReaderStandard.cxx @@ -32,6 +32,7 @@ #include "AliAODMCParticle.h" #include "AliVertexerTracks.h" +#include "assert.h" ClassImp(AliFemtoEventReaderStandard) @@ -288,8 +289,10 @@ AliFemtoEvent* AliFemtoEventReaderStandard::ReturnHbtEvent() if (!mcP) { cout << "AOD MC information requested, but no particle array found!" << endl; } + AliAODHeader * header = dynamic_cast(fAODEvent->GetHeader()); + assert(header&&"Not a standard AOD"); - hbtEvent->SetReactionPlaneAngle(fAODEvent->GetHeader()->GetQTheta(0)/2.0); + hbtEvent->SetReactionPlaneAngle(header->GetQTheta(0)/2.0); if (mcP) { motherids = new Int_t[((AliAODMCParticle *) mcP->At(mcP->GetEntries()-1))->GetLabel()]; diff --git a/PWGCF/FEMTOSCOPY/AliFemtoUser/AliAnalysisTaskParticleEfficiency.cxx b/PWGCF/FEMTOSCOPY/AliFemtoUser/AliAnalysisTaskParticleEfficiency.cxx index fb92409f698..dabeed70df1 100644 --- a/PWGCF/FEMTOSCOPY/AliFemtoUser/AliAnalysisTaskParticleEfficiency.cxx +++ b/PWGCF/FEMTOSCOPY/AliFemtoUser/AliAnalysisTaskParticleEfficiency.cxx @@ -366,7 +366,8 @@ void AliAnalysisTaskParticleEfficiency::UserExec(Option_t *) //AliESDEvent *esdEvent = dynamic_cast(InputEvent()); AliAODEvent* aodEvent = dynamic_cast(InputEvent()); if (!aodEvent) return; - AliAODHeader *fAODheader = aodEvent->GetHeader(); + AliAODHeader *fAODheader = dynamic_cast(aodEvent->GetHeader()); + if(!fAODheader) AliFatal("Not a standard AOD"); Double_t mult = fAODheader->GetRefMultiplicity(); // AliCentrality* alicent= aodEvent->GetCentrality(); //in PbPb and pPb // Double_t mult = alicent->GetCentralityPercentile("V0A"); //in pPb diff --git a/PWGCF/FEMTOSCOPYAOD/AliAnalysisTaskPWG2ESDfilter.cxx b/PWGCF/FEMTOSCOPYAOD/AliAnalysisTaskPWG2ESDfilter.cxx index 1d28b093a44..acf3ef97e5a 100644 --- a/PWGCF/FEMTOSCOPYAOD/AliAnalysisTaskPWG2ESDfilter.cxx +++ b/PWGCF/FEMTOSCOPYAOD/AliAnalysisTaskPWG2ESDfilter.cxx @@ -191,7 +191,8 @@ void AliAnalysisTaskPWG2ESDfilter::Exec(Option_t */*option*/) // Update the header - AliAODHeader* header = fAOD->GetHeader(); + AliAODHeader* header = dynamic_cast(fAOD->GetHeader()); + if(!header) AliFatal("Not a standard AOD"); header->SetRunNumber(fESD->GetRunNumber()); if (old) { header->SetBunchCrossNumber(0); diff --git a/PWGDQ/dielectron/AliAnalysisTaskDielectronFilter.cxx b/PWGDQ/dielectron/AliAnalysisTaskDielectronFilter.cxx index 9c435c8ae67..8d9b7ab2b0e 100644 --- a/PWGDQ/dielectron/AliAnalysisTaskDielectronFilter.cxx +++ b/PWGDQ/dielectron/AliAnalysisTaskDielectronFilter.cxx @@ -326,11 +326,15 @@ void AliAnalysisTaskDielectronFilter::UserExec(Option_t *) nanoEv->GetVertex(0)->SetNContributors((static_cast(InputEvent()))->GetPrimaryVertex()->GetNContributors()); nanoEv->GetVertex(1)->SetNContributors((static_cast(InputEvent()))->GetPrimaryVertexSPD()->GetNContributors()); // set event plane - nanoEv->GetHeader()->SetEventplane((static_cast(InputEvent()))->GetHeader()->GetEventplaneP()); - nanoEv->GetHeader()->ResetEventplanePointer(); + AliAODHeader * header = dynamic_cast(nanoEv->GetHeader()); + if(!header) AliFatal("Not a standard AOD"); + + + header->SetEventplane(((AliVAODHeader*)static_cast(InputEvent())->GetHeader())->GetEventplaneP()); + header->ResetEventplanePointer(); // set multiplicity - nanoEv->GetHeader()->SetRefMultiplicity((Int_t)values[AliDielectronVarManager::kNTrk]); - nanoEv->GetHeader()->SetRefMultiplicityPos((Int_t)values[AliDielectronVarManager::kNacc]); + header->SetRefMultiplicity((Int_t)values[AliDielectronVarManager::kNTrk]); + header->SetRefMultiplicityPos((Int_t)values[AliDielectronVarManager::kNacc]); //nanoEv->GetHeader()->SetRefMultiplicityNeg(values[AliDielectronVarManager::kMatchEffITSTPC]); for(int kj=0; kj<(fDielectron->GetTrackArray(0))->GetEntries(); kj++){ @@ -384,10 +388,12 @@ void AliAnalysisTaskDielectronFilter::UserExec(Option_t *) if(fCreateNanoAOD && isAOD && (!hasCand) && fStoreHeader) { - // set event plane - extDielectron->GetAOD()->GetHeader()->SetEventplane((static_cast(InputEvent()))->GetHeader()->GetEventplaneP()); - extDielectron->GetAOD()->GetHeader()->ResetEventplanePointer(); - extDielectron->GetTree()->Fill(); // fill header for all events without tracks + // set event plane + AliAODHeader * header = dynamic_cast(extDielectron->GetAOD()->GetHeader()); + if(!header) AliFatal("Not a standard AOD"); + header->SetEventplane(((AliAODHeader*)(static_cast(InputEvent()))->GetHeader())->GetEventplaneP()); + header->ResetEventplanePointer(); + extDielectron->GetTree()->Fill(); // fill header for all events without tracks } PostData(1, const_cast(fDielectron->GetHistogramList())); diff --git a/PWGDQ/dielectron/AliAnalysisTaskReducedTree.cxx b/PWGDQ/dielectron/AliAnalysisTaskReducedTree.cxx index a8d4fa9adf5..2ccc92652fd 100644 --- a/PWGDQ/dielectron/AliAnalysisTaskReducedTree.cxx +++ b/PWGDQ/dielectron/AliAnalysisTaskReducedTree.cxx @@ -420,12 +420,16 @@ void AliAnalysisTaskReducedTree::FillEventInfo() } } if(isAOD) { - fReducedEvent->fIRIntClosestIntMap[0] = aodEvent->GetHeader()->GetIRInt1ClosestInteractionMap(); - fReducedEvent->fIRIntClosestIntMap[1] = aodEvent->GetHeader()->GetIRInt2ClosestInteractionMap(); - fReducedEvent->fEventNumberInFile = aodEvent->GetHeader()->GetEventNumberESDFile(); - fReducedEvent->fL0TriggerInputs = aodEvent->GetHeader()->GetL0TriggerInputs(); - fReducedEvent->fL1TriggerInputs = aodEvent->GetHeader()->GetL1TriggerInputs(); - fReducedEvent->fL2TriggerInputs = aodEvent->GetHeader()->GetL2TriggerInputs(); + AliAODHeader * header = dynamic_cast(aodEvent->GetHeader()); + if(!header) AliFatal("Not a standard AOD"); + + + fReducedEvent->fIRIntClosestIntMap[0] = header->GetIRInt1ClosestInteractionMap(); + fReducedEvent->fIRIntClosestIntMap[1] = header->GetIRInt2ClosestInteractionMap(); + fReducedEvent->fEventNumberInFile = header->GetEventNumberESDFile(); + fReducedEvent->fL0TriggerInputs = header->GetL0TriggerInputs(); + fReducedEvent->fL1TriggerInputs = header->GetL1TriggerInputs(); + fReducedEvent->fL2TriggerInputs = header->GetL2TriggerInputs(); fReducedEvent->fTimeStamp = 0; fReducedEvent->fNpileupSPD = aodEvent->GetNumberOfPileupVerticesSPD(); fReducedEvent->fNpileupTracks = aodEvent->GetNumberOfPileupVerticesTracks(); diff --git a/PWGDQ/dielectron/AliDielectronVarManager.h b/PWGDQ/dielectron/AliDielectronVarManager.h index a75e1c3d912..4f4e150186a 100644 --- a/PWGDQ/dielectron/AliDielectronVarManager.h +++ b/PWGDQ/dielectron/AliDielectronVarManager.h @@ -73,6 +73,7 @@ #include "AliVZEROEPSelectionTask.h" #include "AliAODMCHeader.h" +#include "assert.h" class AliVEvent; @@ -2013,7 +2014,8 @@ inline void AliDielectronVarManager::FillVarAODEvent(const AliAODEvent *event, D FillVarVEvent(event, values); // Fill AliAODEvent interface specific information - AliAODHeader *header = event->GetHeader(); + AliAODHeader *header = dynamic_cast(event->GetHeader()); + assert(header&&"Not a standard AOD"); Double_t centralityF=-1; Double_t centralitySPD=-1; AliCentrality *aodCentrality = header->GetCentralityP(); @@ -2587,7 +2589,8 @@ inline void AliDielectronVarManager::GetVzeroRP(const AliVEvent* event, Double_t } if(event->IsA() == AliAODEvent::Class()) { const AliAODEvent* aodEv = static_cast(event); - AliAODHeader *header = aodEv->GetHeader(); + AliAODHeader *header = dynamic_cast(aodEv->GetHeader()); + assert(header&&"Not a standard AOD"); AliCentrality *aodCentrality = header->GetCentralityP(); if(aodCentrality) centralitySPD = aodCentrality->GetCentralityPercentile("CL1"); } @@ -2716,7 +2719,7 @@ inline void AliDielectronVarManager::GetZDCRP(const AliVEvent* event, Double_t q if(fgZDCRecentering[0][0]){ const AliAODEvent* aodEv = static_cast(event); - AliAODHeader *header = aodEv->GetHeader(); + AliAODHeader *header = dynamic_cast(aodEv->GetHeader()); if(!header) return; TPCRefMulti = header -> GetTPConlyRefMultiplicity(); diff --git a/PWGGA/CaloTasks/AliAnalysisTaskCaloFilter.cxx b/PWGGA/CaloTasks/AliAnalysisTaskCaloFilter.cxx index e39501c84ad..cadc787ef53 100644 --- a/PWGGA/CaloTasks/AliAnalysisTaskCaloFilter.cxx +++ b/PWGGA/CaloTasks/AliAnalysisTaskCaloFilter.cxx @@ -623,12 +623,13 @@ void AliAnalysisTaskCaloFilter::FillAODHeader() { // AOD header copy - AliAODHeader* header = AODEvent()->GetHeader(); + AliAODHeader* header = dynamic_cast(AODEvent()->GetHeader()); + if(!header) AliFatal("Not a standard AOD"); // Copy from AODs if(fAODEvent) { - *header = *(fAODEvent->GetHeader()); + *header = *((AliAODHeader*)fAODEvent->GetHeader()); return; } diff --git a/PWGGA/EMCALTasks/AliAnalysisTaskEMCALClusterize.cxx b/PWGGA/EMCALTasks/AliAnalysisTaskEMCALClusterize.cxx index 11f67ef4a93..843d62333e3 100644 --- a/PWGGA/EMCALTasks/AliAnalysisTaskEMCALClusterize.cxx +++ b/PWGGA/EMCALTasks/AliAnalysisTaskEMCALClusterize.cxx @@ -849,7 +849,8 @@ void AliAnalysisTaskEMCALClusterize::FillAODHeader() Double_t covVtx[6]; for (Int_t i = 0; i < 6; i++) covVtx[i] = 0.; - AliAODHeader* header = AODEvent()->GetHeader(); + AliAODHeader* header = dynamic_cast(AODEvent()->GetHeader()); + if(!header) AliFatal("Not a standard AOD"); header->SetRunNumber(fEvent->GetRunNumber()); if(esdevent) @@ -861,7 +862,11 @@ void AliAnalysisTaskEMCALClusterize::FillAODHeader() if (file) header->SetESDFileName(file->GetName()); } } - else if (aodevent) header->SetESDFileName(aodevent->GetHeader()->GetESDFileName()); + else if (aodevent) { + AliAODHeader * aodheader = dynamic_cast(aodevent->GetHeader()); + if(!aodheader) AliFatal("Not a standard AOD"); + header->SetESDFileName(aodheader->GetESDFileName()); + } header->SetBunchCrossNumber(fEvent->GetBunchCrossNumber()); header->SetOrbitNumber(fEvent->GetOrbitNumber()); diff --git a/PWGGA/EMCALTasks/AliAnalysisTaskEMCALPi0PbPb.cxx b/PWGGA/EMCALTasks/AliAnalysisTaskEMCALPi0PbPb.cxx index 7d5b93abb74..1dd7c2a3466 100644 --- a/PWGGA/EMCALTasks/AliAnalysisTaskEMCALPi0PbPb.cxx +++ b/PWGGA/EMCALTasks/AliAnalysisTaskEMCALPi0PbPb.cxx @@ -603,7 +603,7 @@ void AliAnalysisTaskEMCALPi0PbPb::UserExec(Option_t *) return; } am->LoadBranch("header"); - offtrigger = fAodEv->GetHeader()->GetOfflineTrigger(); + offtrigger = ((AliVAODHeader*)fAodEv->GetHeader())->GetOfflineTrigger(); } if (!fMcMode && (offtrigger & AliVEvent::kFastOnly)) { AliWarning(Form("EMCAL not in fast only partition")); @@ -626,8 +626,11 @@ void AliAnalysisTaskEMCALPi0PbPb::UserExec(Option_t *) const TGeoHMatrix *geom = 0; if (fEsdEv) geom = fEsdEv->GetESDRun()->GetEMCALMatrix(i); - else - geom = fAodEv->GetHeader()->GetEMCALMatrix(i); + else { + AliAODHeader * aodheader = dynamic_cast(fAodEv->GetHeader()); + if(!aodheader) AliFatal("Not a standard AOD"); + geom = aodheader->GetEMCALMatrix(i); + } if (!geom) continue; geom->Print(); @@ -1614,17 +1617,20 @@ void AliAnalysisTaskEMCALPi0PbPb::FillNtuple() AliAnalysisManager *am = AliAnalysisManager::GetAnalysisManager(); if (fAodEv) { + AliAODHeader * aodheader = dynamic_cast(fAodEv->GetHeader()); + if(!aodheader) AliFatal("Not a standard AOD"); + fHeader->fRun = fAodEv->GetRunNumber(); - fHeader->fOrbit = fAodEv->GetHeader()->GetOrbitNumber(); - fHeader->fPeriod = fAodEv->GetHeader()->GetPeriodNumber(); - fHeader->fBx = fAodEv->GetHeader()->GetBunchCrossNumber(); - fHeader->fL0 = fAodEv->GetHeader()->GetL0TriggerInputs(); - fHeader->fL1 = fAodEv->GetHeader()->GetL1TriggerInputs(); - fHeader->fL2 = fAodEv->GetHeader()->GetL2TriggerInputs(); - fHeader->fTrClassMask = fAodEv->GetHeader()->GetTriggerMask(); - fHeader->fTrCluster = fAodEv->GetHeader()->GetTriggerCluster(); - fHeader->fOffTriggers = fAodEv->GetHeader()->GetOfflineTrigger(); - fHeader->fFiredTriggers = fAodEv->GetHeader()->GetFiredTriggerClasses(); + fHeader->fOrbit = aodheader->GetOrbitNumber(); + fHeader->fPeriod = aodheader->GetPeriodNumber(); + fHeader->fBx = aodheader->GetBunchCrossNumber(); + fHeader->fL0 = aodheader->GetL0TriggerInputs(); + fHeader->fL1 = aodheader->GetL1TriggerInputs(); + fHeader->fL2 = aodheader->GetL2TriggerInputs(); + fHeader->fTrClassMask = aodheader->GetTriggerMask(); + fHeader->fTrCluster = aodheader->GetTriggerCluster(); + fHeader->fOffTriggers = aodheader->GetOfflineTrigger(); + fHeader->fFiredTriggers = aodheader->GetFiredTriggerClasses(); } else { fHeader->fRun = fEsdEv->GetRunNumber(); fHeader->fOrbit = fEsdEv->GetHeader()->GetOrbitNumber(); diff --git a/PWGGA/GammaConv/AliConvEventCuts.cxx b/PWGGA/GammaConv/AliConvEventCuts.cxx index 0998f168f90..40301699ae8 100644 --- a/PWGGA/GammaConv/AliConvEventCuts.cxx +++ b/PWGGA/GammaConv/AliConvEventCuts.cxx @@ -1310,7 +1310,7 @@ Double_t AliConvEventCuts::GetCentrality(AliVEvent *event) AliAODEvent *aodEvent=dynamic_cast(event); if(aodEvent){ - if(aodEvent->GetHeader()){return aodEvent->GetHeader()->GetCentrality();} + if(aodEvent->GetHeader()){return ((AliVAODHeader*)aodEvent->GetHeader())->GetCentrality();} } return -1; diff --git a/PWGGA/GammaConv/AliConversionCuts.cxx b/PWGGA/GammaConv/AliConversionCuts.cxx index 44f229cd1c9..18db1a93440 100644 --- a/PWGGA/GammaConv/AliConversionCuts.cxx +++ b/PWGGA/GammaConv/AliConversionCuts.cxx @@ -3542,7 +3542,7 @@ Double_t AliConversionCuts::GetCentrality(AliVEvent *event) AliAODEvent *aodEvent=dynamic_cast(event); if(aodEvent){ - if(aodEvent->GetHeader()){return aodEvent->GetHeader()->GetCentrality();} + if(aodEvent->GetHeader()){return ((AliVAODHeader*)aodEvent->GetHeader())->GetCentrality();} } return -1; diff --git a/PWGGA/PHOSTasks/PHOS_embedding/AliAnalysisTaskPi0DiffEfficiency.cxx b/PWGGA/PHOSTasks/PHOS_embedding/AliAnalysisTaskPi0DiffEfficiency.cxx index 08b3b93eda0..2d427903b34 100644 --- a/PWGGA/PHOSTasks/PHOS_embedding/AliAnalysisTaskPi0DiffEfficiency.cxx +++ b/PWGGA/PHOSTasks/PHOS_embedding/AliAnalysisTaskPi0DiffEfficiency.cxx @@ -714,7 +714,8 @@ void AliAnalysisTaskPi0DiffEfficiency::UserExec(Option_t *) } FillHistogram("hSelEvents",1.5) ; - AliAODHeader *header = event->GetHeader() ; + AliAODHeader *header = dynamic_cast(event->GetHeader()) ; + if(!header) AliFatal("Not a standard AOD"); // Checks if we have a primary vertex // Get primary vertices form ESD diff --git a/PWGGA/PHOSTasks/PHOS_embedding/AliAnalysisTaskPi0Efficiency.cxx b/PWGGA/PHOSTasks/PHOS_embedding/AliAnalysisTaskPi0Efficiency.cxx index 23575f7fa3c..94bc6360a41 100644 --- a/PWGGA/PHOSTasks/PHOS_embedding/AliAnalysisTaskPi0Efficiency.cxx +++ b/PWGGA/PHOSTasks/PHOS_embedding/AliAnalysisTaskPi0Efficiency.cxx @@ -418,7 +418,8 @@ void AliAnalysisTaskPi0Efficiency::UserExec(Option_t *) } FillHistogram("hSelEvents",1.5) ; - AliAODHeader *header = event->GetHeader() ; + AliAODHeader *header = dynamic_cast(event->GetHeader()) ; + if(!header) AliFatal("Not a standard AOD"); // Checks if we have a primary vertex // Get primary vertices form ESD diff --git a/PWGGA/PHOSTasks/PHOS_embedding/AliPHOSEmbedding.cxx b/PWGGA/PHOSTasks/PHOS_embedding/AliPHOSEmbedding.cxx index e7acd7af0d5..da7710954f8 100644 --- a/PWGGA/PHOSTasks/PHOS_embedding/AliPHOSEmbedding.cxx +++ b/PWGGA/PHOSTasks/PHOS_embedding/AliPHOSEmbedding.cxx @@ -589,7 +589,8 @@ void AliPHOSEmbedding::ConvertMCParticles(const AliAODEvent* aod) //______________________________________________________________________________ void AliPHOSEmbedding::ConvertHeader(AliESDEvent & esd){ - AliAODHeader* header = AODEvent()->GetHeader(); + AliAODHeader* header = dynamic_cast(AODEvent()->GetHeader()); + if(!header) AliFatal("Not a standard AOD"); header->SetRunNumber(esd.GetRunNumber()); header->SetOfflineTrigger(fInputHandler->IsEventSelected()); // propagate the decision of the physics selection diff --git a/PWGHF/correlationHF/AliAnalysisTaskDStarCorrelations.cxx b/PWGHF/correlationHF/AliAnalysisTaskDStarCorrelations.cxx index 66eb93073ff..ca4cd4b2510 100644 --- a/PWGHF/correlationHF/AliAnalysisTaskDStarCorrelations.cxx +++ b/PWGHF/correlationHF/AliAnalysisTaskDStarCorrelations.cxx @@ -1454,7 +1454,7 @@ void AliAnalysisTaskDStarCorrelations::EventMixingChecks(AliAODEvent* AOD){ } if(fSystem == AA){ // PbPb - centralityObj = AOD->GetHeader()->GetCentralityP(); + centralityObj = ((AliVAODHeader*)AOD->GetHeader())->GetCentralityP(); MultipOrCent = centralityObj->GetCentralityPercentileUnchecked("V0M"); AliInfo(Form("Centrality is %f", MultipOrCent)); } diff --git a/PWGHF/correlationHF/AliAnalysisTaskDxHFECorrelation.cxx b/PWGHF/correlationHF/AliAnalysisTaskDxHFECorrelation.cxx index 519210ce303..2ebac97dd3d 100644 --- a/PWGHF/correlationHF/AliAnalysisTaskDxHFECorrelation.cxx +++ b/PWGHF/correlationHF/AliAnalysisTaskDxHFECorrelation.cxx @@ -283,7 +283,7 @@ void AliAnalysisTaskDxHFECorrelation::UserExec(Option_t* /*option*/) Double_t MultipOrCent = -1; AliAODEvent* aodEvent=dynamic_cast(pEvent); if (aodEvent) { - centralityObj = aodEvent->GetHeader()->GetCentralityP(); + centralityObj = ((AliVAODHeader*)aodEvent->GetHeader())->GetCentralityP(); if (centralityObj) { MultipOrCent = centralityObj->GetCentralityPercentileUnchecked("V0M"); } diff --git a/PWGHF/correlationHF/AliDxHFECorrelation.cxx b/PWGHF/correlationHF/AliDxHFECorrelation.cxx index ab09c6778a7..fbdd4219d75 100644 --- a/PWGHF/correlationHF/AliDxHFECorrelation.cxx +++ b/PWGHF/correlationHF/AliDxHFECorrelation.cxx @@ -597,7 +597,7 @@ void AliDxHFECorrelation::EventMixingChecks(const AliVEvent* pEvent){ MultipOrCent = multiplicity; // convert from Int_t to Double_t } if(fSystem){ // PbPb - centralityObj = AOD->GetHeader()->GetCentralityP(); + centralityObj = ((AliVAODHeader*)AOD->GetHeader())->GetCentralityP(); MultipOrCent = centralityObj->GetCentralityPercentileUnchecked("V0M"); AliInfo(Form("Centrality is %f", MultipOrCent)); } diff --git a/PWGHF/correlationHF/AliHFCorrelator.cxx b/PWGHF/correlationHF/AliHFCorrelator.cxx index 54d3b7b7ee3..5881603ad36 100644 --- a/PWGHF/correlationHF/AliHFCorrelator.cxx +++ b/PWGHF/correlationHF/AliHFCorrelator.cxx @@ -261,7 +261,7 @@ Bool_t AliHFCorrelator::Initialize(){ if(fUseCentrality){ // PbPb if(!fDMesonCutObject){ - centralityObj = fAODEvent->GetHeader()->GetCentralityP(); + centralityObj = ((AliVAODHeader*)fAODEvent->GetHeader())->GetCentralityP(); fMultCentr = centralityObj->GetCentralityPercentileUnchecked("V0M"); } else fMultCentr = fDMesonCutObject->GetCentrality(fAODEvent); diff --git a/PWGHF/hfe/AliAnalysisTaskEMCalHFEpA.cxx b/PWGHF/hfe/AliAnalysisTaskEMCalHFEpA.cxx index 243f4df8c4d..026042b07a3 100644 --- a/PWGHF/hfe/AliAnalysisTaskEMCalHFEpA.cxx +++ b/PWGHF/hfe/AliAnalysisTaskEMCalHFEpA.cxx @@ -1835,7 +1835,7 @@ void AliAnalysisTaskEMCalHFEpA::UserExec(Option_t *) if(fIsAOD) { - fCentrality = fAOD->GetHeader()->GetCentralityP(); + fCentrality = ((AliVAODHeader*)fAOD->GetHeader())->GetCentralityP(); } else { @@ -2159,7 +2159,8 @@ void AliAnalysisTaskEMCalHFEpA::UserExec(Option_t *) if(fIsAOD){ - //AliAODHeader * aodh = fAOD->GetHeader(); + //AliAODHeader * aodh = dynamic_cast(fAOD->GetHeader()); + // if(!aodh) AliFatal("Not a standard AOD"); //Int_t bc= aodh->GetBunchCrossNumber(); diff --git a/PWGHF/hfe/AliAnalysisTaskHFEemcQA.cxx b/PWGHF/hfe/AliAnalysisTaskHFEemcQA.cxx index 59a33a4b57d..ee859a72b79 100644 --- a/PWGHF/hfe/AliAnalysisTaskHFEemcQA.cxx +++ b/PWGHF/hfe/AliAnalysisTaskHFEemcQA.cxx @@ -407,7 +407,9 @@ void AliAnalysisTaskHFEemcQA::UserExec(Option_t *) Int_t trigger = -1; if (fAOD){ - Double_t multiplicity=fAOD->GetHeader()->GetRefMultiplicity(); + AliAODHeader * header = dynamic_cast(fAOD->GetHeader()); + if(!header) AliFatal("Not a standard AOD"); + Double_t multiplicity=header->GetRefMultiplicity(); fTrigMulti->Fill(-0.5, multiplicity); if(evSelMask & AliVEvent::kAny) fTrigMulti->Fill(0.5, multiplicity); if(evSelMask & AliVEvent::kMB) fTrigMulti->Fill(1.5, multiplicity); diff --git a/PWGHF/vertexingHF/AliAnalysisTaskSEHFQA.cxx b/PWGHF/vertexingHF/AliAnalysisTaskSEHFQA.cxx index 6fe3ffa0ce2..116dbf68dcf 100644 --- a/PWGHF/vertexingHF/AliAnalysisTaskSEHFQA.cxx +++ b/PWGHF/vertexingHF/AliAnalysisTaskSEHFQA.cxx @@ -1317,7 +1317,10 @@ void AliAnalysisTaskSEHFQA::UserExec(Option_t */*option*/) UInt_t evSelMask=((AliInputEventHandler*)(AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler()))->IsEventSelected(); Double_t centrality=fCuts->GetCentrality(aod); - Double_t multiplicity=aod->GetHeader()->GetRefMultiplicity(); + AliAODHeader * header = dynamic_cast(aod->GetHeader()); + if(!header) AliFatal("Not a standard AOD"); + + Double_t multiplicity=header->GetRefMultiplicity(); Int_t runNumber = aod->GetRunNumber(); TString trigClass=aod->GetFiredTriggerClasses(); Int_t nAODtracks=aod->GetNumberOfTracks(); @@ -1742,14 +1745,18 @@ void AliAnalysisTaskSEHFQA::UserExec(Option_t */*option*/) if(mincent==100)mincent--; ((AliCounterCollection*)fOutputCounters->FindObject("secondEstimator"))->Count(Form("centralityclass:%d_%d/Run:%d",mincent,mincent+10,runNumber)); + AliAODHeader * header = dynamic_cast(aod->GetHeader()); + if(!header) AliFatal("Not a standard AOD"); + + if(stdCentGetMinCentrality() || stdCent>fCuts->GetMaxCentrality()){ ((TH1F*)fOutputCheckCentrality->FindObject("hNtrackletsOut"))->Fill(aod->GetTracklets()->GetNumberOfTracklets()); - ((TH1F*)fOutputCheckCentrality->FindObject("hMultOut"))->Fill(aod->GetHeader()->GetRefMultiplicity()); + ((TH1F*)fOutputCheckCentrality->FindObject("hMultOut"))->Fill(header->GetRefMultiplicity()); }else{ ((TH1F*)fOutputCheckCentrality->FindObject("hNtrackletsIn"))->Fill(aod->GetTracklets()->GetNumberOfTracklets()); - ((TH1F*)fOutputCheckCentrality->FindObject("hMultIn"))->Fill(aod->GetHeader()->GetRefMultiplicity()); + ((TH1F*)fOutputCheckCentrality->FindObject("hMultIn"))->Fill(header->GetRefMultiplicity()); } - ((TH2F*)fOutputCheckCentrality->FindObject("hMultvsPercentile"))->Fill(aod->GetHeader()->GetRefMultiplicity(),stdCentf); + ((TH2F*)fOutputCheckCentrality->FindObject("hMultvsPercentile"))->Fill(header->GetRefMultiplicity(),stdCentf); ((TH2F*)fOutputCheckCentrality->FindObject("hntrklvsPercentile"))->Fill(aod->GetTracklets()->GetNumberOfTracklets(),stdCentf); ((TH2F*)fOutputCheckCentrality->FindObject("hntrklvsPercentile01"))->Fill(AliVertexingHFUtils::GetNumberOfTrackletsInEtaRange(aod,-1.,1.),stdCentf); ((TH2F*)fOutputCheckCentrality->FindObject("hnTPCTracksvsPercentile"))->Fill(nSelTracksTPCOnly,stdCentf); @@ -1767,10 +1774,10 @@ void AliAnalysisTaskSEHFQA::UserExec(Option_t */*option*/) if(fOnOff[0]){ ((TH1F*)fOutputTrack->FindObject("hNtracklets"))->Fill(aod->GetTracklets()->GetNumberOfTracklets()); ((TH1F*)fOutputTrack->FindObject("hNtracklets01"))->Fill(AliVertexingHFUtils::GetNumberOfTrackletsInEtaRange(aod,-1.,1.)); - ((TH1F*)fOutputTrack->FindObject("hMult"))->Fill(aod->GetHeader()->GetRefMultiplicity()); + ((TH1F*)fOutputTrack->FindObject("hMult"))->Fill(header->GetRefMultiplicity()); ((TH1F*)fOutputTrack->FindObject("hMultFBit4"))->Fill(ntracksFBit4); - ((TH1F*)fOutputTrack->FindObject("hMultComb05"))->Fill(aod->GetHeader()->GetRefMultiplicityComb05()); - ((TH1F*)fOutputTrack->FindObject("hMultComb08"))->Fill(aod->GetHeader()->GetRefMultiplicityComb08()); + ((TH1F*)fOutputTrack->FindObject("hMultComb05"))->Fill(header->GetRefMultiplicityComb05()); + ((TH1F*)fOutputTrack->FindObject("hMultComb08"))->Fill(header->GetRefMultiplicityComb08()); } } } diff --git a/PWGHF/vertexingHF/AliRDHFCuts.cxx b/PWGHF/vertexingHF/AliRDHFCuts.cxx index 8810918d0ed..800fe492faa 100644 --- a/PWGHF/vertexingHF/AliRDHFCuts.cxx +++ b/PWGHF/vertexingHF/AliRDHFCuts.cxx @@ -1178,7 +1178,8 @@ Float_t AliRDHFCuts::GetCentrality(AliAODEvent* aodEvent,AliRDHFCuts::ECentralit TClonesArray *mcArray = (TClonesArray*)((AliAODEvent*)aodEvent)->GetList()->FindObject(AliAODMCParticle::StdBranchName()); if(mcArray) {fUseAOD049=kFALSE;} - AliAODHeader *header=aodEvent->GetHeader(); + AliAODHeader *header=dynamic_cast(aodEvent->GetHeader()); + if(!header) AliFatal("Not a standard AOD"); AliCentrality *centrality=header->GetCentralityP(); Float_t cent=-999.; Bool_t isSelRun=kFALSE; diff --git a/PWGJE/AliAnalysisTaskAj.cxx b/PWGJE/AliAnalysisTaskAj.cxx index 94b28883d55..5f613e18325 100644 --- a/PWGJE/AliAnalysisTaskAj.cxx +++ b/PWGJE/AliAnalysisTaskAj.cxx @@ -479,7 +479,7 @@ void AliAnalysisTaskAj::UserExec(Option_t *) Double_t centValue = 0.; if(fESD) {cent = fESD->GetCentrality(); if(cent) centValue = cent->GetCentralityPercentile("V0M");} - else centValue=aod->GetHeader()->GetCentrality(); + else centValue=((AliVAODHeader*)aod->GetHeader())->GetCentrality(); if(fDebug) printf("centrality: %f\n", centValue); // if (centValue < fCentMin || centValue > fCentMax){ diff --git a/PWGJE/AliAnalysisTaskFragmentationFunction.cxx b/PWGJE/AliAnalysisTaskFragmentationFunction.cxx index 3b75fe8966a..ccd5ac70845 100644 --- a/PWGJE/AliAnalysisTaskFragmentationFunction.cxx +++ b/PWGJE/AliAnalysisTaskFragmentationFunction.cxx @@ -2140,7 +2140,7 @@ void AliAnalysisTaskFragmentationFunction::UserExec(Option_t *) Int_t cl = 0; if(handler->InheritsFrom("AliAODInputHandler")){ // since it is not supported by the helper task define own classes - centPercent = fAOD->GetHeader()->GetCentrality(); + centPercent = ((AliVAODHeader*)fAOD->GetHeader())->GetCentrality(); cl = 1; if(centPercent>10) cl = 2; if(centPercent>30) cl = 3; diff --git a/PWGJE/AliAnalysisTaskJetCore.cxx b/PWGJE/AliAnalysisTaskJetCore.cxx index 66e1c556bc2..3562be28467 100644 --- a/PWGJE/AliAnalysisTaskJetCore.cxx +++ b/PWGJE/AliAnalysisTaskJetCore.cxx @@ -665,7 +665,7 @@ void AliAnalysisTaskJetCore::UserExec(Option_t *) if(fIsPbPb){ if(fESD) {cent = fESD->GetCentrality(); if(cent) centValue = cent->GetCentralityPercentile("V0M");} - else centValue=aod->GetHeader()->GetCentrality(); + else centValue=((AliVAODHeader*)aod->GetHeader())->GetCentrality(); if(fDebug) printf("centrality: %f\n", centValue); if (centValue < fCentMin || centValue > fCentMax){ @@ -765,8 +765,8 @@ void AliAnalysisTaskJetCore::UserExec(Option_t *) Int_t iCount=0; Int_t trigJet=-1; Int_t trigBBTrack=-1; - // Int_t trigInTrack=-1; - fRPAngle = aod->GetHeader()->GetEventplane(); + // Int_t trigInTrack=-1; + fRPAngle = ((AliVAODHeader*)aod->GetHeader())->GetEventplane(); if(fHardest==0 || fHardest==1){ AliVParticle *partback = (AliVParticle*)ParticleList.At(nT); @@ -1166,7 +1166,8 @@ Int_t AliAnalysisTaskJetCore::GetListOfTracks(TList *list){ for(int it = 0;it < aod->GetNumberOfTracks();++it){ - AliAODTrack *tr = aod->GetTrack(it); + AliAODTrack *tr = dynamic_cast(aod->GetTrack(it)); + if(!tr) AliFatal("Not a standard AOD"); Bool_t bGood = false; if(fFilterType == 0)bGood = true; else if(fFilterType == 1)bGood = tr->IsHybridTPCConstrainedGlobal(); @@ -1230,7 +1231,8 @@ Int_t AliAnalysisTaskJetCore::SelectTrigger(TList *list,Double_t minT,Double_t for(Int_t cr=0;cr<100;cr++){triggers[cr]=-1;} Int_t im=0; for(int it = 0;it < aod->GetNumberOfTracks();++it){ - AliAODTrack *tr = aod->GetTrack(it); + AliAODTrack *tr = dynamic_cast(aod->GetTrack(it)); + if(!tr) AliFatal("Not a standard AOD"); Bool_t bGood = false; if(fFilterType == 0)bGood = true; else if(fFilterType == 1)bGood = tr->IsHybridTPCConstrainedGlobal(); @@ -1292,7 +1294,8 @@ Int_t AliAnalysisTaskJetCore::SelectTrigger(TList *list,Double_t minT,Double_t // Double_t dif=0; Int_t iCount=0; for(int it = 0;it < aod->GetNumberOfTracks();++it){ - AliAODTrack *tr = aod->GetTrack(it); + AliAODTrack *tr = dynamic_cast(aod->GetTrack(it)); + if(!tr) AliFatal("Not a standard AOD"); if((fFilterMask>0)&&!(tr->TestFilterBit(fFilterMask)))continue; if(TMath::Abs(tr->Eta())>0.9)continue; if(tr->Pt()<0.15)continue; @@ -1324,7 +1327,8 @@ Int_t AliAnalysisTaskJetCore::SelectTrigger(TList *list,Double_t minT,Double_t else aod = fAODOut; for(int it = 0;it < aod->GetNumberOfTracks();++it){ - AliAODTrack *tr = aod->GetTrack(it); + AliAODTrack *tr = dynamic_cast(aod->GetTrack(it)); + if(!tr) AliFatal("Not a standard AOD"); if((fFilterMask>0)&&!(tr->TestFilterBit(fFilterMask)))continue; if(TMath::Abs(tr->Eta())>0.9)continue; if(tr->Pt()<0.15)continue; diff --git a/PWGJE/AliAnalysisTaskJetCorePP.cxx b/PWGJE/AliAnalysisTaskJetCorePP.cxx index 173e610550d..678856c00eb 100644 --- a/PWGJE/AliAnalysisTaskJetCorePP.cxx +++ b/PWGJE/AliAnalysisTaskJetCorePP.cxx @@ -1102,7 +1102,7 @@ void AliAnalysisTaskJetCorePP::UserExec(Option_t *) cent = fESD->GetCentrality(); if(cent) centValue = cent->GetCentralityPercentile("V0M"); }else{ - centValue = aod->GetHeader()->GetCentrality(); + centValue = ((AliVAODHeader*)aod->GetHeader())->GetCentrality(); } if(fDebug) printf("centrality: %f\n", centValue); //Input events @@ -1120,7 +1120,10 @@ void AliAnalysisTaskJetCorePP::UserExec(Option_t *) //-----------------select disjunct event subsamples ---------------- if(!fIsKine){ //reconstructed data - Int_t eventnum = aod->GetHeader()->GetEventNumberESDFile(); + AliAODHeader * header = dynamic_cast(aod->GetHeader()); + if(!header) AliFatal("Not a standard AOD"); + + Int_t eventnum = header->GetEventNumberESDFile(); Int_t lastdigit = eventnum % 10; if(!(fEventNumberRangeLow<=lastdigit && lastdigit<=fEventNumberRangeHigh)){ fHistEvtSelection->Fill(5); diff --git a/PWGJE/AliAnalysisTaskJetServices.cxx b/PWGJE/AliAnalysisTaskJetServices.cxx index 868f8c98458..34d2db6037c 100644 --- a/PWGJE/AliAnalysisTaskJetServices.cxx +++ b/PWGJE/AliAnalysisTaskJetServices.cxx @@ -535,7 +535,7 @@ void AliAnalysisTaskJetServices::UserExec(Option_t */*option*/) Float_t cent = 0; if(fCollisionType==kPbPb){ - if(aod)cent = aod->GetHeader()->GetCentrality(); + if(aod)cent = ((AliVAODHeader*)aod->GetHeader())->GetCentrality(); if(fDebug)Printf("%s:%d %3.3f",(char*)__FILE__,__LINE__,cent); if(cent<0)cent = 101; } @@ -624,8 +624,8 @@ void AliAnalysisTaskJetServices::UserExec(Option_t */*option*/) Float_t zvtx = vtxAOD->GetZ(); Int_t iCl = GetEventClass(aod); AliAnalysisHelperJetTasks::EventClass(kTRUE,iCl); - Bool_t cand = aod->GetHeader()->GetOfflineTrigger()&fPhysicsSelectionFlag; - if(fDebug)Printf("%s:%d AOD selection %d %d",(char*)__FILE__,__LINE__,cand,aod->GetHeader()->GetOfflineTrigger()); + Bool_t cand = ((AliVAODHeader*)aod->GetHeader())->GetOfflineTrigger()&fPhysicsSelectionFlag; + if(fDebug)Printf("%s:%d AOD selection %d %d",(char*)__FILE__,__LINE__,cand,((AliVAODHeader*)aod->GetHeader())->GetOfflineTrigger()); fh2TriggerCount->Fill(0.,kAllTriggered); fh2TriggerCount->Fill(iCl,kAllTriggered); if(cand){ @@ -671,7 +671,7 @@ void AliAnalysisTaskJetServices::UserExec(Option_t */*option*/) TList recTracks; GetListOfTracks(&recTracks); CalculateReactionPlaneAngleVZERO(aod); - fRPAngle = aod->GetHeader()->GetEventplane(); + fRPAngle = ((AliVAODHeader*)aod->GetHeader())->GetEventplane(); fh1RP->Fill(fRPAngle); fh2RPCentrality->Fill(fCentrality,fRPAngle); fh2RPACentrality->Fill(fCentrality,fPsiVZEROA); @@ -1069,7 +1069,7 @@ Int_t AliAnalysisTaskJetServices::GetEventClass(AliESDEvent *esd){ Int_t AliAnalysisTaskJetServices::GetEventClass(AliAODEvent *aod){ if(fCollisionType==kPbPb){ - Float_t cent = aod->GetHeader()->GetCentrality(); + Float_t cent = ((AliVAODHeader*)aod->GetHeader())->GetCentrality(); if(cent>80||cent<0)return 5; if(cent>50)return 4; if(cent>30)return 3; diff --git a/PWGJE/AliAnalysisTaskJetSpectrum2.cxx b/PWGJE/AliAnalysisTaskJetSpectrum2.cxx index dd811e2a674..1863b73ef23 100644 --- a/PWGJE/AliAnalysisTaskJetSpectrum2.cxx +++ b/PWGJE/AliAnalysisTaskJetSpectrum2.cxx @@ -964,9 +964,11 @@ void AliAnalysisTaskJetSpectrum2::UserExec(Option_t */*option*/){ // CalculateReactionPlaneAngle(&recParticles); fRPAngle = 0; - if(fRPMethod==0)fRPAngle = aod->GetHeader()->GetEventplane(); + if(fRPMethod==0)fRPAngle = ((AliVAODHeader*)aod->GetHeader())->GetEventplane(); else if(fRPMethod==1||fRPMethod==2){ - fRPAngle = aod->GetHeader()->GetQTheta(fRPMethod); + AliAODHeader * aodheader = dynamic_cast(aod->GetHeader()); + if(!aodheader) AliFatal("Not a standard AOD"); + fRPAngle = aodheader->GetQTheta(fRPMethod); } fh1RP->Fill(fRPAngle); fh2RPCentrality->Fill(fCentrality,fRPAngle); @@ -1671,7 +1673,7 @@ Float_t AliAnalysisTaskJetSpectrum2::GetCentrality(){ if(!aod){ return 101; } - return aod->GetHeader()->GetCentrality(); + return ((AliVAODHeader*)aod->GetHeader())->GetCentrality(); } diff --git a/PWGJE/AliAnalysisTaskPartonDisc.cxx b/PWGJE/AliAnalysisTaskPartonDisc.cxx index 03b99064952..a16378504b8 100644 --- a/PWGJE/AliAnalysisTaskPartonDisc.cxx +++ b/PWGJE/AliAnalysisTaskPartonDisc.cxx @@ -1174,7 +1174,8 @@ void AliAnalysisTaskPartonDisc::UserExec(Option_t *) if(fIsHIevent) { - AliAODHeader *aodHeader = fAOD->GetHeader(); + AliAODHeader *aodHeader = dynamic_cast(fAOD->GetHeader()); + if(!aodHeader) AliFatal("Not a standard AOD"); fEventCent = aodHeader->GetCentrality(); } @@ -1190,7 +1191,10 @@ void AliAnalysisTaskPartonDisc::UserExec(Option_t *) fJetAcceptance = 0.5 - fIncExcR; // if the increase is 0.1 -> only jets within |eta|<0.4 // First test of reference multiplicity - Int_t refMultiplicity = fAOD->GetHeader()->GetRefMultiplicity(); + AliAODHeader * header = dynamic_cast(fAOD->GetHeader()); + if(!header) AliFatal("Not a standard AOD"); + + Int_t refMultiplicity = header->GetRefMultiplicity(); fRefMult->Fill(refMultiplicity); // Multiplicity from V0 (V0A+V0C) diff --git a/PWGJE/AliPWG4HighPtTrackQA.cxx b/PWGJE/AliPWG4HighPtTrackQA.cxx index 2a99def4b6e..7f195f2d82a 100644 --- a/PWGJE/AliPWG4HighPtTrackQA.cxx +++ b/PWGJE/AliPWG4HighPtTrackQA.cxx @@ -842,8 +842,8 @@ Bool_t AliPWG4HighPtTrackQA::SelectEvent() } } else if(fDataType==kAOD) { - if(dynamic_cast(fEvent)->GetHeader()->GetCentrality()) - cent = dynamic_cast(fEvent)->GetHeader()->GetCentrality(); + if(((AliVAODHeader*)dynamic_cast(fEvent)->GetHeader())->GetCentrality()) + cent = ((AliVAODHeader*)dynamic_cast(fEvent)->GetHeader())->GetCentrality(); } if(cent>90.) { fNEventReject->Fill("cent>90",1); @@ -901,7 +901,7 @@ Int_t AliPWG4HighPtTrackQA::CalculateCentrality(const AliAODEvent *aod) // if(!aod) return 5; - Float_t cent = aod->GetHeader()->GetCentrality(); + Float_t cent = ((AliVAODHeader*)aod->GetHeader())->GetCentrality(); if(fDebug>3) printf("centrality: %f\n",cent); return GetCentralityClass(cent); diff --git a/PWGJE/StrangenessInJets/AliAnalysisTaskJetChem.cxx b/PWGJE/StrangenessInJets/AliAnalysisTaskJetChem.cxx index 0b436825d3d..d40f61e543d 100644 --- a/PWGJE/StrangenessInJets/AliAnalysisTaskJetChem.cxx +++ b/PWGJE/StrangenessInJets/AliAnalysisTaskJetChem.cxx @@ -2178,7 +2178,7 @@ void AliAnalysisTaskJetChem::UserExec(Option_t *) if(handler && handler->InheritsFrom("AliAODInputHandler")){ - centPercent = fAOD->GetHeader()->GetCentrality(); + centPercent = ((AliVAODHeader*)fAOD->GetHeader())->GetCentrality(); cl = 1; //std::cout<<"centPercent: "<GetHeader()->GetCentrality(); // event centrality - fdCentrality = fAODIn->GetHeader()->GetCentralityP()->GetCentralityPercentile("V0M"); // event centrality +// fdCentrality = ((AliVAODHeader*)fAODIn->GetHeader())->GetCentrality(); // event centrality + fdCentrality = ((AliVAODHeader*)fAODIn->GetHeader())->GetCentralityP()->GetCentralityPercentile("V0M"); // event centrality if(!fbIsPbPb) fdCentrality = 0.; Int_t iCentIndex = GetCentralityBinIndex(fdCentrality); // get index of centrality bin @@ -3118,7 +3118,7 @@ Bool_t AliAnalysisTaskV0sInJets::IsSelectedForJets(AliAODEvent* fAOD, Double_t d return kFALSE; Double_t centrality; // centrality = fAOD->GetHeader()->GetCentrality(); - centrality = fAOD->GetHeader()->GetCentralityP()->GetCentralityPercentile("V0M"); + centrality = ((AliVAODHeader*)fAOD->GetHeader())->GetCentralityP()->GetCentralityPercentile("V0M"); if(fbIsPbPb) { if(centrality < 0) diff --git a/PWGJE/UserTasks/AliAnalysisTaskIDFFTCF.cxx b/PWGJE/UserTasks/AliAnalysisTaskIDFFTCF.cxx index 50c8167970d..2183c9132a5 100644 --- a/PWGJE/UserTasks/AliAnalysisTaskIDFFTCF.cxx +++ b/PWGJE/UserTasks/AliAnalysisTaskIDFFTCF.cxx @@ -1650,7 +1650,7 @@ void AliAnalysisTaskIDFFTCF::UserExec(Option_t *) Int_t cl = 0; if(handler->InheritsFrom("AliAODInputHandler")){ // since it is not supported by the helper task define own classes - centPercent = fAOD->GetHeader()->GetCentrality(); + centPercent = ((AliVAODHeader*)fAOD->GetHeader())->GetCentrality(); cl = 1; if(centPercent>10) cl = 2; if(centPercent>30) cl = 3; diff --git a/PWGJE/UserTasks/AliAnalysisTaskIDFragmentationFunction.cxx b/PWGJE/UserTasks/AliAnalysisTaskIDFragmentationFunction.cxx index 477e3d07559..672c81376c1 100644 --- a/PWGJE/UserTasks/AliAnalysisTaskIDFragmentationFunction.cxx +++ b/PWGJE/UserTasks/AliAnalysisTaskIDFragmentationFunction.cxx @@ -2580,7 +2580,7 @@ void AliAnalysisTaskIDFragmentationFunction::UserExec(Option_t *) Int_t cl = 0; if(handler->InheritsFrom("AliAODInputHandler")){ // since it is not supported by the helper task define own classes - centPercent = fAOD->GetHeader()->GetCentrality(); + centPercent = ((AliVAODHeader*)fAOD->GetHeader())->GetCentrality(); cl = 1; if(centPercent>10) cl = 2; if(centPercent>30) cl = 3; @@ -2601,8 +2601,12 @@ void AliAnalysisTaskIDFragmentationFunction::UserExec(Option_t *) } // Retrieve reference multiplicities in |eta|<0.8 and <0.5 - const Int_t refMult5 = fAOD->GetHeader()->GetRefMultiplicityComb05(); - const Int_t refMult8 = fAOD->GetHeader()->GetRefMultiplicityComb08(); + + AliAODHeader * header = dynamic_cast(fAOD->GetHeader()); + if(!header) AliFatal("Not a standard AOD"); + + const Int_t refMult5 = header->GetRefMultiplicityComb05(); + const Int_t refMult8 = header->GetRefMultiplicityComb08(); const Double_t centPercentPP = fAnaUtils->GetMultiplicityPercentile(fAOD, "V0M"); @@ -4005,7 +4009,8 @@ Int_t AliAnalysisTaskIDFragmentationFunction::GetListOfTracks(TList *list, Int_t // all rec. tracks, esd filter mask, eta range for(Int_t it=0; itGetNumberOfTracks(); ++it){ - AliAODTrack *tr = fAOD->GetTrack(it); + AliAODTrack *tr = dynamic_cast(fAOD->GetTrack(it)); + if(!tr) AliFatal("Not a standard AOD"); if(type == kTrackAODCuts || type==kTrackAODQualityCuts || type==kTrackAODExtraCuts){ diff --git a/PWGJE/UserTasks/AliAnalysisTaskJetAntenna.cxx b/PWGJE/UserTasks/AliAnalysisTaskJetAntenna.cxx index cfb9ba107f4..03830cdcf82 100644 --- a/PWGJE/UserTasks/AliAnalysisTaskJetAntenna.cxx +++ b/PWGJE/UserTasks/AliAnalysisTaskJetAntenna.cxx @@ -396,7 +396,7 @@ void AliAnalysisTaskJetAntenna::UserExec(Option_t *) if(fIsPbPb){ if(fESD) {cent = fESD->GetCentrality(); if(cent) centValue = cent->GetCentralityPercentile("V0M");} - else centValue=aod->GetHeader()->GetCentrality(); + else centValue=((AliVAODHeader*)aod->GetHeader())->GetCentrality(); if(fDebug) printf("centrality: %f\n", centValue); if (centValue < fCentMin || centValue > fCentMax){ diff --git a/PWGJE/UserTasks/AliAnalysisTaskJetShape.cxx b/PWGJE/UserTasks/AliAnalysisTaskJetShape.cxx index 537e6d505d3..77829e3dc06 100644 --- a/PWGJE/UserTasks/AliAnalysisTaskJetShape.cxx +++ b/PWGJE/UserTasks/AliAnalysisTaskJetShape.cxx @@ -1929,7 +1929,7 @@ if(!aodH){ if(fESD) {cent = fESD->GetCentrality(); if(cent) centrality = cent->GetCentralityPercentile("V0M");} - else centrality=aodE->GetHeader()->GetCentrality(); + else centrality=((AliVAODHeader*)aodE->GetHeader())->GetCentrality(); if(!fkIsPbPb) { diff --git a/PWGJE/UserTasks/AliAnalysisTaskPPJetSpectra.cxx b/PWGJE/UserTasks/AliAnalysisTaskPPJetSpectra.cxx index 2796b36f849..bbabf17abd0 100755 --- a/PWGJE/UserTasks/AliAnalysisTaskPPJetSpectra.cxx +++ b/PWGJE/UserTasks/AliAnalysisTaskPPJetSpectra.cxx @@ -501,7 +501,7 @@ Bool_t AliAnalysisTaskPPJetSpectra::EventSelection(Double_t evtContainer[6]) { Float_t centrality = -1; if(fEventClass > 0) { - if(handler->InheritsFrom("AliAODHandler")) centrality = fAODIn->GetHeader()->GetCentrality(); + if(handler->InheritsFrom("AliAODHandler")) centrality = ((AliVAODHeader*)fAODIn->GetHeader())->GetCentrality(); else if(fESD) centrality = fESD->GetCentrality()->GetCentralityPercentile("V0M"); else centrality = AliAnalysisHelperJetTasks::EventClass(); } diff --git a/PWGLF/FORWARD/analysis2/AliBasedNdetaTask.cxx b/PWGLF/FORWARD/analysis2/AliBasedNdetaTask.cxx index 8d22ebffad3..e6044d6d927 100644 --- a/PWGLF/FORWARD/analysis2/AliBasedNdetaTask.cxx +++ b/PWGLF/FORWARD/analysis2/AliBasedNdetaTask.cxx @@ -447,7 +447,8 @@ AliBasedNdetaTask::Event(AliAODEvent& aod) if (HasCentrality()) { Double_t cent = forward->GetCentrality(); if (!fCentMethod.IsNull()) { - AliAODHeader* hdr = aod.GetHeader(); + AliAODHeader* hdr = dynamic_cast(aod.GetHeader()); + if(!hdr) AliFatal("Not a standard AOD"); if (hdr) { AliCentrality* cP = hdr->GetCentralityP(); if (cP) { diff --git a/PWGLF/FORWARD/analysis2/AliCopyHeaderTask.cxx b/PWGLF/FORWARD/analysis2/AliCopyHeaderTask.cxx index ffc92a5118f..f476b7a5ba2 100644 --- a/PWGLF/FORWARD/analysis2/AliCopyHeaderTask.cxx +++ b/PWGLF/FORWARD/analysis2/AliCopyHeaderTask.cxx @@ -46,7 +46,8 @@ AliCopyHeaderTask::UserExec(Option_t*) LoadBranches(); - AliAODHeader* aodHeader = aod->GetHeader(); + AliAODHeader* aodHeader = dynamic_cast(aod->GetHeader()); + if(!aodHeader) AliFatal("Not a standard AOD"); if (!aodHeader) { AliWarning("Missing AOD header"); aodHeader = new AliAODHeader(esd->GetRunNumber(), diff --git a/PWGLF/RESONANCES/AliRsnMiniAnalysisTask.cxx b/PWGLF/RESONANCES/AliRsnMiniAnalysisTask.cxx index 4caaf37de35..7f5118313c8 100644 --- a/PWGLF/RESONANCES/AliRsnMiniAnalysisTask.cxx +++ b/PWGLF/RESONANCES/AliRsnMiniAnalysisTask.cxx @@ -842,7 +842,7 @@ Double_t AliRsnMiniAnalysisTask::ComputeAngle() plane = fInputEvent->GetEventplane(); else if (fInputEvent->InheritsFrom(AliAODEvent::Class())) { AliAODEvent *aodEvent = (AliAODEvent *)fInputEvent; - plane = aodEvent->GetHeader()->GetEventplaneP(); + plane = ((AliVAODHeader*)aodEvent->GetHeader())->GetEventplaneP(); } if (plane) diff --git a/PWGLF/SPECTRA/Nuclei/B2/AliAnalysisTaskB2AOD.cxx b/PWGLF/SPECTRA/Nuclei/B2/AliAnalysisTaskB2AOD.cxx index b76c59ed82e..4e705913ecb 100644 --- a/PWGLF/SPECTRA/Nuclei/B2/AliAnalysisTaskB2AOD.cxx +++ b/PWGLF/SPECTRA/Nuclei/B2/AliAnalysisTaskB2AOD.cxx @@ -212,8 +212,10 @@ void AliAnalysisTaskB2AOD::UserExec(Option_t* ) if(fLnID == 0) AliFatal("PID not set"); // multiplicity and centrality - - fNtrk = (fMaxEta > 0.5) ? fAODevent->GetHeader()->GetRefMultiplicityComb08() : fAODevent->GetHeader()->GetRefMultiplicityComb05(); + AliAODHeader * header = dynamic_cast(fAODevent->GetHeader()); + if(!header) AliFatal("Not a standard AOD"); + + fNtrk = (fMaxEta > 0.5) ? header->GetRefMultiplicityComb08() : header->GetRefMultiplicityComb05(); if(fSimulation) fNch = this->GetChargedMultiplicity(fMaxEta); @@ -226,7 +228,7 @@ void AliAnalysisTaskB2AOD::UserExec(Option_t* ) if(fHeavyIons) { - Double_t centrality = fAODevent->GetHeader()->GetCentrality(); + Double_t centrality = ((AliVAODHeader*)fAODevent->GetHeader())->GetCentrality(); fCentTriggerFired = (centrality >= fMinCentrality) && (centrality < fMaxCentrality); } diff --git a/PWGLF/SPECTRA/PiKaPr/TestAOD/AliSpectraAODEventCuts.cxx b/PWGLF/SPECTRA/PiKaPr/TestAOD/AliSpectraAODEventCuts.cxx index 8b8ccb8b57a..6bce8620bc6 100644 --- a/PWGLF/SPECTRA/PiKaPr/TestAOD/AliSpectraAODEventCuts.cxx +++ b/PWGLF/SPECTRA/PiKaPr/TestAOD/AliSpectraAODEventCuts.cxx @@ -186,8 +186,9 @@ Bool_t AliSpectraAODEventCuts::IsSelected(AliAODEvent * aod,AliSpectraAODTrackCu } Int_t Nch=0; for (Int_t iTracks = 0; iTracks < fAOD->GetNumberOfTracks(); iTracks++) { - AliAODTrack* track = fAOD->GetTrack(iTracks); - if (!fTrackCuts->IsSelected(track,kFALSE)) continue; + AliAODTrack* track = dynamic_cast(fAOD->GetTrack(iTracks)); + if(!track) AliFatal("Not a standard AOD"); + if (!fTrackCuts->IsSelected(track,kFALSE)) continue; ((TH1F*)fOutput->FindObject("fHistoEtaBefSel"))->Fill(track->Eta()); if(fIsSelected){ ((TH1F*)fOutput->FindObject("fHistoEtaAftSel"))->Fill(track->Eta()); @@ -237,7 +238,8 @@ Bool_t AliSpectraAODEventCuts::CheckMultiplicityCut() // FIXME: why this is not tracket in the track stats histos? Int_t Ncharged=0; for (Int_t iTracks = 0; iTracks < fAOD->GetNumberOfTracks(); iTracks++){ - AliAODTrack* track = fAOD->GetTrack(iTracks); + AliAODTrack* track = dynamic_cast(fAOD->GetTrack(iTracks)); + if(!track) AliFatal("Not a standard AOD"); if (!fTrackCuts->IsSelected(track,kFALSE)) continue; Ncharged++; } @@ -401,7 +403,8 @@ Double_t AliSpectraAODEventCuts::CalculateQVectorTPC(Double_t etaMin,Double_t et Double_t Qx2 = 0, Qy2 = 0; Int_t mult = 0; for(Int_t iT = 0; iT < fAOD->GetNumberOfTracks(); iT++) { - AliAODTrack* aodTrack = fAOD->GetTrack(iT); + AliAODTrack* aodTrack = dynamic_cast(fAOD->GetTrack(iT)); + if(!aodTrack) AliFatal("Not a standard AOD"); if (!aodTrack->TestFilterBit(128)) continue; //FIXME track type hard coded -> TPC only constrained to the vertex if (aodTrack->Eta() Eta() > etaMax)continue; mult++; @@ -458,7 +461,7 @@ Short_t AliSpectraAODEventCuts::GetCentrCode(AliVEvent* ev) AliCentrality* centrality = 0; AliAODEvent* aod = (AliAODEvent*)ev; - centrality = aod->GetHeader()->GetCentralityP(); + centrality = ((AliVAODHeader*)aod->GetHeader())->GetCentralityP(); Float_t centV0 = centrality->GetCentralityPercentile("V0M"); Float_t centTrk = centrality->GetCentralityPercentile("TRK"); diff --git a/PWGLF/SPECTRA/PiKaPr/TestAOD/AliSpectraBothEventCuts.cxx b/PWGLF/SPECTRA/PiKaPr/TestAOD/AliSpectraBothEventCuts.cxx index fa3e54c02ad..87b75e53400 100644 --- a/PWGLF/SPECTRA/PiKaPr/TestAOD/AliSpectraBothEventCuts.cxx +++ b/PWGLF/SPECTRA/PiKaPr/TestAOD/AliSpectraBothEventCuts.cxx @@ -386,10 +386,13 @@ Bool_t AliSpectraBothEventCuts::CheckMultiplicityCut() { AliAODEvent* aodevent=0x0; aodevent=dynamic_cast(fAOD); + AliAODHeader * header = dynamic_cast(aodevent->GetHeader()); + if(!header) AliFatal("Not a standard AOD"); + if(TMath::Abs(0.8-fetarangeofmultiplicitycut)<0.1) - Ncharged=aodevent->GetHeader()->GetRefMultiplicityComb08(); + Ncharged=header->GetRefMultiplicityComb08(); else if (TMath::Abs(0.5-fetarangeofmultiplicitycut)<0.1) - Ncharged=aodevent->GetHeader()->GetRefMultiplicityComb05(); + Ncharged=header->GetRefMultiplicityComb05(); else Ncharged=-1; } diff --git a/PWGLF/STRANGENESS/Correlations/AliAnalysisTaskV0ChCorrelations.cxx b/PWGLF/STRANGENESS/Correlations/AliAnalysisTaskV0ChCorrelations.cxx index 5cf673d1c07..d7ae9652654 100644 --- a/PWGLF/STRANGENESS/Correlations/AliAnalysisTaskV0ChCorrelations.cxx +++ b/PWGLF/STRANGENESS/Correlations/AliAnalysisTaskV0ChCorrelations.cxx @@ -300,7 +300,7 @@ void AliAnalysisTaskV0ChCorrelations::UserExec(Option_t *) // Centrality definition Double_t lCent = 0.0; AliCentrality *centralityObj = 0; - centralityObj = aod->GetHeader()->GetCentralityP(); + centralityObj = ((AliVAODHeader*)aod->GetHeader())->GetCentralityP(); lCent = centralityObj->GetCentralityPercentile("V0M"); if ((lCent < 0.)||(lCent > 90.)) return; fHistCentVtx->Fill(lCent,lPVz); diff --git a/PWGLF/STRANGENESS/Hypernuclei/AliAODMCNuclExReplicator.cxx b/PWGLF/STRANGENESS/Hypernuclei/AliAODMCNuclExReplicator.cxx index 5dea00dc673..31caac331ea 100644 --- a/PWGLF/STRANGENESS/Hypernuclei/AliAODMCNuclExReplicator.cxx +++ b/PWGLF/STRANGENESS/Hypernuclei/AliAODMCNuclExReplicator.cxx @@ -266,7 +266,9 @@ void AliAODMCNuclExReplicator::ReplicateAndFilter(const AliAODEvent& source) if (fReplicateHeader) { - *fHeader = *(source.GetHeader()); + AliAODHeader * header = dynamic_cast(source.GetHeader()); + if(!header) AliFatal("Not a standard AOD"); + *fHeader = *(header); } fVertices->Clear("C"); diff --git a/PWGLF/STRANGENESS/Hypernuclei/AliAODNuclExReplicator.cxx b/PWGLF/STRANGENESS/Hypernuclei/AliAODNuclExReplicator.cxx index cd338d9062c..bfc853d6ed2 100644 --- a/PWGLF/STRANGENESS/Hypernuclei/AliAODNuclExReplicator.cxx +++ b/PWGLF/STRANGENESS/Hypernuclei/AliAODNuclExReplicator.cxx @@ -438,7 +438,9 @@ void AliAODNuclExReplicator::ReplicateAndFilter(const AliAODEvent& source) if (fReplicateHeader) { - *fHeader = *(source.GetHeader()); + AliAODHeader * header = dynamic_cast(source.GetHeader()); + if(!header) AliFatal("Not a standard AOD"); + *fHeader = *(header); } fVertices->Clear("C"); diff --git a/PWGLF/STRANGENESS/Hypernuclei/AliAnalysisTaskReadNuclexAOD.cxx b/PWGLF/STRANGENESS/Hypernuclei/AliAnalysisTaskReadNuclexAOD.cxx index 61411412b70..2529d11baca 100644 --- a/PWGLF/STRANGENESS/Hypernuclei/AliAnalysisTaskReadNuclexAOD.cxx +++ b/PWGLF/STRANGENESS/Hypernuclei/AliAnalysisTaskReadNuclexAOD.cxx @@ -638,7 +638,7 @@ void AliAnalysisTaskReadNuclexAOD::UserExec(Option_t *) centrality = lAODevent->GetCentrality(); Float_t percentile=centrality->GetCentralityPercentile("V0M"); - Float_t refMult = lAODevent->GetHeader()->GetRefMultiplicity(); + Float_t refMult = ((AliVAODHeader*)lAODevent->GetHeader())->GetRefMultiplicity(); fHistTrackMultiplicity->Fill(refMult,percentile); //tracce per evento diff --git a/PWGPP/EvTrkSelection/AliCFSingleTrackEfficiencyTask.cxx b/PWGPP/EvTrkSelection/AliCFSingleTrackEfficiencyTask.cxx index eebed12500e..3dbf229a59b 100644 --- a/PWGPP/EvTrkSelection/AliCFSingleTrackEfficiencyTask.cxx +++ b/PWGPP/EvTrkSelection/AliCFSingleTrackEfficiencyTask.cxx @@ -710,7 +710,8 @@ Double_t AliCFSingleTrackEfficiencyTask::GetCentrality() if(isAOD) { AliAODEvent* aodEvent = dynamic_cast(fInputEvent); if(!aodEvent) return cent; - AliAODHeader* header = aodEvent->GetHeader(); + AliAODHeader* header = dynamic_cast(aodEvent->GetHeader()); + if(!header) AliFatal("Not a standard AOD"); if(!header) return cent; AliCentrality *centrality = header->GetCentralityP(); if(!centrality) return cent; diff --git a/PWGPP/ITS/AliAnalysisTaskSEImpParRes.cxx b/PWGPP/ITS/AliAnalysisTaskSEImpParRes.cxx index 5262ad51f67..e40658645e8 100644 --- a/PWGPP/ITS/AliAnalysisTaskSEImpParRes.cxx +++ b/PWGPP/ITS/AliAnalysisTaskSEImpParRes.cxx @@ -1163,7 +1163,7 @@ void AliAnalysisTaskSEImpParRes::UserExec(Option_t */*option*/) if(fIsAOD){ Int_t nclsITS = 0; runNumber=((AliAODEvent*)event)->GetRunNumber(); - nclsITS = ((AliAODEvent*)event)->GetHeader()->GetNumberOfITSClusters(1); + nclsITS = ((AliVAODHeader*)((AliAODEvent*)event)->GetHeader())->GetNumberOfITSClusters(1); if(nclsITSfMaxMult) return; firedTriggerClasses=((AliAODEvent*)event)->GetFiredTriggerClasses(); if(!firedTriggerClasses.Contains(fTriggerClass.Data())) return; diff --git a/PWGUD/UPC/AliAnalysisTaskUpcK0sK0s.cxx b/PWGUD/UPC/AliAnalysisTaskUpcK0sK0s.cxx index 0bbf707c2bb..6b985b546e3 100644 --- a/PWGUD/UPC/AliAnalysisTaskUpcK0sK0s.cxx +++ b/PWGUD/UPC/AliAnalysisTaskUpcK0sK0s.cxx @@ -314,8 +314,11 @@ void AliAnalysisTaskUpcK0sK0s::RunAODtree() if(!aod) return; //input data - fDataFilnam = aod->GetHeader()->GetESDFileName(); - fEvtNum = aod->GetHeader()->GetEventNumberESDFile(); + AliAODHeader * header = dynamic_cast(aod->GetHeader()); +if(!header) AliFatal("Not a standard AOD"); + + fDataFilnam = header->GetESDFileName(); + fEvtNum = header->GetEventNumberESDFile(); fRunNum = aod->GetRunNumber(); @@ -334,8 +337,8 @@ void AliAnalysisTaskUpcK0sK0s::RunAODtree() if(!isTriggered ) return; //trigger inputs - fL0inputs = aod->GetHeader()->GetL0TriggerInputs(); - fL1inputs = aod->GetHeader()->GetL1TriggerInputs(); + fL0inputs = ((AliVAODHeader*)aod->GetHeader())->GetL0TriggerInputs(); + fL1inputs = ((AliVAODHeader*)aod->GetHeader())->GetL1TriggerInputs(); //Event identification fPerNum = aod ->GetPeriodNumber(); diff --git a/STEER/AOD/AliNanoAODHeader.h b/STEER/AOD/AliNanoAODHeader.h index 60e618df599..dd396f67a32 100644 --- a/STEER/AOD/AliNanoAODHeader.h +++ b/STEER/AOD/AliNanoAODHeader.h @@ -63,6 +63,7 @@ public: // FIXME: use dynamic cast in AliAODEVent? virtual AliCentrality* GetCentralityP() const {AliError("Not Implemented"); return 0;}; virtual AliEventplane* GetEventplaneP() const {AliError("Not Implemented"); return 0;}; + virtual Double_t GetEventplane() const {AliError("Not Implemented"); return 0;}; virtual const Float_t* GetVZEROEqFactors() const {AliError("Not Implemented"); return 0;}; virtual Float_t GetVZEROEqFactors(Int_t /* i */) const {AliError("Not Implemented"); return 0;}; virtual void SetVZEROEqFactors(const Float_t* /*factors*/) {AliError("Not Implemented"); } @@ -74,9 +75,14 @@ public: virtual void Clear(Option_t * opt) ; + virtual Int_t GetIRInt2ClosestInteractionMap() const {AliError("Not Implemented"); return 0;}; + virtual Int_t GetIRInt1ClosestInteractionMap(Int_t /*gap = 3*/) const {AliError("Not Implemented"); return 0;}; + + + virtual Int_t GetRefMultiplicity() const { AliError("Not Impletented"); return 0; } Double_t GetMagneticField() const { return GetVar(1); } - Double_t GetCentrality (const char *estimator = "V0M") { return GetVar(0);} + Double_t GetCentrality (/*estimator = "V0M"*/) const { return GetVar(0);} ClassDef(AliNanoAODHeader, 1) diff --git a/STEER/STEERBase/AliVAODHeader.h b/STEER/STEERBase/AliVAODHeader.h index 3ccce648e15..64949fb738c 100644 --- a/STEER/STEERBase/AliVAODHeader.h +++ b/STEER/STEERBase/AliVAODHeader.h @@ -71,10 +71,13 @@ class AliVAODHeader : public AliVHeader { // FIXME: use dynamic cast in AliAODEVent? virtual AliCentrality* GetCentralityP() const = 0; virtual AliEventplane* GetEventplaneP() const = 0; + virtual Double_t GetCentrality () const = 0; virtual const Float_t* GetVZEROEqFactors() const = 0; virtual Float_t GetVZEROEqFactors(Int_t /* i */) const = 0; virtual void SetVZEROEqFactors(const Float_t* /*factors*/) = 0; - virtual UInt_t GetOfflineTrigger() = 0; + virtual UInt_t GetOfflineTrigger() = 0; + virtual Int_t GetRefMultiplicity() const =0; + virtual Double_t GetEventplane() const =0; }; diff --git a/STEER/STEERBase/AliVHeader.h b/STEER/STEERBase/AliVHeader.h index 20c2c6bacef..6b8a50a4643 100644 --- a/STEER/STEERBase/AliVHeader.h +++ b/STEER/STEERBase/AliVHeader.h @@ -42,6 +42,9 @@ class AliVHeader : public TNamed { virtual TBits GetIRInt2InteractionMap() const { static TBits dummy; return dummy; } virtual TBits GetIRInt1InteractionMap() const { static TBits dummy; return dummy; } + virtual Int_t GetIRInt2ClosestInteractionMap() const {return 0;} + virtual Int_t GetIRInt1ClosestInteractionMap(Int_t =3/*gap = 3*/) const {return 0;} + virtual void Print(Option_t* option = "") const = 0;