RPCGeomServ.cc

Go to the documentation of this file.

#include "Geometry/RPCGeometry/interface/RPCGeomServ.h"
#include "DataFormats/MuonDetId/interface/RPCDetId.h"
#include <sstream>
#include <iomanip>

RPCGeomServ::RPCGeomServ::RPCGeomServ( const RPCDetId& id )
  : _id( &id ),
    _n( "" ),
    _sn( "" ),
    _cn( "" ),
    _t( -99 ),
    _z( true ),
    _a( true )
{}

RPCGeomServ::~RPCGeomServ( void )
{}
 
std::string 
RPCGeomServ::name( void )
{
  if( _n.empty() )
  {
    int station = _id->station();
    int region = _id->region();
    int roll = _id->roll();
    int ring = _id->ring();
    int layer = _id->layer();
    int sector = _id->sector();
    int subsector = _id->subsector();

    std::stringstream os;
    
    if( region == 0 )
    {
      os << "W";
      os << std::setw(2) << std::setfill('+') << ring
     << std::setfill(' ') << "_";

      os << "RB" << station;
      if ( station <= 2){
    (layer == 1 ) ? os << "in" : os << "out";
    
      }else if( station > 2 )
      {
    if( sector == 4 && station == 4 )
    {
      if( subsector == 1 )
      {
        os << "--";
      }
      else if( subsector == 2 )
      {
        os << "-";
      }
      else if( subsector == 3 )
      {
        os << "+";
      }
      else if( subsector == 4 )
      {
        os << "++";
      }
    }
      
    if( station == 3 )
    {
      if( subsector == 1 )
        os << "-";
      else
        os << "+";
    }   
    else if( station == 4
         && sector != 9
         && sector !=11
         && sector != 4 )
    {
      if( subsector == 1 )
        os << "-";
      else
        os << "+";
    }
      }
    
      os << "_";
      os << "S" << std::setw(2) << std::setfill('0')
     << sector << std::setfill(' ');
      if( roll == 1 )
    os << "_Backward";
      else if( roll == 3 )
    os << "_Forward";
      else if( roll == 2 )
    os << "_Middle";
    }
    else
    {
      os << "RE";

      os << std::setw(2) << std::setfill('+') << station * region
     << std::setfill(' ') << "_";

      os << "R" << ring;
      os << "_CH" << std::setw(2) << std::setfill('0') << this->segment();

      if( roll == 1 )
    os << "_A";
      else if( roll == 2 )
    os << "_B";
      else if( roll == 3 )
    os << "_C";
      else if( roll == 4 )
    os << "_D";
      else if( roll == 5 )
    os << "_E";
    }
    _n = os.str();
  }
  return _n;
}

std::string 
RPCGeomServ::chambername()
{
  if( _cn.empty() )
  {
    int station = _id->station();
    int region = _id->region();
    int ring = _id->ring();
    int layer = _id->layer();
    int sector = _id->sector();
    int subsector = _id->subsector();

    std::stringstream os;
    
    if( region == 0 )
    {
      os << "W";

      os << std::setw(2) << std::setfill('+') << ring
     << std::setfill(' ') << "_";

      os << "RB" << station;
      if ( station <= 2 ) {
    
    (layer == 1 ) ? os << "in" : os << "out";
    
      }else if( station > 2 )
      {
    if( sector == 4 && station == 4 )
    {
      if( subsector == 1 )
      {
        os << "--";
      }
      else if( subsector == 2 )
      {
        os << "-";
      }
      else if( subsector == 3 )
      {
        os << "+";
      }
      else if( subsector == 4 )
      {
        os <<"++";
      }
    }
      
    if( station == 3 )
    {
      if( subsector == 1 )
        os << "-";
      else
        os << "+";
    }
    else if( station == 4
         && sector != 9
         && sector != 11
         && sector != 4 )
    {
      if( subsector == 1 )
        os << "-";
      else
        os << "+";
    }
      }
    
      os << "_";
      os << "S" << std::setw(2) << std::setfill('0')
     << sector << std::setfill(' ');
    }
    else
    {
      os << "RE"; 
    
      os << std::setw(2) << std::setfill('+') << station * region
     << std::setfill(' ') << "_";

      os << "R" << ring;
      os << "_CH" << std::setw(2) << std::setfill('0') << this->segment();
    }
    _cn = os.str();
  }
  return _cn;
}

std::string 
RPCGeomServ::shortname( void )
{
  if( _sn.empty() )
  {
    int station = _id->station();
    int region = _id->region();
    int roll = _id->roll();
    int ring = _id->ring();
    int layer = _id->layer();
    int sector = _id->sector();
    int subsector = _id->subsector();

    std::stringstream os;

    if( region == 0 )
    {
      os << "RB" << station;
      if( station <= 2 ){

    (layer == 1 ) ? os << "in" : os << "out";

      }
      else
      {
    if( sector == 4 && station == 4 )
    {
      if( subsector == 1 )
      {
        os << "--";
      }
      else if( subsector == 2 )
      {
        os << ",-";
      }
      else if( subsector == 3 )
      {
        os << "+";
      }
      else if( subsector == 4 )
      {
        os << "++";
      }
    }
    else
    {
      if( subsector == 1 )
        os << ",-";
      else
        os << "+";
    }
      }
      if( roll == 1 )
    os << " B";
      else if( roll == 3 )
    os << " F";
      else if( roll == 2 )
    os << " M";
    }
    else
    {
      os << "Ri" << ring << " Su" << subsector;
    }
    _sn = os.str();
  }
  return _sn;
}

// returns a vector with number of channels for each chip in each FEB
std::vector<int>
RPCGeomServ::channelInChip( void )
{
  std::vector<int> chipCh(4,8);//Endcap
  
  if(_id->region()==0){//Barrel
    chipCh.clear();

    int station = _id->station();
  
    if (station<3 && _id->layer()==1){ // i.e. RB1in ||RB2in  
      chipCh.emplace_back(7);
      chipCh.emplace_back(8);
    }else if (station == 1 || station == 3){//i.e. RB1out || RB3 
      chipCh.emplace_back(7);
      chipCh.emplace_back(7);
    }else if (station == 2){// i.e. RB2out
      chipCh.emplace_back(6);
      chipCh.emplace_back(8);
    }else if (_id->sector() == 4 || _id->sector()==10 ||(_id->sector() == 8 &&  _id->subsector()!=1) || (_id->sector() == 12  &&  _id->subsector()==1)){
      chipCh.emplace_back(6);//i.e. Sector 4 &  10 RB4 and Sector 8 &12 RB4+
      chipCh.emplace_back(6);
    }else {
      chipCh.emplace_back(8);
      chipCh.emplace_back(8);
    }   
  }

  return chipCh;
}


int 
RPCGeomServ::eta_partition()
{
  if (_t<-90){
    if (_id->region() == 0 ){
      if (this->inverted()) {
    _t = 3*(_id->ring())+ (3-_id->roll())-1;
      }else{
    _t = 3*(_id->ring())+ _id->roll()-2;
      }
    }else{
      _t = _id->region() * (3*(3-_id->ring()) + _id->roll() + 7);
    }
  }
  return _t;
} 

int
RPCGeomServ::chambernr()
{

  // Station1
  if( _id->station() ==1) {
    
    // in
    if(_id->layer() ==1) { 
      
      if(_id->roll()==1) 
    _cnr = 1;
      else 
    _cnr = 2;
    }
    //out
    else 
      {
    if(_id->roll()==1) 
      _cnr = 3;
    else 
      _cnr = 4;
    
      }
  }


  //Station 2
  if (_id->station()==2) {
   
    //in
    if(_id->layer()==1) {
      
      if(_id->roll()==1)//backward
    _cnr = 5;
      if(_id->roll()==3)//forward
    _cnr=6;
      if(_id->roll()==2)//middle
    _cnr=7;
    }
    //out
    else {
      
      if(_id->roll()==2)
    
    _cnr=7;

      if(_id->roll()==1)
    _cnr=8;
      if(_id->roll()==3)
    _cnr=9;
    
    }
  }
    
  //RB3- RB3+
  if(_id->station()==3)
    {
      if(_id->subsector()==1) {
    
    if(_id->roll()==1)
      _cnr=10;
    else 
      _cnr=11;
      }
      else {
    
    if(_id->roll()==1)
      _cnr=12;
    else
      _cnr=13;
      }
      
    }

  //RB4
  if(_id->station()==4) {
    
    if (_id->sector()== 4) {
      
      if ( _id->subsector()==2){//RB4-
    
    if(_id->roll()==1)
      _cnr=14;
    else
      _cnr=15;
    
      }
      
      if ( _id->subsector()==3){//RB4+
    
    if(_id->roll()==1)
      _cnr=16;
    else
      _cnr=17;
    
      }
      
      if ( _id->subsector()==1) {//RB4--
    
    if(_id->roll()==1)
      _cnr=18;
    else
      _cnr=19;
      }
      
      if ( _id->subsector()==4){//RB4++
    
    if(_id->roll()==1)
      _cnr=20;
    else
      _cnr=21;
    
      }
      
    }  
    
    else 
      
      {
    if(_id->subsector()==1) {
      
      if(_id->roll()==1)
        _cnr=14;
      else 
        _cnr=15;
    }
    else {
      
      if(_id->roll()==1)
        _cnr=16;
      else
        _cnr=17;
    } 
      } 
  }
  

  // _cnr=10;
  return _cnr;
  
}

int 
RPCGeomServ::segment( void )
{
  int nsub = 6;
  int station = _id->station();
  int ring = _id->ring();
  ( ring == 1 && station > 1 ) ? nsub = 3 : nsub = 6;

  return( _id->subsector() + nsub * ( _id->sector() - 1 ));
}

bool
RPCGeomServ::inverted()
{
  // return !(this->zpositive() && this->aclockwise());
  return !(this->zpositive());
}


bool
RPCGeomServ::zpositive()
{
  if (_id->region()==0 && _t<-90 ){
    if (_id->ring()<0){
      _z=false;
    }
    if (_id->ring()==0){
      if (_id->sector() == 1 || _id->sector() == 4 ||
      _id->sector() == 5 || _id->sector() == 8 ||
      _id->sector() == 9 || _id->sector() == 12){
    _z=false;
      }
    } 
  }
 
  return _z;
}

bool
RPCGeomServ::aclockwise()
{
  if (_id->region()==0 && _t<-90 ){
    if (_id->ring() > 0){
      if (_id->layer()==2){
    _a=false;
      }
    }else if(_id->ring() <0){
      if (_id->layer()==1){
    _a=false;
      }
    }else if(_id->ring() ==0) {
      if ((_id->sector() == 1 || _id->sector() == 4 ||
       _id->sector() == 5 || _id->sector() == 8 ||
       _id->sector() == 9 || _id->sector() == 12) && _id->layer()==1)
    _a=false;
      else if ((_id->sector() == 2 || _id->sector() == 3 ||
        _id->sector() == 6 || _id->sector() == 7 ||
        _id->sector() == 10|| _id->sector() == 11) && _id->layer()==2)
    _a=false;
    }
  }
  return _a;
}

RPCGeomServ::RPCGeomServ() : _id(nullptr), _n(""), _sn(""), _cn(""), _t (-99), _z(false), _a(false)
{}