AliMCEventHandler *mcH = (AliMCEventHandler*) ((AliAnalysisManager::GetAnalysisManager())->GetMCtruthEventHandler());
if (!mcH) {
Printf("ERROR: Could not get MCInputHandler");
+ fSim = kFALSE;
+ }
+ else{
+ fMC = mcH->MCEvent();
+ fSim = kTRUE;
}
- else fMC = mcH->MCEvent();
}
}
extPi->SetParticle("kPiPlus"); //set the particle type
extPi->ApplyEff(); //correct the pt distribution !!HAS TO RUN BEFORE extrapolation!!
extPi->BoltzmannFit(); //fit and extrapolate using Boltzmann-Gibbs Blast wave model
- extPi->TsallisFit(); //fit and extrapolate using Tsallis Blast wave model
+// extPi->TsallisFit(); //fit and extrapolate using Tsallis Blast wave model
TList *extOutListPi = extPi->GetOutputList();
AliAnalysisCentralExtrapolate *extK = new AliAnalysisCentralExtrapolate("extrapolationK");
extK->SetParticle("kKPlus");
extK->ApplyEff();
extK->BoltzmannFit();
- extK->TsallisFit();
+// extK->TsallisFit();
TList *extOutListK = extK->GetOutputList();
AliAnalysisCentralExtrapolate *extP = new AliAnalysisCentralExtrapolate("extrapolationP");
extP->SetParticle("kProton");
extP->ApplyEff();
extP->BoltzmannFit();
- extP->TsallisFit();
+// extP->TsallisFit();
TList *extOutListP = extP->GetOutputList();
Bool_t CheckCuts(Int_t no, TObject *obj) const; //used to check if a track/particle is selected
- void SetSimulation(Bool_t type) {fSim = type;} // set to kTRUE if running on simulated data
+// void SetSimulation(Bool_t type) {fSim = type;} // set to kTRUE if running on simulated data
virtual void ConnectInputData(Option_t *);
virtual void CreateOutputObjects();
// Create and configure the task
AliAnalysisTaskCentral *taskcentral = new AliAnalysisTaskCentral("TaskCentral");
- taskcentral->SetSimulation(kTRUE); //kTRUE if we are running on simulated data
mgr->AddTask(taskcentral);
// Create ONLY the output containers for the data produced by the task.