From 74b0bf8547e39fcb52b6ad978f5f02a48ef28f3f Mon Sep 17 00:00:00 2001 From: cvetan Date: Thu, 21 Feb 2013 20:52:36 +0000 Subject: [PATCH] Fix for true centrality estimator with hijing npart (Alberica) --- ANALYSIS/AliCentralitySelectionTask.cxx | 12 +++++------- 1 file changed, 5 insertions(+), 7 deletions(-) diff --git a/ANALYSIS/AliCentralitySelectionTask.cxx b/ANALYSIS/AliCentralitySelectionTask.cxx index 8b2101f571f..24ff26d55be 100644 --- a/ANALYSIS/AliCentralitySelectionTask.cxx +++ b/ANALYSIS/AliCentralitySelectionTask.cxx @@ -1222,24 +1222,22 @@ void AliCentralitySelectionTask::UserExec(Option_t */*option*/) AliStack* stack=0; AliMCEvent* mcEvent=0; if (fIsMCInput && eventHandler && (mcEvent=eventHandler->MCEvent()) && (stack=mcEvent->Stack())) { - AliGenHijingEventHeader* hHijing=0; AliGenDPMjetEventHeader* dpmHeader=0; AliGenEventHeader* mcGenH = mcEvent->GenEventHeader(); - if (mcGenH->InheritsFrom(AliGenHijingEventHeader::Class())) { - hHijing = (AliGenHijingEventHeader*)mcGenH; - if(hHijing) Npart = hHijing->ProjectileParticipants()+hHijing->TargetParticipants(); - } + if (mcGenH->InheritsFrom(AliGenHijingEventHeader::Class())) + hHijing = (AliGenHijingEventHeader*)mcGenH; else if (mcGenH->InheritsFrom(AliGenCocktailEventHeader::Class())) { TList* headers = ((AliGenCocktailEventHeader*)mcGenH)->GetHeaders(); hHijing = dynamic_cast(headers->FindObject("Hijing")); - if(hHijing) Npart = hHijing->ProjectileParticipants()+hHijing->TargetParticipants(); + if (!hHijing) hHijing = dynamic_cast(headers->FindObject("Hijing pPb_0")); } else if (mcGenH->InheritsFrom(AliGenDPMjetEventHeader::Class())) { dpmHeader = (AliGenDPMjetEventHeader*)mcGenH; - if(dpmHeader) Npart = dpmHeader->ProjectileParticipants()+ dpmHeader->TargetParticipants(); } + if(hHijing) Npart = hHijing->ProjectileParticipants()+hHijing->TargetParticipants(); + if(dpmHeader) Npart = dpmHeader->ProjectileParticipants()+ dpmHeader->TargetParticipants(); } } -- 2.39.3