#include <MeasurementTrackerImpl.h>

Inheritance diagram for MeasurementTrackerImpl:

Public Types
typedef std::map< DetId, MeasurementDet * >	DetContainer
enum	QualityFlags { BadModules = 1, BadAPVFibers = 2, BadStrips = 4, MaskBad128StripBlocks = 8, BadROCs = 2 }
Public Member Functions
const DetContainer &	allDets () const
	For debug only.
const GeometricSearchTracker *	geometricSearchTracker () const
const TrackingGeometry *	geomTracker () const
const std::vector < TkGluedMeasurementDet * > &	gluedDets () const
virtual const MeasurementDet *	idToDet (const DetId &id) const
	MeasurementDetSystem interface.
	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)
const std::vector < TkPixelMeasurementDet * > &	pixelDets () const
void	setClusterToSkip (const edm::InputTag &cluster, const edm::Event &event) const
const std::vector < TkStripMeasurementDet * > &	stripDets () const
void	unsetClusterToSkip () const
virtual void	update (const edm::Event &) const
void	updatePixels (const edm::Event &) const
void	updateStrips (const edm::Event &) const
virtual	~MeasurementTrackerImpl ()
Protected Member Functions
void	addGluedDet (const GluedGeomDet gd, const SiStripRecHitMatcher matcher) const
void	addPixelDet (const GeomDet gd, const PixelClusterParameterEstimator cpe) const
void	addPixelDets (const TrackingGeometry::DetContainer &dets) const
void	addStripDet (const GeomDet gd, const StripClusterParameterEstimator cpe) const
void	addStripDets (const TrackingGeometry::DetContainer &dets) const
void	getInactiveStrips (const edm::Event &event, std::vector< uint32_t > &rawInactiveDetIds) const
void	initialize () const
void	initializePixelStatus (const SiPixelQuality stripQuality, const SiPixelFedCabling pixelCabling, int qualityFlags, int qualityDebugFlags) const
void	initializeStripStatus (const SiStripQuality *stripQuality, int qualityFlags, int qualityDebugFlags) const
Protected Attributes
bool	isRegional_
const std::string	name_
const edm::ParameterSet &	pset_
bool	selfUpdateSkipClusters_
DetContainer	theDetMap
std::vector < TkGluedMeasurementDet * >	theGluedDets
const SiStripRecHitMatcher *	theHitMatcher
const std::vector< edm::InputTag >	theInactivePixelDetectorLabels
const std::vector< edm::InputTag >	theInactiveStripDetectorLabels
const SiPixelFedCabling *	thePixelCabling
const PixelClusterParameterEstimator *	thePixelCPE
std::vector < TkPixelMeasurementDet * >	thePixelDets
const StripClusterParameterEstimator *	theStripCPE
std::vector < TkStripMeasurementDet * >	theStripDets

Detailed Description

Definition at line 40 of file MeasurementTrackerImpl.h.

Member Typedef Documentation

typedef std::map<DetId,MeasurementDet*> MeasurementTrackerImpl::DetContainer

Definition at line 74 of file MeasurementTrackerImpl.h.

Member Enumeration Documentation

enum MeasurementTrackerImpl::QualityFlags

Enumerator:

BadModules
BadAPVFibers
BadStrips
MaskBad128StripBlocks
BadROCs

Reimplemented from MeasurementTracker.

Definition at line 42 of file MeasurementTrackerImpl.h.

                     { BadModules=1, // for everybody
                       /* Strips: */ BadAPVFibers=2, BadStrips=4, MaskBad128StripBlocks=8, 
                       /* Pixels: */ BadROCs=2 };

Constructor & Destructor Documentation

MeasurementTrackerImpl::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`
	)

Definition at line 86 of file MeasurementTrackerImpl.cc.

References edm::ParameterSet::exists(), edm::ParameterSet::getParameter(), initialize(), initializePixelStatus(), initializeStripStatus(), LogDebug, selfUpdateSkipClusters_, and createPayload::skip.

                                                        :
  MeasurementTracker(trackerGeom,geometricSearchTracker),
  pset_(conf),
  name_(conf.getParameter<std::string>("ComponentName")),
  thePixelCPE(pixelCPE),theStripCPE(stripCPE),theHitMatcher(hitMatcher),
  theInactivePixelDetectorLabels(conf.getParameter<std::vector<edm::InputTag> >("inactivePixelDetectorLabels")),
  theInactiveStripDetectorLabels(conf.getParameter<std::vector<edm::InputTag> >("inactiveStripDetectorLabels")),
  isRegional_(isRegional)
{
  this->initialize();
  this->initializeStripStatus(stripQuality, stripQualityFlags, stripQualityDebugFlags);
  this->initializePixelStatus(pixelQuality, pixelCabling, pixelQualityFlags, pixelQualityDebugFlags);
  //the measurement tracking is set to skip clusters, the other option is set from outside
  selfUpdateSkipClusters_=conf.exists("skipClusters");
  if (selfUpdateSkipClusters_)
    {
             edm::InputTag skip=conf.getParameter<edm::InputTag>("skipClusters");
             if (skip==edm::InputTag("")) selfUpdateSkipClusters_=false;
    }


  LogDebug("MeasurementTracker")<<"skipping clusters: "<<selfUpdateSkipClusters_;
}

MeasurementTrackerImpl::~MeasurementTrackerImpl ( ) [virtual]

Definition at line 123 of file MeasurementTrackerImpl.cc.

References theGluedDets, thePixelDets, and theStripDets.

{
  for(vector<TkPixelMeasurementDet*>::const_iterator it=thePixelDets.begin(); it!=thePixelDets.end(); ++it){
    delete *it;
  }

  for(vector<TkStripMeasurementDet*>::const_iterator it=theStripDets.begin(); it!=theStripDets.end(); ++it){
    delete *it;
  }

  for(vector<TkGluedMeasurementDet*>::const_iterator it=theGluedDets.begin(); it!=theGluedDets.end(); ++it){
    delete *it;
  }
}

Member Function Documentation

void MeasurementTrackerImpl::addGluedDet	(	const GluedGeomDet *	gd,
		const SiStripRecHitMatcher *	matcher
	)		const `[protected]`

Definition at line 209 of file MeasurementTrackerImpl.cc.

References addStripDet(), GeomDet::geographicalId(), idToDet(), GluedGeomDet::monoDet(), GluedGeomDet::stereoDet(), theDetMap, theGluedDets, theHitMatcher, and theStripCPE.

Referenced by addStripDets().

{
  const MeasurementDet* monoDet = idToDet( gd->monoDet()->geographicalId());
  if (monoDet == 0) {
    addStripDet(gd->monoDet(), theStripCPE);  // in case glued det comes before components
    monoDet = idToDet( gd->monoDet()->geographicalId());
  }

  const MeasurementDet* stereoDet = idToDet( gd->stereoDet()->geographicalId());
  if (stereoDet == 0) {
    addStripDet(gd->stereoDet(), theStripCPE);  // in case glued det comes before components
    stereoDet = idToDet( gd->stereoDet()->geographicalId());
  }

  if (monoDet == 0 || stereoDet == 0) {
    edm::LogError("MeasurementDet") << "MeasurementTracker ERROR: GluedDet components not found as MeasurementDets ";
    throw MeasurementDetException("MeasurementTracker ERROR: GluedDet components not found as MeasurementDets");
  }

  TkGluedMeasurementDet* det = new TkGluedMeasurementDet( gd, theHitMatcher,
                                                          monoDet, stereoDet);
  theGluedDets.push_back( det);
  theDetMap[gd->geographicalId()] = det;
}

void MeasurementTrackerImpl::addPixelDet	(	const GeomDet *	gd,
		const PixelClusterParameterEstimator *	cpe
	)		const `[protected]`

Definition at line 201 of file MeasurementTrackerImpl.cc.

References GeomDet::geographicalId(), theDetMap, and thePixelDets.

Referenced by addPixelDets().

{
  TkPixelMeasurementDet* det = new TkPixelMeasurementDet( gd, cpe);
  thePixelDets.push_back(det);
  theDetMap[gd->geographicalId()] = det;
}

void MeasurementTrackerImpl::addPixelDets ( const TrackingGeometry::DetContainer & dets ) const [protected]

Definition at line 155 of file MeasurementTrackerImpl.cc.

References addPixelDet(), and thePixelCPE.

Referenced by initialize().

{
  for (TrackerGeometry::DetContainer::const_iterator gd=dets.begin();
       gd != dets.end(); gd++) {
    addPixelDet(*gd, thePixelCPE);
  }  
}

void MeasurementTrackerImpl::addStripDet	(	const GeomDet *	gd,
		const StripClusterParameterEstimator *	cpe
	)		const `[protected]`

Definition at line 188 of file MeasurementTrackerImpl.cc.

References GeomDet::geographicalId(), isRegional_, theDetMap, theStripDets, and cms::Exception::what().

Referenced by addGluedDet(), and addStripDets().

{
  try {
    TkStripMeasurementDet* det = new TkStripMeasurementDet( gd, cpe,isRegional_);
    theStripDets.push_back(det);
    theDetMap[gd->geographicalId()] = det;
  }
  catch(MeasurementDetException& err){
    edm::LogError("MeasurementDet") << "Oops, got a MeasurementDetException: " << err.what() ;
  }
}

void MeasurementTrackerImpl::addStripDets ( const TrackingGeometry::DetContainer & dets ) const [protected]

Definition at line 163 of file MeasurementTrackerImpl.cc.

References addGluedDet(), addStripDet(), theHitMatcher, and theStripCPE.

Referenced by initialize().

{
  for (TrackerGeometry::DetContainer::const_iterator gd=dets.begin();
       gd != dets.end(); gd++) {

    const GeomDetUnit* gdu = dynamic_cast<const GeomDetUnit*>(*gd);

    //    StripSubdetector stripId( (**gd).geographicalId());
    //     bool isDetUnit( gdu != 0);
    //     cout << "StripSubdetector glued? " << stripId.glued() 
    //   << " is DetUnit? " << isDetUnit << endl;

    if (gdu != 0) {
      addStripDet(*gd, theStripCPE);
    }
    else {
      const GluedGeomDet* gluedDet = dynamic_cast<const GluedGeomDet*>(*gd);
      if (gluedDet == 0) {
        throw MeasurementDetException("MeasurementTracker ERROR: GeomDet neither DetUnit nor GluedDet");
      }
      addGluedDet(gluedDet, theHitMatcher);
    }  
  }
}

const DetContainer& MeasurementTrackerImpl::allDets ( ) const [inline]

For debug only.

Definition at line 77 of file MeasurementTrackerImpl.h.

References theDetMap.

{return theDetMap;}

const GeometricSearchTracker* MeasurementTrackerImpl::geometricSearchTracker ( ) const [inline]

Reimplemented from MeasurementTracker.

Definition at line 69 of file MeasurementTrackerImpl.h.

References MeasurementTracker::theGeometricSearchTracker.

{return theGeometricSearchTracker;}

const TrackingGeometry* MeasurementTrackerImpl::geomTracker ( ) const [inline]

Reimplemented from MeasurementTracker.

Definition at line 67 of file MeasurementTrackerImpl.h.

References MeasurementTracker::theTrackerGeom.

{ return theTrackerGeom;}

void MeasurementTrackerImpl::getInactiveStrips	(	const edm::Event &	event,
		std::vector< uint32_t > &	rawInactiveDetIds
	)		const `[protected]`

Definition at line 390 of file MeasurementTrackerImpl.cc.

References python::multivaluedict::sort(), and theInactiveStripDetectorLabels.

Referenced by OnDemandMeasurementTracker::updateStrips(), and updateStrips().

                                                                                              {
  if (!theInactiveStripDetectorLabels.empty()) {
    edm::Handle<DetIdCollection> detIds;
    for (std::vector<edm::InputTag>::const_iterator itt = theInactiveStripDetectorLabels.begin(), edt = theInactiveStripDetectorLabels.end(); 
         itt != edt; ++itt) {
      event.getByLabel(*itt, detIds);
      rawInactiveDetIds.insert(rawInactiveDetIds.end(), detIds->begin(), detIds->end());
    }
    if (!rawInactiveDetIds.empty()) std::sort(rawInactiveDetIds.begin(), rawInactiveDetIds.end());
  }
}

const std::vector<TkGluedMeasurementDet*>& MeasurementTrackerImpl::gluedDets ( ) const [inline]

Definition at line 80 of file MeasurementTrackerImpl.h.

References theGluedDets.

{return theGluedDets;}

const MeasurementDet * MeasurementTrackerImpl::idToDet ( const DetId & id ) const [virtual]

MeasurementDetSystem interface.

Implements MeasurementTracker.

Reimplemented in OnDemandMeasurementTracker.

Definition at line 577 of file MeasurementTrackerImpl.cc.

References theDetMap.

Referenced by addGluedDet(), and updateStrips().

{
  std::map<DetId,MeasurementDet*>::const_iterator it = theDetMap.find(id);
  if(it !=theDetMap.end()) {
    return it->second;
  }else{
    //throw exception;
  }
  
  return 0; //to avoid compile warning
}

void MeasurementTrackerImpl::initialize ( ) const [protected]

Definition at line 139 of file MeasurementTrackerImpl.cc.

References addPixelDets(), addStripDets(), TrackerGeometry::detsPXB(), TrackerGeometry::detsPXF(), TrackerGeometry::detsTEC(), TrackerGeometry::detsTIB(), TrackerGeometry::detsTID(), TrackerGeometry::detsTOB(), theGluedDets, thePixelDets, theStripDets, and MeasurementTracker::theTrackerGeom.

Referenced by MeasurementTrackerImpl().

{  
  addPixelDets( theTrackerGeom->detsPXB());
  addPixelDets( theTrackerGeom->detsPXF());
  addStripDets( theTrackerGeom->detsTIB());
  addStripDets( theTrackerGeom->detsTID());
  addStripDets( theTrackerGeom->detsTOB());
  addStripDets( theTrackerGeom->detsTEC());  

  sortTKD(theStripDets);
  sortTKD(thePixelDets);
  sortTKD(theGluedDets);

}

void MeasurementTrackerImpl::initializePixelStatus	(	const SiPixelQuality *	stripQuality,
		const SiPixelFedCabling *	pixelCabling,
		int	qualityFlags,
		int	qualityDebugFlags
	)		const `[protected]`

Definition at line 667 of file MeasurementTrackerImpl.cc.

References BadModules, BadROCs, cond::rpcobgas::detid, SiPixelQuality::getBadRocPositions(), SiPixelQuality::getBadRocs(), i, SiPixelQuality::IsModuleUsable(), thePixelDets, and MeasurementTracker::theTrackerGeom.

Referenced by MeasurementTrackerImpl().

                                                                                                                                                                      {
  if ((quality != 0) && (qualityFlags != 0))  {
    edm::LogInfo("MeasurementTracker") << "qualityFlags = " << qualityFlags;
    unsigned int on = 0, tot = 0, badrocs = 0; 
    for (std::vector<TkPixelMeasurementDet*>::const_iterator i=thePixelDets.begin();
         i!=thePixelDets.end(); i++) {
      uint32_t detid = ((**i).geomDet().geographicalId()).rawId();
      if (qualityFlags & BadModules) {
          bool isOn = quality->IsModuleUsable(detid);
          (*i)->setActive(isOn);
          tot++; on += (unsigned int) isOn;
          if (qualityDebugFlags & BadModules) {
            edm::LogInfo("MeasurementTracker")<< "MeasurementTrackerImpl::initializePixelStatus : detid " << detid << " is " << (isOn ?  "on" : "off");
          }
       } else {
          (*i)->setActive(true);
       }
       if ((qualityFlags & BadROCs) && (quality->getBadRocs(detid) != 0)) {
          std::vector<LocalPoint> badROCs = quality->getBadRocPositions(detid, *theTrackerGeom, pixelCabling);
          badrocs += badROCs.size();
          (*i)->setBadRocPositions(badROCs);
       } else {
          (*i)->clearBadRocPositions();  
       }
    }
    if (qualityDebugFlags & BadModules) {
        edm::LogInfo("MeasurementTracker PixelModuleStatus") << 
            " Total modules: " << tot << ", active " << on <<", inactive " << (tot - on);
    }
    if (qualityDebugFlags & BadROCs) {
        edm::LogInfo("MeasurementTracker PixelROCStatus") << " Total of bad ROCs: " << badrocs ;
    }
  } else {
    for (std::vector<TkPixelMeasurementDet*>::const_iterator i=thePixelDets.begin();
         i!=thePixelDets.end(); i++) {
      (*i)->setActive(true);          // module ON
    }
  }
}

void MeasurementTrackerImpl::initializeStripStatus	(	const SiStripQuality *	stripQuality,
		int	qualityFlags,
		int	qualityDebugFlags
	)		const `[protected]`

Definition at line 589 of file MeasurementTrackerImpl.cc.

References BadAPVFibers, BadModules, BadStrips, SiStripBadStrip::decode(), cond::rpcobgas::detid, SiStripQuality::getBadApvs(), SiStripQuality::getBadFibers(), edm::ParameterSet::getParameter(), SiStripBadStrip::getRange(), i, SiStripQuality::IsModuleUsable(), j, MaskBad128StripBlocks, pset_, SiStripDetId::TEC, theStripDets, SiStripDetId::TIB, SiStripDetId::TID, and SiStripDetId::TOB.

Referenced by MeasurementTrackerImpl().

                                                                                                                               {
  TkStripMeasurementDet::BadStripCuts badStripCuts[4];
  if (qualityFlags & BadStrips) {
     edm::ParameterSet cutPset = pset_.getParameter<edm::ParameterSet>("badStripCuts");
     badStripCuts[SiStripDetId::TIB-3] = TkStripMeasurementDet::BadStripCuts(cutPset.getParameter<edm::ParameterSet>("TIB"));
     badStripCuts[SiStripDetId::TOB-3] = TkStripMeasurementDet::BadStripCuts(cutPset.getParameter<edm::ParameterSet>("TOB"));
     badStripCuts[SiStripDetId::TID-3] = TkStripMeasurementDet::BadStripCuts(cutPset.getParameter<edm::ParameterSet>("TID"));
     badStripCuts[SiStripDetId::TEC-3] = TkStripMeasurementDet::BadStripCuts(cutPset.getParameter<edm::ParameterSet>("TEC"));
  }

  if ((quality != 0) && (qualityFlags != 0))  {
    edm::LogInfo("MeasurementTracker") << "qualityFlags = " << qualityFlags;
    unsigned int on = 0, tot = 0; 
    unsigned int foff = 0, ftot = 0, aoff = 0, atot = 0; 
    for (std::vector<TkStripMeasurementDet*>::const_iterator i=theStripDets.begin();
         i!=theStripDets.end(); i++) {
      uint32_t detid = ((**i).geomDet().geographicalId()).rawId();
      if (qualityFlags & BadModules) {
          bool isOn = quality->IsModuleUsable(detid);
          (*i)->setActive(isOn);
          tot++; on += (unsigned int) isOn;
          if (qualityDebugFlags & BadModules) {
            edm::LogInfo("MeasurementTracker")<< "MeasurementTrackerImpl::initializeStripStatus : detid " << detid << " is " << (isOn ?  "on" : "off");
          }
       } else {
          (*i)->setActive(true);
       }
       // first turn all APVs and fibers ON
       (*i)->set128StripStatus(true); 
       if (qualityFlags & BadAPVFibers) {
          short badApvs   = quality->getBadApvs(detid);
          short badFibers = quality->getBadFibers(detid);
          for (int j = 0; j < 6; j++) {
             atot++;
             if (badApvs & (1 << j)) {
                (*i)->set128StripStatus(false, j);
                aoff++;
             }
          }
          for (int j = 0; j < 3; j++) {
             ftot++;
             if (badFibers & (1 << j)) {
                (*i)->set128StripStatus(false, 2*j);
                (*i)->set128StripStatus(false, 2*j+1);
                foff++;
             }
          }
          (*i)->setMaskBad128StripBlocks((qualityFlags & MaskBad128StripBlocks) != 0);
       } 
       std::vector<TkStripMeasurementDet::BadStripBlock> &badStrips = (*i)->getBadStripBlocks();
       badStrips.clear();
       if (qualityFlags & BadStrips) {
            SiStripBadStrip::Range range = quality->getRange(detid);
            for (SiStripBadStrip::ContainerIterator bit = range.first; bit != range.second; ++bit) {
                badStrips.push_back(quality->decode(*bit));
            }
            (*i)->setBadStripCuts(badStripCuts[SiStripDetId(detid).subdetId()-3]);
       }
    }
    if (qualityDebugFlags & BadModules) {
        edm::LogInfo("MeasurementTracker StripModuleStatus") << 
            " Total modules: " << tot << ", active " << on <<", inactive " << (tot - on);
    }
    if (qualityDebugFlags & BadAPVFibers) {
        edm::LogInfo("MeasurementTracker StripAPVStatus") << 
            " Total APVs: " << atot << ", active " << (atot-aoff) <<", inactive " << (aoff);
        edm::LogInfo("MeasurementTracker StripFiberStatus") << 
            " Total Fibers: " << ftot << ", active " << (ftot-foff) <<", inactive " << (foff);
    }
  } else {
    for (std::vector<TkStripMeasurementDet*>::const_iterator i=theStripDets.begin();
         i!=theStripDets.end(); i++) {
      (*i)->setActive(true);          // module ON
      (*i)->set128StripStatus(true);  // all APVs and fibers ON
    }
  }
}

const std::vector<TkPixelMeasurementDet*>& MeasurementTrackerImpl::pixelDets ( ) const [inline]

Definition at line 79 of file MeasurementTrackerImpl.h.

References thePixelDets.

{return thePixelDets;}

void MeasurementTrackerImpl::setClusterToSkip	(	const edm::InputTag &	cluster,
		const edm::Event &	event
	)		const `[virtual]`

Implements MeasurementTracker.

Definition at line 255 of file MeasurementTrackerImpl.cc.

References edmNew::DetSetVector< T >::begin(), edmNew::DetSetVector< T >::end(), f, i, LogDebug, selfUpdateSkipClusters_, thePixelDets, and theStripDets.

{
  LogDebug("MeasurementTracker")<<"setClusterToSkip";
  if (selfUpdateSkipClusters_)
    edm::LogError("MeasurementTracker")<<"this mode of operation is not supported, either the measurement tracker is set to skip clusters, or is being told to skip clusters. not both";

  edm::Handle<edmNew::DetSetVector<SiPixelClusterRefNew> > pixelClusterRefs;
  event.getByLabel(cluster,pixelClusterRefs);
  for (std::vector<TkPixelMeasurementDet*>::const_iterator i=thePixelDets.begin();
       i!=thePixelDets.end(); i++) {
    edmNew::DetSetVector<SiPixelClusterRefNew>::const_iterator f=pixelClusterRefs->find((**i).geomDet().geographicalId().rawId());
    if (f!=pixelClusterRefs->end())
      (**i).setClusterToSkip(f->begin(),f->end());
    else
      (**i).unset();
  }

  edm::Handle<edmNew::DetSetVector<TkStripMeasurementDet::SiStripClusterRef> > stripClusterRefs;
  event.getByLabel(cluster,stripClusterRefs);
  for (std::vector<TkStripMeasurementDet*>::const_iterator i=theStripDets.begin();
       i!=theStripDets.end(); i++){
    edmNew::DetSetVector<TkStripMeasurementDet::SiStripClusterRef>::const_iterator f=stripClusterRefs->find((**i).geomDet().geographicalId().rawId());
    if (f!=stripClusterRefs->end())
      (**i).setClusterToSkip(f->begin(),f->end());
    else
      (**i).unset();
  }
}

const std::vector<TkStripMeasurementDet*>& MeasurementTrackerImpl::stripDets ( ) const [inline]

Definition at line 78 of file MeasurementTrackerImpl.h.

References theStripDets.

{return theStripDets;}

void MeasurementTrackerImpl::unsetClusterToSkip ( ) const [virtual]

Implements MeasurementTracker.

Definition at line 284 of file MeasurementTrackerImpl.cc.

References i, LogDebug, selfUpdateSkipClusters_, thePixelDets, and theStripDets.

                                                      {
  LogDebug("MeasurementTracker")<<"unsetClusterToSkip";
  if (selfUpdateSkipClusters_)
    edm::LogError("MeasurementTracker")<<"this mode of operation is not supported, either the measurement tracker is set to skip clusters, or is being told to skip clusters. not both";

  for (std::vector<TkPixelMeasurementDet*>::const_iterator i=thePixelDets.begin();
       i!=thePixelDets.end(); i++) {
    (**i).unset();
  }
  for (std::vector<TkStripMeasurementDet*>::const_iterator i=theStripDets.begin();
       i!=theStripDets.end(); i++){
    (**i).unset();
  }
}

void MeasurementTrackerImpl::update ( const edm::Event & event ) const [virtual]

Implements MeasurementTracker.

Reimplemented in OnDemandMeasurementTracker.

Definition at line 235 of file MeasurementTrackerImpl.cc.

References updatePixels(), and updateStrips().

{
  updatePixels(event);
  updateStrips(event);
  
  /*
  for (std::vector<TkStripMeasurementDet*>::const_iterator i=theStripDets.begin();
       i!=theStripDets.end(); i++) {
    if( (*i)->isEmpty()){
      std::cout << "stripDet id, #hits: " 
                <<  (*i)->geomDet().geographicalId().rawId() << " , "
                << 0 << std::endl;
    }else{
      std::cout << "stripDet id, #hits: " 
                <<  (*i)->geomDet().geographicalId().rawId() << " , "
                << (*i)->size() << std::endl;
    }
  }
  */
}

void MeasurementTrackerImpl::updatePixels ( const edm::Event & event ) const [virtual]

Implements MeasurementTracker.

Definition at line 299 of file MeasurementTrackerImpl.cc.

References edmNew::DetSetVector< T >::begin(), edmNew::DetSetVector< T >::empty(), edmNew::DetSetVector< T >::end(), edm::ParameterSet::existsAs(), f, edm::HandleBase::failedToGet(), edmNew::DetSetVector< T >::find(), edm::Event::getByLabel(), edm::ParameterSet::getParameter(), i, LogDebug, name_, edm::Handle< T >::product(), pset_, selfUpdateSkipClusters_, python::multivaluedict::sort(), theInactivePixelDetectorLabels, and thePixelDets.

Referenced by update(), and OnDemandMeasurementTracker::update().

{
  // avoid to update twice from the same event
  if (!edm::Service<UpdaterService>()->checkOnce("MeasurementTrackerImpl::updatePixels::"+name_)) return;

  typedef edmNew::DetSet<SiPixelCluster> PixelDetSet;

  bool switchOffPixelsIfEmpty = (!pset_.existsAs<bool>("switchOffPixelsIfEmpty")) ||
                                (pset_.getParameter<bool>("switchOffPixelsIfEmpty"));

  std::vector<uint32_t> rawInactiveDetIds; 
  if (!theInactivePixelDetectorLabels.empty()) {
    edm::Handle<DetIdCollection> detIds;
    for (std::vector<edm::InputTag>::const_iterator itt = theInactivePixelDetectorLabels.begin(), edt = theInactivePixelDetectorLabels.end(); 
            itt != edt; ++itt) {
      if (event.getByLabel(*itt, detIds)){
        rawInactiveDetIds.insert(rawInactiveDetIds.end(), detIds->begin(), detIds->end());
      }else{
        static bool iFailedAlready=false;
        if (!iFailedAlready){
          edm::LogError("MissingProduct")<<"I fail to get the list of inactive pixel modules, because of 4.2/4.4 event content change.";
          iFailedAlready=true;
        }
      }
    }
    if (!rawInactiveDetIds.empty()) std::sort(rawInactiveDetIds.begin(), rawInactiveDetIds.end());
  }
  // Pixel Clusters
  std::string pixelClusterProducer = pset_.getParameter<std::string>("pixelClusterProducer");
  if( pixelClusterProducer.empty() ) { //clusters have not been produced
    for (std::vector<TkPixelMeasurementDet*>::const_iterator i=thePixelDets.begin();
         i!=thePixelDets.end(); i++) {
      if (switchOffPixelsIfEmpty) {
        (**i).setActiveThisEvent(false);
      }else{
        (**i).setEmpty();
      }
    }
  }else{  

    edm::Handle<edmNew::DetSetVector<SiPixelClusterRefNew> > pixelClusterRefs;
    if (selfUpdateSkipClusters_){
      //and get the collection of pixel ref to skip
      event.getByLabel(pset_.getParameter<edm::InputTag>("skipClusters"),pixelClusterRefs);
      LogDebug("MeasurementTracker")<<"getting pxl refs to skip";
      if (pixelClusterRefs.failedToGet())edm::LogError("MeasurementTracker")<<"not getting the pixel clusters to skip";
    }

    edm::Handle<edmNew::DetSetVector<SiPixelCluster> > pixelClusters;
    event.getByLabel(pixelClusterProducer, pixelClusters);
    const  edmNew::DetSetVector<SiPixelCluster>* pixelCollection = pixelClusters.product();
   
    if (switchOffPixelsIfEmpty && pixelCollection->empty()) {
        for (std::vector<TkPixelMeasurementDet*>::const_iterator i=thePixelDets.begin();
             i!=thePixelDets.end(); i++) {
              (**i).setActiveThisEvent(false);
        }
    } else { 

    for (std::vector<TkPixelMeasurementDet*>::const_iterator i=thePixelDets.begin();
         i!=thePixelDets.end(); i++) {

      // foreach det get cluster range
      unsigned int id = (**i).geomDet().geographicalId().rawId();
      if (!rawInactiveDetIds.empty() && std::binary_search(rawInactiveDetIds.begin(), rawInactiveDetIds.end(), id)) {
        (**i).setActiveThisEvent(false); continue;
      }
      //FIXME
      //fill the set with what needs to be skipped
      edmNew::DetSetVector<SiPixelCluster>::const_iterator it = pixelCollection->find( id );
      if ( it != pixelCollection->end() ){            
        // push cluster range in det
        (**i).update( *it, pixelClusters, id );
        if (selfUpdateSkipClusters_){
          edmNew::DetSetVector<SiPixelClusterRefNew>::const_iterator f=pixelClusterRefs->find(id);
          if (f!=pixelClusterRefs->end())
            (**i).setClusterToSkip(f->begin(),f->end());
          else
            (**i).unset();
        }
      }else{
        (**i).setEmpty();
      }
    }
    }
  }


  
}

void MeasurementTrackerImpl::updateStrips ( const edm::Event & event ) const [virtual]

Implements MeasurementTracker.

Reimplemented in OnDemandMeasurementTracker.

Definition at line 403 of file MeasurementTrackerImpl.cc.

References edmNew::DetSetVector< T >::begin(), edm::RegionIndex< T >::begin(), edmNew::DetSetVector< T >::end(), edm::RegionIndex< T >::end(), end, f, edm::HandleBase::failedToGet(), getInactiveStrips(), edm::ParameterSet::getParameter(), i, edmNew::DetSetVector< T >::id(), idToDet(), isRegional_, LogDebug, name_, edm::Handle< T >::product(), pset_, selfUpdateSkipClusters_, TkStripMeasurementDet::setRegionalClustersToSkip(), theStripDets, TkStripMeasurementDet::unset(), and TkStripMeasurementDet::update().

Referenced by update().

{
  // avoid to update twice from the same event
  if (!edm::Service<UpdaterService>()->checkOnce("MeasurementTrackerImpl::updateStrips::"+name_)) return;

  typedef edmNew::DetSet<SiStripCluster>   StripDetSet;

  std::vector<uint32_t> rawInactiveDetIds;
  getInactiveStrips(event,rawInactiveDetIds);

  // Strip Clusters
  std::string stripClusterProducer = pset_.getParameter<std::string>("stripClusterProducer");
  //first clear all of them

  {
    std::vector<TkStripMeasurementDet*>::const_iterator end = theStripDets.end()-200;
    for (std::vector<TkStripMeasurementDet*>::const_iterator i=theStripDets.begin();
         i!=end; i++) {
      (**i).setEmpty();
#ifdef __SSE2__
      _mm_prefetch(((char *)(*(i+200))),_MM_HINT_T0); 
#endif
    }
   for (std::vector<TkStripMeasurementDet*>::const_iterator i=end;
         i!=theStripDets.end(); i++)
      (**i).setEmpty();
  }
  if( !stripClusterProducer.compare("") ) { //clusters have not been produced
  }else{
    //=========  actually load cluster =============
    if(!isRegional_){
        edm::Handle<edmNew::DetSetVector<TkStripMeasurementDet::SiStripClusterRef> > stripClusterRefs;
        if (selfUpdateSkipClusters_){
          //and get the collection of pixel ref to skip
          LogDebug("MeasurementTracker")<<"getting strp refs to skip";
          event.getByLabel(pset_.getParameter<edm::InputTag>("skipClusters"),stripClusterRefs);
          if (stripClusterRefs.failedToGet())edm::LogError("MeasurementTracker")<<"not getting the strip clusters to skip";
        }

      edm::Handle<edmNew::DetSetVector<SiStripCluster> > clusterHandle;
      event.getByLabel(stripClusterProducer, clusterHandle);
      const edmNew::DetSetVector<SiStripCluster>* clusterCollection = clusterHandle.product();

      std::vector<TkStripMeasurementDet*>::const_iterator i=theStripDets.begin();
      std::vector<TkStripMeasurementDet*>::const_iterator endDet=theStripDets.end();
      edmNew::DetSetVector<SiStripCluster>::const_iterator it = (*clusterCollection).begin();
      edmNew::DetSetVector<SiStripCluster>::const_iterator endColl = (*clusterCollection).end();
      // cluster and det and in order (both) and unique so let's use set intersection
      for (;it!=endColl; ++it) {
        StripDetSet detSet = *it;
        unsigned int id = detSet.id();
        while ( id != (**i).rawId()) { // eventually change to lower_range
          ++i;
          if (i==endDet) throw "we have a problem!!!!";
        }

        if (!rawInactiveDetIds.empty() && std::binary_search(rawInactiveDetIds.begin(), rawInactiveDetIds.end(), id)) {
            (**i).setActiveThisEvent(false); continue;
        }
        // push cluster range in det

        (**i).update( detSet, clusterHandle, id );
        if (selfUpdateSkipClusters_){
          edmNew::DetSetVector<TkStripMeasurementDet::SiStripClusterRef>::const_iterator f=stripClusterRefs->find(id);
          if (f!=stripClusterRefs->end())
            (**i).setClusterToSkip(f->begin(),f->end());
          else
            (**i).unset();
        }
      }
    }else{

      edm::Handle<edmNew::DetSetVector<TkStripMeasurementDet::SiStripRegionalClusterRef> > stripClusterRefs;
      if(selfUpdateSkipClusters_){
        LogDebug("MeasurementTracker")<<"getting reg strp refs to skip";
        event.getByLabel(pset_.getParameter<edm::InputTag>("skipClusters"),stripClusterRefs);
        if (stripClusterRefs.failedToGet())edm::LogError("MeasurementTracker")<<"not getting the strip clusters to skip";
      }

      //then set the not-empty ones only
      edm::Handle<edm::RefGetter<SiStripCluster> > refClusterHandle;
      event.getByLabel(stripClusterProducer, refClusterHandle);
      
      std::string lazyGetter = pset_.getParameter<std::string>("stripLazyGetterProducer");
      edm::Handle<edm::LazyGetter<SiStripCluster> > lazyClusterHandle;
      event.getByLabel(lazyGetter,lazyClusterHandle);

      uint32_t tmpId=0;
      vector<SiStripCluster>::const_iterator beginIterator;
      edm::RefGetter<SiStripCluster>::const_iterator iregion = refClusterHandle->begin();
      for(;iregion!=refClusterHandle->end();++iregion) {
        const edm::RegionIndex<SiStripCluster>& region = *iregion;
        vector<SiStripCluster>::const_iterator icluster = region.begin();
        const vector<SiStripCluster>::const_iterator endIterator = region.end();
        tmpId = icluster->geographicalId();
        beginIterator = icluster;

        //std::cout << "== tmpId ad inizio loop dentro region: " << tmpId << std::endl;

        for (;icluster!=endIterator;icluster++) {
          //std::cout << "===== cluster id,pos " 
          //  << icluster->geographicalId() << " , " << icluster->barycenter()
          //  << std::endl;
          //std::cout << "=====making ref in recHits() " << std::endl;
          if( icluster->geographicalId() != tmpId){ 
            //std::cout << "geo!=tmpId" << std::endl;
            //we should find a way to avoid this casting. it is slow
            //create also another map for TkStripMeasurementDet ??

            // the following castings are really ugly. To be corrected ASAP
            const TkStripMeasurementDet* theConcreteDet = 
              dynamic_cast<const TkStripMeasurementDet*>(idToDet(DetId(tmpId)));
            
            if(theConcreteDet == 0)
              throw MeasurementDetException("failed casting to TkStripMeasurementDet*");            
            TkStripMeasurementDet*  theConcreteDetUpdatable = 
              const_cast<TkStripMeasurementDet*>(theConcreteDet);
            theConcreteDetUpdatable->update(beginIterator,icluster,lazyClusterHandle,tmpId);
            if (selfUpdateSkipClusters_){
              edmNew::DetSetVector<TkStripMeasurementDet::SiStripRegionalClusterRef>::const_iterator f=stripClusterRefs->find(tmpId);
              if (f!=stripClusterRefs->end())
                theConcreteDetUpdatable->setRegionalClustersToSkip(f->begin(),f->end());
              else
                theConcreteDetUpdatable->unset();
            }
            //cannot we avoid to update the det with detId of itself??

            tmpId = icluster->geographicalId();
            beginIterator = icluster;
            if( icluster == (endIterator-1)){
              const TkStripMeasurementDet* theConcreteDet = 
              dynamic_cast<const TkStripMeasurementDet*>(idToDet(DetId(tmpId)));
              
              if(theConcreteDet == 0)
              throw MeasurementDetException("failed casting to TkStripMeasurementDet*");            
              TkStripMeasurementDet*  theConcreteDetUpdatable = 
              const_cast<TkStripMeasurementDet*>(theConcreteDet);
              theConcreteDetUpdatable->update(icluster,endIterator,lazyClusterHandle,tmpId);
              if (selfUpdateSkipClusters_){
                edmNew::DetSetVector<TkStripMeasurementDet::SiStripRegionalClusterRef>::const_iterator f=stripClusterRefs->find(tmpId);
                if (f!=stripClusterRefs->end())
                  theConcreteDetUpdatable->setRegionalClustersToSkip(f->begin(),f->end());
                else
                  theConcreteDetUpdatable->unset();
              }
            }   
          }else if( icluster == (endIterator-1)){          
            const TkStripMeasurementDet* theConcreteDet = 
              dynamic_cast<const TkStripMeasurementDet*>(idToDet(DetId(tmpId)));
            
            if(theConcreteDet == 0)
              throw MeasurementDetException("failed casting to TkStripMeasurementDet*");            
            TkStripMeasurementDet*  theConcreteDetUpdatable = 
              const_cast<TkStripMeasurementDet*>(theConcreteDet);
            //std::cout << "=== option3. fill det with id,#clust: " << tmpId  << " , " 
            //      << iregion->end() - beginIterator << std::endl;
            theConcreteDetUpdatable->update(beginIterator,endIterator,lazyClusterHandle,tmpId);  
            if (selfUpdateSkipClusters_){
              edmNew::DetSetVector<TkStripMeasurementDet::SiStripRegionalClusterRef>::const_iterator f=stripClusterRefs->find(tmpId);
              if (f!=stripClusterRefs->end())
                theConcreteDetUpdatable->setRegionalClustersToSkip(f->begin(),f->end());
              else
                theConcreteDetUpdatable->unset();
            }    
          }
        }//end loop cluster in one ragion
      }
    }//end of block for updating with regional clusters 
  }
}