added plotting versus refmult and modified some classes to not stream
[u/mrichter/AliRoot.git] / PWG2 / FLOW / Tools / macros / plotCumulants.C
1 // add comment
2  
3 // Set how many output analysis files in total you want to access:
4 const Int_t nFiles = 1;
5  
6 // Set how many of those output analysis files you want to represent with a mesh (usually used to represent results of simulations):
7 const Int_t nSim = 0;
8  
9 // Set paths of all output analysis files (first the ones to be represented with mesh (simulations), then the ones to be represented with markers (real data))
10 TString files[nFiles] = {"data"};
11  
12 // Set analysis types for all output analysis files (can be "ESD","AOD","MC",""):
13 TString type[nFiles] = {"ESD"};
14  
15 // Set mesh color:
16 Int_t meshColor[nSim] = {};
17
18 // Set marker styles:
19 Int_t markerStyle[nFiles-nSim] = {kStar};
20
21 // Set marker colors:
22 Int_t markerColor[nFiles-nSim] = {kBlack};
23
24 // Set legend entries:
25 TString legendEntry[nFiles] = {"pp example run"};
26  
27 // Set if you whish to plot cumulants versus <reference multiplicity> (by default they are plotted versus # of RPs):
28 Bool_t plotCumulantsVsReferenceMultiplicity = kTRUE;
29
30 // Set flow values whose theoretical contribution to cumulants will be shown on the plots with the straight coloured lines: 
31 Bool_t showTheoreticalLines = kFALSE;
32 const Int_t nFlowValues = 2;
33 Double_t v[nFlowValues] = {0.1,0.05};
34 Int_t lineColor[nFlowValues] = {kRed,kBlue}; 
35
36 // If the statistical error of 6th and 8th order cumulant is huge you may prefer not to show them:
37 Bool_t plotOnly2ndAnd4thOrderCumulant = kTRUE;
38
39 // For comparison sake show also GFC results with dotted line:
40 Bool_t showAlsoGFCResults = kFALSE;
41 Int_t gfcLineStyle = 3;
42
43 // Set method names which calculate cumulants vs multiplicity:
44 const Int_t nMethods = 2;
45 TString method[nMethods] = {"QC","GFC"}; 
46
47 TFile *commonOutputFiles[nFiles] = {NULL}; // common output files "AnalysisResults.root"
48 TList *lists[nFiles][nMethods] = {{NULL}}; // lists cobj<method> holding objects with results for each method
49 TH1D *cumulantsVsM[nFiles][nMethods][4] = {{{NULL}}}; // histograms with results for cumulants vs multiplicity (4 stands for 4 cumulant orders)
50 TLine *lines[nFlowValues][4] = {{NULL}}; // lines denoting theoretical flow contribution to cumulants
51 TProfile *refMultVsNoOfRPs = NULL; // <reference multipicity> versus # of RPs
52
53 // Ranges for plots:
54 Double_t xMin[4]={0.};
55 Double_t xMax[4]={0.};
56 Double_t yMin[4]={0.};
57 Double_t yMax[4]={0.};
58
59 enum libModes {mLocal,mLocalSource};
60
61 void plotCumulants(Int_t analysisMode=mLocal)
62 {
63  // analysisMode: if analysisMode = mLocal -> analyze data on your computer using aliroot
64  //               if analysisMode = mLocalSource -> analyze data on your computer using root + source files 
65   
66  // Load needed libraries:
67  LoadLibrariesPC(analysisMode);
68  
69  // Access all common output files:
70  TString commonOutputFileName = "AnalysisResults.root"; 
71  AccessCommonOutputFiles(commonOutputFileName);
72  
73  // Get from common output files the lists holding histograms for each method:
74  GetLists();
75  
76  // Get histograms with results for cumulants vs multiplicity:
77  GetHistograms();
78  
79  // Determine ranges for plots:
80  DetermineMinMax();
81  
82  // Make lines which indicate theoretical contributions of flow to cumulants:
83  Lines();
84  
85  // Print number of events and average multiplicities for each common output file:
86  Print();  
87  
88  // Make plots:
89  Plot(); 
90  
91  // Global settings which will affect all plots:
92  GlobalSettings();
93  
94 } // end of void plotCumulants(Int_t analysisMode=mLocal) 
95  
96 // =====================================================================================
97
98 void Plot()
99 {
100  // Make all plots.
101  
102  TCanvas *c = NULL;
103  Int_t coMax = 0;
104  if(!plotOnly2ndAnd4thOrderCumulant)
105  {
106   c = new TCanvas("c","cumulants");
107   c->Divide(2,2);
108   coMax = 4; 
109  } else 
110    {
111     c = new TCanvas("c","cumulants",1200,500);
112     c->Divide(2,1);  
113     coMax = 2; 
114    } 
115    
116  TLegend *legend = new TLegend(0.1,0.7,0.33,0.9);
117  legend->SetFillStyle(0);
118
119  TString qcFlag[4] = {"QC{2}","QC{4}","QC{6}","QC{8}"};
120  
121  for(Int_t co=0;co<coMax;co++) // cumulant order
122  {
123   c->cd(co+1);
124   StyleHist(qcFlag[co].Data(),co)->Draw(); 
125   // simulations:
126   for(Int_t s=0;s<nSim;s++) 
127   {
128    TGraph *errorMesh = GetErrorMesh(cumulantsVsM[s][0][co]);
129    if(errorMesh)
130    {
131     errorMesh->SetFillColor(meshColor[s]);
132     errorMesh->Draw("lfsame");
133    }
134    if(co==0){legend->AddEntry(errorMesh,legendEntry[s].Data(),"f");}
135   } // end of if(Int_t s=0;s<nSim;s++) 
136   // data:
137   for(Int_t f=nSim;f<nFiles;f++)
138   {
139    // Theoretical lines:
140    if(showTheoreticalLines)
141    {
142     if(f==nSim) // plot them only once
143     {
144      for(Int_t fv=0;fv<nFlowValues;fv++)
145      { 
146       lines[fv][co]->Draw("same");
147       if(co==0){legend->AddEntry(lines[fv][co],Form("v_{2} = %g",v[fv]),"l");}  
148      } 
149     }
150    } 
151    // QC results:
152    if(cumulantsVsM[f][0][co])
153    {
154     cumulantsVsM[f][0][co]->Draw("e1same"); 
155     cumulantsVsM[f][0][co]->SetMarkerStyle(markerStyle[f-nSim]);  
156     cumulantsVsM[f][0][co]->SetMarkerColor(markerColor[f-nSim]); 
157     if(co==0)
158     {
159      if(showAlsoGFCResults)
160      {
161       legend->AddEntry(cumulantsVsM[f][0][co],Form("%s (QC)",legendEntry[f].Data()),"p");
162      } else
163        {
164         legend->AddEntry(cumulantsVsM[f][0][co],legendEntry[f].Data(),"p");     
165        }
166     }
167    }
168    // GFC results:
169    if(showAlsoGFCResults && cumulantsVsM[f][1][co])
170    {
171     cumulantsVsM[f][1][co]->Draw("lsame");  
172     cumulantsVsM[f][1][co]->SetLineStyle(gfcLineStyle);    
173     if(co==0){legend->AddEntry(cumulantsVsM[f][1][co],Form("%s (GFC)",legendEntry[f].Data()),"l");}
174    }
175   } 
176   // Draw legend:
177   if(co==0){legend->Draw("same");}
178  } // end of for(Int_t co=0;co<4;co++) // cumulant order
179  
180  // Plot also <reference multiplicity> vs # of RPs:
181  if(plotCumulantsVsReferenceMultiplicity)
182  {
183   if(refMultVsNoOfRPs)
184   {
185    TCanvas *cRefMultVsNoOfRPs = new TCanvas("cRefMultVsNoOfRPs","#LTreference multiplicity#GT vs # of RPs");
186    refMultVsNoOfRPs->SetTitle("");
187    refMultVsNoOfRPs->GetXaxis()->SetRangeUser(0,refMultVsNoOfRPs->FindLastBinAbove());
188    refMultVsNoOfRPs->Draw();   
189   } else
190     {
191      cout<<endl;
192      cout<<"WARNING: refMultVsNoOfRPs is NULL in Plot() !!!!"<<endl;
193      cout<<endl;
194     }
195  }   
196  
197 } // end of void Plot()
198  
199 // =====================================================================================
200
201 void Lines()
202 {
203  // Make lines denoting theoretical contribution of flow to cumulants.
204
205  for(Int_t co=0;co<4;co++)
206  {
207   xMin[co] = 0.;
208   //xMax[co] = 59.5;
209  }
210  
211  for(Int_t fv=0;fv<nFlowValues;fv++)
212  {
213   lines[fv][0] = new TLine(xMin[0],pow(v[fv],2),xMax[0]+0.5,pow(v[fv],2));
214   lines[fv][0]->SetLineColor(lineColor[fv]);
215   lines[fv][1] = new TLine(xMin[1],-pow(v[fv],4),xMax[1]+0.5,-pow(v[fv],4));
216   lines[fv][1]->SetLineColor(lineColor[fv]);
217   lines[fv][2] = new TLine(xMin[2],4.*pow(v[fv],6),xMax[2]+0.5,4.*pow(v[fv],6));
218   lines[fv][2]->SetLineColor(lineColor[fv]);
219   lines[fv][3] = new TLine(xMin[3],-33.*pow(v[fv],8),xMax[3]+0.5,-33.*pow(v[fv],8));
220   lines[fv][3]->SetLineColor(lineColor[fv]);
221  }
222
223 } // end of void Lines()
224
225 // =====================================================================================
226
227 void Print()
228 {
229  // Print number of events and average multiplicities for each common output file.
230  
231  cout<<endl;
232  cout<<"Accessed files:"<<endl;
233  cout<<endl;
234  for(Int_t f=0;f<nFiles;f++)
235  {
236   cout<<commonOutputFiles[f]->GetName()<<endl;
237   for(Int_t m=0;m<nMethods;m++)
238   {
239    AliFlowCommonHist *commonHist = dynamic_cast<AliFlowCommonHist*> (lists[f][m]->FindObject(Form("AliFlowCommonHist%s",method[m].Data())));
240    Double_t nEvts = -1.;
241    Double_t AvM = -1.;
242    if(commonHist && commonHist->GetHistMultRP())
243    {
244     nEvts = commonHist->GetHistMultRP()->GetEntries();
245     AvM = commonHist->GetHistMultRP()->GetMean();
246    }
247    if(!(strcmp(method[m].Data(),"QC")))
248    {
249     cout<<Form("%s:",method[m].Data())<<"  <M> = "<<AvM<<", N = "<<nEvts<<endl;
250    }
251    if(!(strcmp(method[m].Data(),"GFC")) && showAlsoGFCResults)
252    {
253     cout<<Form("%s:",method[m].Data())<<"  <M> = "<<AvM<<", N = "<<nEvts<<endl;
254    }
255   }
256   cout<<endl;
257  } // end of for(Int_t f=0;f<nFiles;f++) 
258
259 } // end of void Print()
260
261 // =====================================================================================
262
263 void DetermineMinMax()
264 {
265  // Determine ranges for plots.
266  
267  for(Int_t co=0;co<4;co++)
268  {
269   xMin[co] = 0.; yMin[co] = 44.;
270   xMax[co] = -440000.; yMax[co] = -44.;
271  }
272  
273  Double_t tfc[nFlowValues][4] = {{0.}}; // theoretical flow contribution
274  for(Int_t fv=0;fv<nFlowValues;fv++)
275  {
276   tfc[fv][0] = pow(v[fv],2);
277   tfc[fv][1] = -pow(v[fv],4);
278   tfc[fv][2] = 4.*pow(v[fv],6);
279   tfc[fv][3] = -33.*pow(v[fv],8);
280  }
281   
282  for(Int_t f=0;f<nFiles;f++)
283  {
284   for(Int_t m=0;m<nMethods;m++)
285   { 
286    for(Int_t co=0;co<4;co++)
287    { 
288     if(cumulantsVsM[f][m][co]) 
289     {
290      for(Int_t b=1;b<=cumulantsVsM[f][m][co]->GetXaxis()->GetNbins();b++)
291      {
292       Double_t result = cumulantsVsM[f][m][co]->GetBinContent(b);
293       Double_t error = cumulantsVsM[f][m][co]->GetBinError(b);
294       if(TMath::Abs(result)>1.e-44 && TMath::Abs(error)>1.e-44)
295       {
296        // y-axis:
297        if(yMin[co] > result-error){yMin[co] = result-error;} // min value
298        if(yMax[co] < result+error) {yMax[co] = result+error;} // max value    
299        // x-axis:
300        xMax[co] = b; 
301       }
302      } // end of for(Int_t b=1;b<=cumulantsVsM[f][m][co]->GetXaxis()->GetNbins();b++) 
303      // theoretical contributions:
304      for(Int_t fv=0;fv<nFlowValues;fv++)
305      {
306       //if(yMin[co] > tfc[fv][0]) {yMin[co] = tfc[fv][0];} // min value
307       //if(yMax[co] < tfc[fv][0]) {yMax[co] = tfc[fv][0];} // max value      
308      } // end of for(Int_t fv=0;fv<nFlowValues;fv++)
309     } // end of if(cumulantsVsM[f][m][co])
310    } // end of for(Int_t co=0;co<4;co++)
311   } // end of for(Int_t m=0;m<nMethods;m++)
312  } // end of for(Int_t f=0;f<nFiles;f++)
313  
314 } // end of void DetermineMinMax()
315
316 // =====================================================================================
317
318 void GetHistograms()
319 {
320  // Get histograms with results for cumulants vs multiplicity.
321  
322  if(plotCumulantsVsReferenceMultiplicity) {GetProfileForRefMultVsNoOfRPs();}
323   
324  TString qcFlag[4] = {"QC{2}","QC{4}","QC{6}","QC{8}"};
325  TString gfcFlag[4] = {"GFC{2}","GFC{4}","GFC{6}","GFC{8}"};
326  for(Int_t f=0;f<nFiles;f++)
327  {
328   for(Int_t m=0;m<nMethods;m++)
329   { 
330    TList *temp = NULL;
331    if(!(strcmp(method[m].Data(),"QC")))
332    {
333     temp = dynamic_cast<TList*> (lists[f][m]->FindObject("Integrated Flow"));
334     if(temp) {temp = dynamic_cast<TList*> (temp->FindObject("Results"));}
335     if(temp) 
336     {
337      for(Int_t co=0;co<4;co++)
338      {
339       cumulantsVsM[f][m][co] = dynamic_cast<TH1D*> (temp->FindObject(Form("fIntFlowQcumulantsVsM, %s",qcFlag[co].Data())));
340       if(plotCumulantsVsReferenceMultiplicity)
341       {
342        Map(cumulantsVsM[f][m][co]);
343       }    
344      } 
345     } 
346    } // end of if(!(strcmp(method[m].Data(),"QC")))
347    else if(!(strcmp(method[m].Data(),"GFC")))
348    {
349     temp = dynamic_cast<TList*> (lists[f][m]->FindObject("Reference Flow"));
350     if(temp) {temp = dynamic_cast<TList*> (temp->FindObject("Results"));}
351     if(temp) 
352     {
353      for(Int_t co=0;co<4;co++)
354      {
355       cumulantsVsM[f][m][co] = dynamic_cast<TH1D*> (temp->FindObject(Form("fReferenceFlowCumulantsVsM, %s",gfcFlag[co].Data())));
356       if(plotCumulantsVsReferenceMultiplicity)
357       {
358        Map(cumulantsVsM[f][m][co]);
359       }
360      }
361     } 
362    } // end of else if(!(strcmp(method[m].Data(),"QC")))
363   } // end of  for(Int_t m=0;m<nMethods;m++)
364  } // end of for(Int_t f=0;f<nFiles;f++)
365
366 } // end of void GetHistograms()
367
368 // =====================================================================================
369
370 void GetProfileForRefMultVsNoOfRPs()
371 {
372  // Get profile holding <reference multiplicity> versus # of RPs.
373  
374  AliFlowCommonHist *commonHist = dynamic_cast<AliFlowCommonHist*> (lists[0][0]->FindObject("AliFlowCommonHistQC"));
375  if(commonHist && commonHist->GetRefMultVsNoOfRPs())
376  {
377   refMultVsNoOfRPs = commonHist->GetRefMultVsNoOfRPs();
378  } else
379    {
380     cout<<endl;
381     cout<<"WARNING: commonHist && commonHist->GetRefMultVsNoOfRPs() is NULL in method"<<endl;
382     cout<<"         GetProfileForRefMultVsNoOfRPs() !!!! But was searching only in the file"<<endl;
383     cout<<"         "<<commonOutputFiles[0]->GetName()<<", though."<<endl;
384     cout<<"         Do you have this file available?"<<endl; 
385     cout<<endl;
386    }
387  
388 } // end of void GetProfileForRefMultVsNoOfRPs()
389
390 // =====================================================================================
391
392 void Map(TH1D *hist)
393 {
394  // Map cumulant versus # of RPs into cumulants versus <reference multiplicity>.
395  
396  TH1D *temp = NULL;
397  if(hist)
398  {
399   temp = (TH1D*) hist->Clone();
400   temp->Reset();
401   Int_t nBins = hist->GetNbinsX();
402   for(Int_t b=1;b<=nBins;b++)
403   {
404    temp->SetBinContent((Int_t)refMultVsNoOfRPs->GetBinContent(b),hist->GetBinContent(b));
405    temp->SetBinError((Int_t)refMultVsNoOfRPs->GetBinContent(b),hist->GetBinError(b));
406   }
407  }
408   
409  hist = temp;
410    
411  return;
412   
413 } // end of Map()
414
415 // =====================================================================================
416
417 TGraphErrors* GetGraphErrors(Int_t bin, Int_t nFiles, TH1D** qc)
418 {
419  TGraphErrors *ge = new TGraphErrors(nFiles);
420  for(Int_t f=0;f<nFiles;f++)
421  {
422   ge->SetPoint(f,f+0.5,qc[f]->GetBinContent(bin+1));
423   ge->SetPointError(f,0,qc[f]->GetBinError(bin+1));
424  }
425
426  return ge;
427
428 } // end of TGraphErrors* GetGraphErrors(Int_t bin, Int_t nFiles, TH1D** qc)
429
430 // =====================================================================================
431
432 TH1D* StyleHist(TString yAxisTitle, Int_t co) 
433 {
434  // Style histogram.
435  
436  TH1D *styleHist = new TH1D("","",10000,0,10000); // to be improved (hardwired 10000)
437  // x-axis:
438  styleHist->GetXaxis()->SetRangeUser(xMin[co],xMax[co]);
439  // y-axis:
440  styleHist->GetYaxis()->SetRangeUser(yMin[co],yMax[co]);
441    
442  if(plotCumulantsVsReferenceMultiplicity)
443  {
444   styleHist->GetXaxis()->SetTitle("#LTreference multiplicity#GT");     
445  } else
446    {
447     styleHist->GetXaxis()->SetTitle("# of RPs");         
448    }
449  
450  styleHist->GetYaxis()->SetTitle(yAxisTitle.Data());
451  
452  return styleHist;
453
454 } // end of TH1D* StyleHist(TString yAxisTitle, Int_t co)  
455   
456 // ===========================================================================================
457
458 TGraph* GetErrorMesh(TH1D *hist)
459 {
460  // Error mesh.
461  
462  TGraph* errorMesh = NULL;
463  if(hist)
464  {
465   Int_t nBins = hist->GetNbinsX();
466   Double_t binWidth = hist->GetBinWidth(1); // assuming that all bins have the same width
467   // Counting the non-empty bins: 
468   Int_t nNonEmptyBins=0;
469   for(Int_t i=1;i<nBins+1;i++)
470   {
471    if(!(hist)->GetBinError(i)==0.0))
472    {
473     nNonEmptyBins++;
474    }
475   }   
476   errorMesh = new TGraph(2*nNonEmptyBins+1); 
477   Double_t value=0.,error=0.;
478   Int_t count=1;
479   Double_t xFirst=0.,yUpFirst=0.; // needed to close up the mesh
480   for(Int_t i=1;i<nBins+1;i++)
481   {
482    // Setting up the upper limit of the mesh:
483    value = hist->GetBinContent(i);
484    error = hist->GetBinError(i);   
485    if(!(error==0.0))
486    {    
487     errorMesh->SetPoint(count++,(i-0.5)*binWidth,value+error);
488     if(xFirst==0.)
489     {
490      xFirst=(i-0.5)*binWidth;
491      yUpFirst=value+error;
492     }
493    } 
494   }   
495   for(Int_t i=nBins+1;i<2*nBins+1;i++)
496   {
497    // Setting up the lower limit of the mesh:
498    value = hist->GetBinContent(2*nBins+1-i);
499    error = hist->GetBinError(2*nBins+1-i); 
500    if(!(error==0.0))
501    {      
502     errorMesh->SetPoint(count++,(2*nBins-i+0.5)*binWidth,value-error);
503    }  
504   }
505   // Closing the mesh area:
506   errorMesh->SetPoint(2*nNonEmptyBins+1,xFirst,yUpFirst);   
507  } // end if(hist)
508  
509  errorMesh->SetFillStyle(1001);
510  
511  return errorMesh;
512  
513 } // end of TGraph* GetErrorMesh(TH1D *hist)
514
515 // ===========================================================================================
516
517 void GlobalSettings()
518 {
519  // Settings which will affect all plots.
520  
521  gROOT->SetStyle("Plain"); // default color is white instead of gray
522  gStyle->SetOptStat(0); // remove stat. box from all histos
523  
524 } // end of void GlobalSettings()
525
526 // ===========================================================================================
527
528 void GetLists() 
529 {
530  // Get from common output files the lists holding histograms for each method.
531
532  TString fileName[nFiles][nMethods]; 
533  TDirectoryFile *dirFile[nFiles][nMethods] = {{NULL}}; 
534  TString listName[nFiles][nMethods]; 
535  for(Int_t f=0;f<nFiles;f++)
536  { 
537   for(Int_t i=0;i<nMethods;i++)
538   {
539    // Form a file name for each method:
540    fileName[f][i]+="output";
541    fileName[f][i]+=method[i].Data();
542    fileName[f][i]+="analysis";
543    fileName[f][i]+=type[f].Data();
544    // Access this file:
545    if(commonOutputFiles[f]){dirFile[f][i] = (TDirectoryFile*)commonOutputFiles[f]->FindObjectAny(fileName[f][i].Data());}
546    // Form a list name for each method:
547    listName[f][i]+="cobj";
548    listName[f][i]+=method[i].Data();
549    // Access this lists:
550    if(dirFile[f][i])
551    {
552     dirFile[f][i]->GetObject(listName[f][i].Data(),lists[f][i]); 
553    } else 
554      {
555       cout<<"WARNING: Couldn't find a file "<<fileName[f][i].Data()<<".root in "<<commonOutputFiles[f]->GetName()<<" !!!!"<<endl;exit(0);
556      }
557   } // end of for(Int_t i=0;i<nMethods;i++)   
558  } // end of for(Int_t f=0;f<nFiles;f++)
559
560 } // end of void GetLists() 
561
562 // ===========================================================================================
563
564 void AccessCommonOutputFiles(TString commonOutputFileName)
565 {
566  // Access all output files.
567  for(Int_t f=0;f<nFiles;f++)
568  { 
569   if(!(gSystem->AccessPathName(Form("%s/%s/%s",gSystem->pwd(),files[f].Data(),commonOutputFileName.Data()),kFileExists)))
570   {
571    commonOutputFiles[f] = TFile::Open(Form("%s/%s/%s",gSystem->pwd(),files[f].Data(),commonOutputFileName.Data()),"READ");
572   } else
573     { 
574      cout<<endl;
575      cout<<"WARNING: Couldn't find the file "<<Form("%s/%s/%s",gSystem->pwd(),files[f].Data(),commonOutputFileName.Data())<<" !!!!"<<endl;
576      cout<<endl;
577      exit(0);
578     }
579  } // end of for(Int_t f=0;f<nFiles;f++)
580  
581 } // void AccessCommonOutputFiles(TString commonOutputFileName);
582
583 // ===========================================================================================
584
585 void LoadLibrariesPC(const libModes analysisMode) {
586   
587   //--------------------------------------
588   // Load the needed libraries most of them already loaded by aliroot
589   //--------------------------------------
590   //gSystem->Load("libTree");
591   gSystem->Load("libGeom");
592   gSystem->Load("libVMC");
593   gSystem->Load("libXMLIO");
594   gSystem->Load("libPhysics");
595   
596   //----------------------------------------------------------
597   // >>>>>>>>>>> Local mode <<<<<<<<<<<<<< 
598   //----------------------------------------------------------
599   if (analysisMode==mLocal) {
600     //--------------------------------------------------------
601     // If you want to use already compiled libraries 
602     // in the aliroot distribution
603     //--------------------------------------------------------
604     
605     //==================================================================================  
606     //load needed libraries:
607     gSystem->AddIncludePath("-I$ROOTSYS/include");
608     //gSystem->Load("libTree");
609     
610     // for AliRoot
611     gSystem->AddIncludePath("-I$ALICE_ROOT/include");
612     //gSystem->Load("libANALYSIS");
613     gSystem->Load("libPWG2flowCommon");
614     //cerr<<"libPWG2flowCommon loaded ..."<<endl;
615     
616   }
617   
618   else if (analysisMode==mLocalSource) {
619     
620     // In root inline compile
621    
622    // Constants  
623     gROOT->LoadMacro("AliFlowCommon/AliFlowCommonConstants.cxx+");
624     gROOT->LoadMacro("AliFlowCommon/AliFlowLYZConstants.cxx+");
625     gROOT->LoadMacro("AliFlowCommon/AliFlowCumuConstants.cxx+");
626     
627     // Flow event
628     gROOT->LoadMacro("AliFlowCommon/AliFlowVector.cxx+"); 
629     gROOT->LoadMacro("AliFlowCommon/AliFlowTrackSimple.cxx+");    
630     gROOT->LoadMacro("AliFlowCommon/AliFlowEventSimple.cxx+");
631     
632     // Cuts
633     gROOT->LoadMacro("AliFlowCommon/AliFlowTrackSimpleCuts.cxx+");    
634     
635     // Output histosgrams
636     gROOT->LoadMacro("AliFlowCommon/AliFlowCommonHist.cxx+");
637     gROOT->LoadMacro("AliFlowCommon/AliFlowCommonHistResults.cxx+");
638     gROOT->LoadMacro("AliFlowCommon/AliFlowLYZHist1.cxx+");
639     gROOT->LoadMacro("AliFlowCommon/AliFlowLYZHist2.cxx+");
640     
641     cout << "finished loading macros!" << endl;  
642     
643   }  
644   
645 } // end of void LoadLibraries(const libModes analysisMode)