#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	nBarrelModules () const

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

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 35 of file PixelUnpackingRegions.cc.

References beamSpotTag_, dPhi_, edm::ParameterSet::getParameter(), inputs_, maxZ_, lumiQTWidget::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.

References MillePedeFileConverter_cfg::e, mayUnpackFED(), mayUnpackModule(), and run().

55 {}

Member Function Documentation

void PixelUnpackingRegions::addRegion ( Region & r )

private

Definition at line 159 of file PixelUnpackingRegions.cc.

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

Referenced by run().

 {
   ++nreg_;
 
   float phi = r.v.phi();
 
   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 180 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 210 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 87 of file PixelUnpackingRegions.cc.

References cabling_, SiPixelFedCablingMap::cablingTree(), edm::ESWatcher< T >::check(), TrackerGeometry::dets(), PixelUnpackingRegions::Module::fed, relativeConstraints::geom, edm::EventSetup::get(), PixelUnpackingRegions::Module::id, TrackerGeometry::isThere(), funct::m, FEDNumbering::MAXSiPixelFEDID, FEDNumbering::MAXSiPixeluTCAFEDID, FEDNumbering::MINSiPixelFEDID, FEDNumbering::MINSiPixeluTCAFEDID, GeomDetEnumerators::P1PXB, GeomDetEnumerators::P1PXEC, PixelUnpackingRegions::Module::phi, phiBPIX_, phiFPIXm_, phiFPIXp_, PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, 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 );
 
     // switch on the phase1 
     unsigned int fedMin = FEDNumbering::MINSiPixelFEDID; // phase0
     unsigned int fedMax = FEDNumbering::MAXSiPixelFEDID;
     if( (geom->isThere(GeomDetEnumerators::P1PXB)) && 
     (geom->isThere(GeomDetEnumerators::P1PXEC)) ) {
       fedMin = FEDNumbering::MINSiPixeluTCAFEDID; // phase1
       fedMax = FEDNumbering::MAXSiPixeluTCAFEDID;
     }
 
 
     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 = (*it)->position().phi();
 
       m.id = (*it)->geographicalId().rawId();
       const std::vector<sipixelobjects::CablingPathToDetUnit> path2det = cabling_->pathToDetUnit(m.id);
 
       m.fed = path2det[0].fed;
       assert( (m.fed<=fedMax) && (m.fed>=fedMin) );
 
       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 226 of file PixelUnpackingRegions.cc.

References feds_.

Referenced by SiPixelRawToDigi::produce(), and ~PixelUnpackingRegions().

 {
   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 233 of file PixelUnpackingRegions.cc.

References modules_.

Referenced by ~PixelUnpackingRegions().

 {
   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:106

unsigned int PixelUnpackingRegions::nBarrelModules ( ) const

Definition at line 239 of file PixelUnpackingRegions.cc.

References modules_.

Referenced by nModules(), and 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:105

unsigned int PixelUnpackingRegions::nForwardModules ( ) const

Definition at line 244 of file PixelUnpackingRegions.cc.

References modules_.

Referenced by nModules(), and SiPixelRawToDigi::produce().

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

unsigned int PixelUnpackingRegions::nModules ( ) const

inline

Definition at line 71 of file PixelUnpackingRegions.h.

References modules_, nBarrelModules(), and nForwardModules().

Referenced by SiPixelRawToDigi::produce().

71 { return modules_.size(); }

PixelUnpackingRegions::modules_

std::set< unsigned int > modules_

Definition: PixelUnpackingRegions.h:106

unsigned int PixelUnpackingRegions::nRegions ( ) const

inline

Definition at line 74 of file PixelUnpackingRegions.h.

References nreg_.

74 { return nreg_; }

PixelUnpackingRegions::nreg_

unsigned int nreg_

Definition: PixelUnpackingRegions.h:107

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

has to be run during each event

Definition at line 55 of file PixelUnpackingRegions.cc.

References addRegion(), ecalDrivenElectronSeedsParameters_cff::beamSpot, beamSpot_, EnergyCorrector::c, dPhi_, feds_, edm::Event::getByToken(), h, mps_fire::i, initialize(), input, inputs_, maxZ_, modules_, reco::Candidate::momentum(), gen::n, nreg_, reco::BeamSpot::position(), alignCSCRings::r, edm::View< T >::size(), tBeamSpot, and tCandidateView.

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

 {
   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);
     }
   }
 }