combined set of macros to merge the output

[u/mrichter/AliRoot.git] / PWG2 / FLOW / macros / showSpread.C

enum libModes {mLocal,mLocalSource};

void showSpread(TString type="", const Int_t nRuns=-1, Int_t mode=mLocal)
{
 // type:  type of analysis can be ESD, AOD, MC, ESDMC0, ESDMC1
 //        (if type="" output files are from MC simulation (default))
 // nRuns: specify here from how many small statistics output .root files          
 //        flow estimates will be shown (if nRuns = -1 all of them will be shown)
 // mode:  if mode=mLocal analyze data on your computer using aliroot
 //        if mode=mLocalSource analyze data on your computer using root + source files
 
 // load needed libraries:                       
 LoadSpreadLibraries(mode);  
 
 // access the path of current directory:
 TString pwd(gSystem->pwd());
 pwd+="/"; 
 
 // standard magic:
 TString execDir(gSystem->pwd());  
 TSystemDirectory* baseDir = new TSystemDirectory(".",execDir.Data());          
 TList* dirList = baseDir->GetListOfFiles();
 Int_t nDirs = dirList->GetEntries();
 gSystem->cd(execDir);
 
 Int_t countDirectories = 0;
 if(dirList)
 {
  for(Int_t i=0;i<nDirs;i++)
  {
   if((dirList->At(i))->InheritsFrom("TSystemDirectory") &&  
      !(strcmp((dirList->At(i))->GetName(),".") == 0 ) &&  
      !(strcmp((dirList->At(i))->GetName(),"..") == 0 )) 
   {
    countDirectories++;
   } 
  }
 }   
 
 // assigning bin numbers to methods:
 Int_t binMCEP = 1; 
 Int_t binSP = 2;
 Int_t binGFC2 = 3; 
 Int_t binGFC4 = 5; 
 Int_t binGFC6 = 7; 
 Int_t binGFC8 = 9; 
 Int_t binQC2 = 4; 
 Int_t binQC4 = 6; 
 Int_t binQC6 = 8; 
 Int_t binQC8 = 10; 
 Int_t binFQD = 11; 
 Int_t binLYZ1 = 12; 
 Int_t binLYZEP = 13;
 
 // one subdirectory gives one estimate for each method: 
 const Int_t nEstimates = countDirectories;
 
 // arrays to store estimates and errors of each method from different small statistics runs:
 // MCEP:
 Double_t mcepValueNONAME[nEstimates] = {0.}; 
 Double_t mcepErrorNONAME[nEstimates] = {0.}; 
 Double_t mcepMaxValueNONAME = 0.; // to be improved      
 Double_t mcepMinValueNONAME = 1000.; // to be improved
 
 // SP:
 Double_t spValueNONAME[nEstimates] = {0.}; 
 Double_t spErrorNONAME[nEstimates] = {0.}; 
 Double_t spMaxValueNONAME = 0.; // to be improved      
 Double_t spMinValueNONAME = 1000.; // to be improved
 
 // GFC{2}:
 Double_t gfc2ValueNONAME[nEstimates] = {0.}; 
 Double_t gfc2ErrorNONAME[nEstimates] = {0.}; 
 Double_t gfc2MaxValueNONAME = 0.; // to be improved      
 Double_t gfc2MinValueNONAME = 1000.; // to be improved
 
 // GFC{4}:
 Double_t gfc4ValueNONAME[nEstimates] = {0.}; 
 Double_t gfc4ErrorNONAME[nEstimates] = {0.}; 
 Double_t gfc4MaxValueNONAME = 0.; // to be improved      
 Double_t gfc4MinValueNONAME = 1000.; // to be improved
    
 // GFC{6}:
 Double_t gfc6ValueNONAME[nEstimates] = {0.}; 
 Double_t gfc6ErrorNONAME[nEstimates] = {0.}; 
 Double_t gfc6MaxValueNONAME = 0.; // to be improved      
 Double_t gfc6MinValueNONAME = 1000.; // to be improved

 // GFC{8}:
 Double_t gfc8ValueNONAME[nEstimates] = {0.}; 
 Double_t gfc8ErrorNONAME[nEstimates] = {0.}; 
 Double_t gfc8MaxValueNONAME = 0.; // to be improved      
 Double_t gfc8MinValueNONAME = 1000.; // to be improved
 
 // QC{2}:
 Double_t qc2ValueNONAME[nEstimates] = {0.}; 
 Double_t qc2ErrorNONAME[nEstimates] = {0.}; 
 Double_t qc2MaxValueNONAME = 0.; // to be improved      
 Double_t qc2MinValueNONAME = 1000.; // to be improved
 
 // QC{4}:
 Double_t qc4ValueNONAME[nEstimates] = {0.}; 
 Double_t qc4ErrorNONAME[nEstimates] = {0.}; 
 Double_t qc4MaxValueNONAME = 0.; // to be improved      
 Double_t qc4MinValueNONAME = 1000.; // to be improved
    
 // QC{6}:
 Double_t qc6ValueNONAME[nEstimates] = {0.}; 
 Double_t qc6ErrorNONAME[nEstimates] = {0.}; 
 Double_t qc6MaxValueNONAME = 0.; // to be improved      
 Double_t qc6MinValueNONAME = 1000.; // to be improved

 // QC{8}:
 Double_t qc8ValueNONAME[nEstimates] = {0.}; 
 Double_t qc8ErrorNONAME[nEstimates] = {0.}; 
 Double_t qc8MaxValueNONAME = 0.; // to be improved      
 Double_t qc8MinValueNONAME = 1000.; // to be improved
 
 // FQD:
 Double_t fqdValueNONAME[nEstimates] = {0.}; 
 Double_t fqdErrorNONAME[nEstimates] = {0.}; 
 Double_t fqdMaxValueNONAME = 0.; // to be improved      
 Double_t fqdMinValueNONAME = 1000.; // to be improved
       
 // LYZ1:
 Double_t lyz1ValueNONAME[nEstimates] = {0.}; 
 Double_t lyz1ErrorNONAME[nEstimates] = {0.}; 
 Double_t lyz1MaxValueNONAME = 0.; // to be improved      
 Double_t lyz1MinValueNONAME = 1000.; // to be improved
             
 // LYZEP:
 Double_t lyzepValueNONAME[nEstimates] = {0.}; 
 Double_t lyzepErrorNONAME[nEstimates] = {0.}; 
 Double_t lyzepMaxValueNONAME = 0.; // to be improved      
 Double_t lyzepMinValueNONAME = 1000.; // to be improved        
             
 Int_t counter = 0;
  
 for(Int_t iDir=0;iDir<nDirs;++iDir)
 {
  TSystemFile* presentDir = (TSystemFile*)dirList->At(iDir);
  if(!presentDir || !presentDir->IsDirectory() || strcmp(presentDir->GetName(), ".") == 0 || 
     strcmp(presentDir->GetName(), "..") == 0) continue; 
  
  if(!(nRuns == -1))
  {                
   if (counter >= nRuns) break;       
  }
                                            
  TString presentDirName(gSystem->pwd()); 
  presentDirName += "/";
  presentDirName += presentDir->GetName();
  presentDirName += "/";
   
  // accessing the small statistics output .root files for each method:
  // MCEP:     
  TString fileNameMCEP = presentDirName;   
  fileNameMCEP+="outputMCEPanalysis";
  (fileNameMCEP+=type.Data())+=".root";
  TFile *fileMCEP = TFile::Open(fileNameMCEP.Data(), "READ");      
  if(fileMCEP) 
  {
   TList *listMCEP = NULL;
   AliFlowCommonHistResults *mcepCommonHistRes = NULL; 
   fileMCEP->GetObject("cobjMCEP",listMCEP); 
   fileMCEP->Close();
   if(listMCEP) 
   {
    mcepCommonHistRes = dynamic_cast<AliFlowCommonHistResults*> (listMCEP->FindObject("AliFlowCommonHistResultsMCEP")); 
    if(mcepCommonHistRes && mcepCommonHistRes->GetHistIntFlow())
    {
     mcepValueNONAME[counter] = (mcepCommonHistRes->GetHistIntFlow())->GetBinContent(1);
     mcepErrorNONAME[counter] = (mcepCommonHistRes->GetHistIntFlow())->GetBinError(1);
     if(mcepValueNONAME[counter]>0.) // to be improved 
     {
      if(mcepMaxValueNONAME < mcepValueNONAME[counter]) mcepMaxValueNONAME = mcepValueNONAME[counter]; 
      if(mcepMinValueNONAME > mcepValueNONAME[counter]) mcepMinValueNONAME = mcepValueNONAME[counter]; 
     } 
    }
   } // end of if(listMCEP)
  } // end of if(fileMCEP) 
  
  // SP:     
  TString fileNameSP = presentDirName;   
  fileNameSP+="outputSPanalysis";
  (fileNameSP+=type.Data())+=".root";
  TFile *fileSP = TFile::Open(fileNameSP.Data(), "READ");      
  if(fileSP) 
  {
   TList *listSP = NULL;
   AliFlowCommonHistResults *spCommonHistRes = NULL; 
   fileSP->GetObject("cobjSP",listSP); 
   fileSP->Close();
   if(listSP) 
   {
    spCommonHistRes = dynamic_cast<AliFlowCommonHistResults*> (listSP->FindObject("AliFlowCommonHistResultsSP")); 
    if(spCommonHistRes && spCommonHistRes->GetHistIntFlow())
    {
     spValueNONAME[counter] = (spCommonHistRes->GetHistIntFlow())->GetBinContent(1);
     spErrorNONAME[counter] = (spCommonHistRes->GetHistIntFlow())->GetBinError(1);
     if(spValueNONAME[counter]>0.) // to be improved 
     {
      if(spMaxValueNONAME < spValueNONAME[counter]) spMaxValueNONAME = spValueNONAME[counter]; 
      if(spMinValueNONAME > spValueNONAME[counter]) spMinValueNONAME = spValueNONAME[counter]; 
     } 
    }
   } // end of if(listSP)
  } // end of if(fileSP) 
  
  // GFC:     
  TString fileNameGFC = presentDirName;   
  fileNameGFC+="outputGFCanalysis";
  (fileNameGFC+=type.Data())+=".root";
  TFile *fileGFC = TFile::Open(fileNameGFC.Data(), "READ");      
  if(fileGFC) 
  {
   TList *listGFC = NULL;
   AliFlowCommonHistResults *gfc2CommonHistRes = NULL; 
   AliFlowCommonHistResults *gfc4CommonHistRes = NULL; 
   AliFlowCommonHistResults *gfc6CommonHistRes = NULL; 
   AliFlowCommonHistResults *gfc8CommonHistRes = NULL; 
   fileGFC->GetObject("cobjGFC",listGFC); 
   fileGFC->Close();
   if(listGFC) 
   {
    gfc2CommonHistRes = dynamic_cast<AliFlowCommonHistResults*> (listGFC->FindObject("AliFlowCommonHistResults2ndOrderGFC"));
    gfc4CommonHistRes = dynamic_cast<AliFlowCommonHistResults*> (listGFC->FindObject("AliFlowCommonHistResults4thOrderGFC"));
    gfc6CommonHistRes = dynamic_cast<AliFlowCommonHistResults*> (listGFC->FindObject("AliFlowCommonHistResults6thOrderGFC"));
    gfc8CommonHistRes = dynamic_cast<AliFlowCommonHistResults*> (listGFC->FindObject("AliFlowCommonHistResults8thOrderGFC")); 
    if(gfc2CommonHistRes && gfc2CommonHistRes->GetHistIntFlow())
    {
     gfc2ValueNONAME[counter] = (gfc2CommonHistRes->GetHistIntFlow())->GetBinContent(1);
     gfc2ErrorNONAME[counter] = (gfc2CommonHistRes->GetHistIntFlow())->GetBinError(1);
     if(gfc2ValueNONAME[counter]>0.) // to be improved 
     {
      if(gfc2MaxValueNONAME < gfc2ValueNONAME[counter]) gfc2MaxValueNONAME = gfc2ValueNONAME[counter]; 
      if(gfc2MinValueNONAME > gfc2ValueNONAME[counter]) gfc2MinValueNONAME = gfc2ValueNONAME[counter]; 
     } 
    }
    if(gfc4CommonHistRes && gfc4CommonHistRes->GetHistIntFlow())
    {
     gfc4ValueNONAME[counter] = (gfc4CommonHistRes->GetHistIntFlow())->GetBinContent(1);
     gfc4ErrorNONAME[counter] = (gfc4CommonHistRes->GetHistIntFlow())->GetBinError(1);
     if(gfc4ValueNONAME[counter]>0.) // to be improved 
     {
      if(gfc4MaxValueNONAME < gfc4ValueNONAME[counter]) gfc4MaxValueNONAME = gfc4ValueNONAME[counter]; 
      if(gfc4MinValueNONAME > gfc4ValueNONAME[counter]) gfc4MinValueNONAME = gfc4ValueNONAME[counter]; 
     } 
    }
    if(gfc6CommonHistRes && gfc6CommonHistRes->GetHistIntFlow())
    {
     gfc6ValueNONAME[counter] = (gfc6CommonHistRes->GetHistIntFlow())->GetBinContent(1);
     gfc6ErrorNONAME[counter] = (gfc6CommonHistRes->GetHistIntFlow())->GetBinError(1);
     if(gfc6ValueNONAME[counter]>0.) // to be improved 
     {
      if(gfc6MaxValueNONAME < gfc6ValueNONAME[counter]) gfc6MaxValueNONAME = gfc6ValueNONAME[counter]; 
      if(gfc6MinValueNONAME > gfc6ValueNONAME[counter]) gfc6MinValueNONAME = gfc6ValueNONAME[counter]; 
     } 
    }
    if(gfc8CommonHistRes && gfc8CommonHistRes->GetHistIntFlow())
    {
     gfc8ValueNONAME[counter] = (gfc8CommonHistRes->GetHistIntFlow())->GetBinContent(1);
     gfc8ErrorNONAME[counter] = (gfc8CommonHistRes->GetHistIntFlow())->GetBinError(1);
     if(gfc8ValueNONAME[counter]>0.) // to be improved 
     {
      if(gfc8MaxValueNONAME < gfc8ValueNONAME[counter]) gfc8MaxValueNONAME = gfc8ValueNONAME[counter]; 
      if(gfc8MinValueNONAME > gfc8ValueNONAME[counter]) gfc8MinValueNONAME = gfc8ValueNONAME[counter]; 
     } 
    }
   } // end of if(listGFC)
  } // end of if(fileGFC) 
  
  // QC:     
  TString fileNameQC = presentDirName;   
  fileNameQC+="outputQCanalysis";
  (fileNameQC+=type.Data())+=".root";
  TFile *fileQC = TFile::Open(fileNameQC.Data(), "READ");      
  if(fileQC) 
  {
   TList *listQC = NULL;
   AliFlowCommonHistResults *qc2CommonHistRes = NULL; 
   AliFlowCommonHistResults *qc4CommonHistRes = NULL; 
   AliFlowCommonHistResults *qc6CommonHistRes = NULL; 
   AliFlowCommonHistResults *qc8CommonHistRes = NULL; 
   fileQC->GetObject("cobjQC",listQC); 
   fileQC->Close();
   if(listQC) 
   {
    qc2CommonHistRes = dynamic_cast<AliFlowCommonHistResults*> (listQC->FindObject("AliFlowCommonHistResults2ndOrderQC"));
    qc4CommonHistRes = dynamic_cast<AliFlowCommonHistResults*> (listQC->FindObject("AliFlowCommonHistResults4thOrderQC"));
    qc6CommonHistRes = dynamic_cast<AliFlowCommonHistResults*> (listQC->FindObject("AliFlowCommonHistResults6thOrderQC"));
    qc8CommonHistRes = dynamic_cast<AliFlowCommonHistResults*> (listQC->FindObject("AliFlowCommonHistResults8thOrderQC")); 
    if(qc2CommonHistRes && qc2CommonHistRes->GetHistIntFlow())
    {
     qc2ValueNONAME[counter] = (qc2CommonHistRes->GetHistIntFlow())->GetBinContent(1);
     qc2ErrorNONAME[counter] = (qc2CommonHistRes->GetHistIntFlow())->GetBinError(1);
     if(qc2ValueNONAME[counter]>0.) // to be improved 
     {
      if(qc2MaxValueNONAME < qc2ValueNONAME[counter]) qc2MaxValueNONAME = qc2ValueNONAME[counter]; 
      if(qc2MinValueNONAME > qc2ValueNONAME[counter]) qc2MinValueNONAME = qc2ValueNONAME[counter]; 
     } 
    }
    if(qc4CommonHistRes && qc4CommonHistRes->GetHistIntFlow())
    {
     qc4ValueNONAME[counter] = (qc4CommonHistRes->GetHistIntFlow())->GetBinContent(1);
     qc4ErrorNONAME[counter] = (qc4CommonHistRes->GetHistIntFlow())->GetBinError(1);
     if(qc4ValueNONAME[counter]>0.) // to be improved 
     {
      if(qc4MaxValueNONAME < qc4ValueNONAME[counter]) qc4MaxValueNONAME = qc4ValueNONAME[counter]; 
      if(qc4MinValueNONAME > qc4ValueNONAME[counter]) qc4MinValueNONAME = qc4ValueNONAME[counter]; 
     } 
    }
    if(qc6CommonHistRes && qc6CommonHistRes->GetHistIntFlow())
    {
     qc6ValueNONAME[counter] = (qc6CommonHistRes->GetHistIntFlow())->GetBinContent(1);
     qc6ErrorNONAME[counter] = (qc6CommonHistRes->GetHistIntFlow())->GetBinError(1);
     if(qc6ValueNONAME[counter]>0.) // to be improved 
     {
      if(qc6MaxValueNONAME < qc6ValueNONAME[counter]) qc6MaxValueNONAME = qc6ValueNONAME[counter]; 
      if(qc6MinValueNONAME > qc6ValueNONAME[counter]) qc6MinValueNONAME = qc6ValueNONAME[counter]; 
     } 
    }
    if(qc8CommonHistRes && qc8CommonHistRes->GetHistIntFlow())
    {
     qc8ValueNONAME[counter] = (qc8CommonHistRes->GetHistIntFlow())->GetBinContent(1);
     qc8ErrorNONAME[counter] = (qc8CommonHistRes->GetHistIntFlow())->GetBinError(1);
     if(qc8ValueNONAME[counter]>0.) // to be improved 
     {
      if(qc8MaxValueNONAME < qc8ValueNONAME[counter]) qc8MaxValueNONAME = qc8ValueNONAME[counter]; 
      if(qc8MinValueNONAME > qc8ValueNONAME[counter]) qc8MinValueNONAME = qc8ValueNONAME[counter]; 
     } 
    }
   } // end of if(listQC)
  } // end of if(fileQC) 
 
  // FQD:     
  TString fileNameFQD = presentDirName;   
  fileNameFQD+="outputFQDanalysis";
  (fileNameFQD+=type.Data())+=".root";
  TFile *fileFQD = TFile::Open(fileNameFQD.Data(), "READ");      
  if(fileFQD) 
  {
   TList *listFQD = NULL;
   AliFlowCommonHistResults *fqdCommonHistRes = NULL; 
   fileFQD->GetObject("cobjFQD",listFQD); 
   fileFQD->Close();
   if(listFQD) 
   {
    fqdCommonHistRes = dynamic_cast<AliFlowCommonHistResults*> (listFQD->FindObject("AliFlowCommonHistResultsFQD")); 
    if(fqdCommonHistRes && fqdCommonHistRes->GetHistIntFlow())
    {
     fqdValueNONAME[counter] = (fqdCommonHistRes->GetHistIntFlow())->GetBinContent(1);
     fqdErrorNONAME[counter] = (fqdCommonHistRes->GetHistIntFlow())->GetBinError(1);
     if(fqdValueNONAME[counter]>0.) // to be improved 
     {
      if(fqdMaxValueNONAME < fqdValueNONAME[counter]) fqdMaxValueNONAME = fqdValueNONAME[counter]; 
      if(fqdMinValueNONAME > fqdValueNONAME[counter]) fqdMinValueNONAME = fqdValueNONAME[counter]; 
     } 
    }
   } // end of if(listFQD)
  } // end of if(fileFQD)   
  
 // ... LYZ* to be added ... 
 
  counter++;
  
 } // end of for(Int_t iDir=0;iDir<nDirs;++iDir)
 
 // accessing the large statistics merged output .root files for each method:
 // MCEP:
 TString mergedOutputFileNameMCEP(pwd.Data());
 ((mergedOutputFileNameMCEP+="outputMCEPanalysis")+=type.Data())+=".root";
 TFile *mergedOutputFileMCEP = NULL;
 TList *mergedOutputListMCEP = NULL;
 Double_t mergedValueMCEP = 0.;
 Double_t mergedErrorMCEP = 0.; 
 if(gSystem->AccessPathName(mergedOutputFileNameMCEP.Data(),kFileExists))
 {
  cout<<"WARNING: You do not have a merged output file "<<mergedOutputFileNameMCEP.Data()<<endl;
 } else 
   {     
    mergedOutputFileMCEP = TFile::Open(mergedOutputFileNameMCEP.Data(),"READ");      
    if(mergedOutputFileMCEP) 
    { 
     mergedOutputFileMCEP->GetObject("cobjMCEP",mergedOutputListMCEP); 
     mergedOutputFileMCEP->Close();
     if(mergedOutputListMCEP) 
     {
      AliFlowCommonHistResults *mcepCommonHistRes = dynamic_cast<AliFlowCommonHistResults*> 
                                                    (mergedOutputListMCEP->FindObject("AliFlowCommonHistResultsMCEP")); 
      if(mcepCommonHistRes && mcepCommonHistRes->GetHistIntFlow())
      {
       mergedValueMCEP = (mcepCommonHistRes->GetHistIntFlow())->GetBinContent(1);
       mergedErrorMCEP = (mcepCommonHistRes->GetHistIntFlow())->GetBinError(1);
      } 
     }
    } // end of if(mergedOutputFileMCEP)
   } // end of else 
   
 // SP:
 TString mergedOutputFileNameSP(pwd.Data());
 ((mergedOutputFileNameSP+="outputSPanalysis")+=type.Data())+=".root";
 TFile *mergedOutputFileSP = NULL;
 TList *mergedOutputListSP = NULL;
 Double_t mergedValueSP = 0.;
 Double_t mergedErrorSP = 0.; 
 if(gSystem->AccessPathName(mergedOutputFileNameSP.Data(),kFileExists))
 {
  cout<<"WARNING: You do not have a merged output file "<<mergedOutputFileNameSP.Data()<<endl;
 } else 
   {     
    mergedOutputFileSP = TFile::Open(mergedOutputFileNameSP.Data(),"READ");      
    if(mergedOutputFileSP) 
    { 
     mergedOutputFileSP->GetObject("cobjSP",mergedOutputListSP); 
     mergedOutputFileSP->Close();
     if(mergedOutputListSP) 
     {
      AliFlowCommonHistResults *spCommonHistRes = dynamic_cast<AliFlowCommonHistResults*> 
                                                    (mergedOutputListSP->FindObject("AliFlowCommonHistResultsSP")); 
      if(spCommonHistRes && spCommonHistRes->GetHistIntFlow())
      {
       mergedValueSP = (spCommonHistRes->GetHistIntFlow())->GetBinContent(1);
       mergedErrorSP = (spCommonHistRes->GetHistIntFlow())->GetBinError(1);
      } 
     }
    } // end of if(mergedOutputFileSP)
   } // end of else 
 
 // GFC:
 TString mergedOutputFileNameGFC(pwd.Data());
 ((mergedOutputFileNameGFC+="outputGFCanalysis")+=type.Data())+=".root";
 TFile *mergedOutputFileGFC = NULL;
 TList *mergedOutputListGFC = NULL;
 Double_t mergedValueGFC2 = 0.;
 Double_t mergedErrorGFC2 = 0.; 
 Double_t mergedValueGFC4 = 0.;
 Double_t mergedErrorGFC4 = 0.; 
 Double_t mergedValueGFC6 = 0.;
 Double_t mergedErrorGFC6 = 0.; 
 Double_t mergedValueGFC8 = 0.;
 Double_t mergedErrorGFC8 = 0.; 
 if(gSystem->AccessPathName(mergedOutputFileNameGFC.Data(),kFileExists))
 {
  cout<<"WARNING: You do not have a merged output file "<<mergedOutputFileNameGFC.Data()<<endl;
 } else 
   {     
    mergedOutputFileGFC = TFile::Open(mergedOutputFileNameGFC.Data(),"READ");      
    if(mergedOutputFileGFC) 
    { 
     mergedOutputFileGFC->GetObject("cobjGFC",mergedOutputListGFC); 
     mergedOutputFileGFC->Close();
     if(mergedOutputListGFC) 
     {
      AliFlowCommonHistResults *gfc2CommonHistRes = dynamic_cast<AliFlowCommonHistResults*> 
                                                    (mergedOutputListGFC->FindObject("AliFlowCommonHistResults2ndOrderGFC")); 
      AliFlowCommonHistResults *gfc4CommonHistRes = dynamic_cast<AliFlowCommonHistResults*> 
                                                    (mergedOutputListGFC->FindObject("AliFlowCommonHistResults4thOrderGFC")); 
      AliFlowCommonHistResults *gfc6CommonHistRes = dynamic_cast<AliFlowCommonHistResults*> 
                                                    (mergedOutputListGFC->FindObject("AliFlowCommonHistResults6thOrderGFC")); 
      AliFlowCommonHistResults *gfc8CommonHistRes = dynamic_cast<AliFlowCommonHistResults*> 
                                                    (mergedOutputListGFC->FindObject("AliFlowCommonHistResults8thOrderGFC")); 
      if(gfc2CommonHistRes && gfc2CommonHistRes->GetHistIntFlow())
      {
       mergedValueGFC2 = (gfc2CommonHistRes->GetHistIntFlow())->GetBinContent(1);
       mergedErrorGFC2 = (gfc2CommonHistRes->GetHistIntFlow())->GetBinError(1);
      } 
      if(gfc4CommonHistRes && gfc4CommonHistRes->GetHistIntFlow())
      {
       mergedValueGFC4 = (gfc4CommonHistRes->GetHistIntFlow())->GetBinContent(1);
       mergedErrorGFC4 = (gfc4CommonHistRes->GetHistIntFlow())->GetBinError(1);
      }
      if(gfc6CommonHistRes && gfc6CommonHistRes->GetHistIntFlow())
      {
       mergedValueGFC6 = (gfc6CommonHistRes->GetHistIntFlow())->GetBinContent(1);
       mergedErrorGFC6 = (gfc6CommonHistRes->GetHistIntFlow())->GetBinError(1);
      } 
      if(gfc8CommonHistRes && gfc8CommonHistRes->GetHistIntFlow())
      {
       mergedValueGFC8 = (gfc8CommonHistRes->GetHistIntFlow())->GetBinContent(1);
       mergedErrorGFC8 = (gfc8CommonHistRes->GetHistIntFlow())->GetBinError(1);
      }  
     }
    } // end of if(mergedOutputFileGFC)
   } // end of else  

 // QC:
 TString mergedOutputFileNameQC(pwd.Data());
 ((mergedOutputFileNameQC+="outputQCanalysis")+=type.Data())+=".root";
 TFile *mergedOutputFileQC = NULL;
 TList *mergedOutputListQC = NULL;
 Double_t mergedValueQC2 = 0.;
 Double_t mergedErrorQC2 = 0.; 
 Double_t mergedValueQC4 = 0.;
 Double_t mergedErrorQC4 = 0.; 
 Double_t mergedValueQC6 = 0.;
 Double_t mergedErrorQC6 = 0.; 
 Double_t mergedValueQC8 = 0.;
 Double_t mergedErrorQC8 = 0.; 
 if(gSystem->AccessPathName(mergedOutputFileNameQC.Data(),kFileExists))
 {
  cout<<"WARNING: You do not have a merged output file "<<mergedOutputFileNameQC.Data()<<endl;
 } else 
   {     
    mergedOutputFileQC = TFile::Open(mergedOutputFileNameQC.Data(),"READ");      
    if(mergedOutputFileQC) 
    { 
     mergedOutputFileQC->GetObject("cobjQC",mergedOutputListQC); 
     mergedOutputFileQC->Close();
     if(mergedOutputListQC) 
     {
      AliFlowCommonHistResults *qc2CommonHistRes = dynamic_cast<AliFlowCommonHistResults*> 
                                                    (mergedOutputListQC->FindObject("AliFlowCommonHistResults2ndOrderQC")); 
      AliFlowCommonHistResults *qc4CommonHistRes = dynamic_cast<AliFlowCommonHistResults*> 
                                                    (mergedOutputListQC->FindObject("AliFlowCommonHistResults4thOrderQC")); 
      AliFlowCommonHistResults *qc6CommonHistRes = dynamic_cast<AliFlowCommonHistResults*> 
                                                    (mergedOutputListQC->FindObject("AliFlowCommonHistResults6thOrderQC")); 
      AliFlowCommonHistResults *qc8CommonHistRes = dynamic_cast<AliFlowCommonHistResults*> 
                                                    (mergedOutputListQC->FindObject("AliFlowCommonHistResults8thOrderQC")); 
      if(qc2CommonHistRes && qc2CommonHistRes->GetHistIntFlow())
      {
       mergedValueQC2 = (qc2CommonHistRes->GetHistIntFlow())->GetBinContent(1);
       mergedErrorQC2 = (qc2CommonHistRes->GetHistIntFlow())->GetBinError(1);
      } 
      if(qc4CommonHistRes && qc4CommonHistRes->GetHistIntFlow())
      {
       mergedValueQC4 = (qc4CommonHistRes->GetHistIntFlow())->GetBinContent(1);
       mergedErrorQC4 = (qc4CommonHistRes->GetHistIntFlow())->GetBinError(1);
      }
      if(qc6CommonHistRes && qc6CommonHistRes->GetHistIntFlow())
      {
       mergedValueQC6 = (qc6CommonHistRes->GetHistIntFlow())->GetBinContent(1);
       mergedErrorQC6 = (qc6CommonHistRes->GetHistIntFlow())->GetBinError(1);
      } 
      if(qc8CommonHistRes && qc8CommonHistRes->GetHistIntFlow())
      {
       mergedValueQC8 = (qc8CommonHistRes->GetHistIntFlow())->GetBinContent(1);
       mergedErrorQC8 = (qc8CommonHistRes->GetHistIntFlow())->GetBinError(1);
      }  
     }
    } // end of if(mergedOutputFileGFC)
   } // end of else  

 // FQD:
 TString mergedOutputFileNameFQD(pwd.Data());
 ((mergedOutputFileNameFQD+="outputFQDanalysis")+=type.Data())+=".root";
 TFile *mergedOutputFileFQD = NULL;
 TList *mergedOutputListFQD = NULL;
 Double_t mergedValueFQD = 0.;
 Double_t mergedErrorFQD = 0.; 
 if(gSystem->AccessPathName(mergedOutputFileNameFQD.Data(),kFileExists))
 {
  cout<<"WARNING: You do not have a merged output file "<<mergedOutputFileNameFQD.Data()<<endl;
 } else 
   {     
    mergedOutputFileFQD = TFile::Open(mergedOutputFileNameFQD.Data(),"READ");      
    if(mergedOutputFileFQD) 
    { 
     mergedOutputFileFQD->GetObject("cobjFQD",mergedOutputListFQD); 
     mergedOutputFileFQD->Close();
     if(mergedOutputListFQD) 
     {
      AliFlowCommonHistResults *fqdCommonHistRes = dynamic_cast<AliFlowCommonHistResults*> 
                                                    (mergedOutputListFQD->FindObject("AliFlowCommonHistResultsFQD")); 
      if(fqdCommonHistRes && fqdCommonHistRes->GetHistIntFlow())
      {
       mergedValueFQD = (fqdCommonHistRes->GetHistIntFlow())->GetBinContent(1);
       mergedErrorFQD = (fqdCommonHistRes->GetHistIntFlow())->GetBinError(1);
      } 
     }
    } // end of if(mergedOutputFileFQD)
   } // end of else  
   
 // ... LYZ* to be added ... 
                   
 // removing the title and stat. box from all histograms:
 // gStyle->SetOptTitle(0);
 gStyle->SetOptStat(0);  
 
 // box width:
 const Double_t boxWidth = 0.25;
  
 // the number of different methods:
 const Int_t nMethods = 13;
 
 // the number of small statistics runs:
 const Int_t nPoints = counter;
   
 // booking the style histogram for the integrated flow results from all methods for NONAME, RPs and POIs:
 TH1D* intFlowAll = new TH1D("intFlowAll","Integrated Flow",nMethods,0,nMethods);      
 //intFlowAll->SetLabelSize(0.036,"X");
 //intFlowAll->SetLabelSize(0.036,"Y");
 //intFlowAll->SetMarkerStyle(markerStyle);
 //intFlowAll->SetMarkerColor(markerColor);
 (intFlowAll->GetYaxis())->SetTickLength(0.01);
 (intFlowAll->GetXaxis())->SetBinLabel(binMCEP,"v_{2}{MC}");
 (intFlowAll->GetXaxis())->SetBinLabel(binSP,"v_{2}{SP}");
 (intFlowAll->GetXaxis())->SetBinLabel(binGFC2,"v_{2}{2,GFC}");
 (intFlowAll->GetXaxis())->SetBinLabel(binQC2,"v_{2}{2,QC}");
 (intFlowAll->GetXaxis())->SetBinLabel(binGFC4,"v_{2}{4,GFC}");
 (intFlowAll->GetXaxis())->SetBinLabel(binQC4,"v_{2}{4,QC}");
 (intFlowAll->GetXaxis())->SetBinLabel(binGFC6,"v_{2}{6,GFC}");
 (intFlowAll->GetXaxis())->SetBinLabel(binQC6,"v_{2}{6,QC}");
 (intFlowAll->GetXaxis())->SetBinLabel(binGFC8,"v_{2}{8,GFC}");
 (intFlowAll->GetXaxis())->SetBinLabel(binQC8,"v_{2}{8,QC}");
 (intFlowAll->GetXaxis())->SetBinLabel(binFQD,"v_{2}{FQD}");
 (intFlowAll->GetXaxis())->SetBinLabel(binLYZ1,"v_{2}{LYZ}");
 (intFlowAll->GetXaxis())->SetBinLabel(binLYZEP,"v_{2}{LYZEP}");
 
 //=============================================================================
 //                             NONAME histogram
 //=============================================================================
   
 TH1D *intFlowNONAME = new TH1D(*intFlowAll); 
 
 Double_t mcepNONAME[nPoints] = {0.};  
 Double_t spNONAME[nPoints] = {0.}; 
 Double_t gfc2NONAME[nPoints] = {0.};
 Double_t gfc4NONAME[nPoints] = {0.};
 Double_t gfc6NONAME[nPoints] = {0.};
 Double_t gfc8NONAME[nPoints] = {0.};
 Double_t qc2NONAME[nPoints] = {0.};
 Double_t qc4NONAME[nPoints] = {0.};
 Double_t qc6NONAME[nPoints] = {0.};
 Double_t qc8NONAME[nPoints] = {0.};
 Double_t fqdNONAME[nPoints] = {0.}; 
 Double_t lyz1NONAME[nPoints] = {0.}; 
 Double_t lyzepNONAME[nPoints] = {0.}; 
 
 for(Int_t i=0;i<nPoints;i++)
 {
  mcepNONAME[i]=binMCEP-0.5;
  spNONAME[i]=binSP-0.5;
  gfc2NONAME[i]=binGFC2-0.5;
  gfc4NONAME[i]=binGFC4-0.5;
  gfc6NONAME[i]=binGFC6-0.5;
  gfc8NONAME[i]=binGFC8-0.5;
  qc2NONAME[i]=binQC2-0.5;
  qc4NONAME[i]=binQC4-0.5;
  qc6NONAME[i]=binQC6-0.5;
  qc8NONAME[i]=binQC8-0.5;
  fqdNONAME[i]=binFQD-0.5;
  lyz1NONAME[i]=binLYZ1-0.5;
  lyzepNONAME[i]=binLYZEP-0.5;
 }
 
 // MCEP:
 TGraph *mcepMeanNONAME = new TGraph(1);    
 mcepMeanNONAME->SetPoint(0,binMCEP-0.5,mergedValueMCEP);  
 mcepMeanNONAME->SetMarkerStyle(25);
 mcepMeanNONAME->SetMarkerColor(kBlack); 
 mcepMeanNONAME->SetMarkerSize(1.25); 
 
 TGraph* mcepTGraphNONAME = new TGraph(nPoints, mcepNONAME, mcepValueNONAME);
 mcepTGraphNONAME->SetMarkerStyle(21);
 mcepTGraphNONAME->SetMarkerColor(kBlack); 
 mcepTGraphNONAME->SetMarkerSize(0.75); 
 
 TGraph *mcepBoxNONAME = new TGraph(5);
 mcepBoxNONAME->SetPoint(0,(binMCEP-0.5)-boxWidth,mcepMinValueNONAME);
 mcepBoxNONAME->SetPoint(1,(binMCEP-0.5)+boxWidth,mcepMinValueNONAME);
 mcepBoxNONAME->SetPoint(2,(binMCEP-0.5)+boxWidth,mcepMaxValueNONAME);
 mcepBoxNONAME->SetPoint(3,(binMCEP-0.5)-boxWidth,mcepMaxValueNONAME);
 mcepBoxNONAME->SetPoint(4,(binMCEP-0.5)-boxWidth,mcepMinValueNONAME);    
 mcepBoxNONAME->SetFillStyle(1001);
 mcepBoxNONAME->SetFillColor(kGray);
 
 // SP:
 TGraph *spMeanNONAME = new TGraph(1);    
 spMeanNONAME->SetPoint(0,binSP-0.5,mergedValueSP);  
 spMeanNONAME->SetMarkerStyle(25);
 spMeanNONAME->SetMarkerColor(kViolet+3); 
 spMeanNONAME->SetMarkerSize(1.25); 
 
 TGraph* spTGraphNONAME = new TGraph(nPoints, spNONAME, spValueNONAME);
 spTGraphNONAME->SetMarkerStyle(21);
 spTGraphNONAME->SetMarkerColor(kViolet+3); 
 spTGraphNONAME->SetMarkerSize(0.75); 
 
 TGraph *spBoxNONAME = new TGraph(5);
 spBoxNONAME->SetPoint(0,(binSP-0.5)-boxWidth,spMinValueNONAME);
 spBoxNONAME->SetPoint(1,(binSP-0.5)+boxWidth,spMinValueNONAME);
 spBoxNONAME->SetPoint(2,(binSP-0.5)+boxWidth,spMaxValueNONAME);
 spBoxNONAME->SetPoint(3,(binSP-0.5)-boxWidth,spMaxValueNONAME);
 spBoxNONAME->SetPoint(4,(binSP-0.5)-boxWidth,spMinValueNONAME);    
 spBoxNONAME->SetFillStyle(1001);
 spBoxNONAME->SetFillColor(kViolet-9);
 
 // GFC{2}:
 TGraph *gfc2MeanNONAME = new TGraph(1);    
 gfc2MeanNONAME->SetPoint(0,binGFC2-0.5,mergedValueGFC2);  
 gfc2MeanNONAME->SetMarkerStyle(25);
 gfc2MeanNONAME->SetMarkerColor(kBlue); 
 gfc2MeanNONAME->SetMarkerSize(1.25); 
 
 TGraph* gfc2TGraphNONAME = new TGraph(nPoints, gfc2NONAME, gfc2ValueNONAME);
 gfc2TGraphNONAME->SetMarkerStyle(21);
 gfc2TGraphNONAME->SetMarkerColor(kBlue); 
 gfc2TGraphNONAME->SetMarkerSize(0.75); 
 
 TGraph *gfc2BoxNONAME = new TGraph(5);
 gfc2BoxNONAME->SetPoint(0,(binGFC2-0.5)-boxWidth,gfc2MinValueNONAME);
 gfc2BoxNONAME->SetPoint(1,(binGFC2-0.5)+boxWidth,gfc2MinValueNONAME);
 gfc2BoxNONAME->SetPoint(2,(binGFC2-0.5)+boxWidth,gfc2MaxValueNONAME);
 gfc2BoxNONAME->SetPoint(3,(binGFC2-0.5)-boxWidth,gfc2MaxValueNONAME);
 gfc2BoxNONAME->SetPoint(4,(binGFC2-0.5)-boxWidth,gfc2MinValueNONAME);    
 gfc2BoxNONAME->SetFillStyle(1001);
 gfc2BoxNONAME->SetFillColor(kBlue-10);
 
 // GFC{4}:
 TGraph *gfc4MeanNONAME = new TGraph(1);    
 gfc4MeanNONAME->SetPoint(0,binGFC4-0.5,mergedValueGFC4);  
 gfc4MeanNONAME->SetMarkerStyle(25);
 gfc4MeanNONAME->SetMarkerColor(kBlue); 
 gfc4MeanNONAME->SetMarkerSize(1.25); 
 
 TGraph* gfc4TGraphNONAME = new TGraph(nPoints, gfc4NONAME, gfc4ValueNONAME);
 gfc4TGraphNONAME->SetMarkerStyle(21);
 gfc4TGraphNONAME->SetMarkerColor(kBlue); 
 gfc4TGraphNONAME->SetMarkerSize(0.75); 
 
 TGraph *gfc4BoxNONAME = new TGraph(5);
 gfc4BoxNONAME->SetPoint(0,(binGFC4-0.5)-boxWidth,gfc4MinValueNONAME);
 gfc4BoxNONAME->SetPoint(1,(binGFC4-0.5)+boxWidth,gfc4MinValueNONAME);
 gfc4BoxNONAME->SetPoint(2,(binGFC4-0.5)+boxWidth,gfc4MaxValueNONAME);
 gfc4BoxNONAME->SetPoint(3,(binGFC4-0.5)-boxWidth,gfc4MaxValueNONAME);
 gfc4BoxNONAME->SetPoint(4,(binGFC4-0.5)-boxWidth,gfc4MinValueNONAME);    
 gfc4BoxNONAME->SetFillStyle(1001);
 gfc4BoxNONAME->SetFillColor(kBlue-10);
 
 // GFC{6}:
 TGraph *gfc6MeanNONAME = new TGraph(1);    
 gfc6MeanNONAME->SetPoint(0,binGFC6-0.5,mergedValueGFC6);  
 gfc6MeanNONAME->SetMarkerStyle(25);
 gfc6MeanNONAME->SetMarkerColor(kBlue); 
 gfc6MeanNONAME->SetMarkerSize(1.25); 
 
 TGraph* gfc6TGraphNONAME = new TGraph(nPoints, gfc6NONAME, gfc6ValueNONAME);
 gfc6TGraphNONAME->SetMarkerStyle(21);
 gfc6TGraphNONAME->SetMarkerColor(kBlue); 
 gfc6TGraphNONAME->SetMarkerSize(0.75); 
 
 TGraph *gfc6BoxNONAME = new TGraph(5);
 gfc6BoxNONAME->SetPoint(0,(binGFC6-0.5)-boxWidth,gfc6MinValueNONAME);
 gfc6BoxNONAME->SetPoint(1,(binGFC6-0.5)+boxWidth,gfc6MinValueNONAME);
 gfc6BoxNONAME->SetPoint(2,(binGFC6-0.5)+boxWidth,gfc6MaxValueNONAME);
 gfc6BoxNONAME->SetPoint(3,(binGFC6-0.5)-boxWidth,gfc6MaxValueNONAME);
 gfc6BoxNONAME->SetPoint(4,(binGFC6-0.5)-boxWidth,gfc6MinValueNONAME);    
 gfc6BoxNONAME->SetFillStyle(1001);
 gfc6BoxNONAME->SetFillColor(kBlue-10);
 
 // GFC{8}:
 TGraph *gfc8MeanNONAME = new TGraph(1);    
 gfc8MeanNONAME->SetPoint(0,binGFC8-0.5,mergedValueGFC8);  
 gfc8MeanNONAME->SetMarkerStyle(25);
 gfc8MeanNONAME->SetMarkerColor(kBlue); 
 gfc8MeanNONAME->SetMarkerSize(1.25); 
 
 TGraph* gfc8TGraphNONAME = new TGraph(nPoints, gfc8NONAME, gfc8ValueNONAME);
 gfc8TGraphNONAME->SetMarkerStyle(21);
 gfc8TGraphNONAME->SetMarkerColor(kBlue); 
 gfc8TGraphNONAME->SetMarkerSize(0.75); 
 
 TGraph *gfc8BoxNONAME = new TGraph(5);
 gfc8BoxNONAME->SetPoint(0,(binGFC8-0.5)-boxWidth,gfc8MinValueNONAME);
 gfc8BoxNONAME->SetPoint(1,(binGFC8-0.5)+boxWidth,gfc8MinValueNONAME);
 gfc8BoxNONAME->SetPoint(2,(binGFC8-0.5)+boxWidth,gfc8MaxValueNONAME);
 gfc8BoxNONAME->SetPoint(3,(binGFC8-0.5)-boxWidth,gfc8MaxValueNONAME);
 gfc8BoxNONAME->SetPoint(4,(binGFC8-0.5)-boxWidth,gfc8MinValueNONAME);    
 gfc8BoxNONAME->SetFillStyle(1001);
 gfc8BoxNONAME->SetFillColor(kBlue-10);
 
 // QC{2}:
 TGraph *qc2MeanNONAME = new TGraph(1);    
 qc2MeanNONAME->SetPoint(0,binQC2-0.5,mergedValueQC2);  
 qc2MeanNONAME->SetMarkerStyle(25);
 qc2MeanNONAME->SetMarkerColor(kRed); 
 qc2MeanNONAME->SetMarkerSize(1.25); 
 
 TGraph* qc2TGraphNONAME = new TGraph(nPoints, qc2NONAME, qc2ValueNONAME);
 qc2TGraphNONAME->SetMarkerStyle(21);
 qc2TGraphNONAME->SetMarkerColor(kRed); 
 qc2TGraphNONAME->SetMarkerSize(0.75); 
 
 TGraph *qc2BoxNONAME = new TGraph(5);
 qc2BoxNONAME->SetPoint(0,(binQC2-0.5)-boxWidth,qc2MinValueNONAME);
 qc2BoxNONAME->SetPoint(1,(binQC2-0.5)+boxWidth,qc2MinValueNONAME);
 qc2BoxNONAME->SetPoint(2,(binQC2-0.5)+boxWidth,qc2MaxValueNONAME);
 qc2BoxNONAME->SetPoint(3,(binQC2-0.5)-boxWidth,qc2MaxValueNONAME);
 qc2BoxNONAME->SetPoint(4,(binQC2-0.5)-boxWidth,qc2MinValueNONAME);    
 qc2BoxNONAME->SetFillStyle(1001);
 qc2BoxNONAME->SetFillColor(kRed-10);
 
 // QC{4}:
 TGraph *qc4MeanNONAME = new TGraph(1);    
 qc4MeanNONAME->SetPoint(0,binQC4-0.5,mergedValueQC4);  
 qc4MeanNONAME->SetMarkerStyle(25);
 qc4MeanNONAME->SetMarkerColor(kRed); 
 qc4MeanNONAME->SetMarkerSize(1.25); 
 
 TGraph* qc4TGraphNONAME = new TGraph(nPoints, qc4NONAME, qc4ValueNONAME);
 qc4TGraphNONAME->SetMarkerStyle(21);
 qc4TGraphNONAME->SetMarkerColor(kRed); 
 qc4TGraphNONAME->SetMarkerSize(0.75); 
 
 TGraph *qc4BoxNONAME = new TGraph(5);
 qc4BoxNONAME->SetPoint(0,(binQC4-0.5)-boxWidth,qc4MinValueNONAME);
 qc4BoxNONAME->SetPoint(1,(binQC4-0.5)+boxWidth,qc4MinValueNONAME);
 qc4BoxNONAME->SetPoint(2,(binQC4-0.5)+boxWidth,qc4MaxValueNONAME);
 qc4BoxNONAME->SetPoint(3,(binQC4-0.5)-boxWidth,qc4MaxValueNONAME);
 qc4BoxNONAME->SetPoint(4,(binQC4-0.5)-boxWidth,qc4MinValueNONAME);    
 qc4BoxNONAME->SetFillStyle(1001);
 qc4BoxNONAME->SetFillColor(kRed-10);
 
 // QC{6}:
 TGraph *qc6MeanNONAME = new TGraph(1);    
 qc6MeanNONAME->SetPoint(0,binQC6-0.5,mergedValueQC6);  
 qc6MeanNONAME->SetMarkerStyle(25);
 qc6MeanNONAME->SetMarkerColor(kRed); 
 qc6MeanNONAME->SetMarkerSize(1.25); 
 
 TGraph* qc6TGraphNONAME = new TGraph(nPoints, qc6NONAME, qc6ValueNONAME);
 qc6TGraphNONAME->SetMarkerStyle(21);
 qc6TGraphNONAME->SetMarkerColor(kRed); 
 qc6TGraphNONAME->SetMarkerSize(0.75); 
 
 TGraph *qc6BoxNONAME = new TGraph(5);
 qc6BoxNONAME->SetPoint(0,(binQC6-0.5)-boxWidth,qc6MinValueNONAME);
 qc6BoxNONAME->SetPoint(1,(binQC6-0.5)+boxWidth,qc6MinValueNONAME);
 qc6BoxNONAME->SetPoint(2,(binQC6-0.5)+boxWidth,qc6MaxValueNONAME);
 qc6BoxNONAME->SetPoint(3,(binQC6-0.5)-boxWidth,qc6MaxValueNONAME);
 qc6BoxNONAME->SetPoint(4,(binQC6-0.5)-boxWidth,qc6MinValueNONAME);    
 qc6BoxNONAME->SetFillStyle(1001);
 qc6BoxNONAME->SetFillColor(kRed-10);
 
 // QC{8}:
 TGraph *qc8MeanNONAME = new TGraph(1);    
 qc8MeanNONAME->SetPoint(0,binQC8-0.5,mergedValueQC8);  
 qc8MeanNONAME->SetMarkerStyle(25);
 qc8MeanNONAME->SetMarkerColor(kRed); 
 qc8MeanNONAME->SetMarkerSize(1.25); 
 
 TGraph* qc8TGraphNONAME = new TGraph(nPoints, qc8NONAME, qc8ValueNONAME);
 qc8TGraphNONAME->SetMarkerStyle(21);
 qc8TGraphNONAME->SetMarkerColor(kRed); 
 qc8TGraphNONAME->SetMarkerSize(0.75); 
 
 TGraph *qc8BoxNONAME = new TGraph(5);
 qc8BoxNONAME->SetPoint(0,(binQC8-0.5)-boxWidth,qc8MinValueNONAME);
 qc8BoxNONAME->SetPoint(1,(binQC8-0.5)+boxWidth,qc8MinValueNONAME);
 qc8BoxNONAME->SetPoint(2,(binQC8-0.5)+boxWidth,qc8MaxValueNONAME);
 qc8BoxNONAME->SetPoint(3,(binQC8-0.5)-boxWidth,qc8MaxValueNONAME);
 qc8BoxNONAME->SetPoint(4,(binQC8-0.5)-boxWidth,qc8MinValueNONAME);    
 qc8BoxNONAME->SetFillStyle(1001);
 qc8BoxNONAME->SetFillColor(kRed-10);
 
 // FQD:
 TGraph *fqdMeanNONAME = new TGraph(1);    
 fqdMeanNONAME->SetPoint(0,binFQD-0.5,mergedValueFQD);  
 fqdMeanNONAME->SetMarkerStyle(25);
 fqdMeanNONAME->SetMarkerColor(kViolet+3); 
 fqdMeanNONAME->SetMarkerSize(1.25); 
 
 TGraph* fqdTGraphNONAME = new TGraph(nPoints, fqdNONAME, fqdValueNONAME);
 fqdTGraphNONAME->SetMarkerStyle(21);
 fqdTGraphNONAME->SetMarkerColor(kOrange+7); 
 fqdTGraphNONAME->SetMarkerSize(0.75); 
 
 TGraph *fqdBoxNONAME = new TGraph(5);
 fqdBoxNONAME->SetPoint(0,(binFQD-0.5)-boxWidth,fqdMinValueNONAME);
 fqdBoxNONAME->SetPoint(1,(binFQD-0.5)+boxWidth,fqdMinValueNONAME);
 fqdBoxNONAME->SetPoint(2,(binFQD-0.5)+boxWidth,fqdMaxValueNONAME);
 fqdBoxNONAME->SetPoint(3,(binFQD-0.5)-boxWidth,fqdMaxValueNONAME);
 fqdBoxNONAME->SetPoint(4,(binFQD-0.5)-boxWidth,fqdMinValueNONAME);    
 fqdBoxNONAME->SetFillStyle(1001);
 fqdBoxNONAME->SetFillColor(kOrange-9);
 
 // ... LYZ* to be added ...
   
 TCanvas* intFlowCanvasNONAME = new TCanvas("Integrated Flow NONAME","Integrated Flow NONAME",1000,600);
 
 intFlowCanvasNONAME->Divide(2,1);
 
 // 1st pad is for plot:
 (intFlowCanvasNONAME->cd(1))->SetPad(0.0,0.0,0.75,1.0);
 
 
 if(intFlowNONAME) 
 {
  TString intFlowNameNONAME("Superimposing ");
  intFlowNameNONAME+=counter;
  intFlowNameNONAME+=" independent runs";
  intFlowNONAME->SetTitle(intFlowNameNONAME.Data());
  (intFlowNONAME->GetYaxis())->SetRangeUser(0,0.20);
  intFlowNONAME->Draw();
 }
 
 // MCEP: 
 if(mcepBoxNONAME && mcepMinValueNONAME < 1000.) mcepBoxNONAME->Draw("LFSAME");
 if(mcepTGraphNONAME) mcepTGraphNONAME->Draw("PSAME");
 if(mcepMeanNONAME) mcepMeanNONAME->Draw("PSAME");
 
 // SP: 
 if(spBoxNONAME && spMinValueNONAME < 1000.) spBoxNONAME->Draw("LFSAME");
 if(spTGraphNONAME) spTGraphNONAME->Draw("PSAME");
 if(spMeanNONAME) spMeanNONAME->Draw("PSAME");
 
 // GFC{2}: 
 if(gfc2BoxNONAME && gfc2MinValueNONAME < 1000.) gfc2BoxNONAME->Draw("LFSAME");
 if(gfc2TGraphNONAME) gfc2TGraphNONAME->Draw("PSAME");
 if(gfc2MeanNONAME) gfc2MeanNONAME->Draw("PSAME");
 
 // GFC{4}: 
 if(gfc4BoxNONAME && gfc4MinValueNONAME < 1000.) gfc4BoxNONAME->Draw("LFSAME");
 if(gfc4TGraphNONAME) gfc4TGraphNONAME->Draw("PSAME");
 if(gfc4MeanNONAME) gfc4MeanNONAME->Draw("PSAME");
 
 // GFC{6}: 
 if(gfc6BoxNONAME && gfc6MinValueNONAME < 1000.) gfc6BoxNONAME->Draw("LFSAME");
 if(gfc6TGraphNONAME) gfc6TGraphNONAME->Draw("PSAME");
 if(gfc6MeanNONAME) gfc6MeanNONAME->Draw("PSAME");
 
 // GFC{8}: 
 if(gfc8BoxNONAME && gfc8MinValueNONAME < 1000.) gfc8BoxNONAME->Draw("LFSAME");
 if(gfc8TGraphNONAME) gfc8TGraphNONAME->Draw("PSAME");
 if(gfc8MeanNONAME) gfc8MeanNONAME->Draw("PSAME");
 
 // QC{2}: 
 if(qc2BoxNONAME && qc2MinValueNONAME < 1000.) qc2BoxNONAME->Draw("LFSAME");
 if(qc2TGraphNONAME) qc2TGraphNONAME->Draw("PSAME");
 if(qc2MeanNONAME) qc2MeanNONAME->Draw("PSAME");
 
 // QC{4}: 
 if(qc4BoxNONAME && qc4MinValueNONAME < 1000.) qc4BoxNONAME->Draw("LFSAME");
 if(qc4TGraphNONAME) qc4TGraphNONAME->Draw("PSAME");
 if(qc4MeanNONAME) qc4MeanNONAME->Draw("PSAME");
 
 // QC{6}: 
 if(qc6BoxNONAME && qc6MinValueNONAME < 1000.) qc6BoxNONAME->Draw("LFSAME");
 if(qc6TGraphNONAME) qc6TGraphNONAME->Draw("PSAME");
 if(qc6MeanNONAME) qc6MeanNONAME->Draw("PSAME");
 
 // QC{8}: 
 if(qc8BoxNONAME && qc8MinValueNONAME < 1000.) qc8BoxNONAME->Draw("LFSAME");
 if(qc8TGraphNONAME) qc8TGraphNONAME->Draw("PSAME");
 if(qc8MeanNONAME) qc8MeanNONAME->Draw("PSAME");
 
 // FQD: 
 if(fqdBoxNONAME && fqdMinValueNONAME < 1000.) fqdBoxNONAME->Draw("LFSAME");
 if(fqdTGraphNONAME) fqdTGraphNONAME->Draw("PSAME");
 if(fqdMeanNONAME) fqdMeanNONAME->Draw("PSAME");

 // ... LYZ* to be added ...
 
 // 2nd pad is for legend:   
 (intFlowCanvasNONAME->cd(2))->SetPad(0.75,0.0,1.0,1.0);

 // count real estimates:
 Int_t mcepCountRealNONAME = 0;
 Int_t spCountRealNONAME = 0;
 Int_t gfc2CountRealNONAME = 0;
 Int_t gfc4CountRealNONAME = 0;
 Int_t gfc6CountRealNONAME = 0;
 Int_t gfc8CountRealNONAME = 0;
 Int_t qc2CountRealNONAME = 0;
 Int_t qc4CountRealNONAME = 0;
 Int_t qc6CountRealNONAME = 0;
 Int_t qc8CountRealNONAME = 0;
 Int_t fqdCountRealNONAME = 0;
 // ... LYZ* to be added ...
 
 for(Int_t i=0;i<nEstimates;i++)
 {
  if(mcepValueNONAME[i]>0.) mcepCountRealNONAME++; 
  if(spValueNONAME[i]>0.) spCountRealNONAME++; 
  if(gfc2ValueNONAME[i]>0.) gfc2CountRealNONAME++; 
  if(gfc4ValueNONAME[i]>0.) gfc4CountRealNONAME++; 
  if(gfc6ValueNONAME[i]>0.) gfc6CountRealNONAME++; 
  if(gfc8ValueNONAME[i]>0.) gfc8CountRealNONAME++; 
  if(qc2ValueNONAME[i]>0.) qc2CountRealNONAME++; 
  if(qc4ValueNONAME[i]>0.) qc4CountRealNONAME++; 
  if(qc6ValueNONAME[i]>0.) qc6CountRealNONAME++; 
  if(qc8ValueNONAME[i]>0.) qc8CountRealNONAME++;
  if(fqdValueNONAME[i]>0.) fqdCountRealNONAME++; 
  // ... LYZ* to be added ... 
 }

 TPaveText *textDefaultNONAME = new TPaveText(0.05,0.67,0.95,0.90,"NDC");
 textDefaultNONAME->SetTextFont(72);
 textDefaultNONAME->SetTextSize(0.08);
  
 TString *entryDefaultRealNONAME = new TString("Real estimates");
 //TString *entryDefaultOutOfNONAME = new TString("out of"); 
 TString *entryDefaultTotalNumberNONAME = new TString(" out of total number of ");
 TString *entryDefaultTotalIndNONAME = new TString("independent");
 TString *entryDefaultTotalSimNONAME = new TString("simulations:");
  
 textDefaultNONAME->AddText(entryDefaultRealNONAME->Data());
 //textDefaultNONAME->AddText(entryDefaultOutOfNONAME->Data());
 textDefaultNONAME->AddText(entryDefaultTotalNumberNONAME->Data());
 textDefaultNONAME->AddText(entryDefaultTotalIndNONAME->Data());
 textDefaultNONAME->AddText(entryDefaultTotalSimNONAME->Data());
 
 // results:
 TPaveText *textResultsNONAME = new TPaveText(0.05,0.12,0.95,0.60,"NDC");
 textResultsNONAME->SetTextFont(72);
 textResultsNONAME->SetTextSize(0.06);
 
 // entries:
 TString *entryIntFlowMCNONAME    = new TString("MC ................ ");
 TString *entryIntFlowSPNONAME    = new TString("SP ................ ");  
 TString *entryIntFlowGFC2NONAME  = new TString("GFC{2} ........ ");
 TString *entryIntFlowGFC4NONAME  = new TString("GFC{4} ........ ");
 TString *entryIntFlowGFC6NONAME  = new TString("GFC{6} ........ ");
 TString *entryIntFlowGFC8NONAME  = new TString("GFC{8} ........ "); 
 TString *entryIntFlowQC2NONAME   = new TString("QC{2} .......... ");
 TString *entryIntFlowQC4NONAME   = new TString("QC{4} .......... ");
 TString *entryIntFlowQC6NONAME   = new TString("QC{6} .......... ");
 TString *entryIntFlowQC8NONAME   = new TString("QC{8} .......... ");
 TString *entryIntFlowFQDNONAME   = new TString("FQD ............. ");
 TString *entryIntFlowLYZ1NONAME  = new TString("LYZ ............. ");
 TString *entryIntFlowLYZEPNONAME = new TString("LYZEP ........ ");

 if(entryIntFlowMCNONAME)
 { 
  (*entryIntFlowMCNONAME)+=(Long_t)mcepCountRealNONAME;
  entryIntFlowMCNONAME->Append(" out of ");
  (*entryIntFlowMCNONAME)+=(Long_t)counter;
  if(textResultsNONAME)textResultsNONAME->AddText(entryIntFlowMCNONAME->Data());
 }
 
 if(entryIntFlowSPNONAME)
 { 
  (*entryIntFlowSPNONAME)+=(Long_t)spCountRealNONAME;
  entryIntFlowSPNONAME->Append(" out of ");
  (*entryIntFlowSPNONAME)+=(Long_t)counter;
  if(textResultsNONAME)textResultsNONAME->AddText(entryIntFlowSPNONAME->Data());
 }
 
 if(entryIntFlowGFC2NONAME)
 { 
  (*entryIntFlowGFC2NONAME)+=(Long_t)gfc2CountRealNONAME;
  entryIntFlowGFC2NONAME->Append(" out of ");
  (*entryIntFlowGFC2NONAME)+=(Long_t)counter;
  if(textResultsNONAME)textResultsNONAME->AddText(entryIntFlowGFC2NONAME->Data());
 }
 
 if(entryIntFlowGFC4NONAME)
 { 
  (*entryIntFlowGFC4NONAME)+=(Long_t)gfc4CountRealNONAME;
  entryIntFlowGFC4NONAME->Append(" out of ");
  (*entryIntFlowGFC4NONAME)+=(Long_t)counter;
  if(textResultsNONAME)textResultsNONAME->AddText(entryIntFlowGFC4NONAME->Data());
 }

 if(entryIntFlowGFC6NONAME)
 { 
  (*entryIntFlowGFC6NONAME)+=(Long_t)gfc6CountRealNONAME;
  entryIntFlowGFC6NONAME->Append(" out of ");
  (*entryIntFlowGFC6NONAME)+=(Long_t)counter;
  if(textResultsNONAME)textResultsNONAME->AddText(entryIntFlowGFC6NONAME->Data());
 }

 if(entryIntFlowGFC8NONAME)
 { 
  (*entryIntFlowGFC8NONAME)+=(Long_t)gfc8CountRealNONAME;
  entryIntFlowGFC8NONAME->Append(" out of ");
  (*entryIntFlowGFC8NONAME)+=(Long_t)counter;
  if(textResultsNONAME)textResultsNONAME->AddText(entryIntFlowGFC8NONAME->Data());
 }
 
 if(entryIntFlowQC2NONAME)
 { 
  (*entryIntFlowQC2NONAME)+=(Long_t)qc2CountRealNONAME;
  entryIntFlowQC2NONAME->Append(" out of ");
  (*entryIntFlowQC2NONAME)+=(Long_t)counter;
  if(textResultsNONAME)textResultsNONAME->AddText(entryIntFlowQC2NONAME->Data());
 }
 
 if(entryIntFlowQC4NONAME)
 { 
  (*entryIntFlowQC4NONAME)+=(Long_t)qc4CountRealNONAME;
  entryIntFlowQC4NONAME->Append(" out of ");
  (*entryIntFlowQC4NONAME)+=(Long_t)counter;
  if(textResultsNONAME)textResultsNONAME->AddText(entryIntFlowQC4NONAME->Data());
 }

 if(entryIntFlowQC2NONAME)
 { 
  (*entryIntFlowQC6NONAME)+=(Long_t)qc6CountRealNONAME;
  entryIntFlowQC6NONAME->Append(" out of ");
  (*entryIntFlowQC6NONAME)+=(Long_t)counter;
  if(textResultsNONAME)textResultsNONAME->AddText(entryIntFlowQC6NONAME->Data());
 }

 if(entryIntFlowQC8NONAME)
 { 
  (*entryIntFlowQC8NONAME)+=(Long_t)qc8CountRealNONAME;
  entryIntFlowQC8NONAME->Append(" out of ");
  (*entryIntFlowQC8NONAME)+=(Long_t)counter;
  if(textResultsNONAME)textResultsNONAME->AddText(entryIntFlowQC8NONAME->Data());
 }
 
 if(entryIntFlowFQDNONAME)
 { 
  (*entryIntFlowFQDNONAME)+=(Long_t)fqdCountRealNONAME;
  entryIntFlowFQDNONAME->Append(" out of ");
  (*entryIntFlowFQDNONAME)+=(Long_t)counter;
  if(textResultsNONAME)textResultsNONAME->AddText(entryIntFlowFQDNONAME->Data());
 }
 
 /*
 if(entryIntFlowLYZ1NONAME)
 { 
  (*entryIntFlowLYZ1NONAME)+=(Long_t)lyz1CountRealNONAME;
  entryIntFlowLYZ1NONAME->Append(" out of ");
  (*entryIntFlowLYZ1NONAME)+=(Long_t)counter;
  if(textResultsNONAME)textResultsNONAME->AddText(entryIntFlowLYZ1NONAME->Data());
 }
 
 if(entryIntFlowLYZEPNONAME)
 { 
  (*entryIntFlowLYZEPNONAME)+=(Long_t)lyzepCountRealNONAME;
  entryIntFlowLYZEPNONAME->Append(" out of ");
  (*entryIntFlowLYZEPNONAME)+=(Long_t)counter;
  if(textResultsNONAME)textResultsNONAME->AddText(entryIntFlowLYZEPNONAME->Data());
 }
 */
 
 if(textDefaultNONAME) textDefaultNONAME->Draw();
 if(textResultsNONAME) textResultsNONAME->Draw();
 
}

void LoadSpreadLibraries(const libModes mode) {
  
  //--------------------------------------
  // Load the needed libraries most of them already loaded by aliroot
  //--------------------------------------
  gSystem->Load("libTree.so");
  gSystem->Load("libGeom.so");
  gSystem->Load("libVMC.so");
  gSystem->Load("libXMLIO.so");
  gSystem->Load("libPhysics.so");
  
  //----------------------------------------------------------
  // >>>>>>>>>>> Local mode <<<<<<<<<<<<<< 
  //----------------------------------------------------------
  if (mode==mLocal) {
    //--------------------------------------------------------
    // If you want to use already compiled libraries 
    // in the aliroot distribution
    //--------------------------------------------------------

  //==================================================================================  
  //load needed libraries:
  gSystem->AddIncludePath("-I$ROOTSYS/include");
  gSystem->Load("libTree.so");

  // for AliRoot
  gSystem->AddIncludePath("-I$ALICE_ROOT/include");
  gSystem->Load("libANALYSIS.so");
  gSystem->Load("libPWG2flowCommon.so");
  cerr<<"libPWG2flowCommon.so loaded ..."<<endl;
  
  }
  
  else if (mode==mLocalSource) {
 
    // In root inline compile
  
    // Constants  
    gROOT->LoadMacro("AliFlowCommon/AliFlowCommonConstants.cxx+");
    gROOT->LoadMacro("AliFlowCommon/AliFlowLYZConstants.cxx+");
    gROOT->LoadMacro("AliFlowCommon/AliFlowCumuConstants.cxx+");
    
    // Flow event
    gROOT->LoadMacro("AliFlowCommon/AliFlowVector.cxx+"); 
    gROOT->LoadMacro("AliFlowCommon/AliFlowTrackSimple.cxx+");    
    gROOT->LoadMacro("AliFlowCommon/AliFlowEventSimple.cxx+");
    
    // Cuts
    gROOT->LoadMacro("AliFlowCommon/AliFlowTrackSimpleCuts.cxx+");    
    
    // Output histosgrams
    gROOT->LoadMacro("AliFlowCommon/AliFlowCommonHist.cxx+");
    gROOT->LoadMacro("AliFlowCommon/AliFlowCommonHistResults.cxx+");
    gROOT->LoadMacro("AliFlowCommon/AliFlowLYZHist1.cxx+");
    gROOT->LoadMacro("AliFlowCommon/AliFlowLYZHist2.cxx+");
       
    cout << "finished loading macros!" << endl;  
    
  }  
  
  
}