#include "AliAODEvent.h"\r
#include "AliAODTrack.h"\r
#include "AliAODInputHandler.h"\r
+#include "AliCollisionGeometry.h"\r
#include "AliGenEventHeader.h"\r
-#include "AliGenHijingEventHeader.h"\r
#include "AliMCEventHandler.h"\r
#include "AliMCEvent.h"\r
#include "AliStack.h"\r
\r
// take only events inside centrality class\r
if((fImpactParameterMin < gCentrality) && (fImpactParameterMax > gCentrality)){\r
+ fHistEventStats->Fill(5,gCentrality); //events with correct centrality\r
return gCentrality; \r
}//centrality class\r
}//Vz cut\r
}//ESD\r
else if(gAnalysisLevel == "MC"){\r
Double_t gImpactParameter = 0.;\r
- AliGenHijingEventHeader* headerH = dynamic_cast<AliGenHijingEventHeader*>(dynamic_cast<AliMCEvent*>(event)->GenEventHeader());\r
- if(headerH){\r
- gImpactParameter = headerH->ImpactParameter();\r
- gCentrality = gImpactParameter;\r
- }//MC header\r
+ if(dynamic_cast<AliMCEvent*>(event)){\r
+ AliCollisionGeometry* headerH = dynamic_cast<AliCollisionGeometry*>(dynamic_cast<AliMCEvent*>(event)->GenEventHeader()); \r
+ if(headerH){\r
+ gImpactParameter = headerH->ImpactParameter();\r
+ gCentrality = gImpactParameter;\r
+ }//MC header\r
+ }//MC event cast\r
}//MC\r
else{\r
gCentrality = -1.;\r
\r
//MC: from reaction plane\r
if(gAnalysisLevel == "MC"){\r
- \r
- AliGenHijingEventHeader* headerH = dynamic_cast<AliGenHijingEventHeader*>(dynamic_cast<AliMCEvent*>(event)->GenEventHeader());\r
- if (headerH) {\r
- gReactionPlane = headerH->ReactionPlaneAngle();\r
- //gReactionPlane *= TMath::RadToDeg();\r
- }\r
+ if(dynamic_cast<AliMCEvent*>(event)){\r
+ AliCollisionGeometry* headerH = dynamic_cast<AliCollisionGeometry*>(dynamic_cast<AliMCEvent*>(event)->GenEventHeader()); \r
+ if (headerH) {\r
+ gReactionPlane = headerH->ReactionPlaneAngle();\r
+ //gReactionPlane *= TMath::RadToDeg();\r
+ }//MC header\r
+ }//MC event cast\r
}//MC\r
\r
// AOD,ESD,ESDMC: from VZERO Event Plane\r