EcalTPGPhysicsConst_PayloadInspector.cc

Go to the documentation of this file.

#include "CondCore/Utilities/interface/PayloadInspectorModule.h"
#include "CondCore/Utilities/interface/PayloadInspector.h"
#include "CondCore/CondDB/interface/Time.h"
#include "DataFormats/EcalDetId/interface/EBDetId.h"
#include "DataFormats/EcalDetId/interface/EEDetId.h"
#include "CondCore/EcalPlugins/plugins/EcalDrawUtils.h"

// the data format of the condition to be inspected
#include "CondFormats/EcalObjects/interface/EcalTPGPhysicsConst.h"

#include "TH2F.h"  // a 2-D histogram with four bytes per cell (float)
#include "TCanvas.h"
#include "TLine.h"
#include "TStyle.h"
#include "TLatex.h"  //write mathematical equations.
#include "TPave.h"
#include "TPaveStats.h"
#include <string>
#include <fstream>

namespace {
  /*****************************************
 2d plot of Ecal TPG Physics Const of 1 IOV
 ******************************************/
  class EcalTPGPhysicsConstPlot : public cond::payloadInspector::PlotImage<EcalTPGPhysicsConst> {
  public:
    EcalTPGPhysicsConstPlot()
        : cond::payloadInspector::PlotImage<EcalTPGPhysicsConst>("ECAL TPG Physics Constant - map ") {
      setSingleIov(true);
    }

    bool fill(const std::vector<std::tuple<cond::Time_t, cond::Hash> >& iovs) override {
      auto iov = iovs.front();
      std::shared_ptr<EcalTPGPhysicsConst> payload = fetchPayload(std::get<1>(iov));
      unsigned int run = std::get<0>(iov);
      TH2F* align;
      int NbRows;

      if (payload.get()) {
        EcalTPGPhysicsConstMap map = (*payload).getMap();
        NbRows = map.size();

        align = new TH2F("TPGPhysicsConstant",
                         "mapKey   EtSat         ttf_threshold_Low         ttf_threshold_High         FG_lowThreshold  "
                         "       FG_highThreshold         FG_lowRatio         FG_highRatio",
                         8,
                         0,
                         8,
                         NbRows,
                         0,
                         NbRows);

        double row = NbRows - 0.5;
        for (std::map<uint32_t, EcalTPGPhysicsConst::Item>::const_iterator it = map.begin(); it != map.end(); it++) {
          uint32_t mapKey = it->first;
          EcalTPGPhysicsConst::Item item = it->second;

          align->Fill(0.5, row, mapKey);
          align->Fill(1.5, row, item.EtSat);
          align->Fill(2.5, row, item.ttf_threshold_Low);
          align->Fill(3.5, row, item.ttf_threshold_High);
          align->Fill(4.5, row, item.FG_lowThreshold);
          align->Fill(5.5, row, item.FG_highThreshold);
          align->Fill(6.5, row, item.FG_lowRatio);
          align->Fill(7.5, row, item.FG_highRatio);

          row = row - 1.;
        }

      } else
        return false;

      gStyle->SetPalette(1);
      gStyle->SetOptStat(0);
      TCanvas canvas("CC map", "CC map", 1000, 1000);
      TLatex t1;
      t1.SetNDC();
      t1.SetTextAlign(26);
      t1.SetTextSize(0.05);
      t1.SetTextColor(2);
      t1.DrawLatex(0.5, 0.96, Form("ECAL TPG Physics Constant, IOV %i", run));

      TPad* pad = new TPad("pad", "pad", 0.0, 0.0, 1.0, 0.94);
      pad->Draw();
      pad->cd();
      align->Draw("TEXT");

      drawTable(NbRows, 8);

      align->GetXaxis()->SetTickLength(0.);
      align->GetXaxis()->SetLabelSize(0.);
      align->GetYaxis()->SetTickLength(0.);
      align->GetYaxis()->SetLabelSize(0.);

      std::string ImageName(m_imageFileName);
      canvas.SaveAs(ImageName.c_str());

      return true;
    }
  };

}  // namespace
// Register the classes as boost python plugin
PAYLOAD_INSPECTOR_MODULE(EcalTPGPhysicsConst) { PAYLOAD_INSPECTOR_CLASS(EcalTPGPhysicsConstPlot); }