.so cleanup: removed from gSystem->Load()
[u/mrichter/AliRoot.git] / PWGGA / PHOSTasks / PHOS_pp_pi0 / macros / Pi0Spectrum_LHC11a.C
1 void Pi0Spectrum(const char* dataset="")
2 {
3     
4   gSystem->Load("libTree");
5   gSystem->Load("libGeom");
6   gSystem->Load("libVMC");
7   gSystem->Load("libPhysics");
8   
9   //load analysis framework
10   gSystem->Load("libANALYSIS");
11   gSystem->Load("libANALYSISalice"); //AliAnalysisTaskSE
12   
13   gSystem->AddIncludePath("-I$ALICE_ROOT/include -I$ALICE_ROOT/PHOS");
14
15   // A task can be compiled dynamically with AClic
16   gROOT->LoadMacro("AliCaloPhoton.cxx+g");
17   gROOT->LoadMacro("AliAnalysisTaskPi0.cxx+g");
18   
19   // Connect to alien
20   TString token = gSystem->Getenv("GRID_TOKEN") ;
21   if ( token == "OK" ) 
22     TGrid::Connect("alien://");
23   else 
24     AliInfo("You are not connected to the GRID") ; 
25
26   cout << "Pi0Analysis: processing collection " << dataset << endl;
27   
28   // Create the chain
29   TChain* chain = new TChain("esdTree");
30   TGridCollection * collection = dynamic_cast<TGridCollection*>(TAlienCollection::Open(dataset));
31   
32   TAlienResult* result = collection->GetGridResult("",0 ,0);
33   TList* rawFileList = result->GetFileInfoList();
34   
35   for (Int_t counter=0 ; counter < rawFileList->GetEntries() ; counter++) {
36     TFileInfo * fi =  static_cast<TFileInfo*>(rawFileList->At(counter)) ; 
37     const char * rawFile = fi->GetCurrentUrl()->GetUrl() ;  
38     printf("Processing %s\n", rawFile) ;
39     chain->Add(rawFile);
40     printf("Chain: %d entries.\n",chain->GetEntries()); 
41   }
42   
43   // Make the analysis manager
44   AliAnalysisManager *mgr = new AliAnalysisManager("Pi0Spectrum");
45   
46   // ESD input handler
47   AliESDInputHandler* esdH = new AliESDInputHandler();
48   esdH->SetReadFriends(kFALSE);
49   mgr->SetInputEventHandler(esdH);
50   
51   // Debug level
52   mgr->SetDebugLevel(0);
53
54   // Add physics selection
55   gROOT->LoadMacro("$ALICE_ROOT/ANALYSIS/macros/AddTaskPhysicsSelection.C");
56   AliPhysicsSelectionTask* physSelTask = AddTaskPhysicsSelection();
57
58   // Add my task
59   AliAnalysisTaskPi0 *task1 = new AliAnalysisTaskPi0("Pi0Spectrum");
60
61   TFile *fBadMap = TFile::Open("BadMap_LHC11a_pp2760_20130902.root");
62   if(fBadMap->IsOpen()){
63     printf("\n\n...Adding PHOS bad channel map \n") ;
64     gROOT->cd();
65     char key[55] ;
66     for(Int_t mod=1;mod<4; mod++){
67       sprintf(key,"PHOS_BadMap_mod%d",mod) ;
68       TH2I * h = (TH2I*)fBadMap->Get(key) ;
69       if(h)
70         task1->SetPHOSBadMap(mod,h) ;
71     }
72     fBadMap->Close() ;
73   }
74   // Set BC gap for LHC10e in seconds
75   task1->SetBCgap(525.e-09);
76
77   // Set abs.recalibration for LHC11a
78   task1->SetRecalib(1, 0.9942);
79   task1->SetRecalib(2, 0.9822);
80   task1->SetRecalib(3, 1.0072);
81
82   task1->SelectCollisionCandidates(AliVEvent::kMB);
83   mgr->AddTask(task1);
84   
85   // Create containers for input/output
86   AliAnalysisDataContainer *cinput   = mgr->GetCommonInputContainer(); 
87   AliAnalysisDataContainer *coutput1 = mgr->CreateContainer("histESD",TList::Class(),AliAnalysisManager::kOutputContainer,"histos.root");
88   
89   // Connect input/output
90   mgr->ConnectInput(task1 , 0, cinput);
91   mgr->ConnectOutput(task1, 1, coutput1);
92   
93   if (mgr->InitAnalysis()) {
94     mgr->PrintStatus();
95     mgr->StartAnalysis("local", chain);
96   }
97   
98 }