#include <ClusterShapeHitFilter.h>

Classes
struct	PixelData

Public Types
typedef TrajectoryFilter::Record	Record

Public Member Functions
	ClusterShapeHitFilter (const TrackerGeometry theTracker_, const MagneticField theMagneticField_, const SiPixelLorentzAngle theSiPixelLorentzAngle_, const SiStripLorentzAngle theSiStripLorentzAngle_)

bool	getSizes (const SiPixelRecHit &recHit, const LocalVector &ldir, int &part, std::vector< std::pair< int, int > > &meas, std::pair< float, float > &predr, PixelData const *pd=0) const

bool	getSizes (DetId detId, const SiStripCluster &cluster, const LocalVector &ldir, int &meas, float &pred) const

bool	getSizes (const SiStripRecHit2D &recHit, const LocalVector &ldir, int &meas, float &pred) const

bool	isCompatible (const SiPixelRecHit &recHit, const LocalVector &ldir, PixelData const *pd=0) const

bool	isCompatible (const SiPixelRecHit &recHit, const GlobalVector &gdir, PixelData const *pd=0) const

bool	isCompatible (DetId detId, const SiStripCluster &cluster, const LocalVector &ldir) const

bool	isCompatible (DetId detId, const SiStripCluster &cluster, const GlobalVector &gdir) const

bool	isCompatible (const SiStripRecHit2D &recHit, const LocalVector &ldir) const

bool	isCompatible (const SiStripRecHit2D &recHit, const GlobalVector &gdir) const

	~ClusterShapeHitFilter ()

Private Member Functions
	ClusterShapeHitFilter ()

void	fillPixelData ()

std::pair< float, float >	getCotangent (const PixelGeomDetUnit *pixelDet) const

float	getCotangent (const StripGeomDetUnit *stripDet) const

std::pair< float, float >	getDrift (const PixelGeomDetUnit *pixelDet) const

float	getDrift (const StripGeomDetUnit *stripDet) const

const PixelData &	getpd (const SiPixelRecHit &recHit, PixelData const *pd=0) const

bool	isNormalOriented (const GeomDetUnit *geomDet) const

void	loadPixelLimits ()

void	loadStripLimits ()

Private Attributes
std::unordered_map< unsigned int, PixelData >	pixelData

PixelLimits	pixelLimits [PixelKeys::N+1]

StripLimits	stripLimits [StripKeys::N+1]

float	theAngle [6]

const MagneticField *	theMagneticField

const SiPixelLorentzAngle *	theSiPixelLorentzAngle

const SiStripLorentzAngle *	theSiStripLorentzAngle

const TrackerGeometry *	theTracker

Detailed Description

Definition at line 143 of file ClusterShapeHitFilter.h.

Member Typedef Documentation

typedef TrajectoryFilter::Record ClusterShapeHitFilter::Record

Definition at line 155 of file ClusterShapeHitFilter.h.

Constructor & Destructor Documentation

ClusterShapeHitFilter::ClusterShapeHitFilter	(	const TrackerGeometry *	theTracker_,
		const MagneticField *	theMagneticField_,
		const SiPixelLorentzAngle *	theSiPixelLorentzAngle_,
		const SiStripLorentzAngle *	theSiStripLorentzAngle_
	)

Definition at line 44 of file ClusterShapeHitFilter.cc.

    : theTracker(theTracker_),
      theMagneticField(theMagneticField_),
      theSiPixelLorentzAngle(theSiPixelLorentzAngle_),
      theSiStripLorentzAngle(theSiStripLorentzAngle_)
 
 {
   // Load pixel limits
   loadPixelLimits();
   fillPixelData();
 
   // Load strip limits
   loadStripLimits();
 }

ClusterShapeHitFilter::~ClusterShapeHitFilter ( )

Definition at line 63 of file ClusterShapeHitFilter.cc.

64 {

65 }

ClusterShapeHitFilter::ClusterShapeHitFilter ( )

inlineprivate

Definition at line 202 of file ClusterShapeHitFilter.h.

202 {}

Member Function Documentation

void ClusterShapeHitFilter::fillPixelData ( )

private

Definition at line 138 of file ClusterShapeHitFilter.cc.

References ClusterShapeHitFilter::PixelData::cotangent, ClusterShapeHitFilter::PixelData::det, ClusterShapeHitFilter::PixelData::drift, GeomDet::geographicalId(), and ClusterShapeHitFilter::PixelData::part.

                                           {
 
   //barrel
   for (auto det : theTracker->detsPXB()) {
     // better not to fail..
     const PixelGeomDetUnit * pixelDet =
       dynamic_cast<const PixelGeomDetUnit*>(det);
     assert(pixelDet);
     PixelData & pd = pixelData[pixelDet->geographicalId()];
     pd.det = pixelDet;
     pd.part=0;
     pd.cotangent=getCotangent(pixelDet);
     pd.drift=getDrift(pixelDet);
   }
 
   //endcap
   for (auto det : theTracker->detsPXF()) {
     // better not to fail..
     const PixelGeomDetUnit * pixelDet =
       dynamic_cast<const PixelGeomDetUnit*>(det);
     assert(pixelDet);
     PixelData & pd = pixelData[pixelDet->geographicalId()];
     pd.det = pixelDet;
     pd.part=1;
     pd.cotangent=getCotangent(pixelDet);
     pd.drift=getDrift(pixelDet);
   }
 
 }

pair< float, float > ClusterShapeHitFilter::getCotangent ( const PixelGeomDetUnit * pixelDet ) const

private

Definition at line 171 of file ClusterShapeHitFilter.cc.

References BoundSurface::bounds(), PixelTopology::pitch(), PixelGeomDetUnit::specificTopology(), GeomDet::surface(), and Bounds::thickness().

 {
   pair<float,float> cotangent;
 
   cotangent.first  = pixelDet->surface().bounds().thickness() /
                      pixelDet->specificTopology().pitch().first;
   cotangent.second = pixelDet->surface().bounds().thickness() /
                      pixelDet->specificTopology().pitch().second;
 
   return cotangent;
 }

float ClusterShapeHitFilter::getCotangent ( const StripGeomDetUnit * stripDet ) const

private

Definition at line 185 of file ClusterShapeHitFilter.cc.

References BoundSurface::bounds(), StripTopology::localPitch(), StripGeomDetUnit::specificTopology(), GeomDet::surface(), and Bounds::thickness().

 {
   // FIXME may be problematic in case of RadialStriptolopgy
   return stripDet->surface().bounds().thickness() /
          stripDet->specificTopology().localPitch(LocalPoint(0,0,0));
 }

pair< float, float > ClusterShapeHitFilter::getDrift ( const PixelGeomDetUnit * pixelDet ) const

private

Definition at line 194 of file ClusterShapeHitFilter.cc.

References dir, GeomDet::geographicalId(), GloballyPositioned< T >::position(), DetId::rawId(), GeomDet::surface(), toLocal(), PV3DBase< T, PVType, FrameType >::x(), and PV3DBase< T, PVType, FrameType >::y().

 {
   LocalVector lBfield =
      (pixelDet->surface()).toLocal(
       theMagneticField->inTesla(
       pixelDet->surface().position()));
 
   double theTanLorentzAnglePerTesla =
          theSiPixelLorentzAngle->getLorentzAngle(
            pixelDet->geographicalId().rawId());
 
   pair<float,float> dir;
   dir.first  = - theTanLorentzAnglePerTesla * lBfield.y();
   dir.second =   theTanLorentzAnglePerTesla * lBfield.x();
 
   return dir;
 }

float ClusterShapeHitFilter::getDrift ( const StripGeomDetUnit * stripDet ) const

private

Definition at line 213 of file ClusterShapeHitFilter.cc.

References dir, GeomDet::geographicalId(), GloballyPositioned< T >::position(), DetId::rawId(), GeomDet::surface(), toLocal(), and PV3DBase< T, PVType, FrameType >::y().

 {
   LocalVector lBfield =
      (stripDet->surface()).toLocal(
       theMagneticField->inTesla(
       stripDet->surface().position()));
     
   double theTanLorentzAnglePerTesla =
          theSiStripLorentzAngle->getLorentzAngle(
            stripDet->geographicalId().rawId());
     
   float dir = theTanLorentzAnglePerTesla * lBfield.y();
 
   return dir;
 }

const PixelData& ClusterShapeHitFilter::getpd	(	const SiPixelRecHit &	recHit,
		PixelData const *	pd = `0`
	)		const

inlineprivate

Definition at line 204 of file ClusterShapeHitFilter.h.

References TrackingRecHit::geographicalId(), AlCaHLTBitMon_ParallelJobs::p, and pixelData.

                                                                                              {
     if (pd) return *pd;
     // Get detector
     DetId id = recHit.geographicalId();
     auto p = pixelData.find(id);
     return (*p).second;
   }

bool ClusterShapeHitFilter::getSizes	(	const SiPixelRecHit &	recHit,
		const LocalVector &	ldir,
		int &	part,
		std::vector< std::pair< int, int > > &	meas,
		std::pair< float, float > &	predr,
		PixelData const *	pd = `0`
	)		const

Definition at line 252 of file ClusterShapeHitFilter.cc.

References ClusterShapeHitFilter::PixelData::cotangent, data, ClusterShapeHitFilter::PixelData::det, ClusterShape::determineShape(), shallow::drift(), ClusterShapeHitFilter::PixelData::drift, ClusterData::isComplete, ClusterData::isStraight, ClusterShapeHitFilter::PixelData::part, alignCSCRings::s, ClusterData::size, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by getSizes().

 {
   // Get detector
   const PixelData & pd = getpd(recHit,ipd);
 
   // Get shape information
   ClusterData data;
   ClusterShape theClusterShape;
   theClusterShape.determineShape(*pd.det, recHit, data);
   bool usable = (data.isStraight && data.isComplete);
  
   // Usable?
   //if(usable)
   {
     part = pd.part;
 
     // Predicted size
     pred.first  = ldir.x() / ldir.z();
     pred.second = ldir.y() / ldir.z();
 
     if(data.size.front().second < 0)
       pred.second = - pred.second;
 
     meas.reserve(data.size.size());
     for(vector<pair<int,int> >::const_iterator s = data.size.begin();
     s!= data.size.end(); s++)
       {
     meas.push_back(*s);
     
     if(data.size.front().second < 0)
       meas.back().second = - meas.back().second;
       }
 
     // Take out drift 
     std::pair<float,float> const & drift = pd.drift;
     pred.first  += drift.first;
     pred.second += drift.second;
 
     // Apply cotangent
     std::pair<float,float> const & cotangent = pd.cotangent;
     pred.first  *= cotangent.first;
     pred.second *= cotangent.second;
   }
 
   // Usable?
   return usable;
 }

bool ClusterShapeHitFilter::getSizes	(	DetId	detId,
		const SiStripCluster &	cluster,
		const LocalVector &	ldir,
		int &	meas,
		float &	pred
	)		const

Definition at line 345 of file ClusterShapeHitFilter.cc.

References SiStripCluster::amplitudes(), shallow::drift(), SiStripCluster::firstStrip(), StripTopology::nstrips(), StripGeomDetUnit::specificTopology(), PV3DBase< T, PVType, FrameType >::x(), and PV3DBase< T, PVType, FrameType >::z().

 {
   // Get detector
   const StripGeomDetUnit* stripDet =
     dynamic_cast<const StripGeomDetUnit*> (theTracker->idToDet(id));
 
   // Measured width
   meas   = cluster.amplitudes().size();
 
   // Usable?
   int fs = cluster.firstStrip();
   int ns = stripDet->specificTopology().nstrips();
   // bool usable = (fs > 1 && fs + meas - 1 < ns);
   bool usable = (fs >= 1 && fs + meas - 1 <= ns);
 
   // Usable?
   //if(usable)
   {
     // Predicted width
     pred = ldir.x() / ldir.z();
   
     // Take out drift
     float drift = getDrift(stripDet);
     pred += drift;
   
     // Apply cotangent
     pred *= getCotangent(stripDet);
   }
 
   return usable;
 }   

bool ClusterShapeHitFilter::getSizes	(	const SiStripRecHit2D &	recHit,
		const LocalVector &	ldir,
		int &	meas,
		float &	pred
	)		const

inline

Definition at line 180 of file ClusterShapeHitFilter.h.

References TrackingRecHit::geographicalId(), getSizes(), and TrackerSingleRecHit::stripCluster().

                                      {
     return getSizes(recHit.geographicalId(), recHit.stripCluster(), ldir, meas, pred);
   }

bool ClusterShapeHitFilter::isCompatible	(	const SiPixelRecHit &	recHit,
		const LocalVector &	ldir,
		PixelData const *	pd = `0`
	)		const

Definition at line 303 of file ClusterShapeHitFilter.cc.

References PixelKeys::isValid(), combine::key, and m.

Referenced by isCompatible().

 {
  // Get detector
   const PixelData & pd = getpd(recHit,ipd);
 
   int part;
   vector<pair<int,int> > meas;
   pair<float,float> pred;
 
   if(getSizes(recHit, ldir, part,meas, pred,&pd))
   {
     for(vector<pair<int,int> >::const_iterator m = meas.begin();
                                                m!= meas.end(); m++)
     {
       PixelKeys key(part, (*m).first, (*m).second);
       if (!key.isValid()) return true; // FIXME original logic
       if (pixelLimits[key].isInside(pred)) return true;
     }
     // none of the choices worked
     return false;
   }
   // not usable
   return true;
 }

bool ClusterShapeHitFilter::isCompatible	(	const SiPixelRecHit &	recHit,
		const GlobalVector &	gdir,
		PixelData const *	pd = `0`
	)		const

Definition at line 330 of file ClusterShapeHitFilter.cc.

References ClusterShapeHitFilter::PixelData::det, and GeomDet::toLocal().

 {
  // Get detector
   const PixelData & pd = getpd(recHit,ipd);
 
   LocalVector ldir =pd.det->toLocal(gdir);
 
   return isCompatible(recHit, ldir,&pd);
 }

bool ClusterShapeHitFilter::isCompatible	(	DetId	detId,
		const SiStripCluster &	cluster,
		const LocalVector &	ldir
	)		const

Definition at line 381 of file ClusterShapeHitFilter.cc.

References StripKeys::isValid(), and combine::key.

 {
   int meas;
   float pred;
 
   if(getSizes(detId, cluster, ldir, meas, pred))
   {
     StripKeys key(meas);
     if (key.isValid())
       return stripLimits[key].isInside(pred);
   }
 
   // Not usable or no limits
   return true;
 }

bool ClusterShapeHitFilter::isCompatible	(	DetId	detId,
		const SiStripCluster &	cluster,
		const GlobalVector &	gdir
	)		const

Definition at line 400 of file ClusterShapeHitFilter.cc.

 {
   LocalVector ldir = theTracker->idToDet(detId)->toLocal(gdir);
   return isCompatible(detId, cluster, ldir);
 }

bool ClusterShapeHitFilter::isCompatible	(	const SiStripRecHit2D &	recHit,
		const LocalVector &	ldir
	)		const

inline

Definition at line 190 of file ClusterShapeHitFilter.h.

References TrackingRecHit::geographicalId(), isCompatible(), and TrackerSingleRecHit::stripCluster().

                                                     {
             return isCompatible(recHit.geographicalId(), recHit.stripCluster(), ldir);
   }

bool ClusterShapeHitFilter::isCompatible	(	const SiStripRecHit2D &	recHit,
		const GlobalVector &	gdir
	)		const

inline

Definition at line 194 of file ClusterShapeHitFilter.h.

References TrackingRecHit::geographicalId(), isCompatible(), and TrackerSingleRecHit::stripCluster().

                                                       {
             return isCompatible(recHit.geographicalId(), recHit.stripCluster(), gdir);
   } 

bool ClusterShapeHitFilter::isNormalOriented ( const GeomDetUnit * geomDet ) const

private

Definition at line 232 of file ClusterShapeHitFilter.cc.

References GeomDetType::isBarrel(), PV3DBase< T, PVType, FrameType >::perp(), pos, makeMuonMisalignmentScenario::rot, GeomDet::toGlobal(), GeomDetUnit::type(), and PV3DBase< T, PVType, FrameType >::z().

 {
   if(geomDet->type().isBarrel())
   { // barrel
     float perp0 = geomDet->toGlobal( Local3DPoint(0.,0.,0.) ).perp();
     float perp1 = geomDet->toGlobal( Local3DPoint(0.,0.,1.) ).perp();
     return (perp1 > perp0);
   }
   else
   { // endcap
     float rot = geomDet->toGlobal( LocalVector (0.,0.,1.) ).z();
     float pos = geomDet->toGlobal( Local3DPoint(0.,0.,0.) ).z();
     return (rot * pos > 0);
   }
 }

void ClusterShapeHitFilter::loadPixelLimits ( )

private

Definition at line 68 of file ClusterShapeHitFilter.cc.

References b, f, edm::FileInPath::fullPath(), gen::k, combine::key, and LogTrace.

 {
   edm::FileInPath
     fileInPath("RecoPixelVertexing/PixelLowPtUtilities/data/pixelShape.par");
   ifstream inFile(fileInPath.fullPath().c_str());
 
 
   while(inFile.eof() == false)
   {
     int part,dx,dy;
 
     inFile >> part; // 0or 1
     inFile >> dx;   // 0 to 10
     inFile >> dy;   // 0 to 15 ...
 
     const PixelKeys key(part,dx,dy);
     auto & pl = pixelLimits[key];
 
     for(int b = 0; b<2 ; b++) // branch
     for(int d = 0; d<2 ; d++) // direction
     for(int k = 0; k<2 ; k++) // lower and upper
       inFile >> pl.data[b][d][k];
 
 
     double f;
     int d;
 
     inFile >> f; // density
     inFile >> d; // points
     inFile >> f; // density
     inFile >> d; // points
   }
   
   inFile.close();
   
   LogTrace("MinBiasTracking|ClusterShapeHitFilter")
     << " [ClusterShapeHitFilter] pixel-cluster-shape filter loaded";
  }

void ClusterShapeHitFilter::loadStripLimits ( )

private

Definition at line 108 of file ClusterShapeHitFilter.cc.

References b, edm::FileInPath::fullPath(), gen::k, combine::key, and LogTrace.

 {
   // Load strip
   edm::FileInPath
     fileInPath("RecoPixelVertexing/PixelLowPtUtilities/data/stripShape.par");
   ifstream inFile(fileInPath.fullPath().c_str());
 
   
   while(inFile.eof() == false)
   {
     int dx;
     inFile >> dx;
 
     StripKeys key(dx);
     auto & sl = stripLimits[key];
 
     for(int b = 0; b<2 ; b++) // branch
     for(int k = 0; k<2 ; k++) // lower and upper
       inFile >> sl.data[b][k];
 
   } 
   
   inFile.close();
   
   LogTrace("MinBiasTracking|ClusterShapeHitFilter")
     << " [ClusterShapeHitFilter] strip-cluster-width filter loaded";
 }

Member Data Documentation

std::unordered_map<unsigned int, PixelData> ClusterShapeHitFilter::pixelData

private

Definition at line 230 of file ClusterShapeHitFilter.h.

Referenced by getpd().

PixelLimits ClusterShapeHitFilter::pixelLimits[PixelKeys::N+1]

private

Definition at line 232 of file ClusterShapeHitFilter.h.

StripLimits ClusterShapeHitFilter::stripLimits[StripKeys::N+1]

private

Definition at line 234 of file ClusterShapeHitFilter.h.

float ClusterShapeHitFilter::theAngle[6]

private

Definition at line 236 of file ClusterShapeHitFilter.h.

const MagneticField* ClusterShapeHitFilter::theMagneticField

private

Definition at line 225 of file ClusterShapeHitFilter.h.

const SiPixelLorentzAngle* ClusterShapeHitFilter::theSiPixelLorentzAngle

private

Definition at line 227 of file ClusterShapeHitFilter.h.

const SiStripLorentzAngle* ClusterShapeHitFilter::theSiStripLorentzAngle

private

Definition at line 228 of file ClusterShapeHitFilter.h.

const TrackerGeometry* ClusterShapeHitFilter::theTracker

private

Definition at line 224 of file ClusterShapeHitFilter.h.

Classes

Public Types

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation