#include <RPCPointVsRecHit.h>

Inheritance diagram for RPCPointVsRecHit:

Public Member Functions
void	analyze (const edm::Event &event, const edm::EventSetup &eventSetup) override

void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override

	RPCPointVsRecHit (const edm::ParameterSet &pset)

	~RPCPointVsRecHit ()

Public Member Functions inherited from DQMEDAnalyzer
virtual void	beginRun (edm::Run const &, edm::EventSetup const &) final

virtual void	beginStream (edm::StreamID id) final

virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)

	DQMEDAnalyzer (void)

virtual void	endLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, dqmDetails::NoCache *) const final

virtual void	endRunSummary (edm::Run const &, edm::EventSetup const &, dqmDetails::NoCache *) const final

uint32_t	streamId () const

Public Member Functions inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
	EDAnalyzer ()=default

Public Member Functions inherited from edm::stream::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

virtual	~EDAnalyzerBase ()

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

virtual	~EDConsumerBase () noexcept(false)

Private Attributes
RPCValidHistograms	h_

edm::EDGetTokenT< RPCRecHitCollection >	recHitToken_

edm::EDGetTokenT< RPCRecHitCollection >	refHitToken_

std::string	subDir_

Additional Inherited Members
Public Types inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
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

Public Types inherited from edm::stream::EDAnalyzerBase
typedef EDAnalyzerAdaptorBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from DQMEDAnalyzer
static std::shared_ptr< dqmDetails::NoCache >	globalBeginLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const *)

static std::shared_ptr< dqmDetails::NoCache >	globalBeginRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const *)

static void	globalEndLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const , dqmDetails::NoCache )

static void	globalEndRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const , dqmDetails::NoCache )

Static Public Member Functions inherited from edm::stream::EDAnalyzerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Definition at line 19 of file RPCPointVsRecHit.h.

Constructor & Destructor Documentation

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

Definition at line 15 of file RPCPointVsRecHit.cc.

References edm::ParameterSet::getParameter(), and AlCaHLTBitMon_QueryRunRegistry::string.

 {
   refHitToken_ = consumes<RPCRecHitCollection>(pset.getParameter<edm::InputTag>("refHit"));
   recHitToken_ = consumes<RPCRecHitCollection>(pset.getParameter<edm::InputTag>("recHit"));
 
   subDir_ = pset.getParameter<std::string>("subDir");
 }

RPCPointVsRecHit::~RPCPointVsRecHit ( )

inline

Definition at line 23 of file RPCPointVsRecHit.h.

References analyze(), and bookHistograms().

23 {};

Member Function Documentation

void RPCPointVsRecHit::analyze	(	const edm::Event &	event,
		const edm::EventSetup &	eventSetup
	)

override

Definition at line 23 of file RPCPointVsRecHit.cc.

References fftjetvertexadder_cfi::errX, edm::EventSetup::get(), edm::Event::getByToken(), RPCRoll::id(), match(), RPCDetId::region(), relativeConstraints::ring, RPCDetId::ring(), RPCGeometry::roll(), mathSSE::sqrt(), relativeConstraints::station, and RPCDetId::station().

Referenced by ~RPCPointVsRecHit().

 {
   // Get the RPC Geometry
   edm::ESHandle<RPCGeometry> rpcGeom;
   eventSetup.get<MuonGeometryRecord>().get(rpcGeom);
 
   // Retrieve RefHits from the event
   edm::Handle<RPCRecHitCollection> refHitHandle;
   if ( !event.getByToken(refHitToken_, refHitHandle) )
   {
     edm::LogInfo("RPCPointVsRecHit") << "Cannot find reference hit collection\n";
     return;
   }
 
   // Retrieve RecHits from the event
   edm::Handle<RPCRecHitCollection> recHitHandle;
   if ( !event.getByToken(recHitToken_, recHitHandle) )
   {
     edm::LogInfo("RPCPointVsRecHit") << "Cannot find recHit collection\n";
     return;
   }
 
   typedef RPCRecHitCollection::const_iterator RecHitIter;
 
   // Loop over refHits, fill histograms which does not need associations
   int nRefHitBarrel = 0, nRefHitEndcap = 0;
   for ( RecHitIter refHitIter = refHitHandle->begin();
         refHitIter != refHitHandle->end(); ++refHitIter )
   {
     const RPCDetId detId = static_cast<const RPCDetId>(refHitIter->rpcId());
     const RPCRoll* roll = dynamic_cast<const RPCRoll*>(rpcGeom->roll(detId()));
     if ( !roll ) continue;
 
     const int region = roll->id().region();
     const int ring = roll->id().ring();
     //const int sector = roll->id().sector();
     const int station = roll->id().station();
     //const int layer = roll->id().layer();
     //const int subSector = roll->id().subsector();
 
     if ( region == 0 )
     {
       h_.refHitOccupancyBarrel_wheel->Fill(ring);
       h_.refHitOccupancyBarrel_station->Fill(station);
       h_.refHitOccupancyBarrel_wheel_station->Fill(ring, station);
     }
     else
     {
       h_.refHitOccupancyEndcap_disk->Fill(region*station);
       h_.refHitOccupancyEndcap_disk_ring->Fill(region*station, ring);
     }
 
   }
   h_.nRefHitBarrel->Fill(nRefHitBarrel);
   h_.nRefHitEndcap->Fill(nRefHitEndcap);
 
   // Loop over recHits, fill histograms which does not need associations
   int sumClusterSizeBarrel = 0, sumClusterSizeEndcap = 0;
   int nRecHitBarrel = 0, nRecHitEndcap = 0;
   for ( RecHitIter recHitIter = recHitHandle->begin();
         recHitIter != recHitHandle->end(); ++recHitIter )
   {
     const RPCDetId detId = static_cast<const RPCDetId>(recHitIter->rpcId());
     const RPCRoll* roll = dynamic_cast<const RPCRoll*>(rpcGeom->roll(detId()));
     if ( !roll ) continue;
 
     const int region = roll->id().region();
     const int ring = roll->id().ring();
     //const int sector = roll->id().sector();
     const int station = (roll->id().station());
     //const int layer = roll->id().layer();
     //const int subSector = roll->id().subsector();
 
     h_.clusterSize->Fill(recHitIter->clusterSize());
 
     if ( region == 0 )
     {
       ++nRecHitBarrel;
       sumClusterSizeBarrel += recHitIter->clusterSize();
       h_.clusterSizeBarrel->Fill(recHitIter->clusterSize());
       h_.recHitOccupancyBarrel_wheel->Fill(ring);
       h_.recHitOccupancyBarrel_station->Fill(station);
       h_.recHitOccupancyBarrel_wheel_station->Fill(ring, station);
     }
     else
     {
       ++nRecHitEndcap;
       sumClusterSizeEndcap += recHitIter->clusterSize();
       h_.clusterSizeEndcap->Fill(recHitIter->clusterSize());
       h_.recHitOccupancyEndcap_disk->Fill(region*station);
       h_.recHitOccupancyEndcap_disk_ring->Fill(ring, region*station);
     }
 
   }
   const double nRecHit = nRecHitBarrel+nRecHitEndcap;
   h_.nRecHitBarrel->Fill(nRecHitBarrel);
   h_.nRecHitEndcap->Fill(nRecHitEndcap);
   if ( nRecHit > 0 )
   {
     const int sumClusterSize = sumClusterSizeBarrel+sumClusterSizeEndcap;
     h_.avgClusterSize->Fill(double(sumClusterSize)/nRecHit);
 
     if ( nRecHitBarrel > 0 )
     {
       h_.avgClusterSizeBarrel->Fill(double(sumClusterSizeBarrel)/nRecHitBarrel);
     }
     if ( nRecHitEndcap > 0 )
     {
       h_.avgClusterSizeEndcap->Fill(double(sumClusterSizeEndcap)/nRecHitEndcap);
     }
   }
 
   // Start matching RefHits to RecHits
   typedef std::map<RecHitIter, RecHitIter> RecToRecHitMap;
   RecToRecHitMap refToRecHitMap;
 
   for ( RecHitIter refHitIter = refHitHandle->begin();
         refHitIter != refHitHandle->end(); ++refHitIter )
   {
     const RPCDetId refDetId = static_cast<const RPCDetId>(refHitIter->rpcId());
     const RPCRoll* refRoll = dynamic_cast<const RPCRoll*>(rpcGeom->roll(refDetId));
     if ( !refRoll ) continue;
 
     const double refX = refHitIter->localPosition().x();
 
     for ( RecHitIter recHitIter = recHitHandle->begin();
           recHitIter != recHitHandle->end(); ++recHitIter )
     {
       const RPCDetId recDetId = static_cast<const RPCDetId>(recHitIter->rpcId());
       const RPCRoll* recRoll = dynamic_cast<const RPCRoll*>(rpcGeom->roll(recDetId));
       if ( !recRoll ) continue;
 
       if ( refDetId != recDetId ) continue;
 
       const double recX = recHitIter->localPosition().x();
       const double newDx = fabs(recX - refX);
 
       // Associate RefHit to RecHit
       RecToRecHitMap::const_iterator prevRefToReco = refToRecHitMap.find(refHitIter);
       if ( prevRefToReco == refToRecHitMap.end() )
       {
         refToRecHitMap.insert(std::make_pair(refHitIter, recHitIter));
       }
       else
       {
         const double oldDx = fabs(prevRefToReco->second->localPosition().x() - refX);
 
         if ( newDx < oldDx )
         {
           refToRecHitMap[refHitIter] = recHitIter;
         }
       }
     }
   }
 
   // Now we have refHit-recHit mapping
   // So we can fill up relavant histograms
   for ( RecToRecHitMap::const_iterator match = refToRecHitMap.begin();
         match != refToRecHitMap.end(); ++match )
   {
     RecHitIter refHitIter = match->first;
     RecHitIter recHitIter = match->second;
 
     const RPCDetId detId = static_cast<const RPCDetId>(refHitIter->rpcId());
     const RPCRoll* roll = dynamic_cast<const RPCRoll*>(rpcGeom->roll(detId));
 
     const int region = roll->id().region();
     const int ring = roll->id().ring();
     //const int sector = roll->id().sector();
     const int station = roll->id().station();
     //const int layer = roll->id().layer();
     //const int subsector = roll->id().subsector();
 
     const double refX = refHitIter->localPosition().x();
     const double recX = recHitIter->localPosition().x();
     const double errX = sqrt(recHitIter->localPositionError().xx());
     const double dX = recX - refX;
     const double pull = errX == 0 ? -999 : dX/errX;
 
     //const GlobalPoint refPos = roll->toGlobal(refHitIter->localPosition());
     //const GlobalPoint recPos = roll->toGlobal(recHitIter->localPosition());
 
     if ( region == 0 )
     {
       h_.resBarrel->Fill(dX);
       h_.pullBarrel->Fill(pull);
       h_.matchOccupancyBarrel_wheel->Fill(ring);
       h_.matchOccupancyBarrel_station->Fill(station);
       h_.matchOccupancyBarrel_wheel_station->Fill(ring, station);
 
       h_.res_wheel_res->Fill(ring, dX);
       h_.res_station_res->Fill(station, dX);
       h_.pull_wheel_pull->Fill(ring, pull);
       h_.pull_station_pull->Fill(station, pull);
     }
     else
     {
       h_.resEndcap->Fill(dX);
       h_.pullEndcap->Fill(pull);
       h_.matchOccupancyEndcap_disk->Fill(region*station);
       h_.matchOccupancyEndcap_disk_ring->Fill(region*station, ring);
 
       h_.res_disk_res->Fill(region*station, dX);
       h_.res_ring_res->Fill(ring, dX);
       h_.pull_disk_pull->Fill(region*station, pull);
       h_.pull_ring_pull->Fill(ring, pull);
     }
   }
 
 }

void RPCPointVsRecHit::bookHistograms	(	DQMStore::IBooker &	booker,
		edm::Run const &	run,
		edm::EventSetup const &	eventSetup
	)

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 234 of file RPCPointVsRecHit.cc.

References DEFINE_FWK_MODULE.

Referenced by ~RPCPointVsRecHit().

 {
   // Book MonitorElements
   h_.bookHistograms(booker, subDir_);
 }

Member Data Documentation

RPCValidHistograms RPCPointVsRecHit::h_

private

Definition at line 32 of file RPCPointVsRecHit.h.

edm::EDGetTokenT<RPCRecHitCollection> RPCPointVsRecHit::recHitToken_

private

Definition at line 29 of file RPCPointVsRecHit.h.

edm::EDGetTokenT<RPCRecHitCollection> RPCPointVsRecHit::refHitToken_

private

Definition at line 29 of file RPCPointVsRecHit.h.

std::string RPCPointVsRecHit::subDir_

private

Definition at line 31 of file RPCPointVsRecHit.h.

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation