#include <SiPixelPhase1RecHitsV.h>

Inheritance diagram for SiPixelPhase1RecHitsV:

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

	SiPixelPhase1RecHitsV (const edm::ParameterSet &conf)

Public Member Functions inherited from SiPixelPhase1Base
void	analyze (edm::Event const &e, edm::EventSetup const &eSetup) override=0

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

	SiPixelPhase1Base (const edm::ParameterSet &iConfig)

	~SiPixelPhase1Base () override

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

void	beginStream (edm::StreamID id) final

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

	DQMEDAnalyzer (void)

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

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

	~EDAnalyzerBase () override

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)

Public Member Functions inherited from HistogramManagerHolder
	HistogramManagerHolder (const edm::ParameterSet &iConfig)

Private Types
enum	{ IN_TIME_BUNCH, OUT_TIME_BUNCH, NSIMHITS, RECHIT_X, RECHIT_Y, RES_X, RES_Y, ERROR_X, ERROR_Y, PULL_X, PULL_Y }

Private Attributes
edm::EDGetTokenT< SiPixelRecHitCollection >	srcToken_

TrackerHitAssociator::Config	trackerHitAssociatorConfig_

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 SiPixelPhase1Base
bool	checktrigger (const edm::Event &iEvent, const edm::EventSetup &iSetup, const unsigned trgidx) const

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)

Protected Attributes inherited from HistogramManagerHolder
GeometryInterface	geometryInterface

std::vector< HistogramManager >	histo

Detailed Description

Definition at line 16 of file SiPixelPhase1RecHitsV.h.

Member Enumeration Documentation

anonymous enum

private

Enumerator
IN_TIME_BUNCH
OUT_TIME_BUNCH
NSIMHITS
RECHIT_X
RECHIT_Y
RES_X
RES_Y
ERROR_X
ERROR_Y
PULL_X
PULL_Y

Definition at line 17 of file SiPixelPhase1RecHitsV.h.

        {
     IN_TIME_BUNCH,
     OUT_TIME_BUNCH,
     NSIMHITS,
     RECHIT_X,
     RECHIT_Y,
     RES_X,
     RES_Y,
     ERROR_X,
     ERROR_Y,
     PULL_X,
     PULL_Y,
   };

Constructor & Destructor Documentation

SiPixelPhase1RecHitsV::SiPixelPhase1RecHitsV ( const edm::ParameterSet & conf )

explicit

Definition at line 19 of file SiPixelPhase1RecHitsV.cc.

                                                                            :
   SiPixelPhase1Base(iConfig),
   trackerHitAssociatorConfig_( iConfig, consumesCollector() ),
   srcToken_ ( consumes<SiPixelRecHitCollection>(iConfig.getParameter<edm::InputTag>("src")) )
 {}

Member Function Documentation

void SiPixelPhase1RecHitsV::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

Definition at line 25 of file SiPixelPhase1RecHitsV.cc.

References funct::abs(), TrackerHitAssociator::associateHit(), edmNew::DetSetVector< T >::begin(), EncodedEventId::bunchCrossing(), DEFINE_FWK_MODULE, edmNew::DetSetVector< T >::end(), PSimHit::entryPoint(), ERROR_X, ERROR_Y, EncodedEventId::event(), PSimHit::eventId(), PSimHit::exitPoint(), edm::Event::getByToken(), HistogramManagerHolder::histo, IN_TIME_BUNCH, input, edm::HandleBase::isValid(), funct::m, NSIMHITS, OUT_TIME_BUNCH, PULL_X, PULL_Y, RECHIT_X, RECHIT_Y, RES_X, RES_Y, rpcPointValidation_cfi::simHit, mathSSE::sqrt(), srcToken_, trackerHitAssociatorConfig_, PV3DBase< T, PVType, FrameType >::x(), LocalError::xx(), PV3DBase< T, PVType, FrameType >::y(), and LocalError::yy().

                                                                                        {
   edm::Handle<SiPixelRecHitCollection> input;
   iEvent.getByToken(srcToken_, input);
   if (!input.isValid()) return;
 
   TrackerHitAssociator associate(iEvent, trackerHitAssociatorConfig_);
 
   SiPixelRecHitCollection::const_iterator it;
   for (it = input->begin(); it != input->end(); ++it) {
     auto id = DetId(it->detId());
 
     for(SiPixelRecHit const& rechit : *it) {
       SiPixelRecHit::ClusterRef const& clust = rechit.cluster();
 
       std::vector<PSimHit> associateSimHit;
       associateSimHit = associate.associateHit(rechit);
       std::vector<PSimHit>::const_iterator closestIt = associateSimHit.begin();
 
       LocalPoint lp = rechit.localPosition();
       float rechit_x = lp.x();
       float rechit_y = lp.y();
       
       LocalError lerr = rechit.localPositionError();
       float lerr_x = sqrt(lerr.xx());
       float lerr_y = sqrt(lerr.yy());
 
       // loop over associated sim hits and find the closest
       if ( !associateSimHit.empty() ) {
       float closestSimHit = 9999.9;
 
        for (std::vector<PSimHit>::const_iterator m = associateSimHit.begin(); m < associateSimHit.end(); m++) {
       float sim_x1 ( (*m).entryPoint().x() ), sim_x2 ( (*m).exitPoint().x() ), sim_xpos ( 0.5*(sim_x1+sim_x2) );
       float sim_y1 ( (*m).entryPoint().y() ), sim_y2 ( (*m).exitPoint().y() ), sim_ypos ( 0.5*(sim_y1+sim_y2) );
 
           float xres ( std::abs(sim_xpos - rechit_x) ), yres ( std::abs(sim_ypos - rechit_y) );
           float dist = std::sqrt( xres*xres + yres*yres );
 
           if ( dist < closestSimHit ) {
             closestSimHit = dist;
             closestIt = m;
           }
         }
       }
 
       // Sim Hit stuff
       
       if ( !associateSimHit.empty() ) {
 
         const PSimHit& simHit = *closestIt;
       
         int bunch = simHit.eventId().bunchCrossing();
         int event = simHit.eventId().event();
 
         float sim_x1 ( simHit.entryPoint().x() ), sim_x2 ( simHit.exitPoint().x() ), sim_xpos ( 0.5*(sim_x1 + sim_x2) );
         float sim_y1 ( simHit.entryPoint().y() ), sim_y2 ( simHit.exitPoint().y() ), sim_ypos ( 0.5*(sim_y1 + sim_y2) );
 
         float res_x = (rechit_x - sim_xpos) * 10000.0;
         float res_y = (rechit_y - sim_ypos) * 10000.0;
 
         float pull_x = ( rechit_x - sim_xpos ) / lerr_x;
         float pull_y = ( rechit_y - sim_ypos ) / lerr_y;
   
         // Now Plotting stuff
         if ( bunch == 0 ) histo[IN_TIME_BUNCH].fill(bunch, id, &iEvent);
         if ( bunch != 0 ) histo[OUT_TIME_BUNCH].fill(bunch, id, &iEvent);
  
         histo[NSIMHITS].fill(event, id, &iEvent);
 
         histo[RECHIT_X].fill(rechit_x, id, &iEvent);
         histo[RECHIT_Y].fill(rechit_y, id, &iEvent);
 
         histo[RES_X].fill(res_x, id, &iEvent);
         histo[RES_Y].fill(res_y, id, &iEvent);
 
         histo[ERROR_X].fill(lerr_x, id, &iEvent);
         histo[ERROR_Y].fill(lerr_y, id, &iEvent);
 
         histo[PULL_X].fill(pull_x, id, &iEvent);
         histo[PULL_Y].fill(pull_y, id, &iEvent); 
       }
     }
   }
 }

Member Data Documentation

edm::EDGetTokenT<SiPixelRecHitCollection> SiPixelPhase1RecHitsV::srcToken_

private

Definition at line 37 of file SiPixelPhase1RecHitsV.h.

Referenced by analyze().

TrackerHitAssociator::Config SiPixelPhase1RecHitsV::trackerHitAssociatorConfig_

private

Definition at line 36 of file SiPixelPhase1RecHitsV.h.

Referenced by analyze().

Public Member Functions

Private Types

Private Attributes

Additional Inherited Members

Detailed Description

Member Enumeration Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation