//Formats a given input string such that it follows the iEta, iPhi, iSlice naming convention
//Will take text written in quotes without requiring a variable
string AnalyzeResponseUniformity::getNameByIndex(int iEta, int iPhi, int iSlice, const std::string & strInputPrefix, const std::string & strInputName){
//Variable Declaration
string ret_Name = "";
if (iSlice > -1) {
ret_Name = strInputPrefix + "_iEta" + getString(iEta) + "iPhi" + getString(iPhi) + "Slice" + getString(iSlice) + "_" + strInputName;
}
else if (iPhi > -1){ //Case: Specific (iEta,iPhi) sector
ret_Name = strInputPrefix + "_iEta" + getString(iEta) + "iPhi" + getString(iPhi) + "_" + strInputName;
} //End Case: Specific (iEta,iPhi) sector
else{ //Case: iEta Sector, sum over sector's iPhi
ret_Name = strInputPrefix + "_iEta" + getString(iEta) + "_" + strInputName;
} //End Case: iEta Sector, sum over sector's iPhi
return ret_Name;
} //End AnalyzeResponseUniformity::getNameByIndex()
void ParameterLoaderAnaysis::loadAnalysisParametersHistograms(ifstream &inputFileStream, HistoSetup &hSetup){
//Variable Declaration
bool bExitSuccess;
pair<string,string> pair_strParam;
string strLine;
vector<string> vec_strList;
//Loop through the section
while ( getlineNoSpaces(inputFileStream, strLine) ) {
bExitSuccess = false;
//Debugging
//cout<<"loadAnalysisParametersHistograms (Base Level); strLine = " << strLine << endl;
//Does the user want to comment out this line?
if ( 0 == strLine.compare(0,1,"#") ) continue;
//Has this section ended?
if ( 0 == strLine.compare(strSecEnd_Uniformity_Hiso) ) break;
//Get Parameter Pairing
pair_strParam = Timing::getParsedLine(strLine, bExitSuccess);
//transform field name to upper case
transform(pair_strParam.first.begin(),pair_strParam.first.end(),pair_strParam.first.begin(),toupper);
//Parse pair
if (bExitSuccess) { //Case: Parameter Fetched Successfully
if( 0 == pair_strParam.first.compare("HISTO_BINRANGE") ){
//Get comma separated list
vec_strList = Timing::getCharSeparatedList(pair_strParam.second,',');
//Debugging
//for(int i=0; i<vec_strList.size(); ++i){
// cout<<"vec_strList["<<i<<"] = " << vec_strList[i] << endl;
//}
if (vec_strList.size() >= 2) { //Case: at least 2 numbers
//Fetch
hSetup.fHisto_xLower = Timing::stofSafe(pair_strParam.first, vec_strList[0]);
hSetup.fHisto_xUpper = Timing::stofSafe(pair_strParam.first, vec_strList[1]);
//Reset to ensure they are both correctly lower & upper values
hSetup.fHisto_xLower = std::min(hSetup.fHisto_xLower, hSetup.fHisto_xUpper);
hSetup.fHisto_xUpper = std::max(hSetup.fHisto_xLower, hSetup.fHisto_xUpper);
//Tell the user they entered more than what was expected
if (vec_strList.size() > 2) { //Case: 3 or more numbers
printClassMethodMsg("ParameterLoaderAnaysis","loadAnalysisParametersHistograms", ( "Sorry you entered 3 or more numbers for " + pair_strParam.first + "\n" ).c_str() );
printClassMethodMsg("ParameterLoaderAnaysis","loadAnalysisParametersHistograms", "\tI have only used the first two and ignored the rest:\n" );
printClassMethodMsg("ParameterLoaderAnaysis","loadAnalysisParametersHistograms", ("\t" + getString( hSetup.fHisto_xLower ) ).c_str() );
printClassMethodMsg("ParameterLoaderAnaysis","loadAnalysisParametersHistograms", ("\t" + getString( hSetup.fHisto_xUpper ) ).c_str() );
} //End Case: 3 or more numbers
} //End Case: at least 2 numbers
else{ //Case: Not enough numbers
if (vec_strList.size() == 1) { //Case: only 1 number entered
hSetup.fHisto_xUpper = Timing::stofSafe(pair_strParam.first, vec_strList[0]);
} //End Case: only 1 number entered
//Reset to ensure they are both correctly lower & upper values
hSetup.fHisto_xLower = std::min(hSetup.fHisto_xLower, hSetup.fHisto_xUpper);
hSetup.fHisto_xUpper = std::max(hSetup.fHisto_xLower, hSetup.fHisto_xUpper);
//Output to User
printClassMethodMsg("ParameterLoaderAnaysis","loadAnalysisParametersHistograms", ( "Sorry I was expecting a comma separated list of 2 numbers for: " + pair_strParam.first + "\n" ).c_str() );
printClassMethodMsg("ParameterLoaderAnaysis","loadAnalysisParametersHistograms", "\tRight now I have set:\n" );
printClassMethodMsg("ParameterLoaderAnaysis","loadAnalysisParametersHistograms", ("\t" + getString( hSetup.fHisto_xLower ) ).c_str() );
printClassMethodMsg("ParameterLoaderAnaysis","loadAnalysisParametersHistograms", ("\t" + getString( hSetup.fHisto_xUpper ) ).c_str() );
} //End Case: Not enough numbers
} //End Case: Assign Histo Bin Range
/*else if( 0 == pair_strParam.first.compare("HISTO_NAME") ) {
hSetup.strHisto_Name = pair_strParam.second;
}*/
else if( 0 == pair_strParam.first.compare("HISTO_NUMBINS") ){
hSetup.iHisto_nBins = Timing::stoiSafe(pair_strParam.first, pair_strParam.second);
}
else if( 0 == pair_strParam.first.compare("HISTO_XTITLE") ){
hSetup.strHisto_Title_X = pair_strParam.second;
}
else if( 0 == pair_strParam.first.compare("HISTO_YTITLE") ){
hSetup.strHisto_Title_Y = pair_strParam.second;
}
else if( 0 == pair_strParam.first.compare("PERFORM_FIT") ){
hSetup.bFit = convert2bool(pair_strParam.second, bExitSuccess);
}
else{ //Case: Parameter not recognized
printClassMethodMsg("ParameterLoaderAnaysis","loadAnalysisParametersHistograms","Error!!! Parameter Not Recognized:\n");
printClassMethodMsg("ParameterLoaderAnaysis","loadAnalysisParametersHistograms",( "\tField = " + pair_strParam.first + "\n" ).c_str() );
printClassMethodMsg("ParameterLoaderAnaysis","loadAnalysisParametersHistograms",( "\tValue = " + pair_strParam.second + "\n" ).c_str() );
} //End Case: Parameter not recognized
} //End Case: Parameter Fetched Successfully
else{ //Case: Parameter was NOT fetched Successfully
printClassMethodMsg("ParameterLoaderAnaysis","loadAnalysisParametersHistograms","Error!!! I didn't parse parameter correctly\n");
printClassMethodMsg("ParameterLoaderAnaysis","loadAnalysisParametersHistograms",("\tCurrent line: " + strLine).c_str() );
} //End Case: Parameter was NOT fetched Successfully
} //End Loop through Section
return;
} //End ParameterLoaderAnalysis::loadAnalysisParametersHistograms() - Histogram specific
//Loops over all stored clusters in detMPGD and Book histograms for the full detector
void AnalyzeResponseUniformity::fillHistos(){
//Variable Declaration
//vector<Cluster> vec_clust;
//Initialize Summary Histograms
//hEta_v_SliceNum_Occupancy = std::make_shared<TH2F>( TH2F("hEta_v_SliceNum_Occupancy","",3. * aSetup.iUniformityGranularity, 1, 3. * aSetup.iUniformityGranularity + 1, detMPGD.map_sectorsEta.size(), 1, detMPGD.map_sectorsEta.size() + 1 ) );
//Loop Over Stored iEta Sectors
for (auto iterEta = detMPGD.map_sectorsEta.begin(); iterEta != detMPGD.map_sectorsEta.end(); ++iterEta) { //Loop Over iEta Sectors
//Grab Eta Sector width (for ClustPos Histo)
aSetup.histoSetup_clustPos.fHisto_xLower = -0.5*(*iterEta).second.fWidth;
aSetup.histoSetup_clustPos.fHisto_xUpper = 0.5*(*iterEta).second.fWidth;
//Initialize iEta Histograms - 1D
(*iterEta).second.hEta_ClustADC = make_shared<TH1F>(getHistogram((*iterEta).first, -1, aSetup.histoSetup_clustADC ) );
(*iterEta).second.hEta_ClustMulti = make_shared<TH1F>(getHistogram((*iterEta).first, -1, aSetup.histoSetup_clustMulti ) );
(*iterEta).second.hEta_ClustPos = make_shared<TH1F>(getHistogram((*iterEta).first, -1, aSetup.histoSetup_clustPos ) );
(*iterEta).second.hEta_ClustSize = make_shared<TH1F>(getHistogram((*iterEta).first, -1, aSetup.histoSetup_clustSize ) );
(*iterEta).second.hEta_ClustTime = make_shared<TH1F>(getHistogram((*iterEta).first, -1, aSetup.histoSetup_clustTime ) );
//Initialize iEta Histograms - 2D
(*iterEta).second.hEta_ClustADC_v_ClustPos = make_shared<TH2F>( TH2F( ("hiEta" + getString( (*iterEta).first ) + "_ClustADC_v_ClustPos").c_str(),"Response Uniformity", 3. * aSetup.iUniformityGranularity,-0.5*(*iterEta).second.fWidth,0.5*(*iterEta).second.fWidth,300,0,15000) );
(*iterEta).second.hEta_ClustADC_v_ClustPos->Sumw2();
//Debugging
//cout<<"(*iterEta).second.hEta_ClustADC->GetName() = " << (*iterEta).second.hEta_ClustADC->GetName() << endl;
//Loop Over Stored iPhi Sectors
for (auto iterPhi = (*iterEta).second.map_sectorsPhi.begin(); iterPhi != (*iterEta).second.map_sectorsPhi.end(); ++iterPhi) { //Loop Over iPhi Sectors
//Initialize iPhi Histograms - 1D
(*iterPhi).second.hPhi_ClustADC = make_shared<TH1F>(getHistogram( (*iterEta).first, (*iterPhi).first, aSetup.histoSetup_clustADC ) );
(*iterPhi).second.hPhi_ClustMulti = make_shared<TH1F>(getHistogram( (*iterEta).first, (*iterPhi).first, aSetup.histoSetup_clustMulti ) );
(*iterPhi).second.hPhi_ClustSize = make_shared<TH1F>(getHistogram( (*iterEta).first, (*iterPhi).first, aSetup.histoSetup_clustSize ) );
(*iterPhi).second.hPhi_ClustTime = make_shared<TH1F>(getHistogram( (*iterEta).first, (*iterPhi).first, aSetup.histoSetup_clustTime ) );
//Initialize iPhi Histograms - 2D
(*iterPhi).second.hPhi_ClustADC_v_ClustPos = make_shared<TH2F>( TH2F( ("hiEta" + getString( (*iterEta).first ) + "iPhi" + getString( (*iterPhi).first ) + "_ClustADC_v_ClustPos").c_str(),"Response Uniformity", aSetup.iUniformityGranularity, (*iterPhi).second.fPos_Xlow, (*iterPhi).second.fPos_Xhigh,300,0,15000) );
(*iterPhi).second.hPhi_ClustADC_v_ClustPos->Sumw2();
//Loop Over Stored Clusters
for (auto iterClust = (*iterPhi).second.vec_clusters.begin(); iterClust != (*iterPhi).second.vec_clusters.end(); ++iterClust) { //Loop Over Stored Clusters
//Fill iEta Histograms
(*iterEta).second.hEta_ClustADC->Fill( (*iterClust).fADC );
//(*iterEta).second.hEta_ClustMulti->Fill( (*iterClust). );
(*iterEta).second.hEta_ClustPos->Fill( (*iterClust).fPos_X );
(*iterEta).second.hEta_ClustSize->Fill( (*iterClust).iSize );
(*iterEta).second.hEta_ClustTime->Fill( (*iterClust).iTimeBin );
(*iterEta).second.hEta_ClustADC_v_ClustPos->Fill( (*iterClust).fPos_X, (*iterClust).fADC );
//Fill iPhi Histograms
(*iterPhi).second.hPhi_ClustADC->Fill( (*iterClust).fADC );
//(*iterPhi).second.hPhi_ClustMulti
(*iterPhi).second.hPhi_ClustSize->Fill( (*iterClust).iSize);
(*iterPhi).second.hPhi_ClustTime->Fill( (*iterClust).iTimeBin);
(*iterPhi).second.hPhi_ClustADC_v_ClustPos->Fill( (*iterClust).fPos_X, (*iterClust).fADC );
} //End Loop Over Stored Clusters
//Slices
//Now that all clusters have been analyzed we extract the slices
for (int i=1; i<= aSetup.iUniformityGranularity; ++i) { //Loop Over Slices
//Create the slice
SectorSlice slice;
//string strSliceName = "hiEta" + getString( (*iterEta).first ) + "iPhi" + getString( (*iterPhi).first ) + "Slice" + getString(i) + "_ClustADC";
//Grab ADC spectrum for this slice
slice.hSlice_ClustADC = make_shared<TH1F>( *( (TH1F*) (*iterPhi).second.hPhi_ClustADC_v_ClustPos->ProjectionY( ("hiEta" + getString( (*iterEta).first ) + "iPhi" + getString( (*iterPhi).first ) + "Slice" + getString(i) + "_ClustADC").c_str(),i,i,"") ) );
//Store position information for this slice
slice.fPos_Center = (*iterPhi).second.hPhi_ClustADC_v_ClustPos->GetXaxis()->GetBinCenter(i);
slice.fWidth = (*iterPhi).second.hPhi_ClustADC_v_ClustPos->GetXaxis()->GetBinWidth(i);
//Store the slice
(*iterPhi).second.map_slices[i] = slice;
} //End Loop Over Slices
} //End Loop Over iPhi Sectors
std::cout<<"(*iterEta).second.hEta_ClustADC->Integral() = " << (*iterEta).second.hEta_ClustADC->Integral() << std::endl;
} //End Loop Over iEta Sectors
} //End AnalyzeResponseUniformity::fillHistos() - Full Detector
TF1 AnalyzeResponseUniformity::getFit(int iEta, int iPhi, int iSlice, HistoSetup & setupHisto, shared_ptr<TH1F> hInput, TSpectrum &specInput ){
//Variable Declaration
float fLimit_Max = setupHisto.fHisto_xUpper, fLimit_Min = setupHisto.fHisto_xLower;
vector<string>::const_iterator iterVec_IGuess; //Iterator to use for setting initial guess of fit
vector<float> vec_fFitRange;
for (auto iterRange = aSetup.histoSetup_clustADC.vec_strFit_Range.begin(); iterRange != aSetup.histoSetup_clustADC.vec_strFit_Range.end(); ++iterRange) { //Loop Over Fit Range
vec_fFitRange.push_back( getFitBoundary( (*iterRange), hInput, specInput ) );
} //End Loop Over Fit Range
if (vec_fFitRange.size() > 1) {
fLimit_Min = (*std::min_element(vec_fFitRange.begin(), vec_fFitRange.end() ) );
fLimit_Max = (*std::max_element(vec_fFitRange.begin(), vec_fFitRange.end() ) );
}
TF1 ret_Func( getNameByIndex(iEta, iPhi, iSlice, "fit", setupHisto.strHisto_Name).c_str(), setupHisto.strFit_Formula.c_str(), fLimit_Min, fLimit_Max);
//Check to see if the number of parameters in the TF1 meets the expectation
if ( ret_Func.GetNpar() < setupHisto.vec_strFit_ParamIGuess.size() || ret_Func.GetNpar() < setupHisto.vec_strFit_ParamLimit_Min.size() || ret_Func.GetNpar() < setupHisto.vec_strFit_ParamLimit_Max.size() ) { //Case: Set points for initial parameters do not meet expectations
printClassMethodMsg("AnalyzeResponseUniformity","getFit","Error! Number of Parameters in Function Less Than Requested Initial Guess Parameters!");
printClassMethodMsg("AnalyzeResponseUniformity","getFit", ("\tNum Parameter: " + getString( ret_Func.GetNpar() ) ).c_str() );
printClassMethodMsg("AnalyzeResponseUniformity","getFit", ("\tNum Initial Guesses: " + getString( setupHisto.vec_strFit_ParamIGuess.size() ) ).c_str() );
printClassMethodMsg("AnalyzeResponseUniformity","getFit", ("\tNum Initial Guess Limits (Min): " + getString( setupHisto.vec_strFit_ParamLimit_Min.size() ) ).c_str() );
printClassMethodMsg("AnalyzeResponseUniformity","getFit", ("\tNum Initial Guess Limits (Max): " + getString( setupHisto.vec_strFit_ParamLimit_Max.size() ) ).c_str() );
printClassMethodMsg("AnalyzeResponseUniformity","getFit", "No Initial Parameters Have Been Set! Please Cross-Check Input Analysis Config File" );
return ret_Func;
} //End Case: Set points for initial parameters do not meet expectations
//Set Fit Parameters - Initial Value
//------------------------------------------------------
//Keywords are AMPLITUDE, MEAN, PEAK, SIGMA
for (int i=0; i<setupHisto.vec_strFit_ParamIGuess.size(); ++i) { //Loop over parameters - Initial Guess
iterVec_IGuess = std::find(vec_strSupportedKeywords.begin(), vec_strSupportedKeywords.end(), setupHisto.vec_strFit_ParamIGuess[i]);
if ( iterVec_IGuess == vec_strSupportedKeywords.end() ) { //Case: No Keyword Found; Try to set a Numeric Value
ret_Func.SetParameter(i, stofSafe( setupHisto.vec_strFit_ParamIGuess[i] ) );
} //End Case: No Keyword Found; Try to set a Numeric Value
else{ //Case: Keyword Found; Set Value based on Keyword
ret_Func.SetParameter(i, getValByKeyword( (*iterVec_IGuess), hInput, specInput ) );
} //End Case: Keyword Found; Set Value based on Keyword
} //End Loop over parameters - Initial Guess
//Set Fit Parameters - Boundaries
//------------------------------------------------------
if (setupHisto.vec_strFit_ParamLimit_Min.size() == setupHisto.vec_strFit_ParamLimit_Max.size() ) { //Check: Stored Parameter Limits Match
//Here we use vec_strFit_ParamLimit_Min but we know it has the same number of parameters as vec_strFit_ParamLimit_Max
//For each fit parameter, set the boundary
for (int i=0; i<setupHisto.vec_strFit_ParamLimit_Min.size(); ++i) { //Loop over boundary parameters
fLimit_Min = getFitBoundary(setupHisto.vec_strFit_ParamLimit_Min[i], hInput, specInput);
fLimit_Max = getFitBoundary(setupHisto.vec_strFit_ParamLimit_Max[i], hInput, specInput);
//cout<<"(fLimit_Min, fLimit_Max) = (" << fLimit_Min << ", " << fLimit_Max << ")\n";
(fLimit_Max > fLimit_Min) ? ret_Func.SetParLimits(i, fLimit_Min, fLimit_Max ) : ret_Func.SetParLimits(i, fLimit_Max, fLimit_Min );
} //End Loop over boundary parameters
} //End Check: Stored Parameter Limits Match
//Set Fit Parameters - Fixed?
//------------------------------------------------------
//Placeholder; maybe we add functionality in the future
//Set Other Fit Data Members
//------------------------------------------------------
ret_Func.SetLineColor(kRed);
ret_Func.SetLineWidth(3);
//Delete Pointers
//delete iterVec_IGuess;
//Return fit
//------------------------------------------------------
return ret_Func;
} //End AnalyzeResponseUniformity::getFit()
void AnalyzeResponseUniformity::storeFits( string strOutputROOTFileName, std::string strOption ){
//Variable Declaration
TFile * ptr_fileOutput = new TFile(strOutputROOTFileName.c_str(), strOption.c_str(),"",1);
//Check if File Failed to Open Correctly
if ( !ptr_fileOutput->IsOpen() || ptr_fileOutput->IsZombie() ) {
printClassMethodMsg("AnalyzeResponseUniformity","storeFits","Error: File I/O");
printROOTFileStatus(ptr_fileOutput);
printClassMethodMsg("AnalyzeResponseUniformity","storeFits", "\tPlease cross check input file name, option, and the execution directory\n" );
printClassMethodMsg("AnalyzeResponseUniformity","storeFits", "\tExiting; No Fits have been stored!\n" );
return;
} //End Check if File Failed to Open Correctly
//Loop Over Stored iEta Sectors
for (auto iterEta = detMPGD.map_sectorsEta.begin(); iterEta != detMPGD.map_sectorsEta.end(); ++iterEta) { //Loop Over iEta Sectors
//Get Directory
//-------------------------------------
//Check to see if the directory exists already
TDirectory *dir_SectorEta = ptr_fileOutput->GetDirectory( ( "SectorEta" + getString( (*iterEta).first ) ).c_str(), false, "GetDirectory" );
//If the above pointer is null the directory does NOT exist, create it
if (dir_SectorEta == nullptr) { //Case: Directory did not exist in file, CREATE
dir_SectorEta = ptr_fileOutput->mkdir( ( "SectorEta" + getString( (*iterEta).first ) ).c_str() );
} //End Case: Directory did not exist in file, CREATE
//Debugging
//cout<<"dir_SectorEta->GetName() = " << dir_SectorEta->GetName()<<endl;
//Store Fits - SectorEta Level
//-------------------------------------
dir_SectorEta->cd();
(*iterEta).second.gEta_ClustADC_Fit_NormChi2->Write();
(*iterEta).second.gEta_ClustADC_Fit_PkPos->Write();
(*iterEta).second.gEta_ClustADC_Fit_Failures->Write();
(*iterEta).second.gEta_ClustADC_Spec_NumPks->Write();
(*iterEta).second.gEta_ClustADC_Spec_PkPos->Write();
//Loop Over Stored iPhi Sectors within this iEta Sector
for (auto iterPhi = (*iterEta).second.map_sectorsPhi.begin(); iterPhi != (*iterEta).second.map_sectorsPhi.end(); ++iterPhi) { //Loop Over Stored iPhi Sectors
//Get Directory
//-------------------------------------
//Check to see if the directory exists already
TDirectory *dir_SectorPhi = dir_SectorEta->GetDirectory( ( "SectorPhi" + getString( (*iterPhi).first ) ).c_str(), false, "GetDirectory" );
//If the above pointer is null the directory does NOT exist, create it
if (dir_SectorPhi == nullptr) { //Case: Directory did not exist in file, CREATE
dir_SectorPhi = dir_SectorEta->mkdir( ( "SectorPhi" + getString( (*iterPhi).first ) ).c_str() );
} //End Case: Directory did not exist in file, CREATE
//Debugging
//cout<<"dir_SectorPhi->GetName() = " << dir_SectorPhi->GetName()<<endl;
//Store Fits - SectorPhi Level
//-------------------------------------
dir_SectorPhi->cd();
//No Fits defined at this level - yet
//Slices
//Now that all clusters have been analyzed we extract the slices
for (auto iterSlice = (*iterPhi).second.map_slices.begin(); iterSlice != (*iterPhi).second.map_slices.end(); ++iterSlice ) { //Loop Over Slices
//int iSliceCount = std::distance( (*iterPhi).second.map_slices.begin(), iterSlice );
//Get Directory
//-------------------------------------
//Check to see if the directory exists already
TDirectory *dir_Slice = dir_SectorPhi->GetDirectory( ( "Slice" + getString( (*iterSlice).first ) ).c_str(), false, "GetDirectory" );
//If the above pointer is null the directory does NOT exist, create it
if (dir_Slice == nullptr) { //Case: Directory did not exist in file, CREATE
dir_Slice = dir_SectorPhi->mkdir( ( "Slice" + getString( (*iterSlice).first ) ).c_str() );
} //End Case: Directory did not exist in file, CREATE
//Store Fits - Slice Level
//-------------------------------------
dir_Slice->cd();
//(*iterSlice).second.pmrkSlice_ClustADC->Write( getNameByIndex( (*iterEta).first, (*iterPhi).first, (*iterSlice).first, "PeakMrk", "clustADC" ).c_str() );
(*iterSlice).second.fitSlice_ClustADC->Write();
} //End Loop Over Slices
} //End Loop Over Stored iPhi Sectors
} //End Loop Over Stored iEta Sectors
//Close the ROOT file
ptr_fileOutput->Close();
return;
} //End storeHistos()
//Stores booked histograms (for those histograms that are non-null)
void AnalyzeResponseUniformity::storeHistos( string strOutputROOTFileName, std::string strOption ){
//Variable Declaration
//std::shared_ptr<TFile> ptr_fileOutput;
TFile * ptr_fileOutput = new TFile(strOutputROOTFileName.c_str(), strOption.c_str(),"",1);
//Assign the TFile to the ptr_fileOutput
//ptr_fileOutput = make_shared<TFile>(TFile(strOutputROOTFileName.c_str(), strOption.c_str(),"",1) );
//Check if File Failed to Open Correctly
if ( !ptr_fileOutput->IsOpen() || ptr_fileOutput->IsZombie() ) {
printClassMethodMsg("AnalyzeResponseUniformity","storeHistos","Error: File I/O");
printROOTFileStatus(ptr_fileOutput);
printClassMethodMsg("AnalyzeResponseUniformity","storeHistos", "\tPlease cross check input file name, option, and the execution directory\n" );
printClassMethodMsg("AnalyzeResponseUniformity","storeHistos", "\tExiting; No Histograms have been stored!\n" );
return;
} //End Check if File Failed to Open Correctly
//Loop over ieta's
//Create/Load file structure
//Store ieta level histograms
//Loop over iphi's within ieta's
//Create/Load file structure
//Store iphi level histograms
//Loop over slices
//Create/Load file structure
//store slice level histograms
//Close File
//Loop Over Stored iEta Sectors
for (auto iterEta = detMPGD.map_sectorsEta.begin(); iterEta != detMPGD.map_sectorsEta.end(); ++iterEta) { //Loop Over iEta Sectors
//Get Directory
//-------------------------------------
//Check to see if the directory exists already
TDirectory *dir_SectorEta = ptr_fileOutput->GetDirectory( ( "SectorEta" + getString( (*iterEta).first ) ).c_str(), false, "GetDirectory" );
//If the above pointer is null the directory does NOT exist, create it
if (dir_SectorEta == nullptr) { //Case: Directory did not exist in file, CREATE
dir_SectorEta = ptr_fileOutput->mkdir( ( "SectorEta" + getString( (*iterEta).first ) ).c_str() );
} //End Case: Directory did not exist in file, CREATE
//Debugging
cout<<"dir_SectorEta->GetName() = " << dir_SectorEta->GetName()<<endl;
//Store Histograms - SectorEta Level
//-------------------------------------
dir_SectorEta->cd();
(*iterEta).second.hEta_ClustADC->Write();
(*iterEta).second.hEta_ClustPos->Write();
(*iterEta).second.hEta_ClustSize->Write();
(*iterEta).second.hEta_ClustTime->Write();
(*iterEta).second.hEta_ClustADC_v_ClustPos->Write();
//Loop Over Stored iPhi Sectors within this iEta Sector
for (auto iterPhi = (*iterEta).second.map_sectorsPhi.begin(); iterPhi != (*iterEta).second.map_sectorsPhi.end(); ++iterPhi) { //Loop Over Stored iPhi Sectors
//Get Directory
//-------------------------------------
//Check to see if the directory exists already
TDirectory *dir_SectorPhi = dir_SectorEta->GetDirectory( ( "SectorPhi" + getString( (*iterPhi).first ) ).c_str(), false, "GetDirectory" );
//If the above pointer is null the directory does NOT exist, create it
if (dir_SectorPhi == nullptr) { //Case: Directory did not exist in file, CREATE
dir_SectorPhi = dir_SectorEta->mkdir( ( "SectorPhi" + getString( (*iterPhi).first ) ).c_str() );
} //End Case: Directory did not exist in file, CREATE
//Debugging
cout<<"dir_SectorPhi->GetName() = " << dir_SectorPhi->GetName()<<endl;
//Store Histograms - SectorPhi Level
//-------------------------------------
dir_SectorPhi->cd();
(*iterPhi).second.hPhi_ClustADC->Write();
(*iterPhi).second.hPhi_ClustSize->Write();
(*iterPhi).second.hPhi_ClustTime->Write();
(*iterPhi).second.hPhi_ClustADC_v_ClustPos->Write();
//Loop through Slices
//To be implemented
//Slices
//Now that all clusters have been analyzed we extract the slices
for (auto iterSlice = (*iterPhi).second.map_slices.begin(); iterSlice != (*iterPhi).second.map_slices.end(); ++iterSlice ) { //Loop Over Slices
//int iSliceCount = std::distance( (*iterPhi).second.map_slices.begin(), iterSlice ) + 1;
//Get Directory
//-------------------------------------
//Check to see if the directory exists already
TDirectory *dir_Slice = dir_SectorPhi->GetDirectory( ( "Slice" + getString( (*iterSlice).first ) ).c_str(), false, "GetDirectory" );
//If the above pointer is null the directory does NOT exist, create it
if (dir_Slice == nullptr) { //Case: Directory did not exist in file, CREATE
dir_Slice = dir_SectorPhi->mkdir( ( "Slice" + getString( (*iterSlice).first ) ).c_str() );
} //End Case: Directory did not exist in file, CREATE
//Store Histograms - Slice Level
//-------------------------------------
dir_Slice->cd();
(*iterSlice).second.hSlice_ClustADC->Write();
} //End Loop Over Slices
} //End Loop Over Stored iPhi Sectors
} //End Loop Over Stored iEta Sectors
//Close the ROOT file
ptr_fileOutput->Close();
return;
} //End storeHistos()
