#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

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

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

Private Attributes
edm::ESGetToken< GEMGeometry, MuonGeometryRecord >	geomToken_

edm::ESGetToken< GEMGeometry, MuonGeometryRecord >	geomTokenBeginRun_

MonitorElement *	me_cls_

MEMap3Ids	me_detail_cls_

MEMap3Ids	me_detail_occ_polar_

MEMap3Ids	me_detail_occ_xy_

MEMap3Ids	me_detail_pull_x_

MEMap3Ids	me_detail_pull_y_

MEMap3Ids	me_detail_residual_x_

MEMap3Ids	me_detail_residual_y_

MEMap1Ids	me_occ_zr_

MEMap1Ids	me_pull_x_

MEMap1Ids	me_pull_y_

MEMap2Ids	me_rechit_occ_det_

MEMap1Ids	me_rechit_occ_eta_

MEMap2Ids	me_rechit_occ_phi_

MEMap1Ids	me_residual_x_

MEMap1Ids	me_residual_y_

MEMap2Ids	me_simhit_occ_det_

MEMap1Ids	me_simhit_occ_eta_

MEMap2Ids	me_simhit_occ_phi_

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 >
typedef CacheContexts< T... >	CacheTypes

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T... >	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	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 *	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)

template<typename T >
dqm::impl::MonitorElement *	bookZROccupancy (DQMStore::IBooker &booker, const T &key, const char name_prfix, const char title_prefix)

dqm::impl::MonitorElement *	bookZROccupancy (DQMStore::IBooker &booker, Int_t region_id, const char name_prfix, const char title_prefix)

Int_t	getDetOccBinX (Int_t chamber_id, Int_t layer_id)

Bool_t	isMuonSimHit (const PSimHit &)

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

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 7 of file GEMRecHitValidation.h.

Constructor & Destructor Documentation

◆ GEMRecHitValidation()

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

explicit

Definition at line 6 of file GEMRecHitValidation.cc.

     : 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");
   simhit_token_ = consumes<edm::PSimHitContainer>(simhit_tag);
   geomToken_ = esConsumes<GEMGeometry, MuonGeometryRecord>();
   geomTokenBeginRun_ = esConsumes<GEMGeometry, MuonGeometryRecord, edm::Transition::BeginRun>();
 }

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

◆ ~GEMRecHitValidation()

GEMRecHitValidation::~GEMRecHitValidation ( )

override

Definition at line 19 of file GEMRecHitValidation.cc.

19 {}

Member Function Documentation

◆ analyze()

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

overridevirtual

Implements GEMBaseValidation.

Definition at line 193 of file GEMRecHitValidation.cc.

                                                                                    {
   const GEMGeometry* gem = &setup.getData(geomToken_);
  
   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;
   }
  
   for (auto rechit = rechit_collection->begin(); rechit != rechit_collection->end(); rechit++) {
     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();
  
     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 = rechit_global_pos.phi();
  
     Int_t cls = rechit->clusterSize();
  
     me_cls_->Fill(cls);
     me_occ_zr_[region_id]->Fill(rechit_g_abs_z, rechit_g_r);
  
     if (detail_plot_) {
       me_detail_cls_[key3]->Fill(cls);
       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
   }
  
   // 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;
     }
  
     if (not isMuonSimHit(simhit))
       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 roll_id = simhit_gemid.roll();
  
     ME2IdsKey key2{region_id, station_id};
     ME3IdsKey key3{region_id, station_id, layer_id};
  
     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 = simhit_global_pos.phi();
  
     Int_t simhit_strip = gem->etaPartition(simhit_gemid)->strip(simhit_local_pos);
     Int_t det_occ_bin_x = getDetOccBinX(chamber_id, layer_id);
  
     me_simhit_occ_eta_[region_id]->Fill(simhit_g_abs_eta);
     me_simhit_occ_phi_[key2]->Fill(simhit_g_phi);
     me_simhit_occ_det_[key2]->Fill(det_occ_bin_x, roll_id);
  
     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 (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 pull_x = residual_x / resolution_x;
         Float_t pull_y = residual_y / resolution_y;
  
         me_residual_x_[region_id]->Fill(residual_x);
         me_residual_y_[region_id]->Fill(residual_y);
         me_pull_x_[region_id]->Fill(pull_x);
         me_pull_y_[region_id]->Fill(pull_y);
  
         me_rechit_occ_eta_[region_id]->Fill(simhit_g_abs_eta);
         me_rechit_occ_phi_[key2]->Fill(simhit_g_phi);
         me_rechit_occ_det_[key2]->Fill(det_occ_bin_x, roll_id);
  
         if (detail_plot_) {
           me_detail_residual_x_[key3]->Fill(residual_x);
           me_detail_residual_y_[key3]->Fill(residual_y);
  
           me_detail_pull_x_[key3]->Fill(pull_x);
           me_detail_pull_y_[key3]->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
 }

References GEMBaseValidation::detail_plot_, PV3DBase< T, PVType, FrameType >::eta(), dqm::impl::MonitorElement::Fill(), geomToken_, GEMBaseValidation::getDetOccBinX(), GEMBaseValidation::isMuonSimHit(), edm::HandleBase::isValid(), GEMBaseValidation::kLogCategory_, matchRecHitAgainstSimHit(), me_cls_, me_detail_cls_, me_detail_occ_polar_, me_detail_occ_xy_, me_detail_pull_x_, me_detail_pull_y_, me_detail_residual_x_, me_detail_residual_y_, me_occ_zr_, me_pull_x_, me_pull_y_, me_rechit_occ_det_, me_rechit_occ_eta_, me_rechit_occ_phi_, me_residual_x_, me_residual_y_, me_simhit_occ_det_, me_simhit_occ_eta_, me_simhit_occ_phi_, PV3DBase< T, PVType, FrameType >::perp(), PV3DBase< T, PVType, FrameType >::phi(), edm::Handle< T >::product(), FastTimerService_cff::range, rechit_token_, GEMDetId::region(), singleTopDQM_cfi::setup, simhit_token_, mathSSE::sqrt(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

◆ bookHistograms()

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

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 21 of file GEMRecHitValidation.cc.

                                                                                                                {
   const GEMGeometry* gem = &setup.getData(geomTokenBeginRun_);
  
   // NOTE Cluster Size
   booker.setCurrentFolder("MuonGEMRecHitsV/GEMRecHitsTask/ClusterSize");
  
   TString cls_title = "Cluster Size Distribution";
   TString cls_x_title = "Cluster size";
  
   me_cls_ = booker.book1D("cls", cls_title + ";" + cls_x_title + ";" + "Entries", 11, -0.5, 10.5);
  
   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) {
           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();
             ME3IdsKey key3{region_id, station_id, layer_id};
  
             me_detail_cls_[key3] = bookHist1D(booker, key3, "cls", cls_title, 11, -0.5, 10.5, cls_x_title);
           }  // chamber loop
         }    // end else
       }      // station loop
     }        // region loop
   }          // detail plot
  
   // NOTE Residual
   booker.setCurrentFolder("MuonGEMRecHitsV/GEMRecHitsTask/Residual");
  
   for (const auto& region : gem->regions()) {
     Int_t region_id = region->region();
  
     me_residual_x_[region_id] =
         bookHist1D(booker, region_id, "residual_x", "Residual in X", 120, -3, 3, "Residual in X [cm]");
  
     me_residual_y_[region_id] =
         bookHist1D(booker, region_id, "residual_y", "Residual in Y", 600, -15, 15, "Residual in Y [cm]");
  
     if (detail_plot_) {
       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();
             ME3IdsKey key3{region_id, station_id, layer_id};
  
             // Occupancy histograms of SimHits and RecHits for Efficiency
             me_detail_residual_x_[key3] =
                 bookHist1D(booker, key3, "residual_x", "Residual in X", 120, -3, 3, "Residual in X [cm]");
  
             me_detail_residual_y_[key3] =
                 bookHist1D(booker, key3, "residual_y", "Residual in Y", 600, -15, 15, "Residual in Y [cm]");
  
           }  // chamber loop
         }    // end if
       }      // station loop
     }        // detail_plot
   }          // region loop
  
   // NOTE Pull
   booker.setCurrentFolder("MuonGEMRecHitsV/GEMRecHitsTask/Pull");
  
   for (const auto& region : gem->regions()) {
     Int_t region_id = region->region();
  
     me_pull_x_[region_id] = bookHist1D(booker, region_id, "pull_x", "Pull in X", 100, -3, 3);
  
     me_pull_y_[region_id] = bookHist1D(booker, region_id, "pull_y", "Pull in Y", 100, -3, 3);
  
     if (detail_plot_) {
       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();
             ME3IdsKey key3{region_id, station_id, layer_id};
  
             me_detail_pull_x_[key3] = bookHist1D(booker, key3, "pull_x", "Pull in X", 100, -3, 3);
  
             me_detail_pull_y_[key3] = bookHist1D(booker, key3, "pull_y", "Pull in Y", 100, -3, 3);
  
           }  // chamber loop
         }    // end if
       }      // station loop
     }        // detail plot
   }          // region loop
  
   // NOTE Occupancy
   booker.setCurrentFolder("MuonGEMRecHitsV/GEMRecHitsTask/Occupancy");
   for (const auto& region : gem->regions()) {
     Int_t region_id = region->region();
  
     me_occ_zr_[region_id] = bookZROccupancy(booker, region_id, "rechit", "RecHit");
  
     me_simhit_occ_eta_[region_id] = bookHist1D(booker,
                                                region_id,
                                                "muon_simhit_occ_eta",
                                                "Muon SimHit Eta Occupancy",
                                                50,
                                                eta_range_[0],
                                                eta_range_[1],
                                                "|#eta|");
  
     me_rechit_occ_eta_[region_id] = bookHist1D(booker,
                                                region_id,
                                                "matched_rechit_occ_eta",
                                                "Matched RecHit Eta Occupancy",
                                                50,
                                                eta_range_[0],
                                                eta_range_[1],
                                                "|#eta|");
  
     for (const auto& station : region->stations()) {
       Int_t station_id = station->station();
       ME2IdsKey key2{region_id, station_id};
  
       me_simhit_occ_phi_[key2] =
           bookHist1D(booker, key2, "muon_simhit_occ_phi", "Muon SimHit Phi Occupancy", 51, -M_PI, M_PI, "#phi");
  
       me_rechit_occ_phi_[key2] =
           bookHist1D(booker, key2, "matched_rechit_occ_phi", "Matched RecHit Phi Occupancy", 51, -M_PI, M_PI, "#phi");
  
       me_simhit_occ_det_[key2] = bookDetectorOccupancy(booker, key2, station, "muon_simhit", "Muon SimHit");
  
       me_rechit_occ_det_[key2] = bookDetectorOccupancy(booker, key2, station, "matched_rechit", "Matched RecHit");
  
       if (detail_plot_) {
         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};
  
             me_detail_occ_xy_[key3] = bookXYOccupancy(booker, key3, "rechit", "RecHit");
  
             me_detail_occ_polar_[key3] = bookPolarOccupancy(booker, key3, "rechit", "RecHit");
  
           }  // chamber loop
         }    // end if
       }      // detail plot
     }        // station loop
   }          // region_loop
 }

◆ matchRecHitAgainstSimHit()

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

private

Definition at line 177 of file GEMRecHitValidation.cc.

                                                                                                                  {
   Bool_t matched = false;
  
   Int_t cls = rechit->clusterSize();
   Int_t rechit_first_strip = rechit->firstClusterStrip();
  
   if (cls == 1) {
     matched = simhit_strip == rechit_first_strip;
   } else {
     Int_t rechit_last_strip = rechit_first_strip + cls - 1;
     matched = (simhit_strip >= rechit_first_strip) and (simhit_strip <= rechit_last_strip);
   }
  
   return matched;
 }