#include <PixelUnpackingRegions.h>

Classes
struct	Module

struct	Region

Public Member Functions
bool	mayUnpackFED (unsigned int fed_n) const
	check whether a FED has to be unpacked More...

bool	mayUnpackModule (unsigned int id) const
	check whether a module has to be unpacked More...

const std::set< unsigned int > *	modulesToUnpack () const
	full set of module ids to unpack More...

unsigned int	nBarrelFEDs () const

unsigned int	nBarrelModules () const

unsigned int	nFEDs () const
	various informational accessors: More...

unsigned int	nForwardFEDs () const

unsigned int	nForwardModules () const

unsigned int	nModules () const

unsigned int	nRegions () const

	PixelUnpackingRegions (const edm::ParameterSet &, edm::ConsumesCollector &&iC)

void	run (const edm::Event &e, const edm::EventSetup &es)
	has to be run during each event More...

	~PixelUnpackingRegions ()

Private Member Functions
void	addRegion (Region &r)

void	addRegionLocal (Region &r, std::vector< Module > &container, const Module &lo, const Module &hi)

void	gatherFromRange (Region &r, std::vector< Module >::const_iterator, std::vector< Module >::const_iterator)

void	initialize (const edm::EventSetup &es)
	run by the run method: (re)initialize the cabling data when it's necessary More...

Private Attributes
math::XYZPoint	beamSpot_

edm::InputTag	beamSpotTag_

std::unique_ptr < SiPixelFedCablingTree >	cabling_

std::vector< double >	dPhi_

std::set< unsigned int >	feds_

std::vector< edm::InputTag >	inputs_

std::vector< double >	maxZ_

std::set< unsigned int >	modules_

unsigned int	nreg_

std::vector< Module >	phiBPIX_

std::vector< Module >	phiFPIXm_

std::vector< Module >	phiFPIXp_

edm::EDGetTokenT< reco::BeamSpot >	tBeamSpot

std::vector< edm::EDGetTokenT < reco::CandidateView > >	tCandidateView

edm::ESWatcher < SiPixelFedCablingMapRcd >	watcherSiPixelFedCablingMap_

Detailed Description

Input: One or several collections of Candidate-based seeds with their objects defining the directions of unpacking regions; separate deltaPhi and maxZ tolerances could be given to each input collection. Output: FED ids and module detIds that need to be unpacked

Definition at line 30 of file PixelUnpackingRegions.h.

Constructor & Destructor Documentation

PixelUnpackingRegions::PixelUnpackingRegions	(	const edm::ParameterSet &	conf,
		edm::ConsumesCollector &&	iC
	)

Definition at line 41 of file PixelUnpackingRegions.cc.

References beamSpotTag_, dPhi_, edm::ParameterSet::getParameter(), inputs_, maxZ_, tree::t, tBeamSpot, and tCandidateView.

 {
   edm::ParameterSet regPSet = conf.getParameter<edm::ParameterSet>("Regions");
   beamSpotTag_ = regPSet.getParameter<edm::InputTag>("beamSpot");
   inputs_      = regPSet.getParameter<std::vector<edm::InputTag> >("inputs");
   dPhi_ = regPSet.getParameter<std::vector<double> >("deltaPhi");
   maxZ_ = regPSet.getParameter<std::vector<double> >("maxZ");
 
   tBeamSpot = iC.consumes<reco::BeamSpot>(beamSpotTag_);
   for (unsigned int t=0; t<inputs_.size(); t++ ) tCandidateView.push_back(iC.consumes< reco::CandidateView >(inputs_[t]));
 
   if (inputs_.size() != dPhi_.size() || dPhi_.size() != maxZ_.size() )
   {
     edm::LogError("PixelUnpackingRegions")<<"Not the same size of config parameters vectors!\n"
         <<"   inputs "<<inputs_.size()<<"  deltaPhi "<<dPhi_.size() <<"  maxZ "<< maxZ_.size();
   }
 }

PixelUnpackingRegions::~PixelUnpackingRegions ( )

inline

Definition at line 55 of file PixelUnpackingRegions.h.

55 {}

Member Function Documentation

void PixelUnpackingRegions::addRegion ( Region & r )

private

Definition at line 154 of file PixelUnpackingRegions.cc.

References addRegionLocal(), PixelUnpackingRegions::Region::dPhi, normalizedPhi(), nreg_, phi, phiBPIX_, phiFPIXm_, phiFPIXp_, and PixelUnpackingRegions::Region::v.

Referenced by run().

 {
   ++nreg_;
 
   float phi = normalizedPhi(r.v.phi());  // ensure [-pi,+pi]
 
   Module lo(phi - r.dPhi);
   Module hi(phi + r.dPhi);
 
   addRegionLocal(r, phiBPIX_, lo, hi);
   if (r.v.eta() >  1.)
   {
     addRegionLocal(r, phiFPIXp_, lo, hi);
   }
   if (r.v.eta() < -1.)
   {
     addRegionLocal(r, phiFPIXm_, lo, hi);
   }
 }

void PixelUnpackingRegions::addRegionLocal	(	Region &	r,
		std::vector< Module > &	container,
		const Module &	lo,
		const Module &	hi
	)

private

Definition at line 175 of file PixelUnpackingRegions.cc.

References a, b, gatherFromRange(), M_PI, and PixelUnpackingRegions::Module::phi.

Referenced by addRegion().

 {
   Module lo = _lo;
   Module hi = _hi;
   Module pi_m(-M_PI);
   Module pi_p( M_PI);
 
   std::vector<Module>::const_iterator a, b;
 
   if (lo.phi >= -M_PI && hi.phi <= M_PI) // interval doesn't cross the +-pi overlap
   {
     a = lower_bound(container.begin(), container.end(), lo);
     b = upper_bound(container.begin(), container.end(), hi);
     gatherFromRange(r, a, b);
   }
   else // interval is torn by the +-pi overlap
   {
     if (hi.phi >  M_PI) hi.phi -= 2.*M_PI;
     a = lower_bound(container.begin(), container.end(), pi_m);
     b = upper_bound(container.begin(), container.end(), hi);
     gatherFromRange(r, a, b);
 
     if (lo.phi < -M_PI) lo.phi += 2.*M_PI;
     a = lower_bound(container.begin(), container.end(), lo);
     b = upper_bound(container.begin(), container.end(), pi_p);
     gatherFromRange(r, a, b);
   }
 }

void PixelUnpackingRegions::gatherFromRange	(	Region &	r,
		std::vector< Module >::const_iterator	a,
		std::vector< Module >::const_iterator	b
	)

private

Definition at line 205 of file PixelUnpackingRegions.cc.

References a, funct::abs(), PixelUnpackingRegions::Region::atantheta, b, beamSpot_, PixelUnpackingRegions::Region::cosphi, feds_, PixelUnpackingRegions::Region::maxZ, modules_, and PixelUnpackingRegions::Region::sinphi.

Referenced by addRegionLocal().

 {
   for(; a != b; ++a)
   {
     // projection in r's direction onto beam's z
     float zmodule = a->z - (  (a->x - beamSpot_.x())*r.cosphi + (a->y - beamSpot_.y())*r.sinphi ) * r.atantheta;
 
     // do not include modules that project too far in z
     if ( std::abs(zmodule) > r.maxZ ) continue;
 
     feds_.insert(a->fed);
     modules_.insert(a->id);
   }
 }

void PixelUnpackingRegions::initialize ( const edm::EventSetup & es )

private

run by the run method: (re)initialize the cabling data when it's necessary

Definition at line 92 of file PixelUnpackingRegions.cc.

References assert(), cabling_, edm::ESWatcher< T >::check(), PixelUnpackingRegions::Module::fed, relativeConstraints::geom, edm::EventSetup::get(), PixelUnpackingRegions::Module::id, visualization-live-secondInstance_cfg::m, normalizedPhi(), PixelUnpackingRegions::Module::phi, phiBPIX_, phiFPIXm_, phiFPIXp_, PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, python.multivaluedict::sort(), watcherSiPixelFedCablingMap_, PixelUnpackingRegions::Module::x, PixelUnpackingRegions::Module::y, and PixelUnpackingRegions::Module::z.

Referenced by run().

 {
   // initialize cabling map or update it if necessary
   // and re-cache modules information
   if (watcherSiPixelFedCablingMap_.check( es ))
   {
     edm::ESTransientHandle<SiPixelFedCablingMap> cablingMap;
     es.get<SiPixelFedCablingMapRcd>().get( cablingMap );
     cabling_ = cablingMap->cablingTree();
 
     edm::ESHandle<TrackerGeometry> geom;
     // get the TrackerGeom
     es.get<TrackerDigiGeometryRecord>().get( geom );
 
     phiBPIX_.clear();
     phiFPIXp_.clear();
     phiFPIXm_.clear();
 
     phiBPIX_.reserve(1024);
     phiFPIXp_.reserve(512);
     phiFPIXm_.reserve(512);
 
     auto it = geom->dets().begin();
     for ( ; it != geom->dets().end(); ++it)
     {
       int subdet = (*it)->geographicalId().subdetId();
       if (! (subdet == PixelSubdetector::PixelBarrel ||
              subdet == PixelSubdetector::PixelEndcap) ) continue;
 
       Module m;
 
       m.x = (*it)->position().x();
       m.y = (*it)->position().y();
       m.z = (*it)->position().z();
 
       m.phi = normalizedPhi( (*it)->position().phi() ); // ensure [-pi,+pi]
 
       m.id = (*it)->geographicalId().rawId();
       const std::vector<sipixelobjects::CablingPathToDetUnit> path2det = cabling_->pathToDetUnit(m.id);
 
       m.fed = path2det[0].fed;
       assert(m.fed<40);
 
       if (subdet == PixelSubdetector::PixelBarrel)
       {
         phiBPIX_.push_back(m);
       }
       else if (subdet == PixelSubdetector::PixelEndcap)
       {
         if (m.z > 0.) phiFPIXp_.push_back(m);
         else phiFPIXm_.push_back(m);
       }
     }
 
     // pre-sort by phi
     std::sort(phiBPIX_.begin(),  phiBPIX_.end());
     std::sort(phiFPIXp_.begin(), phiFPIXp_.end());
     std::sort(phiFPIXm_.begin(), phiFPIXm_.end());
   }
 }

bool PixelUnpackingRegions::mayUnpackFED ( unsigned int fed_n ) const

check whether a FED has to be unpacked

Definition at line 221 of file PixelUnpackingRegions.cc.

References feds_.

Referenced by SiPixelRawToDigi::produce().

 {
   if (feds_.count(fed_n)) return true;
   return false;
 }

bool PixelUnpackingRegions::mayUnpackModule ( unsigned int id ) const

check whether a module has to be unpacked

Definition at line 238 of file PixelUnpackingRegions.cc.

References modules_.

 {
   if (modules_.count(id)) return true;
   return false;
 }

const std::set<unsigned int>* PixelUnpackingRegions::modulesToUnpack ( ) const

inline

full set of module ids to unpack

Definition at line 67 of file PixelUnpackingRegions.h.

References modules_.

Referenced by SiPixelRawToDigi::produce().

67 {return &modules_;}

PixelUnpackingRegions::modules_

std::set< unsigned int > modules_

Definition: PixelUnpackingRegions.h:108

unsigned int PixelUnpackingRegions::nBarrelFEDs ( ) const

Definition at line 227 of file PixelUnpackingRegions.cc.

References feds_.

Referenced by SiPixelRawToDigi::produce().

 {
   return std::count_if(feds_.begin(), feds_.end(), isBPIXFED );
 }

unsigned int PixelUnpackingRegions::nBarrelModules ( ) const

Definition at line 244 of file PixelUnpackingRegions.cc.

References modules_.

Referenced by SiPixelRawToDigi::produce().

 {
   return std::count_if(modules_.begin(), modules_.end(), isBPIXModule );
 }

unsigned int PixelUnpackingRegions::nFEDs ( ) const

inline

various informational accessors:

Definition at line 70 of file PixelUnpackingRegions.h.

References feds_.

Referenced by SiPixelRawToDigi::produce().

70 { return feds_.size(); }

PixelUnpackingRegions::feds_

std::set< unsigned int > feds_

Definition: PixelUnpackingRegions.h:107

unsigned int PixelUnpackingRegions::nForwardFEDs ( ) const

Definition at line 232 of file PixelUnpackingRegions.cc.

References feds_.

Referenced by SiPixelRawToDigi::produce().

 {
   return std::count_if(feds_.begin(), feds_.end(), isFPIXFED );
 }

unsigned int PixelUnpackingRegions::nForwardModules ( ) const

Definition at line 249 of file PixelUnpackingRegions.cc.

References modules_.

Referenced by SiPixelRawToDigi::produce().

 {
   return std::count_if(modules_.begin(), modules_.end(), isFPIXModule );
 }

unsigned int PixelUnpackingRegions::nModules ( ) const

inline

Definition at line 73 of file PixelUnpackingRegions.h.

References modules_.

Referenced by SiPixelRawToDigi::produce().

73 { return modules_.size(); }

PixelUnpackingRegions::modules_

std::set< unsigned int > modules_

Definition: PixelUnpackingRegions.h:108

unsigned int PixelUnpackingRegions::nRegions ( ) const

inline

Definition at line 76 of file PixelUnpackingRegions.h.

References nreg_.

76 { return nreg_; }

PixelUnpackingRegions::nreg_

unsigned int nreg_

Definition: PixelUnpackingRegions.h:109

void PixelUnpackingRegions::run	(	const edm::Event &	e,
		const edm::EventSetup &	es
	)

has to be run during each event

Definition at line 60 of file PixelUnpackingRegions.cc.

References addRegion(), SiPixelRawToDigiRegional_cfi::beamSpot, beamSpot_, EnergyCorrector::c, dPhi_, feds_, edm::Event::getByToken(), h, i, initialize(), input, inputs_, maxZ_, modules_, reco::Candidate::momentum(), gen::n, nreg_, alignCSCRings::r, tBeamSpot, and tCandidateView.

Referenced by Types.EventID::cppID(), Types.LuminosityBlockID::cppID(), and SiPixelRawToDigi::produce().

 {
   feds_.clear();
   modules_.clear();
   nreg_ = 0;
 
   initialize(es);
 
   edm::Handle<reco::BeamSpot> beamSpot;
   e.getByToken(tBeamSpot, beamSpot);
   beamSpot_ = beamSpot->position();
   //beamSpot_ = math::XYZPoint(0.,0.,0.);
 
   size_t ninputs = inputs_.size();
   for(size_t input = 0; input < ninputs; ++input)
   {
     edm::Handle< reco::CandidateView > h;
     e.getByToken(tCandidateView[input], h);
 
     size_t n = h->size();
     for(size_t i = 0; i < n; ++i )
     {
       const reco::Candidate & c = (*h)[i];
 
       // different input collections can have different dPhi and maxZ
       Region r(c.momentum(), dPhi_[input], maxZ_[input]);
       addRegion(r);
     }
   }
 }