CMS 3D CMS Logo

 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
List of all members | Public Types | Public Member Functions | Private Attributes
CSCTFTrackBuilder Class Reference

#include <CSCTFTrackBuilder.h>

Public Types

enum  { nEndcaps = 2, nSectors = 6 }
 

Public Member Functions

void buildTracks (const CSCCorrelatedLCTDigiCollection *, const CSCTriggerContainer< csctf::TrackStub > *, L1CSCTrackCollection *, CSCTriggerContainer< csctf::TrackStub > *)
 
 CSCTFTrackBuilder (const edm::ParameterSet &pset, bool TMB07, const L1MuTriggerScales *scales, const L1MuTriggerPtScale *ptScale)
 
void initialize (const edm::EventSetup &c)
 
 ~CSCTFTrackBuilder ()
 

Private Attributes

int m_maxBX
 
int m_minBX
 
CSCTFSectorProcessormy_SPs [nEndcaps][nSectors]
 

Detailed Description

Definition at line 19 of file CSCTFTrackBuilder.h.

Member Enumeration Documentation

anonymous enum
Enumerator
nEndcaps 
nSectors 

Definition at line 25 of file CSCTFTrackBuilder.h.

Constructor & Destructor Documentation

CSCTFTrackBuilder::CSCTFTrackBuilder ( const edm::ParameterSet pset,
bool  TMB07,
const L1MuTriggerScales scales,
const L1MuTriggerPtScale ptScale 
)

Definition at line 18 of file CSCTFTrackBuilder.cc.

References edm::ParameterSet::getParameter(), m_maxBX, m_minBX, CSCDetId::maxEndcapId(), CSCTriggerNumbering::maxTriggerSectorId(), CSCDetId::minEndcapId(), CSCTriggerNumbering::minTriggerSectorId(), my_SPs, and asciidump::s.

20  {
21  m_minBX = pset.getParameter<int>("MinBX");
22  m_maxBX = pset.getParameter<int>("MaxBX");
23 
24  for(int e = CSCDetId::minEndcapId(); e <= CSCDetId::maxEndcapId(); ++e)
25  {
28  {
29  // All SPs work with the same configuration (impossible to make it more exclusive in this framework)
30  my_SPs[e-1][s-1] = new CSCTFSectorProcessor(e, s, pset, TMB07,
31  scales, ptScale);
32  }
33  }
34 }
T getParameter(std::string const &) const
static int minEndcapId()
Definition: CSCDetId.h:238
static int maxTriggerSectorId()
static int minTriggerSectorId()
static int maxEndcapId()
Definition: CSCDetId.h:239
string s
Definition: asciidump.py:422
CSCTFSectorProcessor * my_SPs[nEndcaps][nSectors]
CSCTFTrackBuilder::~CSCTFTrackBuilder ( )

Definition at line 47 of file CSCTFTrackBuilder.cc.

References CSCDetId::maxEndcapId(), CSCTriggerNumbering::maxTriggerSectorId(), CSCDetId::minEndcapId(), CSCTriggerNumbering::minTriggerSectorId(), my_SPs, NULL, and asciidump::s.

48 {
49  for(int e = CSCDetId::minEndcapId(); e <= CSCDetId::maxEndcapId(); ++e)
50  {
53  {
54  delete my_SPs[e-1][s-1];
55  my_SPs[e-1][s-1] = NULL;
56  }
57  }
58 }
static int minEndcapId()
Definition: CSCDetId.h:238
static int maxTriggerSectorId()
#define NULL
Definition: scimark2.h:8
static int minTriggerSectorId()
static int maxEndcapId()
Definition: CSCDetId.h:239
string s
Definition: asciidump.py:422
CSCTFSectorProcessor * my_SPs[nEndcaps][nSectors]

Member Function Documentation

void CSCTFTrackBuilder::buildTracks ( const CSCCorrelatedLCTDigiCollection lcts,
const CSCTriggerContainer< csctf::TrackStub > *  dtstubss,
L1CSCTrackCollection trkcoll,
CSCTriggerContainer< csctf::TrackStub > *  stubs_to_dt 
)

Hmmm how should I implement this??? Maybe change the L1Track to use stubs not LCTs?

Definition at line 60 of file CSCTFTrackBuilder.cc.

References CSCDetId, end, CSCTFSectorProcessor::filteredStubs(), CSCTriggerContainer< T >::get(), m_maxBX, m_minBX, CSCDetId::maxEndcapId(), CSCTriggerNumbering::maxTriggerSectorId(), CSCDetId::minEndcapId(), CSCTriggerNumbering::minTriggerSectorId(), my_SPs, CSCTriggerContainer< T >::push_back(), CSCTriggerContainer< T >::push_many(), DTTTrigCorrFirst::run, asciidump::s, edm::shift, CSCTFSectorProcessor::tracks(), and CSCTriggerNumbering::triggerSubSectorFromLabels().

Referenced by CSCTFTrackProducer::produce().

63 {
64  std::vector<csc::L1Track> trks;
66 
68 
69  for(Citer = lcts->begin(); Citer != lcts->end(); Citer++)
70  {
71  CSCCorrelatedLCTDigiCollection::const_iterator Diter = (*Citer).second.first;
72  CSCCorrelatedLCTDigiCollection::const_iterator Dend = (*Citer).second.second;
73 
74  for(; Diter != Dend; Diter++)
75  {
76  csctf::TrackStub theStub((*Diter),(*Citer).first);
77  stub_list.push_back(theStub);
78  }
79  }
80 
81  // Now we append the track stubs the the DT Sector Collector
82  // after processing from the DT Receiver.
83 
84 // CSCTriggerContainer<csctf::TrackStub> dtstubs = my_dtrc->process(dttrig);
85 // stub_list.push_many(dtstubs);
86  stub_list.push_many(*dtstubss);
87 
88  // run each sector processor in the TF
89  for(int e = CSCDetId::minEndcapId(); e <= CSCDetId::maxEndcapId(); ++e)
90  {
93  {
94  CSCTriggerContainer<csctf::TrackStub> current_e_s = stub_list.get(e, s);
95  if(my_SPs[e-1][s-1]->run(current_e_s))
96  {
97  std::vector<csc::L1Track> theTracks = my_SPs[e-1][s-1]->tracks().get();
98  trks.insert(trks.end(), theTracks.begin(), theTracks.end());
99  }
100  stubs_to_dt->push_many(my_SPs[e-1][s-1]->dtStubs()); // send stubs whether or not we find a track!!!
101  }
102  }
103 
104  // Now to combine tracks with their track stubs and send them off.
105  trkcoll->resize(trks.size());
106  std::vector<csc::L1Track>::const_iterator titr = trks.begin();
107  L1CSCTrackCollection::iterator tcitr = trkcoll->begin();
108 
109  for(; titr != trks.end(); titr++)
110  {
111  tcitr->first = (*titr);
112  std::vector<csctf::TrackStub> possible_stubs = my_SPs[titr->endcap()-1][titr->sector()-1]->filteredStubs();
113  std::vector<csctf::TrackStub>::const_iterator tkstbs = possible_stubs.begin();
114 
115  int me1ID = titr->me1ID();
116  int me2ID = titr->me2ID();
117  int me3ID = titr->me3ID();
118  int me4ID = titr->me4ID();
119  int mb1ID = titr->mb1ID();
120  int me1delay = titr->me1Tbin();
121  int me2delay = titr->me2Tbin();
122  int me3delay = titr->me3Tbin();
123  int me4delay = titr->me4Tbin();
124  int mb1delay = titr->mb1Tbin();
125  // BX analyzer: some stub could be delayed by BXA so that all the stubs will run through the core at the same BX;
126  // then there is a rule of "second earlies LCT": resulting track will be placed at BX of the "second earliest LCT";
127  // in the end there are two parameters in place: the delay by BXA w.r.t to the last LCT and track tbin assignment
128  std::map< int, std::list<int> > timeline;
129  if( me1ID ) timeline[me1delay].push_back(1);
130  if( me2ID ) timeline[me2delay].push_back(2);
131  if( me3ID ) timeline[me3delay].push_back(3);
132  if( me4ID ) timeline[me4delay].push_back(4);
133  int earliest_tbin=0, second_earliest_tbin=0;
134  for(int bx=7; bx>=0; bx--){
135  std::list<int>::const_iterator iter = timeline[bx].begin();
136  while( iter != timeline[bx].end() ){
137  if( earliest_tbin==0 ) earliest_tbin=bx;
138  else if( second_earliest_tbin==0 ) second_earliest_tbin=bx;
139  iter++;
140  }
141  }
142  // Core's input was loaded in a relative time window BX=[0-7)
143  // To relate it to time window of tracks (centred at BX=0) we introduce a shift:
144  int shift = (m_maxBX + m_minBX)/2 - m_minBX + m_minBX;
145  int me1Tbin = titr->bx() - me1delay + second_earliest_tbin + shift;
146  int me2Tbin = titr->bx() - me2delay + second_earliest_tbin + shift;
147  int me3Tbin = titr->bx() - me3delay + second_earliest_tbin + shift;
148  int me4Tbin = titr->bx() - me4delay + second_earliest_tbin + shift;
149  int mb1Tbin = titr->bx() - mb1delay + second_earliest_tbin + shift;
150 
151  for(; tkstbs != possible_stubs.end(); tkstbs++)
152  {
153  switch(tkstbs->station())
154  {
155  case 1:
156  if((tkstbs->getMPCLink()
157  +(3*(CSCTriggerNumbering::triggerSubSectorFromLabels(CSCDetId(tkstbs->getDetId().rawId())) - 1))) == me1ID && me1ID != 0 && me1Tbin == tkstbs->BX() )
158  {
159  tcitr->second.insertDigi(CSCDetId(tkstbs->getDetId().rawId()), *(tkstbs->getDigi()));
160  }
161  break;
162  case 2:
163  if(tkstbs->getMPCLink() == me2ID && me2ID != 0 && me2Tbin == tkstbs->BX() )
164  {
165  tcitr->second.insertDigi(CSCDetId(tkstbs->getDetId().rawId()), *(tkstbs->getDigi()));
166  }
167  break;
168  case 3:
169  if(tkstbs->getMPCLink() == me3ID && me3ID != 0 && me3Tbin == tkstbs->BX() )
170  {
171  tcitr->second.insertDigi(CSCDetId(tkstbs->getDetId().rawId()), *(tkstbs->getDigi()));
172  }
173  break;
174  case 4:
175  if(tkstbs->getMPCLink() == me4ID && me4ID != 0 && me4Tbin == tkstbs->BX() )
176  {
177  tcitr->second.insertDigi(CSCDetId(tkstbs->getDetId().rawId()), *(tkstbs->getDigi()));
178  }
179  break;
180  case 5:
181  if(tkstbs->getMPCLink() == mb1ID && mb1ID != 0 && mb1Tbin == tkstbs->BX() )
182  {
184  }
185  break;
186  default:
187  edm::LogWarning("CSCTFTrackBuilder::buildTracks()") << "SERIOUS ERROR: STATION " << tkstbs->station() << " NOT IN RANGE [1,5]\n";
188  };
189  }
190  tcitr++; // increment to next track in the collection
191  }
192 
193 }
std::vector< csctf::TrackStub > filteredStubs() const
std::vector< T > get() const
static int minEndcapId()
Definition: CSCDetId.h:238
static int maxTriggerSectorId()
void push_back(const T data)
void push_many(const std::vector< T > data)
static int minTriggerSectorId()
static int maxEndcapId()
Definition: CSCDetId.h:239
#define end
Definition: vmac.h:38
std::vector< DigiType >::const_iterator const_iterator
static int triggerSubSectorFromLabels(int station, int chamber)
CSCTriggerContainer< csc::L1Track > tracks() const
static unsigned int const shift
string s
Definition: asciidump.py:422
CSCTFSectorProcessor * my_SPs[nEndcaps][nSectors]
void CSCTFTrackBuilder::initialize ( const edm::EventSetup c)

Definition at line 36 of file CSCTFTrackBuilder.cc.

References CSCTFSectorProcessor::initialize(), CSCDetId::maxEndcapId(), CSCTriggerNumbering::maxTriggerSectorId(), CSCDetId::minEndcapId(), CSCTriggerNumbering::minTriggerSectorId(), my_SPs, and asciidump::s.

Referenced by CSCTFTrackProducer::produce().

36  {
37  //my_dtrc->initialize(c);
38  for(int e = CSCDetId::minEndcapId(); e <= CSCDetId::maxEndcapId(); ++e)
39  {
41  {
42  my_SPs[e-1][s-1]->initialize(c);
43  }
44  }
45 }
static int minEndcapId()
Definition: CSCDetId.h:238
static int maxTriggerSectorId()
void initialize(const edm::EventSetup &c)
KK.
static int minTriggerSectorId()
static int maxEndcapId()
Definition: CSCDetId.h:239
string s
Definition: asciidump.py:422
CSCTFSectorProcessor * my_SPs[nEndcaps][nSectors]

Member Data Documentation

int CSCTFTrackBuilder::m_maxBX
private

Definition at line 39 of file CSCTFTrackBuilder.h.

Referenced by buildTracks(), and CSCTFTrackBuilder().

int CSCTFTrackBuilder::m_minBX
private

Definition at line 39 of file CSCTFTrackBuilder.h.

Referenced by buildTracks(), and CSCTFTrackBuilder().

CSCTFSectorProcessor* CSCTFTrackBuilder::my_SPs[nEndcaps][nSectors]
private

Definition at line 38 of file CSCTFTrackBuilder.h.

Referenced by buildTracks(), CSCTFTrackBuilder(), initialize(), and ~CSCTFTrackBuilder().