SiPixelDigiValid.cc

Go to the documentation of this file.

#include "DQMServices/Core/interface/DQMStore.h"
#include "DQMServices/Core/interface/MonitorElement.h"
#include "DataFormats/Common/interface/DetSetVector.h"
#include "DataFormats/SiPixelDetId/interface/PixelSubdetector.h"
#include "DataFormats/SiPixelDigi/interface/PixelDigi.h"
#include "DataFormats/TrackerCommon/interface/TrackerTopology.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "Geometry/Records/interface/TrackerDigiGeometryRecord.h"
#include "Geometry/Records/interface/TrackerTopologyRcd.h"
#include "Geometry/TrackerGeometryBuilder/interface/PixelGeomDetType.h"
#include "Geometry/TrackerGeometryBuilder/interface/PixelGeomDetUnit.h"
#include "Geometry/TrackerGeometryBuilder/interface/RectangularPixelTopology.h"
#include "Geometry/TrackerGeometryBuilder/interface/TrackerGeometry.h"
#include "Geometry/TrackerNumberingBuilder/interface/GeometricDet.h"
#include "Validation/TrackerDigis/interface/SiPixelDigiValid.h"

// using namespace std;
// using namespace edm;

SiPixelDigiValid::SiPixelDigiValid(const edm::ParameterSet &ps)
    : outputFile_(ps.getUntrackedParameter<std::string>("outputFile", "pixeldigihisto.root")),
      runStandalone(ps.getParameter<bool>("runStandalone")),
      dbe_(nullptr),
      edmDetSetVector_PixelDigi_Token_(consumes<edm::DetSetVector<PixelDigi>>(ps.getParameter<edm::InputTag>("src"))) {}

SiPixelDigiValid::~SiPixelDigiValid() {}

void SiPixelDigiValid::bookHistograms(DQMStore::IBooker &ibooker, const edm::Run &run, const edm::EventSetup &es) {
  dbe_ = edm::Service<DQMStore>().operator->();
  es.get<TrackerRecoGeometryRecord>().get(tracker);

  if (dbe_) {
    ibooker.setCurrentFolder("TrackerDigisV/TrackerDigis/Pixel");

    meDigiMultiLayer1Ring1_ = ibooker.book1D("digimulti_layer1ring1", "Digi Multiplicity ", 30, 0., 30.);
    meDigiMultiLayer1Ring2_ = ibooker.book1D("digimulti_layer1ring2", "Digi Multiplicity ", 30, 0., 30.);
    meDigiMultiLayer1Ring3_ = ibooker.book1D("digimulti_layer1ring3", "Digi Multiplicity ", 30, 0., 30.);
    meDigiMultiLayer1Ring4_ = ibooker.book1D("digimulti_layer1ring4", "Digi Multiplicity ", 30, 0., 30.);
    meDigiMultiLayer1Ring5_ = ibooker.book1D("digimulti_layer1ring5", "Digi Multiplicity ", 30, 0., 30.);
    meDigiMultiLayer1Ring6_ = ibooker.book1D("digimulti_layer1ring6", "Digi Multiplicity ", 30, 0., 30.);
    meDigiMultiLayer1Ring7_ = ibooker.book1D("digimulti_layer1ring7", "Digi Multiplicity ", 30, 0., 30.);
    meDigiMultiLayer1Ring8_ = ibooker.book1D("digimulti_layer1ring8", "Digi Multiplicity ", 30, 0., 30.);

    meDigiMultiLayer2Ring1_ = ibooker.book1D("digimulti_layer2ring1", "Digi Multiplicity ", 30, 0., 30.);
    meDigiMultiLayer2Ring2_ = ibooker.book1D("digimulti_layer2ring2", "Digi Multiplicity ", 30, 0., 30.);
    meDigiMultiLayer2Ring3_ = ibooker.book1D("digimulti_layer2ring3", "Digi Multiplicity ", 30, 0., 30.);
    meDigiMultiLayer2Ring4_ = ibooker.book1D("digimulti_layer2ring4", "Digi Multiplicity ", 30, 0., 30.);
    meDigiMultiLayer2Ring5_ = ibooker.book1D("digimulti_layer2ring5", "Digi Multiplicity ", 30, 0., 30.);
    meDigiMultiLayer2Ring6_ = ibooker.book1D("digimulti_layer2ring6", "Digi Multiplicity ", 30, 0., 30.);
    meDigiMultiLayer2Ring7_ = ibooker.book1D("digimulti_layer2ring7", "Digi Multiplicity ", 30, 0., 30.);
    meDigiMultiLayer2Ring8_ = ibooker.book1D("digimulti_layer2ring8", "Digi Multiplicity ", 30, 0., 30.);

    meDigiMultiLayer3Ring1_ = ibooker.book1D("digimulti_layer3ring1", "Digi Multiplicity ", 30, 0., 30.);
    meDigiMultiLayer3Ring2_ = ibooker.book1D("digimulti_layer3ring2", "Digi Multiplicity ", 30, 0., 30.);
    meDigiMultiLayer3Ring3_ = ibooker.book1D("digimulti_layer3ring3", "Digi Multiplicity ", 30, 0., 30.);
    meDigiMultiLayer3Ring4_ = ibooker.book1D("digimulti_layer3ring4", "Digi Multiplicity ", 30, 0., 30.);
    meDigiMultiLayer3Ring5_ = ibooker.book1D("digimulti_layer3ring5", "Digi Multiplicity ", 30, 0., 30.);
    meDigiMultiLayer3Ring6_ = ibooker.book1D("digimulti_layer3ring6", "Digi Multiplicity ", 30, 0., 30.);
    meDigiMultiLayer3Ring7_ = ibooker.book1D("digimulti_layer3ring7", "Digi Multiplicity ", 30, 0., 30.);
    meDigiMultiLayer3Ring8_ = ibooker.book1D("digimulti_layer3ring8", "Digi Multiplicity ", 30, 0., 30.);

    meAdcLayer1Ring1_ = ibooker.book1D("adc_layer1ring1", "Digi charge", 50, 0., 300.);
    meAdcLayer1Ring2_ = ibooker.book1D("adc_layer1ring2", "Digi charge", 50, 0., 300.);
    meAdcLayer1Ring3_ = ibooker.book1D("adc_layer1ring3", "Digi charge", 50, 0., 300.);
    meAdcLayer1Ring4_ = ibooker.book1D("adc_layer1ring4", "Digi charge", 50, 0., 300.);
    meAdcLayer1Ring5_ = ibooker.book1D("adc_layer1ring5", "Digi charge", 50, 0., 300.);
    meAdcLayer1Ring6_ = ibooker.book1D("adc_layer1ring6", "Digi charge", 50, 0., 300.);
    meAdcLayer1Ring7_ = ibooker.book1D("adc_layer1ring7", "Digi charge", 50, 0., 300.);
    meAdcLayer1Ring8_ = ibooker.book1D("adc_layer1ring8", "Digi charge", 50, 0., 300.);

    meRowLayer1Ring1_ = ibooker.book1D("row_layer1ring1", "Digi row", 50, 0., 200.);
    meRowLayer1Ring2_ = ibooker.book1D("row_layer1ring2", "Digi row", 50, 0., 200.);
    meRowLayer1Ring3_ = ibooker.book1D("row_layer1ring3", "Digi row", 50, 0., 200.);
    meRowLayer1Ring4_ = ibooker.book1D("row_layer1ring4", "Digi row", 50, 0., 200.);
    meRowLayer1Ring5_ = ibooker.book1D("row_layer1ring5", "Digi row", 50, 0., 200.);
    meRowLayer1Ring6_ = ibooker.book1D("row_layer1ring6", "Digi row", 50, 0., 200.);
    meRowLayer1Ring7_ = ibooker.book1D("row_layer1ring7", "Digi row", 50, 0., 200.);
    meRowLayer1Ring8_ = ibooker.book1D("row_layer1ring8", "Digi row", 50, 0., 200.);

    meColLayer1Ring1_ = ibooker.book1D("col_layer1ring1", "Digi column", 50, 0., 500.);
    meColLayer1Ring2_ = ibooker.book1D("col_layer1ring2", "Digi column", 50, 0., 500.);
    meColLayer1Ring3_ = ibooker.book1D("col_layer1ring3", "Digi column", 50, 0., 500.);
    meColLayer1Ring4_ = ibooker.book1D("col_layer1ring4", "Digi column", 50, 0., 500.);
    meColLayer1Ring5_ = ibooker.book1D("col_layer1ring5", "Digi column", 50, 0., 500.);
    meColLayer1Ring6_ = ibooker.book1D("col_layer1ring6", "Digi column", 50, 0., 500.);
    meColLayer1Ring7_ = ibooker.book1D("col_layer1ring7", "Digi column", 50, 0., 500.);
    meColLayer1Ring8_ = ibooker.book1D("col_layer1ring8", "Digi column", 50, 0., 500.);

    meAdcLayer2Ring1_ = ibooker.book1D("adc_layer2ring1", "Digi charge", 50, 0., 300.);
    meAdcLayer2Ring2_ = ibooker.book1D("adc_layer2ring2", "Digi charge", 50, 0., 300.);
    meAdcLayer2Ring3_ = ibooker.book1D("adc_layer2ring3", "Digi charge", 50, 0., 300.);
    meAdcLayer2Ring4_ = ibooker.book1D("adc_layer2ring4", "Digi charge", 50, 0., 300.);
    meAdcLayer2Ring5_ = ibooker.book1D("adc_layer2ring5", "Digi charge", 50, 0., 300.);
    meAdcLayer2Ring6_ = ibooker.book1D("adc_layer2ring6", "Digi charge", 50, 0., 300.);
    meAdcLayer2Ring7_ = ibooker.book1D("adc_layer2ring7", "Digi charge", 50, 0., 300.);
    meAdcLayer2Ring8_ = ibooker.book1D("adc_layer2ring8", "Digi charge", 50, 0., 300.);

    meRowLayer2Ring1_ = ibooker.book1D("row_layer2ring1", "Digi row", 50, 0., 200.);
    meRowLayer2Ring2_ = ibooker.book1D("row_layer2ring2", "Digi row", 50, 0., 200.);
    meRowLayer2Ring3_ = ibooker.book1D("row_layer2ring3", "Digi row", 50, 0., 200.);
    meRowLayer2Ring4_ = ibooker.book1D("row_layer2ring4", "Digi row", 50, 0., 200.);
    meRowLayer2Ring5_ = ibooker.book1D("row_layer2ring5", "Digi row", 50, 0., 200.);
    meRowLayer2Ring6_ = ibooker.book1D("row_layer2ring6", "Digi row", 50, 0., 200.);
    meRowLayer2Ring7_ = ibooker.book1D("row_layer2ring7", "Digi row", 50, 0., 200.);
    meRowLayer2Ring8_ = ibooker.book1D("row_layer2ring8", "Digi row", 50, 0., 200.);

    meColLayer2Ring1_ = ibooker.book1D("col_layer2ring1", "Digi column", 50, 0., 500.);
    meColLayer2Ring2_ = ibooker.book1D("col_layer2ring2", "Digi column", 50, 0., 500.);
    meColLayer2Ring3_ = ibooker.book1D("col_layer2ring3", "Digi column", 50, 0., 500.);
    meColLayer2Ring4_ = ibooker.book1D("col_layer2ring4", "Digi column", 50, 0., 500.);
    meColLayer2Ring5_ = ibooker.book1D("col_layer2ring5", "Digi column", 50, 0., 500.);
    meColLayer2Ring6_ = ibooker.book1D("col_layer2ring6", "Digi column", 50, 0., 500.);
    meColLayer2Ring7_ = ibooker.book1D("col_layer2ring7", "Digi column", 50, 0., 500.);
    meColLayer2Ring8_ = ibooker.book1D("col_layer2ring8", "Digi column", 50, 0., 500.);

    meAdcLayer3Ring1_ = ibooker.book1D("adc_layer3ring1", "Digi charge", 50, 0., 300.);
    meAdcLayer3Ring2_ = ibooker.book1D("adc_layer3ring2", "Digi charge", 50, 0., 300.);
    meAdcLayer3Ring3_ = ibooker.book1D("adc_layer3ring3", "Digi charge", 50, 0., 300.);
    meAdcLayer3Ring4_ = ibooker.book1D("adc_layer3ring4", "Digi charge", 50, 0., 300.);
    meAdcLayer3Ring5_ = ibooker.book1D("adc_layer3ring5", "Digi charge", 50, 0., 300.);
    meAdcLayer3Ring6_ = ibooker.book1D("adc_layer3ring6", "Digi charge", 50, 0., 300.);
    meAdcLayer3Ring7_ = ibooker.book1D("adc_layer3ring7", "Digi charge", 50, 0., 300.);
    meAdcLayer3Ring8_ = ibooker.book1D("adc_layer3ring8", "Digi charge", 50, 0., 300.);

    meRowLayer3Ring1_ = ibooker.book1D("row_layer3ring1", "Digi row", 50, 0., 200.);
    meRowLayer3Ring2_ = ibooker.book1D("row_layer3ring2", "Digi row", 50, 0., 200.);
    meRowLayer3Ring3_ = ibooker.book1D("row_layer3ring3", "Digi row", 50, 0., 200.);
    meRowLayer3Ring4_ = ibooker.book1D("row_layer3ring4", "Digi row", 50, 0., 200.);
    meRowLayer3Ring5_ = ibooker.book1D("row_layer3ring5", "Digi row", 50, 0., 200.);
    meRowLayer3Ring6_ = ibooker.book1D("row_layer3ring6", "Digi row", 50, 0., 200.);
    meRowLayer3Ring7_ = ibooker.book1D("row_layer3ring7", "Digi row", 50, 0., 200.);
    meRowLayer3Ring8_ = ibooker.book1D("row_layer3ring8", "Digi row", 50, 0., 200.);

    meColLayer3Ring1_ = ibooker.book1D("col_layer3ring1", "Digi column", 50, 0., 500.);
    meColLayer3Ring2_ = ibooker.book1D("col_layer3ring2", "Digi column", 50, 0., 500.);
    meColLayer3Ring3_ = ibooker.book1D("col_layer3ring3", "Digi column", 50, 0., 500.);
    meColLayer3Ring4_ = ibooker.book1D("col_layer3ring4", "Digi column", 50, 0., 500.);
    meColLayer3Ring5_ = ibooker.book1D("col_layer3ring5", "Digi column", 50, 0., 500.);
    meColLayer3Ring6_ = ibooker.book1D("col_layer3ring6", "Digi column", 50, 0., 500.);
    meColLayer3Ring7_ = ibooker.book1D("col_layer3ring7", "Digi column", 50, 0., 500.);
    meColLayer3Ring8_ = ibooker.book1D("col_layer3ring8", "Digi column", 50, 0., 500.);

    meDigiMultiLayer1Ladders_ =
        ibooker.bookProfile("digi_layer1_ladders", "Digi Num. per ladder", 22, 0.0, 21.0, 100, 0.0, 100);
    meDigiMultiLayer2Ladders_ =
        ibooker.bookProfile("digi_layer2_ladders", "Digi Num. per ladder", 34, 0.0, 32.0, 100, 0.0, 100);
    meDigiMultiLayer3Ladders_ =
        ibooker.bookProfile("digi_layer3_ladders", "Digi Num. per ladder", 46, 0.0, 45.0, 100, 0.0, 100);

    // Forward Pixel
    /* ZMinus Side 1st Disk */
    meAdcZmDisk1Panel1Plaq1_ = ibooker.book1D("adc_zm_disk1_panel1_plaq1", "Digi charge", 50, 0., 300.);
    meAdcZmDisk1Panel1Plaq2_ = ibooker.book1D("adc_zm_disk1_panel1_plaq2", "Digi charge", 50, 0., 300.);
    meAdcZmDisk1Panel1Plaq3_ = ibooker.book1D("adc_zm_disk1_panel1_plaq3", "Digi charge", 50, 0., 300.);
    meAdcZmDisk1Panel1Plaq4_ = ibooker.book1D("adc_zm_disk1_panel1_plaq4", "Digi charge", 50, 0., 300.);
    meAdcZmDisk1Panel2Plaq1_ = ibooker.book1D("adc_zm_disk1_panel2_plaq1", "Digi charge", 50, 0., 300.);
    meAdcZmDisk1Panel2Plaq2_ = ibooker.book1D("adc_zm_disk1_panel2_plaq2", "Digi charge", 50, 0., 300.);
    meAdcZmDisk1Panel2Plaq3_ = ibooker.book1D("adc_zm_disk1_panel2_plaq3", "Digi charge", 50, 0., 300.);

    meRowZmDisk1Panel1Plaq1_ = ibooker.book1D("row_zm_disk1_panel1_plaq1", "Digi row", 50, 0., 100.);
    meRowZmDisk1Panel1Plaq2_ = ibooker.book1D("row_zm_disk1_panel1_plaq2", "Digi row", 50, 0., 200.);
    meRowZmDisk1Panel1Plaq3_ = ibooker.book1D("row_zm_disk1_panel1_plaq3", "Digi row", 50, 0., 200.);
    meRowZmDisk1Panel1Plaq4_ = ibooker.book1D("row_zm_disk1_panel1_plaq4", "Digi row", 50, 0., 100.);
    meRowZmDisk1Panel2Plaq1_ = ibooker.book1D("row_zm_disk1_panel2_plaq1", "Digi row", 50, 0., 200.);
    meRowZmDisk1Panel2Plaq2_ = ibooker.book1D("row_zm_disk1_panel2_plaq2", "Digi row", 50, 0., 200.);
    meRowZmDisk1Panel2Plaq3_ = ibooker.book1D("row_zm_disk1_panel2_plaq3", "Digi row", 50, 0., 200.);

    meColZmDisk1Panel1Plaq1_ = ibooker.book1D("col_zm_disk1_panel1_plaq1", "Digi column", 50, 0., 150.);
    meColZmDisk1Panel1Plaq2_ = ibooker.book1D("col_zm_disk1_panel1_plaq2", "Digi column", 50, 0., 200.);
    meColZmDisk1Panel1Plaq3_ = ibooker.book1D("col_zm_disk1_panel1_plaq3", "Digi column", 50, 0., 250.);
    meColZmDisk1Panel1Plaq4_ = ibooker.book1D("col_zm_disk1_panel1_plaq4", "Digi column", 50, 0., 300.);
    meColZmDisk1Panel2Plaq1_ = ibooker.book1D("col_zm_disk1_panel2_plaq1", "Digi column", 50, 0., 200.);
    meColZmDisk1Panel2Plaq2_ = ibooker.book1D("col_zm_disk1_panel2_plaq2", "Digi column", 50, 0., 250.);
    meColZmDisk1Panel2Plaq3_ = ibooker.book1D("col_zm_disk1_panel2_plaq3", "Digi column", 50, 0., 300.);
    meNdigiZmDisk1PerPanel1_ =
        ibooker.book1D("digi_zm_disk1_panel1", "Digi Num. Panel1 Of 1st Disk In ZMinus Side ", 30, 0., 30.);
    meNdigiZmDisk1PerPanel2_ =
        ibooker.book1D("digi_zm_disk1_panel2", "Digi Num. Panel2 Of 1st Disk In ZMinus Side ", 30, 0., 30.);

    /* ZMius Side 2nd disk */
    meAdcZmDisk2Panel1Plaq1_ = ibooker.book1D("adc_zm_disk2_panel1_plaq1", "Digi charge", 50, 0., 300.);
    meAdcZmDisk2Panel1Plaq2_ = ibooker.book1D("adc_zm_disk2_panel1_plaq2", "Digi charge", 50, 0., 300.);
    meAdcZmDisk2Panel1Plaq3_ = ibooker.book1D("adc_zm_disk2_panel1_plaq3", "Digi charge", 50, 0., 300.);
    meAdcZmDisk2Panel1Plaq4_ = ibooker.book1D("adc_zm_disk2_panel1_plaq4", "Digi charge", 50, 0., 300.);
    meAdcZmDisk2Panel2Plaq1_ = ibooker.book1D("adc_zm_disk2_panel2_plaq1", "Digi charge", 50, 0., 300.);
    meAdcZmDisk2Panel2Plaq2_ = ibooker.book1D("adc_zm_disk2_panel2_plaq2", "Digi charge", 50, 0., 300.);
    meAdcZmDisk2Panel2Plaq3_ = ibooker.book1D("adc_zm_disk2_panel2_plaq3", "Digi charge", 50, 0., 300.);

    meRowZmDisk2Panel1Plaq1_ = ibooker.book1D("row_zm_disk2_panel1_plaq1", "Digi row", 50, 0., 100.);
    meRowZmDisk2Panel1Plaq2_ = ibooker.book1D("row_zm_disk2_panel1_plaq2", "Digi row", 50, 0., 200.);
    meRowZmDisk2Panel1Plaq3_ = ibooker.book1D("row_zm_disk2_panel1_plaq3", "Digi row", 50, 0., 200.);
    meRowZmDisk2Panel1Plaq4_ = ibooker.book1D("row_zm_disk2_panel1_plaq4", "Digi row", 50, 0., 100.);
    meRowZmDisk2Panel2Plaq1_ = ibooker.book1D("row_zm_disk2_panel2_plaq1", "Digi row", 50, 0., 200.);
    meRowZmDisk2Panel2Plaq2_ = ibooker.book1D("row_zm_disk2_panel2_plaq2", "Digi row", 50, 0., 200.);
    meRowZmDisk2Panel2Plaq3_ = ibooker.book1D("row_zm_disk2_panel2_plaq3", "Digi row", 50, 0., 200.);

    meColZmDisk2Panel1Plaq1_ = ibooker.book1D("col_zm_disk2_panel1_plaq1", "Digi Column", 50, 0., 150.);
    meColZmDisk2Panel1Plaq2_ = ibooker.book1D("col_zm_disk2_panel1_plaq2", "Digi Column", 50, 0., 200.);
    meColZmDisk2Panel1Plaq3_ = ibooker.book1D("col_zm_disk2_panel1_plaq3", "Digi Column", 50, 0., 250.);
    meColZmDisk2Panel1Plaq4_ = ibooker.book1D("col_zm_disk2_panel1_plaq4", "Digi Column", 50, 0., 300.);
    meColZmDisk2Panel2Plaq1_ = ibooker.book1D("col_zm_disk2_panel2_plaq1", "Digi Column", 50, 0., 200.);
    meColZmDisk2Panel2Plaq2_ = ibooker.book1D("col_zm_disk2_panel2_plaq2", "Digi Column", 50, 0., 250.);
    meColZmDisk2Panel2Plaq3_ = ibooker.book1D("col_zm_disk2_panel2_plaq3", "Digi Column", 50, 0., 300.);
    meNdigiZmDisk2PerPanel1_ =
        ibooker.book1D("digi_zm_disk2_panel1", "Digi Num. Panel1 Of 2nd Disk In ZMinus Side ", 30, 0., 30.);
    meNdigiZmDisk2PerPanel2_ =
        ibooker.book1D("digi_zm_disk2_panel2", "Digi Num. Panel2 Of 2nd Disk In ZMinus Side ", 30, 0., 30.);

    /* ZPlus Side 1st Disk */
    meAdcZpDisk1Panel1Plaq1_ = ibooker.book1D("adc_zp_disk1_panel1_plaq1", "Digi charge", 50, 0., 300.);
    meAdcZpDisk1Panel1Plaq2_ = ibooker.book1D("adc_zp_disk1_panel1_plaq2", "Digi charge", 50, 0., 300.);
    meAdcZpDisk1Panel1Plaq3_ = ibooker.book1D("adc_zp_disk1_panel1_plaq3", "Digi charge", 50, 0., 300.);
    meAdcZpDisk1Panel1Plaq4_ = ibooker.book1D("adc_zp_disk1_panel1_plaq4", "Digi charge", 50, 0., 300.);
    meAdcZpDisk1Panel2Plaq1_ = ibooker.book1D("adc_zp_disk1_panel2_plaq1", "Digi charge", 50, 0., 300.);
    meAdcZpDisk1Panel2Plaq2_ = ibooker.book1D("adc_zp_disk1_panel2_plaq2", "Digi charge", 50, 0., 300.);
    meAdcZpDisk1Panel2Plaq3_ = ibooker.book1D("adc_zp_disk1_panel2_plaq3", "Digi charge", 50, 0., 300.);

    meRowZpDisk1Panel1Plaq1_ = ibooker.book1D("row_zp_disk1_panel1_plaq1", "Digi row", 50, 0., 100.);
    meRowZpDisk1Panel1Plaq2_ = ibooker.book1D("row_zp_disk1_panel1_plaq2", "Digi row", 50, 0., 200.);
    meRowZpDisk1Panel1Plaq3_ = ibooker.book1D("row_zp_disk1_panel1_plaq3", "Digi row", 50, 0., 200.);
    meRowZpDisk1Panel1Plaq4_ = ibooker.book1D("row_zp_disk1_panel1_plaq4", "Digi row", 50, 0., 100.);
    meRowZpDisk1Panel2Plaq1_ = ibooker.book1D("row_zp_disk1_panel2_plaq1", "Digi row", 50, 0., 200.);
    meRowZpDisk1Panel2Plaq2_ = ibooker.book1D("row_zp_disk1_panel2_plaq2", "Digi row", 50, 0., 200.);
    meRowZpDisk1Panel2Plaq3_ = ibooker.book1D("row_zp_disk1_panel2_plaq3", "Digi row", 50, 0., 200.);

    meColZpDisk1Panel1Plaq1_ = ibooker.book1D("col_zp_disk1_panel1_plaq1", "Digi Column", 50, 0., 150.);
    meColZpDisk1Panel1Plaq2_ = ibooker.book1D("col_zp_disk1_panel1_plaq2", "Digi column", 50, 0., 200.);
    meColZpDisk1Panel1Plaq3_ = ibooker.book1D("col_zp_disk1_panel1_plaq3", "Digi column", 50, 0., 250.);
    meColZpDisk1Panel1Plaq4_ = ibooker.book1D("col_zp_disk1_panel1_plaq4", "Digi column", 50, 0., 300.);
    meColZpDisk1Panel2Plaq1_ = ibooker.book1D("col_zp_disk1_panel2_plaq1", "Digi column", 50, 0., 200.);
    meColZpDisk1Panel2Plaq2_ = ibooker.book1D("col_zp_disk1_panel2_plaq2", "Digi column", 50, 0., 250.);
    meColZpDisk1Panel2Plaq3_ = ibooker.book1D("col_zp_disk1_panel2_plaq3", "Digi column", 50, 0., 300.);
    meNdigiZpDisk1PerPanel1_ =
        ibooker.book1D("digi_zp_disk1_panel1", "Digi Num. Panel1 Of 1st Disk In ZPlus Side ", 30, 0., 30.);
    meNdigiZpDisk1PerPanel2_ =
        ibooker.book1D("digi_zp_disk1_panel2", "Digi Num. Panel2 Of 1st Disk In ZPlus Side ", 30, 0., 30.);

    /* ZPlus Side 2nd disk */
    meAdcZpDisk2Panel1Plaq1_ = ibooker.book1D("adc_zp_disk2_panel1_plaq1", "Digi charge", 50, 0., 300.);
    meAdcZpDisk2Panel1Plaq2_ = ibooker.book1D("adc_zp_disk2_panel1_plaq2", "Digi charge", 50, 0., 300.);
    meAdcZpDisk2Panel1Plaq3_ = ibooker.book1D("adc_zp_disk2_panel1_plaq3", "Digi charge", 50, 0., 300.);
    meAdcZpDisk2Panel1Plaq4_ = ibooker.book1D("adc_zp_disk2_panel1_plaq4", "Digi charge", 50, 0., 300.);
    meAdcZpDisk2Panel2Plaq1_ = ibooker.book1D("adc_zp_disk2_panel2_plaq1", "Digi charge", 50, 0., 300.);
    meAdcZpDisk2Panel2Plaq2_ = ibooker.book1D("adc_zp_disk2_panel2_plaq2", "Digi charge", 50, 0., 300.);
    meAdcZpDisk2Panel2Plaq3_ = ibooker.book1D("adc_zp_disk2_panel2_plaq3", "Digi charge", 50, 0., 300.);

    meRowZpDisk2Panel1Plaq1_ = ibooker.book1D("row_zp_disk2_panel1_plaq1", "Digi row", 10, 0., 100.);
    meRowZpDisk2Panel1Plaq2_ = ibooker.book1D("row_zp_disk2_panel1_plaq2", "Digi row", 10, 0., 200.);
    meRowZpDisk2Panel1Plaq3_ = ibooker.book1D("row_zp_disk2_panel1_plaq3", "Digi row", 10, 0., 200.);
    meRowZpDisk2Panel1Plaq4_ = ibooker.book1D("row_zp_disk2_panel1_plaq4", "Digi row", 10, 0., 100.);
    meRowZpDisk2Panel2Plaq1_ = ibooker.book1D("row_zp_disk2_panel2_plaq1", "Digi row", 10, 0., 200.);
    meRowZpDisk2Panel2Plaq2_ = ibooker.book1D("row_zp_disk2_panel2_plaq2", "Digi row", 10, 0., 200.);
    meRowZpDisk2Panel2Plaq3_ = ibooker.book1D("row_zp_disk2_panel2_plaq3", "Digi row", 10, 0., 200.);

    meColZpDisk2Panel1Plaq1_ = ibooker.book1D("col_zp_disk2_panel1_plaq1", "Digi column", 50, 0., 150.);
    meColZpDisk2Panel1Plaq2_ = ibooker.book1D("col_zp_disk2_panel1_plaq2", "Digi column", 50, 0., 200.);
    meColZpDisk2Panel1Plaq3_ = ibooker.book1D("col_zp_disk2_panel1_plaq3", "Digi column", 50, 0., 250.);
    meColZpDisk2Panel1Plaq4_ = ibooker.book1D("col_zp_disk2_panel1_plaq4", "Digi column", 50, 0., 300.);
    meColZpDisk2Panel2Plaq1_ = ibooker.book1D("col_zp_disk2_panel2_plaq1", "Digi column", 50, 0., 200.);
    meColZpDisk2Panel2Plaq2_ = ibooker.book1D("col_zp_disk2_panel2_plaq2", "Digi column", 50, 0., 250.);
    meColZpDisk2Panel2Plaq3_ = ibooker.book1D("col_zp_disk2_panel2_plaq3", "Digi column", 50, 0., 300.);
    meNdigiZpDisk2PerPanel1_ =
        ibooker.book1D("digi_zp_disk2_panel1", "Digi Num. Panel1 Of 2nd Disk In ZPlus Side ", 30, 0., 30.);
    meNdigiZpDisk2PerPanel2_ =
        ibooker.book1D("digi_zp_disk2_panel2", "Digi Num. Panel2 Of 2nd Disk In ZPlus Side ", 30, 0., 30.);
  }
}

void SiPixelDigiValid::endJob() {
  // Save histos in a file only in standalone mode
  if (runStandalone && !outputFile_.empty() && dbe_) {
    dbe_->save(outputFile_);
  }
}

void SiPixelDigiValid::analyze(const edm::Event &e, const edm::EventSetup &c) {
  // Retrieve tracker topology from geometry
  edm::ESHandle<TrackerTopology> tTopoHandle;
  c.get<TrackerTopologyRcd>().get(tTopoHandle);
  const TrackerTopology *const tTopo = tTopoHandle.product();
  // Number of blades and ladders.
  // TODO: other Geometry-Dependent quantities, e.g. num layers.
  int nblades = tracker->posPixelForwardLayers()[0]->components().size();
  int nladders1 = tracker->pixelBarrelLayers()[0]->components().size();
  int nladders2 = tracker->pixelBarrelLayers()[1]->components().size();
  int nladders3 = tracker->pixelBarrelLayers()[2]->components().size();

  int ndigiperRingLayer1[8];
  int ndigiperRingLayer2[8];
  int ndigiperRingLayer3[8];
  for (int i = 0; i < 8; i++) {
    ndigiperRingLayer1[i] = 0;
    ndigiperRingLayer2[i] = 0;
    ndigiperRingLayer3[i] = 0;
  }

  int ndigiZpDisk1PerPanel1[nblades];
  int ndigiZpDisk1PerPanel2[nblades];
  int ndigiZpDisk2PerPanel1[nblades];
  int ndigiZpDisk2PerPanel2[nblades];
  int ndigiZmDisk1PerPanel1[nblades];
  int ndigiZmDisk1PerPanel2[nblades];
  int ndigiZmDisk2PerPanel1[nblades];
  int ndigiZmDisk2PerPanel2[nblades];

  for (int i = 0; i < nblades; i++) {
    ndigiZpDisk1PerPanel1[i] = 0;
    ndigiZpDisk1PerPanel2[i] = 0;
    ndigiZpDisk2PerPanel1[i] = 0;
    ndigiZpDisk2PerPanel2[i] = 0;
    ndigiZmDisk1PerPanel1[i] = 0;
    ndigiZmDisk1PerPanel2[i] = 0;
    ndigiZmDisk2PerPanel1[i] = 0;
    ndigiZmDisk2PerPanel2[i] = 0;
  }

  int ndigilayer1ladders[nladders1];
  int ndigilayer2ladders[nladders2];
  int ndigilayer3ladders[nladders3];

  for (int i = 0; i < nladders1; i++) {
    ndigilayer1ladders[i] = 0;
  }

  for (int i = 0; i < nladders2; i++) {
    ndigilayer2ladders[i] = 0;
  }

  for (int i = 0; i < nladders3; i++) {
    ndigilayer3ladders[i] = 0;
  }

  // LogInfo("EventInfo") << " Run = " << e.id().run() << " Event = " <<
  // e.id().event();

  edm::ESHandle<TrackerGeometry> tracker;
  c.get<TrackerDigiGeometryRecord>().get(tracker);

  // string digiProducer = "siPixelDigis";
  edm::Handle<edm::DetSetVector<PixelDigi>> pixelDigis;
  e.getByToken(edmDetSetVector_PixelDigi_Token_, pixelDigis);

  edm::DetSetVector<PixelDigi>::const_iterator DSViter = pixelDigis->begin();
  for (; DSViter != pixelDigis->end(); DSViter++) {
    unsigned int id = DSViter->id;
    DetId detId(id);
    edm::DetSet<PixelDigi>::const_iterator begin = DSViter->data.begin();
    edm::DetSet<PixelDigi>::const_iterator end = DSViter->data.end();
    edm::DetSet<PixelDigi>::const_iterator iter;

    if (detId.subdetId() == PixelSubdetector::PixelBarrel) {
      unsigned int layer = tTopo->pxbLayer(id);    // Layer:1,2,3.
      unsigned int ladder = tTopo->pxbLadder(id);  // Ladeer: 1-20, 32, 44.
      unsigned int zindex = tTopo->pxbModule(id);  // Z-index: 1-8.
      // LogInfo("SiPixelDigiValid")<<"Barrel:: Layer="<<layer<<"
      // Ladder="<<ladder<<" zindex="<<zindex;
      for (iter = begin; iter != end; iter++) {
        if (layer == 1) {
          ++ndigilayer1ladders[ladder - 1];
          ++ndigiperRingLayer1[zindex - 1];
          if (zindex == 1) {
            meAdcLayer1Ring1_->Fill((*iter).adc());
            meRowLayer1Ring1_->Fill((*iter).row());
            meColLayer1Ring1_->Fill((*iter).column());
          }
          if (zindex == 2) {
            meAdcLayer1Ring2_->Fill((*iter).adc());
            meRowLayer1Ring2_->Fill((*iter).row());
            meColLayer1Ring2_->Fill((*iter).column());
          }

          if (zindex == 3) {
            meAdcLayer1Ring3_->Fill((*iter).adc());
            meRowLayer1Ring3_->Fill((*iter).row());
            meColLayer1Ring3_->Fill((*iter).column());
          }

          if (zindex == 4) {
            meAdcLayer1Ring4_->Fill((*iter).adc());
            meRowLayer1Ring4_->Fill((*iter).row());
            meColLayer1Ring4_->Fill((*iter).column());
          }

          if (zindex == 5) {
            meAdcLayer1Ring5_->Fill((*iter).adc());
            meRowLayer1Ring5_->Fill((*iter).row());
            meColLayer1Ring5_->Fill((*iter).column());
          }

          if (zindex == 6) {
            meAdcLayer1Ring6_->Fill((*iter).adc());
            meRowLayer1Ring6_->Fill((*iter).row());
            meColLayer1Ring6_->Fill((*iter).column());
          }

          if (zindex == 7) {
            meAdcLayer1Ring7_->Fill((*iter).adc());
            meRowLayer1Ring7_->Fill((*iter).row());
            meColLayer1Ring7_->Fill((*iter).column());
          }
          if (zindex == 8) {
            meAdcLayer1Ring8_->Fill((*iter).adc());
            meRowLayer1Ring8_->Fill((*iter).row());
            meColLayer1Ring8_->Fill((*iter).column());
          }
        }
        if (layer == 2) {
          ++ndigilayer2ladders[ladder - 1];
          ++ndigiperRingLayer2[zindex - 1];
          if (zindex == 1) {
            meAdcLayer2Ring1_->Fill((*iter).adc());
            meRowLayer2Ring1_->Fill((*iter).row());
            meColLayer2Ring1_->Fill((*iter).column());
          }
          if (zindex == 2) {
            meAdcLayer2Ring2_->Fill((*iter).adc());
            meRowLayer2Ring2_->Fill((*iter).row());
            meColLayer2Ring2_->Fill((*iter).column());
          }

          if (zindex == 3) {
            meAdcLayer2Ring3_->Fill((*iter).adc());
            meRowLayer2Ring3_->Fill((*iter).row());
            meColLayer2Ring3_->Fill((*iter).column());
          }

          if (zindex == 4) {
            meAdcLayer2Ring4_->Fill((*iter).adc());
            meRowLayer2Ring4_->Fill((*iter).row());
            meColLayer2Ring4_->Fill((*iter).column());
          }

          if (zindex == 5) {
            meAdcLayer2Ring5_->Fill((*iter).adc());
            meRowLayer2Ring5_->Fill((*iter).row());
            meColLayer2Ring5_->Fill((*iter).column());
          }

          if (zindex == 6) {
            meAdcLayer2Ring6_->Fill((*iter).adc());
            meRowLayer2Ring6_->Fill((*iter).row());
            meColLayer2Ring6_->Fill((*iter).column());
          }

          if (zindex == 7) {
            meAdcLayer2Ring7_->Fill((*iter).adc());
            meRowLayer2Ring7_->Fill((*iter).row());
            meColLayer2Ring7_->Fill((*iter).column());
          }
          if (zindex == 8) {
            meAdcLayer2Ring8_->Fill((*iter).adc());
            meRowLayer2Ring8_->Fill((*iter).row());
            meColLayer2Ring8_->Fill((*iter).column());
          }
        }
        if (layer == 3) {
          ++ndigilayer3ladders[ladder - 1];
          ++ndigiperRingLayer3[zindex - 1];
          if (zindex == 1) {
            meAdcLayer3Ring1_->Fill((*iter).adc());
            meRowLayer3Ring1_->Fill((*iter).row());
            meColLayer3Ring1_->Fill((*iter).column());
          }
          if (zindex == 2) {
            meAdcLayer3Ring2_->Fill((*iter).adc());
            meRowLayer3Ring2_->Fill((*iter).row());
            meColLayer3Ring2_->Fill((*iter).column());
          }

          if (zindex == 3) {
            meAdcLayer3Ring3_->Fill((*iter).adc());
            meRowLayer3Ring3_->Fill((*iter).row());
            meColLayer3Ring3_->Fill((*iter).column());
          }

          if (zindex == 4) {
            meAdcLayer3Ring4_->Fill((*iter).adc());
            meRowLayer3Ring4_->Fill((*iter).row());
            meColLayer3Ring4_->Fill((*iter).column());
          }

          if (zindex == 5) {
            meAdcLayer3Ring5_->Fill((*iter).adc());
            meRowLayer3Ring5_->Fill((*iter).row());
            meColLayer3Ring5_->Fill((*iter).column());
          }

          if (zindex == 6) {
            meAdcLayer3Ring6_->Fill((*iter).adc());
            meRowLayer3Ring6_->Fill((*iter).row());
            meColLayer3Ring6_->Fill((*iter).column());
          }

          if (zindex == 7) {
            meAdcLayer3Ring7_->Fill((*iter).adc());
            meRowLayer3Ring7_->Fill((*iter).row());
            meColLayer3Ring7_->Fill((*iter).column());
          }
          if (zindex == 8) {
            meAdcLayer3Ring8_->Fill((*iter).adc());
            meRowLayer3Ring8_->Fill((*iter).row());
            meColLayer3Ring8_->Fill((*iter).column());
          }
        }
      }
    }

    if (detId.subdetId() == PixelSubdetector::PixelEndcap) {  // Endcap

      unsigned int side = tTopo->pxfSide(id);
      unsigned int disk = tTopo->pxfDisk(id);
      unsigned int blade = tTopo->pxfBlade(id);
      unsigned int panel = tTopo->pxfPanel(id);
      unsigned int mod = tTopo->pxfModule(id);
      // LogInfo("SiPixelDigiValid")<<"EndcaP="<<side<<" Disk="<<disk<<"
      // Blade="<<blade<<" Panel="<<panel<<" Module="<<mod;
      for (iter = begin; iter != end; iter++) {
        if (side == 1 && disk == 1 && panel == 1) {
          if (mod == 1) {
            meAdcZmDisk1Panel1Plaq1_->Fill((*iter).adc());
            meRowZmDisk1Panel1Plaq1_->Fill((*iter).row());
            meColZmDisk1Panel1Plaq1_->Fill((*iter).column());
          } else if (mod == 2) {
            meAdcZmDisk1Panel1Plaq2_->Fill((*iter).adc());
            meRowZmDisk1Panel1Plaq2_->Fill((*iter).row());
            meColZmDisk1Panel1Plaq2_->Fill((*iter).column());
          } else if (mod == 3) {
            meAdcZmDisk1Panel1Plaq3_->Fill((*iter).adc());
            meRowZmDisk1Panel1Plaq3_->Fill((*iter).row());
            meColZmDisk1Panel1Plaq3_->Fill((*iter).column());
          } else if (mod == 4) {
            meAdcZmDisk1Panel1Plaq4_->Fill((*iter).adc());
            meRowZmDisk1Panel1Plaq4_->Fill((*iter).row());
            meColZmDisk1Panel1Plaq4_->Fill((*iter).column());
          } else {
            // LogError("SiPixelDigiValid")<<" The number of module is Wrong";
          }
          ++ndigiZmDisk1PerPanel1[blade - 1];
        }

        if (side == 1 && disk == 1 && panel == 2) {
          if (mod == 1) {
            meAdcZmDisk1Panel2Plaq1_->Fill((*iter).adc());
            meRowZmDisk1Panel2Plaq1_->Fill((*iter).row());
            meColZmDisk1Panel2Plaq1_->Fill((*iter).column());
          } else if (mod == 2) {
            meAdcZmDisk1Panel2Plaq2_->Fill((*iter).adc());
            meRowZmDisk1Panel2Plaq2_->Fill((*iter).row());
            meColZmDisk1Panel2Plaq2_->Fill((*iter).column());
          } else if (mod == 3) {
            meAdcZmDisk1Panel2Plaq3_->Fill((*iter).adc());
            meRowZmDisk1Panel2Plaq3_->Fill((*iter).row());
            meColZmDisk1Panel2Plaq3_->Fill((*iter).column());
          } else {
            // LogError("SiPixelDigiValid")<<" The number of module is Wrong";
          }
          ++ndigiZmDisk1PerPanel2[blade - 1];
        }

        if (side == 1 && disk == 2 && panel == 1) {
          if (mod == 1) {
            meAdcZmDisk2Panel1Plaq1_->Fill((*iter).adc());
            meRowZmDisk2Panel1Plaq1_->Fill((*iter).row());
            meColZmDisk2Panel1Plaq1_->Fill((*iter).column());
          } else if (mod == 2) {
            meAdcZmDisk2Panel1Plaq2_->Fill((*iter).adc());
            meRowZmDisk2Panel1Plaq2_->Fill((*iter).row());
            meColZmDisk2Panel1Plaq2_->Fill((*iter).column());
          } else if (mod == 3) {
            meAdcZmDisk2Panel1Plaq3_->Fill((*iter).adc());
            meRowZmDisk2Panel1Plaq3_->Fill((*iter).row());
            meColZmDisk2Panel1Plaq3_->Fill((*iter).column());
          } else if (mod == 4) {
            meAdcZmDisk2Panel1Plaq4_->Fill((*iter).adc());
            meRowZmDisk2Panel1Plaq4_->Fill((*iter).row());
            meColZmDisk2Panel1Plaq4_->Fill((*iter).column());
          } else {
            // LogError("SiPixelDigiValid")<<" The number of module is Wrong";
          }
          ++ndigiZmDisk2PerPanel1[blade - 1];
        }

        if (side == 1 && disk == 2 && panel == 2) {
          if (mod == 1) {
            meAdcZmDisk2Panel2Plaq1_->Fill((*iter).adc());
            meRowZmDisk2Panel2Plaq1_->Fill((*iter).row());
            meColZmDisk2Panel2Plaq1_->Fill((*iter).column());
          } else if (mod == 2) {
            meAdcZmDisk2Panel2Plaq2_->Fill((*iter).adc());
            meRowZmDisk2Panel2Plaq2_->Fill((*iter).row());
            meColZmDisk2Panel2Plaq2_->Fill((*iter).column());
          } else if (mod == 3) {
            meAdcZmDisk2Panel2Plaq3_->Fill((*iter).adc());
            meRowZmDisk2Panel2Plaq3_->Fill((*iter).row());
            meColZmDisk2Panel2Plaq3_->Fill((*iter).column());
          } else {
            // LogError("SiPixelDigiValid")<<" The number of module is Wrong";
          }
          ++ndigiZmDisk2PerPanel2[blade - 1];
        }

        if (side == 2 && disk == 1 && panel == 1) {
          if (mod == 1) {
            meAdcZpDisk1Panel1Plaq1_->Fill((*iter).adc());
            meRowZpDisk1Panel1Plaq1_->Fill((*iter).row());
            meColZpDisk1Panel1Plaq1_->Fill((*iter).column());
          } else if (mod == 2) {
            meAdcZpDisk1Panel1Plaq2_->Fill((*iter).adc());
            meRowZpDisk1Panel1Plaq2_->Fill((*iter).row());
            meColZpDisk1Panel1Plaq2_->Fill((*iter).column());
          } else if (mod == 3) {
            meAdcZpDisk1Panel1Plaq3_->Fill((*iter).adc());
            meRowZpDisk1Panel1Plaq3_->Fill((*iter).row());
            meColZpDisk1Panel1Plaq3_->Fill((*iter).column());
          } else if (mod == 4) {
            meAdcZpDisk1Panel1Plaq4_->Fill((*iter).adc());
            meRowZpDisk1Panel1Plaq4_->Fill((*iter).row());
            meColZpDisk1Panel1Plaq4_->Fill((*iter).column());
          } else {
            // LogError("SiPixelDigiValid")<<" The number of module is Wrong";
          }
          ++ndigiZpDisk1PerPanel1[blade - 1];
        }

        if (side == 2 && disk == 1 && panel == 2) {
          if (mod == 1) {
            meAdcZpDisk1Panel2Plaq1_->Fill((*iter).adc());
            meRowZpDisk1Panel2Plaq1_->Fill((*iter).row());
            meColZpDisk1Panel2Plaq1_->Fill((*iter).column());
          } else if (mod == 2) {
            meAdcZpDisk1Panel2Plaq2_->Fill((*iter).adc());
            meRowZpDisk1Panel2Plaq2_->Fill((*iter).row());
            meColZpDisk1Panel2Plaq2_->Fill((*iter).column());
          } else if (mod == 3) {
            meAdcZpDisk1Panel2Plaq3_->Fill((*iter).adc());
            meRowZpDisk1Panel2Plaq3_->Fill((*iter).row());
            meColZpDisk1Panel2Plaq3_->Fill((*iter).column());
          } else {
            // LogError("SiPixelDigiValid")<<" The number of module is Wrong";
          }
          ++ndigiZpDisk1PerPanel2[blade - 1];
        }

        if (side == 2 && disk == 2 && panel == 1) {
          if (mod == 1) {
            meAdcZpDisk2Panel1Plaq1_->Fill((*iter).adc());
            meRowZpDisk2Panel1Plaq1_->Fill((*iter).row());
            meColZpDisk2Panel1Plaq1_->Fill((*iter).column());
          } else if (mod == 2) {
            meAdcZpDisk2Panel1Plaq2_->Fill((*iter).adc());
            meRowZpDisk2Panel1Plaq2_->Fill((*iter).row());
            meColZpDisk2Panel1Plaq2_->Fill((*iter).column());
          } else if (mod == 3) {
            meAdcZpDisk2Panel1Plaq3_->Fill((*iter).adc());
            meRowZpDisk2Panel1Plaq3_->Fill((*iter).row());
            meColZpDisk2Panel1Plaq3_->Fill((*iter).column());
          } else if (mod == 4) {
            meAdcZpDisk2Panel1Plaq4_->Fill((*iter).adc());
            meRowZpDisk2Panel1Plaq4_->Fill((*iter).row());
            meColZpDisk2Panel1Plaq4_->Fill((*iter).column());
          } else {
            // LogError("SiPixelDigiValid")<<" The number of module is Wrong";
          }
          ++ndigiZpDisk2PerPanel1[blade - 1];
        }

        if (side == 2 && disk == 2 && panel == 2) {
          if (mod == 1) {
            meAdcZpDisk2Panel2Plaq1_->Fill((*iter).adc());
            meRowZpDisk2Panel2Plaq1_->Fill((*iter).row());
            meColZpDisk2Panel2Plaq1_->Fill((*iter).column());
          } else if (mod == 2) {
            meAdcZpDisk2Panel2Plaq2_->Fill((*iter).adc());
            meRowZpDisk2Panel2Plaq2_->Fill((*iter).row());
            meColZpDisk2Panel2Plaq2_->Fill((*iter).column());
          } else if (mod == 3) {
            meAdcZpDisk2Panel2Plaq3_->Fill((*iter).adc());
            meRowZpDisk2Panel2Plaq3_->Fill((*iter).row());
            meColZpDisk2Panel2Plaq3_->Fill((*iter).column());
          } else {
            // LogError("SiPixelDigiValid")<<" The number of module is Wrong";
          }
          ++ndigiZpDisk2PerPanel2[blade - 1];
        }
      }  // iterating the digi

    }  // Endcap

  }  // end for loop

  meDigiMultiLayer1Ring1_->Fill(ndigiperRingLayer1[0]);
  meDigiMultiLayer1Ring2_->Fill(ndigiperRingLayer1[1]);
  meDigiMultiLayer1Ring3_->Fill(ndigiperRingLayer1[2]);
  meDigiMultiLayer1Ring4_->Fill(ndigiperRingLayer1[3]);
  meDigiMultiLayer1Ring5_->Fill(ndigiperRingLayer1[4]);
  meDigiMultiLayer1Ring6_->Fill(ndigiperRingLayer1[5]);
  meDigiMultiLayer1Ring7_->Fill(ndigiperRingLayer1[6]);
  meDigiMultiLayer1Ring8_->Fill(ndigiperRingLayer1[7]);

  meDigiMultiLayer2Ring1_->Fill(ndigiperRingLayer2[0]);
  meDigiMultiLayer2Ring2_->Fill(ndigiperRingLayer2[1]);
  meDigiMultiLayer2Ring3_->Fill(ndigiperRingLayer2[2]);
  meDigiMultiLayer2Ring4_->Fill(ndigiperRingLayer2[3]);
  meDigiMultiLayer2Ring5_->Fill(ndigiperRingLayer2[4]);
  meDigiMultiLayer2Ring6_->Fill(ndigiperRingLayer2[5]);
  meDigiMultiLayer2Ring7_->Fill(ndigiperRingLayer2[6]);
  meDigiMultiLayer2Ring8_->Fill(ndigiperRingLayer2[7]);

  meDigiMultiLayer3Ring1_->Fill(ndigiperRingLayer3[0]);
  meDigiMultiLayer3Ring2_->Fill(ndigiperRingLayer3[1]);
  meDigiMultiLayer3Ring3_->Fill(ndigiperRingLayer3[2]);
  meDigiMultiLayer3Ring4_->Fill(ndigiperRingLayer3[3]);
  meDigiMultiLayer3Ring5_->Fill(ndigiperRingLayer3[4]);
  meDigiMultiLayer3Ring6_->Fill(ndigiperRingLayer3[5]);
  meDigiMultiLayer3Ring7_->Fill(ndigiperRingLayer3[6]);
  meDigiMultiLayer3Ring8_->Fill(ndigiperRingLayer3[7]);

  for (int i = 0; i < nblades; i++) {
    meNdigiZmDisk1PerPanel1_->Fill(ndigiZmDisk1PerPanel1[i]);
    meNdigiZmDisk1PerPanel2_->Fill(ndigiZmDisk1PerPanel2[i]);
    meNdigiZmDisk2PerPanel1_->Fill(ndigiZmDisk2PerPanel1[i]);
    meNdigiZmDisk2PerPanel2_->Fill(ndigiZmDisk2PerPanel2[i]);
    meNdigiZpDisk1PerPanel1_->Fill(ndigiZpDisk1PerPanel1[i]);
    meNdigiZpDisk1PerPanel2_->Fill(ndigiZpDisk1PerPanel2[i]);
    meNdigiZpDisk2PerPanel1_->Fill(ndigiZpDisk2PerPanel1[i]);
    meNdigiZpDisk2PerPanel2_->Fill(ndigiZpDisk2PerPanel2[i]);
  }

  for (int i = 0; i < nladders1; i++) {
    meDigiMultiLayer1Ladders_->Fill(i + 1, ndigilayer1ladders[i]);
  }

  for (int i = 0; i < nladders2; i++) {
    meDigiMultiLayer2Ladders_->Fill(i + 1, ndigilayer2ladders[i]);
  }

  for (int i = 0; i < nladders3; i++) {
    meDigiMultiLayer3Ladders_->Fill(i + 1, ndigilayer3ladders[i]);
  }
}