82 for (
unsigned idet=0; idet<whichDet.size(); ++idet ) {
84 switch ( whichDet[idet] ) {
108 if(seedingAlgo ==
"TobTecLayerPairs"){
117 if(seedingAlgo ==
"PixelLessPairs"){
119 }
else if (seedingAlgo ==
"TobTecLayerPairs"){
122 }
else if (seedingAlgo ==
"MixedTriplets"){
127 std::cout <<
"DEBUG - this should never happen" << std::endl;
139 if ( isOnDet )
break;
149 bool isAProperPair =
false;
151 unsigned int anotherSubDetId = anotherHit.
subDetId();
152 unsigned int anotherLayerNumber = anotherHit.
layerNumber();
156 ( anotherSubDetId == 1 && anotherLayerNumber == 2) ||
157 ( anotherSubDetId == 1 && anotherLayerNumber == 3) ||
158 ( anotherSubDetId == 2 && anotherLayerNumber == 1) ||
159 ( anotherSubDetId == 2 && anotherLayerNumber == 2) ) ) ||
162 ( anotherSubDetId == 1 && anotherLayerNumber == 3) ||
163 ( anotherSubDetId == 2 && anotherLayerNumber == 1) ||
164 ( anotherSubDetId == 2 && anotherLayerNumber == 2) ) ) ||
167 ( anotherSubDetId == 2 && anotherLayerNumber == 2 ) ) ||
170 ( anotherSubDetId == 6 && anotherLayerNumber == 1 ) ||
171 ( anotherSubDetId == 6 && anotherLayerNumber == 2 ) ) ) ||
174 ( anotherSubDetId == 6 && anotherLayerNumber == 2 ) ) ||
177 ( anotherSubDetId == 6 && anotherLayerNumber == 3 ) ) ||
182 (( anotherSubDetId == 3 && anotherLayerNumber == 2 ) ||
183 ( anotherSubDetId == 4 && anotherLayerNumber == 1 )) ) ||
186 ( anotherSubDetId == 4 && anotherLayerNumber == 2 ) ) ||
189 ( anotherSubDetId == 4 && anotherLayerNumber == 3 ) ) ||
192 ( anotherSubDetId == 6 && anotherLayerNumber == 1 ) ) ||
195 ( ( anotherSubDetId == 6 && anotherLayerNumber == 4 ) ||
196 ( anotherSubDetId == 6 && anotherLayerNumber == 5 )) ) ||
199 ( anotherSubDetId == 6 && anotherLayerNumber == 5 ) ) ||
204 (( anotherSubDetId == 5 && anotherLayerNumber == 2 ) ||
205 ( anotherSubDetId == 6 && anotherLayerNumber == 1 )) ) ||
208 ( anotherSubDetId == 6 && anotherLayerNumber == 2 ) ) ||
211 ( anotherSubDetId == 6 && anotherLayerNumber == 3 ) ) ||
214 ( anotherSubDetId == 6 && anotherLayerNumber == 4 ) ) ||
218 ( anotherSubDetId == 6 && anotherLayerNumber == 5 ) ) ||
221 ( anotherSubDetId == 6 && anotherLayerNumber == 6 ) ) ||
224 ( anotherSubDetId == 6 && anotherLayerNumber == 7 ) ) ;
226 return isAProperPair;
233 bool isAProperPair =
false;
235 unsigned int anotherSubDetId = anotherHit.
subDetId();
236 unsigned int anotherLayerNumber = anotherHit.
layerNumber();
240 ( anotherSubDetId == 1 && anotherLayerNumber == 2) ||
241 ( anotherSubDetId == 2 && anotherLayerNumber == 1) ) ) ||
244 ( anotherSubDetId == 1 && anotherLayerNumber == 3) ) ||
247 ( anotherSubDetId == 2 && anotherLayerNumber == 2) ) ||
250 ( anotherSubDetId == 6 && anotherLayerNumber == 2 ) );
252 return isAProperPair;
260 bool isAProperTriplet =
false;
262 unsigned int anotherSubDetId = anotherHit.
subDetId();
263 unsigned int anotherLayerNumber = anotherHit.
layerNumber();
264 unsigned int yetAnotherSubDetId = yetAnotherHit.
subDetId();
265 unsigned int yetAnotherLayerNumber = yetAnotherHit.
layerNumber();
269 ( anotherSubDetId == 1 && anotherLayerNumber == 2) && (
270 ( yetAnotherSubDetId == 1 && yetAnotherLayerNumber == 3) ||
271 ( yetAnotherSubDetId == 2 && yetAnotherLayerNumber == 1) ||
272 ( yetAnotherSubDetId == 3 && yetAnotherLayerNumber == 1) ) ) ||
275 ( anotherSubDetId == 1 && anotherLayerNumber == 3) &&
276 ( yetAnotherSubDetId == 3 && yetAnotherLayerNumber == 1) ) ||
279 ( anotherSubDetId == 1 && anotherLayerNumber == 3) &&
280 ( yetAnotherSubDetId == 3 && yetAnotherLayerNumber == 1) ) ||
283 ( anotherSubDetId == 2 && anotherLayerNumber == 1) && (
284 ( yetAnotherSubDetId == 2 && yetAnotherLayerNumber == 2) ||
285 ( yetAnotherSubDetId == 4 && yetAnotherLayerNumber == 1) ||
286 ( yetAnotherSubDetId == 4 && yetAnotherLayerNumber == 2) ) ) ||
289 ( anotherSubDetId == 2 && anotherLayerNumber == 2 ) && (
290 ( yetAnotherSubDetId == 6 && yetAnotherLayerNumber == 1 ) ||
291 ( yetAnotherSubDetId == 6 && yetAnotherLayerNumber == 2 ) ) ) ||
294 ( anotherSubDetId == 3 && anotherLayerNumber == 2 ) &&
295 ( yetAnotherSubDetId == 3 && yetAnotherLayerNumber == 3 ) );
297 return isAProperTriplet;
unsigned int theLayerNumber
unsigned int theCylinderNumber
TrackerRecHit()
Default Constructor.
unsigned int theRingNumber
bool makesATripletWith(const TrackerRecHit &anotherHit, const TrackerRecHit &yetAnotherHit) const
Check if a triplet is on the proper combination of detectors.
uint32_t rawId() const
get the raw id
const GeomDet * theGeomDet
bool makesAPairWith3rd(const TrackerRecHit &anotherHit) const
int subdetId() const
get the contents of the subdetector field (not cast into any detector's numbering enum) ...
virtual const GeomDet * idToDet(DetId) const
unsigned int layerNumber() const
The Layer Number.
const GSSiTrackerRecHit2DLocalPos * hit() const
bool makesAPairWith(const TrackerRecHit &anotherHit) const
Check if a pair is on the proper combination of detectors.
bool isOnRequestedDet(const std::vector< unsigned int > &whichDet, const std::string &seedingAlgo) const
Check if the hit is on one of the requested detector.
unsigned int subDetId() const
The subdet Id.
void init(const TrackerGeometry *theGeometry)
Initialization at construction time.
DetId geographicalId() const