RBCProcessRPCSimDigis Class Reference

#include <interface/RBCProcessRPCSimDigis.h>

Inheritance diagram for RBCProcessRPCSimDigis:

Public Member Functions
void	builddata ()
int	next () override
void	print_output ()
	RBCProcessRPCSimDigis ()
	Standard constructor. More...
	RBCProcessRPCSimDigis (const edm::ESHandle< RPCGeometry > &, const edm::Handle< edm::DetSetVector< RPCDigiSimLink > > &)
void	reset ()
RPCInputSignal *	retrievedata () override
void	rewind ()
void	showfirst ()
	~RBCProcessRPCSimDigis () override
	Destructor. More...
Public Member Functions inherited from ProcessInputSignal
virtual	~ProcessInputSignal ()

Private Member Functions
int	getBarrelLayer (const int &, const int &)
void	initialize (std::vector< RPCData > &)
void	setDigiAt (int, int, RPCData &)
void	setInputBit (std::bitset< 15 > &, int)

Private Attributes
std::map< int, RBCInput * >	m_data
bool	m_debug
edm::DetSet< RPCDigiSimLink >::const_iterator	m_digiItr
std::unique_ptr< RPCInputSignal >	m_lbin
edm::DetSetVector< RPCDigiSimLink >::const_iterator	m_linkItr
int	m_maxBxWindow
const edm::Handle< edm::DetSetVector< RPCDigiSimLink > > *	m_ptr_digiSimLink
const edm::ESHandle< RPCGeometry > *	m_ptr_rpcGeom
std::map< int, std::vector< RPCData > >	m_vecDataperBx

Detailed Description

Author

Andres Felipe Osorio Oliveros

Date

2009-09-20

Definition at line 40 of file RBCProcessRPCSimDigis.h.

Constructor & Destructor Documentation

RBCProcessRPCSimDigis::RBCProcessRPCSimDigis

(

)

Standard constructor.

RBCProcessRPCSimDigis::RBCProcessRPCSimDigis	(	const edm::ESHandle< RPCGeometry > &	rpcGeom,
		const edm::Handle< edm::DetSetVector< RPCDigiSimLink > > &	digiSimLink
	)

Definition at line 17 of file RBCProcessRPCSimDigis.cc.

{
  
  m_ptr_rpcGeom  = & rpcGeom;
  m_ptr_digiSimLink = & digiSimLink;
  
  m_lbin = std::make_unique<RBCLinkBoardGLSignal>( &m_data ) ;
  
  m_debug = false;
  m_maxBxWindow = 3;
  
}

RBCProcessRPCSimDigis::~RBCProcessRPCSimDigis ( )

override

Destructor.

Definition at line 34 of file RBCProcessRPCSimDigis.cc.

References reset().

Referenced by showfirst().

                                              {
  reset();
} 

Member Function Documentation

void RBCProcessRPCSimDigis::builddata

(

)

Definition at line 180 of file RBCProcessRPCSimDigis.cc.

References funct::abs(), gather_cfg::cout, gen::k, and RBCInput::needmapping.

{
  
  int code(0);
  
  for(auto& dataPerBx: m_vecDataperBx) {
    
    int bx = dataPerBx.first;
    
    int bxsign;
    if ( bx != 0 ) bxsign = ( bx / abs(bx) );
    else bxsign = 1;
    
    for(auto& item : dataPerBx.second) {
      
      for(int k=0; k < 6; ++k) {
        
        code = bxsign * ( 1000000*abs(bx)
                          + 10000*item.wheelIdx()
                          + 100  *item.m_sec1[k]
                          + 1    *item.m_sec2[k] );

        
        RBCInput * signal = & item.m_orsignals[k];
        signal->needmapping = false;
        
        if ( signal->hasData )
          m_data.insert( std::make_pair( code , signal) );
        
      }
    }
  }
  
  if ( m_debug and not m_vecDataperBx.empty()) std::cout << "builddata: completed. size of data: " << m_data.size() << std::endl;
  
}

int RBCProcessRPCSimDigis::getBarrelLayer ( const int &  _layer, const int &  _station  )

private

Definition at line 217 of file RBCProcessRPCSimDigis.cc.

Referenced by showfirst().

{
  
  //... Calculates the generic Barrel Layer (1 to 6)
  int blayer(0);
  
  if ( _station < 3 ) {
    blayer = ( (_station-1) * 2 ) + _layer;
  }
  else {
    blayer = _station + 2;
  }
  
  return blayer;
  
}

void RBCProcessRPCSimDigis::initialize ( std::vector< RPCData > &  dataVec )

private

Definition at line 151 of file RBCProcessRPCSimDigis.cc.

References groupFilesInBlocks::block, constexpr, gather_cfg::cout, mps_fire::i, rpctechnicaltrigger::s_sec1id, rpctechnicaltrigger::s_sec2id, and rpctechnicaltrigger::s_wheelid.

Referenced by showfirst().

{
  
  if ( m_debug ) std::cout << "initialize" << std::endl;
  
  constexpr int maxWheels = 5;
  constexpr int maxRbcBrds = 6;
  
  dataVec.reserve(maxWheels);
  for(int i=0; i < maxWheels; ++i) {
    
    auto& block = dataVec.emplace_back();
    
    block.m_wheel = s_wheelid[i];
    
    for(int j=0; j < maxRbcBrds; ++j) {
      block.m_sec1[j] = s_sec1id[j];
      block.m_sec2[j] = s_sec2id[j];
      block.m_orsignals[j].input_sec[0].reset();
      block.m_orsignals[j].input_sec[1].reset();
      block.m_orsignals[j].needmapping = false;
      block.m_orsignals[j].hasData = false;
    }
  }
  
  if ( m_debug ) std::cout << "initialize: completed" << std::endl;
  
}

int RBCProcessRPCSimDigis::next ( void  )

overridevirtual

Implements ProcessInputSignal.

Definition at line 39 of file RBCProcessRPCSimDigis.cc.

References funct::abs(), groupFilesInBlocks::block, gather_cfg::cout, RPCRoll::id(), triggerObjects_cff::id, initialize(), RPCRoll::isForward(), RPCDetId::layer(), reset(), RPCDetId::ring(), RPCDetId::sector(), relativeConstraints::station, RPCDetId::station(), and makeMuonMisalignmentScenario::wheel.

                                {
  
  //...clean up previous data contents
  
  reset();
  
  int ndigis(0);

  for( m_linkItr = (*m_ptr_digiSimLink)->begin();
       m_linkItr != (*m_ptr_digiSimLink)->end();
       ++m_linkItr ) {
    
    for ( m_digiItr = m_linkItr->data.begin();
          m_digiItr != m_linkItr->data.end();
          ++m_digiItr ) {
      
      if ( m_debug ) std::cout << "looping over digis 1 ..." << std::endl;
      
      int bx = (*m_digiItr).getBx();
      
      if ( abs(bx) >= m_maxBxWindow ) {
        if ( m_debug )  std::cout << "RBCProcessRPCSimDigis> found a bx bigger than max allowed: "
                                  << bx << std::endl;
        continue;
      }
      
      uint32_t detid = m_digiItr->getDetUnitId();
      const RPCDetId id( detid );
      const RPCRoll * roll = dynamic_cast<const RPCRoll* >( (*m_ptr_rpcGeom)->roll(id));
      
      if((roll->isForward())) {
        if( m_debug ) std::cout << "RBCProcessRPCSimDigis: roll is forward" << std::endl;
        continue;
      }
      
      int wheel   = roll->id().ring();                    // -2,-1,0,+1,+2
      int sector  = roll->id().sector();                  // 1 to 12 
      int layer   = roll->id().layer();                   // 1,2
      int station = roll->id().station();                 // 1-4
      int blayer  = getBarrelLayer( layer, station );     // 1 to 6
      int rollid  = id.roll();
      
      int digipos = (station * 100) + (layer * 10) + rollid;
      
      if ( (wheel == -1 || wheel == 0 || wheel == 1) && station == 2 && layer == 1 )
        digipos = 30000 + digipos;
      if ( (wheel == -2 || wheel == 2) && station == 2 && layer == 2 )
        digipos = 30000 + digipos;
      
      if ( (wheel == -1 || wheel == 0 || wheel == 1) && station == 2 && layer == 2 )
        digipos = 20000 + digipos;
      if ( (wheel == -2 || wheel == 2) && station == 2 && layer == 1 )
        digipos = 20000 + digipos;
      
      if ( m_debug ) std::cout << "Bx: "      << bx      << '\t'
                               << "Wheel: "   << wheel   << '\t'
                               << "Sector: "  << sector  << '\t'
                               << "Station: " << station << '\t'
                               << "Layer: "   << layer   << '\t'
                               << "B-Layer: " << blayer  << '\t'
                               << "Roll id: " << rollid  << '\t'
                               << "Digi at: " << digipos << '\n';
      
      //... Construct the RBCinput objects
      auto itr = m_vecDataperBx.find( bx );
      
      if ( itr == m_vecDataperBx.end() ) {
        if ( m_debug ) std::cout << "Found a new Bx: " << bx << std::endl;
        auto& wheelData = m_vecDataperBx[bx];
        initialize(wheelData);
        auto& block = wheelData[ (wheel + 2) ];
        setDigiAt( sector, digipos, block );
      }
      else{
        auto& block = (*itr).second[ (wheel + 2) ];
        setDigiAt( sector, digipos, block );
      }
      
      if ( m_debug ) std::cout << "looping over digis 2 ..." << std::endl;
      
      ++ndigis;
      
    }
    
  }
  
  if ( m_debug ) std::cout << "size of data vectors: " << m_vecDataperBx.size() << std::endl;
  
  builddata();
  
  if ( m_debug ) {
    std::cout << "after reset" << std::endl;
    print_output();
  }
  
  if ( m_debug ) std::cout << "RBCProcessRPCSimDigis: DataSize: " << m_data.size() 
                           << " ndigis " << ndigis << std::endl;
  
  if ( m_data.empty() ) return 0;
  
  return 1;
  
}

void RBCProcessRPCSimDigis::print_output

(

)

Definition at line 276 of file RBCProcessRPCSimDigis.cc.

References gather_cfg::cout.

{

  std::cout << "RBCProcessRPCSimDigis> Output starts" << std::endl;
  
  std::map<int,RBCInput*>::const_iterator itr;
  for( itr = m_data.begin(); itr != m_data.end(); ++itr) {
    std::cout << (*itr).first << '\t' << (* (*itr).second ) << '\n';
  }

  std::cout << "RBCProcessRPCSimDigis> Output ends" << std::endl;
  
}

void RBCProcessRPCSimDigis::reset

(

void

)

Definition at line 143 of file RBCProcessRPCSimDigis.cc.

{
  
  m_vecDataperBx.clear();
  
}

RPCInputSignal* RBCProcessRPCSimDigis::retrievedata ( )

inlineoverridevirtual

Implements ProcessInputSignal.

Definition at line 56 of file RBCProcessRPCSimDigis.h.

References m_lbin.

                                           {
    return  m_lbin.get();
  };

void RBCProcessRPCSimDigis::rewind ( )

inline

Definition at line 60 of file RBCProcessRPCSimDigis.h.

{};

void RBCProcessRPCSimDigis::setDigiAt ( int  sector, int  digipos, RPCData &  block  )

private

Definition at line 235 of file RBCProcessRPCSimDigis.cc.

References gather_cfg::cout, spr::find(), RPCData::m_orsignals, RPCData::m_sec1, RPCData::m_sec2, rpctechnicaltrigger::s_sec1id, and rpctechnicaltrigger::s_sec2id.

Referenced by showfirst().

{
  
  int pos   = 0;
  int isAoB = 0;

  if ( m_debug ) std::cout << "setDigiAt" << std::endl;
  
  auto itr = std::find( s_sec1id.begin(), s_sec1id.end(), sector );
  
  if ( itr == s_sec1id.end()) {
    itr = std::find( s_sec2id.begin(), s_sec2id.end(), sector );
    isAoB = 1;
  } 
  
  for ( pos = 0; pos < 6; ++pos ) {
    if (block.m_sec1[pos] == sector || block.m_sec2[pos] == sector )
      break;
  }
  
  if ( m_debug ) std::cout << block.m_orsignals[pos];
  
  setInputBit( block.m_orsignals[pos].input_sec[ isAoB ] , digipos );
  
  block.m_orsignals[pos].hasData = true;
  
  if ( m_debug ) std::cout << block.m_orsignals[pos];
  
  if ( m_debug ) std::cout << "setDigiAt completed" << std::endl;
  
}

void RBCProcessRPCSimDigis::setInputBit ( std::bitset< 15 > &  signals, int  digipos  )

private

Definition at line 267 of file RBCProcessRPCSimDigis.cc.

References gather_cfg::cout, and rpctechnicaltrigger::s_layermap.

Referenced by showfirst().

{
  
  int bitpos = s_layermap.at(digipos);
  if( m_debug ) std::cout << "Bitpos: " << bitpos << std::endl;
  signals.set( bitpos , true );
  
}

void RBCProcessRPCSimDigis::showfirst ( )

inline

Definition at line 61 of file RBCProcessRPCSimDigis.h.

References getBarrelLayer(), initialize(), setDigiAt(), setInputBit(), and ~RBCProcessRPCSimDigis().

{};

Member Data Documentation

std::map<int, RBCInput*> RBCProcessRPCSimDigis::m_data

private

Definition at line 85 of file RBCProcessRPCSimDigis.h.

bool RBCProcessRPCSimDigis::m_debug

private

Definition at line 89 of file RBCProcessRPCSimDigis.h.

edm::DetSet<RPCDigiSimLink>::const_iterator RBCProcessRPCSimDigis::m_digiItr

private

Definition at line 81 of file RBCProcessRPCSimDigis.h.

std::unique_ptr<RPCInputSignal> RBCProcessRPCSimDigis::m_lbin

private

Definition at line 83 of file RBCProcessRPCSimDigis.h.

Referenced by retrievedata().

edm::DetSetVector<RPCDigiSimLink>::const_iterator RBCProcessRPCSimDigis::m_linkItr

private

Definition at line 80 of file RBCProcessRPCSimDigis.h.

int RBCProcessRPCSimDigis::m_maxBxWindow

private

Definition at line 90 of file RBCProcessRPCSimDigis.h.

const edm::Handle<edm::DetSetVector<RPCDigiSimLink> >* RBCProcessRPCSimDigis::m_ptr_digiSimLink

private

Definition at line 78 of file RBCProcessRPCSimDigis.h.

const edm::ESHandle<RPCGeometry>* RBCProcessRPCSimDigis::m_ptr_rpcGeom

private

Definition at line 77 of file RBCProcessRPCSimDigis.h.

std::map<int, std::vector<RPCData> > RBCProcessRPCSimDigis::m_vecDataperBx

private

Definition at line 87 of file RBCProcessRPCSimDigis.h.