//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()