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 ( 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 ( )

override

Definition at line 57 of file SiPixelPhase1EfficiencyExtras.cc.

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

Member Function Documentation

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

overrideprotected

Definition at line 61 of file SiPixelPhase1EfficiencyExtras.cc.

61 {}

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_, 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;
 }