JME::JetScaleFactorSummary Class Reference

Inheritance diagram for JME::JetScaleFactorSummary:

Public Member Functions
bool	fill () override
	JetScaleFactorSummary ()
Public Member Functions inherited from cond::payloadInspector::PlotImage< JetResolutionObject, SINGLE_IOV >
std::shared_ptr< JetResolutionObject >	fetchPayload (const cond::Hash &payloadHash)
	PlotImage (const std::string &title)
std::string	serializeData () override
Public Member Functions inherited from cond::payloadInspector::PlotImpl< IOV_M, 0 >
	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

Additional Inherited Members
Public Types inherited from cond::payloadInspector::PlotImage< JetResolutionObject, SINGLE_IOV >
typedef PlotImpl< IOV_M, 0 >	Base
Protected Attributes inherited from cond::payloadInspector::PlotImage< JetResolutionObject, SINGLE_IOV >
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

Detailed Description

Definition at line 697 of file JetResolution_PayloadInspector.cc.

Constructor & Destructor Documentation

JetScaleFactorSummary()

JME::JetScaleFactorSummary::JetScaleFactorSummary ( )

inline

Definition at line 699 of file JetResolution_PayloadInspector.cc.

References cond::payloadInspector::PlotBase::addInputParam().

        : cond::payloadInspector::PlotImage<JetResolutionObject, SINGLE_IOV>("Jet ScaleFactor Summary") {
      cond::payloadInspector::PlotBase::addInputParam("Jet_Pt");
      cond::payloadInspector::PlotBase::addInputParam("Jet_Eta");
    }

Member Function Documentation

fill()

bool JME::JetScaleFactorSummary::fill ( )

inlineoverridevirtual

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

Definition at line 705 of file JetResolution_PayloadInspector.cc.

References svgfig::canvas(), cond::payloadInspector::PlotBase::inputParamValues(), MAX_ETA, MAX_PT, MIN_ETA, MIN_PT, NBIN_ETA, NBIN_PT, jetsAK4_Puppi_cff::payload, AlCaHarvesting_cff::record, writedatasetfile::run, AlCaHLTBitMon_QueryRunRegistry::string, makeGlobalPositionRcd_cfg::tag, createPayload::tagname, and runGCPTkAlMap::title.

                         {
      double par_Pt = 100.;
      double par_Eta = 1.;

      auto paramValues = cond::payloadInspector::PlotBase::inputParamValues();
      auto ip = paramValues.find("Jet_Pt");
      if (ip != paramValues.end()) {
        par_Pt = std::stod(ip->second);
      }
      ip = paramValues.find("Jet_Eta");
      if (ip != paramValues.end()) {
        par_Eta = std::stod(ip->second);
      }

      TH1D* sf_eta_norm = new TH1D("Jet SF vs #eta NORM", "", NBIN_ETA, MIN_ETA, MAX_ETA);
      TH1D* sf_eta_down = new TH1D("Jet SF vs #eta DOWN", "", NBIN_ETA, MIN_ETA, MAX_ETA);
      TH1D* sf_eta_up = new TH1D("Jet SF vs #eta UP", "", NBIN_ETA, MIN_ETA, MAX_ETA);
      TH1D* sf_pt_norm = new TH1D("Jet SF vs p_T NORM", "", NBIN_PT, MIN_PT, MAX_PT);
      TH1D* sf_pt_down = new TH1D("Jet SF vs p_T DOWN", "", NBIN_PT, MIN_PT, MAX_PT);
      TH1D* sf_pt_up = new TH1D("Jet SF vs p_T UP", "", NBIN_PT, MIN_PT, MAX_PT);

      TLegend* leg_eta = new TLegend(0.26, 0.73, 0.935, 0.90);
      TLegend* leg_pt = new TLegend(0.26, 0.73, 0.935, 0.90);

      leg_eta->SetBorderSize(0);
      leg_eta->SetLineStyle(0);
      leg_eta->SetFillStyle(0);

      leg_eta->SetTextFont(42);
      leg_pt->SetBorderSize(0);
      leg_pt->SetLineStyle(0);
      leg_pt->SetFillStyle(0);

      auto tag = PlotBase::getTag<0>();
      auto iov = tag.iovs.front();
      std::shared_ptr<JetResolutionObject> payload = fetchPayload(std::get<1>(iov));
      unsigned int run = std::get<0>(iov);
      std::string tagname = tag.name;
      std::stringstream ss_tagname(tag.name);
      std::string stmp;

      std::string tag_ver;
      std::string tag_res;
      std::string tag_jet;

      getline(ss_tagname, stmp, '_');  // drop first
      getline(ss_tagname, stmp, '_');  // year
      tag_ver = stmp;
      getline(ss_tagname, stmp, '_');  // ver
      tag_ver += '_' + stmp;
      getline(ss_tagname, stmp, '_');  // cmssw
      tag_ver += '_' + stmp;
      getline(ss_tagname, stmp, '_');  // data/mc
      tag_ver += '_' + stmp;
      getline(ss_tagname, stmp, '_');  // bin
      tag_res = stmp;
      getline(ss_tagname, stmp, '_');  // jet algorithm
      tag_jet = stmp;

      bool is_2bin = false;

      if (payload.get()) {
        if (!payload->getRecords().empty() &&  // No formula for SF
            payload->getDefinition().getFormulaString().compare("") != 0)
          return false;

        is_2bin = false;
        for (const auto& record : payload->getRecords()) {
          if (record.getBinsRange().size() > 1)
            is_2bin = true;

          if (!record.getBinsRange().empty() && payload->getDefinition().getBinName(0) == "JetEta" &&
              record.getParametersValues().size() == 3) {  // norm, down, up

            for (size_t it = 0; it <= NBIN_ETA; it++) {
              double x_axis = (it + 0.5) * (MAX_ETA - MIN_ETA) / NBIN_ETA + MIN_ETA;
              if (((is_2bin == false) || (is_2bin == true && record.getBinsRange()[1].is_inside(par_Pt))) &&
                  record.getBinsRange()[0].is_inside(x_axis)) {
                sf_eta_norm->SetBinContent(it + 1, record.getParametersValues()[0]);
                sf_eta_down->SetBinContent(it + 1, record.getParametersValues()[1]);
                sf_eta_up->SetBinContent(it + 1, record.getParametersValues()[2]);
              }
            }
          }

          if (record.getBinsRange().size() > 1 && payload->getDefinition().getBinName(0) == "JetEta" &&
              record.getBinsRange()[0].is_inside(par_Eta) &&  // take jeteta=2.5
              payload->getDefinition().getBinName(1) == "JetPt" &&
              record.getParametersValues().size() == 3) {  // norm, down, up

            is_2bin = true;

            for (size_t it = 0; it <= NBIN_PT; it++) {
              double x_axis = (it + 0.5) * (MAX_PT - MIN_PT) / NBIN_PT + MIN_PT;
              if (record.getBinsRange()[1].is_inside(x_axis)) {
                sf_pt_norm->SetBinContent(it + 1, record.getParametersValues()[0]);
                sf_pt_down->SetBinContent(it + 1, record.getParametersValues()[1]);
                sf_pt_up->SetBinContent(it + 1, record.getParametersValues()[2]);
              }
            }
          }  // 2-bin
        }    // records

        gStyle->SetOptStat(0);
        gStyle->SetLabelFont(42, "XYZ");
        gStyle->SetLabelSize(0.05, "XYZ");
        gStyle->SetFrameLineWidth(3);

        std::string title = Form("Summary Run %i", run);
        TCanvas canvas("Jet ScaleFactor Summary", title.c_str(), 800, 1200);
        canvas.Divide(1, 2);

        canvas.cd(1);
        sf_eta_up->SetTitle("ScaleFactor vs. #eta");
        sf_eta_up->SetXTitle("#eta");
        sf_eta_up->SetYTitle("Scale Factor");
        sf_eta_up->SetLineStyle(7);
        sf_eta_up->SetLineWidth(3);
        sf_eta_up->SetFillColorAlpha(kGray, 0.5);
        sf_eta_up->SetMaximum(sf_eta_up->GetMaximum() * 1.25);
        sf_eta_up->SetMinimum(0.);
        sf_eta_up->Draw("][");

        sf_eta_down->SetLineStyle(7);
        sf_eta_down->SetLineWidth(3);
        sf_eta_down->SetFillColorAlpha(kWhite, 1);
        sf_eta_down->Draw("][ same");

        sf_eta_norm->SetLineStyle(1);
        sf_eta_norm->SetLineWidth(5);
        sf_eta_norm->SetFillColor(0);
        sf_eta_norm->Draw("][ same");
        sf_eta_norm->Draw("axis same");

        leg_eta->AddEntry(sf_eta_norm, (tag_ver + '_' + tag_jet).c_str(), "l");
        leg_eta->AddEntry((TObject*)nullptr, Form("JetPt=%.2f", par_Pt), "");
        leg_eta->Draw();

        if (is_2bin == true) {
          canvas.cd(2);
          sf_pt_up->SetTitle("ScaleFactor vs. p_{T}");
          sf_pt_up->SetXTitle("p_{T} [GeV]");
          sf_pt_up->SetYTitle("Scale Factor");
          sf_pt_up->SetLineStyle(7);
          sf_pt_up->SetLineWidth(3);
          sf_pt_up->SetFillColorAlpha(kGray, 0.5);
          sf_pt_up->SetMaximum(sf_pt_up->GetMaximum() * 1.25);
          sf_pt_up->SetMinimum(0.);
          sf_pt_up->Draw("][");

          sf_pt_down->SetLineStyle(7);
          sf_pt_down->SetLineWidth(3);
          sf_pt_down->SetFillColorAlpha(kWhite, 1);
          sf_pt_down->Draw("][ same");

          sf_pt_norm->SetLineStyle(1);
          sf_pt_norm->SetLineWidth(5);
          sf_pt_norm->SetFillColor(0);
          sf_pt_norm->Draw("][ same");
          sf_pt_norm->Draw("axis same");

          leg_pt->AddEntry(sf_pt_norm, (tag_ver + '_' + tag_jet).c_str(), "l");
          leg_pt->AddEntry((TObject*)nullptr, Form("JetEta=%.2f", par_Eta), "");
          leg_pt->Draw();
        }

        canvas.SaveAs(m_imageFileName.c_str());

        return true;
      } else  // no payload.get()
        return false;
    }  // fill