1 //Create by Christine Nattrass, Rebecca Scott, Irakli Martashvili
2 //University of Tennessee at Knoxville
4 //by default this runs locally
5 //With the argument true this submits jobs to the grid
6 //As written this requires an xml script tag.xml in the ~/et directory on the grid to submit jobs
7 void runHadEt(bool submit = false, bool data = false, bool PbPb = false) {
10 gSystem->Load("libTree.so");
11 gSystem->Load("libGeom.so");
12 gSystem->Load("libVMC.so");
13 gSystem->Load("libXMLIO.so");
15 gSystem->Load("libSTEERBase.so");
16 gSystem->Load("libESD.so");
17 gSystem->Load("libAOD.so");
19 gSystem->Load("libANALYSIS");
20 gSystem->Load("libANALYSISalice");
22 gSystem->AddIncludePath("-I$ALICE_ROOT/include");
23 gROOT->ProcessLine(".L AliAnalysisEtCuts.cxx+g");
24 gROOT->ProcessLine(".L AliAnalysisHadEtCorrections.cxx+g");
25 gROOT->ProcessLine(".L AliAnalysisEtCommon.cxx+g");
26 gROOT->ProcessLine(".L AliAnalysisHadEt.cxx+g");
27 gROOT->ProcessLine(".L AliAnalysisHadEtMonteCarlo.cxx+g");
28 gROOT->ProcessLine(".L AliAnalysisHadEtReconstructed.cxx+g");
29 gROOT->ProcessLine(".L AliAnalysisEtSelectionContainer.cxx+g");
30 gROOT->ProcessLine(".L AliAnalysisEtSelectionHandler.cxx+g");
31 gROOT->ProcessLine(".L AliAnalysisTaskTransverseEnergy.cxx+g");
32 gROOT->ProcessLine(".L AliAnalysisTaskHadEt.cxx+g");
35 char *kTreeName = "esdTree" ;
36 TChain * chain = new TChain(kTreeName,"myESDTree") ;
38 gSystem->Load("libNetx.so") ;
39 gSystem->Load("libgapiUI.so");
40 gSystem->Load("libRAliEn.so");
41 TGrid::Connect("alien://") ;
43 //chain->Add("/data/LHC10h12/999/AliESDs.root");//Hijing Pb+Pb
44 chain->Add("/data/LHC10d15/1821/AliESDs.root");//simulation p+p
45 //chain->Add("/data/LHC10dpass2/10000126403050.70/AliESDs.root");//data
47 // Make the analysis manager
48 AliAnalysisManager *mgr = new AliAnalysisManager("TotEtManager");
51 gROOT->LoadMacro("CreateAlienHandlerHadEtSimPbPb.C");
52 AliAnalysisGrid *alienHandler = CreateAlienHandlerHadEtSimPbPb();
55 gROOT->LoadMacro("CreateAlienHandlerHadEtSim.C");
56 AliAnalysisGrid *alienHandler = CreateAlienHandlerHadEtSim();
58 if (!alienHandler) return;
59 mgr->SetGridHandler(alienHandler);
62 AliVEventHandler* esdH = new AliESDInputHandler;
63 mgr->SetInputEventHandler(esdH);
64 AliMCEventHandler* handler = new AliMCEventHandler;
66 handler->SetReadTR(kFALSE);
67 mgr->SetMCtruthEventHandler(handler);
69 AliAnalysisTaskHadEt *task2 = new AliAnalysisTaskHadEt("TaskHadEt");
71 AliAnalysisDataContainer *cinput1 = mgr->GetCommonInputContainer();
72 AliAnalysisDataContainer *coutput2 = mgr->CreateContainer("out2", TList::Class(), AliAnalysisManager::kOutputContainer,"Et.ESD.new.sim.root");
73 mgr->ConnectInput(task2,0,cinput1);
74 mgr->ConnectOutput(task2,1,coutput2);
76 mgr->SetDebugLevel(0);
78 if (!mgr->InitAnalysis()) return;
81 mgr->StartAnalysis("grid");
84 mgr->StartAnalysis("local",chain);