#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. More...

MsgStream &	msg () const
	The standard message stream. More...

MsgStream &	msg (const MSG::Level lvl) const
	The standard message stream. More...

void	setLevel (MSG::Level lvl)
	Change the current logging level. More...

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. More...

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

TGCArguments *	m_tgcArgs

std::string	m_nm
	Message source name. More...

boost::thread_specific_ptr< MsgStream >	m_msg_tls
	MsgStream instance (a std::cout like with print-out levels) More...

std::atomic< IMessageSvc * >	m_imsg { nullptr }
	MessageSvc pointer. More...

std::atomic< MSG::Level >	m_lvl { MSG::NIL }
	Current logging level. More...

std::atomic_flag m_initialized	ATLAS_THREAD_SAFE = ATOMIC_FLAG_INIT
	Messaging initialized (initMessaging) More...

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.

28 {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.

29 {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();
   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 164 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 179 of file AthMessaging.h.

180 { 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 (!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;
         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 195 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.

51 { return m_tgcArgs;}

◆ tgcArgs() [2/2]

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

inline

Definition at line 52 of file TGCNSWCoincidenceMap.h.

52 { 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.

◆ m_lvl

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

mutableprivateinherited

Current logging level.

Definition at line 138 of file AthMessaging.h.

◆ m_module

int LVL1TGCTrigger::TGCNSWCoincidenceMap::m_module

private

Definition at line 64 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.

◆ 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.

◆ m_octant

int LVL1TGCTrigger::TGCNSWCoincidenceMap::m_octant

private

Definition at line 63 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

private

Definition at line 66 of file TGCNSWCoincidenceMap.h.

◆ m_sector

int LVL1TGCTrigger::TGCNSWCoincidenceMap::m_sector

private

Definition at line 65 of file TGCNSWCoincidenceMap.h.

◆ m_side

int LVL1TGCTrigger::TGCNSWCoincidenceMap::m_side

private

Definition at line 62 of file TGCNSWCoincidenceMap.h.

◆ m_tgcArgs

TGCArguments* LVL1TGCTrigger::TGCNSWCoincidenceMap::m_tgcArgs

private

Definition at line 68 of file TGCNSWCoincidenceMap.h.

◆ m_verName

std::string LVL1TGCTrigger::TGCNSWCoincidenceMap::m_verName

private

Definition at line 61 of file TGCNSWCoincidenceMap.h.

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

Public Member Functions

Private Types

Private Member Functions

Private Attributes

Detailed Description

Member Enumeration Documentation

◆ anonymous enum

◆ ReadCW_Type

Constructor & Destructor Documentation

◆ TGCNSWCoincidenceMap() [1/3]

◆ TGCNSWCoincidenceMap() [2/3]

◆ ~TGCNSWCoincidenceMap()

◆ TGCNSWCoincidenceMap() [3/3]

Member Function Documentation

◆ getOctantId()

◆ getSideId()

◆ getVersion()

◆ initMessaging()

◆ isForward()

◆ msg() [1/2]

◆ msg() [2/2]

◆ msgLvl()

◆ readMap()

◆ readShift()

◆ setLevel()

◆ tgcArgs() [1/2]

◆ tgcArgs() [2/2]

◆ TGCNSW_pTcalcu_EtaDtheta()

◆ TGCNSW_pTcalcu_EtaPhi()

Member Data Documentation

◆ ATLAS_THREAD_SAFE

◆ m_EtaDtheta_CW

◆ m_EtaPhi_CW

◆ m_imsg

◆ m_lvl

◆ m_module

◆ m_msg_tls

◆ m_nm

◆ m_NumberOfEtaRaw

◆ m_NumberOfRoI

◆ m_octant

◆ m_Offset_Eta

◆ m_Offset_Phi

◆ m_region

◆ m_sector

◆ m_side

◆ m_tgcArgs

◆ m_verName