GEMRecHitValidation Class Reference

#include <GEMRecHitValidation.h>

Inheritance diagram for GEMRecHitValidation:

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override
void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
	GEMRecHitValidation (const edm::ParameterSet &)
	~GEMRecHitValidation () override
Public Member Functions inherited from GEMBaseValidation
	GEMBaseValidation (const edm::ParameterSet &, std::string)
	~GEMBaseValidation () override=0
Public Member Functions inherited from DQMEDAnalyzer
void	accumulate (edm::Event const &event, edm::EventSetup const &setup) final
void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final
void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final
void	beginStream (edm::StreamID id) final
virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)
	DQMEDAnalyzer ()
void	endLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final
void	endRun (edm::Run const &run, edm::EventSetup const &setup) final
virtual bool	getCanSaveByLumi ()
Public Member Functions inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
	EDProducer ()=default
	EDProducer (const EDProducer &)=delete
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () const final
bool	hasAbilityToProduceInEndRuns () const final
const EDProducer &	operator= (const EDProducer &)=delete

Private Member Functions
Bool_t	matchRecHitAgainstSimHit (GEMRecHitCollection::const_iterator, Int_t)

Private Attributes
edm::EDGetTokenT< edm::DetSetVector< GEMDigiSimLink > >	digisimlink_token_
edm::ESGetToken< GEMGeometry, MuonGeometryRecord >	geomToken_
edm::ESGetToken< GEMGeometry, MuonGeometryRecord >	geomTokenBeginRun_
MEMap4Ids	me_detail_cls_
MEMap2Ids	me_detail_cls_roll_
MonitorElement *	me_detail_cls_total_
MEMap3Ids	me_detail_occ_ieta_
MEMap3Ids	me_detail_occ_phi_
MEMap3Ids	me_detail_occ_pid_
MEMap3Ids	me_detail_occ_polar_
MEMap3Ids	me_detail_occ_xy_
MEMap1Ids	me_detail_occ_zr_
MEMap2Ids	me_detail_pull_x_
MEMap4Ids	me_detail_pull_x_la_
MEMap2Ids	me_detail_pull_y_
MEMap4Ids	me_detail_pull_y_la_
MEMap2Ids	me_detail_rechit_occ_det_
MEMap4Ids	me_detail_residual_rphi_
MEMap4Ids	me_detail_residual_y_
MEMap3Ids	me_detail_total_rechit_
MEMap3Ids	me_rechit_occ_eta_
MEMap3Ids	me_rechit_occ_phi_
MEMap2Ids	me_residual_rphi_
MEMap2Ids	me_residual_y_
edm::EDGetTokenT< GEMRecHitCollection >	rechit_token_
edm::EDGetTokenT< edm::PSimHitContainer >	simhit_token_

Additional Inherited Members
Public Types inherited from DQMEDAnalyzer
typedef dqm::reco::DQMStore	DQMStore
typedef dqm::reco::MonitorElement	MonitorElement
Public Types inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
using	CacheTypes = CacheContexts< T... >
using	GlobalCache = typename CacheTypes::GlobalCache
using	HasAbility = AbilityChecker< T... >
using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache
using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache
using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >
using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache
using	RunCache = typename CacheTypes::RunCache
using	RunContext = RunContextT< RunCache, GlobalCache >
using	RunSummaryCache = typename CacheTypes::RunSummaryCache
Static Public Member Functions inherited from DQMEDAnalyzer
static void	globalEndJob (DQMEDAnalyzerGlobalCache const *)
static void	globalEndLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup, LuminosityBlockContext const *context)
static void	globalEndRunProduce (edm::Run &run, edm::EventSetup const &setup, RunContext const *context)
static std::unique_ptr< DQMEDAnalyzerGlobalCache >	initializeGlobalCache (edm::ParameterSet const &)
Protected Member Functions inherited from GEMBaseValidation
template<typename T >
dqm::impl::MonitorElement *	bookDetectorOccupancy (DQMStore::IBooker &booker, const T &key, const GEMStation *station, const char *name_prfix, const char *title_prefix)
template<typename T >
dqm::impl::MonitorElement *	bookHist1D (DQMStore::IBooker &booker, const T &key, const char *name, const char *title, Int_t nbinsx, Double_t xlow, Double_t xup, const char *x_title="", const char *y_title="Entries")
template<typename T >
dqm::impl::MonitorElement *	bookHist2D (DQMStore::IBooker &booker, const T &key, const char *name, const char *title, Int_t nbinsx, Double_t xlow, Double_t xup, Int_t nbinsy, Double_t ylow, Double_t yup, const char *x_title="", const char *y_title="")
template<typename T >
dqm::impl::MonitorElement *	bookPIDHist (DQMStore::IBooker &booker, const T &key, const char *name, const char *title)
template<typename T >
dqm::impl::MonitorElement *	bookPIDHist (DQMStore::IBooker &booker, const T &key, Int_t ieta, const char *name, const char *title)
template<typename T >
dqm::impl::MonitorElement *	bookPolarOccupancy (DQMStore::IBooker &booker, const T &key, const char *name_prefix, const char *title_prefix)
template<typename T >
dqm::impl::MonitorElement *	bookXYOccupancy (DQMStore::IBooker &booker, const T &key, const char *name_prefix, const char *title_prefix)
dqm::impl::MonitorElement *	bookZROccupancy (DQMStore::IBooker &booker, Int_t region_id, const char *name_prfix, const char *title_prefix)
template<typename T >
dqm::impl::MonitorElement *	bookZROccupancy (DQMStore::IBooker &booker, const T &key, const char *name_prfix, const char *title_prefix)
Int_t	getDetOccBinX (Int_t num_layers, Int_t chamber_id, Int_t layer_id)
Int_t	getPidIdx (Int_t pid)
Bool_t	isMuonSimHit (const PSimHit &)
Float_t	toDegree (Float_t radian)
Protected Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const
Protected Attributes inherited from GEMBaseValidation
Bool_t	detail_plot_
std::vector< Double_t >	eta_range_
edm::ESGetToken< GEMGeometry, MuonGeometryRecord >	geomToken_
edm::ESGetToken< GEMGeometry, MuonGeometryRecord >	geomTokenBeginRun_
const std::string	kLogCategory_
const Int_t	kMuonPDGId_ = 13
std::vector< Int_t >	pid_list_
Int_t	xy_occ_num_bins_
std::vector< Int_t >	zr_occ_num_bins_
std::vector< Double_t >	zr_occ_range_
Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_
edm::EDPutTokenT< DQMToken >	runToken_
unsigned int	streamId_

Detailed Description

Definition at line 10 of file GEMRecHitValidation.h.

Constructor & Destructor Documentation

GEMRecHitValidation()

GEMRecHitValidation::GEMRecHitValidation ( const edm::ParameterSet &  pset )

explicit

Definition at line 6 of file GEMRecHitValidation.cc.

References digisimlink_token_, geomToken_, geomTokenBeginRun_, muonDTDigis_cfi::pset, rechit_token_, and simhit_token_.

    : GEMBaseValidation(pset, "GEMRecHitsValidation") {
  const auto& rechit_pset = pset.getParameterSet("gemRecHit");
  const auto& rechit_tag = rechit_pset.getParameter<edm::InputTag>("inputTag");
  rechit_token_ = consumes<GEMRecHitCollection>(rechit_tag);

  const auto& simhit_pset = pset.getParameterSet("gemSimHit");
  const auto& simhit_tag = simhit_pset.getParameter<edm::InputTag>("inputTag");

  const auto& digisimlink_tag = pset.getParameter<edm::InputTag>("gemDigiSimLink");
  digisimlink_token_ = consumes<edm::DetSetVector<GEMDigiSimLink>>(digisimlink_tag);

  simhit_token_ = consumes<edm::PSimHitContainer>(simhit_tag);
  geomToken_ = esConsumes<GEMGeometry, MuonGeometryRecord>();
  geomTokenBeginRun_ = esConsumes<GEMGeometry, MuonGeometryRecord, edm::Transition::BeginRun>();
}

~GEMRecHitValidation()

GEMRecHitValidation::~GEMRecHitValidation ( )

override

Definition at line 23 of file GEMRecHitValidation.cc.

{}

Member Function Documentation

analyze()

void GEMRecHitValidation::analyze ( const edm::Event &  event, const edm::EventSetup &  setup  )

overridevirtual

Implements GEMBaseValidation.

Definition at line 245 of file GEMRecHitValidation.cc.

References submitPVResolutionJobs::count, GEMBaseValidation::detail_plot_, digisimlink_token_, PV3DBase< T, PVType, FrameType >::eta(), dqm::impl::MonitorElement::Fill(), mixOne_premix_on_sim_cfi::gem, geomToken_, GEMBaseValidation::getDetOccBinX(), GEMBaseValidation::getPidIdx(), hcalRecHitTable_cff::ieta, GEMBaseValidation::isMuonSimHit(), edm::HandleBase::isValid(), crabWrapper::key, GEMBaseValidation::kLogCategory_, electronStore::links, matchRecHitAgainstSimHit(), me_detail_cls_, me_detail_cls_roll_, me_detail_cls_total_, me_detail_occ_ieta_, me_detail_occ_phi_, me_detail_occ_pid_, me_detail_occ_polar_, me_detail_occ_xy_, me_detail_occ_zr_, me_detail_pull_x_, me_detail_pull_x_la_, me_detail_pull_y_, me_detail_pull_y_la_, me_detail_rechit_occ_det_, me_detail_residual_rphi_, me_detail_residual_y_, me_detail_total_rechit_, me_rechit_occ_eta_, me_rechit_occ_phi_, me_residual_rphi_, me_residual_y_, PV3DBase< T, PVType, FrameType >::perp(), PV3DBase< T, PVType, FrameType >::phi(), conifer::pow(), edm::Handle< T >::product(), FastTimerService_cff::range, rechit_token_, GEMDetId::region(), singleTopDQM_cfi::setup, simhit_token_, mathSSE::sqrt(), nano_mu_digi_cff::strip, GEMBaseValidation::toDegree(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

                                                                                   {
  const GEMGeometry* gem = &setup.getData(geomToken_);

  edm::Handle<edm::DetSetVector<GEMDigiSimLink>> digiSimLink;
  event.getByToken(digisimlink_token_, digiSimLink);
  if (not digiSimLink.isValid()) {
    edm::LogError(kLogCategory_) << "Failed to get GEMDigiSimLink." << std::endl;
    return;
  }

  edm::Handle<edm::PSimHitContainer> simhit_container;
  event.getByToken(simhit_token_, simhit_container);
  if (not simhit_container.isValid()) {
    edm::LogError(kLogCategory_) << "Failed to get PSimHitContainer." << std::endl;
    return;
  }

  edm::Handle<GEMRecHitCollection> rechit_collection;
  event.getByToken(rechit_token_, rechit_collection);
  if (not rechit_collection.isValid()) {
    edm::LogError(kLogCategory_) << "Failed to get GEMRecHitCollection" << std::endl;
    return;
  }

  std::map<ME3IdsKey, Int_t> total_rechit;
  for (const auto& rechit : *rechit_collection) {
    GEMDetId gem_id{rechit.gemId()};
    Int_t region_id = gem_id.region();
    Int_t station_id = gem_id.station();
    Int_t layer_id = gem_id.layer();
    Int_t ieta = gem_id.ieta();
    ME4IdsKey key4{region_id, station_id, layer_id, ieta};

    ME2IdsKey key{station_id, ieta};
    ME2IdsKey key2{region_id, station_id};
    ME3IdsKey key3{region_id, station_id, layer_id};

    const BoundPlane& surface = gem->idToDet(gem_id)->surface();
    GlobalPoint&& rechit_global_pos = surface.toGlobal(rechit.localPosition());

    Float_t rechit_g_x = rechit_global_pos.x();
    Float_t rechit_g_y = rechit_global_pos.y();
    Float_t rechit_g_abs_z = std::fabs(rechit_global_pos.z());
    Float_t rechit_g_r = rechit_global_pos.perp();
    Float_t rechit_g_phi = toDegree(rechit_global_pos.phi());

    Int_t first_strip = rechit.firstClusterStrip();
    Int_t cls = rechit.clusterSize();
    cls = cls > 10 ? 10 : cls;

    total_rechit[key3]++;

    if (detail_plot_) {
      me_detail_cls_total_->Fill(cls);
      me_detail_cls_[key4]->Fill(cls);
      me_detail_cls_roll_[key]->Fill(cls);
      me_detail_occ_ieta_[key3]->Fill(ieta);
      me_detail_occ_phi_[key3]->Fill(rechit_g_phi);
      me_detail_occ_zr_[region_id]->Fill(rechit_g_abs_z, rechit_g_r);
      me_detail_occ_xy_[key3]->Fill(rechit_g_x, rechit_g_y);
      me_detail_occ_polar_[key3]->Fill(rechit_g_phi, rechit_g_r);
    }  // detail plot

    auto links = digiSimLink->find(gem_id);

    if (links == digiSimLink->end())
      continue;
    std::map<Int_t, Int_t> pid_count;

    for (Int_t strip = first_strip; strip < first_strip + cls; strip++) {
      for (const auto& link : *links) {
        Int_t link_strip = link.getStrip();
        if (link_strip == strip) {
          Int_t pid = link.getParticleType();
          pid_count[pid]++;
          break;
        }
      }
    }
    if (detail_plot_) {
      Int_t max_pid = 0;
      Int_t max_count = 0;
      for (auto& [pid, count] : pid_count) {
        if (max_count < count) {
          max_pid = pid;
          max_count = count;
        }
      }
      Int_t pid_idx = getPidIdx(max_pid);
      me_detail_occ_pid_[key3]->Fill(pid_idx);
    }
  }

  if (detail_plot_) {
    for (auto [key, num_total_rechit] : total_rechit) {
      me_detail_total_rechit_[key]->Fill(num_total_rechit);
    }
  }

  // NOTE
  for (const auto& simhit : *simhit_container.product()) {
    if (gem->idToDet(simhit.detUnitId()) == nullptr) {
      edm::LogError(kLogCategory_) << "MuonGEMHit didn't matched with GEMGeometry." << std::endl;
      continue;
    }

    GEMDetId simhit_gemid{simhit.detUnitId()};
    const BoundPlane& surface = gem->idToDet(simhit_gemid)->surface();

    Int_t region_id = simhit_gemid.region();
    Int_t station_id = simhit_gemid.station();
    Int_t layer_id = simhit_gemid.layer();
    Int_t chamber_id = simhit_gemid.chamber();
    Int_t ieta = simhit_gemid.ieta();
    Int_t num_layers = simhit_gemid.nlayers();

    ME2IdsKey key{station_id, ieta};
    ME2IdsKey key2{region_id, station_id};
    ME3IdsKey key3{region_id, station_id, layer_id};
    ME4IdsKey key4{region_id, station_id, layer_id, ieta};

    const LocalPoint& simhit_local_pos = simhit.localPosition();
    const GlobalPoint& simhit_global_pos = surface.toGlobal(simhit_local_pos);

    Float_t simhit_g_abs_eta = std::fabs(simhit_global_pos.eta());
    Float_t simhit_g_phi = toDegree(simhit_global_pos.phi());

    Int_t det_occ_bin_x = getDetOccBinX(num_layers, chamber_id, layer_id);

    auto simhit_trackId = simhit.trackId();

    auto links = digiSimLink->find(simhit_gemid);
    if (links == digiSimLink->end())
      continue;

    Int_t simhit_strip = -1;
    for (const auto& link : *links) {
      if (simhit_trackId == link.getTrackId()) {
        simhit_strip = link.getStrip();
        break;
      }
    }

    GEMRecHitCollection::range range = rechit_collection->get(simhit_gemid);
    for (auto rechit = range.first; rechit != range.second; ++rechit) {
      if (gem->idToDet(rechit->gemId()) == nullptr) {
        edm::LogError(kLogCategory_) << "GEMRecHit didn't matched with GEMGeometry." << std::endl;
        continue;
      }

      if (not isMuonSimHit(simhit))
        continue;

      if (matchRecHitAgainstSimHit(rechit, simhit_strip)) {
        const LocalPoint& rechit_local_pos = rechit->localPosition();

        Float_t resolution_x = std::sqrt(rechit->localPositionError().xx());
        Float_t resolution_y = std::sqrt(rechit->localPositionError().yy());

        Float_t residual_x = rechit_local_pos.x() - simhit_local_pos.x();
        Float_t residual_y = rechit_local_pos.y() - simhit_local_pos.y();
        Float_t residual_r = sqrt(pow(residual_x, 2) + pow(residual_y, 2));
        Float_t residual_phi = rechit_local_pos.phi() - simhit_local_pos.phi();
        Float_t residual_rphi = residual_r * residual_phi;

        Float_t pull_x = residual_x / resolution_x;
        Float_t pull_y = residual_y / resolution_y;

        me_residual_y_[key]->Fill(residual_y);
        me_residual_rphi_[key]->Fill(residual_rphi);

        me_rechit_occ_eta_[key3]->Fill(simhit_g_abs_eta);
        me_rechit_occ_phi_[key3]->Fill(simhit_g_phi);

        if (detail_plot_) {
          me_detail_rechit_occ_det_[key2]->Fill(det_occ_bin_x, ieta);

          me_detail_residual_y_[key4]->Fill(residual_y);
          me_detail_residual_rphi_[key4]->Fill(residual_rphi);

          me_detail_pull_x_[key]->Fill(pull_x);
          me_detail_pull_y_[key]->Fill(pull_y);
          me_detail_pull_x_la_[key4]->Fill(pull_x);
          me_detail_pull_y_la_[key4]->Fill(pull_y);
        }  // detail_plot

        // If we find GEMRecHit that matches PSimHit, then exit
        // GEMRecHitCollection loop.
        break;

      }  // if rechit matches against simhit
    }    // rechit loop
  }      // simhit loop
}

bookHistograms()

void GEMRecHitValidation::bookHistograms ( DQMStore::IBooker &  booker, edm::Run const &  run, edm::EventSetup const &  setup  )

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 25 of file GEMRecHitValidation.cc.

References dqm::implementation::IBooker::book1D(), GEMBaseValidation::bookDetectorOccupancy(), GEMBaseValidation::bookHist1D(), GEMBaseValidation::bookPIDHist(), GEMBaseValidation::bookPolarOccupancy(), GEMBaseValidation::bookXYOccupancy(), GEMBaseValidation::bookZROccupancy(), relativeConstraints::chamber, chambers, GEMBaseValidation::detail_plot_, GEMBaseValidation::eta_range_, mixOne_premix_on_sim_cfi::gem, geomTokenBeginRun_, hcalRecHitTable_cff::ieta, crabWrapper::key, GEMBaseValidation::kLogCategory_, me_detail_cls_, me_detail_cls_roll_, me_detail_cls_total_, me_detail_occ_ieta_, me_detail_occ_phi_, me_detail_occ_pid_, me_detail_occ_polar_, me_detail_occ_xy_, me_detail_occ_zr_, me_detail_pull_x_, me_detail_pull_x_la_, me_detail_pull_y_, me_detail_pull_y_la_, me_detail_rechit_occ_det_, me_detail_residual_rphi_, me_detail_residual_y_, me_detail_total_rechit_, me_rechit_occ_eta_, me_rechit_occ_phi_, me_residual_rphi_, me_residual_y_, nano_mu_digi_cff::region, nano_mu_digi_cff::roll, dqm::implementation::NavigatorBase::setCurrentFolder(), singleTopDQM_cfi::setup, and relativeConstraints::station.

                                                                                                               {
  const GEMGeometry* gem = &setup.getData(geomTokenBeginRun_);

  // NOTE Cluster Size
  booker.setCurrentFolder("GEM/RecHits");

  TString cls_title = "Cluster Size Distribution";
  TString cls_x_title = "Cluster size";

  if (detail_plot_) {
    me_detail_cls_total_ = booker.book1D("cls", cls_title + ";" + cls_x_title + ";" + "Entries", 10, 0.5, 10.5);

    for (const auto& station : gem->regions()[0]->stations()) {
      Int_t station_id = station->station();
      for (const auto& roll : station->superChambers()[0]->chambers()[0]->etaPartitions()) {
        Int_t ieta = roll->id().ieta();
        ME2IdsKey key{station_id, ieta};
        me_detail_cls_roll_[key] = booker.book1D(Form("cls_GE%d1-E%d", station_id, ieta),
                                                 Form("Cluster Size Distribution : GE%d1-E%d", station_id, ieta),
                                                 10,
                                                 0.5,
                                                 10.5);
      }
    }

    for (const auto& region : gem->regions()) {
      Int_t region_id = region->region();

      for (const auto& station : region->stations()) {
        Int_t station_id = station->station();

        const auto& superChamberVec = station->superChambers();
        if (superChamberVec.empty() || superChamberVec[0] == nullptr) {
          edm::LogError(kLogCategory_) << "Super chambers missing or null for region = " << region_id
                                       << " and station = " << station_id;
        } else {
          for (const auto& chamber : superChamberVec[0]->chambers()) {
            Int_t layer_id = chamber->id().layer();

            for (const auto& roll : chamber->etaPartitions()) {
              Int_t ieta = roll->id().ieta();
              ME4IdsKey key4{region_id, station_id, layer_id, ieta};

              me_detail_cls_[key4] = bookHist1D(booker, key4, "cls", "Cluster Size Distribution", 11, -0.5, 10.5);
            }  // roll loop
          }    // chamber loop
        }      // end else
      }        // station loop
    }          // region loop
  }            // detail plot

  // NOTE Residual
  for (const auto& station : gem->regions()[0]->stations()) {
    Int_t station_id = station->station();
    for (const auto& roll : station->superChambers()[0]->chambers()[0]->etaPartitions()) {
      Int_t ieta = roll->id().ieta();
      ME2IdsKey key{station_id, ieta};

      me_residual_y_[key] = booker.book1D(Form("residual_y_GE%d1-E%d", station_id, ieta),
                                          Form("Residual in Y : GE%d1-E%d; Residual in Y [cm]", station_id, ieta),
                                          60,
                                          -15,
                                          15);

      me_residual_rphi_[key] =
          booker.book1D(Form("residual_rphi_GE%d1-E%d", station_id, ieta),
                        Form("Residual in R #times #phi : GE%d1-E%d; Residual in r #times #phi [cm]", station_id, ieta),
                        60,
                        -15,
                        15);
    }
  }

  if (detail_plot_) {
    for (const auto& region : gem->regions()) {
      Int_t region_id = region->region();

      for (const auto& station : region->stations()) {
        Int_t station_id = station->station();

        const auto& superChamberVec = station->superChambers();
        if (!superChamberVec.empty() && superChamberVec[0] != nullptr) {
          for (const auto& chamber : superChamberVec[0]->chambers()) {
            Int_t layer_id = chamber->id().layer();

            for (const auto& roll : chamber->etaPartitions()) {
              Int_t ieta = roll->id().ieta();
              ME4IdsKey key4{region_id, station_id, layer_id, ieta};

              me_detail_residual_y_[key4] =
                  bookHist1D(booker, key4, "residual_y", "Residual in y", 60, -15, 15, "Residual in y [cm]");

              me_detail_residual_rphi_[key4] = bookHist1D(booker,
                                                          key4,
                                                          "residual_rphi",
                                                          "Residual in r #times #phi",
                                                          60,
                                                          -15,
                                                          15,
                                                          "Residual in r #times #phi [cm]");
            }  // roll loop
          }    // chamber loop
        }      // end if
      }        // station loop
    }          // region loop
  }            // detail plot

  // NOTE Pull
  if (detail_plot_) {
    for (const auto& station : gem->regions()[0]->stations()) {
      Int_t station_id = station->station();
      for (const auto& roll : station->superChambers()[0]->chambers()[0]->etaPartitions()) {
        Int_t ieta = roll->id().ieta();
        ME2IdsKey key{station_id, ieta};

        me_detail_pull_x_[key] = booker.book1D(
            Form("pull_x_GE%d1-E%d", station_id, ieta), Form("Pull in X : GE%d1-E%d", station_id, ieta), 60, -3, 3);

        me_detail_pull_y_[key] = booker.book1D(
            Form("pull_y_GE%d1-E%d", station_id, ieta), Form("Pull in Y : GE%d1-E%d", station_id, ieta), 60, -3, 3);
      }
    }

    for (const auto& region : gem->regions()) {
      Int_t region_id = region->region();

      for (const auto& station : region->stations()) {
        Int_t station_id = station->station();

        const auto& superChamberVec = station->superChambers();
        if (!superChamberVec.empty() && superChamberVec[0] != nullptr) {
          for (const auto& chamber : superChamberVec[0]->chambers()) {
            Int_t layer_id = chamber->id().layer();

            for (const auto& roll : chamber->etaPartitions()) {
              Int_t ieta = roll->id().ieta();
              ME4IdsKey key4{region_id, station_id, layer_id, ieta};

              me_detail_pull_x_la_[key4] = bookHist1D(booker, key4, "pull_x", "Pull in x", 60, -3, 3);

              me_detail_pull_y_la_[key4] = bookHist1D(booker, key4, "pull_y", "Pull in y", 60, -3, 3);
            }  // roll loop
          }    // chamber loop
        }      // end if
      }        // station loop
    }          // region loop
  }            // detail plot

  // NOTE Occupancy
  for (const auto& region : gem->regions()) {
    Int_t region_id = region->region();

    if (detail_plot_)
      me_detail_occ_zr_[region_id] = bookZROccupancy(booker, region_id, "rechit", "RecHit");

    for (const auto& station : region->stations()) {
      Int_t station_id = station->station();
      ME2IdsKey key2{region_id, station_id};

      if (detail_plot_)
        me_detail_rechit_occ_det_[key2] = bookDetectorOccupancy(booker, key2, station, "sim_matched", "Matched RecHit");

      const auto& superChamberVec = station->superChambers();
      if (!superChamberVec.empty() && superChamberVec[0] != nullptr) {
        for (const auto& chamber : superChamberVec[0]->chambers()) {
          Int_t layer_id = chamber->id().layer();
          ME3IdsKey key3{region_id, station_id, layer_id};

          Int_t num_eta_partitions = chamber->nEtaPartitions();

          me_rechit_occ_eta_[key3] = bookHist1D(booker,
                                                key3,
                                                "sim_matched_occ_eta",
                                                "Matched RecHit Eta Occupancy",
                                                16,
                                                eta_range_[station_id * 2 + 0],
                                                eta_range_[station_id * 2 + 1],
                                                "|#eta|");

          me_rechit_occ_phi_[key3] =
              bookHist1D(booker, key3, "sim_matched_occ_phi", "Matched RecHit Phi Occupancy", 36, -5, 355, "#phi");

          if (detail_plot_) {
            me_detail_total_rechit_[key3] =
                bookHist1D(booker, key3, "total_rechit", "Number of rec hits per event", 25, -0.5, 24.5);

            me_detail_occ_pid_[key3] = bookPIDHist(booker, key3, "sim_occ_pid", "Number of entreis for each particle");

            me_detail_occ_ieta_[key3] = bookHist1D(booker,
                                                   key3,
                                                   "occ_ieta",
                                                   "Rechit Occupancy per eta partition",
                                                   num_eta_partitions,
                                                   0.5,
                                                   num_eta_partitions + 0.5);

            me_detail_occ_phi_[key3] = bookHist1D(booker, key3, "occ_phi", "Rechit Phi Occupancy", 108, -5, 355);

            me_detail_occ_xy_[key3] = bookXYOccupancy(booker, key3, "rechit", "RecHit");

            me_detail_occ_polar_[key3] = bookPolarOccupancy(booker, key3, "rechit", "RecHit");
          }
        }  // chamber loop
      }    // end if
    }      // station loop
  }        // region_loop
}

matchRecHitAgainstSimHit()

Bool_t GEMRecHitValidation::matchRecHitAgainstSimHit ( GEMRecHitCollection::const_iterator  rechit, Int_t  simhit_strip  )

private

Definition at line 233 of file GEMRecHitValidation.cc.

Referenced by analyze().

                                                                                                                 {
  Int_t cls = rechit->clusterSize();
  Int_t rechit_first_strip = rechit->firstClusterStrip();

  if (cls == 1) {
    return simhit_strip == rechit_first_strip;
  } else {
    Int_t rechit_last_strip = rechit_first_strip + cls - 1;
    return (simhit_strip >= rechit_first_strip) and (simhit_strip <= rechit_last_strip);
  }
}

Member Data Documentation

digisimlink_token_

edm::EDGetTokenT<edm::DetSetVector<GEMDigiSimLink> > GEMRecHitValidation::digisimlink_token_

private

Definition at line 23 of file GEMRecHitValidation.h.

Referenced by analyze(), and GEMRecHitValidation().

geomToken_

edm::ESGetToken<GEMGeometry, MuonGeometryRecord> GEMRecHitValidation::geomToken_

private

Definition at line 24 of file GEMRecHitValidation.h.

Referenced by analyze(), and GEMRecHitValidation().

geomTokenBeginRun_

edm::ESGetToken<GEMGeometry, MuonGeometryRecord> GEMRecHitValidation::geomTokenBeginRun_

private

Definition at line 25 of file GEMRecHitValidation.h.

Referenced by bookHistograms(), and GEMRecHitValidation().

me_detail_cls_

MEMap4Ids GEMRecHitValidation::me_detail_cls_

private

Definition at line 30 of file GEMRecHitValidation.h.

Referenced by analyze(), and bookHistograms().

me_detail_cls_roll_

MEMap2Ids GEMRecHitValidation::me_detail_cls_roll_

private

Definition at line 29 of file GEMRecHitValidation.h.

Referenced by analyze(), and bookHistograms().

me_detail_cls_total_

MonitorElement* GEMRecHitValidation::me_detail_cls_total_

private

Definition at line 28 of file GEMRecHitValidation.h.

Referenced by analyze(), and bookHistograms().

me_detail_occ_ieta_

MEMap3Ids GEMRecHitValidation::me_detail_occ_ieta_

private

Definition at line 46 of file GEMRecHitValidation.h.

Referenced by analyze(), and bookHistograms().

me_detail_occ_phi_

MEMap3Ids GEMRecHitValidation::me_detail_occ_phi_

private

Definition at line 47 of file GEMRecHitValidation.h.

Referenced by analyze(), and bookHistograms().

me_detail_occ_pid_

MEMap3Ids GEMRecHitValidation::me_detail_occ_pid_

private

Definition at line 45 of file GEMRecHitValidation.h.

Referenced by analyze(), and bookHistograms().

me_detail_occ_polar_

MEMap3Ids GEMRecHitValidation::me_detail_occ_polar_

private

Definition at line 50 of file GEMRecHitValidation.h.

Referenced by analyze(), and bookHistograms().

me_detail_occ_xy_

MEMap3Ids GEMRecHitValidation::me_detail_occ_xy_

private

Definition at line 49 of file GEMRecHitValidation.h.

Referenced by analyze(), and bookHistograms().

me_detail_occ_zr_

MEMap1Ids GEMRecHitValidation::me_detail_occ_zr_

private

Definition at line 48 of file GEMRecHitValidation.h.

Referenced by analyze(), and bookHistograms().

me_detail_pull_x_

MEMap2Ids GEMRecHitValidation::me_detail_pull_x_

private

Definition at line 38 of file GEMRecHitValidation.h.

Referenced by analyze(), and bookHistograms().

me_detail_pull_x_la_

MEMap4Ids GEMRecHitValidation::me_detail_pull_x_la_

private

Definition at line 40 of file GEMRecHitValidation.h.

Referenced by analyze(), and bookHistograms().

me_detail_pull_y_

MEMap2Ids GEMRecHitValidation::me_detail_pull_y_

private

Definition at line 39 of file GEMRecHitValidation.h.

Referenced by analyze(), and bookHistograms().

me_detail_pull_y_la_

MEMap4Ids GEMRecHitValidation::me_detail_pull_y_la_

private

Definition at line 41 of file GEMRecHitValidation.h.

Referenced by analyze(), and bookHistograms().

me_detail_rechit_occ_det_

MEMap2Ids GEMRecHitValidation::me_detail_rechit_occ_det_

private

Definition at line 55 of file GEMRecHitValidation.h.

Referenced by analyze(), and bookHistograms().

me_detail_residual_rphi_

MEMap4Ids GEMRecHitValidation::me_detail_residual_rphi_

private

Definition at line 36 of file GEMRecHitValidation.h.

Referenced by analyze(), and bookHistograms().

me_detail_residual_y_

MEMap4Ids GEMRecHitValidation::me_detail_residual_y_

private

Definition at line 35 of file GEMRecHitValidation.h.

Referenced by analyze(), and bookHistograms().

me_detail_total_rechit_

MEMap3Ids GEMRecHitValidation::me_detail_total_rechit_

private

Definition at line 44 of file GEMRecHitValidation.h.

Referenced by analyze(), and bookHistograms().

me_rechit_occ_eta_

MEMap3Ids GEMRecHitValidation::me_rechit_occ_eta_

private

Definition at line 53 of file GEMRecHitValidation.h.

Referenced by analyze(), and bookHistograms().

me_rechit_occ_phi_

MEMap3Ids GEMRecHitValidation::me_rechit_occ_phi_

private

Definition at line 54 of file GEMRecHitValidation.h.

Referenced by analyze(), and bookHistograms().

me_residual_rphi_

MEMap2Ids GEMRecHitValidation::me_residual_rphi_

private

Definition at line 34 of file GEMRecHitValidation.h.

Referenced by analyze(), and bookHistograms().

me_residual_y_

MEMap2Ids GEMRecHitValidation::me_residual_y_

private

Definition at line 33 of file GEMRecHitValidation.h.

Referenced by analyze(), and bookHistograms().

rechit_token_

edm::EDGetTokenT<GEMRecHitCollection> GEMRecHitValidation::rechit_token_

private

Definition at line 21 of file GEMRecHitValidation.h.

Referenced by analyze(), and GEMRecHitValidation().

simhit_token_

edm::EDGetTokenT<edm::PSimHitContainer> GEMRecHitValidation::simhit_token_

private

Definition at line 22 of file GEMRecHitValidation.h.

Referenced by analyze(), and GEMRecHitValidation().