2 #include <TStopwatch.h>
6 #include <TTimeStamp.h>
8 #include <Base/AliFlowCommonConstants.h>
9 #include <Base/AliFlowLYZConstants.h>
10 #include <Base/AliFlowCumuConstants.h>
11 #include <Base/AliFlowVector.h>
12 #include <Base/AliFlowTrackSimple.h>
13 #include <Base/AliFlowEvent.h>
14 #include <Base/AliFlowEventSimple.h>
15 #include <Base/AliFlowTrackSimpleCuts.h>
16 #include <Base/AliFlowCommonHist.h>
17 #include <Base/AliFlowCommonHistResults.h>
18 #include <Base/AliFlowLYZHist1.h>
19 #include <Base/AliFlowLYZHist2.h>
20 #include <Base/AliCumulantsFunctions.h>
21 #include <Base/AliFlowLYZEventPlane.h>
22 #include <Base/AliFlowAnalysisWithMCEventPlane.h>
23 #include <Base/AliFlowAnalysisWithScalarProduct.h>
24 #include <Base/AliFlowAnalysisWithLYZEventPlane.h>
25 #include <Base/AliFlowAnalysisWithLeeYangZeros.h>
26 #include <Base/AliFlowAnalysisWithCumulants.h>
27 #include <Base/AliFlowAnalysisWithQCumulants.h>
28 #include <Base/AliFlowAnalysisWithFittingQDistribution.h>
29 #include <Base/AliFlowAnalysisWithMixedHarmonics.h>
30 #include <Base/AliFlowAnalysisWithNestedLoops.h>
32 //--------------------------------------------------------------------------------------
33 // Run flow analysis on local data with custom FlowEvent maker
35 //flow analysis method can be: (set to kTRUE or kFALSE)
37 Bool_t LYZ1SUM = kTRUE;
38 Bool_t LYZ1PROD = kTRUE;
39 Bool_t LYZ2SUM = kFALSE;
40 Bool_t LYZ2PROD = kFALSE;
41 Bool_t LYZEP = kFALSE;
47 Bool_t MCEP = kFALSE; //does not work yet 24/12/08
48 //--------------------------------------------------------------------------------------
51 // Use weights for Q vector
52 Bool_t usePhiWeights = kFALSE; //Phi (correction for non-uniform azimuthal acceptance)
53 Bool_t usePtWeights = kFALSE; //v'(pt) (differential flow in pt)
54 Bool_t useEtaWeights = kFALSE; //v'(eta) (differential flow in eta)
56 //--------------------------------------------------------------------------------------
58 //integrated selection
59 Double_t ptMaxInt = 10.;
60 Double_t ptMinInt = 0.;
61 Double_t etaMaxInt = 1.;
62 Double_t etaMinInt = -1.;
63 Double_t phiMaxInt = 7.5;
64 Double_t phiMinInt = 0.;
67 //differential selection
68 Double_t ptMaxDiff = 10.;
69 Double_t ptMinDiff = 0.;
70 Double_t etaMaxDiff = 1.;
71 Double_t etaMinDiff = -1.;
72 Double_t phiMaxDiff = 7.5;
73 Double_t phiMinDiff = 0.;
76 //--------------------------------------------------------------------------------------
77 // FLOW SETTINGS (R.Rietkerk)
78 Int_t nLoops=1; // Number of times to use the same particle (nonflow).
79 Double_t xEllipticFlowValue=0.1;// Add Elliptic Flow. Must be in range [0,1].
80 Int_t nMultiplicityOfEvent=500; // Set Average Multiplicity.
81 Double_t xSigmaFlow=0.00; // Add Elliptic Flow. Must be in range [0,1].
82 Int_t nSigmaMult=50; // Set Average Multiplicity.
83 //--------------------------------------------------------------------------------------
86 enum anaModes {mLocal,mLocalSource,mLocalPAR,};
87 //mLocal: Analyze data on your computer using aliroot
88 //mLocalPAR: Analyze data on your computer using root + PAR files
89 //mLocalSource: Analyze data on your computer using root + source files
91 void LoadLibraries(const anaModes mode);
95 int runFlowAnalysis(const anaModes mode=mLocal, Int_t aRuns = 100, const char*
96 dir="/data/alice1/kolk/KineOnly3/")
97 // dir="/Users/snelling/alice_data/KineOnly3/")
98 // dir="/Users/snelling/alice_data/stoomboot/5b/")
103 if (LYZ1SUM && LYZ2SUM) {cout<<"WARNING: you cannot run LYZ1 and LYZ2 at the same time! LYZ2 needs the output from LYZ1."<<endl; exit(1); }
104 if (LYZ1PROD && LYZ2PROD) {cout<<"WARNING: you cannot run LYZ1 and LYZ2 at the same time! LYZ2 needs the output from LYZ1."<<endl; exit(1); }
105 if (LYZ2SUM && LYZEP) {cout<<"WARNING: you cannot run LYZ2 and LYZEP at the same time! LYZEP needs the output from LYZ2."<<endl; exit(1); }
106 if (LYZ1SUM && LYZEP) {cout<<"WARNING: you cannot run LYZ1 and LYZEP at the same time! LYZEP needs the output from LYZ2."<<endl; exit(1); }
110 cout<<" ---- BEGIN ANALYSIS ---- "<<endl;
115 TRandom3 random3Temp; //init for manual settings (R.Rietkerk)
117 Int_t sseed = dt.GetNanoSec()/1000;
118 random3Temp.SetSeed(sseed);
120 if (mode == mLocal || mode == mLocalPAR) {
122 else if (mode == mLocalSource) {
125 cout << "No supported running mode selected!" << endl;
130 //------------------------------------------------------------------------
132 AliFlowTrackSimpleCuts* cutsInt = new AliFlowTrackSimpleCuts();
133 cutsInt->SetPtMax(ptMaxInt);
134 cutsInt->SetPtMin(ptMinInt);
135 cutsInt->SetEtaMax(etaMaxInt);
136 cutsInt->SetEtaMin(etaMinInt);
137 cutsInt->SetPhiMax(phiMaxInt);
138 cutsInt->SetPhiMin(phiMinInt);
139 cutsInt->SetPID(PIDInt);
141 AliFlowTrackSimpleCuts* cutsDiff = new AliFlowTrackSimpleCuts();
142 cutsDiff->SetPtMax(ptMaxDiff);
143 cutsDiff->SetPtMin(ptMinDiff);
144 cutsDiff->SetEtaMax(etaMaxDiff);
145 cutsDiff->SetEtaMin(etaMinDiff);
146 cutsDiff->SetPhiMax(phiMaxDiff);
147 cutsDiff->SetPhiMin(phiMinDiff);
148 cutsDiff->SetPID(PIDDiff);
150 //if the weights are used:
151 TFile *fileWithWeights = NULL;
152 TList *listWithWeights = NULL;
154 if(usePhiWeights||usePtWeights||useEtaWeights) {
155 fileWithWeights = TFile::Open("weights.root","READ");
156 if(fileWithWeights) {
157 listWithWeights = (TList*)fileWithWeights->Get("weights");
160 {cout << " WARNING: the file <weights.root> with weights from the previous run was not found."<<endl;
166 AliFlowAnalysisWithQCumulants *qc = NULL;
167 AliFlowAnalysisWithCumulants *gfc = NULL;
168 AliFlowAnalysisWithFittingQDistribution *fqd = NULL;
169 AliFlowAnalysisWithLeeYangZeros *lyz1sum = NULL;
170 AliFlowAnalysisWithLeeYangZeros *lyz1prod = NULL;
171 AliFlowAnalysisWithLeeYangZeros *lyz2sum = NULL;
172 AliFlowAnalysisWithLeeYangZeros *lyz2prod = NULL;
173 AliFlowAnalysisWithLYZEventPlane *lyzep = NULL;
174 AliFlowAnalysisWithScalarProduct *sp = NULL;
175 AliFlowAnalysisWithMCEventPlane *mcep = NULL;
176 AliFlowAnalysisWithMixedHarmonics *mh = NULL;
177 AliFlowAnalysisWithNestedLoops *nl = NULL;
179 //MCEP = monte carlo event plane
181 AliFlowAnalysisWithMCEventPlane *mcep = new AliFlowAnalysisWithMCEventPlane();
187 AliFlowAnalysisWithQCumulants* qc = new AliFlowAnalysisWithQCumulants();
188 if(listWithWeights) qc->SetWeightsList(listWithWeights);
189 if(usePhiWeights) qc->SetUsePhiWeights(usePhiWeights);
190 if(usePtWeights) qc->SetUsePtWeights(usePtWeights);
191 if(useEtaWeights) qc->SetUseEtaWeights(useEtaWeights);
195 //GFC = Generating Function Cumulants
197 AliFlowAnalysisWithCumulants* gfc = new AliFlowAnalysisWithCumulants();
198 if(listWithWeights) gfc->SetWeightsList(listWithWeights);
199 if(usePhiWeights) gfc->SetUsePhiWeights(usePhiWeights);
200 if(usePtWeights) gfc->SetUsePtWeights(usePtWeights);
201 if(useEtaWeights) gfc->SetUseEtaWeights(useEtaWeights);
205 //FQD = Fitting q-distribution
207 AliFlowAnalysisWithFittingQDistribution* fqd = new AliFlowAnalysisWithFittingQDistribution();
208 if(listWithWeights) fqd->SetWeightsList(listWithWeights);
209 if(usePhiWeights) fqd->SetUsePhiWeights(usePhiWeights);
213 //SP = Scalar Product
215 AliFlowAnalysisWithScalarProduct* sp = new AliFlowAnalysisWithScalarProduct();
216 if(usePhiWeights) sp->SetUsePhiWeights(usePhiWeights);
220 //LYZ1 = Lee-Yang Zeroes first run
222 AliFlowAnalysisWithLeeYangZeros* lyz1sum = new AliFlowAnalysisWithLeeYangZeros();
223 lyz1sum->SetFirstRun(kTRUE);
224 lyz1sum->SetUseSum(kTRUE);
228 AliFlowAnalysisWithLeeYangZeros* lyz1prod = new AliFlowAnalysisWithLeeYangZeros();
229 lyz1prod->SetFirstRun(kTRUE);
230 lyz1prod->SetUseSum(kFALSE);
233 //LYZ2 = Lee-Yang Zeroes second run
235 AliFlowAnalysisWithLeeYangZeros* lyz2sum = new AliFlowAnalysisWithLeeYangZeros();
236 // read the input file from the first run
237 TString inputFileNameLYZ2SUM = "outputLYZ1SUManalysis.root" ;
238 TFile* inputFileLYZ2SUM = new TFile(inputFileNameLYZ2SUM.Data(),"READ");
239 if(!inputFileLYZ2SUM || inputFileLYZ2SUM->IsZombie()) {
240 cerr << " ERROR: To run LYZ2SUM you need the output file from LYZ1SUM. This file is not there! Please run LYZ1SUM first." << endl ;
244 TList* inputListLYZ2SUM = (TList*)inputFileLYZ2SUM->Get("cobjLYZ1SUM");
245 if (!inputListLYZ2SUM) {cout<<"SUM Input list is NULL pointer!"<<endl; break;}
247 cout<<"LYZ2SUM input file/list read..."<<endl;
248 lyz2sum->SetFirstRunList(inputListLYZ2SUM);
249 lyz2sum->SetFirstRun(kFALSE);
250 lyz2sum->SetUseSum(kTRUE);
256 AliFlowAnalysisWithLeeYangZeros* lyz2prod = new AliFlowAnalysisWithLeeYangZeros();
257 // read the input file from the first run
258 TString inputFileNameLYZ2PROD = "outputLYZ1PRODanalysis.root" ;
259 TFile* inputFileLYZ2PROD = new TFile(inputFileNameLYZ2PROD.Data(),"READ");
260 if(!inputFileLYZ2PROD || inputFileLYZ2PROD->IsZombie()) {
261 cerr << " ERROR: To run LYZ2PROD you need the output file from LYZ1PROD. This file is not there! Please run LYZ1PROD first." << endl ;
265 TList* inputListLYZ2PROD = (TList*)inputFileLYZ2PROD->Get("cobjLYZ1PROD");
266 if (!inputListLYZ2PROD) {cout<<"PROD Input list is NULL pointer!"<<endl; break;}
268 cout<<"LYZ2PROD input file/list read..."<<endl;
269 lyz2prod->SetFirstRunList(inputListLYZ2PROD);
270 lyz2prod->SetFirstRun(kFALSE);
271 lyz2prod->SetUseSum(kTRUE);
276 //LYZEP = Lee-Yang Zeroes event plane
278 AliFlowLYZEventPlane* ep = new AliFlowLYZEventPlane() ;
279 AliFlowAnalysisWithLYZEventPlane* lyzep = new AliFlowAnalysisWithLYZEventPlane();
280 // read the input file from the second lyz run
281 TString inputFileNameLYZEP = "outputLYZ2SUManalysis.root" ;
282 TFile* inputFileLYZEP = new TFile(inputFileNameLYZEP.Data(),"READ");
283 if(!inputFileLYZEP || inputFileLYZEP->IsZombie()) {
284 cerr << " ERROR: To run LYZEP you need the output file from LYZ2SUM. This file is not there! Please run LYZ2SUM first." << endl ;
288 TList* inputListLYZEP = (TList*)inputFileLYZEP->Get("cobjLYZ2SUM");
289 if (!inputListLYZEP) {cout<<"Input list is NULL pointer!"<<endl; break;}
291 cout<<"LYZEP input file/list read..."<<endl;
292 ep ->SetSecondRunList(inputListLYZEP);
293 lyzep->SetSecondRunList(inputListLYZEP);
299 // MH = Mixed Harmonics:
301 AliFlowAnalysisWithMixedHarmonics* mh = new AliFlowAnalysisWithMixedHarmonics();
302 if(listWithWeights) mh->SetWeightsList(listWithWeights);
303 //if(usePhiWeights) mh->SetUsePhiWeights(usePhiWeights); // to be improved (enabled)
304 //if(usePtWeights) mh->SetUsePtWeights(usePtWeights); // to be improved (enabled)
305 //if(useEtaWeights) mh->SetUseEtaWeights(useEtaWeights); // to be improved (enabled)
308 // NL = Nested Loops:
310 AliFlowAnalysisWithNestedLoops* nl = new AliFlowAnalysisWithNestedLoops();
311 if(listWithWeights) nl->SetWeightsList(listWithWeights);
312 //if(usePhiWeights) nl->SetUsePhiWeights(usePhiWeights); // to be improved (enabled)
313 //if(usePtWeights) nl->SetUsePtWeights(usePtWeights); // to be improved (enabled)
314 //if(useEtaWeights) nl->SetUseEtaWeights(useEtaWeights); // to be improved (enabled)
318 //------------------------------------------------------------------------
321 //standard code to read files in directory
323 TString execDir(gSystem->pwd());
324 TString targetDir(dir);
325 TSystemDirectory* baseDir = new TSystemDirectory(".", dir);
326 TList* dirList = baseDir->GetListOfFiles();
328 cout << endl << "No input files in: " << targetDir.Data() << endl;
331 Int_t nDirs = dirList->GetEntries();
333 cout<<"Int_t nDirs = "<<nDirs<<endl;
334 gSystem->cd(execDir);
336 for(Int_t iDir=0;iDir<nDirs;++iDir)
338 TSystemFile* presentDir = (TSystemFile*)dirList->At(iDir);
339 if(!presentDir || !presentDir->IsDirectory() ||
340 strcmp(presentDir->GetName(), ".") == 0 ||
341 strcmp(presentDir->GetName(), "..") == 0)
344 cout << "Directory (" << iDir << "): " << presentDir->GetName() <<
345 " - Skipping ... " << endl;
349 if(offset > 0) { --offset ; continue ; }
350 if((aRuns > 0) && (fCount >= aRuns)) { break ; }
352 TString presentDirName(dir); // aDataDir
353 presentDirName += presentDir->GetName();
354 presentDirName += "/";
355 //cerr<<" presentDirName = "<<presentDirName<<endl;
357 TString fileName = presentDirName;
358 fileName += "galice.root";
359 Long_t *id, *size, *flags, *modtime;
360 if(gSystem->GetPathInfo(fileName.Data(),id,size,flags,modtime))
362 cout << " File : " << fileName << " does NOT exist ! - Skipping ... "
366 // cout << endl ; cout << "Directory (" << iDir << "): " << presentDirName << " ... " << endl;
368 //loop (simulations in the present dir)
369 TSystemDirectory* evtsDir = new TSystemDirectory(".",
370 presentDirName.Data());
371 TList* fileList = evtsDir->GetListOfFiles();
372 Int_t nFiles = fileList->GetEntries();
373 //cout<<" Int_t nFiles = "<<nFiles<<endl;
374 gSystem->cd(execDir);
375 for(Int_t iFiles=0; iFiles<nFiles; ++iFiles)
377 TSystemFile* presentFile = (TSystemFile*) fileList->At(iFiles);
378 TString presentFileName(presentDirName);
379 presentFileName += presentFile->GetName();
381 if(!(presentFileName.Contains("Kinematics") &&
382 presentFileName.Contains("root"))) { continue ; }
384 //cout << " found: " << presentFileName.Data() << endl;
386 TFile* kineFile = new TFile(presentFileName.Data(), "READ");
388 Int_t nEvts = kineFile->GetNkeys() ;
389 //cout << " . found: " << nEvts << " KineTree(s) in " << presentFileName.Data() << endl;
390 TList* kineEventsList = (TList*)kineFile->GetListOfKeys();
392 TIter next(kineEventsList);
396 //loop over the events
397 while( key=(TKey *)next() )
399 TDirectory* tDir = (TDirectory*)key->ReadObj();
402 TString evtDir(tDir->GetName());
403 //cout << " . . found: " << tDir->GetName() << endl;
405 kTree = (TTree *)tDir->Get("TreeK");
408 Int_t nPart = kTree->GetEntries();
409 //cout << " . . . kTree " << fCount << " has " << nPart << " particles " << endl;
411 //-----------------------------------------------------------
412 //fill and save the flow event
415 Int_t nNewMultOfEvent = random3Temp.Gaus(nMultiplicityOfEvent,nSigmaMult);
416 cout << "new multiplicity: " << nNewMultOfEvent << endl;
417 Double_t xNewFlowValue = random3Temp.Gaus(xEllipticFlowValue,xSigmaFlow);
418 if ( (fCount % 100) == 0) {
419 cout << "new multiplicity: " << nNewMultOfEvent << endl;
420 cout << "new flow value: " << xNewFlowValue << endl;
424 AliFlowEventSimple *fEvent = new AliFlowEventSimple(kTree, cutsInt, cutsDiff);
425 //xRPAngle=random3Temp.Uniform(0.0,TMath::TwoPi());
426 //fEvent->SetMCReactionPlaneAngle(xRPAngle);
427 //fEvent->SetV2(xNewFlowValue);
428 //fEvent->CloneTracks(nLoops);
430 // do flow analysis for various methods
431 if(MCEP) mcep->Make(fEvent);
432 if(QC) qc->Make(fEvent);
433 if(GFC) gfc->Make(fEvent);
434 if(FQD) fqd->Make(fEvent);
435 if(LYZ1SUM) lyz1sum->Make(fEvent);
436 if(LYZ1PROD)lyz1prod->Make(fEvent);
437 if(LYZ2SUM) lyz2sum->Make(fEvent);
438 if(LYZ2PROD)lyz2prod->Make(fEvent);
439 if(LYZEP) lyzep->Make(fEvent,ep);
440 if(SP) sp->Make(fEvent);
441 if(MH) mh->Make(fEvent);
442 if(NL) nl->Make(fEvent);
443 //-----------------------------------------------------------
445 //cout << "# " << fCount << " events processed" << endl;
454 //---------------------------------------------------------------------------------------
455 // create a new file which will hold the final results of all methods:
456 TString outputFileName = "AnalysisResults.root";
457 TFile *outputFile = new TFile(outputFileName.Data(),"RECREATE");
458 // create a new file for each method wich will hold list with final results:
459 const Int_t nMethods = 12;
460 TString method[nMethods] = {"MCEP","SP","GFC","QC","FQD","LYZ1SUM","LYZ1PROD","LYZ2SUM","LYZ2PROD","LYZEP","MH","NL"};
461 TDirectoryFile *dirFileFinal[nMethods] = {NULL};
462 TString fileNameMethod[nMethods];
463 for(Int_t i=0;i<nMethods;i++)
465 // form a file name for each method:
466 fileNameMethod[i]+="output";
467 fileNameMethod[i]+=method[i].Data();
468 fileNameMethod[i]+="analysis";
469 dirFileFinal[i] = new TDirectoryFile(fileNameMethod[i].Data(),fileNameMethod[i].Data());
472 // calculating and storing the final results of default flow analysis:
473 if(MCEP) {mcep->Finish(); mcep->WriteHistograms(dirFileFinal[0]);}
474 if(SP) {sp->Finish(); sp->WriteHistograms(dirFileFinal[1]);}
475 if(GFC) {gfc->Finish(); gfc->WriteHistograms(dirFileFinal[2]);}
476 if(QC) {qc->Finish(); qc->WriteHistograms(dirFileFinal[3]);}
477 if(FQD) {fqd->Finish(); fqd->WriteHistograms(dirFileFinal[4]);}
478 if(LYZ1SUM) {lyz1sum->Finish(); lyz1sum->WriteHistograms(dirFileFinal[5]);}
479 if(LYZ1PROD){lyz1prod->Finish();lyz1prod->WriteHistograms(dirFileFinal[6]);}
480 if(LYZ2SUM) {lyz2sum->Finish(); lyz2sum->WriteHistograms(dirFileFinal[7]);}
481 if(LYZ2PROD){lyz2prod->Finish();lyz2prod->WriteHistograms(dirFileFinal[8]);}
482 if(LYZEP) {lyzep->Finish(); lyzep->WriteHistograms(dirFileFinal[9]);}
483 if(MH) {mh->Finish(); mh->WriteHistograms(dirFileFinal[10]);}
484 if(NL) {nl->Finish(); nl->WriteHistograms(dirFileFinal[11]);}
485 //---------------------------------------------------------------------------------------
491 cout << " Fini ... " << endl;
499 void SetupPar(char* pararchivename)
501 //Load par files, create analysis libraries
502 //For testing, if par file already decompressed and modified
503 //classes then do not decompress.
505 TString cdir(Form("%s", gSystem->WorkingDirectory() )) ;
506 TString parpar(Form("%s.par", pararchivename)) ;
507 if ( gSystem->AccessPathName(parpar.Data()) ) {
508 gSystem->ChangeDirectory(gSystem->Getenv("ALICE_ROOT")) ;
509 TString processline(Form(".! make %s", parpar.Data())) ;
510 gROOT->ProcessLine(processline.Data()) ;
511 gSystem->ChangeDirectory(cdir) ;
512 processline = Form(".! mv /tmp/%s .", parpar.Data()) ;
513 gROOT->ProcessLine(processline.Data()) ;
515 if ( gSystem->AccessPathName(pararchivename) ) {
516 TString processline = Form(".! tar xvzf %s",parpar.Data()) ;
517 gROOT->ProcessLine(processline.Data());
520 TString ocwd = gSystem->WorkingDirectory();
521 gSystem->ChangeDirectory(pararchivename);
523 // check for BUILD.sh and execute
524 if (!gSystem->AccessPathName("PROOF-INF/BUILD.sh")) {
525 printf("*******************************\n");
526 printf("*** Building PAR archive ***\n");
527 cout<<pararchivename<<endl;
528 printf("*******************************\n");
530 if (gSystem->Exec("PROOF-INF/BUILD.sh")) {
531 Error("runProcess","Cannot Build the PAR Archive! - Abort!");
535 // check for SETUP.C and execute
536 if (!gSystem->AccessPathName("PROOF-INF/SETUP.C")) {
537 printf("*******************************\n");
538 printf("*** Setup PAR archive ***\n");
539 cout<<pararchivename<<endl;
540 printf("*******************************\n");
541 gROOT->Macro("PROOF-INF/SETUP.C");
544 gSystem->ChangeDirectory(ocwd.Data());
545 printf("Current dir: %s\n", ocwd.Data());
548 void LoadLibraries(const anaModes mode) {
550 //--------------------------------------
551 // Load the needed libraries most of them already loaded by aliroot
552 //--------------------------------------
553 //gSystem->Load("libTree");
554 gSystem->Load("libGeom");
555 gSystem->Load("libVMC");
556 gSystem->Load("libXMLIO");
557 gSystem->Load("libPhysics");
559 //----------------------------------------------------------
560 // >>>>>>>>>>> Local mode <<<<<<<<<<<<<<
561 //----------------------------------------------------------
563 //--------------------------------------------------------
564 // If you want to use already compiled libraries
565 // in the aliroot distribution
566 //--------------------------------------------------------
567 gSystem->Load("libSTEERBase");
568 gSystem->Load("libESD");
569 gSystem->Load("libAOD");
570 gSystem->Load("libANALYSIS");
571 gSystem->Load("libANALYSISalice");
572 gSystem->Load("libCORRFW");
573 cerr<<"libCORRFW loaded..."<<endl;
574 gSystem->Load("libPWGflowBase");
575 cerr<<"libPWGflowBase loaded..."<<endl;
576 gSystem->Load("libPWGflowTasks");
577 cerr<<"libPWGflowTasks loaded..."<<endl;
580 else if (mode == mLocalPAR) {
581 //--------------------------------------------------------
582 //If you want to use root and par files from aliroot
583 //--------------------------------------------------------
584 //If you want to use root and par files from aliroot
585 //--------------------------------------------------------
586 SetupPar("STEERBase");
589 SetupPar("ANALYSIS");
590 SetupPar("ANALYSISalice");
593 SetupPar("PWGflowBase");
594 cerr<<"PWGflowBase.par loaded..."<<endl;
595 SetupPar("PWGflowTasks");
596 cerr<<"PWGflowTasks.par loaded..."<<endl;
599 //---------------------------------------------------------
600 // <<<<<<<<<< Source mode >>>>>>>>>>>>
601 //---------------------------------------------------------
602 else if (mode==mLocalSource) {
604 // In root inline compile
608 gROOT->LoadMacro("Base/AliFlowCommonConstants.cxx+");
609 gROOT->LoadMacro("Base/AliFlowLYZConstants.cxx+");
612 gROOT->LoadMacro("Base/AliFlowVector.cxx+");
613 gROOT->LoadMacro("Base/AliFlowTrackSimple.cxx+");
614 gROOT->LoadMacro("Base/AliFlowTrackSimpleCuts.cxx+");
615 gROOT->LoadMacro("Base/AliFlowEventSimple.cxx+");
617 // Output histosgrams
618 gROOT->LoadMacro("Base/AliFlowCommonHist.cxx+");
619 gROOT->LoadMacro("Base/AliFlowCommonHistResults.cxx+");
620 gROOT->LoadMacro("Base/AliFlowLYZHist1.cxx+");
621 gROOT->LoadMacro("Base/AliFlowLYZHist2.cxx+");
623 // Functions needed for various methods
624 gROOT->LoadMacro("Base/AliCumulantsFunctions.cxx+");
625 gROOT->LoadMacro("Base/AliFlowLYZEventPlane.cxx+");
627 // Flow Analysis code for various methods
628 gROOT->LoadMacro("Base/AliFlowAnalysisWithMCEventPlane.cxx+");
629 gROOT->LoadMacro("Base/AliFlowAnalysisWithScalarProduct.cxx+");
630 gROOT->LoadMacro("Base/AliFlowAnalysisWithLYZEventPlane.cxx+");
631 gROOT->LoadMacro("Base/AliFlowAnalysisWithLeeYangZeros.cxx+");
632 gROOT->LoadMacro("Base/AliFlowAnalysisWithCumulants.cxx+");
633 gROOT->LoadMacro("Base/AliFlowAnalysisWithQCumulants.cxx+");
634 gROOT->LoadMacro("Base/AliFlowAnalysisWithFittingQDistribution.cxx+");
635 gROOT->LoadMacro("Base/AliFlowAnalysisWithMixedHarmonics.cxx+");
636 gROOT->LoadMacro("Base/AliFlowAnalysisWithNestedLoops.cxx+");
638 cout << "finished loading macros!" << endl;