ATLAS Offline Software
Public Member Functions | Private Member Functions | Private Attributes | List of all members
G4UA::StoppedParticleAction Class Reference

NEEDS DOCUMENTATION. More...

#include <StoppedParticleAction.h>

Inheritance diagram for G4UA::StoppedParticleAction:
Collaboration diagram for G4UA::StoppedParticleAction:

Public Member Functions

 StoppedParticleAction ()
 
virtual void UserSteppingAction (const G4Step *) override
 
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

bool isSUSYParticle (const int) const
 
void initMessaging () const
 Initialize our message level and MessageSvc. More...
 

Private Attributes

TrackFastSimSD * m_fsSD
 
bool m_init
 
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

NEEDS DOCUMENTATION.

Definition at line 19 of file StoppedParticleAction.h.

Constructor & Destructor Documentation

◆ StoppedParticleAction()

G4UA::StoppedParticleAction::StoppedParticleAction ( )

Definition at line 46 of file StoppedParticleAction.cxx.

47  : AthMessaging(Gaudi::svcLocator()->service< IMessageSvc >( "MessageSvc"),
48  "StoppedParticleAction"),
49  m_fsSD(0), m_init(false)
50  {}

Member Function Documentation

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

40 {
41  m_imsg = Athena::getMessageSvc();
42  m_lvl = m_imsg ?
43  static_cast<MSG::Level>( m_imsg.load()->outputLevel(m_nm) ) :
44  MSG::INFO;
45 }

◆ isSUSYParticle()

bool G4UA::StoppedParticleAction::isSUSYParticle ( const int  id) const
private

Definition at line 31 of file StoppedParticleAction.cxx.

32  {
33  if (id==1000021 || id==1000005 || id==1000006 || id==1000512 || id==1000522 || id==1000991 || id==1000993 ||
34  id==1000612 || id==1000622 || id==1000632 || id==1000642 || id==1000652 || id==1005211 ||
35  id==1006113 || id==1006211 || id==1006213 || id==1006223 || id==1006311 ||
36  id==1006313 || id==1006321 || id==1006323 || id==1006333 ||
37  id==1009111 || id==1009113 || id==1009211 || id==1009213 || id==1009311 ||
38  id==1009313 || id==1009321 || id==1009323 || id==1009223 || id==1009333 ||
39  id==1092112 || id==1091114 || id==1092114 || id==1092212 || id==1092214 || id==1092224 ||
40  id==1093114 || id==1093122 || id==1093214 || id==1093224 || id==1093314 || id==1093324 || id==1093334)
41  return true;
42  return false;
43  }

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

165 {
166  MsgStream* ms = m_msg_tls.get();
167  if (!ms) {
168  if (!m_initialized.test_and_set()) initMessaging();
169  ms = new MsgStream(m_imsg,m_nm);
170  m_msg_tls.reset( ms );
171  }
172 
173  ms->setLevel (m_lvl);
174  return *ms;
175 }

◆ 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
lvlThe message level to test against
Returns
boolean Indicating if messages at given level will be printed
Return values
trueMessages at level "lvl" will be printed

Definition at line 151 of file AthMessaging.h.

152 {
153  if (!m_initialized.test_and_set()) initMessaging();
154  if (m_lvl <= lvl) {
155  msg() << lvl;
156  return true;
157  } else {
158  return false;
159  }
160 }

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

29 {
30  m_lvl = lvl;
31 }

◆ UserSteppingAction()

void G4UA::StoppedParticleAction::UserSteppingAction ( const G4Step *  aStep)
overridevirtual

Definition at line 53 of file StoppedParticleAction.cxx.

54  {
55 
56  // Trigger if the energy is below our threshold or if the R-hadron is decaying
57  int id = std::abs(aStep->GetTrack()->GetDynamicParticle()->GetDefinition()->GetPDGEncoding());
58 
59  // Special treatment for SUSY particles and R-hadrons
60  if (id>=1000000 && id<=1100000 &&
61  isSUSYParticle(id)){
62 
63  G4Material * mat = aStep->GetTrack()->GetMaterial();
64  double minA=1500000.;
65  for (unsigned int i=0;i<mat->GetNumberOfElements();++i){
66  if (mat->GetElement(i) &&
67  minA>mat->GetElement(i)->GetN()){
68  minA=mat->GetElement(i)->GetN();
69  }
70  }
71 
72  // Stopping condition
73  if (aStep->GetPostStepPoint()->GetVelocity()>0.15*std::pow(minA,-2./3.)*CLHEP::c_light && // Stopping condition or...
74  ( !aStep->GetPostStepPoint()->GetProcessDefinedStep() || // null pointer?
75  aStep->GetPostStepPoint()->GetProcessDefinedStep()->GetProcessType()!=fDecay) ) // Decaying particle (does not fire for hadronic interactions)
76  return;
77 
78  if (!m_init){
79  m_init = true;
80 
81  G4SDManager * g4sdm = G4SDManager::GetSDMpointer();
82  if (!g4sdm) {
83  ATH_MSG_ERROR( "StoppedParticleFastSim could not get sensitive detector catalog." );
84  } else {
85  G4VSensitiveDetector * g4sd = g4sdm->FindSensitiveDetector("ToolSvc.SensitiveDetectorMasterTool.TrackFastSimSD");
86  if (!g4sd) {
87  ATH_MSG_ERROR( "StoppedParticleFastSim could not get ToolSvc.SensitiveDetectorMasterTool.TrackFastSimSD sensitive detector." );
88  } else {
89  m_fsSD = dynamic_cast<TrackFastSimSD*>(g4sd);
90  if (!m_fsSD) {
91  ATH_MSG_ERROR( "StoppedParticleFastSim could not cast the SD." );
92  }
93  } // found the SD
94  } // got the catalog
95  }
96 
97  if (m_fsSD) {
98  m_fsSD->WriteTrack( aStep->GetTrack() , false , true );
99  }
100  }
101 
102  aStep->GetTrack()->SetTrackStatus(fStopAndKill);
103  const G4TrackVector *tv = aStep->GetSecondary();
104  for (unsigned int i=0;i<tv->size();i++){
105  (*tv)[i]->SetTrackStatus(fStopAndKill);
106  }
107  }

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_fsSD

TrackFastSimSD* G4UA::StoppedParticleAction::m_fsSD
private

Definition at line 26 of file StoppedParticleAction.h.

◆ m_imsg

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

MessageSvc pointer.

Definition at line 135 of file AthMessaging.h.

◆ m_init

bool G4UA::StoppedParticleAction::m_init
private

Definition at line 27 of file StoppedParticleAction.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_nm

std::string AthMessaging::m_nm
privateinherited

Message source name.

Definition at line 129 of file AthMessaging.h.

The documentation for this class was generated from the following files:
AthMessaging::m_lvl
std::atomic< MSG::Level > m_lvl
Current logging level.
Definition: AthMessaging.h:138
TrackFastSimSD
Definition: TrackFastSimSD.h:24
mat
GeoMaterial * mat
Definition: LArDetectorConstructionTBEC.cxx:53
G4UA::StoppedParticleAction::isSUSYParticle
bool isSUSYParticle(const int) const
Definition: StoppedParticleAction.cxx:31
conifer::pow
constexpr int pow(int x)
Definition: conifer.h:20
G4UA::StoppedParticleAction::m_init
bool m_init
Definition: StoppedParticleAction.h:27
AthMessaging::m_imsg
std::atomic< IMessageSvc * > m_imsg
MessageSvc pointer.
Definition: AthMessaging.h:135
python.SystemOfUnits.ms
int ms
Definition: SystemOfUnits.py:132
Athena::getMessageSvc
IMessageSvc * getMessageSvc(bool quiet=false)
Definition: getMessageSvc.cxx:20
AthMessaging::AthMessaging
AthMessaging()
Default constructor:
TrigConf::MSGTC::Level
Level
Definition: Trigger/TrigConfiguration/TrigConfBase/TrigConfBase/MsgStream.h:21
ATH_MSG_ERROR
#define ATH_MSG_ERROR(x)
Definition: AthMsgStreamMacros.h:33
lumiFormat.i
int i
Definition: lumiFormat.py:92
AthMessaging::msg
MsgStream & msg() const
The standard message stream.
Definition: AthMessaging.h:164
G4UA::StoppedParticleAction::m_fsSD
TrackFastSimSD * m_fsSD
Definition: StoppedParticleAction.h:26
python.PhysicalConstants.c_light
float c_light
Definition: PhysicalConstants.py:63
AthMessaging::m_nm
std::string m_nm
Message source name.
Definition: AthMessaging.h:129
TrackFastSimSD::WriteTrack
void WriteTrack(const G4Track *, const bool, const bool)
Definition: TrackFastSimSD.cxx:97
AthMessaging::initMessaging
void initMessaging() const
Initialize our message level and MessageSvc.
Definition: AthMessaging.cxx:39
AthMessaging::m_msg_tls
boost::thread_specific_ptr< MsgStream > m_msg_tls
MsgStream instance (a std::cout like with print-out levels)
Definition: AthMessaging.h:132
Generated on Tue Jul 2 2024 22:16:40 for ATLAS Offline Software by   doxygen 1.8.18