SiPixelPhase1RecHits Class Reference

#include <SiPixelPhase1RecHits.h>

Inheritance diagram for SiPixelPhase1RecHits:

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &)
	SiPixelPhase1RecHits (const edm::ParameterSet &conf)
Public Member Functions inherited from SiPixelPhase1Base
void	analyze (edm::Event const &e, edm::EventSetup const &eSetup)=0
void	bookHistograms (DQMStore::IBooker &iBooker, edm::Run const &run, edm::EventSetup const &iSetup)
	SiPixelPhase1Base (const edm::ParameterSet &iConfig)
virtual	~SiPixelPhase1Base ()
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
	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 &	itemsToGetFromEvent () 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	{   NRECHITS, CLUST_X, CLUST_Y, ERROR_X,   ERROR_Y, POS, CLUSTER_PROB }

Private Attributes
bool	onlyValid_
edm::EDGetTokenT< reco::TrackCollection >	srcToken_

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

Member Enumeration Documentation

anonymous enum

private

Enumerator
NRECHITS
CLUST_X
CLUST_Y
ERROR_X
ERROR_Y
POS
CLUSTER_PROB

Definition at line 16 of file SiPixelPhase1RecHits.h.

       {
    NRECHITS,
    CLUST_X,
    CLUST_Y,
    ERROR_X,
    ERROR_Y,
    POS,
    CLUSTER_PROB
  };

Constructor & Destructor Documentation

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

explicit

Definition at line 27 of file SiPixelPhase1RecHits.cc.

References edm::ParameterSet::getParameter(), onlyValid_, and srcToken_.

                                                                         :
  SiPixelPhase1Base(iConfig) 
{
  srcToken_ = consumes<reco::TrackCollection>(iConfig.getParameter<edm::InputTag>("src"));
  onlyValid_=iConfig.getParameter<bool>("onlyValidHits");
}

Member Function Documentation

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

Definition at line 34 of file SiPixelPhase1RecHits.cc.

References CLUST_X, CLUST_Y, SiPixelRecHit::cluster(), CLUSTER_PROB, SiPixelRecHit::clusterProbability(), cuy::col, DEFINE_FWK_MODULE, ERROR_X, ERROR_Y, edm::EventSetup::get(), edm::Event::getByToken(), h, hcalSimParameters_cfi::hb, HistogramManagerHolder::histo, TrackerGeometry::idToDet(), createfilelist::int, edm::ESHandleBase::isValid(), Topology::measurementPosition(), NRECHITS, onlyValid_, PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, POS, RecoTauValidation_cfi::sizeX, RecoTauValidation_cfi::sizeY, PixelGeomDetUnit::specificTopology(), mathSSE::sqrt(), srcToken_, HiIsolationCommonParameters_cff::track, mixOne_simraw_on_sim_cfi::tracker, l1t::tracks, TrackingRecHit::valid, PV2DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::x(), LocalError::xx(), PV2DBase< T, PVType, FrameType >::y(), PV3DBase< T, PVType, FrameType >::y(), and LocalError::yy().

                                                                                      {

  edm::ESHandle<TrackerGeometry> tracker;
  iSetup.get<TrackerDigiGeometryRecord>().get(tracker);
  assert(tracker.isValid());

  edm::Handle<reco::TrackCollection> tracks;
  iEvent.getByToken( srcToken_, tracks);
  if (!tracks.isValid()) return;

  for (auto const & track : *tracks) {

    bool isBpixtrack = false, isFpixtrack = false;

    auto const & trajParams = track.extra()->trajParams();
    auto hb = track.recHitsBegin();
    for(unsigned int h=0;h<track.recHitsSize();h++){
     
      auto hit = *(hb+h);
      if(!hit->isValid()) continue;
      
      DetId id = hit->geographicalId();
      uint32_t subdetid = (id.subdetId());

      if (subdetid == PixelSubdetector::PixelBarrel) isBpixtrack = true;
      if (subdetid == PixelSubdetector::PixelEndcap) isFpixtrack = true;
    }

    if (!isBpixtrack && !isFpixtrack) continue;

    // then, look at each hit
    for(unsigned int h=0;h<track.recHitsSize();h++){
      auto rechit = *(hb+h);
    
      if(!rechit->isValid()) continue;

      //continue if not a Pixel recHit
      DetId id = rechit->geographicalId();
      uint32_t subdetid = (id.subdetId());

      if (   subdetid != PixelSubdetector::PixelBarrel 
         && subdetid != PixelSubdetector::PixelEndcap) continue;

      bool isHitValid   = rechit->getType()==TrackingRecHit::valid;
      if (onlyValid_ && !isHitValid) continue; //useful to run on cosmics where the TrackEfficiency plugin is not used

      const SiPixelRecHit* prechit = dynamic_cast<const SiPixelRecHit*>(rechit);//to be used to get the associated cluster and the cluster probability      

      int sizeX=0, sizeY=0;

      if (isHitValid){
    SiPixelRecHit::ClusterRef const& clust = prechit->cluster();
    sizeX = (*clust).sizeX();
    sizeY = (*clust).sizeY();
      }

      const PixelGeomDetUnit* geomdetunit = dynamic_cast<const PixelGeomDetUnit*> ( tracker->idToDet(id) );
      const PixelTopology& topol = geomdetunit->specificTopology();

      LocalPoint lp = trajParams[h].position();
      MeasurementPoint mp = topol.measurementPosition(lp);
      
      int row = (int) mp.x();
      int col = (int) mp.y();
      
      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());

      histo[NRECHITS].fill(id, &iEvent, col, row); //in general a inclusive counter of missing/valid/inactive hits

      if (isHitValid){
    histo[CLUST_X].fill(sizeX, id, &iEvent, col, row);
    histo[CLUST_Y].fill(sizeY, id, &iEvent, col, row);
      }

      histo[ERROR_X].fill(lerr_x, id, &iEvent);
      histo[ERROR_Y].fill(lerr_y, id, &iEvent);

      histo[POS].fill(rechit_x, rechit_y, id, &iEvent);
      
      if (isHitValid){
    double clusterProbability= prechit->clusterProbability(0);
    if (clusterProbability > 0)
      histo[CLUSTER_PROB].fill(log10(clusterProbability), id, &iEvent);
      }
    }
  }

  histo[NRECHITS].executePerEventHarvesting(&iEvent);
}

Member Data Documentation

bool SiPixelPhase1RecHits::onlyValid_

private

Definition at line 26 of file SiPixelPhase1RecHits.h.

Referenced by analyze(), and SiPixelPhase1RecHits().

edm::EDGetTokenT<reco::TrackCollection> SiPixelPhase1RecHits::srcToken_

private

Definition at line 33 of file SiPixelPhase1RecHits.h.

Referenced by analyze(), and SiPixelPhase1RecHits().