CMS 3D CMS Logo
List of all members | Public Member Functions | Private Member Functions | Private Attributes
RPixRoadFinder Class Reference

#include <RPixRoadFinder.h>

Inheritance diagram for RPixRoadFinder:
RPixDetPatternFinder

Public Member Functions

void findPattern (bool *is2planepot) 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_
 
double roadRadiusBadPot_
 
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 34 of file RPixRoadFinder.h.

Constructor & Destructor Documentation

◆ RPixRoadFinder()

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

Definition at line 18 of file RPixRoadFinder.cc.

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

18  : RPixDetPatternFinder(parameterSet) {
19  verbosity_ = parameterSet.getUntrackedParameter<int>("verbosity");
20  roadRadius_ = parameterSet.getParameter<double>("roadRadius");
21  minRoadSize_ = parameterSet.getParameter<int>("minRoadSize");
22  maxRoadSize_ = parameterSet.getParameter<int>("maxRoadSize");
23  roadRadiusBadPot_ = parameterSet.getParameter<double>("roadRadiusBadPot");
24 }
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
Definition: ParameterSet.h:303
RPixRoadFinder::maxRoadSize_
unsigned int maxRoadSize_
Definition: RPixRoadFinder.h:44
RPixDetPatternFinder::RPixDetPatternFinder
RPixDetPatternFinder(edm::ParameterSet const &parameterSet)
Definition: RPixDetPatternFinder.h:28
edm::parameterSet
ParameterSet const & parameterSet(StableProvenance const &provenance, ProcessHistory const &history)
Definition: Provenance.cc:11
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
RPixRoadFinder::roadRadiusBadPot_
double roadRadiusBadPot_
Definition: RPixRoadFinder.h:45
RPixRoadFinder::minRoadSize_
unsigned int minRoadSize_
Definition: RPixRoadFinder.h:43

◆ ~RPixRoadFinder()

RPixRoadFinder::~RPixRoadFinder ( )
override

Definition at line 28 of file RPixRoadFinder.cc.

28 {}

Member Function Documentation

◆ findPattern()

void RPixRoadFinder::findPattern ( bool *  is2planepot)
overridevirtual

1mm

Implements RPixDetPatternFinder.

Definition at line 32 of file RPixRoadFinder.cc.

References mps_fire::end, RPixDetPatternFinder::geometry_, RPixDetPatternFinder::hitVector_, mps_fire::i, CTPPSGeometry::localToGlobal(), maxRoadSize_, minRoadSize_, RPixDetPatternFinder::patternVector_, roadRadius_, roadRadiusBadPot_, DetGeomDesc::rotation(), CTPPSDetId::rpId(), CTPPSGeometry::sensor(), findQualityFiles::size, and verbosity_.

32  {
33  Road temp_all_hits;
34  Road temp_all_hits_2PlanePot[4];
35 
36  // convert local hit sto global and push them to a vector
37  for (const auto& ds_rh2 : *hitVector_) {
38  const auto myid = CTPPSPixelDetId(ds_rh2.id);
39  for (const auto& it_rh : ds_rh2.data) {
40  CTPPSGeometry::Vector localV(it_rh.point().x(), it_rh.point().y(), it_rh.point().z());
41  const auto& globalV = geometry_->localToGlobal(ds_rh2.id, localV);
42  math::Error<3>::type localError;
43  localError[0][0] = it_rh.error().xx();
44  localError[0][1] = it_rh.error().xy();
45  localError[0][2] = 0.;
46  localError[1][0] = it_rh.error().xy();
47  localError[1][1] = it_rh.error().yy();
48  localError[1][2] = 0.;
49  localError[2][0] = 0.;
50  localError[2][1] = 0.;
51  localError[2][2] = 0.;
52  if (verbosity_ > 2)
53  edm::LogInfo("RPixRoadFinder") << "Hits = " << ds_rh2.data.size();
54 
55  DetGeomDesc::RotationMatrix theRotationMatrix = geometry_->sensor(myid)->rotation();
56  AlgebraicMatrix33 theRotationTMatrix;
57  theRotationMatrix.GetComponents(theRotationTMatrix(0, 0),
58  theRotationTMatrix(0, 1),
59  theRotationTMatrix(0, 2),
60  theRotationTMatrix(1, 0),
61  theRotationTMatrix(1, 1),
62  theRotationTMatrix(1, 2),
63  theRotationTMatrix(2, 0),
64  theRotationTMatrix(2, 1),
65  theRotationTMatrix(2, 2));
66 
67  math::Error<3>::type globalError = ROOT::Math::SimilarityT(theRotationTMatrix, localError);
68 
69  // create new collections for bad (2 planes) and good pots
70 
71  if (is2PlanePot[0] == true && myid.arm() == 0 && myid.station() == 2) { // 45-220
72  temp_all_hits_2PlanePot[0].emplace_back(PointInPlane{globalV, globalError, it_rh, myid});
73  } else if (is2PlanePot[1] == true && myid.arm() == 0 && myid.station() == 0) { // 45-210
74  temp_all_hits_2PlanePot[1].emplace_back(PointInPlane{globalV, globalError, it_rh, myid});
75  } else if (is2PlanePot[2] == true && myid.arm() == 1 && myid.station() == 0) { // 56-210
76  temp_all_hits_2PlanePot[2].emplace_back(PointInPlane{globalV, globalError, it_rh, myid});
77  } else if (is2PlanePot[3] == true && myid.arm() == 1 && myid.station() == 2) { // 56-220
78  temp_all_hits_2PlanePot[3].emplace_back(PointInPlane{globalV, globalError, it_rh, myid});
79  } else {
80  temp_all_hits.emplace_back(PointInPlane{globalV, globalError, it_rh, myid});
81  }
82  }
83  }
84 
85  Road::iterator it_gh1 = temp_all_hits.begin();
86  Road::iterator it_gh2;
87 
88  patternVector_.clear();
89 
90  //look for points near wrt each other
91  // starting algorithm
92  while (it_gh1 != temp_all_hits.end() && temp_all_hits.size() >= minRoadSize_) {
93  Road temp_road;
94 
95  it_gh2 = it_gh1;
96 
97  const auto currPoint = it_gh1->globalPoint;
98  CTPPSPixelDetId currDet = CTPPSPixelDetId(it_gh1->detId);
99 
100  while (it_gh2 != temp_all_hits.end()) {
101  bool same_pot = false;
102  CTPPSPixelDetId tmpGh2Id = CTPPSPixelDetId(it_gh2->detId);
103  if (currDet.rpId() == tmpGh2Id.rpId())
104  same_pot = true;
105  const auto subtraction = currPoint - it_gh2->globalPoint;
106 
107  if (subtraction.Rho() < roadRadius_ && same_pot) {
108  temp_road.push_back(*it_gh2);
109  temp_all_hits.erase(it_gh2);
110  } else {
111  ++it_gh2;
112  }
113  }
114 
115  if (temp_road.size() >= minRoadSize_ && temp_road.size() < maxRoadSize_)
116  patternVector_.push_back(temp_road);
117  }
118  // end of algorithm
119 
120  // 2PlanePot algorithm
121 
122  for (unsigned int i = 0; i < 4; i++) {
123  if (is2PlanePot[i]) {
124  Road::iterator it_gh1_bP = temp_all_hits_2PlanePot[i].begin();
125  Road::iterator it_gh2_bP;
126 
127  while (it_gh1_bP != temp_all_hits_2PlanePot[i].end() && temp_all_hits_2PlanePot[i].size() >= 2) {
128  Road temp_road;
129 
130  it_gh2_bP = it_gh1_bP;
131 
132  const auto currPoint = it_gh1_bP->globalPoint;
133 
134  while (it_gh2_bP != temp_all_hits_2PlanePot[i].end()) {
135  const auto subtraction = currPoint - it_gh2_bP->globalPoint;
136 
137  if (subtraction.Rho() < roadRadiusBadPot_) {
138  temp_road.push_back(*it_gh2_bP);
139  temp_all_hits_2PlanePot[i].erase(it_gh2_bP);
140  } else {
141  ++it_gh2_bP;
142  }
143  }
144 
145  if (temp_road.size() == 2) { // look for isolated tracks
146  patternVector_.push_back(temp_road);
147  }
148  }
149  }
150  }
151 }
findQualityFiles.size
size
Write out results.
Definition: findQualityFiles.py:443
mps_fire.i
i
Definition: mps_fire.py:429
AlgebraicMatrix33
ROOT::Math::SMatrix< double, 3, 3, ROOT::Math::MatRepStd< double, 3, 3 > > AlgebraicMatrix33
Definition: AlgebraicROOTObjects.h:41
RPixRoadFinder::maxRoadSize_
unsigned int maxRoadSize_
Definition: RPixRoadFinder.h:44
math::Error::type
ErrorD< N >::type type
Definition: Error.h:32
RPixDetPatternFinder::Road
std::vector< PointInPlane > Road
Definition: RPixDetPatternFinder.h:38
CTPPSGeometry::localToGlobal
Vector localToGlobal(const DetGeomDesc *, const Vector &) const
Definition: CTPPSGeometry.cc:159
CTPPSGeometry::Vector
DetGeomDesc::Translation Vector
Definition: CTPPSGeometry.h:36
RPixDetPatternFinder::geometry_
const CTPPSGeometry * geometry_
Definition: RPixDetPatternFinder.h:49
RPixDetPatternFinder::hitVector_
const edm::DetSetVector< CTPPSPixelRecHit > * hitVector_
Definition: RPixDetPatternFinder.h:47
CTPPSDetId::rpId
CTPPSDetId rpId() const
Definition: CTPPSDetId.h:84
CTPPSGeometry::sensor
const DetGeomDesc * sensor(unsigned int id) const
returns geometry of a detector performs necessary checks, returns NULL if fails
Definition: CTPPSGeometry.cc:92
RPixDetPatternFinder::patternVector_
std::vector< Road > patternVector_
Definition: RPixDetPatternFinder.h:48
PointInPlane
std::pair< double, double > PointInPlane
Definition: HoughGrouping.h:29
edm::LogInfo
Log< level::Info, false > LogInfo
Definition: MessageLogger.h:131
mps_fire.end
end
Definition: mps_fire.py:242
RPixRoadFinder::roadRadiusBadPot_
double roadRadiusBadPot_
Definition: RPixRoadFinder.h:45
DetGeomDesc::rotation
const RotationMatrix & rotation() const
Definition: DetGeomDesc.h:81
RPixRoadFinder::minRoadSize_
unsigned int minRoadSize_
Definition: RPixRoadFinder.h:43
DetGeomDesc::RotationMatrix
ROOT::Math::Rotation3D RotationMatrix
Definition: DetGeomDesc.h:54

◆ 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 44 of file RPixRoadFinder.h.

Referenced by findPattern(), and RPixRoadFinder().

◆ minRoadSize_

unsigned int RPixRoadFinder::minRoadSize_
private

Definition at line 43 of file RPixRoadFinder.h.

Referenced by findPattern(), and RPixRoadFinder().

◆ roadRadius_

double RPixRoadFinder::roadRadius_
private

Definition at line 42 of file RPixRoadFinder.h.

Referenced by findPattern(), and RPixRoadFinder().

◆ roadRadiusBadPot_

double RPixRoadFinder::roadRadiusBadPot_
private

Definition at line 45 of file RPixRoadFinder.h.

Referenced by findPattern(), and RPixRoadFinder().

◆ verbosity_

int RPixRoadFinder::verbosity_
private

Definition at line 41 of file RPixRoadFinder.h.

Referenced by findPattern(), and RPixRoadFinder().

Generated for CMSSW Reference Manual by doxygen 1.8.14