void runLocal(const char *chainlistfile, int dataFromAlien=0) {
TStopwatch timer;
timer.Start();
printf("*** Connect to AliEn ***\n");
if (dataFromAlien)
TGrid::Connect("alien://");
gSystem->Load("libProofPlayer");
//____________________________________________________//
//_____________Setting up STEERBase.par_______________//
//____________________________________________________//
setupPar("STEERBase");
gSystem->Load("libSTEERBase");
//____________________________________________________//
//_____________Setting up ESD.par_____________________//
//____________________________________________________//
setupPar("ESD");
gSystem->Load("libVMC");
gSystem->Load("libESD");
//____________________________________________________//
//_____________Setting up AOD.par_____________________//
//____________________________________________________//
setupPar("AOD");
gSystem->Load("libAOD");
//_________________________________________________________//
//_____________Setting up ANALYSIS.par_____________________//
//_________________________________________________________//
setupPar("ANALYSIS");
gSystem->Load("libANALYSIS");
//_________________________________________________________//
//_____________Setting up ANALYSISalice.par________________//
//_________________________________________________________//
setupPar("ANALYSISalice");
gSystem->Load("libANALYSISalice");
//____________________________________________________//
//_____________Setting up PWG2AOD.par_________________//
//____________________________________________________//
setupPar("PWG2AOD");
gSystem->Load("libPWG2AOD");
//____________________________________________________//
//_____________Setting up PWG2femtoscopy.par__________//
//____________________________________________________//
setupPar("PWG2femtoscopy");
gSystem->Load("libPWG2femtoscopy");
//____________________________________________________//
//_____________Setting up PWG2femtoscopyUser.par______//
//____________________________________________________//
setupPar("PWG2femtoscopyUser");
gSystem->Load("libPWG2femtoscopyUser");
//ANALYSIS PART
gSystem->SetIncludePath("-I$ROOTSYS/include -I\"/usr/local/CERN/root/include\" -I./PWG2femtoscopy/FEMTOSCOPY/AliFemto -I./PWG2femtoscopyUser/FEMTOSCOPY/AliFemtoUser -I./ESD -I./AOD -I./ANALYSIS -I./PWG2AOD/AOD");
gROOT->LoadMacro("ConfigFemtoAnalysis.C++");
//____________________________________________//
//Usage of event tags
AliTagAnalysis *analysis = new AliTagAnalysis();
TChain *chain = 0x0;
// chain = analysis->GetChainFromCollection(collectionfile,"esdTree");
if (dataFromAlien) {
AliTagAnalysis *analysis = new AliTagAnalysis();
chain = analysis->GetChainFromCollection(chainlistfile,"esdTree");
}
else {
gROOT->LoadMacro("CreateESDChain.C");
chain = CreateESDChain(chainlistfile,500);
}
//____________________________________________//
// Make the analysis manager
AliAnalysisManager *mgr = new AliAnalysisManager("TestManager");
AliESDInputHandler* esdH = new AliESDInputHandler;
AliMCEventHandler *mcH = new AliMCEventHandler;
esdH->SetInactiveBranches("FMD CaloCluster");
mgr->SetInputEventHandler(esdH);
mgr->SetMCtruthEventHandler(mcH);
//____________________________________________//
// 1st Pt task
AliAnalysisTaskFemto *task1 = new AliAnalysisTaskFemto("TaskFemto");
mgr->AddTask(task1);
// Create containers for input/output
// AliAnalysisDataContainer *cinput1 = mgr->GetCommonInputContainer();
AliAnalysisDataContainer *cinput1 = mgr->CreateContainer("input0",
TTree::Class(), AliAnalysisManager::kInputContainer);
AliAnalysisDataContainer *coutput1 = mgr->CreateContainer("clist1", TList::Class(),AliAnalysisManager::kOutputContainer,"Femto.ESD.root");
//____________________________________________//
cinput1->SetData(chain);
mgr->ConnectInput(task1,0,cinput1);
mgr->ConnectOutput(task1,0,coutput1);
if (!mgr->InitAnalysis()) return;
mgr->PrintStatus();
mgr->StartAnalysis("local",chain);
timer.Stop();
timer.Print();
}
void runLocal() {
TStopwatch timer;
timer.Start();
gSystem->Load("libTree.so");
//____________________________________________________//
//_____________Setting up STEERBase.par_______________//
//____________________________________________________//
setupPar("STEERBase");
gSystem->Load("libSTEERBase.so");
//____________________________________________________//
//_____________Setting up ESD.par_____________________//
//____________________________________________________//
setupPar("ESD");
gSystem->Load("libVMC.so");
gSystem->Load("libESD.so");
//____________________________________________________//
//_____________Setting up AOD.par_____________________//
//____________________________________________________//
setupPar("AOD");
gSystem->Load("libAOD.so");
//_____________________________________________________________//
//_____________Setting up ANALYSIS_NEW.par_____________________//
//_____________________________________________________________//
setupPar("ANALYSIS");
gSystem->Load("libANALYSIS");
gSystem->Load("libANALYSISalice");
// gROOT->LoadMacro("../tasks/AliAnalysisTaskTotEt.cxx+");
gSystem->AddIncludePath("-I$ALICE_ROOT/include");
gROOT->LoadMacro("AliAnalysisTaskTotEt.cxx+g");
//____________________________________________//
AliTagAnalysis *TagAna = new AliTagAnalysis("ESD");
TagAna->ChainLocalTags("../Tags");
AliRunTagCuts *runCuts = new AliRunTagCuts();
AliLHCTagCuts *lhcCuts = new AliLHCTagCuts();
AliDetectorTagCuts *detCuts = new AliDetectorTagCuts();
AliEventTagCuts *evCuts = new AliEventTagCuts();
// evCuts->SetMultiplicityRange(11,12);
evCuts->SetNPHOSClustersRange(0,100);
TChain* chain = 0x0;
chain = TagAna->QueryTags(runCuts,lhcCuts,detCuts,evCuts);
// chain->SetBranchStatus("*Calo*",0);
// Make the analysis manager
AliAnalysisManager *mgr = new AliAnalysisManager("TotEtManager");
AliVEventHandler* esdH = new AliESDInputHandler;
mgr->SetInputEventHandler(esdH);
AliAnalysisTaskTotEt *task1 = new AliAnalysisTaskTotEt("TaskTotEt");
mgr->AddTask(task1);
// Create containers for input/output
AliAnalysisDataContainer *cinput1 = mgr->CreateContainer("cchain1", TChain::Class(),AliAnalysisManager::kInputContainer);
AliAnalysisDataContainer *coutput1 = mgr->CreateContainer("chist1", TH1::Class(), AliAnalysisManager::kOutputContainer,"Et.ESD.root");
AliAnalysisDataContainer *coutput2 = mgr->CreateContainer("chist2", TH1::Class(), AliAnalysisManager::kOutputContainer,"Et.ESD.root");
AliAnalysisDataContainer *coutput3 = mgr->CreateContainer("chist3", TH1::Class(), AliAnalysisManager::kOutputContainer,"Et.ESD.root");
AliAnalysisDataContainer *coutput4 = mgr->CreateContainer("chist4", TH1::Class(), AliAnalysisManager::kOutputContainer,"Et.ESD.root");
AliAnalysisDataContainer *coutput5 = mgr->CreateContainer("cnt", TNtuple::Class(), AliAnalysisManager::kOutputContainer,"Et.ESD.root");
// AliAnalysisDataContainer *coutput4 = mgr->CreateContainer("carray1", TClonesArray::Class(), AliAnalysisManager::kOutputContainer,"Pt.ESD.root");
//____________________________________________//
mgr->ConnectInput(task1,0,cinput1);
mgr->ConnectOutput(task1,1,coutput1);
mgr->ConnectOutput(task1,2,coutput2);
mgr->ConnectOutput(task1,3,coutput3);
mgr->ConnectOutput(task1,4,coutput4);
mgr->ConnectOutput(task1,5,coutput5);
// mgr->ConnectOutput(task1,3,coutput4);
mgr->SetDebugLevel(2);
if (!mgr->InitAnalysis()) return;
mgr->PrintStatus();
mgr->StartAnalysis("local",chain);
timer.Stop();
timer.Print();
}
