#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 40 of file RPixRoadFinder.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 30 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 42 of file RPixRoadFinder.cc.

42 {

43 }

Member Function Documentation

void RPixRoadFinder::findPattern ( )

overridevirtual

1mm

Implements RPixDetPatternFinder.

Definition at line 47 of file RPixRoadFinder.cc.

References RPixDetPatternFinder::geometry_, CTPPSDetId::getRPId(), CTPPSGeometry::getSensor(), RPixDetPatternFinder::hitVector_, CTPPSGeometry::localToGlobal(), maxRoadSize_, minRoadSize_, RPixDetPatternFinder::patternVector_, roadRadius_, DetGeomDesc::rotation(), 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){
       CLHEP::Hep3Vector localV(it_rh.getPoint().x(),it_rh.getPoint().y(),it_rh.getPoint().z() );
       CLHEP::Hep3Vector globalV = geometry_->localToGlobal(ds_rh2.id,localV);
       math::Error<3>::type localError;
       localError[0][0] = it_rh.getError().xx();
       localError[0][1] = it_rh.getError().xy();
       localError[0][2] =                  0.;
       localError[1][0] = it_rh.getError().xy();
       localError[1][1] = it_rh.getError().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_->getSensor(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);
       PointInPlane thePointAndRecHit = {globalV,globalError,it_rh,myid};
       temp_all_hits.push_back(thePointAndRecHit);
     }
 
   }
   
   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;
 
     CLHEP::Hep3Vector 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.getRPId() == tmpGh2Id.getRPId() ) same_pot = true;
       CLHEP::Hep3Vector subtraction = currPoint - it_gh2->globalPoint;
 
       if(subtraction.perp() < 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
 
 
 
 }