CalibrationHistograms Class Reference

#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 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.

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_.

    : 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_ = make_unique<CalibrationScanSummaryFactory>();
  else
    factory_ = make_unique<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::~CalibrationHistograms ( )

override

Definition at line 46 of file CalibrationHistograms.cc.

References LogTrace, and sistrip::mlDqmClient_.

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

Member Function Documentation

histoAnalysis()

void CalibrationHistograms::histoAnalysis ( bool  debug )

overridevirtual

Reimplemented from CommissioningHistograms.

Definition at line 53 of file CalibrationHistograms.cc.

References conversion_template_cfg::anal, sistrip::CALIBRATION_DECO, sistrip::CALIBRATION_SCAN, sistrip::CALIBRATION_SCAN_DECO, CommissioningHistograms::data(), extract(), CommissioningHistograms::histos(), sistrip::mlDqmClient_, or, CommissioningHistograms::pset(), targetDecayTime_, targetRiseTime_, CommissioningHistograms::task(), 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;
    }
  }
}

printAnalyses()

void CalibrationHistograms::printAnalyses ( )

overridevirtual

Reimplemented from CommissioningHistograms.

Definition at line 114 of file CalibrationHistograms.cc.

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

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

save()

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

Definition at line 133 of file CalibrationHistograms.cc.

References conversion_template_cfg::anal, CommissioningHistograms::bei(), SiStripFecKey::ccuAddr(), SiStripFecKey::ccuChan(), sistrip::collate_, relmon_rootfiles_spy::contents, CommissioningHistograms::data(), DeadROC_duringRun::dir, createBeamHaloJobs::directory, sistrip::dqmClientFileName_, SiStripFecKey::fecCrate(), SiStripFecKey::fecRing(), SiStripFecKey::fecSlot(), dqm::implementation::IGetter::getAllContents(), LogTrace, hlt_dqm_clientPB-live_cfg::me, sistrip::mlDqmClient_, download_sqlite_cfg::outputFile, FrameHeaderEvents_cfi::partitionName, castor_dqm_sourceclient_file_cfg::path, contentValuesCheck::run_number, dqm::implementation::DQMStore::save(), contentValuesCheck::ss, and AlCaHLTBitMon_QueryRunRegistry::string.

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

Member Data Documentation

targetDecayTime_

float CalibrationHistograms::targetDecayTime_

private

Definition at line 21 of file CalibrationHistograms.h.

Referenced by CalibrationHistograms(), and histoAnalysis().

targetRiseTime_

float CalibrationHistograms::targetRiseTime_

private

Definition at line 20 of file CalibrationHistograms.h.

Referenced by CalibrationHistograms(), and histoAnalysis().

tuneSimultaneously_

bool CalibrationHistograms::tuneSimultaneously_

private

Definition at line 22 of file CalibrationHistograms.h.

Referenced by CalibrationHistograms(), and histoAnalysis().