#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 ()

std::vector< Road > const &	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 ( const edm::ParameterSet & param )

explicit

Definition at line 29 of file RPixRoadFinder.cc.

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

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

RPixRoadFinder::~RPixRoadFinder ( )

override

Definition at line 38 of file RPixRoadFinder.cc.

38 {}

Member Function Documentation

void RPixRoadFinder::findPattern ( )

overridevirtual

1mm

Implements RPixDetPatternFinder.

Definition at line 42 of file RPixRoadFinder.cc.

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

                                  {
   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
 }