LVL1TGCTrigger::TGCNSWCoincidenceMap Class Reference

#include <TGCNSWCoincidenceMap.h>

Inheritance diagram for LVL1TGCTrigger::TGCNSWCoincidenceMap:

Collaboration diagram for LVL1TGCTrigger::TGCNSWCoincidenceMap:

Public Member Functions
	TGCNSWCoincidenceMap (TGCArguments *tgcargs, const std::string &version, int side, int oct, int mod)
	~TGCNSWCoincidenceMap ()=default
bool	isForward (int module)
int	TGCNSW_pTcalcu_EtaPhi (const LVL1TGC::NSWTrigOut *nswOut, int RoI) const
int	TGCNSW_pTcalcu_EtaDtheta (const LVL1TGC::NSWTrigOut *nswOut, int RoI) const
const std::string &	getVersion () const
int	getSideId () const
int	getOctantId () const
	TGCNSWCoincidenceMap (const TGCNSWCoincidenceMap &right)
bool	readMap (const std::string &moduleName, ReadCW_Type cw_type)
bool	readShift ()
TGCArguments *	tgcArgs ()
const TGCArguments *	tgcArgs () const
bool	msgLvl (const MSG::Level lvl) const
	Test the output level.
MsgStream &	msg () const
	The standard message stream.
MsgStream &	msg (const MSG::Level lvl) const
	The standard message stream.
void	setLevel (MSG::Level lvl)
	Change the current logging level.

Private Types
enum	{ N_dEta =64 , N_dPhi =16 , N_Dtheta =32 , N_PT_THRESH =15 }
enum	ReadCW_Type { EtaPhi_CW =0 , EtaDtheta_CW }

Private Member Functions
	TGCNSWCoincidenceMap ()=delete
void	initMessaging () const
	Initialize our message level and MessageSvc.

Private Attributes
std::map< TGCRegionType, int >	m_NumberOfEtaRaw = {{TGCRegionType::ENDCAP,37}, {TGCRegionType::FORWARD,16}}
std::map< TGCRegionType, int >	m_NumberOfRoI = {{TGCRegionType::ENDCAP, LVL1TGC::kNumberOfEndcapRoI}, {TGCRegionType::FORWARD,64}}
std::vector< short int >	m_EtaPhi_CW [N_dEta][N_dPhi]
std::vector< short int >	m_EtaDtheta_CW [N_dEta][N_Dtheta]
std::vector< short int >	m_Offset_Eta
std::vector< short int >	m_Offset_Phi
std::string	m_verName
int	m_side {}
int	m_octant {}
int	m_module {}
int	m_sector {}
TGCRegionType	m_region {TGCRegionType::FORWARD}
TGCArguments *	m_tgcArgs {}
std::string	m_nm
	Message source name.
boost::thread_specific_ptr< MsgStream >	m_msg_tls
	MsgStream instance (a std::cout like with print-out levels)
std::atomic< IMessageSvc * >	m_imsg { nullptr }
	MessageSvc pointer.
std::atomic< MSG::Level >	m_lvl { MSG::NIL }
	Current logging level.
std::atomic_flag m_initialized	ATLAS_THREAD_SAFE = ATOMIC_FLAG_INIT
	Messaging initialized (initMessaging)

Detailed Description

Definition at line 24 of file TGCNSWCoincidenceMap.h.

Member Enumeration Documentation

anonymous enum

anonymous enum

private

Enumerator
N_dEta
N_dPhi
N_Dtheta
N_PT_THRESH

Definition at line 28 of file TGCNSWCoincidenceMap.h.

{N_dEta=64,N_dPhi=16,N_Dtheta=32,N_PT_THRESH=15};//6bit eta,4bit phi,5bit Dtheta, 4bit pT

ReadCW_Type

enum LVL1TGCTrigger::TGCNSWCoincidenceMap::ReadCW_Type

private

Enumerator
EtaPhi_CW
EtaDtheta_CW

Definition at line 29 of file TGCNSWCoincidenceMap.h.

{EtaPhi_CW=0,EtaDtheta_CW};

Constructor & Destructor Documentation

TGCNSWCoincidenceMap() [1/3]

LVL1TGCTrigger::TGCNSWCoincidenceMap::TGCNSWCoincidenceMap ( )

privatedelete

TGCNSWCoincidenceMap() [2/3]

LVL1TGCTrigger::TGCNSWCoincidenceMap::TGCNSWCoincidenceMap	(	TGCArguments *	tgcargs,
		const std::string &	version,
		int	side,
		int	oct,
		int	mod )

Definition at line 22 of file TGCNSWCoincidenceMap.cxx.

    :AthMessaging("TGCNSWCoincidenceMap"),
     m_verName(version),
     m_side(side),
     m_octant(oct),
     m_module(mod),
     m_tgcArgs(tgcargs)
  {
 
    setLevel(tgcArgs()->MSGLEVEL());
 
    if(!tgcArgs()->USE_NSW()){return;}
 
    std::string moduleName = std::to_string(mod);
    if( isForward(mod)){
      m_region=TGCRegionType::FORWARD;
      m_sector  = m_module/3; 
      m_sector += 3*m_octant;
    }
    else{
      m_region=TGCRegionType::ENDCAP;
      m_sector  = m_module%3 + 2*(m_module/3);
      m_sector += 6*m_octant;
    }
 
 
    m_Offset_Eta.resize(m_NumberOfRoI[m_region]);
    m_Offset_Phi.resize(m_NumberOfRoI[m_region]);
 
    for(int dr=0;dr!=N_dEta;dr++){
      for(int dphi=0;dphi!=N_dPhi;dphi++){
    m_EtaPhi_CW[dr][dphi].resize(m_NumberOfRoI[m_region]);
    std::fill(m_EtaPhi_CW[dr][dphi].begin(), m_EtaPhi_CW[dr][dphi].end(),0);
      }
      for(int dTheta=0;dTheta!=N_Dtheta;dTheta++){
    m_EtaDtheta_CW[dr][dTheta].resize(m_NumberOfRoI[m_region]);
    std::fill(m_EtaDtheta_CW[dr][dTheta].begin(), m_EtaDtheta_CW[dr][dTheta].end(),0);
      }
    }
 
    if( ( oct%2==0 && mod==5 ) || (oct%2==1 && (mod==2 || mod==8)) ){moduleName+="b";}
    else{moduleName+="a";}
 
 
    //---------Read out CW data---------
    if(!this->readMap( moduleName, ReadCW_Type::EtaPhi_CW) || 
       !this->readMap( moduleName, ReadCW_Type::EtaDtheta_CW)
       ){
      // NSW trigger flag is set to false when the map reading failed.
      tgcArgs()->set_USE_NSW(false);
    }
 
    //---------Fill Shift data---------
    if(!this->readShift()){ 
      // NSW trigger flag is set to false when the map reading failed.
      tgcArgs()->set_USE_NSW(false);
    }
 
  }

~TGCNSWCoincidenceMap()

LVL1TGCTrigger::TGCNSWCoincidenceMap::~TGCNSWCoincidenceMap ( )

default

TGCNSWCoincidenceMap() [3/3]

LVL1TGCTrigger::TGCNSWCoincidenceMap::TGCNSWCoincidenceMap

(

const TGCNSWCoincidenceMap &

right

)

Member Function Documentation

getOctantId()

int LVL1TGCTrigger::TGCNSWCoincidenceMap::getOctantId

(

)

const

getSideId()

int LVL1TGCTrigger::TGCNSWCoincidenceMap::getSideId

(

)

const

getVersion()

const std::string & LVL1TGCTrigger::TGCNSWCoincidenceMap::getVersion

(

)

const

initMessaging()

void AthMessaging::initMessaging ( ) const

privateinherited

Initialize our message level and MessageSvc.

This method should only be called once.

Definition at line 39 of file AthMessaging.cxx.

{
  m_imsg = Athena::getMessageSvc();
  // If user did not set an explicit level, set a default
  if (m_lvl == MSG::NIL) {
    m_lvl = m_imsg ?
      static_cast<MSG::Level>( m_imsg.load()->outputLevel(m_nm) ) :
      MSG::INFO;
  }
}

isForward()

bool LVL1TGCTrigger::TGCNSWCoincidenceMap::isForward

(

int

module

)

Definition at line 83 of file TGCNSWCoincidenceMap.cxx.

                                               {
    return ( mod==2 || mod==5 || mod==8 );
  }

msg() [1/2]

MsgStream & AthMessaging::msg ( ) const

inlineinherited

The standard message stream.

Returns a reference to the default message stream May not be invoked before sysInitialize() has been invoked.

Definition at line 167 of file AthMessaging.h.

{
  MsgStream* ms = m_msg_tls.get();
  if (!ms) {
    if (!m_initialized.test_and_set()) initMessaging();
    ms = new MsgStream(m_imsg,m_nm);
    m_msg_tls.reset( ms );
  }
 
  ms->setLevel (m_lvl);
  return *ms;
}

msg() [2/2]

MsgStream & AthMessaging::msg ( const MSG::Level lvl ) const

inlineinherited

The standard message stream.

Returns a reference to the default message stream May not be invoked before sysInitialize() has been invoked.

Definition at line 182 of file AthMessaging.h.

{ return msg() << lvl; }

msgLvl()

bool AthMessaging::msgLvl ( const MSG::Level lvl ) const

inlineinherited

Test the output level.

Parameters

lvl	The message level to test against

Returns

boolean Indicating if messages at given level will be printed

Return values

true	Messages at level "lvl" will be printed

Definition at line 151 of file AthMessaging.h.

{
  // If user did not set explicit message level we have to initialize
  // the messaging and retrieve the default via the MessageSvc.
  if (m_lvl==MSG::NIL && !m_initialized.test_and_set()) initMessaging();
 
  if (m_lvl <= lvl) {
    msg() << lvl;
    return true;
  } else {
    return false;
  }
}

readMap()

bool LVL1TGCTrigger::TGCNSWCoincidenceMap::readMap	(	const std::string &	moduleName,
		ReadCW_Type	cw_type )

Definition at line 129 of file TGCNSWCoincidenceMap.cxx.

  {
    std::string kSide[2] = {"a", "c"};
    std::string kCWtype[2] = {"EtaPhi","EtaDtheta"};
 
    std::string dbname="";
 
    dbname = "/NSW/cm_" + kSide[m_side] + moduleName +kCWtype[cw_type]+"_Octant_"+m_verName+".db";
    std::string fullName = PathResolver::FindCalibDirectory("dev")+"/TrigT1TGC"+dbname;
 
    std::ifstream data(fullName);
    if(!data.is_open()){
      ATH_MSG_WARNING("Cannot open file dev/TrigT1TGC" << dbname);
      return false;
    }
    char delimiter = '\n';
    std::string field;
    std::string tag;
 
    while (std::getline(data, field, delimiter)) {
      int roi=-1;
      unsigned int n_Etabit=0;
      unsigned int n_Phibit=0;
      std::istringstream header(field); 
      header >> tag;
      if(tag=="#"){ // read header part.     
    header >> roi >> n_Etabit >> n_Phibit;
      }
      // get trigger word
      std::string word;
      unsigned int pT=0;

      for(size_t posR=0; posR<N_dEta ; posR++){
        std::getline(data, field, delimiter);
        std::istringstream cont(field);
        //----Read EtaPhi CW---------
        if(cw_type==ReadCW_Type::EtaPhi_CW){
          for(size_t posPHI=0; posPHI<N_dPhi; posPHI++){
            cont >> word;
            std::istringstream(word) >> std::hex >> pT;
            //values should be checked against sensible limits; do minimal check here:
            if (roi<0) continue;
            //coverity[TAINTED_SCALAR]
            m_EtaPhi_CW[posR][posPHI][roi]=pT;
          }
        }
        //-----Read EtaDtheta CW-----
        if(cw_type==ReadCW_Type::EtaDtheta_CW){
          for(size_t posDTHETA=0; posDTHETA<N_Dtheta; posDTHETA++){
            cont >> word;
            std::istringstream(word) >> std::hex >> pT;
            auto & roiIndexedVector = m_EtaDtheta_CW[posR][posDTHETA];
            if ((roi<0) or (static_cast<size_t>(roi)>=roiIndexedVector.size())){
              throw std::out_of_range("roi outside of vector limits in TGCNSWCoincidenceMap::readMap");
            }
            roiIndexedVector[roi]=pT;
          }
        }
      }      
    }
    data.close();
 
 
    return true;
  } 

readShift()

bool LVL1TGCTrigger::TGCNSWCoincidenceMap::readShift

(

)

Definition at line 197 of file TGCNSWCoincidenceMap.cxx.

                                      {
 
    char delimiter = '\n';
    std::string field;
    std::string tag;
    int side;
    int triggerSector;
 
    //------- Read Endcap Shift
    std::string dbname="/NSW/";
 
    if(m_region==TGCRegionType::ENDCAP){dbname += "RoIpos_Endcap.db";}
    if(m_region==TGCRegionType::FORWARD){dbname += "RoIpos_Forward.db";}
    std::string fullName = PathResolver::FindCalibDirectory("dev")+"/TrigT1TGC"+dbname;
 
    std::ifstream data(fullName);
    if(!data.is_open()){return false;}
    while (std::getline(data, field, delimiter)) {
      std::istringstream header(field); 
      header >> tag;
      if(tag=="#"){ // read header part.
    header >> side >> triggerSector ;
      }
 
      // get trigger word
      std::getline(data, field, delimiter);
      std::istringstream cont(field);
      std::string word;
      int shift=0;
      std::array<int, 4> phi_shift {};
      std::array<int, 37> eta_shift {};
 
      if(side!=m_side || triggerSector!=m_sector){ continue; }
      for(int phiN=0; phiN!=4; phiN++){
    cont >> word;
    std::istringstream(word) >> shift;
    phi_shift[phiN]=shift;
      }
      for(int etaN=0; etaN!=m_NumberOfEtaRaw[m_region]; etaN++){
    cont >> word;
    std::istringstream(word) >> shift;
    eta_shift[etaN]=shift;
      }
      for(int roi=0;roi!=m_NumberOfRoI[m_region];roi++){
    m_Offset_Eta[roi]=eta_shift[roi/4]+N_dEta/2;
    m_Offset_Phi[roi]=phi_shift[roi%4]+N_dPhi/2;
      }
      
      break;
    }
    data.close();
 
    return true;
  }

setLevel()

void AthMessaging::setLevel ( MSG::Level lvl )

inherited

Change the current logging level.

Use this rather than msg().setLevel() for proper operation with MT.

Definition at line 28 of file AthMessaging.cxx.

{
  m_lvl = lvl;
}

tgcArgs() [1/2]

TGCArguments * LVL1TGCTrigger::TGCNSWCoincidenceMap::tgcArgs ( )

inline

Definition at line 51 of file TGCNSWCoincidenceMap.h.

{ return m_tgcArgs;}

tgcArgs() [2/2]

const TGCArguments * LVL1TGCTrigger::TGCNSWCoincidenceMap::tgcArgs ( ) const

inline

Definition at line 52 of file TGCNSWCoincidenceMap.h.

{ return m_tgcArgs;}

TGCNSW_pTcalcu_EtaDtheta()

int LVL1TGCTrigger::TGCNSWCoincidenceMap::TGCNSW_pTcalcu_EtaDtheta	(	const LVL1TGC::NSWTrigOut *	nswOut,
		int	RoI ) const

just for now. before we finalize the strategy, we always return the lowest pT value if we can find the NSW hit

Definition at line 110 of file TGCNSWCoincidenceMap.cxx.

  {
    std::vector<uint8_t> nswEta_vec=nswOut->getNSWeta();
    std::vector<uint8_t> nswDtheta_vec=nswOut->getNSWDtheta();
    int highest_pT=0;
 
    for(unsigned int nswTrk_id=0;nswTrk_id!=nswEta_vec.size();nswTrk_id++){
      int eta_decode=m_Offset_Eta[roi]-nswEta_vec[nswTrk_id];
      int dTheta_decode=nswDtheta_vec[nswTrk_id];
      if(eta_decode<0 || eta_decode>=N_dEta){continue;}
      if(highest_pT<m_EtaDtheta_CW[eta_decode][dTheta_decode][roi]){
    highest_pT=m_EtaDtheta_CW[eta_decode][dTheta_decode][roi];
      }
    }
 
    return 1; 
  }

TGCNSW_pTcalcu_EtaPhi()

int LVL1TGCTrigger::TGCNSWCoincidenceMap::TGCNSW_pTcalcu_EtaPhi	(	const LVL1TGC::NSWTrigOut *	nswOut,
		int	RoI ) const

just for now. before we finalize the strategy, we always return the lowest pT value if we can find the NSW hit

Definition at line 88 of file TGCNSWCoincidenceMap.cxx.

  {
    std::vector<uint8_t> nswEta_vec=nswOut->getNSWeta();
    std::vector<uint8_t> nswPhi_vec=nswOut->getNSWphi();
    int highest_pT=0;
 
    for(unsigned int nswTrk_id=0;nswTrk_id!=nswEta_vec.size();nswTrk_id++){
      int eta_decode=m_Offset_Eta[roi]-nswEta_vec[nswTrk_id];
      int phi_decode=m_Offset_Phi[roi]-nswPhi_vec[nswTrk_id];
      if(eta_decode<0 || eta_decode>=N_dEta){continue;}
      if(phi_decode<0 || phi_decode>=N_dPhi){continue;}
 
      if(highest_pT<m_EtaPhi_CW[eta_decode][phi_decode][roi]){
    highest_pT=m_EtaPhi_CW[eta_decode][phi_decode][roi];
      }
    }
    
    return 1; 
  }

Member Data Documentation

ATLAS_THREAD_SAFE

std::atomic_flag m_initialized AthMessaging::ATLAS_THREAD_SAFE = ATOMIC_FLAG_INIT

mutableprivateinherited

Messaging initialized (initMessaging)

Definition at line 141 of file AthMessaging.h.

m_EtaDtheta_CW

std::vector<short int> LVL1TGCTrigger::TGCNSWCoincidenceMap::m_EtaDtheta_CW[N_dEta][N_Dtheta]

private

Definition at line 56 of file TGCNSWCoincidenceMap.h.

m_EtaPhi_CW

std::vector<short int> LVL1TGCTrigger::TGCNSWCoincidenceMap::m_EtaPhi_CW[N_dEta][N_dPhi]

private

Definition at line 55 of file TGCNSWCoincidenceMap.h.

m_imsg

std::atomic<IMessageSvc*> AthMessaging::m_imsg { nullptr }

mutableprivateinherited

MessageSvc pointer.

Definition at line 135 of file AthMessaging.h.

{ nullptr };

m_lvl

std::atomic<MSG::Level> AthMessaging::m_lvl { MSG::NIL }

mutableprivateinherited

Current logging level.

Definition at line 138 of file AthMessaging.h.

{ MSG::NIL };

m_module

int LVL1TGCTrigger::TGCNSWCoincidenceMap::m_module {}

private

Definition at line 63 of file TGCNSWCoincidenceMap.h.

{};

m_msg_tls

boost::thread_specific_ptr<MsgStream> AthMessaging::m_msg_tls

mutableprivateinherited

MsgStream instance (a std::cout like with print-out levels)

Definition at line 132 of file AthMessaging.h.

m_nm

std::string AthMessaging::m_nm

privateinherited

Message source name.

Definition at line 129 of file AthMessaging.h.

m_NumberOfEtaRaw

std::map<TGCRegionType,int> LVL1TGCTrigger::TGCNSWCoincidenceMap::m_NumberOfEtaRaw = {{TGCRegionType::ENDCAP,37}, {TGCRegionType::FORWARD,16}}

private

Definition at line 30 of file TGCNSWCoincidenceMap.h.

{{TGCRegionType::ENDCAP,37}, {TGCRegionType::FORWARD,16}};

m_NumberOfRoI

std::map<TGCRegionType,int> LVL1TGCTrigger::TGCNSWCoincidenceMap::m_NumberOfRoI = {{TGCRegionType::ENDCAP, LVL1TGC::kNumberOfEndcapRoI}, {TGCRegionType::FORWARD,64}}

private

Definition at line 31 of file TGCNSWCoincidenceMap.h.

{{TGCRegionType::ENDCAP, LVL1TGC::kNumberOfEndcapRoI}, {TGCRegionType::FORWARD,64}};

m_octant

int LVL1TGCTrigger::TGCNSWCoincidenceMap::m_octant {}

private

Definition at line 62 of file TGCNSWCoincidenceMap.h.

{};

m_Offset_Eta

std::vector<short int> LVL1TGCTrigger::TGCNSWCoincidenceMap::m_Offset_Eta

private

Definition at line 57 of file TGCNSWCoincidenceMap.h.

m_Offset_Phi

std::vector<short int> LVL1TGCTrigger::TGCNSWCoincidenceMap::m_Offset_Phi

private

Definition at line 58 of file TGCNSWCoincidenceMap.h.

m_region

TGCRegionType LVL1TGCTrigger::TGCNSWCoincidenceMap::m_region {TGCRegionType::FORWARD}

private

Definition at line 65 of file TGCNSWCoincidenceMap.h.

{TGCRegionType::FORWARD};

m_sector

int LVL1TGCTrigger::TGCNSWCoincidenceMap::m_sector {}

private

Definition at line 64 of file TGCNSWCoincidenceMap.h.

{};

m_side

int LVL1TGCTrigger::TGCNSWCoincidenceMap::m_side {}

private

Definition at line 61 of file TGCNSWCoincidenceMap.h.

{};

m_tgcArgs

TGCArguments* LVL1TGCTrigger::TGCNSWCoincidenceMap::m_tgcArgs {}

private

Definition at line 67 of file TGCNSWCoincidenceMap.h.

{};

m_verName

std::string LVL1TGCTrigger::TGCNSWCoincidenceMap::m_verName

private

Definition at line 60 of file TGCNSWCoincidenceMap.h.

---

The documentation for this class was generated from the following files:

• TGCNSWCoincidenceMap.h
• TGCNSWCoincidenceMap.cxx