}
ClassImp(AliAnalysisTaskPi0DiffEfficiency)
+//________________________________________________________________________
+AliAnalysisTaskPi0DiffEfficiency::AliAnalysisTaskPi0DiffEfficiency()
+: AliAnalysisTaskSE(),
+ fStack(0),
+ fOutputContainer(0),
+ fPHOSEvent(0),
+ fPHOSEvent1(0),
+ fPHOSEvent2(0),
+ fPHOSCalibData(0),
+ fNonLinCorr(0),
+ fRPfull(0),
+ fRPA(0),
+ fRPC(0),
+ fRPFar(0),
+ fRPAFar(0),
+ fRPCFar(0),
+ fCentrality(0),
+ fCenBin(0),
+ fPHOSGeo(0),
+ fEventCounter(0)
+{
+ // Constructor //should not be used
+ for(Int_t i=0;i<1;i++){
+ for(Int_t j=0;j<10;j++)
+ for(Int_t k=0;k<11;k++)
+ fPHOSEvents[i][j][k]=0 ;
+ }
+ for(Int_t mod=0; mod<6; mod++) fPHOSBadMap[mod]=0x0 ;
+}
//________________________________________________________________________
AliAnalysisTaskPi0DiffEfficiency::AliAnalysisTaskPi0DiffEfficiency(const char *name)
fCentrality(0),
fCenBin(0),
fPHOSGeo(0),
- fEventCounter(0)
+ fEventCounter(0)
{
// Constructor
for(Int_t i=0;i<1;i++){
}
-
+//________________________________________________________________________
+AliAnalysisTaskPi0DiffEfficiency::AliAnalysisTaskPi0DiffEfficiency(const AliAnalysisTaskPi0DiffEfficiency& a):
+AliAnalysisTaskSE(a),
+ fStack(0),
+ fOutputContainer(0),
+ fPHOSEvent(0),
+ fPHOSEvent1(0),
+ fPHOSEvent2(0),
+ fPHOSCalibData(0),
+ fNonLinCorr(0),
+ fRPfull(0),
+ fRPA(0),
+ fRPC(0),
+ fRPFar(0),
+ fRPAFar(0),
+ fRPCFar(0),
+ fCentrality(0),
+ fCenBin(0),
+ fPHOSGeo(0),
+ fEventCounter(0)
+{ // not implemented
+ for(Int_t mod=0; mod<6; mod++) fPHOSBadMap[mod]=0x0 ;
+}
//________________________________________________________________________
void AliAnalysisTaskPi0DiffEfficiency::UserCreateOutputObjects()
{
}
FillHistogram("hSelEvents",1.5) ;
- AliAODHeader *header = event->GetHeader() ;
+ AliAODHeader *header = dynamic_cast<AliAODHeader*>(event->GetHeader()) ;
+ if(!header) AliFatal("Not a standard AOD");
// Checks if we have a primary vertex
// Get primary vertices form ESD
6.58365e-01*5.91917e-01*5.91917e-01/((pt-9.61306e-01)*(pt-9.61306e-01)+5.91917e-01*5.91917e-01)+1.59219);
Double_t sz=TMath::Min(2.75,4.90341e+02*1.91456e-02*1.91456e-02/(pt*pt+1.91456e-02*1.91456e-02)+1.60) ;
AliAODEvent *event = dynamic_cast<AliAODEvent*>(InputEvent());
+ if(!event)AliFatal("Can not get ESD event") ;
Double_t mf = event->GetMagneticField(); //Positive for ++ and negative for --
if(mf<0.){ //field --