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

#include <CosmicLayerPairs.h>

Inheritance diagram for CosmicLayerPairs:
SeedLayerPairs

Public Member Functions

 CosmicLayerPairs (std::string geometry)
 
void init (const SiStripRecHit2DCollection &collstereo, const SiStripRecHit2DCollection &collrphi, const SiStripMatchedRecHit2DCollection &collmatched, const edm::EventSetup &iSetup)
 
std::vector< SeedLayerPairs::LayerPair > operator() () override
 
 ~CosmicLayerPairs () override
 
- Public Member Functions inherited from SeedLayerPairs
 SeedLayerPairs ()
 
 SeedLayerPairs ()
 
virtual ~SeedLayerPairs ()
 
virtual ~SeedLayerPairs ()
 

Private Member Functions

std::vector< const TrackingRecHit * > selectTECHit (const SiStripRecHit2DCollection &collrphi, const TrackerTopology &ttopo, int side, int disk)
 
std::vector< const TrackingRecHit * > selectTECHit (const SiStripMatchedRecHit2DCollection &collmatch, const TrackerTopology &ttopo, int side, int disk)
 
std::vector< const TrackingRecHit * > selectTIBHit (const SiStripRecHit2DCollection &collrphi, const TrackerTopology &ttopo, int layer)
 
std::vector< const TrackingRecHit * > selectTIBHit (const SiStripMatchedRecHit2DCollection &collmatch, const TrackerTopology &ttopo, int layer)
 
std::vector< const TrackingRecHit * > selectTOBHit (const SiStripRecHit2DCollection &collrphi, const TrackerTopology &ttopo, int layer)
 
std::vector< const TrackingRecHit * > selectTOBHit (const SiStripMatchedRecHit2DCollection &collmatch, const TrackerTopology &ttopo, int layer)
 

Private Attributes

std::string _geometry
 
std::vector< BarrelDetLayer const * > bl
 
edm::OwnVector< LayerWithHits > CRACKLayerWithHits
 
std::vector< ForwardDetLayer const * > fneg
 
std::vector< ForwardDetLayer const * > fpos
 
edm::OwnVector< LayerWithHits > MTCCLayerWithHits
 
edm::OwnVector< LayerWithHits > TECMinusLayerWithHits
 
edm::OwnVector< LayerWithHits > TECPlusLayerWithHits
 
edm::OwnVector< LayerWithHits > TIBLayerWithHits
 
edm::OwnVector< LayerWithHits > TOBLayerWithHits
 

Additional Inherited Members

- Public Types inherited from SeedLayerPairs
typedef std::pair< const LayerWithHits *, const LayerWithHits * > LayerPair
 
typedef std::pair< const LayerWithHits *, const LayerWithHits * > LayerPair
 

Detailed Description

find all (resonable) pairs of pixel layers

Definition at line 19 of file CosmicLayerPairs.h.

Constructor & Destructor Documentation

CosmicLayerPairs::CosmicLayerPairs ( std::string  geometry)
inline

Definition at line 21 of file CosmicLayerPairs.h.

References init(), operator()(), and ~CosmicLayerPairs().

21 :_geometry(geometry){};//:isFirstCall(true){};
CosmicLayerPairs::_geometry
std::string _geometry
Definition: CosmicLayerPairs.h:38
CosmicLayerPairs::~CosmicLayerPairs ( )
override

Definition at line 252 of file CosmicLayerPairs.cc.

Referenced by CosmicLayerPairs().

252 {}

Member Function Documentation

void CosmicLayerPairs::init ( const SiStripRecHit2DCollection &  collstereo,
const SiStripRecHit2DCollection &  collrphi,
const SiStripMatchedRecHit2DCollection &  collmatched,
const edm::EventSetup &  iSetup 
)

Definition at line 256 of file CosmicLayerPairs.cc.

References _geometry, GeometricSearchTracker::barrelLayers(), bl, CRACKLayerWithHits, fneg, fpos, edm::EventSetup::get(), MTCCLayerWithHits, GeometricSearchTracker::negTecLayers(), GeometricSearchTracker::posTecLayers(), edm::OwnVector< T, P >::push_back(), selectTECHit(), selectTIBHit(), selectTOBHit(), TECMinusLayerWithHits, TECPlusLayerWithHits, TIBLayerWithHits, TOBLayerWithHits, and HiIsolationCommonParameters_cff::track.

Referenced by CosmicLayerPairs(), SeedGeneratorForCRack::init(), and SeedGeneratorForCosmics::init().

260  {
262  //_geometry=geometry;
263  //if(isFirstCall){
264  //std::cout << "in isFirtsCall" << std::endl;
265  edm::ESHandle<GeometricSearchTracker> track;
266  iSetup.get<TrackerRecoGeometryRecord>().get( track );
267  //std::cout << "about to take barrel" << std::endl;
268  bl=track->barrelLayers();
269  //std::cout << "barrel taken" << std::endl;
270  fpos=track->posTecLayers();
271  //std::cout << "pos forw taken" << std::endl;
272  fneg=track->negTecLayers();
273  //std::cout << "neg forw taken" << std::endl;
274  //isFirstCall=false;
275 
276  edm::ESHandle<TrackerTopology> httopo;
277  iSetup.get<TrackerTopologyRcd>().get(httopo);
278  const TrackerTopology& ttopo = *httopo;
279 
280  if (_geometry=="MTCC"){//we have to distinguish the MTCC and CRACK case because they have special geometries with different neumbering of layers
281  MTCCLayerWithHits.push_back(new LayerWithHits(bl[0], selectTIBHit(collrphi, ttopo, 1)));
282  MTCCLayerWithHits.push_back(new LayerWithHits(bl[1], selectTIBHit(collrphi, ttopo, 2)));
283  MTCCLayerWithHits.push_back(new LayerWithHits(bl[2], selectTOBHit(collrphi, ttopo, 1)));
284  MTCCLayerWithHits.push_back(new LayerWithHits(bl[3], selectTOBHit(collrphi, ttopo, 2)));
285  return;
286  }
287  if (_geometry=="CRACK"){
288  CRACKLayerWithHits.push_back(new LayerWithHits(bl[6], selectTOBHit(collmatched, ttopo, 7)));
289  CRACKLayerWithHits.push_back(new LayerWithHits(bl[5], selectTOBHit(collmatched, ttopo, 6)));
290  CRACKLayerWithHits.push_back(new LayerWithHits(bl[3], selectTOBHit(collmatched, ttopo, 4)));
291  CRACKLayerWithHits.push_back(new LayerWithHits(bl[2], selectTOBHit(collmatched, ttopo, 3)));
292  CRACKLayerWithHits.push_back(new LayerWithHits(bl[6], selectTOBHit(collrphi, ttopo, 7)));
293  CRACKLayerWithHits.push_back(new LayerWithHits(bl[5], selectTOBHit(collrphi, ttopo, 6)));
294  CRACKLayerWithHits.push_back(new LayerWithHits(bl[3], selectTOBHit(collrphi, ttopo, 4)));
295  CRACKLayerWithHits.push_back(new LayerWithHits(bl[2], selectTOBHit(collrphi, ttopo, 3)));
296  CRACKLayerWithHits.push_back(new LayerWithHits(bl[4], selectTOBHit(collrphi, ttopo, 5)));
297  CRACKLayerWithHits.push_back(new LayerWithHits(bl[1], selectTOBHit(collrphi, ttopo, 2)));
298  CRACKLayerWithHits.push_back(new LayerWithHits(bl[0], selectTOBHit(collrphi, ttopo, 1)));
299  CRACKLayerWithHits.push_back(new LayerWithHits(bl[4], selectTOBHit(collmatched, ttopo, 5)));
300  CRACKLayerWithHits.push_back(new LayerWithHits(bl[1], selectTOBHit(collmatched, ttopo, 2)));
301  CRACKLayerWithHits.push_back(new LayerWithHits(bl[0], selectTOBHit(collmatched, ttopo, 1)));
302  return;
303  }
304 
305  TIBLayerWithHits.push_back(new LayerWithHits(bl[3], selectTIBHit(collrphi, ttopo, 1))); //layer
306  //std::cout << "TIB 0" << std::endl;
307  TIBLayerWithHits.push_back(new LayerWithHits(bl[4], selectTIBHit(collrphi, ttopo, 2)));
308  //std::cout << "TIB 1" << std::endl;
309  TIBLayerWithHits.push_back(new LayerWithHits(bl[5], selectTIBHit(collrphi, ttopo, 3)));
310  //std::cout << "TIB 2" << std::endl;
311  TIBLayerWithHits.push_back(new LayerWithHits(bl[6], selectTIBHit(collrphi, ttopo, 4)));
312  //std::cout << "TIB 3" << std::endl;
313 
314  TOBLayerWithHits.push_back(new LayerWithHits(bl[7], selectTOBHit(collrphi, ttopo, 1)));
315  //std::cout << "TOB 0" << std::endl;
316  TOBLayerWithHits.push_back(new LayerWithHits(bl[8], selectTOBHit(collrphi, ttopo, 2)));
317  //std::cout << "TOB 1" << std::endl;
318  TOBLayerWithHits.push_back(new LayerWithHits(bl[9], selectTOBHit(collrphi, ttopo, 3)));
319  //std::cout << "TOB 2" << std::endl;
320  TOBLayerWithHits.push_back(new LayerWithHits(bl[10], selectTOBHit(collrphi, ttopo, 4)));
321  //std::cout << "TOB 3" << std::endl;
322  TOBLayerWithHits.push_back(new LayerWithHits(bl[11], selectTOBHit(collrphi, ttopo, 5)));
323  //std::cout << "TOB 4" << std::endl;
324  TOBLayerWithHits.push_back(new LayerWithHits(bl[12], selectTOBHit(collrphi, ttopo, 6)));
325  //std::cout << "TOB 5" << std::endl;
326 
327 
328  TECPlusLayerWithHits.push_back(new LayerWithHits(fpos[0], selectTECHit(collrphi, ttopo, 2, 1))); //side, disk
329  //std::cout << "wheel 0" << std::endl;
330  TECPlusLayerWithHits.push_back(new LayerWithHits(fpos[1], selectTECHit(collrphi, ttopo, 2, 2)));
331  //std::cout << "wheel 1" << std::endl;
332  TECPlusLayerWithHits.push_back(new LayerWithHits(fpos[2], selectTECHit(collrphi, ttopo, 2, 3)));
333  //std::cout << "wheel 2" << std::endl;
334  TECPlusLayerWithHits.push_back(new LayerWithHits(fpos[3], selectTECHit(collrphi, ttopo, 2, 4)));
335  //std::cout << "wheel 3" << std::endl;
336  TECPlusLayerWithHits.push_back(new LayerWithHits(fpos[4], selectTECHit(collrphi, ttopo, 2, 5)));
337  //std::cout << "wheel 4" << std::endl;
338  TECPlusLayerWithHits.push_back(new LayerWithHits(fpos[5], selectTECHit(collrphi, ttopo, 2, 6)));
339  //std::cout << "wheel 5" << std::endl;
340  TECPlusLayerWithHits.push_back(new LayerWithHits(fpos[6], selectTECHit(collrphi, ttopo, 2, 7)));
341  //std::cout << "wheel 6" << std::endl;
342  TECPlusLayerWithHits.push_back(new LayerWithHits(fpos[7], selectTECHit(collrphi, ttopo, 2, 8)));
343  //std::cout << "wheel 7" << std::endl;
344  TECPlusLayerWithHits.push_back(new LayerWithHits(fpos[8], selectTECHit(collrphi, ttopo, 2, 9)));
345  //std::cout << "wheel 8" << std::endl;
346 
347  TECMinusLayerWithHits.push_back(new LayerWithHits(fneg[0], selectTECHit(collrphi, ttopo, 1, 1))); //side, disk
348  //std::cout << "wheel 0" << std::endl;
349  TECMinusLayerWithHits.push_back(new LayerWithHits(fneg[1], selectTECHit(collrphi, ttopo, 1, 2)));
350  //std::cout << "wheel 1" << std::endl;
351  TECMinusLayerWithHits.push_back(new LayerWithHits(fneg[2], selectTECHit(collrphi, ttopo, 1, 3)));
352  //std::cout << "wheel 2" << std::endl;
353  TECMinusLayerWithHits.push_back(new LayerWithHits(fneg[3], selectTECHit(collrphi, ttopo, 1, 4)));
354  //std::cout << "wheel 3" << std::endl;
355  TECMinusLayerWithHits.push_back(new LayerWithHits(fneg[4], selectTECHit(collrphi, ttopo, 1, 5)));
356  //std::cout << "wheel 4" << std::endl;
357  TECMinusLayerWithHits.push_back(new LayerWithHits(fneg[5], selectTECHit(collrphi, ttopo, 1, 6)));
358  //std::cout << "wheel 5" << std::endl;
359  TECMinusLayerWithHits.push_back(new LayerWithHits(fneg[6], selectTECHit(collrphi, ttopo, 1, 7)));
360  //std::cout << "wheel 6" << std::endl;
361  TECMinusLayerWithHits.push_back(new LayerWithHits(fneg[7], selectTECHit(collrphi, ttopo, 1, 8)));
362  //std::cout << "wheel 7" << std::endl;
363  TECMinusLayerWithHits.push_back(new LayerWithHits(fneg[8], selectTECHit(collrphi, ttopo, 1, 9)));
364  //std::cout << "wheel 8" << std::endl;
365 }
CosmicLayerPairs::TIBLayerWithHits
edm::OwnVector< LayerWithHits > TIBLayerWithHits
Definition: CosmicLayerPairs.h:45
CosmicLayerPairs::TOBLayerWithHits
edm::OwnVector< LayerWithHits > TOBLayerWithHits
Definition: CosmicLayerPairs.h:46
edm::OwnVector::push_back
void push_back(D *&d)
Definition: OwnVector.h:290
CosmicLayerPairs::bl
std::vector< BarrelDetLayer const * > bl
Definition: CosmicLayerPairs.h:40
CosmicLayerPairs::selectTECHit
std::vector< const TrackingRecHit * > selectTECHit(const SiStripRecHit2DCollection &collrphi, const TrackerTopology &ttopo, int side, int disk)
Definition: CosmicLayerPairs.cc:367
CosmicLayerPairs::MTCCLayerWithHits
edm::OwnVector< LayerWithHits > MTCCLayerWithHits
Definition: CosmicLayerPairs.h:47
CosmicLayerPairs::TECMinusLayerWithHits
edm::OwnVector< LayerWithHits > TECMinusLayerWithHits
Definition: CosmicLayerPairs.h:44
CosmicLayerPairs::CRACKLayerWithHits
edm::OwnVector< LayerWithHits > CRACKLayerWithHits
Definition: CosmicLayerPairs.h:48
CosmicLayerPairs::selectTOBHit
std::vector< const TrackingRecHit * > selectTOBHit(const SiStripRecHit2DCollection &collrphi, const TrackerTopology &ttopo, int layer)
Definition: CosmicLayerPairs.cc:387
HiIsolationCommonParameters_cff.track
track
Definition: HiIsolationCommonParameters_cff.py:9
GeometricSearchTracker::posTecLayers
std::vector< ForwardDetLayer const * > const & posTecLayers() const
Definition: GeometricSearchTracker.h:48
edm::EventSetup::get
const T & get() const
Definition: EventSetup.h:59
CosmicLayerPairs::fneg
std::vector< ForwardDetLayer const * > fneg
Definition: CosmicLayerPairs.h:42
CosmicLayerPairs::fpos
std::vector< ForwardDetLayer const * > fpos
Definition: CosmicLayerPairs.h:41
CosmicLayerPairs::TECPlusLayerWithHits
edm::OwnVector< LayerWithHits > TECPlusLayerWithHits
Definition: CosmicLayerPairs.h:43
CosmicLayerPairs::selectTIBHit
std::vector< const TrackingRecHit * > selectTIBHit(const SiStripRecHit2DCollection &collrphi, const TrackerTopology &ttopo, int layer)
Definition: CosmicLayerPairs.cc:377
CosmicLayerPairs::_geometry
std::string _geometry
Definition: CosmicLayerPairs.h:38
GeometricSearchTracker::negTecLayers
std::vector< ForwardDetLayer const * > const & negTecLayers() const
Definition: GeometricSearchTracker.h:44
GeometricSearchTracker::barrelLayers
std::vector< BarrelDetLayer const * > const & barrelLayers() const
Definition: GeometricSearchTracker.h:32
std::vector< SeedLayerPairs::LayerPair > CosmicLayerPairs::operator() ( )
overridevirtual

Implements SeedLayerPairs.

Definition at line 15 of file CosmicLayerPairs.cc.

References _geometry, CRACKLayerWithHits, Exception, MTCCLayerWithHits, mps_fire::result, TECMinusLayerWithHits, TECPlusLayerWithHits, TIBLayerWithHits, and TOBLayerWithHits.

Referenced by CosmicLayerPairs().

16 {
17  std::vector<SeedLayerPairs::LayerPair> result;
18 
19  if (_geometry=="STANDARD"){
20 // result.push_back( SeedLayerPairs::LayerPair(&TIBLayerWithHits[1], &TIBLayerWithHits[0]));
21 // result.push_back( SeedLayerPairs::LayerPair(&TIBLayerWithHits[2], &TIBLayerWithHits[0]));
22 
23  result.push_back( SeedLayerPairs::LayerPair(&TOBLayerWithHits[3], &TOBLayerWithHits[5]));
24  result.push_back( SeedLayerPairs::LayerPair(&TOBLayerWithHits[4], &TOBLayerWithHits[5]));
25  result.push_back( SeedLayerPairs::LayerPair(&TOBLayerWithHits[3], &TOBLayerWithHits[4]));
26 
27  result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[7], &TECPlusLayerWithHits[8]));
28  result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[6], &TECPlusLayerWithHits[8]));
29 
30  result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[6], &TECPlusLayerWithHits[7]));
31  result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[5], &TECPlusLayerWithHits[7]));
32 
33  result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[5], &TECPlusLayerWithHits[6]));
34  result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[4], &TECPlusLayerWithHits[6]));
35 
36  result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[4], &TECPlusLayerWithHits[5]));
37  result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[3], &TECPlusLayerWithHits[5]));
38 
39  result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[3], &TECPlusLayerWithHits[4]));
40  result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[2], &TECPlusLayerWithHits[4]));
41 
42  result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[2], &TECPlusLayerWithHits[3]));
43  result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[1], &TECPlusLayerWithHits[3]));
44 
45  result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[1], &TECPlusLayerWithHits[2]));
46  result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[0], &TECPlusLayerWithHits[2]));
47 
48  result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[0], &TECPlusLayerWithHits[1]));
49 
50 
51  result.push_back( SeedLayerPairs::LayerPair(&TECMinusLayerWithHits[7], &TECMinusLayerWithHits[8]));
52  result.push_back( SeedLayerPairs::LayerPair(&TECMinusLayerWithHits[6], &TECMinusLayerWithHits[8]));
53 
54  result.push_back( SeedLayerPairs::LayerPair(&TECMinusLayerWithHits[6], &TECMinusLayerWithHits[7]));
55  result.push_back( SeedLayerPairs::LayerPair(&TECMinusLayerWithHits[5], &TECMinusLayerWithHits[7]));
56 
57  result.push_back( SeedLayerPairs::LayerPair(&TECMinusLayerWithHits[5], &TECMinusLayerWithHits[6]));
58  result.push_back( SeedLayerPairs::LayerPair(&TECMinusLayerWithHits[4], &TECMinusLayerWithHits[6]));
59 
60  result.push_back( SeedLayerPairs::LayerPair(&TECMinusLayerWithHits[4], &TECMinusLayerWithHits[5]));
61  result.push_back( SeedLayerPairs::LayerPair(&TECMinusLayerWithHits[3], &TECMinusLayerWithHits[5]));
62 
63  result.push_back( SeedLayerPairs::LayerPair(&TECMinusLayerWithHits[3], &TECMinusLayerWithHits[4]));
64  result.push_back( SeedLayerPairs::LayerPair(&TECMinusLayerWithHits[2], &TECMinusLayerWithHits[4]));
65 
66  result.push_back( SeedLayerPairs::LayerPair(&TECMinusLayerWithHits[2], &TECMinusLayerWithHits[3]));
67  result.push_back( SeedLayerPairs::LayerPair(&TECMinusLayerWithHits[1], &TECMinusLayerWithHits[3]));
68 
69  result.push_back( SeedLayerPairs::LayerPair(&TECMinusLayerWithHits[1], &TECMinusLayerWithHits[2]));
70  result.push_back( SeedLayerPairs::LayerPair(&TECMinusLayerWithHits[0], &TECMinusLayerWithHits[2]));
71 
72  result.push_back( SeedLayerPairs::LayerPair(&TECMinusLayerWithHits[0], &TECMinusLayerWithHits[1]));
73 
74 
75  }
76  else if(_geometry=="TECPAIRS_TOBTRIPLETS"){
77  result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[7], &TECPlusLayerWithHits[8]));
78  result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[6], &TECPlusLayerWithHits[8]));
79 
80  result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[6], &TECPlusLayerWithHits[7]));
81  result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[5], &TECPlusLayerWithHits[7]));
82 
83  result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[5], &TECPlusLayerWithHits[6]));
84  result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[4], &TECPlusLayerWithHits[6]));
85 
86  result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[4], &TECPlusLayerWithHits[5]));
87  result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[3], &TECPlusLayerWithHits[5]));
88 
89  result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[3], &TECPlusLayerWithHits[4]));
90  result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[2], &TECPlusLayerWithHits[4]));
91 
92  result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[2], &TECPlusLayerWithHits[3]));
93  result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[1], &TECPlusLayerWithHits[3]));
94 
95  result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[1], &TECPlusLayerWithHits[2]));
96  result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[0], &TECPlusLayerWithHits[2]));
97 
98  result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[0], &TECPlusLayerWithHits[1]));
99 
100 
101  result.push_back( SeedLayerPairs::LayerPair(&TECMinusLayerWithHits[7], &TECMinusLayerWithHits[8]));
102  result.push_back( SeedLayerPairs::LayerPair(&TECMinusLayerWithHits[6], &TECMinusLayerWithHits[8]));
103 
104  result.push_back( SeedLayerPairs::LayerPair(&TECMinusLayerWithHits[6], &TECMinusLayerWithHits[7]));
105  result.push_back( SeedLayerPairs::LayerPair(&TECMinusLayerWithHits[5], &TECMinusLayerWithHits[7]));
106 
107  result.push_back( SeedLayerPairs::LayerPair(&TECMinusLayerWithHits[5], &TECMinusLayerWithHits[6]));
108  result.push_back( SeedLayerPairs::LayerPair(&TECMinusLayerWithHits[4], &TECMinusLayerWithHits[6]));
109 
110  result.push_back( SeedLayerPairs::LayerPair(&TECMinusLayerWithHits[4], &TECMinusLayerWithHits[5]));
111  result.push_back( SeedLayerPairs::LayerPair(&TECMinusLayerWithHits[3], &TECMinusLayerWithHits[5]));
112 
113  result.push_back( SeedLayerPairs::LayerPair(&TECMinusLayerWithHits[3], &TECMinusLayerWithHits[4]));
114  result.push_back( SeedLayerPairs::LayerPair(&TECMinusLayerWithHits[2], &TECMinusLayerWithHits[4]));
115 
116  result.push_back( SeedLayerPairs::LayerPair(&TECMinusLayerWithHits[2], &TECMinusLayerWithHits[3]));
117  result.push_back( SeedLayerPairs::LayerPair(&TECMinusLayerWithHits[1], &TECMinusLayerWithHits[3]));
118 
119  result.push_back( SeedLayerPairs::LayerPair(&TECMinusLayerWithHits[1], &TECMinusLayerWithHits[2]));
120  result.push_back( SeedLayerPairs::LayerPair(&TECMinusLayerWithHits[0], &TECMinusLayerWithHits[2]));
121 
122  result.push_back( SeedLayerPairs::LayerPair(&TECMinusLayerWithHits[0], &TECMinusLayerWithHits[1]));
123 
124 
125  }
126  else if (_geometry=="MTCC"){
127  result.push_back( SeedLayerPairs::LayerPair(&MTCCLayerWithHits[1],&MTCCLayerWithHits[0]));
128  result.push_back( SeedLayerPairs::LayerPair(&MTCCLayerWithHits[2],&MTCCLayerWithHits[3]));
129  //IMPORTANT
130  // The seed from overlaps must be at the end
131  result.push_back( SeedLayerPairs::LayerPair(&MTCCLayerWithHits[0],&MTCCLayerWithHits[0]));
132  result.push_back( SeedLayerPairs::LayerPair(&MTCCLayerWithHits[1],&MTCCLayerWithHits[1]));
133  }
134  else if (_geometry=="CRACK"){
135  //TODO: clean all this. Now this is a random choice of layers
136  result.push_back( SeedLayerPairs::LayerPair(&CRACKLayerWithHits[1],&CRACKLayerWithHits[0]));
137  result.push_back( SeedLayerPairs::LayerPair(&CRACKLayerWithHits[2],&CRACKLayerWithHits[0]));
138  result.push_back( SeedLayerPairs::LayerPair(&CRACKLayerWithHits[2],&CRACKLayerWithHits[1]));
139  result.push_back( SeedLayerPairs::LayerPair(&CRACKLayerWithHits[3],&CRACKLayerWithHits[2]));
140  result.push_back( SeedLayerPairs::LayerPair(&CRACKLayerWithHits[5],&CRACKLayerWithHits[4]));
141  result.push_back( SeedLayerPairs::LayerPair(&CRACKLayerWithHits[6],&CRACKLayerWithHits[4]));
142  result.push_back( SeedLayerPairs::LayerPair(&CRACKLayerWithHits[6],&CRACKLayerWithHits[5]));
143  result.push_back( SeedLayerPairs::LayerPair(&CRACKLayerWithHits[7],&CRACKLayerWithHits[6]));
144  result.push_back( SeedLayerPairs::LayerPair(&CRACKLayerWithHits[8],&CRACKLayerWithHits[5]));
145  result.push_back( SeedLayerPairs::LayerPair(&CRACKLayerWithHits[9],&CRACKLayerWithHits[5]));
146  result.push_back( SeedLayerPairs::LayerPair(&CRACKLayerWithHits[10],&CRACKLayerWithHits[5]));
147  result.push_back( SeedLayerPairs::LayerPair(&CRACKLayerWithHits[11],&CRACKLayerWithHits[5]));
148  result.push_back( SeedLayerPairs::LayerPair(&CRACKLayerWithHits[12],&CRACKLayerWithHits[5]));
149  result.push_back( SeedLayerPairs::LayerPair(&CRACKLayerWithHits[13],&CRACKLayerWithHits[5]));
150  }
151  else if (_geometry=="TIBD+"){
152  result.push_back( SeedLayerPairs::LayerPair(&TIBLayerWithHits[1],&TIBLayerWithHits[0]));
153  result.push_back( SeedLayerPairs::LayerPair(&TIBLayerWithHits[2],&TIBLayerWithHits[3]));
154  result.push_back( SeedLayerPairs::LayerPair(&TIBLayerWithHits[0],&TIBLayerWithHits[0]));
155  result.push_back( SeedLayerPairs::LayerPair(&TIBLayerWithHits[1],&TIBLayerWithHits[1]));
156  }
157  else if (_geometry=="TOB") {
158  result.push_back( SeedLayerPairs::LayerPair(&TOBLayerWithHits[3], &TOBLayerWithHits[5]));
159  result.push_back( SeedLayerPairs::LayerPair(&TOBLayerWithHits[4], &TOBLayerWithHits[5]));
160  result.push_back( SeedLayerPairs::LayerPair(&TOBLayerWithHits[3], &TOBLayerWithHits[4]));
161 
162  }
163  else if(_geometry=="TIBTOB") {
164  result.push_back( SeedLayerPairs::LayerPair(&TIBLayerWithHits[1], &TIBLayerWithHits[0]));
165  result.push_back( SeedLayerPairs::LayerPair(&TIBLayerWithHits[2], &TIBLayerWithHits[0]));
166 
167  result.push_back( SeedLayerPairs::LayerPair(&TOBLayerWithHits[3], &TOBLayerWithHits[5]));
168  result.push_back( SeedLayerPairs::LayerPair(&TOBLayerWithHits[4], &TOBLayerWithHits[5]));
169  result.push_back( SeedLayerPairs::LayerPair(&TOBLayerWithHits[3], &TOBLayerWithHits[4]));
170 
171  }
172  else if (_geometry=="TEC+") {
173  result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[7], &TECPlusLayerWithHits[8]));
174  result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[6], &TECPlusLayerWithHits[8]));
175 
176  result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[6], &TECPlusLayerWithHits[7]));
177  result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[5], &TECPlusLayerWithHits[7]));
178 
179  result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[5], &TECPlusLayerWithHits[6]));
180  result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[4], &TECPlusLayerWithHits[6]));
181 
182  result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[4], &TECPlusLayerWithHits[5]));
183  result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[3], &TECPlusLayerWithHits[5]));
184 
185  result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[3], &TECPlusLayerWithHits[4]));
186  result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[2], &TECPlusLayerWithHits[4]));
187 
188  result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[2], &TECPlusLayerWithHits[3]));
189  result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[1], &TECPlusLayerWithHits[3]));
190 
191  result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[2], &TECPlusLayerWithHits[1]));
192  result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[2], &TECPlusLayerWithHits[0]));
193 
194  result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[1], &TECPlusLayerWithHits[0]));
195 
196  }
197  else if (_geometry=="CkfTIBD+"){
198  result.push_back( SeedLayerPairs::LayerPair(&TIBLayerWithHits[0], &TIBLayerWithHits[1]));
199  result.push_back( SeedLayerPairs::LayerPair(&TIBLayerWithHits[2], &TIBLayerWithHits[3]));
200  }
201  else if (_geometry=="CkfTIBTOB"){
202  result.push_back( SeedLayerPairs::LayerPair(&TIBLayerWithHits[0], &TIBLayerWithHits[1]));
203  result.push_back( SeedLayerPairs::LayerPair(&TIBLayerWithHits[0], &TIBLayerWithHits[2]));
204 
205  result.push_back( SeedLayerPairs::LayerPair(&TOBLayerWithHits[3], &TOBLayerWithHits[5]));
206  result.push_back( SeedLayerPairs::LayerPair(&TOBLayerWithHits[4], &TOBLayerWithHits[5]));
207  result.push_back( SeedLayerPairs::LayerPair(&TOBLayerWithHits[3], &TOBLayerWithHits[4]));
208  }
209  else if (_geometry=="CkfTIF3"){
210  result.push_back( SeedLayerPairs::LayerPair(&TIBLayerWithHits[0], &TIBLayerWithHits[1]));
211  result.push_back( SeedLayerPairs::LayerPair(&TIBLayerWithHits[0], &TIBLayerWithHits[2]));
212 
213  result.push_back( SeedLayerPairs::LayerPair(&TOBLayerWithHits[3], &TOBLayerWithHits[5]));
214  result.push_back( SeedLayerPairs::LayerPair(&TOBLayerWithHits[4], &TOBLayerWithHits[5]));
215  result.push_back( SeedLayerPairs::LayerPair(&TOBLayerWithHits[3], &TOBLayerWithHits[4]));
216 
217  result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[7], &TECPlusLayerWithHits[8]));
218  result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[6], &TECPlusLayerWithHits[8]));
219 
220  result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[6], &TECPlusLayerWithHits[7]));
221  result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[5], &TECPlusLayerWithHits[7]));
222 
223  result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[5], &TECPlusLayerWithHits[6]));
224  result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[4], &TECPlusLayerWithHits[6]));
225 
226  result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[4], &TECPlusLayerWithHits[5]));
227  result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[3], &TECPlusLayerWithHits[5]));
228 
229  result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[3], &TECPlusLayerWithHits[4]));
230  result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[2], &TECPlusLayerWithHits[4]));
231 
232  result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[2], &TECPlusLayerWithHits[3]));
233  result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[1], &TECPlusLayerWithHits[3]));
234 
235  result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[1], &TECPlusLayerWithHits[2]));
236  result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[0], &TECPlusLayerWithHits[2]));
237 
238  result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[0], &TECPlusLayerWithHits[1]));
239 
240  }
241  else if (_geometry=="CkfTOB"){
242  result.push_back( SeedLayerPairs::LayerPair(&TOBLayerWithHits[0], &TOBLayerWithHits[1]));
243  result.push_back( SeedLayerPairs::LayerPair(&TOBLayerWithHits[0], &TOBLayerWithHits[2]));
244 
245  result.push_back( SeedLayerPairs::LayerPair(&TOBLayerWithHits[3], &TOBLayerWithHits[4]));
246  result.push_back( SeedLayerPairs::LayerPair(&TOBLayerWithHits[4], &TOBLayerWithHits[5]));
247  result.push_back( SeedLayerPairs::LayerPair(&TOBLayerWithHits[3], &TOBLayerWithHits[5]));
248  }
249  else {throw cms::Exception("CosmicLayerPairs") << "The geometry " << _geometry << " is not implemented ";}
250  return result;
251 }
mps_fire.result
result
Definition: mps_fire.py:226
CosmicLayerPairs::TIBLayerWithHits
edm::OwnVector< LayerWithHits > TIBLayerWithHits
Definition: CosmicLayerPairs.h:45
CosmicLayerPairs::TOBLayerWithHits
edm::OwnVector< LayerWithHits > TOBLayerWithHits
Definition: CosmicLayerPairs.h:46
CosmicLayerPairs::MTCCLayerWithHits
edm::OwnVector< LayerWithHits > MTCCLayerWithHits
Definition: CosmicLayerPairs.h:47
CosmicLayerPairs::TECMinusLayerWithHits
edm::OwnVector< LayerWithHits > TECMinusLayerWithHits
Definition: CosmicLayerPairs.h:44
CosmicLayerPairs::CRACKLayerWithHits
edm::OwnVector< LayerWithHits > CRACKLayerWithHits
Definition: CosmicLayerPairs.h:48
CosmicLayerPairs::TECPlusLayerWithHits
edm::OwnVector< LayerWithHits > TECPlusLayerWithHits
Definition: CosmicLayerPairs.h:43
SeedLayerPairs::LayerPair
std::pair< const LayerWithHits *, const LayerWithHits * > LayerPair
Definition: SeedLayerPairs.h:19
CosmicLayerPairs::_geometry
std::string _geometry
Definition: CosmicLayerPairs.h:38
std::vector< const TrackingRecHit * > CosmicLayerPairs::selectTECHit ( const SiStripRecHit2DCollection &  collrphi,
const TrackerTopology &  ttopo,
int  side,
int  disk 
)
private

Definition at line 367 of file CosmicLayerPairs.cc.

References edmNew::copyDetSetRange(), and TrackerTopology::tecDetIdWheelComparator().

Referenced by init().

370  {
371  std::vector<const TrackingRecHit*> theChoosedHits;
372  edmNew::copyDetSetRange(collrphi, theChoosedHits, ttopo.tecDetIdWheelComparator(side,disk));
373  return theChoosedHits;
374 
375 }
edmNew::copyDetSetRange
void copyDetSetRange(DSTV const &dstv, std::vector< T const * > &v, std::pair< A, B > const &sel)
Definition: DetSetAlgorithm.h:42
TrackerTopology::tecDetIdWheelComparator
std::pair< DetId, SameLayerComparator > tecDetIdWheelComparator(uint32_t side, uint32_t wheel) const
Definition: TrackerTopology.h:567
std::vector< const TrackingRecHit * > CosmicLayerPairs::selectTECHit ( const SiStripMatchedRecHit2DCollection &  collmatch,
const TrackerTopology &  ttopo,
int  side,
int  disk 
)
private

Definition at line 396 of file CosmicLayerPairs.cc.

References edmNew::copyDetSetRange(), and TrackerTopology::tecDetIdWheelComparator().

399  {
400  std::vector<const TrackingRecHit*> theChoosedHits;
401  //std::cout << "in selectTECHit" << std::endl;
402  edmNew::copyDetSetRange(collmatch,theChoosedHits, ttopo.tecDetIdWheelComparator(side,disk));
403  return theChoosedHits;
404 
405 }
edmNew::copyDetSetRange
void copyDetSetRange(DSTV const &dstv, std::vector< T const * > &v, std::pair< A, B > const &sel)
Definition: DetSetAlgorithm.h:42
TrackerTopology::tecDetIdWheelComparator
std::pair< DetId, SameLayerComparator > tecDetIdWheelComparator(uint32_t side, uint32_t wheel) const
Definition: TrackerTopology.h:567
std::vector< const TrackingRecHit * > CosmicLayerPairs::selectTIBHit ( const SiStripRecHit2DCollection &  collrphi,
const TrackerTopology &  ttopo,
int  layer 
)
private

Definition at line 377 of file CosmicLayerPairs.cc.

References edmNew::copyDetSetRange(), and TrackerTopology::tibDetIdLayerComparator().

Referenced by init().

379  {
380  std::vector<const TrackingRecHit*> theChoosedHits;
381  //std::cout << "in selectTIBHit" << std::endl;
382  edmNew::copyDetSetRange(collrphi,theChoosedHits, ttopo.tibDetIdLayerComparator(layer));
383  return theChoosedHits;
384 
385 }
edmNew::copyDetSetRange
void copyDetSetRange(DSTV const &dstv, std::vector< T const * > &v, std::pair< A, B > const &sel)
Definition: DetSetAlgorithm.h:42
TrackerTopology::tibDetIdLayerComparator
std::pair< DetId, SameLayerComparator > tibDetIdLayerComparator(uint32_t layer) const
Definition: TrackerTopology.h:571
std::vector< const TrackingRecHit * > CosmicLayerPairs::selectTIBHit ( const SiStripMatchedRecHit2DCollection &  collmatch,
const TrackerTopology &  ttopo,
int  layer 
)
private

Definition at line 407 of file CosmicLayerPairs.cc.

References edmNew::copyDetSetRange(), and TrackerTopology::tibDetIdLayerComparator().

409  {
410  std::vector<const TrackingRecHit*> theChoosedHits;
411  //std::cout << "in selectTIBHit" << std::endl;
412  edmNew::copyDetSetRange(collmatch,theChoosedHits, ttopo.tibDetIdLayerComparator(layer));
413  return theChoosedHits;
414 
415 }
edmNew::copyDetSetRange
void copyDetSetRange(DSTV const &dstv, std::vector< T const * > &v, std::pair< A, B > const &sel)
Definition: DetSetAlgorithm.h:42
TrackerTopology::tibDetIdLayerComparator
std::pair< DetId, SameLayerComparator > tibDetIdLayerComparator(uint32_t layer) const
Definition: TrackerTopology.h:571
std::vector< const TrackingRecHit * > CosmicLayerPairs::selectTOBHit ( const SiStripRecHit2DCollection &  collrphi,
const TrackerTopology &  ttopo,
int  layer 
)
private

Definition at line 387 of file CosmicLayerPairs.cc.

References edmNew::copyDetSetRange(), and TrackerTopology::tobDetIdLayerComparator().

Referenced by init().

389  {
390  std::vector<const TrackingRecHit*> theChoosedHits;
391  //std::cout << "in selectTOBHit" << std::endl;
392  edmNew::copyDetSetRange(collrphi,theChoosedHits, ttopo.tobDetIdLayerComparator(layer));
393  return theChoosedHits;
394 }
edmNew::copyDetSetRange
void copyDetSetRange(DSTV const &dstv, std::vector< T const * > &v, std::pair< A, B > const &sel)
Definition: DetSetAlgorithm.h:42
TrackerTopology::tobDetIdLayerComparator
std::pair< DetId, SameLayerComparator > tobDetIdLayerComparator(uint32_t layer) const
Definition: TrackerTopology.h:579
std::vector< const TrackingRecHit * > CosmicLayerPairs::selectTOBHit ( const SiStripMatchedRecHit2DCollection &  collmatch,
const TrackerTopology &  ttopo,
int  layer 
)
private

Definition at line 417 of file CosmicLayerPairs.cc.

References edmNew::copyDetSetRange(), and TrackerTopology::tobDetIdLayerComparator().

419  {
420  std::vector<const TrackingRecHit*> theChoosedHits;
421  //std::cout << "in selectTOBHit" << std::endl;
422  edmNew::copyDetSetRange(collmatch,theChoosedHits, ttopo.tobDetIdLayerComparator(layer));
423  return theChoosedHits;
424 }
edmNew::copyDetSetRange
void copyDetSetRange(DSTV const &dstv, std::vector< T const * > &v, std::pair< A, B > const &sel)
Definition: DetSetAlgorithm.h:42
TrackerTopology::tobDetIdLayerComparator
std::pair< DetId, SameLayerComparator > tobDetIdLayerComparator(uint32_t layer) const
Definition: TrackerTopology.h:579

Member Data Documentation

std::string CosmicLayerPairs::_geometry
private

Definition at line 38 of file CosmicLayerPairs.h.

Referenced by init(), and operator()().

std::vector<BarrelDetLayer const*> CosmicLayerPairs::bl
private

Definition at line 40 of file CosmicLayerPairs.h.

Referenced by init().

edm::OwnVector<LayerWithHits> CosmicLayerPairs::CRACKLayerWithHits
private

Definition at line 48 of file CosmicLayerPairs.h.

Referenced by init(), and operator()().

std::vector<ForwardDetLayer const*> CosmicLayerPairs::fneg
private

Definition at line 42 of file CosmicLayerPairs.h.

Referenced by init().

std::vector<ForwardDetLayer const*> CosmicLayerPairs::fpos
private

Definition at line 41 of file CosmicLayerPairs.h.

Referenced by init().

edm::OwnVector<LayerWithHits> CosmicLayerPairs::MTCCLayerWithHits
private

Definition at line 47 of file CosmicLayerPairs.h.

Referenced by init(), and operator()().

edm::OwnVector<LayerWithHits> CosmicLayerPairs::TECMinusLayerWithHits
private

Definition at line 44 of file CosmicLayerPairs.h.

Referenced by init(), and operator()().

edm::OwnVector<LayerWithHits> CosmicLayerPairs::TECPlusLayerWithHits
private

Definition at line 43 of file CosmicLayerPairs.h.

Referenced by init(), and operator()().

edm::OwnVector<LayerWithHits> CosmicLayerPairs::TIBLayerWithHits
private

Definition at line 45 of file CosmicLayerPairs.h.

Referenced by init(), and operator()().

edm::OwnVector<LayerWithHits> CosmicLayerPairs::TOBLayerWithHits
private

Definition at line 46 of file CosmicLayerPairs.h.

Referenced by init(), and operator()().

Generated for CMSSW Reference Manual by doxygen 1.8.11