#include <RBCProcessRPCDigis.h>

Inheritance diagram for RBCProcessRPCDigis:

Public Member Functions
void	builddata ()

void	configure ()

int	next () override

void	print_output ()

	RBCProcessRPCDigis (const edm::ESHandle< RPCGeometry > &, const edm::Handle< RPCDigiCollection > &)
	Standard constructor. More...

void	reset ()

RPCInputSignal *	retrievedata () override

void	rewind ()

void	showfirst ()

	~RBCProcessRPCDigis () 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 > &) const

void	setDigiAt (int, int, RPCData &)

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

Private Attributes
std::map< int, RBCInput * >	m_data

const bool	m_debug

std::map< int, l1trigger::Counters >	m_digiCounters

std::unique_ptr< RPCInputSignal >	m_lbin

const int	m_maxBxWindow

const edm::Handle< RPCDigiCollection > *	m_ptr_digiColl

const edm::ESHandle< RPCGeometry > *	m_ptr_rpcGeom

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

Detailed Description

Author: Andres Felipe Osorio Oliveros

Date: 2009-04-15

Definition at line 39 of file RBCProcessRPCDigis.h.

Constructor & Destructor Documentation

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

Standard constructor.

Definition at line 20 of file RBCProcessRPCDigis.cc.

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

                                                                                          :
   m_maxBxWindow{3},
   m_debug{false}
 {
   
   m_ptr_rpcGeom  = & rpcGeom;
   m_ptr_digiColl = & digiColl;
   
   m_lbin =std::make_unique<RBCLinkBoardGLSignal>( &m_data ) ;
     
   configure();
   
 }

RBCProcessRPCDigis::~RBCProcessRPCDigis ( )

override

Destructor.

Definition at line 44 of file RBCProcessRPCDigis.cc.

                                         {
   
 } 

Member Function Documentation

void RBCProcessRPCDigis::builddata ( )

Definition at line 187 of file RBCProcessRPCDigis.cc.

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

Referenced by next().

 {
   
   
   for(auto & vecData: m_vecDataperBx) {
     
     int bx = vecData.first;
     int bxsign(1);
     
     if ( bx != 0 ) bxsign = ( bx / abs(bx) );
     else bxsign = 1;
     
     for(auto & item : vecData.second) {
       
       for(int k=0; k < 6; ++k) {
         
         int 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.emplace( code , signal );
         
       }
     }
   }
   
   if ( m_debug and not m_vecDataperBx.empty() ) std::cout << "builddata: completed. size of data: " << m_data.size() << std::endl;
   
 }

void RBCProcessRPCDigis::configure ( )

Definition at line 35 of file RBCProcessRPCDigis.cc.

References m_digiCounters, rpctechnicaltrigger::s_wheelid, and makeMuonMisalignmentScenario::wheel.

Referenced by RBCProcessRPCDigis().

 {
   for( auto wheel : s_wheelid)
     m_digiCounters.emplace(wheel, wheel);
 }

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

private

Definition at line 222 of file RBCProcessRPCDigis.cc.

Referenced by next(), and 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 RBCProcessRPCDigis::initialize ( std::vector< RPCData > & dataVec ) const

private

Definition at line 159 of file RBCProcessRPCDigis.cc.

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

Referenced by next(), and 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 RBCProcessRPCDigis::next ( void )

overridevirtual

Implements ProcessInputSignal.

Definition at line 49 of file RBCProcessRPCDigis.cc.

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

                              {
   
   //...clean up previous data contents
   
   reset();
   
   int ndigis(0);
   
   for (auto const& detUnit : *(*m_ptr_digiColl) ) {
     
     if ( m_debug ) std::cout << "looping over digis 1 ..." << std::endl;
     
     auto digiItr = detUnit.second.first;
     int bx = (*digiItr).bx();
     
     if ( abs(bx) >= m_maxBxWindow ) {
       if ( m_debug )  std::cout << "RBCProcessRPCDigis> found a bx bigger than max allowed: "
                                 << bx << std::endl;
       continue;
     }
     
     const RPCDetId & id  = detUnit.first;
     const RPCRoll * roll = dynamic_cast<const RPCRoll* >( (*m_ptr_rpcGeom)->roll(id));
     
     if((roll->isForward())) {
       if( m_debug ) std::cout << "RBCProcessRPCDigis: 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);
       setDigiAt( sector, digipos,  wheelData[ (wheel + 2) ] );
     }
     else{
       setDigiAt( sector, digipos, (*itr).second[ (wheel + 2) ] );
     }
     
     auto wheelCounter = m_digiCounters.find( wheel );
     
     if ( wheelCounter != m_digiCounters.end() )
       (*wheelCounter).second.incrementSector( sector );
     
     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 << "RBCProcessRPCDigis: DataSize: " << m_data.size() 
                            << " ndigis " << ndigis << std::endl;
   
   for( auto& wheelCounter : m_digiCounters) {
     wheelCounter.second.evalCounters();
     if ( m_debug ) wheelCounter.second.printSummary();
   }
   
   if ( m_data.empty() ) return 0;
   
   return 1;
   
 }

void RBCProcessRPCDigis::print_output ( )

Definition at line 281 of file RBCProcessRPCDigis.cc.

References gather_cfg::cout, and m_data.

Referenced by next().

 {
 
   std::cout << "RBCProcessRPCDigis> Output starts" << std::endl;
   
   for( auto const& item: m_data) {
     std::cout << item.first << '\t' << (* item.second ) << '\n';
   }
 
   std::cout << "RBCProcessRPCDigis> Output ends" << std::endl;
   
 }

void RBCProcessRPCDigis::reset ( void )

Definition at line 153 of file RBCProcessRPCDigis.cc.

References m_vecDataperBx.

Referenced by next().

 {  
   m_vecDataperBx.clear();
 }

RPCInputSignal* RBCProcessRPCDigis::retrievedata ( )

inlineoverridevirtual

Implements ProcessInputSignal.

Definition at line 58 of file RBCProcessRPCDigis.h.

References m_lbin.

                                            {
     return  m_lbin.get();
   };

void RBCProcessRPCDigis::rewind ( )

inline

Definition at line 62 of file RBCProcessRPCDigis.h.

62 {};

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

private

Definition at line 240 of file RBCProcessRPCDigis.cc.

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

Referenced by next(), and 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 RBCProcessRPCDigis::setInputBit	(	std::bitset< 15 > &	signals,
		int	digipos
	)

private

Definition at line 272 of file RBCProcessRPCDigis.cc.

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

Referenced by setDigiAt(), and showfirst().

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