SiPixelPhase1EfficiencyExtras Class Reference

Inheritance diagram for SiPixelPhase1EfficiencyExtras:

Public Member Functions
	SiPixelPhase1EfficiencyExtras (const edm::ParameterSet &conf)
	~SiPixelPhase1EfficiencyExtras () override

Protected Member Functions
void	beginRun (edm::Run const &run, edm::EventSetup const &eSetup) override
void	dqmEndJob (DQMStore::IBooker &iBooker, DQMStore::IGetter &iGetter) override

Protected Attributes
const std::string	effFolderName_
const std::string	instLumiFolderName_
const std::string	vtxFolderName_

Detailed Description

Definition at line 35 of file SiPixelPhase1EfficiencyExtras.cc.

Constructor & Destructor Documentation

SiPixelPhase1EfficiencyExtras()

SiPixelPhase1EfficiencyExtras::SiPixelPhase1EfficiencyExtras ( const edm::ParameterSet &  conf )

explicit

Definition at line 50 of file SiPixelPhase1EfficiencyExtras.cc.

    : effFolderName_(iConfig.getParameter<std::string>("EffFolderName")),
      vtxFolderName_(iConfig.getParameter<std::string>("VtxFolderName")),
      instLumiFolderName_(iConfig.getParameter<std::string>("InstLumiFolderName")) {
  LogInfo("PixelDQM") << "SiPixelPhase1EfficiencyExtras::SiPixelPhase1EfficiencyExtras: Hello!" << endl;
}

~SiPixelPhase1EfficiencyExtras()

SiPixelPhase1EfficiencyExtras::~SiPixelPhase1EfficiencyExtras ( )

override

Definition at line 57 of file SiPixelPhase1EfficiencyExtras.cc.

                                                              {
  LogInfo("PixelDQM") << "SiPixelPhase1EfficiencyExtras::~SiPixelPhase1EfficiencyExtras: Destructor" << endl;
}

Member Function Documentation

beginRun()

void SiPixelPhase1EfficiencyExtras::beginRun ( edm::Run const &  run, edm::EventSetup const &  eSetup  )

overrideprotected

Definition at line 61 of file SiPixelPhase1EfficiencyExtras.cc.

{}

dqmEndJob()

void SiPixelPhase1EfficiencyExtras::dqmEndJob ( DQMStore::IBooker &  iBooker, DQMStore::IGetter &  iGetter  )

overrideprotected

Definition at line 66 of file SiPixelPhase1EfficiencyExtras.cc.

References dqm::implementation::IBooker::book2D(), effFolderName_, dqm::implementation::IGetter::get(), dqm::impl::MonitorElement::getBinContent(), dqm::impl::MonitorElement::getEntries(), dqm::impl::MonitorElement::getNbinsX(), dqm::legacy::MonitorElement::getTH2F(), dqm::legacy::MonitorElement::getTProfile(), instLumiFolderName_, createfilelist::int, dqm::impl::MonitorElement::setBinContent(), dqm::implementation::NavigatorBase::setCurrentFolder(), and vtxFolderName_.

                                                                                                {
  iBooker.setCurrentFolder(effFolderName_);

  //Get the existing histos
  MonitorElement* vtx_v_lumi = iGetter.get(vtxFolderName_ + "/NumberOfGoodPVtxVsLS_GenTk");

  MonitorElement* scalLumi_v_lumi = iGetter.get(instLumiFolderName_ + "/lumiVsLS");

  MonitorElement* eff_v_lumi_forward =
      iGetter.get(effFolderName_ + "/hitefficiency_per_Lumisection_per_PXDisk_PXForward");

  MonitorElement* eff_v_lumi_barrel =
      iGetter.get(effFolderName_ + "/hitefficiency_per_Lumisection_per_PXLayer_PXBarrel");

  //set up some booleans that will tell us which graphs to create
  bool createNvtx = true;
  bool createInstLumi = true;

  //check which of the MEs exist and respond appropriately
  if (!eff_v_lumi_forward) {
    edm::LogWarning("SiPixelPhase1EfficiencyExtras")
        << "no hitefficiency_per_Lumisection_per_PXDisk_PXForward ME is available in " << effFolderName_ << std::endl;
    return;
  }
  if (!eff_v_lumi_barrel) {
    edm::LogWarning("SiPixelPhase1EfficiencyExtras")
        << "no hitefficiency_per_Lumisection_per_PXLayer_PXBarrel ME is available in " << effFolderName_ << std::endl;
    return;
  }
  if (!vtx_v_lumi) {
    edm::LogWarning("SiPixelPhase1EfficiencyExtras")
        << "no NumberOfGoodPVtxVsLS_GenTK ME is available in " << vtxFolderName_ << std::endl;
    createNvtx = false;
  }
  if (!scalLumi_v_lumi) {
    edm::LogWarning("SiPixelPhase1EfficiencyExtras")
        << "no lumiVsLS ME is available in " << instLumiFolderName_ << std::endl;
    createInstLumi = false;
  }

  //If the existing MEs are empty, set the boolean to skip booking
  if (vtx_v_lumi && vtx_v_lumi->getEntries() == 0)
    createNvtx = false;
  if (scalLumi_v_lumi && scalLumi_v_lumi->getEntries() == 0)
    createInstLumi = false;

  double eff = 0.0;

  //Will pass if nvtx ME exists and is not empty
  if (createNvtx) {
    //Book new histos
    MonitorElement* eff_v_vtx_barrel =
        iBooker.book2D("hitefficiency_per_meanNvtx_per_PXLayer_PXBarrel",
                       "hitefficiency_per_meanNvtx_per_PXLayer_PXBarrel; meanNvtx; PXLayer",
                       500,
                       0,
                       100,
                       3,
                       .5,
                       3.5);

    MonitorElement* eff_v_vtx_forward =
        iBooker.book2D("hitefficiency_per_meanNvtx_per_PXDisk_PXForward",
                       "hitefficiency_per_meanNvtx_per_PXDisk_PXForward; meanNvtx; PXDisk",
                       500,
                       0,
                       100,
                       7,
                       -3.5,
                       3.5);

    //initialize variables
    int numLumiNvtx = int(vtx_v_lumi->getNbinsX());
    double nvtx = 0.0;
    int binNumVtx = 0;

    //For loop to loop through lumisections
    for (int iLumi = 1; iLumi < numLumiNvtx - 1; iLumi++) {
      //get the meanNvtx for each lumi
      nvtx = vtx_v_lumi->getBinContent(iLumi);

      //Filter out useless iterations
      if (nvtx != 0) {
        //Grab the bin number for the nvtx
        binNumVtx = eff_v_vtx_barrel->getTH2F()->FindBin(nvtx);

        //loop through the layers
        for (int iLayer = 1; iLayer < 8; iLayer++) {
          //get the eff at the lumisection and layer
          eff = eff_v_lumi_forward->getBinContent(iLumi - 1, iLayer);

          //set the efficiency in the new histo
          eff_v_vtx_forward->setBinContent(binNumVtx, iLayer, eff);
        }

        //loop through the layers
        for (int iLayer = 1; iLayer < 5; iLayer++) {
          //get the efficiency for each lumi at each layer
          eff = eff_v_lumi_barrel->getBinContent(iLumi - 1, iLayer);

          //set the efficiency
          eff_v_vtx_barrel->setBinContent(binNumVtx, iLayer, eff);
        }
      }
    }
  }
  // Will pass if InstLumi ME exists and is not empty
  if (createInstLumi) {
    //Get the max value of inst lumi for plot
    int yMax2 = scalLumi_v_lumi->getTProfile()->GetMaximum();
    yMax2 = yMax2 + yMax2 * .1;

    //Book new histos
    MonitorElement* eff_v_scalLumi_barrel =
        iBooker.book2D("hitefficiency_per_scalLumi_per_PXLayer_PXBarrel",
                       "hitefficiency_per_scalLumi_per_PXLayer_PXBarrel; scal inst lumi E30; PXLayer",
                       500,
                       0,
                       yMax2,
                       3,
                       .5,
                       3.5);

    MonitorElement* eff_v_scalLumi_forward =
        iBooker.book2D("hitefficiency_per_scalLumi_per_PXDisk_PXForward",
                       "hitefficiency_per_scalLumi_per_PXDisk_PXForward; scal inst lumi E30; PXDisk",
                       500,
                       0,
                       yMax2,
                       7,
                       -3.5,
                       3.5);

    //initialize variables
    int numLumiScal = int(scalLumi_v_lumi->getNbinsX());
    double scalLumi = 0.0;
    int binNumScal = 0;

    //For loop to loop through lumisections
    for (int iLumi = 1; iLumi < numLumiScal - 1; iLumi++) {
      //get the inst lumi for each lumi
      scalLumi = scalLumi_v_lumi->getBinContent(iLumi);

      //Filter out useless iterations
      if (scalLumi != 0) {
        //Grab the bin number for the inst lumi
        binNumScal = eff_v_scalLumi_barrel->getTH2F()->FindBin(scalLumi);

        //loop through the layers
        for (int iLayer = 1; iLayer < 8; iLayer++) {
          //get the eff at the lumisection and layer
          eff = eff_v_lumi_forward->getBinContent(iLumi - 1, iLayer);

          //set the efficiency in the new histo
          eff_v_scalLumi_forward->setBinContent(binNumScal, iLayer, eff);
        }

        //loop through the layers
        for (int iLayer = 1; iLayer < 5; iLayer++) {
          //get the eff at the lumisection and layer
          eff = eff_v_lumi_barrel->getBinContent(iLumi - 1, iLayer);

          //set the efficiency in the new histo
          eff_v_scalLumi_barrel->setBinContent(binNumScal, iLayer, eff);
        }
      }
    }
  } else
    return;
}

Member Data Documentation

effFolderName_

const std::string SiPixelPhase1EfficiencyExtras::effFolderName_

protected

Definition at line 45 of file SiPixelPhase1EfficiencyExtras.cc.

Referenced by dqmEndJob().

instLumiFolderName_

const std::string SiPixelPhase1EfficiencyExtras::instLumiFolderName_

protected

Definition at line 47 of file SiPixelPhase1EfficiencyExtras.cc.

Referenced by dqmEndJob().

vtxFolderName_

const std::string SiPixelPhase1EfficiencyExtras::vtxFolderName_

protected

Definition at line 46 of file SiPixelPhase1EfficiencyExtras.cc.

Referenced by dqmEndJob().