test
CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
RPCGeomServ.cc
Go to the documentation of this file.
1 #include "Geometry/RPCGeometry/interface/RPCGeomServ.h"
2 #include "DataFormats/MuonDetId/interface/RPCDetId.h"
3 #include <sstream>
4 #include <iomanip>
5 
6 RPCGeomServ::RPCGeomServ::RPCGeomServ( const RPCDetId& id )
7  : _id( &id ),
8  _n( "" ),
9  _sn( "" ),
10  _cn( "" ),
11  _t( -99 ),
12  _z( true ),
13  _a( true )
14 {}
15 
16 RPCGeomServ::~RPCGeomServ( void )
17 {}
18 
19 std::string
20 RPCGeomServ::name( void )
21 {
22  if( _n.size() < 1 )
23  {
24  int station = _id->station();
25  int region = _id->region();
26  int roll = _id->roll();
27  int ring = _id->ring();
28  int layer = _id->layer();
29  int sector = _id->sector();
30  int subsector = _id->subsector();
31 
32  std::stringstream os;
33 
34  if( region == 0 )
35  {
36  os << "W";
37  os << std::setw(2) << std::setfill('+') << ring
38  << std::setfill(' ') << "_";
39 
40  os << "RB" << station;
41  if ( station <= 2){
42  (layer == 1 ) ? os << "in" : os << "out";
43 
44  }else if( station > 2 )
45  {
46  if( sector == 4 && station == 4 )
47  {
48  if( subsector == 1 )
49  {
50  os << "--";
51  }
52  else if( subsector == 2 )
53  {
54  os << "-";
55  }
56  else if( subsector == 3 )
57  {
58  os << "+";
59  }
60  else if( subsector == 4 )
61  {
62  os << "++";
63  }
64  }
65 
66  if( station == 3 )
67  {
68  if( subsector == 1 )
69  os << "-";
70  else
71  os << "+";
72  }
73  else if( station == 4
74  && sector != 9
75  && sector !=11
76  && sector != 4 )
77  {
78  if( subsector == 1 )
79  os << "-";
80  else
81  os << "+";
82  }
83  }
84 
85  os << "_";
86  os << "S" << std::setw(2) << std::setfill('0')
87  << sector << std::setfill(' ');
88  if( roll == 1 )
89  os << "_Backward";
90  else if( roll == 3 )
91  os << "_Forward";
92  else if( roll == 2 )
93  os << "_Middle";
94  }
95  else
96  {
97  os << "RE";
98 
99  os << std::setw(2) << std::setfill('+') << station * region
100  << std::setfill(' ') << "_";
101 
102  os << "R" << ring;
103  os << "_CH" << std::setw(2) << std::setfill('0') << this->segment();
104 
105  if( roll == 1 )
106  os << "_A";
107  else if( roll == 2 )
108  os << "_B";
109  else if( roll == 3 )
110  os << "_C";
111  else if( roll == 4 )
112  os << "_D";
113  else if( roll == 5 )
114  os << "_E";
115  }
116  _n = os.str();
117  }
118  return _n;
119 }
120 
121 std::string
122 RPCGeomServ::chambername()
123 {
124  if( _cn.size() < 1 )
125  {
126  int station = _id->station();
127  int region = _id->region();
128  int ring = _id->ring();
129  int layer = _id->layer();
130  int sector = _id->sector();
131  int subsector = _id->subsector();
132 
133  std::stringstream os;
134 
135  if( region == 0 )
136  {
137  os << "W";
138 
139  os << std::setw(2) << std::setfill('+') << ring
140  << std::setfill(' ') << "_";
141 
142  os << "RB" << station;
143  if ( station <= 2 ) {
144 
145  (layer == 1 ) ? os << "in" : os << "out";
146 
147  }else if( station > 2 )
148  {
149  if( sector == 4 && station == 4 )
150  {
151  if( subsector == 1 )
152  {
153  os << "--";
154  }
155  else if( subsector == 2 )
156  {
157  os << "-";
158  }
159  else if( subsector == 3 )
160  {
161  os << "+";
162  }
163  else if( subsector == 4 )
164  {
165  os <<"++";
166  }
167  }
168 
169  if( station == 3 )
170  {
171  if( subsector == 1 )
172  os << "-";
173  else
174  os << "+";
175  }
176  else if( station == 4
177  && sector != 9
178  && sector != 11
179  && sector != 4 )
180  {
181  if( subsector == 1 )
182  os << "-";
183  else
184  os << "+";
185  }
186  }
187 
188  os << "_";
189  os << "S" << std::setw(2) << std::setfill('0')
190  << sector << std::setfill(' ');
191  }
192  else
193  {
194  os << "RE";
195 
196  os << std::setw(2) << std::setfill('+') << station * region
197  << std::setfill(' ') << "_";
198 
199  os << "R" << ring;
200  os << "_CH" << std::setw(2) << std::setfill('0') << this->segment();
201  }
202  _cn = os.str();
203  }
204  return _cn;
205 }
206 
207 std::string
208 RPCGeomServ::shortname( void )
209 {
210  if( _sn.size() < 1 )
211  {
212  int station = _id->station();
213  int region = _id->region();
214  int roll = _id->roll();
215  int ring = _id->ring();
216  int layer = _id->layer();
217  int sector = _id->sector();
218  int subsector = _id->subsector();
219 
220  std::stringstream os;
221 
222  if( region == 0 )
223  {
224  os << "RB" << station;
225  if( station <= 2 ){
226 
227  (layer == 1 ) ? os << "in" : os << "out";
228 
229  }
230  else
231  {
232  if( sector == 4 && station == 4 )
233  {
234  if( subsector == 1 )
235  {
236  os << "--";
237  }
238  else if( subsector == 2 )
239  {
240  os << ",-";
241  }
242  else if( subsector == 3 )
243  {
244  os << "+";
245  }
246  else if( subsector == 4 )
247  {
248  os << "++";
249  }
250  }
251  else
252  {
253  if( subsector == 1 )
254  os << ",-";
255  else
256  os << "+";
257  }
258  }
259  if( roll == 1 )
260  os << " B";
261  else if( roll == 3 )
262  os << " F";
263  else if( roll == 2 )
264  os << " M";
265  }
266  else
267  {
268  os << "Ri" << ring << " Su" << subsector;
269  }
270  _sn = os.str();
271  }
272  return _sn;
273 }
274 
275 // returns a vector with number of channels for each chip in each FEB
276 std::vector<int>
277 RPCGeomServ::channelInChip( void )
278 {
279  std::vector<int> chipCh(4,8);//Endcap
280 
281  if(_id->region()==0){//Barrel
282  chipCh.clear();
283 
284  int station = _id->station();
285 
286  if (station<3 && _id->layer()==1){ // i.e. RB1in ||RB2in
287  chipCh.push_back(7);
288  chipCh.push_back(8);
289  }else if (station == 1 || station == 3){//i.e. RB1out || RB3
290  chipCh.push_back(7);
291  chipCh.push_back(7);
292  }else if (station == 2){// i.e. RB2out
293  chipCh.push_back(6);
294  chipCh.push_back(8);
295  }else if (_id->sector() == 4 || _id->sector()==10 ||(_id->sector() == 8 && _id->subsector()!=1) || (_id->sector() == 12 && _id->subsector()==1)){
296  chipCh.push_back(6);//i.e. Sector 4 & 10 RB4 and Sector 8 &12 RB4+
297  chipCh.push_back(6);
298  }else {
299  chipCh.push_back(8);
300  chipCh.push_back(8);
301  }
302  }
303 
304  return chipCh;
305 }
306 
307 
308 int
309 RPCGeomServ::eta_partition()
310 {
311  if (_t<-90){
312  if (_id->region() == 0 ){
313  if (this->inverted()) {
314  _t = 3*(_id->ring())+ (3-_id->roll())-1;
315  }else{
316  _t = 3*(_id->ring())+ _id->roll()-2;
317  }
318  }else{
319  _t = _id->region() * (3*(3-_id->ring()) + _id->roll() + 7);
320  }
321  }
322  return _t;
323 }
324 
325 int
326 RPCGeomServ::chambernr()
327 {
328 
329  // Station1
330  if( _id->station() ==1) {
331 
332  // in
333  if(_id->layer() ==1) {
334 
335  if(_id->roll()==1)
336  _cnr = 1;
337  else
338  _cnr = 2;
339  }
340  //out
341  else
342  {
343  if(_id->roll()==1)
344  _cnr = 3;
345  else
346  _cnr = 4;
347 
348  }
349  }
350 
351 
352  //Station 2
353  if (_id->station()==2) {
354 
355  //in
356  if(_id->layer()==1) {
357 
358  if(_id->roll()==1)//backward
359  _cnr = 5;
360  if(_id->roll()==3)//forward
361  _cnr=6;
362  if(_id->roll()==2)//middle
363  _cnr=7;
364  }
365  //out
366  else {
367 
368  if(_id->roll()==2)
369 
370  _cnr=7;
371 
372  if(_id->roll()==1)
373  _cnr=8;
374  if(_id->roll()==3)
375  _cnr=9;
376 
377  }
378  }
379 
380  //RB3- RB3+
381  if(_id->station()==3)
382  {
383  if(_id->subsector()==1) {
384 
385  if(_id->roll()==1)
386  _cnr=10;
387  else
388  _cnr=11;
389  }
390  else {
391 
392  if(_id->roll()==1)
393  _cnr=12;
394  else
395  _cnr=13;
396  }
397 
398  }
399 
400  //RB4
401  if(_id->station()==4) {
402 
403  if (_id->sector()== 4) {
404 
405  if ( _id->subsector()==2){//RB4-
406 
407  if(_id->roll()==1)
408  _cnr=14;
409  else
410  _cnr=15;
411 
412  }
413 
414  if ( _id->subsector()==3){//RB4+
415 
416  if(_id->roll()==1)
417  _cnr=16;
418  else
419  _cnr=17;
420 
421  }
422 
423  if ( _id->subsector()==1) {//RB4--
424 
425  if(_id->roll()==1)
426  _cnr=18;
427  else
428  _cnr=19;
429  }
430 
431  if ( _id->subsector()==4){//RB4++
432 
433  if(_id->roll()==1)
434  _cnr=20;
435  else
436  _cnr=21;
437 
438  }
439 
440  }
441 
442  else
443 
444  {
445  if(_id->subsector()==1) {
446 
447  if(_id->roll()==1)
448  _cnr=14;
449  else
450  _cnr=15;
451  }
452  else {
453 
454  if(_id->roll()==1)
455  _cnr=16;
456  else
457  _cnr=17;
458  }
459  }
460  }
461 
462 
463  // _cnr=10;
464  return _cnr;
465 
466 }
467 
468 int
469 RPCGeomServ::segment( void )
470 {
471  int nsub = 6;
472  int station = _id->station();
473  int ring = _id->ring();
474  ( ring == 1 && station > 1 ) ? nsub = 3 : nsub = 6;
475 
476  return( _id->subsector() + nsub * ( _id->sector() - 1 ));
477 }
478 
479 bool
480 RPCGeomServ::inverted()
481 {
482  // return !(this->zpositive() && this->aclockwise());
483  return !(this->zpositive());
484 }
485 
486 
487 bool
488 RPCGeomServ::zpositive()
489 {
490  if (_id->region()==0 && _t<-90 ){
491  if (_id->ring()<0){
492  _z=false;
493  }
494  if (_id->ring()==0){
495  if (_id->sector() == 1 || _id->sector() == 4 ||
496  _id->sector() == 5 || _id->sector() == 8 ||
497  _id->sector() == 9 || _id->sector() == 12){
498  _z=false;
499  }
500  }
501  }
502 
503  return _z;
504 }
505 
506 bool
507 RPCGeomServ::aclockwise()
508 {
509  if (_id->region()==0 && _t<-90 ){
510  if (_id->ring() > 0){
511  if (_id->layer()==2){
512  _a=false;
513  }
514  }else if(_id->ring() <0){
515  if (_id->layer()==1){
516  _a=false;
517  }
518  }else if(_id->ring() ==0) {
519  if ((_id->sector() == 1 || _id->sector() == 4 ||
520  _id->sector() == 5 || _id->sector() == 8 ||
521  _id->sector() == 9 || _id->sector() == 12) && _id->layer()==1)
522  _a=false;
523  else if ((_id->sector() == 2 || _id->sector() == 3 ||
524  _id->sector() == 6 || _id->sector() == 7 ||
525  _id->sector() == 10|| _id->sector() == 11) && _id->layer()==2)
526  _a=false;
527  }
528  }
529  return _a;
530 }
531 
532 RPCGeomServ::RPCGeomServ() : _id(0), _n(""), _sn(""), _cn(""), _t (-99), _z(false), _a(false)
533 {}
534 
535 
RPCGeomServ::chambernr
virtual int chambernr()
Definition: RPCGeomServ.cc:326
RPCGeomServ::chambername
virtual std::string chambername()
Definition: RPCGeomServ.cc:122
RPCGeomServ::_sn
std::string _sn
Definition: RPCGeomServ.h:29
RPCGeomServ::eta_partition
virtual int eta_partition()
Definition: RPCGeomServ.cc:309
RPCGeomServ::channelInChip
std::vector< int > channelInChip()
Definition: RPCGeomServ.cc:277
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
RPCGeomServ::_n
std::string _n
Definition: RPCGeomServ.h:28
RPCGeomServ::name
virtual std::string name()
Definition: RPCGeomServ.cc:20
RPCDetId::roll
int roll() const
Definition: RPCDetId.h:120
RPCDetId::ring
int ring() const
Definition: RPCDetId.h:72
RPCGeomServ::aclockwise
virtual bool aclockwise()
Definition: RPCGeomServ.cc:507
RPCGeomServ::~RPCGeomServ
virtual ~RPCGeomServ()
Definition: RPCGeomServ.cc:16
RPCGeomServ::shortname
virtual std::string shortname()
Definition: RPCGeomServ.cc:208
RPCDetId::layer
int layer() const
Definition: RPCDetId.h:108
RPCGeomServ::segment
virtual int segment()
Definition: RPCGeomServ.cc:469
HLT_25ns10e33_v2_cff.region
tuple region
Definition: HLT_25ns10e33_v2_cff.py:61599
RPCGeomServ::zpositive
virtual bool zpositive()
Definition: RPCGeomServ.cc:488
RPCDetId::sector
int sector() const
Sector id: the group of chambers at same phi (and increasing r)
Definition: RPCDetId.h:102
RPCGeomServ::_cn
std::string _cn
Definition: RPCGeomServ.h:30
RPCDetId::subsector
int subsector() const
SubSector id : some sectors are divided along the phi direction in subsectors (from 1 to 4 in Barrel...
Definition: RPCDetId.h:114
RPCGeomServ::inverted
virtual bool inverted()
Definition: RPCGeomServ.cc:480
edm::false
volatile std::atomic< bool > shutdown_flag false
Definition: UnixSignalHandlers.cc:22
RPCGeomServ::_id
const RPCDetId * _id
Definition: RPCGeomServ.h:27
RPCDetId::region
int region() const
Region id: 0 for Barrel, +/-1 For +/- Endcap.
Definition: RPCDetId.h:63
RPCDetId::station
int station() const
Definition: RPCDetId.h:96
Generated for CMSSW Reference Manual by doxygen 1.8.5