#include <CalibrationHistograms.h>

Inheritance diagram for CalibrationHistograms:

Public Member Functions
	CalibrationHistograms (const edm::ParameterSet &pset, DQMStore *, const sistrip::RunType &task=sistrip::CALIBRATION)

void	histoAnalysis (bool debug) override

void	printAnalyses () override

void	save (std::string &filename, uint32_t run_number=0, std::string partitionName="")

	~CalibrationHistograms () override

Public Member Functions inherited from CommissioningHistograms
	CommissioningHistograms (const edm::ParameterSet &pset, DQMStore *const, const sistrip::RunType &)

	CommissioningHistograms ()

virtual void	configure (const edm::ParameterSet &, const edm::EventSetup &)

void	createCollations (const std::vector< std::string > &)

virtual void	createSummaryHisto (const sistrip::Monitorable &, const sistrip::Presentation &, const std::string &top_level_dir, const sistrip::Granularity &)

void	extractHistograms (const std::vector< std::string > &)

virtual void	printSummary ()

void	remove (std::string pattern="")

void	save (std::string &filename, uint32_t run_number=0, std::string partitionName="")

virtual	~CommissioningHistograms ()

Private Attributes
float	targetDecayTime_

float	targetRiseTime_

bool	tuneSimultaneously_

Additional Inherited Members
Public Types inherited from CommissioningHistograms
typedef std::map< uint32_t, CommissioningAnalysis * >	Analyses

typedef Analyses::iterator	Analysis

typedef SummaryPlotFactory< CommissioningAnalysis * >	Factory

typedef std::map< uint32_t, uint32_t >	FedToFecMap

typedef std::vector< Histo * >	Histos

typedef std::map< uint32_t, Histos >	HistosMap

Static Public Member Functions inherited from CommissioningHistograms
static void	copyCustomInformation (DQMStore *const, const std::vector< std::string > &)

static uint32_t	runNumber (DQMStore *const, const std::vector< std::string > &)

static sistrip::RunType	runType (DQMStore *const, const std::vector< std::string > &)

Protected Member Functions inherited from CommissioningHistograms
DQMStore *const	bei () const

void	clearHistosMap ()

Analyses &	data (bool getMaskedData=false)

Factory *const	factory ()

TH1 *	histogram (const sistrip::Monitorable &, const sistrip::Presentation &, const sistrip::View &, const std::string &directory, const uint32_t &xbins, const float &xlow=1.sistrip::invalid_, const float &xhigh=1.sistrip::invalid_)

const HistosMap &	histos () const

const FedToFecMap &	mapping () const

void	printHistosMap ()

const edm::ParameterSet &	pset () const

const sistrip::RunType &	task () const

Protected Attributes inherited from CommissioningHistograms
std::unique_ptr< Factory >	factory_

Detailed Description

Definition at line 9 of file CalibrationHistograms.h.

Constructor & Destructor Documentation

CalibrationHistograms::CalibrationHistograms	(	const edm::ParameterSet &	pset,
		DQMStore *	bei,
		const sistrip::RunType &	task = `sistrip::CALIBRATION`
	)

Definition at line 26 of file CalibrationHistograms.cc.

References sistrip::CALIBRATION_SCAN, sistrip::CALIBRATION_SCAN_DECO, edm::ParameterSet::existsAs(), CommissioningHistograms::factory_, edm::ParameterSet::getParameter(), LogTrace, sistrip::mlDqmClient_, or, CommissioningHistograms::pset(), targetDecayTime_, targetRiseTime_, and tuneSimultaneously_.

   : CommissioningHistograms(pset.getParameter<edm::ParameterSet>("CalibrationParameters"),bei,task){  
 
   LogTrace(mlDqmClient_) 
     << "[CalibrationHistograms::" << __func__ << "]"
        << " Constructing object...";
   
   if(task == sistrip::CALIBRATION_SCAN or task == sistrip::CALIBRATION_SCAN_DECO)
     factory_ = auto_ptr<CalibrationScanSummaryFactory>( new CalibrationScanSummaryFactory ); 
   else
     factory_ = auto_ptr<CalibrationSummaryFactory>( new CalibrationSummaryFactory ); 
 
   targetRiseTime_     = this->pset().existsAs<double>("targetRiseTime") ? this->pset().getParameter<double>("targetRiseTime") : 50;
   targetDecayTime_    = this->pset().existsAs<double>("targetDecayTime") ? this->pset().getParameter<double>("targetDecayTime") : 125;
   tuneSimultaneously_ = this->pset().existsAs<bool>("tuneSimultaneously") ? this->pset().getParameter<bool>("tuneSimultaneously") : false;
 }

CalibrationHistograms::~CalibrationHistograms ( )

override

Definition at line 47 of file CalibrationHistograms.cc.

References LogTrace, and sistrip::mlDqmClient_.

                                               {
   LogTrace(mlDqmClient_) 
        << "[CalibrationHistograms::" << __func__ << "]"
        << " Deleting object...";
 }

Member Function Documentation

void CalibrationHistograms::histoAnalysis ( bool debug )

overridevirtual

Reimplemented from CommissioningHistograms.

Definition at line 55 of file CalibrationHistograms.cc.

References patPFMETCorrections_cff::algo, CommissioningAlgorithm::analysis(), begin, sistrip::CALIBRATION_DECO, sistrip::CALIBRATION_SCAN, sistrip::CALIBRATION_SCAN_DECO, CommissioningHistograms::data(), extract(), CalibrationScanAlgorithm::fillTunedObservables(), CommissioningHistograms::histos(), sistrip::mlDqmClient_, or, CommissioningHistograms::pset(), targetDecayTime_, targetRiseTime_, CommissioningHistograms::task(), CalibrationScanAlgorithm::tuneIndependently(), CalibrationScanAlgorithm::tuneSimultaneously(), and tuneSimultaneously_.

                                                       {
 
   // Clear map holding analysis objects
   Analyses::iterator ianal;
   for ( ianal = data().begin(); ianal != data().end(); ianal++ ) {
     if ( ianal->second ) { delete ianal->second; }
   }
   data().clear();
   
   // Iterate through map containing vectors of profile histograms
   HistosMap::const_iterator iter = histos().begin();
 
   // One entry for each LLD channel --> differnt thousand entries
   for ( ; iter != histos().end(); iter++ ) {
     if ( iter->second.empty() ) {
       edm::LogWarning(mlDqmClient_)
        << "[CalibrationHistograms::" << __func__ << "]"
        << " Zero collation histograms found!";
       continue;
     }
     
     // Retrieve pointers to 1D histos for this FED channel --> all strips in the fiber = 256
     vector<TH1*> profs;
     Histos::const_iterator ihis = iter->second.begin();
     for ( ; ihis != iter->second.end(); ihis++ ) {
       TH1F* prof = ExtractTObject<TH1F>().extract( (*ihis)->me_ );
       if ( prof ) { profs.push_back(prof); }
     } 
     
     // Perform histo analysis 
     bool isdeconv = false;
     if(task() == sistrip::CALIBRATION_DECO or task() == sistrip::CALIBRATION_SCAN_DECO)
       isdeconv = true;
 
     if(task() == sistrip::CALIBRATION_SCAN or task() == sistrip::CALIBRATION_SCAN_DECO){
       CalibrationScanAnalysis* anal =  new CalibrationScanAnalysis(iter->first,isdeconv);
       CalibrationScanAlgorithm algo( this->pset(), anal );
       algo.analysis( profs );    
       data()[iter->first] = anal;     
 
       // tune the parameters for this a given target
       for(int iapv = 0; iapv < 2; iapv++){
     if(tuneSimultaneously_)
       algo.tuneSimultaneously(iapv,targetRiseTime_,targetDecayTime_);
     else
       algo.tuneIndependently(iapv,targetRiseTime_,targetDecayTime_);
     algo.fillTunedObservables(iapv);
       }      
     }
     else{
       CalibrationAnalysis* anal = new CalibrationAnalysis( iter->first,isdeconv);
       CalibrationAlgorithm algo( this->pset(), anal );
       algo.analysis( profs );          
       data()[iter->first] = anal;     
     }    
   } 
 }

void CalibrationHistograms::printAnalyses ( )

overridevirtual

Reimplemented from CommissioningHistograms.

Definition at line 115 of file CalibrationHistograms.cc.

References CommissioningHistograms::data(), LogTrace, and sistrip::mlDqmClient_.

                                           {
   Analyses::iterator ianal = data().begin();
   Analyses::iterator janal = data().end();
   for ( ; ianal != janal; ++ianal ) {
     if ( ianal->second ) {
       std::stringstream ss;
       ianal->second->print( ss,0 );
       ianal->second->print( ss,1 );
       if ( ianal->second->isValid() ) { LogTrace(mlDqmClient_) << ss.str();
       } else { edm::LogWarning(mlDqmClient_) << ss.str(); }
     }
   }
 }

void CalibrationHistograms::save	(	std::string &	filename,
		uint32_t	run_number = `0`,
		std::string	partitionName = `""`
	)

Definition at line 132 of file CalibrationHistograms.cc.

References CommissioningHistograms::bei(), SiStripFecKey::ccuAddr(), SiStripFecKey::ccuChan(), sistrip::collate_, ws_sso_content_reader::content, relmon_rootfiles_spy::contents, CommissioningHistograms::data(), CalibrationScanAnalysis::decayTimeVsISHAVsVFS(), CalibrationScanAnalysis::decayTimeVsVFS(), dir, createBeamHaloJobs::directory, sistrip::dqmClientFileName_, SiStripFecKey::fecCrate(), CommissioningAnalysis::fecKey(), SiStripFecKey::fecRing(), SiStripFecKey::fecSlot(), LogTrace, sistrip::mlDqmClient_, download_sqlite_cfg::outputFile, callgraph::path, CalibrationScanAnalysis::riseTimeVsISHA(), CalibrationScanAnalysis::riseTimeVsISHAVsVFS(), and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by Vispa.Main.TabController.TabController::allowClose(), SiStripCommissioningOfflineClient::beginRun(), and Vispa.Main.TabController.TabController::checkModificationTimestamp().

                                            {
 
   // Construct path and filename                                                                                                                                                                       
   std::stringstream ss;
   if ( !path.empty() ) { // create with a specific outputName                                                                                                                                        
     ss << path;
     if ( ss.str().find(".root") == std::string::npos ) { ss << ".root"; }
 
   } else {
 
     // Retrieve SCRATCH directory                                                                                                                                                                    
     std::string scratch = "SCRATCH";
     std::string dir = "";
     if ( getenv(scratch.c_str()) != nullptr ) {
       dir = getenv(scratch.c_str());
     }
 
     // Add directory path                                                                                                                                                                             
     if ( !dir.empty() ) { ss << dir << "/"; }
     else { ss << "/tmp/"; }
 
     // Add filename with run number and ".root" extension                                                                                                                                          
     if(partitionName.empty())
       ss << sistrip::dqmClientFileName_ << "_"
          << std::setfill('0') << std::setw(8) << run_number
          << ".root";
     else
       ss << sistrip::dqmClientFileName_ << "_" << partitionName << "_"
          << std::setfill('0') << std::setw(8) << run_number
          << ".root";
   }
   
   // Save file with appropriate filename                                                                                                                                                               
   LogTrace(mlDqmClient_)
     << "[CommissioningHistograms::" << __func__ << "]"
     << " Saving histograms to root file"
     << " (This may take some time!)";
   path = ss.str();
   bei()->save( path, sistrip::collate_ );
   edm::LogVerbatim(mlDqmClient_)
     << "[CommissioningHistograms::" << __func__ << "]"
     << " Saved histograms to root file \""
     << ss.str() << "\"!";
   
   // In case of calibration-scan, add also the TGraphs                                                                                                                                                
   // re-open the file                                                                                                                                                                                
   TFile* outputFile = TFile::Open(path.c_str(),"UPDATE");
   outputFile->cd();
   
   // get all sub-dirs            
   std::vector<std::string> contents;
   bei()->getContents( contents );
 
   TMultiGraph* graph_isha     = new TMultiGraph("riseTime_vs_isha","");
   TMultiGraph* graph_vfs      = new TMultiGraph("decayTime_vs_vfs","");
   
   bool save_graph_isha = false;
   bool save_graph_vfs  = false;
   
   // loop on the analysis objects which are storing all relevant results                                                                                                                             
   Analyses::iterator ianal = data().begin();
   Analyses::iterator janal = data().end();
   for ( ; ianal != janal; ++ianal ) {
     if ( ianal->second ) {
 
       CalibrationScanAnalysis* anal = dynamic_cast<CalibrationScanAnalysis*> (ianal->second);
       SiStripFecKey feckey = anal->fecKey();
 
       TString directory;
       for(auto content : contents){
     std::vector<std::string> tokens;
     std::string token;
     std::istringstream tokenStream(content);
     while (std::getline(tokenStream, token,':')){
       tokens.push_back(token);
     }
     directory = Form("%s",tokens.at(0).c_str());
     if(directory.Contains(Form("FecCrate%d",feckey.fecCrate())) and
        directory.Contains(Form("FecRing%d",feckey.fecRing())) and
        directory.Contains(Form("FecSlot%d",feckey.fecSlot())) and
        directory.Contains(Form("CcuAddr%d",feckey.ccuAddr())) and
        directory.Contains(Form("CcuChan%d",feckey.ccuChan())))
       break;    
       }
       
       outputFile->cd("DQMData/"+directory);
 
       for(size_t igraph = 0 ; igraph < anal->decayTimeVsVFS().size(); igraph++){    
     graph_vfs->Add(anal->decayTimeVsVFS()[igraph]);
     anal->decayTimeVsVFS()[igraph]->Write();
     save_graph_vfs = true;
       }
 
       for(size_t igraph = 0 ; igraph < anal->riseTimeVsISHA().size(); igraph++){    
     graph_isha->Add(anal->riseTimeVsISHA()[igraph]);
     anal->riseTimeVsISHA()[igraph]->Write();
     save_graph_isha = true;
       }      
 
       for(size_t igraph = 0 ; igraph < anal->riseTimeVsISHAVsVFS().size(); igraph++)
     anal->riseTimeVsISHAVsVFS()[igraph]->Write();
 
 
       for(size_t igraph = 0 ; igraph < anal->decayTimeVsISHAVsVFS().size(); igraph++)
     anal->decayTimeVsISHAVsVFS()[igraph]->Write();      
 
       outputFile->cd();
     }
   }
 
   outputFile->cd();
   outputFile->cd("DQMData/Collate/SiStrip/ControlView");
 
 
   if(save_graph_isha)
     graph_isha->Write("riseTime_vs_isha");
   if(save_graph_vfs)
     graph_vfs->Write("decayTime_vs_vfs");
 
   outputFile->Close();
 } 