[u/mrichter/AliRoot.git] / PWGCF / FLOW / macros / fqd.C

//==========================================================================//
// Macro to refit the q-distribution in FQD method. By using this macro the // 
// fitting parameters for q-distribution can be tuned and correspondingly   //
// the estimate for reference flow harmonic obtained from FQD method can be //
// significantly improved.                                                  // 
//==========================================================================//

// Fitting parameters for q-distribution:
Double_t vStart = 0.05; // starting value for v fit
Double_t vMin = 0.0; // lower bound for v fit 
Double_t vMax = 0.25; // upper bound for v fit  
Double_t sigma2Start = 0.75; // starting value for sigma^2 fit 
Double_t sigma2Min = 0.5; // lower bound for sigma^2 fit (according to theorists must be >= 0.5)   
Double_t sigma2Max = 2.5; // upper bound for sigma^2 fit
Double_t treshold = 5.; // first and last bin taken for fitting are determined as the first and last bin with more than 'treshold' number of entries
Bool_t finalResultIsFromSigma2Fitted = kTRUE; // final saved result is obtained with sigma^2 fitted (kTRUE) or sigma^2 fixed (kFALSE)
Bool_t printOnTheScreen = kTRUE; // print or not the final results on the screen
Bool_t plotResults = kTRUE; // plot on the screen q-distribution and resulting fitting functions (for non-fixed sigma^2 and fixed sigma^2)

// Name of the common output file:
TString outputFileName = "AnalysisResults.root";
//TString outputFileName = "outputCentrality0.root";

enum libModes {mLocal,mLocalSource};

TFile *commonOutputFile = NULL; // common output file "AnalysisResults.root"
TDirectoryFile *dirFileFQD = NULL; // FQD's TDirectoryFile in "AnalysisResults.root"
TList *outputListFQD = NULL; // output list holding all FQD objects

void fqd(TString analysisType="ESD", Int_t analysisMode=mLocal)
{
 // 1. analysisType: type of analysis can be ESD, AOD, MC, MK, ....
 //                  (if type="" output files are from simulation 'on the fly')
 // 2. analysisMode: if analysisMode = mLocal -> analyze data on your computer using aliroot
 //                  if analysisMode = mLocalSource -> analyze data on your computer using root + source files 
  
 // Cross-check if the user's settings make sense:
 CrossCheckSettings();
 // Load needed libraries:                       
 LoadLibrariesFQD(analysisMode); 
 // Get output list which holds all FQD objects:
 GetOutputList(analysisType); 
 // Redo fit of q-distribution:
 RedoFit();
 // Plot q-distribution and resulting fitting functions:
 if(plotResults) Plot(); 
 // Save the new results in the common output file: 
 Save();
  
} // end of void fqd(TString type="ESD", Int_t mode=mLocal)   

// =========================================================================================== 

void RedoFit()
{
 // Redo the fit of q-distribution.

 AliFlowAnalysisWithFittingQDistribution* fqd = new AliFlowAnalysisWithFittingQDistribution();
 fqd->GetOutputHistograms(outputListFQD);
 // Set new fitting parameters:
 fqd->SetvStart(vStart);  
 fqd->SetvMin(vMin);
 fqd->SetvMax(vMax);
 fqd->SetSigma2Start(sigma2Start);  
 fqd->SetSigma2Min(sigma2Min);
 fqd->SetSigma2Max(sigma2Max);
 fqd->SetTreshold(treshold);
 fqd->SetFinalResultIsFromSigma2Fitted(finalResultIsFromSigma2Fitted);
 fqd->SetPrintOnTheScreen(printOnTheScreen); 
 // Save new fitting parameters:  
 fqd->StoreFittingParameters();
 // Redo fit:
 fqd->Finish(kTRUE);
 
} // end of void RedoFit(TList *outputList)
 
// =========================================================================================== 
 
void GetOutputList(TString analysisType)
{
 // Get output list which holds all FQD objects.
  
 // Access common output file:
 commonOutputFile = AccessOutputFile(outputFileName);
 
 // Access from common output file the TDirectoryFile for FQD method
 // and from it the output list holding all objects:
 GetListWithHistograms(commonOutputFile,analysisType,"FQD");
   
} // end of TList* GetOutputList(TString analysisType)
 
// ===========================================================================================

void Plot()
{
 // Plot q-distribution and resulting fitting functions.
 
 gROOT->SetStyle("Plain"); // default color is white instead of gray
 gStyle->SetOptStat(0); // remove stat. box from all histos
 legend = new TLegend(0.6,0.55,0.85,0.7); 
 legend->SetFillStyle(0);
 // q-distribution:
 TH1D *qDistribution = dynamic_cast<TH1D*>(outputListFQD->FindObject("fqDistribution"));
 Cosmetics(qDistribution);
 legend->AddEntry(qDistribution,"q-distribution","f");
 qDistribution->Draw();
 // resulting fitting functions:
 TF1 *fittingFun[2] = {NULL};
 TString sigmaFlag[2] = {"#sigma^{2} not fitted","#sigma^{2} fitted"};
 Double_t lineColors[2] = {kBlue,kRed};
 for(Int_t s2F=0;s2F<2;s2F++)
 {
  fittingFun[s2F] = dynamic_cast<TF1*>(outputListFQD->FindObject(Form("fFittingFunction, %s",sigmaFlag[s2F].Data())));
  fittingFun[s2F]->SetLineColor(lineColors[s2F]);
  fittingFun[s2F]->Draw("SAME");
  legend->AddEntry(fittingFun[s2F]->GetName(),sigmaFlag[s2F].Data(),"l");
 }
 legend->Draw("SAME"); 
  
} // end of Plot()

// ===========================================================================================

void Save()
{
 // Save the new results of fitting for FQD method in the common output file.
 
 dirFileFQD->Add(outputListFQD,kTRUE);
 dirFileFQD->Write(dirFileFQD->GetName(),TObject::kSingleKey+TObject::kOverwrite);
 //delete commonOutputFile;

} // end of void Save()

// ===========================================================================================

void Cosmetics(TH1D *hist)
{
 // Set cosmetics for the q-distribution.
 
 Int_t firstNonEmptyBin = hist->FindFirstBinAbove(0);
 Double_t lowRange = hist->GetBinLowEdge(firstNonEmptyBin);
 Int_t lastNonEmptyBin = hist->FindLastBinAbove(0);
 Double_t upperRange = hist->GetBinLowEdge(lastNonEmptyBin+10);
 hist->GetXaxis()->SetRangeUser(lowRange,upperRange);  
 hist->SetFillColor(16);  
 hist->SetTitle("Fitted q-distribution");
 
} // end of void Cosmetics(TH1D *hist)

// ===========================================================================================

TFile* AccessOutputFile(TString outputFileName)
{
 // Access the common output file.
    
 TFile *outputFile = NULL;
 if(!(gSystem->AccessPathName(Form("%s%s%s",gSystem->pwd(),"/",outputFileName.Data()),kFileExists)))
 {
  outputFile = TFile::Open(outputFileName.Data(),"UPDATE");
 } else
   {
    cout<<endl;
    cout<<"WARNING: Couldn't find the file "<<outputFileName.Data()<<" in "<<endl;
    cout<<"         directory "<<gSystem->pwd()<<" !!!!"<<endl;
    cout<<endl;
    exit(0);
   }
  
 return outputFile;
 
} // end of TFile* AccessOutputFile(TString outputFileName)
  
// ===========================================================================================
 
void GetListWithHistograms(TFile *outputFile, TString analysisType, TString methodName)
{
 // Access from common output file the TDirectoryFile for FQD method and from it
 // the output list holding all objects:

 TString fileName = ""; 
 TString listName = ""; 
 // Form a file name: 
 fileName+="output";
 fileName+=methodName.Data();
 fileName+="analysis";
 fileName+=analysisType.Data();
 // Access this file:
 dirFileFQD = (TDirectoryFile*)outputFile->FindObjectAny(fileName.Data());
 // Form a list name for each method:
 if(dirFileFQD)
 {
  TList* listTemp = dirFileFQD->GetListOfKeys();
  if(listTemp && listTemp->GetEntries() == 1)
  {
   listName = listTemp->At(0)->GetName(); // to be improved - implemented better (use dynamic_cast)
   dirFileFQD->GetObject(listName.Data(),outputListFQD);
  } else
    {
     cout<<" WARNING: Accessing TList from TDirectoryFile failed miserably for FQD method !!!!"<<endl;
     cout<<"          Did you actually use FQD in the analysis?"<<endl;
     cout<<endl;
     exit(0);
    }
 } else 
   {
    cout<<" WARNING: Couldn't find a TDirectoryFile "<<fileName.Data()<<" in file "<<outputFileName.Data()<<" !!!!"<<endl;
    cout<<"          Perhaps wrong analysis type? Can be \"ESD\",\"AOD\",\"\", etc."<<endl;
    cout<<"          Or you didn't use FQD in the analysis which have produced "<<outputFileName.Data()<<"?"<<endl;
   } 
          
} // end of void GetListWithHistograms(TFile *outputFile, TString analysisType, TString methodName)
 
//===========================================================================================

void CrossCheckSettings() 
{
 // Check in this method if the settings make sense.
 
 if(sigma2Min<0.5)
 {
  cout<<endl; 
  cout<<"WARNING: According to theorists sigma2Min >= 0.5 !!!!"<<endl;
  cout<<endl; 
  exit(0);
 }
 
} // end of void CrossCheckSettings()

//===========================================================================================

void LoadLibrariesFQD(const libModes mode) {
  
  //--------------------------------------
  // Load the needed libraries most of them already loaded by aliroot
  //--------------------------------------
  gSystem->Load("libTree");
  gSystem->Load("libGeom");
  gSystem->Load("libVMC");
  gSystem->Load("libXMLIO");
  gSystem->Load("libPhysics");
  
  //----------------------------------------------------------
  // >>>>>>>>>>> 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");

  // for AliRoot
  gSystem->AddIncludePath("-I$ALICE_ROOT/include");
  gSystem->Load("libANALYSIS");
  gSystem->Load("libPWGflowBase");
  cerr<<"libPWGflowBase loaded ..."<<endl;
  
  }
  
  else if (mode==mLocalSource) {
 
    // In root inline compile
  
    // Constants  
    gROOT->LoadMacro("BaseAliFlowCommonConstants.cxx+");
    gROOT->LoadMacro("BaseAliFlowLYZConstants.cxx+");
    
    // Flow event
    gROOT->LoadMacro("BaseAliFlowVector.cxx+"); 
    gROOT->LoadMacro("BaseAliFlowTrackSimple.cxx+");    
    gROOT->LoadMacro("BaseAliFlowTrackSimpleCuts.cxx+");    
    gROOT->LoadMacro("BaseAliFlowEventSimple.cxx+");
    
    // Output histosgrams
    gROOT->LoadMacro("BaseAliFlowCommonHist.cxx+");
    gROOT->LoadMacro("BaseAliFlowCommonHistResults.cxx+");
    gROOT->LoadMacro("BaseAliFlowLYZHist1.cxx+");
    gROOT->LoadMacro("BaseAliFlowLYZHist2.cxx+");
       
    cout << "finished loading macros!" << endl;  
    
  } // end of else if (mode==mLocalSource) 
  
} // end of void LoadLibrariesFQD(const libModes mode)
Commit	Line	Data
f0e0fbcd	1	//==========================================================================//
	2	// Macro to refit the q-distribution in FQD method. By using this macro the //
	3	// fitting parameters for q-distribution can be tuned and correspondingly //
	4	// the estimate for reference flow harmonic obtained from FQD method can be //
	5	// significantly improved. //
	6	//==========================================================================//
	7
	8	// Fitting parameters for q-distribution:
	9	Double_t vStart = 0.05; // starting value for v fit
	10	Double_t vMin = 0.0; // lower bound for v fit
	11	Double_t vMax = 0.25; // upper bound for v fit
	12	Double_t sigma2Start = 0.75; // starting value for sigma^2 fit
	13	Double_t sigma2Min = 0.5; // lower bound for sigma^2 fit (according to theorists must be >= 0.5)
	14	Double_t sigma2Max = 2.5; // upper bound for sigma^2 fit
	15	Double_t treshold = 5.; // first and last bin taken for fitting are determined as the first and last bin with more than 'treshold' number of entries
	16	Bool_t finalResultIsFromSigma2Fitted = kTRUE; // final saved result is obtained with sigma^2 fitted (kTRUE) or sigma^2 fixed (kFALSE)
	17	Bool_t printOnTheScreen = kTRUE; // print or not the final results on the screen
	18	Bool_t plotResults = kTRUE; // plot on the screen q-distribution and resulting fitting functions (for non-fixed sigma^2 and fixed sigma^2)
	19
5c09ff70	20	// Name of the common output file:
	21	TString outputFileName = "AnalysisResults.root";
	22	//TString outputFileName = "outputCentrality0.root";
	23
f0e0fbcd	24	enum libModes {mLocal,mLocalSource};
	25
	26	TFile *commonOutputFile = NULL; // common output file "AnalysisResults.root"
	27	TDirectoryFile *dirFileFQD = NULL; // FQD's TDirectoryFile in "AnalysisResults.root"
	28	TList *outputListFQD = NULL; // output list holding all FQD objects
	29
5c09ff70	30	void fqd(TString analysisType="ESD", Int_t analysisMode=mLocal)
f0e0fbcd	31	{
5c09ff70	32	// 1. analysisType: type of analysis can be ESD, AOD, MC, MK, ....
5c09ff70	33	// (if type="" output files are from simulation 'on the fly')
f0e0fbcd	34	// 2. analysisMode: if analysisMode = mLocal -> analyze data on your computer using aliroot
	35	// if analysisMode = mLocalSource -> analyze data on your computer using root + source files
	36
	37	// Cross-check if the user's settings make sense:
	38	CrossCheckSettings();
	39	// Load needed libraries:
	40	LoadLibrariesFQD(analysisMode);
	41	// Get output list which holds all FQD objects:
	42	GetOutputList(analysisType);
	43	// Redo fit of q-distribution:
	44	RedoFit();
	45	// Plot q-distribution and resulting fitting functions:
	46	if(plotResults) Plot();
	47	// Save the new results in the common output file:
	48	Save();
	49
	50	} // end of void fqd(TString type="ESD", Int_t mode=mLocal)
	51
	52	// ===========================================================================================
	53
	54	void RedoFit()
	55	{
	56	// Redo the fit of q-distribution.
	57
	58	AliFlowAnalysisWithFittingQDistribution* fqd = new AliFlowAnalysisWithFittingQDistribution();
	59	fqd->GetOutputHistograms(outputListFQD);
	60	// Set new fitting parameters:
	61	fqd->SetvStart(vStart);
	62	fqd->SetvMin(vMin);
	63	fqd->SetvMax(vMax);
	64	fqd->SetSigma2Start(sigma2Start);
	65	fqd->SetSigma2Min(sigma2Min);
	66	fqd->SetSigma2Max(sigma2Max);
	67	fqd->SetTreshold(treshold);
	68	fqd->SetFinalResultIsFromSigma2Fitted(finalResultIsFromSigma2Fitted);
	69	fqd->SetPrintOnTheScreen(printOnTheScreen);
	70	// Save new fitting parameters:
	71	fqd->StoreFittingParameters();
	72	// Redo fit:
	73	fqd->Finish(kTRUE);
	74
	75	} // end of void RedoFit(TList *outputList)
	76
	77	// ===========================================================================================
	78
	79	void GetOutputList(TString analysisType)
	80	{
	81	// Get output list which holds all FQD objects.
	82
	83	// Access common output file:
f0e0fbcd	84	commonOutputFile = AccessOutputFile(outputFileName);
	85
	86	// Access from common output file the TDirectoryFile for FQD method
	87	// and from it the output list holding all objects:
	88	GetListWithHistograms(commonOutputFile,analysisType,"FQD");
	89
	90	} // end of TList* GetOutputList(TString analysisType)
	91
	92	// ===========================================================================================
	93
	94	void Plot()
	95	{
	96	// Plot q-distribution and resulting fitting functions.
	97
	98	gROOT->SetStyle("Plain"); // default color is white instead of gray
	99	gStyle->SetOptStat(0); // remove stat. box from all histos
5c09ff70	100	legend = new TLegend(0.6,0.55,0.85,0.7);
5c09ff70	101	legend->SetFillStyle(0);
f0e0fbcd	102	// q-distribution:
	103	TH1D qDistribution = dynamic_cast<TH1D>(outputListFQD->FindObject("fqDistribution"));
	104	Cosmetics(qDistribution);
5c09ff70	105	legend->AddEntry(qDistribution,"q-distribution","f");
f0e0fbcd	106	qDistribution->Draw();
	107	// resulting fitting functions:
	108	TF1 *fittingFun[2] = {NULL};
	109	TString sigmaFlag[2] = {"#sigma^{2} not fitted","#sigma^{2} fitted"};
	110	Double_t lineColors[2] = {kBlue,kRed};
	111	for(Int_t s2F=0;s2F<2;s2F++)
	112	{
	113	fittingFun[s2F] = dynamic_cast<TF1*>(outputListFQD->FindObject(Form("fFittingFunction, %s",sigmaFlag[s2F].Data())));
	114	fittingFun[s2F]->SetLineColor(lineColors[s2F]);
	115	fittingFun[s2F]->Draw("SAME");
5c09ff70	116	legend->AddEntry(fittingFun[s2F]->GetName(),sigmaFlag[s2F].Data(),"l");
f0e0fbcd	117	}
5c09ff70	118	legend->Draw("SAME");
f0e0fbcd	119
	120	} // end of Plot()
	121
	122	// ===========================================================================================
	123
	124	void Save()
	125	{
	126	// Save the new results of fitting for FQD method in the common output file.
	127
	128	dirFileFQD->Add(outputListFQD,kTRUE);
	129	dirFileFQD->Write(dirFileFQD->GetName(),TObject::kSingleKey+TObject::kOverwrite);
	130	//delete commonOutputFile;
	131
	132	} // end of void Save()
	133
	134	// ===========================================================================================
	135
	136	void Cosmetics(TH1D *hist)
	137	{
	138	// Set cosmetics for the q-distribution.
	139
	140	Int_t firstNonEmptyBin = hist->FindFirstBinAbove(0);
	141	Double_t lowRange = hist->GetBinLowEdge(firstNonEmptyBin);
	142	Int_t lastNonEmptyBin = hist->FindLastBinAbove(0);
	143	Double_t upperRange = hist->GetBinLowEdge(lastNonEmptyBin+10);
	144	hist->GetXaxis()->SetRangeUser(lowRange,upperRange);
	145	hist->SetFillColor(16);
	146	hist->SetTitle("Fitted q-distribution");
	147
	148	} // end of void Cosmetics(TH1D *hist)
	149
	150	// ===========================================================================================
	151
	152	TFile* AccessOutputFile(TString outputFileName)
	153	{
	154	// Access the common output file.
	155
	156	TFile *outputFile = NULL;
	157	if(!(gSystem->AccessPathName(Form("%s%s%s",gSystem->pwd(),"/",outputFileName.Data()),kFileExists)))
	158	{
	159	outputFile = TFile::Open(outputFileName.Data(),"UPDATE");
	160	} else
	161	{
	162	cout<<endl;
	163	cout<<"WARNING: Couldn't find the file "<<outputFileName.Data()<<" in "<<endl;
	164	cout<<" directory "<<gSystem->pwd()<<" !!!!"<<endl;
	165	cout<<endl;
	166	exit(0);
	167	}
	168
	169	return outputFile;
	170
	171	} // end of TFile* AccessOutputFile(TString outputFileName)
	172
	173	// ===========================================================================================
	174
	175	void GetListWithHistograms(TFile *outputFile, TString analysisType, TString methodName)
	176	{
	177	// Access from common output file the TDirectoryFile for FQD method and from it
	178	// the output list holding all objects:
	179
	180	TString fileName = "";
	181	TString listName = "";
	182	// Form a file name:
183	fileName+="output";
184	fileName+=methodName.Data();
185	fileName+="analysis";
186	fileName+=analysisType.Data();
187	// Access this file:
188	dirFileFQD = (TDirectoryFile*)outputFile->FindObjectAny(fileName.Data());
5c09ff70	189	// Form a list name for each method:
f0e0fbcd	190	if(dirFileFQD)
f0e0fbcd	191	{
5c09ff70	192	TList* listTemp = dirFileFQD->GetListOfKeys();
5c09ff70	193	if(listTemp && listTemp->GetEntries() == 1)
f0e0fbcd	194	{
5c09ff70	195	listName = listTemp->At(0)->GetName(); // to be improved - implemented better (use dynamic_cast)
f0e0fbcd	196	dirFileFQD->GetObject(listName.Data(),outputListFQD);
	197	} else
	198	{
5c09ff70	199	cout<<" WARNING: Accessing TList from TDirectoryFile failed miserably for FQD method !!!!"<<endl;
	200	cout<<" Did you actually use FQD in the analysis?"<<endl;
	201	cout<<endl;
	202	exit(0);
	203	}
f0e0fbcd	204	} else
f0e0fbcd	205	{
5c09ff70	206	cout<<" WARNING: Couldn't find a TDirectoryFile "<<fileName.Data()<<" in file "<<outputFileName.Data()<<" !!!!"<<endl;
	207	cout<<" Perhaps wrong analysis type? Can be \"ESD\",\"AOD\",\"\", etc."<<endl;
	208	cout<<" Or you didn't use FQD in the analysis which have produced "<<outputFileName.Data()<<"?"<<endl;
	209	}
	210
f0e0fbcd	211	} // end of void GetListWithHistograms(TFile *outputFile, TString analysisType, TString methodName)
	212
	213	//===========================================================================================
	214
	215	void CrossCheckSettings()
	216	{
	217	// Check in this method if the settings make sense.
	218
	219	if(sigma2Min<0.5)
	220	{
	221	cout<<endl;
	222	cout<<"WARNING: According to theorists sigma2Min >= 0.5 !!!!"<<endl;
	223	cout<<endl;
	224	exit(0);
	225	}
	226
	227	} // end of void CrossCheckSettings()
	228
	229	//===========================================================================================
	230
	231	void LoadLibrariesFQD(const libModes mode) {
	232
	233	//--------------------------------------
	234	// Load the needed libraries most of them already loaded by aliroot
	235	//--------------------------------------
4070f709	236	gSystem->Load("libTree");
	237	gSystem->Load("libGeom");
	238	gSystem->Load("libVMC");
	239	gSystem->Load("libXMLIO");
	240	gSystem->Load("libPhysics");
f0e0fbcd	241
	242	//----------------------------------------------------------
	243	// >>>>>>>>>>> Local mode <<<<<<<<<<<<<<
	244	//----------------------------------------------------------
	245	if (mode==mLocal) {
	246	//--------------------------------------------------------
	247	// If you want to use already compiled libraries
	248	// in the aliroot distribution
	249	//--------------------------------------------------------
	250
	251	//==================================================================================
	252	//load needed libraries:
	253	gSystem->AddIncludePath("-I$ROOTSYS/include");
4070f709	254	gSystem->Load("libTree");
f0e0fbcd	255
	256	// for AliRoot
	257	gSystem->AddIncludePath("-I$ALICE_ROOT/include");
4070f709	258	gSystem->Load("libANALYSIS");
4070f709	259	gSystem->Load("libPWGflowBase");
4bc96904	260	cerr<<"libPWGflowBase loaded ..."<<endl;
f0e0fbcd	261
	262	}
	263
	264	else if (mode==mLocalSource) {
	265
	266	// In root inline compile
	267
	268	// Constants
2e311896	269	gROOT->LoadMacro("BaseAliFlowCommonConstants.cxx+");
2e311896	270	gROOT->LoadMacro("BaseAliFlowLYZConstants.cxx+");
f0e0fbcd	271
f0e0fbcd	272	// Flow event
2e311896	273	gROOT->LoadMacro("BaseAliFlowVector.cxx+");
	274	gROOT->LoadMacro("BaseAliFlowTrackSimple.cxx+");
	275	gROOT->LoadMacro("BaseAliFlowTrackSimpleCuts.cxx+");
	276	gROOT->LoadMacro("BaseAliFlowEventSimple.cxx+");
f0e0fbcd	277
f0e0fbcd	278	// Output histosgrams
2e311896	279	gROOT->LoadMacro("BaseAliFlowCommonHist.cxx+");
	280	gROOT->LoadMacro("BaseAliFlowCommonHistResults.cxx+");
	281	gROOT->LoadMacro("BaseAliFlowLYZHist1.cxx+");
	282	gROOT->LoadMacro("BaseAliFlowLYZHist2.cxx+");
f0e0fbcd	283
	284	cout << "finished loading macros!" << endl;
	285
	286	} // end of else if (mode==mLocalSource)
	287
	288	} // end of void LoadLibrariesFQD(const libModes mode)