CMS 3D CMS Logo

CosmicLayerPairs Class Reference

find all (resonable) pairs of pixel layers More...

#include <RecoTracker/TkHitPairs/interface/CosmicLayerPairs.h>

Inheritance diagram for CosmicLayerPairs:

SeedLayerPairs

List of all members.

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() ()
 ~CosmicLayerPairs ()

Private Member Functions

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

Private Attributes

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

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.

00021 :_geometry(geometry){};//:isFirstCall(true){};

CosmicLayerPairs::~CosmicLayerPairs (  ) 

Definition at line 251 of file CosmicLayerPairs.cc.

00251 {}

Member Function Documentation

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

Definition at line 255 of file CosmicLayerPairs.cc.

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

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

00259                                                         {
00261     //_geometry=geometry;
00262   //if(isFirstCall){
00263     //std::cout << "in isFirtsCall" << std::endl;
00264     edm::ESHandle<GeometricSearchTracker> track;
00265     iSetup.get<TrackerRecoGeometryRecord>().get( track );
00266         //std::cout << "about to take barrel" << std::endl; 
00267     bl=track->barrelLayers();
00268         //std::cout << "barrel taken" << std::endl;
00269     fpos=track->posTecLayers();
00270         //std::cout << "pos forw taken" << std::endl;
00271     fneg=track->negTecLayers();         
00272         //std::cout << "neg forw taken" << std::endl;
00273   //isFirstCall=false;
00274         
00275     if (_geometry=="MTCC"){//we have to distinguish the MTCC and CRACK case because they have special geometries with different neumbering of layers
00276         MTCCLayerWithHits.push_back(new LayerWithHits(bl[0], selectTIBHit(collrphi, 1)));
00277         MTCCLayerWithHits.push_back(new LayerWithHits(bl[1], selectTIBHit(collrphi, 2)));
00278         MTCCLayerWithHits.push_back(new LayerWithHits(bl[2], selectTOBHit(collrphi, 1)));
00279         MTCCLayerWithHits.push_back(new LayerWithHits(bl[3], selectTOBHit(collrphi, 2)));
00280         return;
00281     }
00282     if (_geometry=="CRACK"){
00283         CRACKLayerWithHits.push_back(new LayerWithHits(bl[6], selectTOBHit(collmatched, 7)));
00284         CRACKLayerWithHits.push_back(new LayerWithHits(bl[5], selectTOBHit(collmatched, 6)));
00285         CRACKLayerWithHits.push_back(new LayerWithHits(bl[3], selectTOBHit(collmatched, 4)));
00286         CRACKLayerWithHits.push_back(new LayerWithHits(bl[2], selectTOBHit(collmatched, 3)));
00287         CRACKLayerWithHits.push_back(new LayerWithHits(bl[6], selectTOBHit(collrphi, 7)));
00288         CRACKLayerWithHits.push_back(new LayerWithHits(bl[5], selectTOBHit(collrphi, 6)));
00289         CRACKLayerWithHits.push_back(new LayerWithHits(bl[3], selectTOBHit(collrphi, 4)));
00290         CRACKLayerWithHits.push_back(new LayerWithHits(bl[2], selectTOBHit(collrphi, 3)));
00291         CRACKLayerWithHits.push_back(new LayerWithHits(bl[4], selectTOBHit(collrphi, 5)));
00292         CRACKLayerWithHits.push_back(new LayerWithHits(bl[1], selectTOBHit(collrphi, 2)));
00293         CRACKLayerWithHits.push_back(new LayerWithHits(bl[0], selectTOBHit(collrphi, 1)));
00294         CRACKLayerWithHits.push_back(new LayerWithHits(bl[4], selectTOBHit(collmatched, 5)));
00295         CRACKLayerWithHits.push_back(new LayerWithHits(bl[1], selectTOBHit(collmatched, 2)));
00296         CRACKLayerWithHits.push_back(new LayerWithHits(bl[0], selectTOBHit(collmatched, 1)));
00297         return;
00298     }   
00299         
00300     TIBLayerWithHits.push_back(new LayerWithHits(bl[3], selectTIBHit(collrphi, 1)));  //layer
00301         //std::cout << "TIB 0" << std::endl;
00302     TIBLayerWithHits.push_back(new LayerWithHits(bl[4], selectTIBHit(collrphi, 2)));
00303         //std::cout << "TIB 1" << std::endl;
00304     TIBLayerWithHits.push_back(new LayerWithHits(bl[5], selectTIBHit(collrphi, 3)));
00305         //std::cout << "TIB 2" << std::endl;
00306     TIBLayerWithHits.push_back(new LayerWithHits(bl[6], selectTIBHit(collrphi, 4)));
00307         //std::cout << "TIB 3" << std::endl;
00308 
00309     TOBLayerWithHits.push_back(new LayerWithHits(bl[7], selectTOBHit(collrphi, 1)));
00310         //std::cout << "TOB 0" << std::endl;
00311     TOBLayerWithHits.push_back(new LayerWithHits(bl[8], selectTOBHit(collrphi, 2)));
00312         //std::cout << "TOB 1" << std::endl;
00313     TOBLayerWithHits.push_back(new LayerWithHits(bl[9], selectTOBHit(collrphi, 3)));
00314         //std::cout << "TOB 2" << std::endl;
00315     TOBLayerWithHits.push_back(new LayerWithHits(bl[10], selectTOBHit(collrphi, 4)));
00316         //std::cout << "TOB 3" << std::endl;
00317     TOBLayerWithHits.push_back(new LayerWithHits(bl[11], selectTOBHit(collrphi, 5)));
00318         //std::cout << "TOB 4" << std::endl;
00319     TOBLayerWithHits.push_back(new LayerWithHits(bl[12], selectTOBHit(collrphi, 6)));
00320         //std::cout << "TOB 5" << std::endl;
00321 
00322 
00323     TECPlusLayerWithHits.push_back(new LayerWithHits(fpos[0], selectTECHit(collrphi, 2, 1)));  //side, disk
00324         //std::cout << "wheel 0" << std::endl;
00325     TECPlusLayerWithHits.push_back(new LayerWithHits(fpos[1], selectTECHit(collrphi, 2, 2)));  
00326         //std::cout << "wheel 1" << std::endl;
00327     TECPlusLayerWithHits.push_back(new LayerWithHits(fpos[2], selectTECHit(collrphi, 2, 3)));  
00328         //std::cout << "wheel 2" << std::endl;
00329     TECPlusLayerWithHits.push_back(new LayerWithHits(fpos[3], selectTECHit(collrphi, 2, 4)));  
00330         //std::cout << "wheel 3" << std::endl;
00331     TECPlusLayerWithHits.push_back(new LayerWithHits(fpos[4], selectTECHit(collrphi, 2, 5)));  
00332         //std::cout << "wheel 4" << std::endl;
00333     TECPlusLayerWithHits.push_back(new LayerWithHits(fpos[5], selectTECHit(collrphi, 2, 6)));  
00334         //std::cout << "wheel 5" << std::endl;
00335     TECPlusLayerWithHits.push_back(new LayerWithHits(fpos[6], selectTECHit(collrphi, 2, 7)));  
00336         //std::cout << "wheel 6" << std::endl;
00337     TECPlusLayerWithHits.push_back(new LayerWithHits(fpos[7], selectTECHit(collrphi, 2, 8)));  
00338         //std::cout << "wheel 7" << std::endl;
00339     TECPlusLayerWithHits.push_back(new LayerWithHits(fpos[8], selectTECHit(collrphi, 2, 9)));  
00340         //std::cout << "wheel 8" << std::endl;
00341 
00342     TECMinusLayerWithHits.push_back(new LayerWithHits(fneg[0], selectTECHit(collrphi, 1, 1)));  //side, disk
00343         //std::cout << "wheel 0" << std::endl;
00344     TECMinusLayerWithHits.push_back(new LayerWithHits(fneg[1], selectTECHit(collrphi, 1, 2)));
00345         //std::cout << "wheel 1" << std::endl;
00346     TECMinusLayerWithHits.push_back(new LayerWithHits(fneg[2], selectTECHit(collrphi, 1, 3)));
00347         //std::cout << "wheel 2" << std::endl;
00348     TECMinusLayerWithHits.push_back(new LayerWithHits(fneg[3], selectTECHit(collrphi, 1, 4)));
00349         //std::cout << "wheel 3" << std::endl;
00350     TECMinusLayerWithHits.push_back(new LayerWithHits(fneg[4], selectTECHit(collrphi, 1, 5)));
00351         //std::cout << "wheel 4" << std::endl;
00352     TECMinusLayerWithHits.push_back(new LayerWithHits(fneg[5], selectTECHit(collrphi, 1, 6)));
00353         //std::cout << "wheel 5" << std::endl;
00354     TECMinusLayerWithHits.push_back(new LayerWithHits(fneg[6], selectTECHit(collrphi, 1, 7)));
00355         //std::cout << "wheel 6" << std::endl;
00356     TECMinusLayerWithHits.push_back(new LayerWithHits(fneg[7], selectTECHit(collrphi, 1, 8)));
00357         //std::cout << "wheel 7" << std::endl;
00358     TECMinusLayerWithHits.push_back(new LayerWithHits(fneg[8], selectTECHit(collrphi, 1, 9)));
00359         //std::cout << "wheel 8" << std::endl;
00360 }

std::vector< SeedLayerPairs::LayerPair > CosmicLayerPairs::operator() ( void   )  [virtual]

Implements SeedLayerPairs.

Definition at line 14 of file CosmicLayerPairs.cc.

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

00015 {
00016   std::vector<SeedLayerPairs::LayerPair> result;
00017 
00018   if (_geometry=="STANDARD"){
00019 //      result.push_back( SeedLayerPairs::LayerPair(&TIBLayerWithHits[1], &TIBLayerWithHits[0]));
00020 //      result.push_back( SeedLayerPairs::LayerPair(&TIBLayerWithHits[2], &TIBLayerWithHits[0]));
00021 
00022     result.push_back( SeedLayerPairs::LayerPair(&TOBLayerWithHits[3], &TOBLayerWithHits[5]));
00023     result.push_back( SeedLayerPairs::LayerPair(&TOBLayerWithHits[4], &TOBLayerWithHits[5]));
00024     result.push_back( SeedLayerPairs::LayerPair(&TOBLayerWithHits[3], &TOBLayerWithHits[4]));
00025 
00026     result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[7], &TECPlusLayerWithHits[8]));
00027     result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[6], &TECPlusLayerWithHits[8]));
00028 
00029     result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[6], &TECPlusLayerWithHits[7]));
00030     result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[5], &TECPlusLayerWithHits[7]));
00031 
00032     result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[5], &TECPlusLayerWithHits[6]));
00033     result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[4], &TECPlusLayerWithHits[6]));
00034 
00035     result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[4], &TECPlusLayerWithHits[5]));
00036     result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[3], &TECPlusLayerWithHits[5]));
00037 
00038     result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[3], &TECPlusLayerWithHits[4]));
00039     result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[2], &TECPlusLayerWithHits[4]));
00040 
00041     result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[2], &TECPlusLayerWithHits[3]));
00042     result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[1], &TECPlusLayerWithHits[3]));
00043 
00044     result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[1], &TECPlusLayerWithHits[2]));
00045     result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[0], &TECPlusLayerWithHits[2]));
00046 
00047     result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[0], &TECPlusLayerWithHits[1]));
00048 
00049 
00050     result.push_back( SeedLayerPairs::LayerPair(&TECMinusLayerWithHits[7], &TECMinusLayerWithHits[8]));
00051     result.push_back( SeedLayerPairs::LayerPair(&TECMinusLayerWithHits[6], &TECMinusLayerWithHits[8]));
00052 
00053     result.push_back( SeedLayerPairs::LayerPair(&TECMinusLayerWithHits[6], &TECMinusLayerWithHits[7]));
00054     result.push_back( SeedLayerPairs::LayerPair(&TECMinusLayerWithHits[5], &TECMinusLayerWithHits[7]));
00055 
00056     result.push_back( SeedLayerPairs::LayerPair(&TECMinusLayerWithHits[5], &TECMinusLayerWithHits[6]));
00057     result.push_back( SeedLayerPairs::LayerPair(&TECMinusLayerWithHits[4], &TECMinusLayerWithHits[6]));
00058 
00059     result.push_back( SeedLayerPairs::LayerPair(&TECMinusLayerWithHits[4], &TECMinusLayerWithHits[5]));
00060     result.push_back( SeedLayerPairs::LayerPair(&TECMinusLayerWithHits[3], &TECMinusLayerWithHits[5]));
00061 
00062     result.push_back( SeedLayerPairs::LayerPair(&TECMinusLayerWithHits[3], &TECMinusLayerWithHits[4]));
00063     result.push_back( SeedLayerPairs::LayerPair(&TECMinusLayerWithHits[2], &TECMinusLayerWithHits[4]));
00064 
00065     result.push_back( SeedLayerPairs::LayerPair(&TECMinusLayerWithHits[2], &TECMinusLayerWithHits[3]));
00066     result.push_back( SeedLayerPairs::LayerPair(&TECMinusLayerWithHits[1], &TECMinusLayerWithHits[3]));
00067 
00068     result.push_back( SeedLayerPairs::LayerPair(&TECMinusLayerWithHits[1], &TECMinusLayerWithHits[2]));
00069     result.push_back( SeedLayerPairs::LayerPair(&TECMinusLayerWithHits[0], &TECMinusLayerWithHits[2]));
00070 
00071     result.push_back( SeedLayerPairs::LayerPair(&TECMinusLayerWithHits[0], &TECMinusLayerWithHits[1]));
00072 
00073 
00074   } 
00075   else if(_geometry=="TECPAIRS_TOBTRIPLETS"){
00076     result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[7], &TECPlusLayerWithHits[8]));
00077     result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[6], &TECPlusLayerWithHits[8]));
00078     
00079     result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[6], &TECPlusLayerWithHits[7]));
00080     result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[5], &TECPlusLayerWithHits[7]));
00081     
00082     result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[5], &TECPlusLayerWithHits[6]));
00083     result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[4], &TECPlusLayerWithHits[6]));
00084     
00085     result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[4], &TECPlusLayerWithHits[5]));
00086     result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[3], &TECPlusLayerWithHits[5]));
00087 
00088     result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[3], &TECPlusLayerWithHits[4]));
00089     result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[2], &TECPlusLayerWithHits[4]));
00090 
00091     result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[2], &TECPlusLayerWithHits[3]));
00092     result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[1], &TECPlusLayerWithHits[3]));
00093 
00094     result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[1], &TECPlusLayerWithHits[2]));
00095     result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[0], &TECPlusLayerWithHits[2]));
00096 
00097     result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[0], &TECPlusLayerWithHits[1]));
00098 
00099 
00100     result.push_back( SeedLayerPairs::LayerPair(&TECMinusLayerWithHits[7], &TECMinusLayerWithHits[8]));
00101     result.push_back( SeedLayerPairs::LayerPair(&TECMinusLayerWithHits[6], &TECMinusLayerWithHits[8]));
00102 
00103     result.push_back( SeedLayerPairs::LayerPair(&TECMinusLayerWithHits[6], &TECMinusLayerWithHits[7]));
00104     result.push_back( SeedLayerPairs::LayerPair(&TECMinusLayerWithHits[5], &TECMinusLayerWithHits[7]));
00105 
00106     result.push_back( SeedLayerPairs::LayerPair(&TECMinusLayerWithHits[5], &TECMinusLayerWithHits[6]));
00107     result.push_back( SeedLayerPairs::LayerPair(&TECMinusLayerWithHits[4], &TECMinusLayerWithHits[6]));
00108 
00109     result.push_back( SeedLayerPairs::LayerPair(&TECMinusLayerWithHits[4], &TECMinusLayerWithHits[5]));
00110     result.push_back( SeedLayerPairs::LayerPair(&TECMinusLayerWithHits[3], &TECMinusLayerWithHits[5]));
00111 
00112     result.push_back( SeedLayerPairs::LayerPair(&TECMinusLayerWithHits[3], &TECMinusLayerWithHits[4]));
00113     result.push_back( SeedLayerPairs::LayerPair(&TECMinusLayerWithHits[2], &TECMinusLayerWithHits[4]));
00114 
00115     result.push_back( SeedLayerPairs::LayerPair(&TECMinusLayerWithHits[2], &TECMinusLayerWithHits[3]));
00116     result.push_back( SeedLayerPairs::LayerPair(&TECMinusLayerWithHits[1], &TECMinusLayerWithHits[3]));
00117 
00118     result.push_back( SeedLayerPairs::LayerPair(&TECMinusLayerWithHits[1], &TECMinusLayerWithHits[2]));
00119     result.push_back( SeedLayerPairs::LayerPair(&TECMinusLayerWithHits[0], &TECMinusLayerWithHits[2]));
00120 
00121     result.push_back( SeedLayerPairs::LayerPair(&TECMinusLayerWithHits[0], &TECMinusLayerWithHits[1]));
00122     
00123 
00124   }
00125   else if (_geometry=="MTCC"){
00126     result.push_back( SeedLayerPairs::LayerPair(&MTCCLayerWithHits[1],&MTCCLayerWithHits[0]));
00127     result.push_back( SeedLayerPairs::LayerPair(&MTCCLayerWithHits[2],&MTCCLayerWithHits[3]));
00128     //IMPORTANT
00129     // The seed from overlaps must be at the end
00130     result.push_back( SeedLayerPairs::LayerPair(&MTCCLayerWithHits[0],&MTCCLayerWithHits[0]));
00131     result.push_back( SeedLayerPairs::LayerPair(&MTCCLayerWithHits[1],&MTCCLayerWithHits[1]));  
00132   } 
00133   else if (_geometry=="CRACK"){
00134     //TODO: clean all this. Now this is a random choice of layers
00135     result.push_back( SeedLayerPairs::LayerPair(&CRACKLayerWithHits[1],&CRACKLayerWithHits[0]));
00136     result.push_back( SeedLayerPairs::LayerPair(&CRACKLayerWithHits[2],&CRACKLayerWithHits[0]));
00137     result.push_back( SeedLayerPairs::LayerPair(&CRACKLayerWithHits[2],&CRACKLayerWithHits[1]));
00138     result.push_back( SeedLayerPairs::LayerPair(&CRACKLayerWithHits[3],&CRACKLayerWithHits[2]));
00139     result.push_back( SeedLayerPairs::LayerPair(&CRACKLayerWithHits[5],&CRACKLayerWithHits[4]));
00140     result.push_back( SeedLayerPairs::LayerPair(&CRACKLayerWithHits[6],&CRACKLayerWithHits[4]));
00141     result.push_back( SeedLayerPairs::LayerPair(&CRACKLayerWithHits[6],&CRACKLayerWithHits[5]));
00142     result.push_back( SeedLayerPairs::LayerPair(&CRACKLayerWithHits[7],&CRACKLayerWithHits[6]));
00143     result.push_back( SeedLayerPairs::LayerPair(&CRACKLayerWithHits[8],&CRACKLayerWithHits[5]));
00144     result.push_back( SeedLayerPairs::LayerPair(&CRACKLayerWithHits[9],&CRACKLayerWithHits[5]));
00145     result.push_back( SeedLayerPairs::LayerPair(&CRACKLayerWithHits[10],&CRACKLayerWithHits[5]));
00146     result.push_back( SeedLayerPairs::LayerPair(&CRACKLayerWithHits[11],&CRACKLayerWithHits[5]));
00147     result.push_back( SeedLayerPairs::LayerPair(&CRACKLayerWithHits[12],&CRACKLayerWithHits[5]));
00148     result.push_back( SeedLayerPairs::LayerPair(&CRACKLayerWithHits[13],&CRACKLayerWithHits[5]));
00149   } 
00150   else if (_geometry=="TIBD+"){
00151     result.push_back( SeedLayerPairs::LayerPair(&TIBLayerWithHits[1],&TIBLayerWithHits[0]));
00152     result.push_back( SeedLayerPairs::LayerPair(&TIBLayerWithHits[2],&TIBLayerWithHits[3]));
00153     result.push_back( SeedLayerPairs::LayerPair(&TIBLayerWithHits[0],&TIBLayerWithHits[0]));
00154     result.push_back( SeedLayerPairs::LayerPair(&TIBLayerWithHits[1],&TIBLayerWithHits[1]));
00155   } 
00156   else if (_geometry=="TOB") {
00157     result.push_back( SeedLayerPairs::LayerPair(&TOBLayerWithHits[3], &TOBLayerWithHits[5]));
00158     result.push_back( SeedLayerPairs::LayerPair(&TOBLayerWithHits[4], &TOBLayerWithHits[5]));
00159     result.push_back( SeedLayerPairs::LayerPair(&TOBLayerWithHits[3], &TOBLayerWithHits[4]));
00160     
00161   } 
00162   else if(_geometry=="TIBTOB") {
00163     result.push_back( SeedLayerPairs::LayerPair(&TIBLayerWithHits[1], &TIBLayerWithHits[0]));
00164     result.push_back( SeedLayerPairs::LayerPair(&TIBLayerWithHits[2], &TIBLayerWithHits[0]));
00165 
00166     result.push_back( SeedLayerPairs::LayerPair(&TOBLayerWithHits[3], &TOBLayerWithHits[5]));
00167     result.push_back( SeedLayerPairs::LayerPair(&TOBLayerWithHits[4], &TOBLayerWithHits[5]));
00168     result.push_back( SeedLayerPairs::LayerPair(&TOBLayerWithHits[3], &TOBLayerWithHits[4]));
00169  
00170   }
00171   else if (_geometry=="TEC+") {
00172     result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[7], &TECPlusLayerWithHits[8]));
00173     result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[6], &TECPlusLayerWithHits[8]));
00174     
00175     result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[6], &TECPlusLayerWithHits[7]));
00176     result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[5], &TECPlusLayerWithHits[7]));
00177     
00178     result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[5], &TECPlusLayerWithHits[6]));
00179     result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[4], &TECPlusLayerWithHits[6]));
00180     
00181     result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[4], &TECPlusLayerWithHits[5]));
00182     result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[3], &TECPlusLayerWithHits[5]));
00183     
00184     result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[3], &TECPlusLayerWithHits[4]));
00185     result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[2], &TECPlusLayerWithHits[4]));
00186     
00187     result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[2], &TECPlusLayerWithHits[3]));
00188     result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[1], &TECPlusLayerWithHits[3]));
00189     
00190     result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[2], &TECPlusLayerWithHits[1]));
00191     result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[2], &TECPlusLayerWithHits[0]));
00192     
00193     result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[1], &TECPlusLayerWithHits[0]));
00194     
00195   }
00196   else if (_geometry=="CkfTIBD+"){
00197     result.push_back( SeedLayerPairs::LayerPair(&TIBLayerWithHits[0], &TIBLayerWithHits[1]));
00198     result.push_back( SeedLayerPairs::LayerPair(&TIBLayerWithHits[2], &TIBLayerWithHits[3]));
00199   } 
00200   else if (_geometry=="CkfTIBTOB"){
00201     result.push_back( SeedLayerPairs::LayerPair(&TIBLayerWithHits[0], &TIBLayerWithHits[1]));
00202     result.push_back( SeedLayerPairs::LayerPair(&TIBLayerWithHits[0], &TIBLayerWithHits[2]));
00203 
00204     result.push_back( SeedLayerPairs::LayerPair(&TOBLayerWithHits[3], &TOBLayerWithHits[5]));
00205     result.push_back( SeedLayerPairs::LayerPair(&TOBLayerWithHits[4], &TOBLayerWithHits[5]));
00206     result.push_back( SeedLayerPairs::LayerPair(&TOBLayerWithHits[3], &TOBLayerWithHits[4]));
00207   }
00208   else if (_geometry=="CkfTIF3"){
00209     result.push_back( SeedLayerPairs::LayerPair(&TIBLayerWithHits[0], &TIBLayerWithHits[1]));
00210     result.push_back( SeedLayerPairs::LayerPair(&TIBLayerWithHits[0], &TIBLayerWithHits[2]));
00211     
00212     result.push_back( SeedLayerPairs::LayerPair(&TOBLayerWithHits[3], &TOBLayerWithHits[5]));
00213     result.push_back( SeedLayerPairs::LayerPair(&TOBLayerWithHits[4], &TOBLayerWithHits[5]));
00214     result.push_back( SeedLayerPairs::LayerPair(&TOBLayerWithHits[3], &TOBLayerWithHits[4]));
00215 
00216     result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[7], &TECPlusLayerWithHits[8]));   
00217     result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[6], &TECPlusLayerWithHits[8]));   
00218 
00219     result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[6], &TECPlusLayerWithHits[7]));   
00220     result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[5], &TECPlusLayerWithHits[7]));
00221         
00222     result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[5], &TECPlusLayerWithHits[6]));   
00223     result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[4], &TECPlusLayerWithHits[6]));   
00224 
00225     result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[4], &TECPlusLayerWithHits[5]));   
00226     result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[3], &TECPlusLayerWithHits[5]));   
00227 
00228     result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[3], &TECPlusLayerWithHits[4]));   
00229     result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[2], &TECPlusLayerWithHits[4]));   
00230 
00231     result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[2], &TECPlusLayerWithHits[3]));   
00232     result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[1], &TECPlusLayerWithHits[3]));   
00233 
00234     result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[1], &TECPlusLayerWithHits[2]));   
00235     result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[0], &TECPlusLayerWithHits[2]));   
00236 
00237     result.push_back( SeedLayerPairs::LayerPair(&TECPlusLayerWithHits[0], &TECPlusLayerWithHits[1]));   
00238 
00239   } 
00240   else if (_geometry=="CkfTOB"){
00241     result.push_back( SeedLayerPairs::LayerPair(&TOBLayerWithHits[0], &TOBLayerWithHits[1]));
00242     result.push_back( SeedLayerPairs::LayerPair(&TOBLayerWithHits[0], &TOBLayerWithHits[2]));
00243 
00244     result.push_back( SeedLayerPairs::LayerPair(&TOBLayerWithHits[3], &TOBLayerWithHits[4]));   
00245     result.push_back( SeedLayerPairs::LayerPair(&TOBLayerWithHits[4], &TOBLayerWithHits[5]));
00246     result.push_back( SeedLayerPairs::LayerPair(&TOBLayerWithHits[3], &TOBLayerWithHits[5]));
00247   } 
00248   else {throw cms::Exception("CosmicLayerPairs") << "The geometry " << _geometry << " is not implemented ";}
00249   return result;
00250 }

std::vector<const TrackingRecHit*> CosmicLayerPairs::selectTECHit ( const SiStripMatchedRecHit2DCollection &  collmatch,
int  side,
int  disk 
) [private]

std::vector< const TrackingRecHit * > CosmicLayerPairs::selectTECHit ( const SiStripRecHit2DCollection &  collrphi,
int  side,
int  disk 
) [private]

Definition at line 362 of file CosmicLayerPairs.cc.

References edm::RangeMap< ID, C, P >::get(), it, range, and TrackerLayerIdAccessor::stripTECDisk().

Referenced by init().

00364                                                                          {
00365         std::vector<const TrackingRecHit*> theChoosedHits;
00366         TrackerLayerIdAccessor acc;
00367         //std::cout << "in selectTECHit" << std::endl;  
00368         SiStripRecHit2DCollection::range range = collrphi.get(acc.stripTECDisk(side,disk));
00369         for(SiStripRecHit2DCollection::const_iterator it = range.first; it != range.second; it++){
00370                 theChoosedHits.push_back( &(*it) );     
00371         }
00372         return theChoosedHits;
00373         
00374 }

std::vector<const TrackingRecHit*> CosmicLayerPairs::selectTIBHit ( const SiStripMatchedRecHit2DCollection &  collmatch,
int  layer 
) [private]

std::vector< const TrackingRecHit * > CosmicLayerPairs::selectTIBHit ( const SiStripRecHit2DCollection &  collrphi,
int  layer 
) [private]

Definition at line 376 of file CosmicLayerPairs.cc.

References edm::RangeMap< ID, C, P >::get(), it, range, and TrackerLayerIdAccessor::stripTIBLayer().

Referenced by init().

00377                                                                           {
00378         std::vector<const TrackingRecHit*> theChoosedHits;
00379         TrackerLayerIdAccessor acc; 
00380         //std::cout << "in selectTIBHit" << std::endl;
00381         SiStripRecHit2DCollection::range range = collrphi.get(acc.stripTIBLayer(layer));
00382         for(SiStripRecHit2DCollection::const_iterator it = range.first; it != range.second; it++){
00383                 theChoosedHits.push_back( &(*it) );
00384         }
00385         return theChoosedHits;
00386 
00387 }

std::vector<const TrackingRecHit*> CosmicLayerPairs::selectTOBHit ( const SiStripMatchedRecHit2DCollection &  collmatch,
int  layer 
) [private]

std::vector< const TrackingRecHit * > CosmicLayerPairs::selectTOBHit ( const SiStripRecHit2DCollection &  collrphi,
int  layer 
) [private]

Definition at line 389 of file CosmicLayerPairs.cc.

References edm::RangeMap< ID, C, P >::get(), it, range, and TrackerLayerIdAccessor::stripTOBLayer().

Referenced by init().

00390                                                                           {
00391         std::vector<const TrackingRecHit*> theChoosedHits;
00392         TrackerLayerIdAccessor acc;
00393         //std::cout << "in selectTOBHit" << std::endl;
00394         SiStripRecHit2DCollection::range range = collrphi.get(acc.stripTOBLayer(layer));
00395         for(SiStripRecHit2DCollection::const_iterator it = range.first; it != range.second; it++){
00396                 theChoosedHits.push_back( &(*it) );
00397         }
00398         return theChoosedHits;
00399 }

Member Data Documentation

std::string CosmicLayerPairs::_geometry [private]

Definition at line 38 of file CosmicLayerPairs.h.

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

std::vector<BarrelDetLayer*> 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*> CosmicLayerPairs::fneg [private]

Definition at line 42 of file CosmicLayerPairs.h.

Referenced by init().

std::vector<ForwardDetLayer*> 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()().

The documentation for this class was generated from the following files:
Generated on Tue Jun 9 18:16:52 2009 for CMSSW by  doxygen 1.5.4