#include <OnDemandMeasurementTracker.h>

Inheritance diagram for OnDemandMeasurementTracker:

Classes
class	DetODStatus

Public Types
typedef edm::LazyGetter < SiStripCluster >	LazyGetter

typedef edm::RefGetter < SiStripCluster >	RefGetter

Public Types inherited from MeasurementTrackerImpl
typedef std::unordered_map < unsigned int, MeasurementDet * >	DetContainer

enum	QualityFlags { BadModules =1, BadAPVFibers =2, BadStrips =4, MaskBad128StripBlocks =8, BadROCs =2 }

Public Types inherited from MeasurementTracker
enum	QualityFlags { BadModules =1, BadAPVFibers =2, BadStrips =4, MaskBad128StripBlocks =8, BadROCs =2 }

Public Member Functions
void	define (const edm::Handle< edm::LazyGetter< SiStripCluster > > &, RefGetter &, StMeasurementDetSet &) const
	OnDemandMeasurementTracker specific function to be called to define the region in the RefGetter according to MeasurementDet content. More...

const MeasurementDet *	idToDetBare (const DetId &id, const MeasurementTrackerEvent &data) const
	MeasurementTrackerImpl interface. More...

	OnDemandMeasurementTracker (const edm::ParameterSet &conf, const PixelClusterParameterEstimator pixelCPE, const StripClusterParameterEstimator stripCPE, const SiStripRecHitMatcher hitMatcher, const TrackerGeometry trackerGeom, const GeometricSearchTracker geometricSearchTracker, const SiStripQuality stripQuality, int stripQualityFlags, int stripQualityDebugFlags, const SiPixelQuality pixelQuality, const SiPixelFedCabling pixelCabling, int pixelQualityFlags, int pixelQualityDebugFlags, const SiStripRegionCabling *stripRegionCabling, bool isRegional=false)
	constructor More...

virtual	~OnDemandMeasurementTracker ()
	destructor More...

Public Member Functions inherited from MeasurementTrackerImpl
const DetContainer &	allDets () const
	For debug only. More...

const MeasurementDet *	findDet (const DetId &id) const

const GeometricSearchTracker *	geometricSearchTracker () const

const TrackingGeometry *	geomTracker () const

const std::vector < TkGluedMeasurementDet > &	gluedDets () const

virtual MeasurementDetWithData	idToDet (const DetId &id, const MeasurementTrackerEvent &data) const
	MeasurementDetSystem interface (won't be overloaded anymore) More...

	MeasurementTrackerImpl (const edm::ParameterSet &conf, const PixelClusterParameterEstimator pixelCPE, const StripClusterParameterEstimator stripCPE, const SiStripRecHitMatcher hitMatcher, const TrackerGeometry trackerGeom, const GeometricSearchTracker geometricSearchTracker, const SiStripQuality stripQuality, int stripQualityFlags, int stripQualityDebugFlags, const SiPixelQuality pixelQuality, const SiPixelFedCabling pixelCabling, int pixelQualityFlags, int pixelQualityDebugFlags, bool isRegional=false)

virtual const PxMeasurementConditionSet &	pixelDetConditions () const

const std::vector < TkPixelMeasurementDet > &	pixelDets () const

virtual const StMeasurementConditionSet &	stripDetConditions () const
	Provide templates to be filled in. More...

const std::vector < TkStripMeasurementDet > &	stripDets () const

virtual	~MeasurementTrackerImpl ()

Public Member Functions inherited from MeasurementTracker
const GeometricSearchTracker *	geometricSearchTracker () const

const TrackingGeometry *	geomTracker () const

	MeasurementTracker (TrackerGeometry const trackerGeom, GeometricSearchTracker const geometricSearchTracker)

virtual	~MeasurementTracker ()

Private Types
typedef std::unordered_map < unsigned int, DetODStatus >	DetODContainer

typedef std::vector< std::pair < SiStripRegionCabling::ElementIndex, std::vector < DetODContainer::const_iterator > > >	RegionalMap
	mapping of elementIndex -> iterator to the DetODMap: to know what are the regions that needs to be defined in the ref getter More...

Private Member Functions
void	assign (const TkStripMeasurementDet *csmdet, const MeasurementTrackerEvent &data) const
	assigne the cluster iterator to the TkStipMeasurementDet (const_cast in the way) More...

std::string	dumpCluster (const std::vector< SiStripCluster >::const_iterator &begin, const std::vector< SiStripCluster >::const_iterator &end) const
	some printouts, exclusively under LogDebug More...

std::string	dumpRegion (std::pair< unsigned int, unsigned int > indexes, const RefGetter &theGetter, bool stayUnpacked=false) const

Private Attributes
std::string	category_
	log category More...

bool	PixelOnDemand_
	internal flag to do pixel on demand (not configurable) false by default More...

RegionalMap	region_mapping

bool	StayPacked_
	internal flag to avoid unpacking things with LogDebug on More...

bool	StripOnDemand_
	internal flag to do strip on demand (not configurable) true by default More...

DetODContainer	theDetODMap
	mapping of detid -> MeasurementDet+flags+region_range More...

const SiStripRegionCabling *	theStripRegionCabling
	the cabling region tool to update a RefGetter More...

Additional Inherited Members
Protected Member Functions inherited from MeasurementTrackerImpl
void	addGluedDet (const GluedGeomDet *gd)

void	addPixelDet (const GeomDet *gd)

void	addPixelDets (const TrackingGeometry::DetContainer &dets)

void	addStripDet (const GeomDet *gd)

void	addStripDets (const TrackingGeometry::DetContainer &dets)

void	initGluedDet (TkGluedMeasurementDet &det)

void	initialize ()

void	initializePixelStatus (const SiPixelQuality stripQuality, const SiPixelFedCabling pixelCabling, int qualityFlags, int qualityDebugFlags)

void	initializeStripStatus (const SiStripQuality *stripQuality, int qualityFlags, int qualityDebugFlags)

void	initPxMeasurementConditionSet (std::vector< TkPixelMeasurementDet > &pixelDets)

void	initStMeasurementConditionSet (std::vector< TkStripMeasurementDet > &stripDets)

Protected Attributes inherited from MeasurementTrackerImpl
const std::string	name_

const edm::ParameterSet &	pset_

DetContainer	theDetMap

std::vector < TkGluedMeasurementDet >	theGluedDets

const SiPixelFedCabling *	thePixelCabling

std::vector < TkPixelMeasurementDet >	thePixelDets

PxMeasurementConditionSet	thePxDetConditions

StMeasurementConditionSet	theStDetConditions

std::vector < TkStripMeasurementDet >	theStripDets

Protected Attributes inherited from MeasurementTracker
const GeometricSearchTracker *	theGeometricSearchTracker

const TrackerGeometry *	theTrackerGeom

Detailed Description

Definition at line 13 of file OnDemandMeasurementTracker.h.

Member Typedef Documentation

typedef std::unordered_map<unsigned int, DetODStatus> OnDemandMeasurementTracker::DetODContainer

private

Definition at line 66 of file OnDemandMeasurementTracker.h.

typedef edm::LazyGetter<SiStripCluster> OnDemandMeasurementTracker::LazyGetter

Definition at line 34 of file OnDemandMeasurementTracker.h.

typedef edm::RefGetter<SiStripCluster> OnDemandMeasurementTracker::RefGetter

Definition at line 35 of file OnDemandMeasurementTracker.h.

typedef std::vector<std::pair<SiStripRegionCabling::ElementIndex, std::vector<DetODContainer::const_iterator> > > OnDemandMeasurementTracker::RegionalMap

private

mapping of elementIndex -> iterator to the DetODMap: to know what are the regions that needs to be defined in the ref getter

Definition at line 72 of file OnDemandMeasurementTracker.h.

Constructor & Destructor Documentation

OnDemandMeasurementTracker::OnDemandMeasurementTracker	(	const edm::ParameterSet &	conf,
		const PixelClusterParameterEstimator *	pixelCPE,
		const StripClusterParameterEstimator *	stripCPE,
		const SiStripRecHitMatcher *	hitMatcher,
		const TrackerGeometry *	trackerGeom,
		const GeometricSearchTracker *	geometricSearchTracker,
		const SiStripQuality *	stripQuality,
		int	stripQualityFlags,
		int	stripQualityDebugFlags,
		const SiPixelQuality *	pixelQuality,
		const SiPixelFedCabling *	pixelCabling,
		int	pixelQualityFlags,
		int	pixelQualityDebugFlags,
		const SiStripRegionCabling *	stripRegionCabling,
		bool	isRegional = `false`
	)

constructor

Definition at line 51 of file OnDemandMeasurementTracker.cc.

References category_, SiStripRegionCabling::elementIndex(), PV3DBase< T, PVType, FrameType >::eta(), eta(), first, OnDemandMeasurementTracker::DetODStatus::Glued, SiStripRegionCabling::layerFromDetId(), LogDebug, phi, PV3DBase< T, PVType, FrameType >::phi(), OnDemandMeasurementTracker::DetODStatus::Pixel, GeomDetEnumerators::PixelBarrel, GeomDetEnumerators::PixelEndcap, SiStripRegionCabling::position(), SiStripRegionCabling::positionIndex(), SiStripRegionCabling::region(), region_mapping, edm::second(), OnDemandMeasurementTracker::DetODStatus::Strip, SiStripRegionCabling::subdetFromDetId(), GeomDetEnumerators::TEC, MeasurementTrackerImpl::theDetMap, theDetODMap, theStripRegionCabling, GeomDetEnumerators::TIB, GeomDetEnumerators::TID, and GeomDetEnumerators::TOB.

                                                :
   MeasurementTrackerImpl(conf,pixelCPE,stripCPE,hitMatcher,trackerGeom,geometricSearchTracker,
         stripQuality,stripQualityFlags,stripQualityDebugFlags,
         pixelQuality,pixelCabling,pixelQualityFlags,pixelQualityDebugFlags,
         isRegional)
   , category_("OnDemandMeasurementTracker")
   , StayPacked_(true)
   , StripOnDemand_(true)
   , PixelOnDemand_(false)
   , theStripRegionCabling(stripRegionCabling)
 {
   //  the constructor does construct the regular MeasurementTracker
   //  then a smart copy of the DetMap is made into DetODMap: this could be avoided with modification to MeasurementDet interface
   //  the elementIndex to be defined in the refgetter is mapped to the detId
   //  flags are set to initialize the DetODMap
   
   std::map<SiStripRegionCabling::ElementIndex, std::vector< DetODContainer::const_iterator> > local_mapping;
   for (DetContainer::iterator it=theDetMap.begin(); it!= theDetMap.end();++it)
     {
       DetODContainer::iterator inserted = theDetODMap.insert(make_pair(it->first,DetODStatus(it->second))).first;
 
 
       GeomDet::SubDetector subdet = it->second->geomDet().subDetector();
       if (subdet == GeomDetEnumerators::PixelBarrel  || subdet == GeomDetEnumerators::PixelEndcap ){
         inserted->second.index = -1;
         inserted->second.kind  = DetODStatus::Pixel;
       }//pixel module
       else if (subdet == GeomDetEnumerators::TIB || subdet == GeomDetEnumerators::TOB ||
       subdet == GeomDetEnumerators::TID || subdet == GeomDetEnumerators::TEC )
     {
       //set flag to false
           if (typeid(*it->second) == typeid(TkStripMeasurementDet)) {
             const TkStripMeasurementDet & sdet = static_cast<const TkStripMeasurementDet &>(*it->second);
             inserted->second.index = sdet.index();
             inserted->second.kind = DetODStatus::Strip;
           } else {
             inserted->second.kind = DetODStatus::Glued;
           }
 
       //what will be the element index in the refgetter
       GlobalPoint center = it->second->geomDet().position();
       double eta = center.eta();
       double phi = center.phi();
       uint32_t id = it->first;
       SiStripRegionCabling::ElementIndex eIndex = theStripRegionCabling->elementIndex(SiStripRegionCabling::Position(eta,phi),
                                               SiStripRegionCabling::subdetFromDetId(id),
                                               SiStripRegionCabling::layerFromDetId(id));
       LogDebug(category_)<<"region selected (from "<<id<<" center) is:\n"
                  <<"position: "<<center
                  <<"\n center absolute index: "<<theStripRegionCabling->region(theStripRegionCabling->positionIndex(SiStripRegionCabling::Position(eta,phi)))
                  <<"\n center position index: "<<theStripRegionCabling->positionIndex(SiStripRegionCabling::Position(eta,phi)).first<<
         " "<<theStripRegionCabling->positionIndex(SiStripRegionCabling::Position(eta,phi)).second
                  <<"\n center postion: "<<theStripRegionCabling->position(theStripRegionCabling->positionIndex(SiStripRegionCabling::Position(eta,phi))).first<<
         " "<<theStripRegionCabling->position(theStripRegionCabling->positionIndex(SiStripRegionCabling::Position(eta,phi))).second
                  <<"\n eta: "<<eta
                  <<"\n phi: "<<phi
                  <<"\n subedet: "<<SiStripRegionCabling::subdetFromDetId(id)
                  <<" layer: "<<SiStripRegionCabling::layerFromDetId(id);
 
       //      register those in a map
       //to be able to know what are the detid in a given elementIndex
       local_mapping[eIndex].push_back(inserted);
     }//strip module
       else{
     //abort
     edm::LogError(category_)<<"not a tracker geomdet in constructor: "<<it->first;
     throw MeasurementDetException("OnDemandMeasurementTracker dealing with a non tracker GeomDet.");
       }//abort
     }//loop over DetMap
   //move into a vector
   region_mapping.reserve(local_mapping.size());
   for( auto eIt= local_mapping.begin();
        eIt!=local_mapping.end();++eIt)
     region_mapping.push_back(std::make_pair((*eIt).first,(*eIt).second));
 }

virtual OnDemandMeasurementTracker::~OnDemandMeasurementTracker ( )

inlinevirtual

destructor

Definition at line 32 of file OnDemandMeasurementTracker.h.

32 {}

Member Function Documentation

void OnDemandMeasurementTracker::assign	(	const TkStripMeasurementDet *	csmdet,
		const MeasurementTrackerEvent &	data
	)		const

private

assigne the cluster iterator to the TkStipMeasurementDet (const_cast in the way)

Definition at line 219 of file OnDemandMeasurementTracker.cc.

References category_, dumpCluster(), dumpRegion(), edm::RefGetter< T >::find(), LogDebug, DetId::rawId(), StMeasurementDetSet::rawInactiveStripDetIds(), StMeasurementDetSet::refGetter(), StMeasurementDetSet::regionRange(), StMeasurementDetSet::setUpdated(), and MeasurementTrackerEvent::stripData().

Referenced by idToDetBare().

                                                                                    {
   //  assign is using the handle to the refgetter and the region index range to update the MeasurementDet with their clusters
  
   StMeasurementDetSet & rwdata = const_cast<StMeasurementDetSet &>(data.stripData());
 
   DetId id = smdet->rawId();
   
   LogDebug(category_)<<"assigning: "<<id.rawId();
 
     rwdata.setUpdated(smdet->index());
 
     if (!data.stripData().rawInactiveStripDetIds().empty() && std::binary_search(data.stripData().rawInactiveStripDetIds().begin(), data.stripData().rawInactiveStripDetIds().end(), id)) {
       smdet->setActiveThisEvent(rwdata, false); 
       return;
     }
 
     //retrieve the region range index for this module
     const std::pair<unsigned int,unsigned int> & indexes = data.stripData().regionRange(smdet->index());
 
     //this printout will trigger unpacking. no problem. it is done on the next regular line (find(id.rawId())
     LogDebug(category_)<<"between index: "<<indexes.first<<" and: "<<indexes.second
                <<"\nretrieved for module: "<<id.rawId()
                <<"\n"<<dumpRegion(indexes,data.stripData().refGetter());
     
     //look for iterator range in the regions defined for that module
     for (unsigned int iRegion = indexes.first; iRegion != indexes.second; ++iRegion){
       RefGetter::record_pair range = data.stripData().refGetter()[iRegion].find(id.rawId());
       if (range.first!=range.second){
     //  found something not empty
     //update the measurementDet
     smdet->update(rwdata, range.first, range.second);
     LogDebug(category_)<<"Valid clusters for: "<<id.rawId()
                <<"\nnumber of regions defined here: "<< indexes.second-indexes.first
                <<"\n"<<dumpCluster(range.first,range.second);
     /* since theStripsToSkip is a "static" pointer of the MT, no need to set it at all time.
       if (selfUpdateSkipClusters_){
       //assign skip clusters
       smdet->setClusterToSkip(&theStripsToSkip);
     }
     */
     //and you are done
     return;}
     }//loop over regions, between indexes
 
     //if reached. no cluster are found. set the TkStripMeasurementDet to be empty
     smdet->setEmpty(rwdata);
 
 }

void OnDemandMeasurementTracker::define	(	const edm::Handle< edm::LazyGetter< SiStripCluster > > &	aLazyGetterH,
		RefGetter &	aGetter,
		StMeasurementDetSet &	stData
	)		const

OnDemandMeasurementTracker specific function to be called to define the region in the RefGetter according to MeasurementDet content.

Definition at line 141 of file OnDemandMeasurementTracker.cc.

References category_, StMeasurementDetSet::defineStrip(), dumpRegion(), HTMLExport::elem(), OnDemandMeasurementTracker::DetODStatus::index, OnDemandMeasurementTracker::DetODStatus::kind, LogDebug, region_mapping, StMeasurementDetSet::resetOnDemandStrips(), edm::RefGetter< T >::size(), StayPacked_, OnDemandMeasurementTracker::DetODStatus::Strip, theStripRegionCabling, and SiStripRegionCabling::updateSiStripRefGetter().

Referenced by MeasurementTrackerSiStripRefGetterProducer::produce().

 {
   //  define is supposed to be call by an EDProducer module, which wil put the RefGetter in the event
   //  so that reference can be made to it.
   //  the lazy getter is retrieved by the calling module and passed along with the event
   //  the map is cleared, except for pixel
   //  then the known elementIndex are defined to the RefGetter. no unpacking is done at this time
   //  the defined region range is registered in the DetODMap for further use.
   stData.resetOnDemandStrips();
  
   //define all the elementindex in the refgetter
   for(auto eIt= region_mapping.begin();
        eIt!=region_mapping.end();++eIt){
     std::pair<unsigned int, unsigned int> region_range; 
     
     //before update of the refgetter
     region_range.first = aGetter.size();
     //update the refegetter with the elementindex
     theStripRegionCabling->updateSiStripRefGetter<SiStripCluster> (aGetter, aLazyGetterH, eIt->first);
     //after update of the refgetter
     region_range.second = aGetter.size();
 
     LogDebug(category_)<<"between index: "<<region_range.first<<" "<<region_range.second
                <<"\n"<<dumpRegion(region_range,aGetter,StayPacked_);
     
     //now assign to each measurement det for that element index
     for (auto dIt=eIt->second.begin();
      dIt!=eIt->second.end();++dIt){
       const DetODStatus & elem = (*dIt)->second;
       if (elem.kind == DetODStatus::Strip) {
           stData.defineStrip(elem.index, region_range);
       }
       LogDebug(category_)<<"detId: "<<(*dIt)->first<<" in region range: "<<region_range.first<<" "<<region_range.second;
     }//loop over MeasurementDet attached to that elementIndex
   }//loop over know elementindex
 }

std::string OnDemandMeasurementTracker::dumpCluster	(	const std::vector< SiStripCluster >::const_iterator &	begin,
		const std::vector< SiStripCluster >::const_iterator &	end
	)		const

private

some printouts, exclusively under LogDebug

Definition at line 181 of file OnDemandMeasurementTracker.cc.

References begin, end, i, AlCaHLTBitMon_QueryRunRegistry::string, and accesses::tab.

Referenced by assign(), and dumpRegion().

 {
   //  dumpCluster is a printout of all the clusters between the iterator. returns a string
   std::string tab="      ";
   std::stringstream ss;
   std::vector<SiStripCluster> ::const_iterator it = begin;
   unsigned int i=0;
   for (;it!=end;++it){
     ss<<tab<<i++<<") center: "<<it->barycenter()<<",id: "<<it->geographicalId()<<" with: "<<it->amplitudes().size()<<" strips\n"<<tab<<tab<<"{";
     for (unsigned int is=0;is!=it->amplitudes().size();++is){
       ss<<it->amplitudes()[is]<<" ";
     }ss<<"}\n";
   }
   return ss.str();
 }

std::string OnDemandMeasurementTracker::dumpRegion	(	std::pair< unsigned int, unsigned int >	indexes,
		const RefGetter &	theGetter,
		bool	stayUnpacked = `false`
	)		const

private

Definition at line 197 of file OnDemandMeasurementTracker.cc.

References begin, dumpCluster(), end, SiStripRegionCabling::position(), SiStripRegionCabling::positionIndex(), SiStripRegionCabling::region(), and theStripRegionCabling.

Referenced by assign(), and define().

 {
   //  dumpRegion is a printout of all the clusters in a region defined on the RefGetter. returns a string
   std::stringstream ss;
   ss<<"cluster between: "<<indexes.first<<" and: "<<indexes.second<<"\n";
   for (unsigned int iRegion = indexes.first; iRegion != indexes.second; ++iRegion){    
     uint32_t reg = SiStripRegionCabling::region((theGetter)[iRegion].region());
     SiStripRegionCabling::Position pos = theStripRegionCabling->position(reg);
     SiStripRegionCabling::PositionIndex posI = theStripRegionCabling->positionIndex(reg);
     
     ss<<"Clusters for region:["<<iRegion<<"]"
       <<"\n element index: "<<(theGetter)[iRegion].region()
       <<"\n region absolute index: "<<reg
       <<"\n region position index: "<<posI.first<<" "<<posI.second
       <<"\n region position: "<<pos.first<<" "<<pos.second
       <<"\n"<< (stayPacked? " hidden to avoid unpacking." : dumpCluster((theGetter)[iRegion].begin(),(theGetter)[iRegion].end()));
   }
   return ss.str();
 }

const MeasurementDet * OnDemandMeasurementTracker::idToDetBare	(	const DetId &	id,
		const MeasurementTrackerEvent &	data
	)		const

virtual

MeasurementTrackerImpl interface.

Reimplemented from MeasurementTrackerImpl.

Definition at line 274 of file OnDemandMeasurementTracker.cc.

References assign(), category_, data, OnDemandMeasurementTracker::DetODStatus::Glued, OnDemandMeasurementTracker::DetODStatus::Pixel, OnDemandMeasurementTracker::DetODStatus::Strip, MeasurementTrackerEvent::stripData(), StMeasurementDetSet::stripDefined(), StMeasurementDetSet::stripUpdated(), and theDetODMap.

 {
   //  overloaded from MeasurementTracker
   //  find the detid. if not found throw exception
   //  if already updated: always for pixel or strip already queried. return it
   DetODContainer::const_iterator it = theDetODMap.find(id);
   if ( it != theDetODMap.end()) {
     switch (it->second.kind) {
         case DetODStatus::Pixel:
                 // nothing to do
             break;
         case DetODStatus::Strip: {
                 const TkStripMeasurementDet*  theConcreteDet = static_cast<const TkStripMeasurementDet*>(it->second.mdet);
                 assert(data.stripData().stripDefined(theConcreteDet->index()));
                 if (!data.stripData().stripUpdated(theConcreteDet->index())) assign(theConcreteDet, data);
             }
             break;
         case DetODStatus::Glued: {
                 const TkGluedMeasurementDet*  theConcreteDet = static_cast<const TkGluedMeasurementDet*>(it->second.mdet);
                 int imono = theConcreteDet->monoDet()->index(), istereo = theConcreteDet->stereoDet()->index();
                 assert(data.stripData().stripDefined(imono) && data.stripData().stripDefined(istereo));
                 if (!data.stripData().stripUpdated(imono))   assign(theConcreteDet->monoDet(), data);
                 if (!data.stripData().stripUpdated(istereo)) assign(theConcreteDet->stereoDet(), data);
             }
             break;
         }
         return it->second.mdet;
   } 
   else{
     //throw excpetion
     edm::LogError(category_)<<"failed to find the MeasurementDet for: "<<id.rawId();
     throw MeasurementDetException("failed to find the MeasurementDet for: <see message logger>");
   }
   return 0;
 }