.so cleanup: removed from gSystem->Load()
[u/mrichter/AliRoot.git] / HLT / JET / macros / tasks / JetAnalysisManagerHLTMC.C
1 void JetAnalysisManagerHLTMC() {
2
3   if ( getenv("FASTJET") ) {
4     gSystem->Load("libCGAL");
5     gSystem->Load("libfastjet");
6     gSystem->Load("libSISConePlugin");
7   }
8
9   gSystem->Load("libTree");
10   gSystem->Load("libPhysics");
11   gSystem->Load("libGeom");
12   gSystem->Load("libVMC");
13   gSystem->Load("libANALYSIS");
14   gSystem->Load("libSTEERBase");
15   gSystem->Load("libAOD");
16   gSystem->Load("libESD");
17   gSystem->Load("libANALYSISalice");
18   gSystem->Load("libJETAN");
19  
20   gSystem->Load("libHLTbase");
21   gSystem->Load("libAliHLTUtil");
22   gSystem->Load("libAliHLTJET");
23
24   Int_t debugLevel = 10;
25   
26   // --------------------------------------------------------------------------------  
27     
28   if (gApplication) gApplication->InitializeGraphics();
29
30   // --------------------------------------------------------------------------------  
31   //
32   // Create the chain
33   //
34   // --------------------------------------------------------------------------------  
35
36   gROOT->LoadMacro("${ALICE_ROOT}/PWG0/CreateESDChain.C");
37   TChain* chain = new TChain("TE");
38
39   chain->Add("/home/jthaeder/jet/data/HEAD_2009-10-26/FastGen/kPythia6Jets125_150_14TeV/JET-ETA=-0.2,0.2_JET-ET=50,1000_R=0.4_200ev/galice.root");
40
41   // --------------------------------------------------------------------------------  
42   //
43   // Create the analysis manager
44   //
45   // --------------------------------------------------------------------------------  
46
47   // Output
48   AliAODHandler* aodHandler = new AliAODHandler();
49   aodHandler->SetOutputFileName("aod.root");
50
51   // MC Truth
52   AliMCEventHandler* mcHandler = new AliMCEventHandler();
53   mcHandler->SetReadTR(kFALSE); // for fastgen
54
55   AliAnalysisManager *mgr  = new AliAnalysisManager("Jet Manager", "Jet Manager");
56   mgr->SetOutputEventHandler (aodHandler);
57   mgr->SetMCtruthEventHandler(mcHandler);
58   mgr->SetDebugLevel(debugLevel);
59
60   // --------------------------------------------------------------------------------
61   //
62   // Set Configfiles
63   //
64   // --------------------------------------------------------------------------------
65
66   // -- HLT FFSC - MC
67 #if FFSC
68   AliAnalysisTaskJets *taskFFSC = new AliAnalysisTaskJets("JetAnalysisHLTMC_FFSC");
69   taskFFSC->SetConfigFile("./tasks/ConfigJetAnalysisHLTMC.C");
70   taskFFSC->SetNonStdBranch("jetsHLTMC_FFSC");
71   taskFFSC->SetDebugLevel(debugLevel);
72   mgr->AddTask(taskFFSC);
73 #endif
74
75   // -- HLT FastJet - Kt
76 #if FASTJET_KT
77   if ( getenv("FASTJET") ) {
78     AliAnalysisTaskJets *taskKt = new AliAnalysisTaskJets("JetAnalysisHLTMC_Kt");
79     taskKt->SetConfigFile("./tasks/ConfigJetAnalysisHLTMCKt.C");
80     taskKt->SetNonStdBranch("jetsHLTMC_Kt");
81     taskKt->SetDebugLevel(debugLevel);
82     mgr->AddTask(taskKt);
83   }
84 #endif
85
86 #if FASTJET_ANTIKT
87   // -- HLT FastJet - AntiKt
88   if ( getenv("FASTJET") ) {
89     AliAnalysisTaskJets *taskAntiKt = new AliAnalysisTaskJets("JetAnalysisHLTMC_AntiKt");
90     taskAntiKt->SetConfigFile("./tasks/ConfigJetAnalysisHLTMCAntiKt.C");
91     taskAntiKt->SetNonStdBranch("jetsHLTMC_AntiKt");
92     taskAntiKt->SetDebugLevel(debugLevel);
93     mgr->AddTask(taskAntiKt);
94   }
95 #endif
96
97   // --------------------------------------------------------------------------------
98   //
99   // Create containers for input/output
100   //
101   // --------------------------------------------------------------------------------  
102
103   AliAnalysisDataContainer *cinput1 = mgr->GetCommonInputContainer();
104   AliAnalysisDataContainer *coutput1 = mgr->GetCommonOutputContainer();
105
106 #if FFSC
107   AliAnalysisDataContainer *coutputFFSC = mgr->CreateContainer("histos", TList::Class(), 
108                                                                AliAnalysisManager::kOutputContainer, "histos_MC_FFSC.root");
109   mgr->ConnectInput  (taskFFSC,     0, cinput1  );
110   mgr->ConnectOutput (taskFFSC,     0, coutput1 );
111   mgr->ConnectOutput (taskFFSC,     1, coutputFFSC );
112 #endif
113   
114 #if FASTJET_KT
115   if ( getenv("FASTJET") ) {
116     AliAnalysisDataContainer *coutputKt = mgr->CreateContainer("histos", TList::Class(), 
117                                                                AliAnalysisManager::kOutputContainer, "histos_MC_KT.root");
118     mgr->ConnectInput  (taskKt,     0, cinput1  );
119     mgr->ConnectOutput (taskKt,     0, coutput1 );
120     mgr->ConnectOutput (taskKt,     1, coutputKt );
121   }
122 #endif
123
124 #if FASTJET_ANTIKT
125   if ( getenv("FASTJET") ) {
126     AliAnalysisDataContainer *coutputAntiKt = mgr->CreateContainer("histos", TList::Class(), 
127                                                                    AliAnalysisManager::kOutputContainer, "histos_MC_ANTIKT.root");
128     mgr->ConnectInput  (taskAntiKt,     0, cinput1  );
129     mgr->ConnectOutput (taskAntiKt,     0, coutput1 );
130     mgr->ConnectOutput (taskAntiKt,     1, coutputAntiKt );
131   }
132 #endif
133
134   //  AliAnalysisTaskKineFilter *kinefilter = new AliAnalysisTaskKineFilter("Kine Filter");
135   //  mgr->AddTask(kinefilter);
136   //  mgr->ConnectInput  (kinefilter,     0, cinput1  );
137   //  mgr->ConnectOutput (kinefilter,     0, coutput1 );
138
139   // --------------------------------------------------------------------------------  
140   //
141   // Run the analysis
142   //    
143   // --------------------------------------------------------------------------------  
144
145   mgr->InitAnalysis();
146   mgr->PrintStatus();
147   mgr->StartAnalysis("local",chain);
148 }