simBeamSpotHLLHCPI::DisplayParametersDiff< PayloadType, nIOVs, ntags > Class Template Reference

#include <BeamSpotPayloadInspectorHelper.h>

Inheritance diagram for simBeamSpotHLLHCPI::DisplayParametersDiff< PayloadType, nIOVs, ntags >:

Public Member Functions
	DisplayParametersDiff ()
bool	fill () override
Public Member Functions inherited from cond::payloadInspector::PlotImage< PayloadType, nIOVs, ntags >
std::shared_ptr< PayloadType >	fetchPayload (const cond::Hash &payloadHash)
	PlotImage (const std::string &title)
std::string	serializeData () override
Public Member Functions inherited from cond::payloadInspector::PlotImpl< IOV_M, NTAGS >
	PlotImpl (const std::string &type, const std::string &title)
std::string	processData () override
	~PlotImpl () override=default
Public Member Functions inherited from cond::payloadInspector::PlotBase
void	addInputParam (const std::string &paramName)
cond::persistency::Session	dbSession ()
template<typename PayloadType >
std::shared_ptr< PayloadType >	fetchPayload (const cond::Hash &payloadHash)
template<int index>
TagReference	getTag ()
cond::Tag_t	getTagInfo (const std::string &tag)
virtual void	init ()
const std::map< std::string, std::string > &	inputParamValues () const
bool	isSingleIov () const
bool	isTwoTags () const
unsigned int	ntags () const
std::string	payloadType () const
	PlotBase ()
std::string	title () const
std::string	type () const
virtual	~PlotBase ()=default

Protected Attributes
std::shared_ptr< PayloadType >	f_payload
std::shared_ptr< PayloadType >	l_payload
Protected Attributes inherited from cond::payloadInspector::PlotImage< PayloadType, nIOVs, ntags >
std::string	m_imageFileName
Protected Attributes inherited from cond::payloadInspector::PlotBase
std::set< std::string >	m_inputParams
std::map< std::string, std::string >	m_inputParamValues
PlotAnnotations	m_plotAnnotations
std::vector< std::pair< cond::Time_t, cond::Time_t > >	m_tagBoundaries
std::vector< std::vector< std::tuple< cond::Time_t, cond::Hash > > >	m_tagIovs
std::vector< std::string >	m_tagNames

Additional Inherited Members
Public Types inherited from cond::payloadInspector::PlotImage< PayloadType, nIOVs, ntags >
typedef PlotImpl< IOV_M, NTAGS >	Base

Detailed Description

template<class PayloadType, cond::payloadInspector::IOVMultiplicity nIOVs, int ntags>
class simBeamSpotHLLHCPI::DisplayParametersDiff< PayloadType, nIOVs, ntags >

Definition at line 1347 of file BeamSpotPayloadInspectorHelper.h.

Constructor & Destructor Documentation

DisplayParametersDiff()

template<class PayloadType , cond::payloadInspector::IOVMultiplicity nIOVs, int ntags>

simBeamSpotHLLHCPI::DisplayParametersDiff< PayloadType, nIOVs, ntags >::DisplayParametersDiff ( )

inline

Definition at line 1349 of file BeamSpotPayloadInspectorHelper.h.

        : cond::payloadInspector::PlotImage<PayloadType, nIOVs, ntags>(
              "Display of SimBeamSpotHLLHC parameters differences") {}

Member Function Documentation

fill()

template<class PayloadType , cond::payloadInspector::IOVMultiplicity nIOVs, int ntags>

bool simBeamSpotHLLHCPI::DisplayParametersDiff< PayloadType, nIOVs, ntags >::fill ( )

inlineoverridevirtual

Implements cond::payloadInspector::PlotImpl< IOV_M, NTAGS >.

Definition at line 1353 of file BeamSpotPayloadInspectorHelper.h.

References cms::cuda::assert(), svgfig::canvas(), simBeamSpotHLLHCPI::crabbingAngleCrossing, simBeamSpotHLLHCPI::crabbingAngleSeparation, simBeamSpotHLLHCPI::SimBSHLLHCParamsHelper< PayloadType >::diffCentralValues(), END_OF_TYPES, simBeamSpotHLLHCPI::DisplayParametersDiff< PayloadType, nIOVs, ntags >::f_payload, cond::payloadInspector::PlotImage< PayloadType, nIOVs, ntags >::fetchPayload(), MillePedeFileConverter_cfg::fileName, custom_jme_cff::foo, simBeamSpotHLLHCPI::getStringFromParamEnum(), simBeamSpotHLLHCPI::DisplayParametersDiff< PayloadType, nIOVs, ntags >::l_payload, cond::payloadInspector::PlotImage< PayloadType, nIOVs, ntags >::m_imageFileName, cond::payloadInspector::PlotBase::m_plotAnnotations, SiStripPI::max, simBeamSpotHLLHCPI::meanX, SiStripPI::min, Skims_PA_cff::name, cond::payloadInspector::PlotAnnotations::ntags, contentValuesFiles::number, simBeamSpotHLLHCPI::SimBSHLLHCParamsHelper< PayloadType >::printDebug(), AlCaHLTBitMon_QueryRunRegistry::string, simBeamSpotHLLHCPI::timeOffset, to_string(), beamSpotPI::unpack(), and photonAnalyzer_cfi::yBin.

                         {
      // trick to deal with the multi-ioved tag and two tag case at the same time
      auto theIOVs = cond::payloadInspector::PlotBase::getTag<0>().iovs;
      auto f_tagname = cond::payloadInspector::PlotBase::getTag<0>().name;
      std::string l_tagname = "";
      auto firstiov = theIOVs.front();
      std::tuple<cond::Time_t, cond::Hash> lastiov;

      // we don't support (yet) comparison with more than 2 tags
      assert(this->m_plotAnnotations.ntags < 3);

      if (this->m_plotAnnotations.ntags == 2) {
        auto tag2iovs = cond::payloadInspector::PlotBase::getTag<1>().iovs;
        l_tagname = cond::payloadInspector::PlotBase::getTag<1>().name;
        lastiov = tag2iovs.front();
      } else {
        lastiov = theIOVs.back();
      }

      l_payload = this->fetchPayload(std::get<1>(lastiov));
      f_payload = this->fetchPayload(std::get<1>(firstiov));

      std::string lastIOVsince = std::to_string(std::get<0>(lastiov));
      std::string firstIOVsince = std::to_string(std::get<0>(firstiov));

      TCanvas canvas("Sim Beam Spot HL-LHC Parameters Difference Summary",
                     "Sim Beam Spot HL-LHC Parameters Difference summary",
                     1500,
                     1000);
      canvas.Divide(2, 1);

      // Histograms with parameters differences
      // for the "text"-filled histogram
      auto h2_SimBSHLLHCParameters = std::make_unique<TH2F>(
          "Parameters", "", 1, 0.0, 1.0, parameters::crabbingAngleCrossing, 0, parameters::crabbingAngleCrossing);
      h2_SimBSHLLHCParameters->SetStats(false);
      h2_SimBSHLLHCParameters->GetXaxis()->SetBinLabel(1, "Value");
      h2_SimBSHLLHCParameters->GetXaxis()->LabelsOption("h");
      h2_SimBSHLLHCParameters->GetYaxis()->SetLabelSize(0.05);
      h2_SimBSHLLHCParameters->GetXaxis()->SetLabelSize(0.05);
      h2_SimBSHLLHCParameters->SetMarkerSize(1.5);

      const auto span = parameters::END_OF_TYPES - parameters::crabbingAngleSeparation;
      auto h2_SimBSHLLHCParameters_Extras =
          std::make_unique<TH2F>("Parameters Extras", "", 1, 0.0, 1.0, span, 0, static_cast<float>(span));
      h2_SimBSHLLHCParameters_Extras->SetStats(false);
      h2_SimBSHLLHCParameters_Extras->GetXaxis()->SetBinLabel(1, "Value");
      h2_SimBSHLLHCParameters_Extras->GetXaxis()->LabelsOption("h");
      h2_SimBSHLLHCParameters_Extras->GetYaxis()->SetLabelSize(0.05);
      h2_SimBSHLLHCParameters_Extras->GetXaxis()->SetLabelSize(0.05);
      h2_SimBSHLLHCParameters_Extras->SetMarkerSize(1.5);

      // prepare the arrays to fill the histogram
      simBeamSpotHLLHCPI::SimBSHLLHCParamsHelper fBS(f_payload);
      simBeamSpotHLLHCPI::SimBSHLLHCParamsHelper lBS(l_payload);

#ifdef MM_DEBUG
      std::stringstream ss1, ss2;
      edm::LogPrint("") << "**** first payload";
      fBS.printDebug(ss1);
      edm::LogPrint("") << ss1.str();
      edm::LogPrint("") << "**** last payload";
      lBS.printDebug(ss2);
      edm::LogPrint("") << ss2.str();
#endif

      const auto diffPars = fBS.diffCentralValues(lBS);

      // fill "text"-filled histograms
      unsigned int yBin = 8;
      for (int foo = parameters::meanX; foo <= parameters::crabbingAngleCrossing; foo++) {
        parameters param = static_cast<parameters>(foo);
        std::string theLabel = simBeamSpotHLLHCPI::getStringFromParamEnum(param, true /*use units*/);
        h2_SimBSHLLHCParameters->GetYaxis()->SetBinLabel(yBin, theLabel.c_str());
        h2_SimBSHLLHCParameters->SetBinContent(1, yBin, diffPars[foo]); /* profiting of the parameters enum indexing */
        yBin--;
      }

      unsigned int yBin_Extras = span;
      for (int foo = parameters::crabbingAngleSeparation; foo <= parameters::timeOffset; foo++) {
        parameters param = static_cast<parameters>(foo);
        std::string theLabel = simBeamSpotHLLHCPI::getStringFromParamEnum(param, true /*use units*/);
        h2_SimBSHLLHCParameters_Extras->GetYaxis()->SetBinLabel(yBin_Extras, theLabel.c_str());
        h2_SimBSHLLHCParameters_Extras->SetBinContent(
            1, yBin_Extras, diffPars[foo]); /* profiting of the parameters enum indexing */
        yBin_Extras--;
      }

      // for the "colz"-filled histogram (clone from the text-based one)
      auto h2_SimBSHLLHCShadow = (TH2F*)(h2_SimBSHLLHCParameters->Clone("shadow"));
      h2_SimBSHLLHCShadow->GetZaxis()->SetTitle("#Delta Parameter(payload A - payload B)");
      h2_SimBSHLLHCShadow->GetZaxis()->CenterTitle();
      h2_SimBSHLLHCShadow->GetZaxis()->SetTitleOffset(1.5);

      auto h2_SimBSHLLHCShadow_Extras = (TH2F*)(h2_SimBSHLLHCParameters_Extras->Clone("shadow"));
      h2_SimBSHLLHCShadow_Extras->GetZaxis()->SetTitle("#Delta Parameter(payload A - payload B)");
      h2_SimBSHLLHCShadow_Extras->GetZaxis()->CenterTitle();
      h2_SimBSHLLHCShadow_Extras->GetZaxis()->SetTitleOffset(1.5);

      // this is the fine gradient palette (blue to red)
      double max = std::max(h2_SimBSHLLHCShadow->GetMaximum(), h2_SimBSHLLHCShadow_Extras->GetMaximum());
      double min = std::min(h2_SimBSHLLHCShadow->GetMinimum(), h2_SimBSHLLHCShadow_Extras->GetMinimum());
      double val_white = 0.;
      double per_white = (max != min) ? ((val_white - min) / (max - min)) : 0.5;

      const int number = 3;
      double Red[number] = {0., 1., 1.};
      double Green[number] = {0., 1., 0.};
      double Blue[number] = {1., 1., 0.};
      double Stops[number] = {0., per_white, 1.};
      int nb = 256;
      h2_SimBSHLLHCShadow->SetContour(nb);
      h2_SimBSHLLHCShadow_Extras->SetContour(nb);
      TColor::CreateGradientColorTable(number, Stops, Red, Green, Blue, nb);

      // Fill canvas (1,1) with first 9 parameters differences
      canvas.cd(1);
      canvas.cd(1)->SetTopMargin(0.10);
      canvas.cd(1)->SetBottomMargin(0.06);
      canvas.cd(1)->SetLeftMargin(0.35);
      canvas.cd(1)->SetRightMargin(0.25);
      canvas.cd(1)->Modified();
      canvas.cd(1)->SetGrid();
      h2_SimBSHLLHCShadow->Draw("colz");
      h2_SimBSHLLHCParameters->Draw("TEXTsame");

      auto ltx = TLatex();
      ltx.SetTextFont(62);
      ltx.SetTextSize(0.025);
      ltx.SetTextAlign(11);

      // compute the (run,LS) pairs
      auto l_runLS = beamSpotPI::unpack(std::get<0>(lastiov));
      std::string l_runLSs = "(" + std::to_string(l_runLS.first) + "," + std::to_string(l_runLS.second) + ")";
      auto f_runLS = beamSpotPI::unpack(std::get<0>(firstiov));
      std::string f_runLSs = "(" + std::to_string(f_runLS.first) + "," + std::to_string(f_runLS.second) + ")";

      if (this->m_plotAnnotations.ntags == 2) {
        ltx.DrawLatexNDC(
            gPad->GetLeftMargin(),
            1 - gPad->GetTopMargin() + 0.025,
            (fmt::sprintf(
                 "#splitline{A = #color[4]{%s}: %s}{B = #color[4]{%s}: %s}", f_tagname, f_runLSs, l_tagname, l_runLSs))
                .c_str());
      } else {
        ltx.DrawLatexNDC(
            gPad->GetLeftMargin(),
            1 - gPad->GetTopMargin() + 0.025,
            (fmt::sprintf("#splitline{#color[4]{%s}}{A = %s | B = %s}", f_tagname, l_runLSs, f_runLSs)).c_str());
      }

      // Fill canvas (2,1) with first 9 parameters differences
      canvas.cd(2);
      canvas.cd(2)->SetTopMargin(0.10);
      canvas.cd(2)->SetBottomMargin(0.06);
      canvas.cd(2)->SetLeftMargin(0.35);
      canvas.cd(2)->SetRightMargin(0.20);
      canvas.cd(2)->Modified();
      canvas.cd(2)->SetGrid();
      h2_SimBSHLLHCShadow_Extras->Draw("colz");
      h2_SimBSHLLHCParameters_Extras->Draw("TEXTsame");

      if (this->m_plotAnnotations.ntags == 2) {
        ltx.DrawLatexNDC(
            gPad->GetLeftMargin(),
            1 - gPad->GetTopMargin() + 0.025,
            (fmt::sprintf(
                 "#splitline{A = #color[4]{%s}: %s}{B = #color[4]{%s}: %s}", f_tagname, f_runLSs, l_tagname, l_runLSs))
                .c_str());
      } else {
        ltx.DrawLatexNDC(
            gPad->GetLeftMargin(),
            1 - gPad->GetTopMargin() + 0.025,
            (fmt::sprintf("#splitline{#color[4]{%s}}{A = %s | B = %s}", f_tagname, l_runLSs, f_runLSs)).c_str());
      }

      std::string fileName(this->m_imageFileName);
      canvas.SaveAs(fileName.c_str());

      return true;
    }

Member Data Documentation

f_payload

template<class PayloadType , cond::payloadInspector::IOVMultiplicity nIOVs, int ntags>

std::shared_ptr<PayloadType> simBeamSpotHLLHCPI::DisplayParametersDiff< PayloadType, nIOVs, ntags >::f_payload

protected

Definition at line 1536 of file BeamSpotPayloadInspectorHelper.h.

Referenced by simBeamSpotHLLHCPI::DisplayParametersDiff< PayloadType, nIOVs, ntags >::fill().

l_payload

template<class PayloadType , cond::payloadInspector::IOVMultiplicity nIOVs, int ntags>

std::shared_ptr<PayloadType> simBeamSpotHLLHCPI::DisplayParametersDiff< PayloadType, nIOVs, ntags >::l_payload

protected

Definition at line 1537 of file BeamSpotPayloadInspectorHelper.h.

Referenced by simBeamSpotHLLHCPI::DisplayParametersDiff< PayloadType, nIOVs, ntags >::fill().