ValidHitPairFilter Class Reference

#include <ValidHitPairFilter.h>

Inheritance diagram for ValidHitPairFilter:

Public Member Functions
virtual bool	operator() (const reco::Track *track, const std::vector< const TrackingRecHit * > &recHits, const TrackerTopology *tTopo) const
void	update (const edm::Event &ev, const edm::EventSetup &es) override
	ValidHitPairFilter (const edm::ParameterSet &ps, edm::ConsumesCollector &iC)
virtual	~ValidHitPairFilter ()
Public Member Functions inherited from PixelTrackFilter
virtual bool	operator() (const reco::Track *) const
virtual bool	operator() (const reco::Track *, const Hits &) const
virtual bool	operator() (const reco::Track *, const Hits &, const TrackerTopology *tTopo) const
virtual	~PixelTrackFilter ()

Private Member Functions
std::vector< const GeomDet * >	getCloseDets (int il, float rz, const std::vector< float > &rzB, float ph, const std::vector< float > &phB, const TrackerTopology *tTopo) const
int	getLayer (const TrackingRecHit &recHit, const TrackerTopology *tTopo) const
std::vector< int >	getMissingLayers (int a, int b) const
FreeTrajectoryState	getTrajectory (const reco::Track &track) const

Private Attributes
const std::vector< DetLayer * >	detLayers
std::vector< float >	phBounds [7]
std::vector< float >	rzBounds [7]
const GeometricSearchTracker *	theGSTracker
const MagneticField *	theMagneticField
const Propagator *	thePropagator
const TrackerGeometry *	theTracker

Additional Inherited Members
Public Types inherited from PixelTrackFilter
typedef std::vector< const TrackingRecHit * >	Hits

Detailed Description

Definition at line 20 of file ValidHitPairFilter.h.

Constructor & Destructor Documentation

ValidHitPairFilter::ValidHitPairFilter	(	const edm::ParameterSet &	ps,
		edm::ConsumesCollector &	iC
	)

Definition at line 52 of file ValidHitPairFilter.cc.

{
}

ValidHitPairFilter::~ValidHitPairFilter ( )

virtual

Definition at line 57 of file ValidHitPairFilter.cc.

{
}

Member Function Documentation

vector< const GeomDet * > ValidHitPairFilter::getCloseDets ( int  il, float  rz, const std::vector< float > &  rzB, float  ph, const std::vector< float > &  phB, const TrackerTopology *  tTopo  ) const

private

Definition at line 212 of file ValidHitPairFilter.cc.

References BPix3, LogTrace, TrackerTopology::pxbDetId(), and TrackerTopology::pxfDetId().

{
  vector<int> rzVec, phVec;

  // Radius or z
  vector<float>::const_iterator irz = lower_bound(rzB.begin(),rzB.end(), rz);
  if(irz > rzB.begin()) rzVec.push_back(irz - rzB.begin() - 1);
  if(irz < rzB.end()  ) rzVec.push_back(irz - rzB.begin()    );

  // Phi, circular
  vector<float>::const_iterator iph = lower_bound(phB.begin(),phB.end(), ph);
  if(iph > phB.begin()) phVec.push_back(iph         - phB.begin() - 1);
                   else phVec.push_back(phB.end()   - phB.begin() - 1);
  if(iph < phB.end()  ) phVec.push_back(iph         - phB.begin()    );
                   else phVec.push_back(phB.begin() - phB.begin()    );

  // Detectors
  vector<const GeomDet *> dets;

  for(vector<int>::const_iterator irz = rzVec.begin(); irz!= rzVec.end(); irz++)
  for(vector<int>::const_iterator iph = phVec.begin(); iph!= phVec.end(); iph++)
  {
    if(il+1 <= BPix3)
    {
      int layer  =  il  + 1;
      int ladder = *iph + 1;
      int module = *irz + 1;

      LogTrace("MinBiasTracking")
    << "  [ValidHitPairFilter]  third ("<<layer<< "|"<<ladder<<"|"<<module<<")";

      DetId id=tTopo->pxbDetId(layer,ladder,module);
      dets.push_back(theTracker->idToDet(id));
    }
    else
    {
      int side = (il - BPix3) / 2 +1;  
      int disk   = (il - BPix3) % 2 + 1;
      int blade  =  *iph + 1;
      int panel  = (*irz) % 2 + 1;
      int module = (*irz) / 2 + 1;

      LogTrace("MinBiasTracking")
    << "  [ValidHitPairFilter]  third ("<<side<<"|"<<disk<<"|"<<blade<<"|"<<panel<<"|"<<module<<")";

     
      DetId id=tTopo->pxfDetId(side,disk,blade,panel,module);
      dets.push_back(theTracker->idToDet(id));
      
    }
  }

  return dets;
}

int ValidHitPairFilter::getLayer ( const TrackingRecHit &  recHit, const TrackerTopology *  tTopo  ) const

private

Definition at line 141 of file ValidHitPairFilter.cc.

References BPix3, TrackingRecHit::geographicalId(), PixelSubdetector::PixelBarrel, TrackerTopology::pxbLayer(), TrackerTopology::pxfDisk(), and TrackerTopology::pxfSide().

{
  DetId id(recHit.geographicalId());

  if(id.subdetId() == int(PixelSubdetector::PixelBarrel))
  {
    
    return tTopo->pxbLayer(id);
  }
  else
  {
    
    return BPix3 + ((tTopo->pxfSide(id)-1) << 1) + tTopo->pxfDisk(id);
  }
}

vector< int > ValidHitPairFilter::getMissingLayers ( int  a, int  b  ) const

private

Definition at line 158 of file ValidHitPairFilter.cc.

References BPix1, BPix2, BPix3, EnergyCorrector::c, FPix1_neg, FPix1_pos, FPix2_neg, FPix2_pos, and prof2calltree::l.

{
  vector<int> l;
  pair<int,int> c(a,b);

  if(c == pair<int,int>(BPix1,BPix2))     { l.push_back(int(BPix3));
                                            l.push_back(int(FPix1_pos));
                                            l.push_back(int(FPix1_neg)); return l; }
  if(c == pair<int,int>(BPix1,BPix3))     { l.push_back(int(BPix2));     return l; }
  if(c == pair<int,int>(BPix2,BPix3))     { l.push_back(int(BPix1));     return l; }
  if(c == pair<int,int>(BPix1,FPix1_pos)) { l.push_back(int(BPix2));
                                            l.push_back(int(FPix2_pos)); return l; }
  if(c == pair<int,int>(BPix1,FPix1_neg)) { l.push_back(int(BPix2)); 
                                            l.push_back(int(FPix2_neg)); return l; }
  if(c == pair<int,int>(BPix1,FPix2_pos)) { l.push_back(int(BPix2)); 
                                            l.push_back(int(FPix1_pos)); return l; }
  if(c == pair<int,int>(BPix1,FPix2_neg)) { l.push_back(int(BPix2)); 
                                            l.push_back(int(FPix1_neg)); return l; }
  if(c == pair<int,int>(BPix2,FPix1_pos)) { l.push_back(int(BPix1)); 
                                            l.push_back(int(FPix2_pos)); return l; }
  if(c == pair<int,int>(BPix2,FPix1_neg)) { l.push_back(int(BPix1)); 
                                            l.push_back(int(FPix2_neg)); return l; }
  if(c == pair<int,int>(BPix2,FPix2_pos)) { l.push_back(int(BPix1));
                                            l.push_back(int(FPix1_pos)); return l; }
  if(c == pair<int,int>(BPix2,FPix2_neg)) { l.push_back(int(BPix1)); 
                                            l.push_back(int(FPix1_neg)); return l; }
  if(c == pair<int,int>(FPix1_pos,FPix2_pos)) { l.push_back(int(BPix1)); return l; }
  if(c == pair<int,int>(FPix1_neg,FPix2_neg)) { l.push_back(int(BPix1)); return l; }

  return l;
}

FreeTrajectoryState ValidHitPairFilter::getTrajectory ( const reco::Track &  track ) const

private

Definition at line 192 of file ValidHitPairFilter.cc.

References reco::TrackBase::charge(), reco::TrackBase::momentum(), position, and reco::TrackBase::vertex().

{
  GlobalPoint position(track.vertex().x(),
                       track.vertex().y(),
                       track.vertex().z());

  GlobalVector momentum(track.momentum().x(),
                        track.momentum().y(),
                        track.momentum().z());
    
  GlobalTrajectoryParameters gtp(position, momentum,
                                 track.charge(), theMagneticField);
  
  FreeTrajectoryState fts(gtp);

  return fts;
}

bool ValidHitPairFilter::operator() ( const reco::Track *  track, const std::vector< const TrackingRecHit * > &  recHits, const TrackerTopology *  tTopo  ) const

virtual

Definition at line 272 of file ValidHitPairFilter.cc.

References BPix3, HitInfo::getInfo(), TrajectoryStateOnSurface::globalPosition(), TrajectoryStateOnSurface::isValid(), TrajectoryStateOnSurface::localPosition(), LogTrace, PV3DBase< T, PVType, FrameType >::perp(), PV3DBase< T, PVType, FrameType >::phi(), phi(), HLT_25ns14e33_v3_cff::recHits, spin(), ForwardDetLayer::surface(), BarrelDetLayer::surface(), and PV3DBase< T, PVType, FrameType >::z().

{
  bool hasGap = true;

  if(recHits.size() == 2)
  {
    LogTrace("MinBiasTracking")
      << "  [ValidHitPairFilter] pair" << HitInfo::getInfo(*(recHits[0]),tTopo)
      << HitInfo::getInfo(*(recHits[1]),tTopo);

    float tol = 0.1; // cm
    float sc  = -1.; // scale, allow 'tol' of edge to count as outside
    LocalError le(tol*tol, tol*tol, tol*tol);

    // determine missing layers
    vector<int> missingLayers = getMissingLayers(getLayer(*(recHits[0]),tTopo),
                                                 getLayer(*(recHits[1]),tTopo));

    for(vector<int>::const_iterator missingLayer = missingLayers.begin();
                                    missingLayer!= missingLayers.end();
                                    missingLayer++)
    {
      int il = *missingLayer - 1;

      // propagate to missing layer
      FreeTrajectoryState fts = getTrajectory(*track);
      TrajectoryStateOnSurface tsos;
      float rz = 0.;
  
      if(il < BPix3)
      { // barrel
        const BarrelDetLayer * layer =
          (theGSTracker->pixelBarrelLayers())[il];
  
        tsos = thePropagator->propagate(fts, layer->surface());

        if(tsos.isValid())
          rz = tsos.globalPosition().z();
      }
      else
      { // endcap
        const ForwardDetLayer * layer;
        if(il - BPix3 < 2)
          layer = (theGSTracker->negPixelForwardLayers())[il - BPix3    ];
        else
          layer = (theGSTracker->posPixelForwardLayers())[il - BPix3 - 2];

        tsos = thePropagator->propagate(fts, layer->surface());

        if(tsos.isValid())
          rz = tsos.globalPosition().perp();
      }

      if(tsos.isValid())
      {
        float phi = spin(tsos.globalPosition().phi());

        // check close dets
        vector<const GeomDet *> closeDets =
          getCloseDets(il, rz ,rzBounds[il], phi,phBounds[il], tTopo);

        for(vector<const GeomDet *>::const_iterator det = closeDets.begin();
                                                    det!= closeDets.end();
                                                    det++)
        {
          TrajectoryStateOnSurface tsos =
            thePropagator->propagate(fts, (*det)->surface());

          if(tsos.isValid())
            if((*det)->surface().bounds().inside(tsos.localPosition(), le, sc))
              hasGap = false;
        }
      }
    }
  }


  if(hasGap)
    LogTrace("MinBiasTracking") << " [ValidHitPairFilter] has gap --> good";
  else
    LogTrace("MinBiasTracking") << " [ValidHitPairFilter] no gap --> rejected";

  return hasGap;
}

void ValidHitPairFilter::update ( const edm::Event &  ev, const edm::EventSetup &  es  )

overridevirtual

Implements PixelTrackFilter.

Definition at line 61 of file ValidHitPairFilter.cc.

References begin, BPix3, end, edm::EventSetup::get(), i, LogTrace, sysUtil::pid, edm::ESHandle< class >::product(), python.multivaluedict::sort(), and spin().

Referenced by progressbar.ProgressBar::__next__(), MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), Vispa.Gui.VispaWidget.VispaWidget::autosize(), Vispa.Views.LineDecayView.LineDecayContainer::createObject(), Vispa.Views.LineDecayView.LineDecayContainer::deselectAllObjects(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::deselectAllWidgets(), Vispa.Gui.VispaWidget.VispaWidget::enableAutosizing(), progressbar.ProgressBar::finish(), Vispa.Gui.MenuWidget.MenuWidget::leaveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseMoveEvent(), Vispa.Gui.MenuWidget.MenuWidget::mouseMoveEvent(), Vispa.Views.LineDecayView.LineDecayContainer::mouseMoveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseReleaseEvent(), Vispa.Views.LineDecayView.LineDecayContainer::objectMoved(), MatrixUtil.Steps::overwrite(), Vispa.Views.LineDecayView.LineDecayContainer::removeObject(), Vispa.Gui.ConnectableWidget.ConnectableWidget::removePorts(), Vispa.Gui.FindDialog.FindDialog::reset(), Vispa.Gui.PortConnection.PointToPointConnection::select(), Vispa.Gui.VispaWidget.VispaWidget::select(), Vispa.Views.LineDecayView.LineDecayContainer::select(), Vispa.Gui.VispaWidget.VispaWidget::setText(), Vispa.Gui.VispaWidget.VispaWidget::setTitle(), Vispa.Gui.ZoomableWidget.ZoomableWidget::setZoom(), Vispa.Views.LineDecayView.LineDecayContainer::setZoom(), and Vispa.Gui.PortConnection.PointToPointConnection::updateConnection().

                                                                           {
  //Retrieve tracker topology from geometry
  edm::ESHandle<TrackerTopology> tTopo;
  es.get<TrackerTopologyRcd>().get(tTopo);


  // Get tracker
  edm::ESHandle<TrackerGeometry> trackerHandle;
  es.get<TrackerDigiGeometryRecord>().get(trackerHandle);
  theTracker = trackerHandle.product();

  // Get geometric search tracker
  edm::ESHandle<GeometricSearchTracker> geometricSearchTrackerHandle;
  es.get<TrackerRecoGeometryRecord>().get(geometricSearchTrackerHandle);
  theGSTracker = geometricSearchTrackerHandle.product();

  // Get magnetic field
  edm::ESHandle<MagneticField> magneticFieldHandle;
  es.get<IdealMagneticFieldRecord>().get(magneticFieldHandle);
  theMagneticField = magneticFieldHandle.product();

  // Get propagator
  edm::ESHandle<Propagator> thePropagatorHandle;
  es.get<TrackingComponentsRecord>().get("AnalyticalPropagator",
                                          thePropagatorHandle);
  thePropagator = thePropagatorHandle.product();

  // Bounds, hardwired FIXME
  rzBounds[0].resize(8); phBounds[0].resize(20);
  rzBounds[1].resize(8); phBounds[1].resize(32);
  rzBounds[2].resize(8); phBounds[2].resize(44);

  rzBounds[3].resize(7); phBounds[3].resize(24);
  rzBounds[4].resize(7); phBounds[4].resize(24);
  rzBounds[5].resize(7); phBounds[5].resize(24);
  rzBounds[6].resize(7); phBounds[6].resize(24);

  LogTrace("MinBiasTracking")
    << " [ValidHitPairFilter] initializing pixel barrel";

  for(TrackerGeometry::DetContainer::const_iterator
      det = theTracker->detsPXB().begin();
      det!= theTracker->detsPXB().end(); det++)
  {
    
    DetId pid=(*det)->geographicalId();
    int il  = tTopo->pxbLayer(pid)  - 1;
    int irz = tTopo->pxbModule(pid) - 1;
    int iph = tTopo->pxbLadder(pid) - 1;
  
    rzBounds[il][irz] = (*det)->position().z();
    phBounds[il][iph] = spin((*det)->position().phi());
  }

  LogTrace("MinBiasTracking")
    << " [ValidHitPairFilter] initializing pixel endcap";

  for(TrackerGeometry::DetContainer::const_iterator
      det = theTracker->detsPXF().begin();
      det!= theTracker->detsPXF().end(); det++)
  {
    

    DetId pid=(*det)->geographicalId();
    int il  = BPix3 + ((tTopo->pxfSide(pid)  -1) << 1) + (tTopo->pxfDisk(pid) -1);
    int irz =         ((tTopo->pxfModule(pid)-1) << 1) + (tTopo->pxfPanel(pid)-1);
    int iph =          (tTopo->pxfBlade(pid) -1);

    rzBounds[il][irz] = (*det)->position().perp();
    phBounds[il][iph] = spin((*det)->position().phi());
  }

  for(int i = 0; i < 7; i++)
  {
    sort(rzBounds[i].begin(), rzBounds[i].end());
    sort(phBounds[i].begin(), phBounds[i].end());
  }
}

Member Data Documentation

const std::vector<DetLayer *> ValidHitPairFilter::detLayers

private

Definition at line 43 of file ValidHitPairFilter.h.

std::vector<float> ValidHitPairFilter::phBounds[7]

private

Definition at line 45 of file ValidHitPairFilter.h.

std::vector<float> ValidHitPairFilter::rzBounds[7]

private

Definition at line 44 of file ValidHitPairFilter.h.

const GeometricSearchTracker* ValidHitPairFilter::theGSTracker

private

Definition at line 39 of file ValidHitPairFilter.h.

const MagneticField* ValidHitPairFilter::theMagneticField

private

Definition at line 40 of file ValidHitPairFilter.h.

const Propagator* ValidHitPairFilter::thePropagator

private

Definition at line 41 of file ValidHitPairFilter.h.

const TrackerGeometry* ValidHitPairFilter::theTracker

private

Definition at line 38 of file ValidHitPairFilter.h.