sTGCSensitiveDetector Class Reference

#include <sTGCSensitiveDetector.h>

Inheritance diagram for sTGCSensitiveDetector:

Collaboration diagram for sTGCSensitiveDetector:

Public Member Functions
	sTGCSensitiveDetector (const std::string &name, const std::string &hitCollectionName, unsigned baseDepth)
	construction/destruction More...
	~sTGCSensitiveDetector ()=default
void	Initialize (G4HCofThisEvent *HCE) override final
	member functions More...
G4bool	ProcessHits (G4Step *aStep, G4TouchableHistory *ROhist) override final
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 Member Functions
	FRIEND_TEST (sTGCSensitiveDetectortest, Initialize)
	FRIEND_TEST (sTGCSensitiveDetectortest, ProcessHits)
void	initMessaging () const
	Initialize our message level and MessageSvc. More...

Private Attributes
SG::WriteHandle< sTGCSimHitCollection >	m_sTGCSimHitCollection
const sTgcHitIdHelper *	m_muonHelper {sTgcHitIdHelper::GetHelper()}
unsigned	m_baseDepth {0}
	basic depth to travel along the G4 history. More...
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 17 of file sTGCSensitiveDetector.h.

Constructor & Destructor Documentation

sTGCSensitiveDetector()

sTGCSensitiveDetector::sTGCSensitiveDetector	(	const std::string &	name,
		const std::string &	hitCollectionName,
		unsigned	baseDepth
	)

construction/destruction

Definition at line 18 of file sTGCSensitiveDetector.cxx.

                                                                : 
    G4VSensitiveDetector( name ),
    AthMessaging{name},
    m_sTGCSimHitCollection( hitCollectionName ),
    m_baseDepth{baseDepth} {}

~sTGCSensitiveDetector()

sTGCSensitiveDetector::~sTGCSensitiveDetector ( )

default

Member Function Documentation

FRIEND_TEST() [1/2]

sTGCSensitiveDetector::FRIEND_TEST ( sTGCSensitiveDetectortest  , Initialize    )

private

FRIEND_TEST() [2/2]

sTGCSensitiveDetector::FRIEND_TEST ( sTGCSensitiveDetectortest  , ProcessHits    )

private

Initialize()

void sTGCSensitiveDetector::Initialize ( G4HCofThisEvent *  HCE )

finaloverride

member functions

Definition at line 27 of file sTGCSensitiveDetector.cxx.

{
  if (!m_sTGCSimHitCollection.isValid()) m_sTGCSimHitCollection = std::make_unique<sTGCSimHitCollection>();
}

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;
}

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.

{ 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;
  }
}

ProcessHits()

G4bool sTGCSensitiveDetector::ProcessHits ( G4Step *  aStep, G4TouchableHistory *  ROhist  )

finaloverride

Definition at line 32 of file sTGCSensitiveDetector.cxx.

{
  G4Track* currentTrack = aStep->GetTrack();
  int charge=currentTrack->GetDefinition()->GetPDGCharge();
 
  bool geantinoHit = (currentTrack->GetDefinition()==G4Geantino::GeantinoDefinition()) ||
                     (currentTrack->GetDefinition()==G4ChargedGeantino::ChargedGeantinoDefinition());
 
  if (!charge && (!geantinoHit)) return false;
  //  G4cout << "\t\t sTGCSD: Hit in a sensitive layer!!!!! " << G4endl;
  G4StepPoint* postStep=aStep->GetPostStepPoint();
  const G4Step* post_Step=aStep->GetTrack()->GetStep();
 
  Amg::Vector3D position = Amg::Hep3VectorToEigen(postStep->GetPosition());
  
  G4StepPoint* preStep = aStep->GetPreStepPoint();
  Amg::Vector3D preposition = Amg::Hep3VectorToEigen(preStep->GetPosition());
 
  int pdgCode=currentTrack->GetDefinition()->GetPDGEncoding();
 
  float globalTime=postStep->GetGlobalTime();
 
  Amg::Vector3D direction= Amg::Hep3VectorToEigen( postStep->GetMomentumDirection() );
  float depositEnergy=post_Step->GetTotalEnergyDeposit();
 
  if (depositEnergy<0.0001 && (!geantinoHit)) return false;
 
  const G4TouchableHistory* touchHist = static_cast<const G4TouchableHistory*>(aStep->GetPreStepPoint()->GetTouchable());
 
  // int iDepth=touchHist->GetHistoryDepth();
  //  G4cout << "\t\t\t\t Touchable history dump "<<G4endl;
  int nLayer=touchHist->GetVolume(m_baseDepth)->GetCopyNo();
  std::string chName=touchHist->GetVolume(1+m_baseDepth)->GetLogicalVolume()->GetName();
  //G4cout << "sTGCSensitiveDetector name: "<<chName<<G4endl;
  std::string subType=chName.substr(chName.find('-')+1);
  //G4cout << "\t\t sType: "<<subType);
  if (subType[0]!='T'&&subType[0]!='Q' ) ATH_MSG_WARNING(" something is wrong, this is no sTGC! "<<chName<<", "<<Amg::toString(preposition));
  int iRing = std::atoi(&subType[2]) -1;
  ATH_MSG_VERBOSE("Volume name: "<<chName<<", nLayer: "<<nLayer<<", subType: "<<subType<<", iRing: "<<iRing);
 
  // identifiers have eta naming 0-2, eta encoded in subtype is 1-3
  // double phiDiff=2*M_PI;
 
  G4ThreeVector posH=postStep->GetPosition(); //posH is equivalent to position - eigen not used to avoid additional dependence on EventPrimitives
  if (subType[1]=='L') posH.rotateZ(M_PI/8.);
  double phiHit=posH.phi();
  if(phiHit<=0) phiHit+=2.*M_PI;
  int iPhi=1+int(phiHit/(M_PI/4.));
  iPhi*=2;
  if (subType[1]=='L') iPhi-=1;
 
  int iSide=1;
  if (position.z()<0) iSide=-1;
 
  int mLayer=0;
  if (subType[1]=='S')
  {
    if (subType[3]=='C') mLayer=1;
    else if (subType[3]=='P') mLayer=2;
  }
  else if (subType[1]=='L')
  {
    if (subType[3]=='P') mLayer=1;
    else if (subType[3]=='C') mLayer=2;
  }
 
  if (mLayer != 1 && mLayer !=2) G4cout << " something is wrong - multilayer index is " << mLayer << G4endl;
 
  int sTgcId = m_muonHelper->BuildsTgcHitId(subType, iPhi, iRing, mLayer,nLayer, iSide);
  TrackHelper trHelp(aStep->GetTrack());
  m_sTGCSimHitCollection->Emplace(sTgcId,globalTime,position,pdgCode,direction,depositEnergy,
                                  trHelp.GenerateParticleLink(),
                                  preStep->GetKineticEnergy(),preposition);
 
  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;
}

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_baseDepth

unsigned sTGCSensitiveDetector::m_baseDepth {0}

private

basic depth to travel along the G4 history.

For jobs run with the legacy geometry database, it's zero. Otherwise, in the new sqlite workflow it's 1

Definition at line 39 of file sTGCSensitiveDetector.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_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_muonHelper

const sTgcHitIdHelper* sTGCSensitiveDetector::m_muonHelper {sTgcHitIdHelper::GetHelper()}

private

Definition at line 36 of file sTGCSensitiveDetector.h.

m_nm

std::string AthMessaging::m_nm

privateinherited

Message source name.

Definition at line 129 of file AthMessaging.h.

m_sTGCSimHitCollection

SG::WriteHandle<sTGCSimHitCollection> sTGCSensitiveDetector::m_sTGCSimHitCollection

private

Definition at line 35 of file sTGCSensitiveDetector.h.

---

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

• sTGCSensitiveDetector.h
• sTGCSensitiveDetector.cxx