CMS 3D CMS Logo

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

#include <RBCProcessRPCDigis.h>

Inheritance diagram for RBCProcessRPCDigis:
ProcessInputSignal

Public Member Functions

void builddata ()
 
void configure ()
 
void initialize (std::vector< RPCData * > &)
 
int next ()
 
void print_output ()
 
 RBCProcessRPCDigis ()
 Standard constructor. More...
 
 RBCProcessRPCDigis (const edm::ESHandle< RPCGeometry > &, const edm::Handle< RPCDigiCollection > &)
 
void reset ()
 
RPCInputSignalretrievedata ()
 
void rewind ()
 
void showfirst ()
 
virtual ~RBCProcessRPCDigis ()
 Destructor. More...
 
- Public Member Functions inherited from ProcessInputSignal
virtual ~ProcessInputSignal ()
 

Private Member Functions

int getBarrelLayer (const int &, const int &)
 
void setDigiAt (int, int)
 
void setInputBit (std::bitset< 15 > &, int)
 

Private Attributes

RPCDatam_block
 
std::map< int, RBCInput * > m_data
 
bool m_debug
 
RPCDigiCollection::DigiRangeIterator m_detUnitItr
 
std::map< int, Counters * > m_digiCounters
 
RPCDigiCollection::const_iterator m_digiItr
 
std::map< int, int > m_layermap
 
RPCInputSignalm_lbin
 
int m_maxBxWindow
 
const edm::Handle
< RPCDigiCollection > * 
m_ptr_digiColl
 
const edm::ESHandle
< RPCGeometry > * 
m_ptr_rpcGeom
 
std::vector< int > m_sec1id
 
std::vector< int > m_sec2id
 
std::map< int, std::vector
< RPCData * > > 
m_vecDataperBx
 
std::vector< int > m_wheelid
 

Detailed Description

Author
Andres Felipe Osorio Oliveros
Date
2009-04-15

Definition at line 39 of file RBCProcessRPCDigis.h.

Constructor & Destructor Documentation

RBCProcessRPCDigis::RBCProcessRPCDigis ( )
inline

Standard constructor.

Definition at line 42 of file RBCProcessRPCDigis.h.

42 {};
RBCProcessRPCDigis::RBCProcessRPCDigis ( const edm::ESHandle< RPCGeometry > &  rpcGeom,
const edm::Handle< RPCDigiCollection > &  digiColl 
)

Definition at line 20 of file RBCProcessRPCDigis.cc.

References configure(), m_data, m_debug, m_lbin, m_ptr_digiColl, and m_ptr_rpcGeom.

22 {
23 
24  m_ptr_rpcGeom = & rpcGeom;
25  m_ptr_digiColl = & digiColl;
26 
27  m_lbin = dynamic_cast<RPCInputSignal*>( new RBCLinkBoardGLSignal( &m_data ) );
28 
29  m_debug = false;
30 
31  configure();
32 
33 }
const edm::ESHandle< RPCGeometry > * m_ptr_rpcGeom
const edm::Handle< RPCDigiCollection > * m_ptr_digiColl
RPCInputSignal * m_lbin
std::map< int, RBCInput * > m_data
RBCProcessRPCDigis::~RBCProcessRPCDigis ( )
virtual

Destructor.

Definition at line 92 of file RBCProcessRPCDigis.cc.

References m_digiCounters, m_layermap, m_lbin, m_sec1id, m_sec2id, m_wheelid, and reset().

92  {
93 
94  if ( m_lbin ) delete m_lbin;
95 
96  std::vector<int>::iterator wheel;
97 
98  for( wheel = m_wheelid.begin(); wheel != m_wheelid.end(); ++wheel)
99  delete m_digiCounters[(*wheel)];
100 
101  m_sec1id.clear();
102  m_sec2id.clear();
103  m_wheelid.clear();
104  m_layermap.clear();
105 
106  reset();
107 
108 }
std::map< int, Counters * > m_digiCounters
std::vector< int > m_wheelid
std::vector< int > m_sec1id
std::vector< int > m_sec2id
RPCInputSignal * m_lbin
std::map< int, int > m_layermap

Member Function Documentation

void RBCProcessRPCDigis::builddata ( )

Definition at line 268 of file RBCProcessRPCDigis.cc.

References abs, gather_cfg::cout, gen::k, m_data, m_debug, m_vecDataperBx, and RBCInput::needmapping.

Referenced by next().

269 {
270 
271  int bx(0);
272  int code(0);
273  int bxsign(1);
274  std::vector<RPCData*>::iterator itr;
275  std::map<int, std::vector<RPCData*> >::iterator itr2;
276 
277  itr2 = m_vecDataperBx.begin();
278  if( itr2 == ( m_vecDataperBx.end() ) ) return;
279 
280  while ( itr2 != m_vecDataperBx.end() ) {
281 
282  bx = (*itr2).first;
283 
284  if ( bx != 0 ) bxsign = ( bx / abs(bx) );
285  else bxsign = 1;
286 
287  for(itr = (*itr2).second.begin(); itr != (*itr2).second.end(); ++itr) {
288 
289  for(int k=0; k < 6; ++k) {
290 
291  code = bxsign * ( 1000000*abs(bx)
292  + 10000*(*itr)->wheelIdx()
293  + 100 *(*itr)->m_sec1[k]
294  + 1 *(*itr)->m_sec2[k] );
295 
296  RBCInput * signal = & (*itr)->m_orsignals[k];
297  signal->needmapping = false;
298 
299  if ( signal->hasData )
300  m_data.insert( std::make_pair( code , signal) );
301 
302  }
303  }
304 
305  ++itr2;
306 
307  }
308 
309  if ( m_debug ) std::cout << "builddata: completed. size of data: " << m_data.size() << std::endl;
310 
311 }
#define abs(x)
Definition: mlp_lapack.h:159
int k[5][pyjets_maxn]
bool needmapping
Definition: RBCInput.h:76
tuple cout
Definition: gather_cfg.py:121
std::map< int, std::vector< RPCData * > > m_vecDataperBx
std::map< int, RBCInput * > m_data
void RBCProcessRPCDigis::configure ( )

Definition at line 35 of file RBCProcessRPCDigis.cc.

References m_digiCounters, m_layermap, m_maxBxWindow, m_sec1id, m_sec2id, and m_wheelid.

Referenced by RBCProcessRPCDigis().

36 {
37 
38  m_wheelid.push_back(-2); //-2
39  m_wheelid.push_back(-1); //-1
40  m_wheelid.push_back(0); // 0
41  m_wheelid.push_back( 1); //+1
42  m_wheelid.push_back( 2); //+2
43 
44  m_sec1id.push_back(12);
45  m_sec2id.push_back(1);
46  m_sec1id.push_back(2);
47  m_sec2id.push_back(3);
48  m_sec1id.push_back(4);
49  m_sec2id.push_back(5);
50  m_sec1id.push_back(6);
51  m_sec2id.push_back(7);
52  m_sec1id.push_back(8);
53  m_sec2id.push_back(9);
54  m_sec1id.push_back(10);
55  m_sec2id.push_back(11);
56 
57  m_layermap[113] = 0; //RB1InFw
58  m_layermap[123] = 1; //RB1OutFw
59 
60  m_layermap[20213] = 2; //RB22Fw
61  m_layermap[20223] = 2; //RB22Fw
62  m_layermap[30223] = 3; //RB23Fw
63  m_layermap[30213] = 3; //RB23Fw
64  m_layermap[30212] = 4; //RB23M
65  m_layermap[30222] = 4; //RB23M
66 
67  m_layermap[313] = 5; //RB3Fw
68  m_layermap[413] = 6; //RB4Fw
69  m_layermap[111] = 7; //RB1InBk
70  m_layermap[121] = 8; //RB1OutBk
71 
72  m_layermap[20211] = 9; //RB22Bw
73  m_layermap[20221] = 9; //RB22Bw
74  m_layermap[30211] = 10; //RB23Bw
75  m_layermap[30221] = 10; //RB23Bw
76 
77  m_layermap[311] = 11; //RB3Bk
78  m_layermap[411] = 12; //RB4Bk
79 
80  m_maxBxWindow = 3;
81 
82  std::vector<int>::iterator wheel;
83 
84  for( wheel = m_wheelid.begin(); wheel != m_wheelid.end(); ++wheel)
85  m_digiCounters[(*wheel)] = new Counters( (*wheel) );
86 
87 }
std::map< int, Counters * > m_digiCounters
std::vector< int > m_wheelid
std::vector< int > m_sec1id
std::vector< int > m_sec2id
std::map< int, int > m_layermap
int RBCProcessRPCDigis::getBarrelLayer ( const int &  _layer,
const int &  _station 
)
private

Definition at line 313 of file RBCProcessRPCDigis.cc.

Referenced by next().

314 {
315 
316  //... Calculates the generic Barrel Layer (1 to 6)
317  int blayer(0);
318 
319  if ( _station < 3 ) {
320  blayer = ( (_station-1) * 2 ) + _layer;
321  }
322  else {
323  blayer = _station + 2;
324  }
325 
326  return blayer;
327 
328 }
void RBCProcessRPCDigis::initialize ( std::vector< RPCData * > &  dataVec)

Definition at line 237 of file RBCProcessRPCDigis.cc.

References gather_cfg::cout, RBCInput::hasData, i, RBCInput::input_sec, j, m_block, m_debug, RPCData::m_orsignals, RPCData::m_sec1, m_sec1id, RPCData::m_sec2, m_sec2id, RPCData::m_wheel, m_wheelid, and RBCInput::needmapping.

Referenced by next().

238 {
239 
240  if ( m_debug ) std::cout << "initialize" << std::endl;
241 
242  int maxWheels = 5;
243  int maxRbcBrds = 6;
244 
245  for(int i=0; i < maxWheels; ++i) {
246 
247  m_block = new RPCData();
248 
250 
251  for(int j=0; j < maxRbcBrds; ++j) {
252  m_block->m_sec1[j] = m_sec1id[j];
253  m_block->m_sec2[j] = m_sec2id[j];
254  m_block->m_orsignals[j].input_sec[0].reset();
255  m_block->m_orsignals[j].input_sec[1].reset();
256  m_block->m_orsignals[j].needmapping = false;
257  m_block->m_orsignals[j].hasData = false;
258  }
259 
260  dataVec.push_back( m_block );
261 
262  }
263 
264  if ( m_debug ) std::cout << "initialize: completed" << std::endl;
265 
266 }
int i
Definition: DBlmapReader.cc:9
RBCInput * m_orsignals
Definition: RPCData.h:55
int * m_sec2
Definition: RPCData.h:54
std::vector< int > m_wheelid
std::vector< int > m_sec1id
int j
Definition: DBlmapReader.cc:9
std::bitset< 15 > * input_sec
Definition: RBCInput.h:65
std::vector< int > m_sec2id
int m_wheel
Definition: RPCData.h:52
bool needmapping
Definition: RBCInput.h:76
bool hasData
Definition: RBCInput.h:75
tuple cout
Definition: gather_cfg.py:121
int * m_sec1
Definition: RPCData.h:53
int RBCProcessRPCDigis::next ( void  )
virtual

Implements ProcessInputSignal.

Definition at line 111 of file RBCProcessRPCDigis.cc.

References abs, builddata(), gather_cfg::cout, getBarrelLayer(), RPCRoll::id(), initialize(), RPCRoll::isForward(), RPCDetId::layer(), m_block, m_data, m_debug, m_detUnitItr, m_digiCounters, m_digiItr, m_maxBxWindow, m_vecDataperBx, print_output(), reset(), RPCDetId::ring(), RPCDetId::sector(), setDigiAt(), relativeConstraints::station, and RPCDetId::station().

Referenced by BeautifulSoup.PageElement::_invert().

111  {
112 
113  //...clean up previous data contents
114 
115  reset();
116 
117  int ndigis(0);
118 
119  for (m_detUnitItr = (*m_ptr_digiColl)->begin();
120  m_detUnitItr != (*m_ptr_digiColl)->end(); ++m_detUnitItr ) {
121 
122  if ( m_debug ) std::cout << "looping over digis 1 ..." << std::endl;
123 
124  m_digiItr = (*m_detUnitItr ).second.first;
125  int bx = (*m_digiItr).bx();
126 
127  if ( abs(bx) >= m_maxBxWindow ) {
128  if ( m_debug ) std::cout << "RBCProcessRPCDigis> found a bx bigger than max allowed: "
129  << bx << std::endl;
130  continue;
131  }
132 
133  const RPCDetId & id = (*m_detUnitItr).first;
134  const RPCRoll * roll = dynamic_cast<const RPCRoll* >( (*m_ptr_rpcGeom)->roll(id));
135 
136  if((roll->isForward())) {
137  if( m_debug ) std::cout << "RBCProcessRPCDigis: roll is forward" << std::endl;
138  continue;
139  }
140 
141  int wheel = roll->id().ring(); // -2,-1,0,+1,+2
142  int sector = roll->id().sector(); // 1 to 12
143  int layer = roll->id().layer(); // 1,2
144  int station = roll->id().station(); // 1-4
145  int blayer = getBarrelLayer( layer, station ); // 1 to 6
146  int rollid = id.roll();
147 
148  int digipos = (station * 100) + (layer * 10) + rollid;
149 
150  if ( (wheel == -1 || wheel == 0 || wheel == 1) && station == 2 && layer == 1 )
151  digipos = 30000 + digipos;
152  if ( (wheel == -2 || wheel == 2) && station == 2 && layer == 2 )
153  digipos = 30000 + digipos;
154 
155  if ( (wheel == -1 || wheel == 0 || wheel == 1) && station == 2 && layer == 2 )
156  digipos = 20000 + digipos;
157  if ( (wheel == -2 || wheel == 2) && station == 2 && layer == 1 )
158  digipos = 20000 + digipos;
159 
160  if ( m_debug ) std::cout << "Bx: " << bx << '\t'
161  << "Wheel: " << wheel << '\t'
162  << "Sector: " << sector << '\t'
163  << "Station: " << station << '\t'
164  << "Layer: " << layer << '\t'
165  << "B-Layer: " << blayer << '\t'
166  << "Roll id: " << rollid << '\t'
167  << "Digi at: " << digipos << '\n';
168 
169  //... Construct the RBCinput objects
170  std::map<int,std::vector<RPCData*> >::iterator itr;
171  itr = m_vecDataperBx.find( bx );
172 
173  if ( itr == m_vecDataperBx.end() ) {
174  if ( m_debug ) std::cout << "Found a new Bx: " << bx << std::endl;
175  std::vector<RPCData*> wheelData;
176  initialize(wheelData);
177  m_vecDataperBx[bx] = wheelData;
178  this->m_block = wheelData[ (wheel + 2) ];
179  setDigiAt( sector, digipos );
180  }
181  else{
182  this->m_block = (*itr).second[ (wheel + 2) ];
183  setDigiAt( sector, digipos );
184  }
185 
186  std::map<int, Counters* >::iterator wheelCounter;
187  wheelCounter = m_digiCounters.find( wheel );
188 
189  if ( wheelCounter != m_digiCounters.end() )
190  (*wheelCounter).second->incrementSector( sector );
191 
192  if ( m_debug ) std::cout << "looping over digis 2 ..." << std::endl;
193 
194  ++ndigis;
195 
196  }
197 
198  if ( m_debug ) std::cout << "size of data vectors: " << m_vecDataperBx.size() << std::endl;
199 
200  builddata();
201 
202  if ( m_debug ) {
203  std::cout << "after reset" << std::endl;
204  print_output();
205  }
206 
207  if ( m_debug ) std::cout << "RBCProcessRPCDigis: DataSize: " << m_data.size()
208  << " ndigis " << ndigis << std::endl;
209 
210  std::map<int, Counters* >::iterator wheelCounter;
211  for( wheelCounter = m_digiCounters.begin(); wheelCounter != m_digiCounters.end(); ++wheelCounter) {
212  (*wheelCounter).second->evalCounters();
213  if ( m_debug ) (*wheelCounter).second->printSummary();
214  }
215 
216  if ( m_data.size() <= 0 ) return 0;
217 
218  return 1;
219 
220 }
RPCDigiCollection::DigiRangeIterator m_detUnitItr
void initialize(std::vector< RPCData * > &)
std::map< int, Counters * > m_digiCounters
#define abs(x)
Definition: mlp_lapack.h:159
RPCDetId id() const
Definition: RPCRoll.cc:24
int ring() const
Definition: RPCDetId.h:76
void setDigiAt(int, int)
int layer() const
Definition: RPCDetId.h:112
RPCDigiCollection::const_iterator m_digiItr
int sector() const
Sector id: the group of chambers at same phi (and increasing r)
Definition: RPCDetId.h:106
int getBarrelLayer(const int &, const int &)
tuple cout
Definition: gather_cfg.py:121
std::map< int, std::vector< RPCData * > > m_vecDataperBx
bool isForward() const
Definition: RPCRoll.cc:98
std::map< int, RBCInput * > m_data
int station() const
Definition: RPCDetId.h:100
void RBCProcessRPCDigis::print_output ( )

Definition at line 373 of file RBCProcessRPCDigis.cc.

References gather_cfg::cout, and m_data.

Referenced by next().

374 {
375 
376  std::cout << "RBCProcessRPCDigis> Output starts" << std::endl;
377 
378  std::map<int,RBCInput*>::const_iterator itr;
379  for( itr = m_data.begin(); itr != m_data.end(); ++itr) {
380  std::cout << (*itr).first << '\t' << (* (*itr).second ) << '\n';
381  }
382 
383  std::cout << "RBCProcessRPCDigis> Output ends" << std::endl;
384 
385 }
tuple cout
Definition: gather_cfg.py:121
std::map< int, RBCInput * > m_data
void RBCProcessRPCDigis::reset ( void  )

Definition at line 222 of file RBCProcessRPCDigis.cc.

References m_vecDataperBx.

Referenced by next(), and ~RBCProcessRPCDigis().

223 {
224 
225  std::map<int,std::vector<RPCData*> >::iterator itr1;
226  for( itr1 = m_vecDataperBx.begin(); itr1 != m_vecDataperBx.end(); ++itr1) {
227  std::vector<RPCData*>::iterator itr2;
228  for(itr2 = (*itr1).second.begin(); itr2 != (*itr1).second.end();++itr2 )
229  if ( (*itr2) ) delete *itr2;
230  (*itr1).second.clear();
231  }
232  m_vecDataperBx.clear();
233 
234 }
std::map< int, std::vector< RPCData * > > m_vecDataperBx
RPCInputSignal* RBCProcessRPCDigis::retrievedata ( )
inlinevirtual

Implements ProcessInputSignal.

Definition at line 61 of file RBCProcessRPCDigis.h.

References m_lbin.

61  {
62  return m_lbin;
63  };
RPCInputSignal * m_lbin
void RBCProcessRPCDigis::rewind ( )
inline

Definition at line 65 of file RBCProcessRPCDigis.h.

65 {};
void RBCProcessRPCDigis::setDigiAt ( int  sector,
int  digipos 
)
private

Definition at line 331 of file RBCProcessRPCDigis.cc.

References gather_cfg::cout, spr::find(), RBCInput::hasData, RBCInput::input_sec, m_block, m_debug, RPCData::m_orsignals, RPCData::m_sec1, m_sec1id, m_sec2id, pos, and setInputBit().

Referenced by next().

332 {
333 
334  int pos = 0;
335  int isAoB = 0;
336 
337  if ( m_debug ) std::cout << "setDigiAt" << std::endl;
338 
339  std::vector<int>::const_iterator itr;
340  itr = std::find( m_sec1id.begin(), m_sec1id.end(), sector );
341 
342  if ( itr == m_sec1id.end()) {
343  itr = std::find( m_sec2id.begin(), m_sec2id.end(), sector );
344  isAoB = 1;
345  }
346 
347  for ( pos = 0; pos < 6; ++pos ) {
348  if (this->m_block->m_sec1[pos] == sector || this->m_block->m_sec2[pos] == sector )
349  break;
350  }
351 
352  if ( m_debug ) std::cout << this->m_block->m_orsignals[pos];
353 
354  setInputBit( this->m_block->m_orsignals[pos].input_sec[ isAoB ] , digipos );
355 
356  this->m_block->m_orsignals[pos].hasData = true;
357 
358  if ( m_debug ) std::cout << this->m_block->m_orsignals[pos];
359 
360  if ( m_debug ) std::cout << "setDigiAt completed" << std::endl;
361 
362 }
RBCInput * m_orsignals
Definition: RPCData.h:55
void setInputBit(std::bitset< 15 > &, int)
void find(edm::Handle< EcalRecHitCollection > &hits, DetId thisDet, std::vector< EcalRecHitCollection::const_iterator > &hit, bool debug=false)
Definition: FindCaloHit.cc:7
std::vector< int > m_sec1id
std::bitset< 15 > * input_sec
Definition: RBCInput.h:65
std::vector< int > m_sec2id
bool hasData
Definition: RBCInput.h:75
tuple cout
Definition: gather_cfg.py:121
int * m_sec1
Definition: RPCData.h:53
void RBCProcessRPCDigis::setInputBit ( std::bitset< 15 > &  signals,
int  digipos 
)
private

Definition at line 364 of file RBCProcessRPCDigis.cc.

References gather_cfg::cout, m_debug, and m_layermap.

Referenced by setDigiAt().

365 {
366 
367  int bitpos = m_layermap[digipos];
368  if( m_debug ) std::cout << "Bitpos: " << bitpos << std::endl;
369  signals.set( bitpos , 1 );
370 
371 }
list signals
Definition: scaleCards.py:71
std::map< int, int > m_layermap
tuple cout
Definition: gather_cfg.py:121
void RBCProcessRPCDigis::showfirst ( )
inline

Definition at line 66 of file RBCProcessRPCDigis.h.

66 {};

Member Data Documentation

RPCData* RBCProcessRPCDigis::m_block
private

Definition at line 84 of file RBCProcessRPCDigis.h.

Referenced by initialize(), next(), and setDigiAt().

std::map<int, RBCInput*> RBCProcessRPCDigis::m_data
private

Definition at line 90 of file RBCProcessRPCDigis.h.

Referenced by builddata(), next(), print_output(), and RBCProcessRPCDigis().

bool RBCProcessRPCDigis::m_debug
private
RPCDigiCollection::DigiRangeIterator RBCProcessRPCDigis::m_detUnitItr
private

Definition at line 82 of file RBCProcessRPCDigis.h.

Referenced by next().

std::map<int, Counters*> RBCProcessRPCDigis::m_digiCounters
private

Definition at line 101 of file RBCProcessRPCDigis.h.

Referenced by configure(), next(), and ~RBCProcessRPCDigis().

RPCDigiCollection::const_iterator RBCProcessRPCDigis::m_digiItr
private

Definition at line 81 of file RBCProcessRPCDigis.h.

Referenced by next().

std::map<int, int> RBCProcessRPCDigis::m_layermap
private

Definition at line 88 of file RBCProcessRPCDigis.h.

Referenced by configure(), setInputBit(), and ~RBCProcessRPCDigis().

RPCInputSignal* RBCProcessRPCDigis::m_lbin
private

Definition at line 86 of file RBCProcessRPCDigis.h.

Referenced by RBCProcessRPCDigis(), retrievedata(), and ~RBCProcessRPCDigis().

int RBCProcessRPCDigis::m_maxBxWindow
private

Definition at line 95 of file RBCProcessRPCDigis.h.

Referenced by configure(), and next().

const edm::Handle<RPCDigiCollection>* RBCProcessRPCDigis::m_ptr_digiColl
private

Definition at line 79 of file RBCProcessRPCDigis.h.

Referenced by RBCProcessRPCDigis().

const edm::ESHandle<RPCGeometry>* RBCProcessRPCDigis::m_ptr_rpcGeom
private

Definition at line 78 of file RBCProcessRPCDigis.h.

Referenced by RBCProcessRPCDigis().

std::vector<int> RBCProcessRPCDigis::m_sec1id
private

Definition at line 98 of file RBCProcessRPCDigis.h.

Referenced by configure(), initialize(), setDigiAt(), and ~RBCProcessRPCDigis().

std::vector<int> RBCProcessRPCDigis::m_sec2id
private

Definition at line 99 of file RBCProcessRPCDigis.h.

Referenced by configure(), initialize(), setDigiAt(), and ~RBCProcessRPCDigis().

std::map<int, std::vector<RPCData*> > RBCProcessRPCDigis::m_vecDataperBx
private

Definition at line 92 of file RBCProcessRPCDigis.h.

Referenced by builddata(), next(), and reset().

std::vector<int> RBCProcessRPCDigis::m_wheelid
private

Definition at line 97 of file RBCProcessRPCDigis.h.

Referenced by configure(), initialize(), and ~RBCProcessRPCDigis().