SiPixelRecHitModule.cc

Go to the documentation of this file.

#include "DQM/SiPixelCommon/interface/SiPixelHistogramId.h"
#include "DQM/SiPixelMonitorRecHit/interface/SiPixelRecHitModule.h"
#include "DQMServices/Core/interface/DQMStore.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
// STL
#include <cstdlib>
#include <iostream>
#include <memory>
#include <string>
#include <vector>

// Data Formats
#include "DataFormats/DetId/interface/DetId.h"
#include "DataFormats/SiPixelDetId/interface/PixelBarrelName.h"
#include "DataFormats/SiPixelDetId/interface/PixelBarrelNameUpgrade.h"
#include "DataFormats/SiPixelDetId/interface/PixelEndcapName.h"
#include "DataFormats/SiPixelDetId/interface/PixelEndcapNameUpgrade.h"
#include "DataFormats/SiPixelDetId/interface/PixelSubdetector.h"
#include "DataFormats/TrackerCommon/interface/TrackerTopology.h"

#include "Geometry/Records/interface/TrackerDigiGeometryRecord.h"
#include "Geometry/Records/interface/TrackerTopologyRcd.h"
#include "Geometry/TrackerGeometryBuilder/interface/PixelGeomDetUnit.h"
#include "Geometry/TrackerGeometryBuilder/interface/TrackerGeometry.h"

#include "DataFormats/Common/interface/Handle.h"
#include "FWCore/Framework/interface/EDProducer.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/Framework/interface/Event.h"

#include "FWCore/Framework/interface/EventSetup.h"
// Constructors
//
SiPixelRecHitModule::SiPixelRecHitModule() : id_(0) {}
SiPixelRecHitModule::SiPixelRecHitModule(const uint32_t &id) : id_(id) {}

//
// Destructor
//
SiPixelRecHitModule::~SiPixelRecHitModule() {}
//
// Book histograms
//
void SiPixelRecHitModule::book(const edm::ParameterSet &iConfig,
                               DQMStore::IBooker &iBooker,
                               const edm::EventSetup &iSetup,
                               int type,
                               bool twoD,
                               bool reducedSet,
                               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) {
    // if (!isUpgrade) {
    isHalfModule = PixelBarrelName(DetId(id_), pTT, isUpgrade).isHalfModule();
    //} else if (isUpgrade) {
    //  isHalfModule = PixelBarrelNameUpgrade(DetId(id_)).isHalfModule();
    //}
  }

  std::string hid;
  // Get collection name and instantiate Histo Id builder
  edm::InputTag src = iConfig.getParameter<edm::InputTag>("src");
  // Get DQM interface

  if (type == 0) {
    SiPixelHistogramId *theHistogramId = new SiPixelHistogramId(src.label());
    if (!reducedSet) {
      if (twoD) {
        // XYPosition
        hid = theHistogramId->setHistoId("xypos", id_);
        meXYPos_ = iBooker.book2D(hid, "XY Position", 100, -1., 1, 100, -4, 4);
        meXYPos_->setAxisTitle("X Position", 1);
        meXYPos_->setAxisTitle("Y Position", 2);
      } else {
        // projections of XYPosition
        hid = theHistogramId->setHistoId("xypos", id_);
        meXYPos_px_ = iBooker.book1D(hid + "_px", "X Position", 100, -1., 1);
        meXYPos_px_->setAxisTitle("X Position", 1);
        meXYPos_py_ = iBooker.book1D(hid + "_py", "Y Position", 100, -4, 4);
        meXYPos_py_->setAxisTitle("Y Position", 1);
      }
    }
    hid = theHistogramId->setHistoId("ClustX", id_);
    meClustX_ = iBooker.book1D(hid, "RecHit X size", 10, 0., 10.);
    meClustX_->setAxisTitle("RecHit size X dimension", 1);
    hid = theHistogramId->setHistoId("ClustY", id_);
    meClustY_ = iBooker.book1D(hid, "RecHit Y size", 15, 0., 15.);
    meClustY_->setAxisTitle("RecHit size Y dimension", 1);

    hid = theHistogramId->setHistoId("ErrorX", id_);
    meErrorX_ = iBooker.book1D(hid, "RecHit error X", 100, 0., 0.02);
    meErrorX_->setAxisTitle("RecHit error X", 1);
    hid = theHistogramId->setHistoId("ErrorY", id_);
    meErrorY_ = iBooker.book1D(hid, "RecHit error Y", 100, 0., 0.02);
    meErrorY_->setAxisTitle("RecHit error Y", 1);

    // Removed to save offline memory
    // hid = theHistogramId->setHistoId("nRecHits",id_);
    // menRecHits_ = iBooker.book1D(hid, "# of rechits in this module", 8, 0,
    // 8); menRecHits_->setAxisTitle("number of rechits",1);
    delete theHistogramId;
  }

  if (type == 1 && barrel) {
    uint32_t DBladder;
    // if (!isUpgrade) { DBladder = PixelBarrelName(DetId(id_)).ladderName(); }
    // else { DBladder = PixelBarrelNameUpgrade(DetId(id_)).ladderName(); }
    DBladder = PixelBarrelName(DetId(id_), pTT, isUpgrade).ladderName();
    char sladder[80];
    sprintf(sladder, "Ladder_%02i", DBladder);
    hid = src.label() + "_" + sladder;
    if (isHalfModule)
      hid += "H";
    else
      hid += "F";
    if (!reducedSet) {
      if (twoD) {
        meXYPosLad_ = iBooker.book2D("xypos_" + hid, "XY Position", 100, -1., 1, 100, -4, 4);
        meXYPosLad_->setAxisTitle("X Position", 1);
        meXYPosLad_->setAxisTitle("Y Position", 2);
      } else {
        // projections of XYPosition
        meXYPosLad_px_ = iBooker.book1D("xypos_" + hid + "_px", "X Position", 100, -1., 1);
        meXYPosLad_px_->setAxisTitle("X Position", 1);
        meXYPosLad_py_ = iBooker.book1D("xypos_" + hid + "_py", "Y Position", 100, -4, 4);
        meXYPosLad_py_->setAxisTitle("Y Position", 1);
      }
    }
    meClustXLad_ = iBooker.book1D("ClustX_" + hid, "RecHit X size", 10, 0., 10.);
    meClustXLad_->setAxisTitle("RecHit size X dimension", 1);
    meClustYLad_ = iBooker.book1D("ClustY_" + hid, "RecHit Y size", 15, 0., 15.);
    meClustYLad_->setAxisTitle("RecHit size Y dimension", 1);
    meErrorXLad_ = iBooker.book1D("ErrorX_" + hid, "RecHit error X", 100, 0., 0.02);
    meErrorXLad_->setAxisTitle("RecHit error X", 1);
    meErrorYLad_ = iBooker.book1D("ErrorY_" + hid, "RecHit error Y", 100, 0., 0.02);
    meErrorYLad_->setAxisTitle("RecHit error Y", 1);
    menRecHitsLad_ = iBooker.book1D("nRecHits_" + hid, "# of rechits in this module", 8, 0, 8);
    menRecHitsLad_->setAxisTitle("number of rechits", 1);
  }

  if (type == 2 && barrel) {
    uint32_t DBlayer;
    // if (!isUpgrade) { DBlayer = PixelBarrelName(DetId(id_)).layerName(); }
    // else { DBlayer = PixelBarrelNameUpgrade(DetId(id_)).layerName(); }
    DBlayer = PixelBarrelName(DetId(id_), pTT, isUpgrade).layerName();
    char slayer[80];
    sprintf(slayer, "Layer_%i", DBlayer);
    hid = src.label() + "_" + slayer;

    if (!reducedSet) {
      if (twoD) {
        meXYPosLay_ = iBooker.book2D("xypos_" + hid, "XY Position", 100, -1., 1, 100, -4, 4);
        meXYPosLay_->setAxisTitle("X Position", 1);
        meXYPosLay_->setAxisTitle("Y Position", 2);
      } else {
        // projections of XYPosition
        meXYPosLay_px_ = iBooker.book1D("xypos_" + hid + "_px", "X Position", 100, -1., 1);
        meXYPosLay_px_->setAxisTitle("X Position", 1);
        meXYPosLay_py_ = iBooker.book1D("xypos_" + hid + "_py", "Y Position", 100, -4, 4);
        meXYPosLay_py_->setAxisTitle("Y Position", 1);
      }
    }

    meClustXLay_ = iBooker.book1D("ClustX_" + hid, "RecHit X size", 10, 0., 10.);
    meClustXLay_->setAxisTitle("RecHit size X dimension", 1);
    meClustYLay_ = iBooker.book1D("ClustY_" + hid, "RecHit Y size", 15, 0., 15.);
    meClustYLay_->setAxisTitle("RecHit size Y dimension", 1);
    meErrorXLay_ = iBooker.book1D("ErrorX_" + hid, "RecHit error X", 100, 0., 0.02);
    meErrorXLay_->setAxisTitle("RecHit error X", 1);
    meErrorYLay_ = iBooker.book1D("ErrorY_" + hid, "RecHit error Y", 100, 0., 0.02);
    meErrorYLay_->setAxisTitle("RecHit error Y", 1);
    menRecHitsLay_ = iBooker.book1D("nRecHits_" + hid, "# of rechits in this module", 8, 0, 8);
    menRecHitsLay_->setAxisTitle("number of rechits", 1);
  }

  if (type == 3 && barrel) {
    uint32_t DBmodule;
    // if (!isUpgrade) { DBmodule = PixelBarrelName(DetId(id_)).moduleName(); }
    // else { DBmodule = PixelBarrelNameUpgrade(DetId(id_)).moduleName(); }
    DBmodule = PixelBarrelName(DetId(id_), pTT, isUpgrade).moduleName();
    char smodule[80];
    sprintf(smodule, "Ring_%i", DBmodule);
    hid = src.label() + "_" + smodule;

    if (!reducedSet) {
      if (twoD) {
        meXYPosPhi_ = iBooker.book2D("xypos_" + hid, "XY Position", 100, -1., 1, 100, -4, 4);
        meXYPosPhi_->setAxisTitle("X Position", 1);
        meXYPosPhi_->setAxisTitle("Y Position", 2);
      } else {
        // projections of XYPosition
        meXYPosPhi_px_ = iBooker.book1D("xypos_" + hid + "_px", "X Position", 100, -1., 1);
        meXYPosPhi_px_->setAxisTitle("X Position", 1);
        meXYPosPhi_py_ = iBooker.book1D("xypos_" + hid + "_py", "Y Position", 100, -4, 4);
        meXYPosPhi_py_->setAxisTitle("Y Position", 1);
      }
    }
    meClustXPhi_ = iBooker.book1D("ClustX_" + hid, "RecHit X size", 10, 0., 10.);
    meClustXPhi_->setAxisTitle("RecHit size X dimension", 1);
    meClustYPhi_ = iBooker.book1D("ClustY_" + hid, "RecHit Y size", 15, 0., 15.);
    meClustYPhi_->setAxisTitle("RecHit size Y dimension", 1);
    meErrorXPhi_ = iBooker.book1D("ErrorX_" + hid, "RecHit error X", 100, 0., 0.02);
    meErrorXPhi_->setAxisTitle("RecHit error X", 1);
    meErrorYPhi_ = iBooker.book1D("ErrorY_" + hid, "RecHit error Y", 100, 0., 0.02);
    meErrorYPhi_->setAxisTitle("RecHit error Y", 1);
    menRecHitsPhi_ = iBooker.book1D("nRecHits_" + hid, "# of rechits in this module", 8, 0, 8);
    menRecHitsPhi_->setAxisTitle("number of rechits", 1);
  }

  if (type == 4 && endcap) {
    uint32_t blade;
    // if (!isUpgrade) { blade= PixelEndcapName(DetId(id_)).bladeName(); }
    // else { blade= PixelEndcapNameUpgrade(DetId(id_)).bladeName(); }
    blade = PixelEndcapName(DetId(id_), pTT, isUpgrade).bladeName();

    char sblade[80];
    sprintf(sblade, "Blade_%02i", blade);
    hid = src.label() + "_" + sblade;
    //     meXYPosBlade_ = iBooker.book2D("xypos_" + hid,"XY
    //     Position",100,-1.,1,100,-4,4); meXYPosBlade_->setAxisTitle("X
    //     Position",1); meXYPosBlade_->setAxisTitle("Y Position",2);

    meClustXBlade_ = iBooker.book1D("ClustX_" + hid, "RecHit X size", 10, 0., 10.);
    meClustXBlade_->setAxisTitle("RecHit size X dimension", 1);
    meClustYBlade_ = iBooker.book1D("ClustY_" + hid, "RecHit Y size", 15, 0., 15.);
    meClustYBlade_->setAxisTitle("RecHit size Y dimension", 1);
    meErrorXBlade_ = iBooker.book1D("ErrorX_" + hid, "RecHit error X", 100, 0., 0.02);
    meErrorXBlade_->setAxisTitle("RecHit error X", 1);
    meErrorYBlade_ = iBooker.book1D("ErrorY_" + hid, "RecHit error Y", 100, 0., 0.02);
    meErrorYBlade_->setAxisTitle("RecHit error Y", 1);
    menRecHitsBlade_ = iBooker.book1D("nRecHits_" + hid, "# of rechits in this module", 8, 0, 8);
    menRecHitsBlade_->setAxisTitle("number of rechits", 1);
  }
  if (type == 5 && endcap) {
    uint32_t disk;
    // if (!isUpgrade) { disk = PixelEndcapName(DetId(id_)).diskName(); }
    // else { disk = PixelEndcapNameUpgrade(DetId(id_)).diskName(); }
    disk = PixelEndcapName(DetId(id_), pTT, isUpgrade).diskName();

    char sdisk[80];
    sprintf(sdisk, "Disk_%i", disk);
    hid = src.label() + "_" + sdisk;
    //     meXYPosDisk_ = iBooker.book2D("xypos_" + hid,"XY
    //     Position",100,-1.,1,100,-4,4); meXYPosDisk_->setAxisTitle("X
    //     Position",1); meXYPosDisk_->setAxisTitle("Y Position",2);

    meClustXDisk_ = iBooker.book1D("ClustX_" + hid, "RecHit X size", 10, 0., 10.);
    meClustXDisk_->setAxisTitle("RecHit size X dimension", 1);
    meClustYDisk_ = iBooker.book1D("ClustY_" + hid, "RecHit Y size", 15, 0., 15.);
    meClustYDisk_->setAxisTitle("RecHit size Y dimension", 1);
    meErrorXDisk_ = iBooker.book1D("ErrorX_" + hid, "RecHit error X", 100, 0., 0.02);
    meErrorXDisk_->setAxisTitle("RecHit error X", 1);
    meErrorYDisk_ = iBooker.book1D("ErrorY_" + hid, "RecHit error Y", 100, 0., 0.02);
    meErrorYDisk_->setAxisTitle("RecHit error Y", 1);
    menRecHitsDisk_ = iBooker.book1D("nRecHits_" + hid, "# of rechits in this module", 8, 0, 8);
    menRecHitsDisk_->setAxisTitle("number of rechits", 1);
  }

  if (type == 6 && endcap) {
    uint32_t panel;
    uint32_t module;
    /*if (!isUpgrade) {
      panel= PixelEndcapName(DetId(id_)).pannelName();
      module= PixelEndcapName(DetId(id_)).plaquetteName();
    } else {
      panel= PixelEndcapNameUpgrade(DetId(id_)).pannelName();
      module= PixelEndcapNameUpgrade(DetId(id_)).plaquetteName();
    }*/
    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);
    hid = src.label() + "_" + slab;

    if (!reducedSet) {
      if (twoD) {
        meXYPosRing_ = iBooker.book2D("xypos_" + hid, "XY Position", 100, -1., 1, 100, -4, 4);
        meXYPosRing_->setAxisTitle("X Position", 1);
        meXYPosRing_->setAxisTitle("Y Position", 2);
      } else {
        // projections of XYPosition
        meXYPosRing_px_ = iBooker.book1D("xypos_" + hid + "_px", "X Position", 100, -1., 1);
        meXYPosRing_px_->setAxisTitle("X Position", 1);
        meXYPosRing_py_ = iBooker.book1D("xypos_" + hid + "_py", "Y Position", 100, -4, 4);
        meXYPosRing_py_->setAxisTitle("Y Position", 1);
      }
    }
    meClustXRing_ = iBooker.book1D("ClustX_" + hid, "RecHit X size", 10, 0., 10.);
    meClustXRing_->setAxisTitle("RecHit size X dimension", 1);
    meClustYRing_ = iBooker.book1D("ClustY_" + hid, "RecHit Y size", 15, 0., 15.);
    meClustYRing_->setAxisTitle("RecHit size Y dimension", 1);
    meErrorXRing_ = iBooker.book1D("ErrorX_" + hid, "RecHit error X", 100, 0., 0.02);
    meErrorXRing_->setAxisTitle("RecHit error X", 1);
    meErrorYRing_ = iBooker.book1D("ErrorY_" + hid, "RecHit error Y", 100, 0., 0.02);
    meErrorYRing_->setAxisTitle("RecHit error Y", 1);
    menRecHitsRing_ = iBooker.book1D("nRecHits_" + hid, "# of rechits in this module", 8, 0, 8);
    menRecHitsRing_->setAxisTitle("number of rechits", 1);
  }
}
//
// Fill histograms
//
void SiPixelRecHitModule::fill(const float &rechit_x,
                               const float &rechit_y,
                               const int &sizeX,
                               const int &sizeY,
                               const float &lerr_x,
                               const float &lerr_y,
                               bool modon,
                               bool ladon,
                               bool layon,
                               bool phion,
                               bool bladeon,
                               bool diskon,
                               bool ringon,
                               bool twoD,
                               bool reducedSet) {
  bool barrel = DetId(id_).subdetId() == static_cast<int>(PixelSubdetector::PixelBarrel);
  bool endcap = DetId(id_).subdetId() == static_cast<int>(PixelSubdetector::PixelEndcap);

  // std::cout << rechit_x << " " << rechit_y << " " << sizeX << " " << sizeY <<
  // std::endl;
  if (modon) {
    meClustX_->Fill(sizeX);
    meClustY_->Fill(sizeY);
    meErrorX_->Fill(lerr_x);
    meErrorY_->Fill(lerr_y);
    if (!reducedSet) {
      if (twoD)
        meXYPos_->Fill(rechit_x, rechit_y);
      else {
        meXYPos_px_->Fill(rechit_x);
        meXYPos_py_->Fill(rechit_y);
      }
    }
  }
  // std::cout<<"number of detector units="<<numberOfDetUnits<<std::endl;

  if (ladon && barrel) {
    meClustXLad_->Fill(sizeX);
    meClustYLad_->Fill(sizeY);
    meErrorXLad_->Fill(lerr_x);
    meErrorYLad_->Fill(lerr_y);
    if (!reducedSet) {
      if (twoD)
        meXYPosLad_->Fill(rechit_x, rechit_y);
      else {
        meXYPosLad_px_->Fill(rechit_x);
        meXYPosLad_py_->Fill(rechit_y);
      }
    }
  }

  if (layon && barrel) {
    meClustXLay_->Fill(sizeX);
    meClustYLay_->Fill(sizeY);
    meErrorXLay_->Fill(lerr_x);
    meErrorYLay_->Fill(lerr_y);
    if (!reducedSet) {
      if (twoD)
        meXYPosLay_->Fill(rechit_x, rechit_y);
      else {
        meXYPosLay_px_->Fill(rechit_x);
        meXYPosLay_py_->Fill(rechit_y);
      }
    }
  }

  if (phion && barrel) {
    meClustXPhi_->Fill(sizeX);
    meClustYPhi_->Fill(sizeY);
    meErrorXPhi_->Fill(lerr_x);
    meErrorYPhi_->Fill(lerr_y);
    if (!reducedSet) {
      if (twoD)
        meXYPosPhi_->Fill(rechit_x, rechit_y);
      else {
        meXYPosPhi_px_->Fill(rechit_x);
        meXYPosPhi_py_->Fill(rechit_y);
      }
    }
  }

  if (bladeon && endcap) {
    // meXYPosBlade_->Fill(rechit_x, rechit_y);
    meClustXBlade_->Fill(sizeX);
    meClustYBlade_->Fill(sizeY);
    meErrorXBlade_->Fill(lerr_x);
    meErrorYBlade_->Fill(lerr_y);
  }

  if (diskon && endcap) {
    // meXYPosDisk_->Fill(rechit_x, rechit_y);
    meClustXDisk_->Fill(sizeX);
    meClustYDisk_->Fill(sizeY);
    meErrorXDisk_->Fill(lerr_x);
    meErrorYDisk_->Fill(lerr_y);
  }

  if (ringon && endcap) {
    meClustXRing_->Fill(sizeX);
    meClustYRing_->Fill(sizeY);
    meErrorXRing_->Fill(lerr_x);
    meErrorYRing_->Fill(lerr_y);
    if (!reducedSet) {
      if (twoD)
        meXYPosRing_->Fill(rechit_x, rechit_y);
      else {
        meXYPosRing_px_->Fill(rechit_x);
        meXYPosRing_py_->Fill(rechit_y);
      }
    }
  }
}

void SiPixelRecHitModule::nfill(
    const int &nrec, 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);

  // if(modon) menRecHits_->Fill(nrec);
  // barrel
  if (ladon && barrel)
    menRecHitsLad_->Fill(nrec);
  if (layon && barrel)
    menRecHitsLay_->Fill(nrec);
  if (phion && barrel)
    menRecHitsPhi_->Fill(nrec);
  // endcap
  if (bladeon && endcap)
    menRecHitsBlade_->Fill(nrec);
  if (diskon && endcap)
    menRecHitsDisk_->Fill(nrec);
  if (ringon && endcap)
    menRecHitsRing_->Fill(nrec);
}