SiPixelHitEfficiencyModule.cc

Go to the documentation of this file.

// Package:    SiPixelMonitorTrack
// Class:      SiPixelHitEfficiencyModule
//
// class SiPixelHitEfficiencyModule SiPixelHitEfficiencyModule.cc
//       DQM/SiPixelMonitorTrack/src/SiPixelHitEfficiencyModule.cc
//
// Description: SiPixel hit efficiency data quality monitoring modules
// Implementation: prototype -> improved -> never final - end of the 1st step
//
// Original Authors: Romain Rougny & Luca Mucibello
//         Created: Mar Nov 10 13:29:00 CET nbinangle9
 
#include <iostream>
#include <string>
 
#include "DQM/SiPixelCommon/interface/SiPixelHistogramId.h"
#include "DQM/SiPixelMonitorTrack/interface/SiPixelHitEfficiencyModule.h"
#include "DQMServices/Core/interface/DQMStore.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "FWCore/Utilities/interface/InputTag.h"
 
// Data Formats
#include "DataFormats/DetId/interface/DetId.h"
#include "DataFormats/TrackerCommon/interface/PixelBarrelName.h"
#include "DataFormats/SiPixelDetId/interface/PixelBarrelNameUpgrade.h"
#include "DataFormats/TrackerCommon/interface/PixelEndcapName.h"
#include "DataFormats/SiPixelDetId/interface/PixelEndcapNameUpgrade.h"
#include "DataFormats/SiPixelDetId/interface/PixelSubdetector.h"
#include "DataFormats/TrackerCommon/interface/TrackerTopology.h"
#include "DataFormats/VertexReco/interface/Vertex.h"
#include "DataFormats/VertexReco/interface/VertexFwd.h"
#include "Geometry/Records/interface/TrackerTopologyRcd.h"
 
using namespace std;
 
SiPixelHitEfficiencyModule::SiPixelHitEfficiencyModule() : id_(0) { bBookTracks = true; }
 
SiPixelHitEfficiencyModule::SiPixelHitEfficiencyModule(uint32_t id) : id_(id) { bBookTracks = true; }
 
SiPixelHitEfficiencyModule::~SiPixelHitEfficiencyModule() {}
 
void SiPixelHitEfficiencyModule::book(const edm::ParameterSet &iConfig,
                                      edm::EventSetup const &iSetup,
                                      DQMStore::IBooker &iBooker,
                                      int type,
                                      bool isUpgrade) {
  edm::ESHandle<TrackerTopology> tTopoHandle;
  iSetup.get<TrackerTopologyRcd>().get(tTopoHandle);
  const TrackerTopology *pTT = tTopoHandle.product();
 
  bool barrel = DetId(id_).subdetId() == static_cast<int>(PixelSubdetector::PixelBarrel);
  bool endcap = DetId(id_).subdetId() == static_cast<int>(PixelSubdetector::PixelEndcap);
  bool isHalfModule = false;
  if (barrel) {
    isHalfModule = PixelBarrelName(DetId(id_), pTT, isUpgrade).isHalfModule();
  }
 
  edm::InputTag src = iConfig.getParameter<edm::InputTag>("src");
  debug_ = iConfig.getUntrackedParameter<bool>("debug", false);
  updateEfficiencies = iConfig.getUntrackedParameter<bool>("updateEfficiencies", false);
  std::string hisID;
 
  int nbinangle = 28;
  int nbinX = 45;
  int nbinY = 40;
 
  if (type == 0) {
    SiPixelHistogramId *theHistogramId = new SiPixelHistogramId(src.label());
 
    if (updateEfficiencies) {
      // EFFICIENCY
      hisID = theHistogramId->setHistoId("efficiency", id_);
      meEfficiency_ = iBooker.book1D(hisID, "Hit efficiency", 1, 0, 1.);
      meEfficiency_->setAxisTitle("Hit efficiency", 1);
 
      hisID = theHistogramId->setHistoId("efficiencyX", id_);
      meEfficiencyX_ = iBooker.book1D(hisID, "Hit efficiency in X", nbinX, -1.5, 1.5);
      meEfficiencyX_->setAxisTitle("Hit efficiency in X", 1);
 
      hisID = theHistogramId->setHistoId("efficiencyY", id_);
      meEfficiencyY_ = iBooker.book1D(hisID, "Hit efficiency in Y", nbinY, -4., 4.);
      meEfficiencyY_->setAxisTitle("Hit efficiency in Y", 1);
 
      hisID = theHistogramId->setHistoId("efficiencyAlpha", id_);
      meEfficiencyAlpha_ = iBooker.book1D(hisID, "Hit efficiency in Alpha", nbinangle, -3.5, 3.5);
      meEfficiencyAlpha_->setAxisTitle("Hit efficiency in Alpha", 1);
 
      hisID = theHistogramId->setHistoId("efficiencyBeta", id_);
      meEfficiencyBeta_ = iBooker.book1D(hisID, "Hit efficiency in Beta", nbinangle, -3.5, 3.5);
      meEfficiencyBeta_->setAxisTitle("Hit efficiency in Beta", 1);
    }
 
    // VALID
    hisID = theHistogramId->setHistoId("valid", id_);
    meValid_ = iBooker.book1D(hisID, "# Valid hits", 1, 0, 1.);
    meValid_->setAxisTitle("# Valid hits", 1);
 
    // MISSING
    hisID = theHistogramId->setHistoId("missing", id_);
    meMissing_ = iBooker.book1D(hisID, "# Missing hits", 1, 0, 1.);
    meMissing_->setAxisTitle("# Missing hits", 1);
 
    delete theHistogramId;
  }
 
  if (type == 1 && barrel) {
    uint32_t DBladder;
    DBladder = PixelBarrelName(DetId(id_), pTT, isUpgrade).ladderName();
    char sladder[80];
    sprintf(sladder, "Ladder_%02i", DBladder);
    hisID = src.label() + "_" + sladder;
    if (isHalfModule)
      hisID += "H";
    else
      hisID += "F";
 
    if (updateEfficiencies) {
      // EFFICIENCY
      meEfficiencyLad_ = iBooker.book1D("efficiency_" + hisID, "Hit efficiency", 1, 0, 1.);
      meEfficiencyLad_->setAxisTitle("Hit efficiency", 1);
 
      meEfficiencyXLad_ = iBooker.book1D("efficiencyX_" + hisID, "Hit efficiency in X", nbinX, -1.5, 1.5);
      meEfficiencyXLad_->setAxisTitle("Hit efficiency in X", 1);
 
      meEfficiencyYLad_ = iBooker.book1D("efficiencyY_" + hisID, "Hit efficiency in Y", nbinY, -4., 4.);
      meEfficiencyYLad_->setAxisTitle("Hit efficiency in Y", 1);
 
      meEfficiencyAlphaLad_ =
          iBooker.book1D("efficiencyAlpha_" + hisID, "Hit efficiency in Alpha", nbinangle, -3.5, 3.5);
      meEfficiencyAlphaLad_->setAxisTitle("Hit efficiency in Alpha", 1);
 
      meEfficiencyBetaLad_ = iBooker.book1D("efficiencyBeta_" + hisID, "Hit efficiency in Beta", nbinangle, -3.5, 3.5);
      meEfficiencyBetaLad_->setAxisTitle("Hit efficiency in Beta", 1);
    }
 
    // VALID
    meValidLad_ = iBooker.book1D("valid_" + hisID, "# Valid hits", 1, 0, 1.);
    meValidLad_->setAxisTitle("# Valid hits", 1);
 
    meValidXLad_ = iBooker.book1D("validX_" + hisID, "# Valid hits in X", nbinX, -1.5, 1.5);
    meValidXLad_->setAxisTitle("# Valid hits in X", 1);
 
    meValidYLad_ = iBooker.book1D("validY_" + hisID, "# Valid hits in Y", nbinY, -4., 4.);
    meValidYLad_->setAxisTitle("# Valid hits in Y", 1);
 
    meValidModLad_ = iBooker.book1D("validMod_" + hisID, "# Valid hits on Module", 4, 0.5, 4.5);
    meValidModLad_->setAxisTitle("# Valid hits on Module", 1);
 
    meValidAlphaLad_ = iBooker.book1D("validAlpha_" + hisID, "# Valid hits in Alpha", nbinangle, -3.5, 3.5);
    meValidAlphaLad_->setAxisTitle("# Valid hits in Alpha", 1);
 
    meValidBetaLad_ = iBooker.book1D("validBeta_" + hisID, "# Valid hits in Beta", nbinangle, -3.5, 3.5);
    meValidBetaLad_->setAxisTitle("# Valid hits in Beta", 1);
 
    // MISSING
    meMissingLad_ = iBooker.book1D("missing_" + hisID, "# Missing hits", 1, 0, 1.);
    meMissingLad_->setAxisTitle("# Missing hits", 1);
 
    meMissingXLad_ = iBooker.book1D("missingX_" + hisID, "# Missing hits in X", nbinX, -1.5, 1.5);
    meMissingXLad_->setAxisTitle("# Missing hits in X", 1);
 
    meMissingYLad_ = iBooker.book1D("missingY_" + hisID, "# Missing hits in Y", nbinY, -4., 4.);
    meMissingYLad_->setAxisTitle("# Missing hits in Y", 1);
 
    meMissingModLad_ = iBooker.book1D("missingMod_" + hisID, "# Missing hits on Module", 4, 0.5, 4.5);
    meMissingModLad_->setAxisTitle("# Missing hits on Module", 1);
 
    meMissingAlphaLad_ = iBooker.book1D("missingAlpha_" + hisID, "# Missing hits in Alpha", nbinangle, -3.5, 3.5);
    meMissingAlphaLad_->setAxisTitle("# Missing hits in Alpha", 1);
 
    meMissingBetaLad_ = iBooker.book1D("missingBeta_" + hisID, "# Missing hits in Beta", nbinangle, -3.5, 3.5);
    meMissingBetaLad_->setAxisTitle("# Missing hits in Beta", 1);
  }
 
  if (type == 2 && barrel) {
    uint32_t DBlayer;
    DBlayer = PixelBarrelName(DetId(id_), pTT, isUpgrade).layerName();
    char slayer[80];
    sprintf(slayer, "Layer_%i", DBlayer);
    hisID = src.label() + "_" + slayer;
 
    if (updateEfficiencies) {
      // EFFICIENCY
      meEfficiencyLay_ = iBooker.book1D("efficiency_" + hisID, "Hit efficiency", 1, 0, 1.);
      meEfficiencyLay_->setAxisTitle("Hit efficiency", 1);
 
      meEfficiencyXLay_ = iBooker.book1D("efficiencyX_" + hisID, "Hit efficiency in X", nbinX, -1.5, 1.5);
      meEfficiencyXLay_->setAxisTitle("Hit efficiency in X", 1);
 
      meEfficiencyYLay_ = iBooker.book1D("efficiencyY_" + hisID, "Hit efficiency in Y", nbinY, -4., 4.);
      meEfficiencyYLay_->setAxisTitle("Hit efficiency in Y", 1);
 
      meEfficiencyAlphaLay_ =
          iBooker.book1D("efficiencyAlpha_" + hisID, "Hit efficiency in Alpha", nbinangle, -3.5, 3.5);
      meEfficiencyAlphaLay_->setAxisTitle("Hit efficiency in Alpha", 1);
 
      meEfficiencyBetaLay_ = iBooker.book1D("efficiencyBeta_" + hisID, "Hit efficiency in Beta", nbinangle, -3.5, 3.5);
      meEfficiencyBetaLay_->setAxisTitle("Hit efficiency in Beta", 1);
    }
 
    // VALID
    meValidLay_ = iBooker.book1D("valid_" + hisID, "# Valid hits", 1, 0, 1.);
    meValidLay_->setAxisTitle("# Valid hits", 1);
 
    meValidXLay_ = iBooker.book1D("validX_" + hisID, "# Valid hits in X", nbinX, -1.5, 1.5);
    meValidXLay_->setAxisTitle("# Valid hits in X", 1);
 
    meValidYLay_ = iBooker.book1D("validY_" + hisID, "# Valid hits in Y", nbinY, -4., 4.);
    meValidYLay_->setAxisTitle("# Valid hits in Y", 1);
 
    meValidAlphaLay_ = iBooker.book1D("validAlpha_" + hisID, "# Valid hits in Alpha", nbinangle, -3.5, 3.5);
    meValidAlphaLay_->setAxisTitle("# Valid hits in Alpha", 1);
 
    meValidBetaLay_ = iBooker.book1D("validBeta_" + hisID, "# Valid hits in Beta", nbinangle, -3.5, 3.5);
    meValidBetaLay_->setAxisTitle("# Valid hits in Beta", 1);
 
    // MISSING
    meMissingLay_ = iBooker.book1D("missing_" + hisID, "# Missing hits", 1, 0, 1.);
    meMissingLay_->setAxisTitle("# Missing hits", 1);
 
    meMissingXLay_ = iBooker.book1D("missingX_" + hisID, "# Missing hits in X", nbinX, -1.5, 1.5);
    meMissingXLay_->setAxisTitle("# Missing hits in X", 1);
 
    meMissingYLay_ = iBooker.book1D("missingY_" + hisID, "# Missing hits in Y", nbinY, -4., 4.);
    meMissingYLay_->setAxisTitle("# Missing hits in Y", 1);
 
    meMissingAlphaLay_ = iBooker.book1D("missingAlpha_" + hisID, "# Missing hits in Alpha", nbinangle, -3.5, 3.5);
    meMissingAlphaLay_->setAxisTitle("# Missing hits in Alpha", 1);
 
    meMissingBetaLay_ = iBooker.book1D("missingBeta_" + hisID, "# Missing hits in Beta", nbinangle, -3.5, 3.5);
    meMissingBetaLay_->setAxisTitle("# Missing hits in Beta", 1);
  }
 
  if (type == 3 && barrel) {
    uint32_t DBmodule;
    DBmodule = PixelBarrelName(DetId(id_), pTT, isUpgrade).moduleName();
    char smodule[80];
    sprintf(smodule, "Ring_%i", DBmodule);
    hisID = src.label() + "_" + smodule;
 
    if (updateEfficiencies) {
      // EFFICIENCY
      meEfficiencyPhi_ = iBooker.book1D("efficiency_" + hisID, "Hit efficiency", 1, 0, 1.);
      meEfficiencyPhi_->setAxisTitle("Hit efficiency", 1);
 
      meEfficiencyXPhi_ = iBooker.book1D("efficiencyX_" + hisID, "Hit efficiency in X", nbinX, -1.5, 1.5);
      meEfficiencyXPhi_->setAxisTitle("Hit efficiency in X", 1);
 
      meEfficiencyYPhi_ = iBooker.book1D("efficiencyY_" + hisID, "Hit efficiency in Y", nbinY, -4., 4.);
      meEfficiencyYPhi_->setAxisTitle("Hit efficiency in Y", 1);
 
      meEfficiencyAlphaPhi_ =
          iBooker.book1D("efficiencyAlpha_" + hisID, "Hit efficiency in Alpha", nbinangle, -3.5, 3.5);
      meEfficiencyAlphaPhi_->setAxisTitle("Hit efficiency in Alpha", 1);
 
      meEfficiencyBetaPhi_ = iBooker.book1D("efficiencyBeta_" + hisID, "Hit efficiency in Beta", nbinangle, -3.5, 3.5);
      meEfficiencyBetaPhi_->setAxisTitle("Hit efficiency in Beta", 1);
    }
 
    // VALID
    meValidPhi_ = iBooker.book1D("valid_" + hisID, "# Valid hits", 1, 0, 1.);
    meValidPhi_->setAxisTitle("# Valid hits", 1);
 
    meValidXPhi_ = iBooker.book1D("validX_" + hisID, "# Valid hits in X", nbinX, -1.5, 1.5);
    meValidXPhi_->setAxisTitle("# Valid hits in X", 1);
 
    meValidYPhi_ = iBooker.book1D("validY_" + hisID, "# Valid hits in Y", nbinY, -4., 4.);
    meValidYPhi_->setAxisTitle("# Valid hits in Y", 1);
 
    meValidAlphaPhi_ = iBooker.book1D("validAlpha_" + hisID, "# Valid hits in Alpha", nbinangle, -3.5, 3.5);
    meValidAlphaPhi_->setAxisTitle("# Valid hits in Alpha", 1);
 
    meValidBetaPhi_ = iBooker.book1D("validBeta_" + hisID, "# Valid hits in Beta", nbinangle, -3.5, 3.5);
    meValidBetaPhi_->setAxisTitle("# Valid hits in Beta", 1);
 
    // MISSING
    meMissingPhi_ = iBooker.book1D("missing_" + hisID, "# Missing hits", 1, 0, 1.);
    meMissingPhi_->setAxisTitle("# Missing hits", 1);
 
    meMissingXPhi_ = iBooker.book1D("missingX_" + hisID, "# Missing hits in X", nbinX, -1.5, 1.5);
    meMissingXPhi_->setAxisTitle("# Missing hits in X", 1);
 
    meMissingYPhi_ = iBooker.book1D("missingY_" + hisID, "# Missing hits in Y", nbinY, -4., 4.);
    meMissingYPhi_->setAxisTitle("# Missing hits in Y", 1);
 
    meMissingAlphaPhi_ = iBooker.book1D("missingAlpha_" + hisID, "# Missing hits in Alpha", nbinangle, -3.5, 3.5);
    meMissingAlphaPhi_->setAxisTitle("# Missing hits in Alpha", 1);
 
    meMissingBetaPhi_ = iBooker.book1D("missingBeta_" + hisID, "# Missing hits in Beta", nbinangle, -3.5, 3.5);
    meMissingBetaPhi_->setAxisTitle("# Missing hits in Beta", 1);
  }
 
  if (type == 4 && endcap) {
    uint32_t blade;
    blade = PixelEndcapName(DetId(id_), pTT, isUpgrade).bladeName();
 
    char sblade[80];
    sprintf(sblade, "Blade_%02i", blade);
    hisID = src.label() + "_" + sblade;
 
    if (updateEfficiencies) {
      // EFFICIENCY
      meEfficiencyBlade_ = iBooker.book1D("efficiency_" + hisID, "Hit efficiency", 1, 0, 1.);
      meEfficiencyBlade_->setAxisTitle("Hit efficiency", 1);
 
      meEfficiencyXBlade_ = iBooker.book1D("efficiencyX_" + hisID, "Hit efficiency in X", nbinX, -1.5, 1.5);
      meEfficiencyXBlade_->setAxisTitle("Hit efficiency in X", 1);
 
      meEfficiencyYBlade_ = iBooker.book1D("efficiencyY_" + hisID, "Hit efficiency in Y", nbinY, -4., 4.);
      meEfficiencyYBlade_->setAxisTitle("Hit efficiency in Y", 1);
 
      meEfficiencyAlphaBlade_ =
          iBooker.book1D("efficiencyAlpha_" + hisID, "Hit efficiency in Alpha", nbinangle, -3.5, 3.5);
      meEfficiencyAlphaBlade_->setAxisTitle("Hit efficiency in Alpha", 1);
 
      meEfficiencyBetaBlade_ =
          iBooker.book1D("efficiencyBeta_" + hisID, "Hit efficiency in Beta", nbinangle, -3.5, 3.5);
      meEfficiencyBetaBlade_->setAxisTitle("Hit efficiency in Beta", 1);
    }
 
    // VALID
    meValidBlade_ = iBooker.book1D("valid_" + hisID, "# Valid hits", 2, 0.5, 2.5);
    meValidBlade_->setAxisTitle("# Valid hits", 1);
 
    meValidXBlade_ = iBooker.book1D("validX_" + hisID, "# Valid hits in X", nbinX, -1.5, 1.5);
    meValidXBlade_->setAxisTitle("# Valid hits in X", 1);
 
    meValidYBlade_ = iBooker.book1D("validY_" + hisID, "# Valid hits in Y", nbinY, -4., 4.);
    meValidYBlade_->setAxisTitle("# Valid hits in Y", 1);
 
    meValidAlphaBlade_ = iBooker.book1D("validAlpha_" + hisID, "# Valid hits in Alpha", nbinangle, -3.5, 3.5);
    meValidAlphaBlade_->setAxisTitle("# Valid hits in Alpha", 1);
 
    meValidBetaBlade_ = iBooker.book1D("validBeta_" + hisID, "# Valid hits in Beta", nbinangle, -3.5, 3.5);
    meValidBetaBlade_->setAxisTitle("# Valid hits in Beta", 1);
 
    // MISSING
    meMissingBlade_ = iBooker.book1D("missing_" + hisID, "# Missing hits", 2, 0.5, 2.5);
    meMissingBlade_->setAxisTitle("# Missing hits", 1);
 
    meMissingXBlade_ = iBooker.book1D("missingX_" + hisID, "# Missing hits in X", nbinX, -1.5, 1.5);
    meMissingXBlade_->setAxisTitle("# Missing hits in X", 1);
 
    meMissingYBlade_ = iBooker.book1D("missingY_" + hisID, "# Missing hits in Y", nbinY, -4., 4.);
    meMissingYBlade_->setAxisTitle("# Missing hits in Y", 1);
 
    meMissingAlphaBlade_ = iBooker.book1D("missingAlpha_" + hisID, "# Missing hits in Alpha", nbinangle, -3.5, 3.5);
    meMissingAlphaBlade_->setAxisTitle("# Missing hits in Alpha", 1);
 
    meMissingBetaBlade_ = iBooker.book1D("missingBeta_" + hisID, "# Missing hits in Beta", nbinangle, -3.5, 3.5);
    meMissingBetaBlade_->setAxisTitle("# Missing hits in Beta", 1);
  }
 
  if (type == 5 && endcap) {
    uint32_t disk;
    disk = PixelEndcapName(DetId(id_), pTT, isUpgrade).diskName();
 
    char sdisk[80];
    sprintf(sdisk, "Disk_%i", disk);
    hisID = src.label() + "_" + sdisk;
 
    if (updateEfficiencies) {
      // EFFICIENCY
      meEfficiencyDisk_ = iBooker.book1D("efficiency_" + hisID, "Hit efficiency", 1, 0, 1.);
      meEfficiencyDisk_->setAxisTitle("Hit efficiency", 1);
 
      meEfficiencyXDisk_ = iBooker.book1D("efficiencyX_" + hisID, "Hit efficiency in X", nbinX, -1.5, 1.5);
      meEfficiencyXDisk_->setAxisTitle("Hit efficiency in X", 1);
 
      meEfficiencyYDisk_ = iBooker.book1D("efficiencyY_" + hisID, "Hit efficiency in Y", nbinY, -4., 4.);
      meEfficiencyYDisk_->setAxisTitle("Hit efficiency in Y", 1);
 
      meEfficiencyAlphaDisk_ =
          iBooker.book1D("efficiencyAlpha_" + hisID, "Hit efficiency in Alpha", nbinangle, -3.5, 3.5);
      meEfficiencyAlphaDisk_->setAxisTitle("Hit efficiency in Alpha", 1);
 
      meEfficiencyBetaDisk_ = iBooker.book1D("efficiencyBeta_" + hisID, "Hit efficiency in Beta", nbinangle, -3.5, 3.5);
      meEfficiencyBetaDisk_->setAxisTitle("Hit efficiency in Beta", 1);
    }
 
    // VALID
    meValidDisk_ = iBooker.book1D("valid_" + hisID, "# Valid hits", 1, 0, 1.);
    meValidDisk_->setAxisTitle("# Valid hits", 1);
 
    meValidXDisk_ = iBooker.book1D("validX_" + hisID, "# Valid hits in X", nbinX, -1.5, 1.5);
    meValidXDisk_->setAxisTitle("# Valid hits in X", 1);
 
    meValidYDisk_ = iBooker.book1D("validY_" + hisID, "# Valid hits in Y", nbinY, -4., 4.);
    meValidYDisk_->setAxisTitle("# Valid hits in Y", 1);
 
    meValidAlphaDisk_ = iBooker.book1D("validAlpha_" + hisID, "# Valid hits in Alpha", nbinangle, -3.5, 3.5);
    meValidAlphaDisk_->setAxisTitle("# Valid hits in Alpha", 1);
 
    meValidBetaDisk_ = iBooker.book1D("validBeta_" + hisID, "# Valid hits in Beta", nbinangle, -3.5, 3.5);
    meValidBetaDisk_->setAxisTitle("# Valid hits in Beta", 1);
 
    // MISSING
    meMissingDisk_ = iBooker.book1D("missing_" + hisID, "# Missing hits", 1, 0, 1.);
    meMissingDisk_->setAxisTitle("# Missing hits", 1);
 
    meMissingXDisk_ = iBooker.book1D("missingX_" + hisID, "# Missing hits in X", nbinX, -1.5, 1.5);
    meMissingXDisk_->setAxisTitle("# Missing hits in X", 1);
 
    meMissingYDisk_ = iBooker.book1D("missingY_" + hisID, "# Missing hits in Y", nbinY, -4., 4.);
    meMissingYDisk_->setAxisTitle("# Missing hits in Y", 1);
 
    meMissingAlphaDisk_ = iBooker.book1D("missingAlpha_" + hisID, "# Missing hits in Alpha", nbinangle, -3.5, 3.5);
    meMissingAlphaDisk_->setAxisTitle("# Missing hits in Alpha", 1);
 
    meMissingBetaDisk_ = iBooker.book1D("missingBeta_" + hisID, "# Missing hits in Beta", nbinangle, -3.5, 3.5);
    meMissingBetaDisk_->setAxisTitle("# Missing hits in Beta", 1);
  }
 
  if (type == 6 && endcap) {
    uint32_t panel;
    uint32_t module;
    panel = PixelEndcapName(DetId(id_), pTT, isUpgrade).pannelName();
    module = PixelEndcapName(DetId(id_), pTT, isUpgrade).plaquetteName();
 
    char slab[80];
    sprintf(slab, "Panel_%i_Ring_%i", panel, module);
    hisID = src.label() + "_" + slab;
 
    if (updateEfficiencies) {
      // EFFICIENCY
      meEfficiencyRing_ = iBooker.book1D("efficiency_" + hisID, "Hit efficiency", 1, 0, 1.);
      meEfficiencyRing_->setAxisTitle("Hit efficiency", 1);
 
      meEfficiencyXRing_ = iBooker.book1D("efficiencyX_" + hisID, "Hit efficiency in X", nbinX, -1.5, 1.5);
      meEfficiencyXRing_->setAxisTitle("Hit efficiency in X", 1);
 
      meEfficiencyYRing_ = iBooker.book1D("efficiencyY_" + hisID, "Hit efficiency in Y", nbinY, -4., 4.);
      meEfficiencyYRing_->setAxisTitle("Hit efficiency in Y", 1);
 
      meEfficiencyAlphaRing_ =
          iBooker.book1D("efficiencyAlpha_" + hisID, "Hit efficiency in Alpha", nbinangle, -3.5, 3.5);
      meEfficiencyAlphaRing_->setAxisTitle("Hit efficiency in Alpha", 1);
 
      meEfficiencyBetaRing_ = iBooker.book1D("efficiencyBeta_" + hisID, "Hit efficiency in Beta", nbinangle, -3.5, 3.5);
      meEfficiencyBetaRing_->setAxisTitle("Hit efficiency in Beta", 1);
    }
 
    // VALID
    meValidRing_ = iBooker.book1D("valid_" + hisID, "# Valid hits", 1, 0, 1.);
    meValidRing_->setAxisTitle("# Valid hits", 1);
 
    meValidXRing_ = iBooker.book1D("validX_" + hisID, "# Valid hits in X", nbinX, -1.5, 1.5);
    meValidXRing_->setAxisTitle("# Valid hits in X", 1);
 
    meValidYRing_ = iBooker.book1D("validY_" + hisID, "# Valid hits in Y", nbinY, -4., 4.);
    meValidYRing_->setAxisTitle("# Valid hits in Y", 1);
 
    meValidAlphaRing_ = iBooker.book1D("validAlpha_" + hisID, "# Valid hits in Alpha", nbinangle, -3.5, 3.5);
    meValidAlphaRing_->setAxisTitle("# Valid hits in Alpha", 1);
 
    meValidBetaRing_ = iBooker.book1D("validBeta_" + hisID, "# Valid hits in Beta", nbinangle, -3.5, 3.5);
    meValidBetaRing_->setAxisTitle("# Valid hits in Beta", 1);
 
    // MISSING
    meMissingRing_ = iBooker.book1D("missing_" + hisID, "# Missing hits", 1, 0, 1.);
    meMissingRing_->setAxisTitle("# Missing hits", 1);
 
    meMissingXRing_ = iBooker.book1D("missingX_" + hisID, "# Missing hits in X", nbinX, -1.5, 1.5);
    meMissingXRing_->setAxisTitle("# Missing hits in X", 1);
 
    meMissingYRing_ = iBooker.book1D("missingY_" + hisID, "# Missing hits in Y", nbinY, -4., 4.);
    meMissingYRing_->setAxisTitle("# Missing hits in Y", 1);
 
    meMissingAlphaRing_ = iBooker.book1D("missingAlpha_" + hisID, "# Missing hits in Alpha", nbinangle, -3.5, 3.5);
    meMissingAlphaRing_->setAxisTitle("# Missing hits in Alpha", 1);
 
    meMissingBetaRing_ = iBooker.book1D("missingBeta_" + hisID, "# Missing hits in Beta", nbinangle, -3.5, 3.5);
    meMissingBetaRing_->setAxisTitle("# Missing hits in Beta", 1);
  }
}
 
void SiPixelHitEfficiencyModule::fill(const TrackerTopology *pTT,
                                      const LocalTrajectoryParameters &ltp,
                                      bool isHitValid,
                                      bool modon,
                                      bool ladon,
                                      bool layon,
                                      bool phion,
                                      bool bladeon,
                                      bool diskon,
                                      bool ringon) {
  bool barrel = DetId(id_).subdetId() == static_cast<int>(PixelSubdetector::PixelBarrel);
  bool endcap = DetId(id_).subdetId() == static_cast<int>(PixelSubdetector::PixelEndcap);
 
  LocalVector localDir = ltp.momentum() / ltp.momentum().mag();
  float prediction_alpha = atan2(localDir.z(), localDir.x());
  float prediction_beta = atan2(localDir.z(), localDir.y());
  float prediction_x = ltp.position().x();
  float prediction_y = ltp.position().y();
 
  PixelBarrelName PBN = PixelBarrelName(DetId(id_), pTT);
  int imod = PBN.moduleName();
  PixelEndcapName PEN = PixelEndcapName(DetId(id_), pTT);
  int ipan = PEN.pannelName();
 
  if (isHitValid) {
    if (modon) {
      meValid_->Fill(0.5);
    }
    if (barrel && ladon) {
      meValidLad_->Fill(0.5);
      meValidXLad_->Fill(prediction_x);
      meValidYLad_->Fill(prediction_y);
      meValidModLad_->Fill(imod);
      meValidAlphaLad_->Fill(prediction_alpha);
      meValidBetaLad_->Fill(prediction_beta);
    }
    if (barrel && layon) {
      meValidLay_->Fill(0.5);
      meValidXLay_->Fill(prediction_x);
      meValidYLay_->Fill(prediction_y);
      meValidAlphaLay_->Fill(prediction_alpha);
      meValidBetaLay_->Fill(prediction_beta);
    }
    if (barrel && phion) {
      meValidPhi_->Fill(0.5);
      meValidXPhi_->Fill(prediction_x);
      meValidYPhi_->Fill(prediction_y);
      meValidAlphaPhi_->Fill(prediction_alpha);
      meValidBetaPhi_->Fill(prediction_beta);
    }
    if (endcap && bladeon) {
      meValidBlade_->Fill(ipan);
      meValidXBlade_->Fill(prediction_x);
      meValidYBlade_->Fill(prediction_y);
      meValidAlphaBlade_->Fill(prediction_alpha);
      meValidBetaBlade_->Fill(prediction_beta);
    }
    if (endcap && diskon) {
      meValidDisk_->Fill(0.5);
      meValidXDisk_->Fill(prediction_x);
      meValidYDisk_->Fill(prediction_y);
      meValidAlphaDisk_->Fill(prediction_alpha);
      meValidBetaDisk_->Fill(prediction_beta);
    }
    if (endcap && ringon) {
      meValidRing_->Fill(0.5);
      meValidXRing_->Fill(prediction_x);
      meValidYRing_->Fill(prediction_y);
      meValidAlphaRing_->Fill(prediction_alpha);
      meValidBetaRing_->Fill(prediction_beta);
    }
  } else {
    if (modon) {
      meMissing_->Fill(0.5);
    }
    if (barrel && ladon) {
      meMissingLad_->Fill(0.5);
      meMissingXLad_->Fill(prediction_x);
      meMissingYLad_->Fill(prediction_y);
      meMissingModLad_->Fill(imod);
      meMissingAlphaLad_->Fill(prediction_alpha);
      meMissingBetaLad_->Fill(prediction_beta);
    }
    if (barrel && layon) {
      meMissingLay_->Fill(0.5);
      meMissingXLay_->Fill(prediction_x);
      meMissingYLay_->Fill(prediction_y);
      meMissingAlphaLay_->Fill(prediction_alpha);
      meMissingBetaLay_->Fill(prediction_beta);
    }
    if (barrel && phion) {
      meMissingPhi_->Fill(0.5);
      meMissingXPhi_->Fill(prediction_x);
      meMissingYPhi_->Fill(prediction_y);
      meMissingAlphaPhi_->Fill(prediction_alpha);
      meMissingBetaPhi_->Fill(prediction_beta);
    }
    if (endcap && bladeon) {
      meMissingBlade_->Fill(ipan);
      meMissingXBlade_->Fill(prediction_x);
      meMissingYBlade_->Fill(prediction_y);
      meMissingAlphaBlade_->Fill(prediction_alpha);
      meMissingBetaBlade_->Fill(prediction_beta);
    }
    if (endcap && diskon) {
      meMissingDisk_->Fill(0.5);
      meMissingXDisk_->Fill(prediction_x);
      meMissingYDisk_->Fill(prediction_y);
      meMissingAlphaDisk_->Fill(prediction_alpha);
      meMissingBetaDisk_->Fill(prediction_beta);
    }
    if (endcap && ringon) {
      meMissingRing_->Fill(0.5);
      meMissingXRing_->Fill(prediction_x);
      meMissingYRing_->Fill(prediction_y);
      meMissingAlphaRing_->Fill(prediction_alpha);
      meMissingBetaRing_->Fill(prediction_beta);
    }
  }
 
  if (updateEfficiencies)
    computeEfficiencies(modon, ladon, layon, phion, bladeon, diskon, ringon);
}
 
void SiPixelHitEfficiencyModule::computeEfficiencies(
    bool modon, bool ladon, bool layon, bool phion, bool bladeon, bool diskon, bool ringon) {
  if (debug_)
    std::cout << "Now Filling histos for detid " << id_ << std::endl;
 
  bool barrel = DetId(id_).subdetId() == static_cast<int>(PixelSubdetector::PixelBarrel);
  bool endcap = DetId(id_).subdetId() == static_cast<int>(PixelSubdetector::PixelEndcap);
 
  if (modon) {
    meEfficiency_->setBinContent(1, (eff(meValid_->getBinContent(1), meMissing_->getBinContent(1))).first);
    meEfficiency_->setBinError(1, (eff(meValid_->getBinContent(1), meMissing_->getBinContent(1))).second);
  }
  if (ladon && barrel) {
    meEfficiencyLad_->setBinContent(1, (eff(meValidLad_->getBinContent(1), meMissingLad_->getBinContent(1))).first);
    meEfficiencyLad_->setBinError(1, (eff(meValidLad_->getBinContent(1), meMissingLad_->getBinContent(1))).second);
    for (int i = 1; i <= meValidXLad_->getNbinsX(); ++i) {
      meEfficiencyXLad_->setBinContent(i,
                                       (eff(meValidXLad_->getBinContent(i), meMissingXLad_->getBinContent(i))).first);
      meEfficiencyXLad_->setBinError(i, (eff(meValidXLad_->getBinContent(i), meMissingXLad_->getBinContent(i))).second);
    }
    for (int i = 1; i <= meValidYLad_->getNbinsX(); ++i) {
      meEfficiencyYLad_->setBinContent(i,
                                       (eff(meValidYLad_->getBinContent(i), meMissingYLad_->getBinContent(i))).first);
      meEfficiencyYLad_->setBinError(i, (eff(meValidYLad_->getBinContent(i), meMissingYLad_->getBinContent(i))).second);
    }
    for (int i = 1; i <= meValidAlphaLad_->getNbinsX(); ++i) {
      meEfficiencyAlphaLad_->setBinContent(
          i, (eff(meValidAlphaLad_->getBinContent(i), meMissingAlphaLad_->getBinContent(i))).first);
      meEfficiencyAlphaLad_->setBinError(
          i, (eff(meValidAlphaLad_->getBinContent(i), meMissingAlphaLad_->getBinContent(i))).second);
    }
    for (int i = 1; i <= meValidBetaLad_->getNbinsX(); ++i) {
      meEfficiencyBetaLad_->setBinContent(
          i, (eff(meValidBetaLad_->getBinContent(i), meMissingBetaLad_->getBinContent(i))).first);
      meEfficiencyBetaLad_->setBinError(
          i, (eff(meValidBetaLad_->getBinContent(i), meMissingBetaLad_->getBinContent(i))).second);
    }
  }
 
  if (layon && barrel) {
    meEfficiencyLay_->setBinContent(1, (eff(meValidLay_->getBinContent(1), meMissingLay_->getBinContent(1))).first);
    meEfficiencyLay_->setBinError(1, (eff(meValidLay_->getBinContent(1), meMissingLay_->getBinContent(1))).second);
    for (int i = 1; i <= meValidXLay_->getNbinsX(); ++i) {
      meEfficiencyXLay_->setBinContent(i,
                                       (eff(meValidXLay_->getBinContent(i), meMissingXLay_->getBinContent(i))).first);
      meEfficiencyXLay_->setBinError(i, (eff(meValidXLay_->getBinContent(i), meMissingXLay_->getBinContent(i))).second);
    }
    for (int i = 1; i <= meValidYLay_->getNbinsX(); ++i) {
      meEfficiencyYLay_->setBinContent(i,
                                       (eff(meValidYLay_->getBinContent(i), meMissingYLay_->getBinContent(i))).first);
      meEfficiencyYLay_->setBinError(i, (eff(meValidYLay_->getBinContent(i), meMissingYLay_->getBinContent(i))).second);
    }
    for (int i = 1; i <= meValidAlphaLay_->getNbinsX(); ++i) {
      meEfficiencyAlphaLay_->setBinContent(
          i, (eff(meValidAlphaLay_->getBinContent(i), meMissingAlphaLay_->getBinContent(i))).first);
      meEfficiencyAlphaLay_->setBinError(
          i, (eff(meValidAlphaLay_->getBinContent(i), meMissingAlphaLay_->getBinContent(i))).second);
    }
    for (int i = 1; i <= meValidBetaLay_->getNbinsX(); ++i) {
      meEfficiencyBetaLay_->setBinContent(
          i, (eff(meValidBetaLay_->getBinContent(i), meMissingBetaLay_->getBinContent(i))).first);
      meEfficiencyBetaLay_->setBinError(
          i, (eff(meValidBetaLay_->getBinContent(i), meMissingBetaLay_->getBinContent(i))).second);
    }
  }
 
  if (phion && barrel) {
    meEfficiencyPhi_->setBinContent(1, (eff(meValidPhi_->getBinContent(1), meMissingPhi_->getBinContent(1))).first);
    meEfficiencyPhi_->setBinError(1, (eff(meValidPhi_->getBinContent(1), meMissingPhi_->getBinContent(1))).second);
    for (int i = 1; i <= meValidXPhi_->getNbinsX(); ++i) {
      meEfficiencyXPhi_->setBinContent(i,
                                       (eff(meValidXPhi_->getBinContent(i), meMissingXPhi_->getBinContent(i))).first);
      meEfficiencyXPhi_->setBinError(i, (eff(meValidXPhi_->getBinContent(i), meMissingXPhi_->getBinContent(i))).second);
    }
    for (int i = 1; i <= meValidYPhi_->getNbinsX(); ++i) {
      meEfficiencyYPhi_->setBinContent(i,
                                       (eff(meValidYPhi_->getBinContent(i), meMissingYPhi_->getBinContent(i))).first);
      meEfficiencyYPhi_->setBinError(i, (eff(meValidYPhi_->getBinContent(i), meMissingYPhi_->getBinContent(i))).second);
    }
    for (int i = 1; i <= meValidAlphaPhi_->getNbinsX(); ++i) {
      meEfficiencyAlphaPhi_->setBinContent(
          i, (eff(meValidAlphaPhi_->getBinContent(i), meMissingAlphaPhi_->getBinContent(i))).first);
      meEfficiencyAlphaPhi_->setBinError(
          i, (eff(meValidAlphaPhi_->getBinContent(i), meMissingAlphaPhi_->getBinContent(i))).second);
    }
    for (int i = 1; i <= meValidBetaPhi_->getNbinsX(); ++i) {
      meEfficiencyBetaPhi_->setBinContent(
          i, (eff(meValidBetaPhi_->getBinContent(i), meMissingBetaPhi_->getBinContent(i))).first);
      meEfficiencyBetaPhi_->setBinError(
          i, (eff(meValidBetaPhi_->getBinContent(i), meMissingBetaPhi_->getBinContent(i))).second);
    }
  }
  if (bladeon && endcap) {
    meEfficiencyBlade_->setBinContent(1,
                                      (eff(meValidBlade_->getBinContent(1), meMissingBlade_->getBinContent(1))).first);
    meEfficiencyBlade_->setBinError(1,
                                    (eff(meValidBlade_->getBinContent(1), meMissingBlade_->getBinContent(1))).second);
    for (int i = 1; i <= meValidXBlade_->getNbinsX(); ++i) {
      meEfficiencyXBlade_->setBinContent(
          i, (eff(meValidXBlade_->getBinContent(i), meMissingXBlade_->getBinContent(i))).first);
      meEfficiencyXBlade_->setBinError(
          i, (eff(meValidXBlade_->getBinContent(i), meMissingXBlade_->getBinContent(i))).second);
    }
    for (int i = 1; i <= meValidYBlade_->getNbinsX(); ++i) {
      meEfficiencyYBlade_->setBinContent(
          i, (eff(meValidYBlade_->getBinContent(i), meMissingYBlade_->getBinContent(i))).first);
      meEfficiencyYBlade_->setBinError(
          i, (eff(meValidYBlade_->getBinContent(i), meMissingYBlade_->getBinContent(i))).second);
    }
    for (int i = 1; i <= meValidAlphaBlade_->getNbinsX(); ++i) {
      meEfficiencyAlphaBlade_->setBinContent(
          i, (eff(meValidAlphaBlade_->getBinContent(i), meMissingAlphaBlade_->getBinContent(i))).first);
      meEfficiencyAlphaBlade_->setBinError(
          i, (eff(meValidAlphaBlade_->getBinContent(i), meMissingAlphaBlade_->getBinContent(i))).second);
    }
    for (int i = 1; i <= meValidBetaBlade_->getNbinsX(); ++i) {
      meEfficiencyBetaBlade_->setBinContent(
          i, (eff(meValidBetaBlade_->getBinContent(i), meMissingBetaBlade_->getBinContent(i))).first);
      meEfficiencyBetaBlade_->setBinError(
          i, (eff(meValidBetaBlade_->getBinContent(i), meMissingBetaBlade_->getBinContent(i))).second);
    }
  }
  if (diskon && endcap) {
    meEfficiencyDisk_->setBinContent(1, (eff(meValidDisk_->getBinContent(1), meMissingDisk_->getBinContent(1))).first);
    meEfficiencyDisk_->setBinError(1, (eff(meValidDisk_->getBinContent(1), meMissingDisk_->getBinContent(1))).second);
    for (int i = 1; i <= meValidXDisk_->getNbinsX(); ++i) {
      meEfficiencyXDisk_->setBinContent(
          i, (eff(meValidXDisk_->getBinContent(i), meMissingXDisk_->getBinContent(i))).first);
      meEfficiencyXDisk_->setBinError(i,
                                      (eff(meValidXDisk_->getBinContent(i), meMissingXDisk_->getBinContent(i))).second);
    }
    for (int i = 1; i <= meValidYDisk_->getNbinsX(); ++i) {
      meEfficiencyYDisk_->setBinContent(
          i, (eff(meValidYDisk_->getBinContent(i), meMissingYDisk_->getBinContent(i))).first);
      meEfficiencyYDisk_->setBinError(i,
                                      (eff(meValidYDisk_->getBinContent(i), meMissingYDisk_->getBinContent(i))).second);
    }
    for (int i = 1; i <= meValidAlphaDisk_->getNbinsX(); ++i) {
      meEfficiencyAlphaDisk_->setBinContent(
          i, (eff(meValidAlphaDisk_->getBinContent(i), meMissingAlphaDisk_->getBinContent(i))).first);
      meEfficiencyAlphaDisk_->setBinError(
          i, (eff(meValidAlphaDisk_->getBinContent(i), meMissingAlphaDisk_->getBinContent(i))).second);
    }
    for (int i = 1; i <= meValidBetaDisk_->getNbinsX(); ++i) {
      meEfficiencyBetaDisk_->setBinContent(
          i, (eff(meValidBetaDisk_->getBinContent(i), meMissingBetaDisk_->getBinContent(i))).first);
      meEfficiencyBetaDisk_->setBinError(
          i, (eff(meValidBetaDisk_->getBinContent(i), meMissingBetaDisk_->getBinContent(i))).second);
    }
  }
  if (ringon && endcap) {
    meEfficiencyRing_->setBinContent(1, (eff(meValidRing_->getBinContent(1), meMissingRing_->getBinContent(1))).first);
    meEfficiencyRing_->setBinError(1, (eff(meValidRing_->getBinContent(1), meMissingRing_->getBinContent(1))).second);
    for (int i = 1; i <= meValidXRing_->getNbinsX(); ++i) {
      meEfficiencyXRing_->setBinContent(
          i, (eff(meValidXRing_->getBinContent(i), meMissingXRing_->getBinContent(i))).first);
      meEfficiencyXRing_->setBinError(i,
                                      (eff(meValidXRing_->getBinContent(i), meMissingXRing_->getBinContent(i))).second);
    }
    for (int i = 1; i <= meValidYRing_->getNbinsX(); ++i) {
      meEfficiencyYRing_->setBinContent(
          i, (eff(meValidYRing_->getBinContent(i), meMissingYRing_->getBinContent(i))).first);
      meEfficiencyYRing_->setBinError(i,
                                      (eff(meValidYRing_->getBinContent(i), meMissingYRing_->getBinContent(i))).second);
    }
    for (int i = 1; i <= meValidAlphaRing_->getNbinsX(); ++i) {
      meEfficiencyAlphaRing_->setBinContent(
          i, (eff(meValidAlphaRing_->getBinContent(i), meMissingAlphaRing_->getBinContent(i))).first);
      meEfficiencyAlphaRing_->setBinError(
          i, (eff(meValidAlphaRing_->getBinContent(i), meMissingAlphaRing_->getBinContent(i))).second);
    }
    for (int i = 1; i <= meValidBetaRing_->getNbinsX(); ++i) {
      meEfficiencyBetaRing_->setBinContent(
          i, (eff(meValidBetaRing_->getBinContent(i), meMissingBetaRing_->getBinContent(i))).first);
      meEfficiencyBetaRing_->setBinError(
          i, (eff(meValidBetaRing_->getBinContent(i), meMissingBetaRing_->getBinContent(i))).second);
    }
  }
}
 
std::pair<double, double> SiPixelHitEfficiencyModule::eff(double nValid, double nMissing) {
  double efficiency = 0, error = 0;
  if (nValid + nMissing != 0) {
    efficiency = nValid / (nValid + nMissing);
    error = sqrt(efficiency * (1. - efficiency) / (nValid + nMissing));
  }
  return make_pair(efficiency, error);
}