Updates to Trains. create a job-script to help
[u/mrichter/AliRoot.git] / PWGLF / SPECTRA / IdentifiedHighPt / grid / macro_HighPtDeDx_lhc10b_Data_ESDs.C
CommitLineData
4ebdd20e 1const char *anatype = "ESD";
2
3void macro_HighPtDeDx_lhc10b_Data_ESDs()
4{
5// Analysis using ESD data
6// Automatically generated analysis steering macro executed in grid subjobs
7
8 TStopwatch timer;
9 timer.Start();
10
11// Set temporary merging directory to current one
12 gSystem->Setenv("TMPDIR", gSystem->pwd());
13
14// Set temporary compilation directory to current one
15 gSystem->SetBuildDir(gSystem->pwd(), kTRUE);
16
17// Reset existing include path and add current directory first in the search
18 gSystem->SetIncludePath("-I.");
19// Load analysis framework libraries
20 gSystem->Load("libANALYSIS");
21 gSystem->Load("libOADB");
22 gSystem->Load("libANALYSISalice");
23 gSystem->Load("libCORRFW");
24
25// include path
26 TString intPath = gInterpreter->GetIncludePath();
27 TObjArray *listpaths = intPath.Tokenize(" ");
28 TIter nextpath(listpaths);
29 TObjString *pname;
30 while ((pname=(TObjString*)nextpath())) {
31 TString current = pname->GetName();
32 if (current.Contains("AliRoot") || current.Contains("ALICE_ROOT")) continue;
33 gSystem->AddIncludePath(current);
34 }
35 if (listpaths) delete listpaths;
36 gSystem->AddIncludePath("-I$ROOTSYS/include -I$ALICE_ROOT/include ");
37 gROOT->ProcessLine(".include $ALICE_ROOT/include");
38 printf("Include path: %s\n", gSystem->GetIncludePath());
39
40// Add aditional AliRoot libraries
41
42// analysis source to be compiled at runtime (if any)
43 gROOT->ProcessLine(".L DebugClasses.C+g");
44 gROOT->ProcessLine(".L AliAnalysisTaskHighPtDeDx.cxx+g");
45 gROOT->ProcessLine(".L AliAnalysisTaskHighPtDeDxV0.cxx+g");
46
47// fast xrootd reading enabled
48 printf("!!! You requested FastRead option. Using xrootd flags to reduce timeouts. Note that this may skip some files that could be accessed !!!");
49 gEnv->SetValue("XNet.ConnectTimeout",50);
50 gEnv->SetValue("XNet.RequestTimeout",50);
51 gEnv->SetValue("XNet.MaxRedirectCount",2);
52 gEnv->SetValue("XNet.ReconnectTimeout",50);
53 gEnv->SetValue("XNet.FirstConnectMaxCnt",1);
54
55// connect to AliEn and make the chain
56 if (!TGrid::Connect("alien://")) return;
57// read the analysis manager from file
58 AliAnalysisManager *mgr = AliAnalysisAlien::LoadAnalysisManager("HighPtDeDx_lhc10b_Data_ESDs.root");
59 if (!mgr) return;
60 mgr->PrintStatus();
61 AliLog::SetGlobalLogLevel(AliLog::kError);
62 TChain *chain = CreateChain("wn.xml", anatype);
63
64 mgr->StartAnalysis("localfile", chain);
65 timer.Stop();
66 timer.Print();
67}
68
69//________________________________________________________________________________
70TChain* CreateChain(const char *xmlfile, const char *type="ESD")
71{
72// Create a chain using url's from xml file
73 TString filename;
74 Int_t run = 0;
75 TString treename = type;
76 treename.ToLower();
77 treename += "Tree";
78 printf("***************************************\n");
79 printf(" Getting chain of trees %s\n", treename.Data());
80 printf("***************************************\n");
81 TAlienCollection *coll = TAlienCollection::Open(xmlfile);
82 if (!coll) {
83 ::Error("CreateChain", "Cannot create an AliEn collection from %s", xmlfile);
84 return NULL;
85 }
86 AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
87 TChain *chain = new TChain(treename);
88 coll->Reset();
89 while (coll->Next()) {
90 filename = coll->GetTURL();
91 if (mgr) {
92 Int_t nrun = AliAnalysisManager::GetRunFromAlienPath(filename);
93 if (nrun && nrun != run) {
94 printf("### Run number detected from chain: %d\n", nrun);
95 mgr->SetRunFromPath(nrun);
96 run = nrun;
97 }
98 }
99 chain->Add(filename);
100 }
101 if (!chain->GetNtrees()) {
102 ::Error("CreateChain", "No tree found from collection %s", xmlfile);
103 return NULL;
104 }
105 return chain;
106}
107