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

	CommissioningHistograms (const edm::ParameterSet &pset, DQMStore *const, const sistrip::RunType &)

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 dqm::harvesting::DQMStore	DQMStore

typedef SummaryPlotFactory< CommissioningAnalysis * >	Factory

typedef std::map< uint32_t, uint32_t >	FedToFecMap

typedef std::vector< Histo * >	Histos

typedef std::map< uint32_t, Histos >	HistosMap

typedef dqm::harvesting::MonitorElement	MonitorElement

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 7 of file CalibrationHistograms.h.

Constructor & Destructor Documentation

◆ CalibrationHistograms()

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

Definition at line 26 of file CalibrationHistograms.cc.

     : 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;
 }

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

◆ ~CalibrationHistograms()

CalibrationHistograms::~CalibrationHistograms ( )

override

Definition at line 46 of file CalibrationHistograms.cc.

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

References LogTrace, and sistrip::mlDqmClient_.

Member Function Documentation

◆ histoAnalysis()

void CalibrationHistograms::histoAnalysis ( bool debug )

overridevirtual

Reimplemented from CommissioningHistograms.

Definition at line 53 of file CalibrationHistograms.cc.

                                                     {
   // 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;
     }
   }
 }

References conversion_template_cfg::anal, begin, sistrip::CALIBRATION_DECO, sistrip::CALIBRATION_SCAN, sistrip::CALIBRATION_SCAN_DECO, CommissioningHistograms::data(), ExtractTObject< T >::extract(), CommissioningHistograms::histos(), sistrip::mlDqmClient_, or, CommissioningHistograms::pset(), targetDecayTime_, targetRiseTime_, CommissioningHistograms::task(), and tuneSimultaneously_.

◆ printAnalyses()

void CalibrationHistograms::printAnalyses ( )

overridevirtual

Reimplemented from CommissioningHistograms.

Definition at line 114 of file CalibrationHistograms.cc.

                                           {
   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();
       }
     }
   }
 }

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

◆ save()

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

Definition at line 133 of file CalibrationHistograms.cc.

                                                                                               {
   // 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 (std::getenv(scratch.c_str()) != nullptr) {
       dir = std::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();
  
   auto contents = bei()->getAllContents("");
  
   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 me : contents) {
         directory = me->getPathname();
         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();
 }