RPixRoadFinder Class Reference

#include <RPixRoadFinder.h>

Inheritance diagram for RPixRoadFinder:

Public Member Functions
void	findPattern () override
	RPixRoadFinder (const edm::ParameterSet &param)
	~RPixRoadFinder () override
Public Member Functions inherited from RPixDetPatternFinder
void	clear ()
const std::vector< Road > &	getPatterns () const
	RPixDetPatternFinder (edm::ParameterSet const &parameterSet)
void	setGeometry (const CTPPSGeometry *geometry)
void	setHits (const edm::DetSetVector< CTPPSPixelRecHit > *hitVector)
virtual	~RPixDetPatternFinder ()

Private Member Functions
void	run (const edm::DetSetVector< CTPPSPixelRecHit > &input, const CTPPSGeometry &geometry, std::vector< Road > &roads)

Private Attributes
unsigned int	maxRoadSize_
unsigned int	minRoadSize_
double	roadRadius_
int	verbosity_

Additional Inherited Members
Public Types inherited from RPixDetPatternFinder
typedef std::vector< PointInPlane >	Road
Protected Attributes inherited from RPixDetPatternFinder
const CTPPSGeometry *	geometry_
const edm::DetSetVector< CTPPSPixelRecHit > *	hitVector_
std::vector< Road >	patternVector_

Detailed Description

Definition at line 36 of file RPixRoadFinder.h.

Constructor & Destructor Documentation

RPixRoadFinder()

RPixRoadFinder::RPixRoadFinder ( const edm::ParameterSet &  param )

explicit

Definition at line 29 of file RPixRoadFinder.cc.

                                                                  : RPixDetPatternFinder(parameterSet) {
  verbosity_ = parameterSet.getUntrackedParameter<int>("verbosity");
  roadRadius_ = parameterSet.getParameter<double>("roadRadius");
  minRoadSize_ = parameterSet.getParameter<int>("minRoadSize");
  maxRoadSize_ = parameterSet.getParameter<int>("maxRoadSize");
}

References edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), maxRoadSize_, minRoadSize_, edm::parameterSet(), roadRadius_, and verbosity_.

~RPixRoadFinder()

RPixRoadFinder::~RPixRoadFinder ( )

override

Definition at line 38 of file RPixRoadFinder.cc.

{}

Member Function Documentation

findPattern()

void RPixRoadFinder::findPattern ( )

overridevirtual

1mm

Implements RPixDetPatternFinder.

Definition at line 42 of file RPixRoadFinder.cc.

                                 {
  Road temp_all_hits;
  temp_all_hits.clear();
 
  // convert local hit sto global and push them to a vector
  for (const auto& ds_rh2 : *hitVector_) {
    const auto myid = CTPPSPixelDetId(ds_rh2.id);
    for (const auto& it_rh : ds_rh2.data) {
      CTPPSGeometry::Vector localV(it_rh.point().x(), it_rh.point().y(), it_rh.point().z());
      const auto& globalV = geometry_->localToGlobal(ds_rh2.id, localV);
      math::Error<3>::type localError;
      localError[0][0] = it_rh.error().xx();
      localError[0][1] = it_rh.error().xy();
      localError[0][2] = 0.;
      localError[1][0] = it_rh.error().xy();
      localError[1][1] = it_rh.error().yy();
      localError[1][2] = 0.;
      localError[2][0] = 0.;
      localError[2][1] = 0.;
      localError[2][2] = 0.;
      if (verbosity_ > 2)
        edm::LogInfo("RPixRoadFinder") << "Hits = " << ds_rh2.data.size();
 
      DetGeomDesc::RotationMatrix theRotationMatrix = geometry_->sensor(myid)->rotation();
      AlgebraicMatrix33 theRotationTMatrix;
      theRotationMatrix.GetComponents(theRotationTMatrix(0, 0),
                                      theRotationTMatrix(0, 1),
                                      theRotationTMatrix(0, 2),
                                      theRotationTMatrix(1, 0),
                                      theRotationTMatrix(1, 1),
                                      theRotationTMatrix(1, 2),
                                      theRotationTMatrix(2, 0),
                                      theRotationTMatrix(2, 1),
                                      theRotationTMatrix(2, 2));
 
      math::Error<3>::type globalError = ROOT::Math::SimilarityT(theRotationTMatrix, localError);
      temp_all_hits.emplace_back(PointInPlane{globalV, globalError, it_rh, myid});
    }
  }
 
  Road::iterator it_gh1 = temp_all_hits.begin();
  Road::iterator it_gh2;
 
  patternVector_.clear();
 
  //look for points near wrt each other
  // starting algorithm
  while (it_gh1 != temp_all_hits.end() && temp_all_hits.size() >= minRoadSize_) {
    Road temp_road;
 
    it_gh2 = it_gh1;
 
    const auto currPoint = it_gh1->globalPoint;
    CTPPSPixelDetId currDet = CTPPSPixelDetId(it_gh1->detId);
 
    while (it_gh2 != temp_all_hits.end()) {
      bool same_pot = false;
      CTPPSPixelDetId tmpGh2Id = CTPPSPixelDetId(it_gh2->detId);
      if (currDet.rpId() == tmpGh2Id.rpId())
        same_pot = true;
      const auto subtraction = currPoint - it_gh2->globalPoint;
 
      if (subtraction.Rho() < roadRadius_ && same_pot) {  
        temp_road.push_back(*it_gh2);
        temp_all_hits.erase(it_gh2);
      } else {
        ++it_gh2;
      }
    }
 
    if (temp_road.size() >= minRoadSize_ && temp_road.size() < maxRoadSize_)
      patternVector_.push_back(temp_road);
  }
  // end of algorithm
}

References RPixDetPatternFinder::geometry_, RPixDetPatternFinder::hitVector_, CTPPSGeometry::localToGlobal(), maxRoadSize_, minRoadSize_, RPixDetPatternFinder::patternVector_, roadRadius_, DetGeomDesc::rotation(), CTPPSDetId::rpId(), CTPPSGeometry::sensor(), and verbosity_.

run()

void RPixRoadFinder::run ( const edm::DetSetVector< CTPPSPixelRecHit > &  input, const CTPPSGeometry &  geometry, std::vector< Road > &  roads  )

private

Member Data Documentation

maxRoadSize_

unsigned int RPixRoadFinder::maxRoadSize_

private

Definition at line 46 of file RPixRoadFinder.h.

Referenced by findPattern(), and RPixRoadFinder().

minRoadSize_

unsigned int RPixRoadFinder::minRoadSize_

private

Definition at line 45 of file RPixRoadFinder.h.

Referenced by findPattern(), and RPixRoadFinder().

roadRadius_

double RPixRoadFinder::roadRadius_

private

Definition at line 44 of file RPixRoadFinder.h.

Referenced by findPattern(), and RPixRoadFinder().

verbosity_

int RPixRoadFinder::verbosity_

private

Definition at line 43 of file RPixRoadFinder.h.

Referenced by findPattern(), and RPixRoadFinder().