#include <ValidHitPairFilter.h>

Inheritance diagram for ValidHitPairFilter:

Public Member Functions
virtual bool	operator() (const reco::Track track, std::vector< const TrackingRecHit > recHits) const

	ValidHitPairFilter (const edm::ParameterSet &ps, const edm::EventSetup &es)

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 void	update (edm::Event &ev)

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

int	getLayer (const TrackingRecHit &recHit) 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 19 of file ValidHitPairFilter.h.

Constructor & Destructor Documentation

ValidHitPairFilter::ValidHitPairFilter	(	const edm::ParameterSet &	ps,
		const edm::EventSetup &	es
	)

Definition at line 52 of file ValidHitPairFilter.cc.

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

 {
   // 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++)
   {
     PXBDetId pid((*det)->geographicalId());
 
     int il  = pid.layer()  - 1;
     int irz = pid.module() - 1;
     int iph = pid.ladder() - 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++)
   {
     PXFDetId pid((*det)->geographicalId());
 
     int il  = BPix3 + ((pid.side()  -1) << 1) + (pid.disk() -1);
     int irz =         ((pid.module()-1) << 1) + (pid.panel()-1);
     int iph =          (pid.blade() -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());
   }
 }

ValidHitPairFilter::~ValidHitPairFilter ( )

virtual

Definition at line 127 of file ValidHitPairFilter.cc.

128 {

129 }

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

private

Definition at line 202 of file ValidHitPairFilter.cc.

References BPix3, errorMatrix2Lands_multiChannel::id, and LogTrace.

 {
   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<<")";
 
       PXBDetId id(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<<")";
 
       PXFDetId id(side,disk,blade,panel,module);
       dets.push_back(theTracker->idToDet(id));
     }
   }
 
   return dets;
 }

int ValidHitPairFilter::getLayer ( const TrackingRecHit & recHit ) const

private

Definition at line 131 of file ValidHitPairFilter.cc.

References BPix3, PXFDetId::disk(), TrackingRecHit::geographicalId(), errorMatrix2Lands_multiChannel::id, PXBDetId::layer(), evf::utils::pid, PixelSubdetector::PixelBarrel, and PXFDetId::side().

 {
   DetId id(recHit.geographicalId());
 
   if(id.subdetId() == int(PixelSubdetector::PixelBarrel))
   {
     PXBDetId pid(id);
     return pid.layer();
   }
   else
   {
     PXFDetId pid(id);
     return BPix3 + ((pid.side()-1) << 1) + pid.disk();
   }
 }

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

private

Definition at line 148 of file ValidHitPairFilter.cc.

References BPix1, BPix2, BPix3, trackerHits::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 182 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,
		std::vector< const TrackingRecHit * >	recHits
	)		const

virtual

Definition at line 259 of file ValidHitPairFilter.cc.

References BPix3, HitInfo::getInfo(), TrajectoryStateOnSurface::globalPosition(), TrajectoryStateOnSurface::isValid(), vdt::le, TrajectoryStateOnSurface::localPosition(), LogTrace, PV3DBase< T, PVType, FrameType >::perp(), phi, PV3DBase< T, PVType, FrameType >::phi(), 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]))
                                        << HitInfo::getInfo(*(recHits[1]));
 
     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])),
                                                  getLayer(*(recHits[1])));
 
     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]);
 
         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;
 }

Member Data Documentation

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

private

Definition at line 40 of file ValidHitPairFilter.h.

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

private

Definition at line 42 of file ValidHitPairFilter.h.

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

private

Definition at line 41 of file ValidHitPairFilter.h.

const GeometricSearchTracker* ValidHitPairFilter::theGSTracker

private

Definition at line 36 of file ValidHitPairFilter.h.

const MagneticField* ValidHitPairFilter::theMagneticField

private

Definition at line 37 of file ValidHitPairFilter.h.

const Propagator* ValidHitPairFilter::thePropagator

private

Definition at line 38 of file ValidHitPairFilter.h.

const TrackerGeometry* ValidHitPairFilter::theTracker

private

Definition at line 35 of file ValidHitPairFilter.h.

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation