1 #if !defined(__CINT__) || defined(__MAKECINT__)
4 #include "AliAnalysisManager.h"
5 #include "AliESDInputHandler.h"
6 #include "AliAODHandler.h"
7 #include "AliAnalysisTaskESDfilter.h"
8 #include "AliAnalysisDataContainer.h"
11 void CreateAODfromESD(const char *inFileName = "AliESDs.root",
12 const char *outFileName = "AliAOD.root") {
14 gSystem->Load("libTree");
15 gSystem->Load("libGeom");
16 gSystem->Load("libPhysics");
17 gSystem->Load("libVMC");
18 gSystem->Load("libSTEERBase");
19 gSystem->Load("libESD");
20 gSystem->Load("libAOD");
22 gSystem->Load("libANALYSIS");
23 gSystem->Load("libANALYSISalice");
24 gSystem->Load("libPWG3muon");
26 TChain *chain = new TChain("esdTree");
27 // Steering input chain
28 chain->Add(inFileName);
29 AliAnalysisManager *mgr = new AliAnalysisManager("ESD to AOD", "Analysis Manager");
32 AliESDInputHandler* inpHandler = new AliESDInputHandler();
33 mgr->SetInputEventHandler (inpHandler);
36 AliAODHandler* aodHandler = new AliAODHandler();
37 aodHandler->SetOutputFileName(outFileName);
38 mgr->SetOutputEventHandler(aodHandler);
42 AliAnalysisTaskESDfilter *filter = new AliAnalysisTaskESDfilter("Filter");
45 AliAnalysisTaskESDMuonFilter *esdmuonfilter = new AliAnalysisTaskESDMuonFilter("ESD Muon Filter");
46 mgr->AddTask(esdmuonfilter);
48 AliESDtrackCuts* esdTrackCutsL = new AliESDtrackCuts("AliESDtrackCuts", "Standard");
49 esdTrackCutsL->SetMinNClustersTPC(50);
50 esdTrackCutsL->SetMaxChi2PerClusterTPC(3.5);
51 esdTrackCutsL->SetMaxCovDiagonalElements(2,2,0.5,0.5,2);
52 esdTrackCutsL->SetRequireTPCRefit(kTRUE);
53 esdTrackCutsL->SetMinNsigmaToVertex(3);
54 esdTrackCutsL->SetDCAToVertexXY(3.);
55 esdTrackCutsL->SetRequireSigmaToVertex(kTRUE);
56 esdTrackCutsL->SetAcceptKingDaughters(kFALSE);
58 AliAnalysisFilter* trackFilter = new AliAnalysisFilter("trackFilter");
59 trackFilter->AddCuts(esdTrackCutsL);
61 filter->SetTrackFilter(trackFilter);
64 AliAnalysisDataContainer *cinput1 = mgr->CreateContainer("cchain", TChain::Class(),
65 AliAnalysisManager::kInputContainer);
67 AliAnalysisDataContainer *coutput1 = mgr->CreateContainer("tree", TTree::Class(),
68 AliAnalysisManager::kOutputContainer,
72 mgr->ConnectInput (filter, 0, cinput1 );
73 mgr->ConnectOutput(filter, 0, coutput1);
75 mgr->ConnectInput (esdmuonfilter, 0, cinput1 );
76 mgr->ConnectOutput(esdmuonfilter, 0, coutput1);
83 mgr->StartAnalysis("local", chain);