ME0RecHitsValidation Class Reference

#include <ME0RecHitsValidation.h>

Inheritance diagram for ME0RecHitsValidation:

Public Member Functions
void	analyze (const edm::Event &e, const edm::EventSetup &) override
void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
	ME0RecHitsValidation (const edm::ParameterSet &)
	~ME0RecHitsValidation () override
Public Member Functions inherited from ME0BaseValidation
MonitorElement *	BookHistXY (DQMStore::IBooker &, const char *name, const char *label, unsigned int region_num, unsigned int layer_num=99)
MonitorElement *	BookHistZR (DQMStore::IBooker &, const char *name, const char *label, unsigned int region_num, unsigned int layer_num=99)
	ME0BaseValidation (const edm::ParameterSet &ps)
	~ME0BaseValidation () override
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	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndRuns () const final

Private Attributes
edm::EDGetToken	InputTagToken_RecHit
MonitorElement *	me0_rh_DeltaX [2][6]
MonitorElement *	me0_rh_DeltaY [2][6]
MonitorElement *	me0_rh_PullX [2][6]
MonitorElement *	me0_rh_PullY [2][6]
MonitorElement *	me0_rh_xy [2][6]
MonitorElement *	me0_rh_zr [2]
Int_t	npart

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 DQMEDAnalyzer
uint64_t	meId () const
Protected Attributes inherited from ME0BaseValidation
edm::EDGetToken	InputTagToken_
std::vector< std::string >	layerLabel
int	nBinXY_
std::vector< double >	nBinZR_
std::vector< double >	RangeZR_
std::vector< std::string >	regionLabel
Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_
edm::EDPutTokenT< DQMToken >	runToken_
unsigned int	streamId_

Detailed Description

Definition at line 11 of file ME0RecHitsValidation.h.

Constructor & Destructor Documentation

ME0RecHitsValidation()

ME0RecHitsValidation::ME0RecHitsValidation ( const edm::ParameterSet &  cfg )

explicit

Definition at line 4 of file ME0RecHitsValidation.cc.

                                                                     : ME0BaseValidation(cfg) {
  InputTagToken_ = consumes<edm::PSimHitContainer>(cfg.getParameter<edm::InputTag>("simInputLabel"));
  InputTagToken_RecHit = consumes<ME0RecHitCollection>(cfg.getParameter<edm::InputTag>("recHitInputLabel"));
}

References looper::cfg, ME0BaseValidation::InputTagToken_, and InputTagToken_RecHit.

~ME0RecHitsValidation()

ME0RecHitsValidation::~ME0RecHitsValidation ( )

override

Definition at line 59 of file ME0RecHitsValidation.cc.

{}

Member Function Documentation

analyze()

void ME0RecHitsValidation::analyze ( const edm::Event &  e, const edm::EventSetup &  iSetup  )

overridevirtual

Implements ME0BaseValidation.

Definition at line 61 of file ME0RecHitsValidation.cc.

                                                                                 {
  edm::ESHandle<ME0Geometry> hGeom;
  iSetup.get<MuonGeometryRecord>().get(hGeom);
  const ME0Geometry *ME0Geometry_ = (&*hGeom);
  edm::Handle<edm::PSimHitContainer> ME0Hits;
  e.getByToken(InputTagToken_, ME0Hits);
 
  edm::Handle<ME0RecHitCollection> ME0RecHits;
  e.getByToken(InputTagToken_RecHit, ME0RecHits);
 
  if (!ME0Hits.isValid() | !ME0RecHits.isValid()) {
    edm::LogError("ME0RecHitsValidation") << "Cannot get ME0Hits/ME0RecHits by Token simInputTagToken";
    return;
  }
 
  for (auto hits = ME0Hits->begin(); hits != ME0Hits->end(); hits++) {
    const ME0DetId id(hits->detUnitId());
    Int_t sh_region = id.region();
    Int_t sh_layer = id.layer();
    Int_t sh_station = 0;
    Int_t sh_chamber = id.chamber();
    Int_t sh_roll = id.roll();
 
    // Int_t even_odd = id.chamber()%2;
    if (ME0Geometry_->idToDet(hits->detUnitId()) == nullptr) {
      std::cout << "simHit did not matched with GEMGeometry." << std::endl;
      continue;
    }
 
    const LocalPoint hitLP(hits->localPosition());
 
    for (ME0RecHitCollection::const_iterator recHit = ME0RecHits->begin(); recHit != ME0RecHits->end(); ++recHit) {
      Float_t x = recHit->localPosition().x();
      Float_t xErr = recHit->localPositionError().xx();
      Float_t y = recHit->localPosition().y();
      Float_t yErr = recHit->localPositionError().yy();
      // Float_t detId = (Short_t) (*recHit).me0Id();
      Float_t bx = 0;
      Float_t clusterSize = 0;
      Float_t firstClusterStrip = 0;
 
      ME0DetId id((*recHit).me0Id());
 
      Short_t region = (Short_t)id.region();
      Short_t station = 0;
      Short_t layer = (Short_t)id.layer();
      Short_t chamber = (Short_t)id.chamber();
      Short_t roll = (Short_t)id.roll();
 
      LocalPoint rhLP = recHit->localPosition();
      GlobalPoint rhGP = ME0Geometry_->idToDet((*recHit).me0Id())->surface().toGlobal(rhLP);
 
      Float_t globalR = rhGP.perp();
      Float_t globalX = rhGP.x();
      Float_t globalY = rhGP.y();
      Float_t globalZ = rhGP.z();
 
      Float_t x_sim = hitLP.x();
      Float_t y_sim = hitLP.y();
      Float_t pullX = (x_sim - x) / sqrt(xErr);
      Float_t pullY = (y_sim - y) / sqrt(yErr);
 
      if (bx != 0)
        continue;
 
      std::vector<int> stripsFired;
      for (int i = firstClusterStrip; i < (firstClusterStrip + clusterSize); i++) {
        stripsFired.push_back(i);
      }
 
      const bool cond1(sh_region == region and sh_layer == layer and sh_station == station);
      const bool cond2(sh_chamber == chamber and sh_roll == roll);
      const bool cond3(
          (std::find(stripsFired.begin(), stripsFired.end(), (firstClusterStrip + 1)) != stripsFired.end()) or
          clusterSize == 0);
 
      int region_num = 0;
      if (region == -1)
        region_num = 0;
      else if (region == 1)
        region_num = 1;
      int layer_num = layer - 1;
 
      if (cond1 and cond2 and cond3) {
        me0_rh_xy[region_num][layer_num]->Fill(globalX, globalY);
        me0_rh_zr[region_num]->Fill(globalZ, globalR);
 
        me0_rh_DeltaX[region_num][layer_num]->Fill(x_sim - x);
        me0_rh_DeltaY[region_num][layer_num]->Fill(y_sim - y);
        me0_rh_PullX[region_num][layer_num]->Fill(pullX);
        me0_rh_PullY[region_num][layer_num]->Fill(pullY);
      }
    }
  }
}

References simKBmtfDigis_cfi::bx, relativeConstraints::chamber, gather_cfg::cout, MillePedeFileConverter_cfg::e, dqm::impl::MonitorElement::Fill(), spr::find(), edm::EventSetup::get(), get, hfClusterShapes_cfi::hits, mps_fire::i, triggerObjects_cff::id, ME0Geometry::idToDet(), ME0BaseValidation::InputTagToken_, InputTagToken_RecHit, edm::HandleBase::isValid(), me0_rh_DeltaX, me0_rh_DeltaY, me0_rh_PullX, me0_rh_PullY, me0_rh_xy, me0_rh_zr, or, PV3DBase< T, PVType, FrameType >::perp(), rpcPointValidation_cfi::recHit, HLT_2018_cff::region, mathSSE::sqrt(), relativeConstraints::station, GeomDet::toGlobal(), x, PV3DBase< T, PVType, FrameType >::x(), y, PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

bookHistograms()

void ME0RecHitsValidation::bookHistograms ( DQMStore::IBooker &  ibooker, edm::Run const &  Run, edm::EventSetup const &  iSetup  )

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 9 of file ME0RecHitsValidation.cc.

                                                                       {
  LogDebug("MuonME0RecHitsValidation") << "Info : Loading Geometry information\n";
  ibooker.setCurrentFolder("MuonME0RecHitsV/ME0RecHitsTask");
 
  unsigned int nregion = 2;
 
  edm::LogInfo("MuonME0RecHitsValidation") << "+++ Info : # of region : " << nregion << std::endl;
 
  LogDebug("MuonME0RecHitsValidation") << "+++ Info : finish to get geometry information from ES.\n";
 
  for (unsigned int region_num = 0; region_num < nregion; region_num++) {
    me0_rh_zr[region_num] = BookHistZR(ibooker, "me0_rh_tot", "Digi", region_num);
    for (unsigned int layer_num = 0; layer_num < 6; layer_num++) {
      me0_rh_xy[region_num][layer_num] = BookHistXY(ibooker, "me0_rh", "RecHit", region_num, layer_num);
 
      std::string histo_name_DeltaX =
          std::string("me0_rh_DeltaX_r") + regionLabel[region_num] + "_l" + layerLabel[layer_num];
      std::string histo_name_DeltaY =
          std::string("me0_rh_DeltaY_r") + regionLabel[region_num] + "_l" + layerLabel[layer_num];
      std::string histo_label_DeltaX = "RecHit Delta X : region" + regionLabel[region_num] + " layer " +
                                       layerLabel[layer_num] + " " + " ; x_{SimHit} - x_{RecHit} ; entries";
      std::string histo_label_DeltaY = "RecHit Delta Y : region" + regionLabel[region_num] + " layer " +
                                       layerLabel[layer_num] + " " + " ; y_{SimHit} - y_{RecHit} ; entries";
 
      me0_rh_DeltaX[region_num][layer_num] =
          ibooker.book1D(histo_name_DeltaX.c_str(), histo_label_DeltaX.c_str(), 100, -10, 10);
      me0_rh_DeltaY[region_num][layer_num] =
          ibooker.book1D(histo_name_DeltaY.c_str(), histo_label_DeltaY.c_str(), 100, -10, 10);
 
      std::string histo_name_PullX =
          std::string("me0_rh_PullX_r") + regionLabel[region_num] + "_l" + layerLabel[layer_num];
      std::string histo_name_PullY =
          std::string("me0_rh_PullY_r") + regionLabel[region_num] + "_l" + layerLabel[layer_num];
      std::string histo_label_PullX = "RecHit Pull X : region" + regionLabel[region_num] + " layer " +
                                      layerLabel[layer_num] + " " +
                                      " ; #frac{x_{SimHit} - x_{RecHit}}{#sigma_{x,RecHit}} ; entries";
      std::string histo_label_PullY = "RecHit Pull Y : region" + regionLabel[region_num] + " layer " +
                                      layerLabel[layer_num] + " " +
                                      " ; #frac{y_{SimHit} - y_{RecHit}}{#sigma_{y,RecHit}} ; entries";
 
      me0_rh_PullX[region_num][layer_num] =
          ibooker.book1D(histo_name_PullX.c_str(), histo_label_DeltaX.c_str(), 100, -10, 10);
      me0_rh_PullY[region_num][layer_num] =
          ibooker.book1D(histo_name_PullY.c_str(), histo_label_DeltaY.c_str(), 100, -10, 10);
    }
  }
}

References dqm::implementation::IBooker::book1D(), ME0BaseValidation::BookHistXY(), ME0BaseValidation::BookHistZR(), ME0BaseValidation::layerLabel, LogDebug, me0_rh_DeltaX, me0_rh_DeltaY, me0_rh_PullX, me0_rh_PullY, me0_rh_xy, me0_rh_zr, ME0BaseValidation::regionLabel, dqm::implementation::NavigatorBase::setCurrentFolder(), and AlCaHLTBitMon_QueryRunRegistry::string.

Member Data Documentation

InputTagToken_RecHit

edm::EDGetToken ME0RecHitsValidation::InputTagToken_RecHit

private

Definition at line 27 of file ME0RecHitsValidation.h.

Referenced by analyze(), and ME0RecHitsValidation().

me0_rh_DeltaX

MonitorElement* ME0RecHitsValidation::me0_rh_DeltaX[2][6]

private

Definition at line 22 of file ME0RecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

me0_rh_DeltaY

MonitorElement* ME0RecHitsValidation::me0_rh_DeltaY[2][6]

private

Definition at line 23 of file ME0RecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

me0_rh_PullX

MonitorElement* ME0RecHitsValidation::me0_rh_PullX[2][6]

private

Definition at line 24 of file ME0RecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

me0_rh_PullY

MonitorElement* ME0RecHitsValidation::me0_rh_PullY[2][6]

private

Definition at line 25 of file ME0RecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

me0_rh_xy

MonitorElement* ME0RecHitsValidation::me0_rh_xy[2][6]

private

Definition at line 19 of file ME0RecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

me0_rh_zr

MonitorElement* ME0RecHitsValidation::me0_rh_zr[2]

private

Definition at line 20 of file ME0RecHitsValidation.h.

Referenced by analyze(), and bookHistograms().

npart

Int_t ME0RecHitsValidation::npart

private

Definition at line 29 of file ME0RecHitsValidation.h.