G4UA::G4AtlantisDumper Class Reference

#include <G4AtlantisDumper.h>

Inheritance diagram for G4UA::G4AtlantisDumper:

Collaboration diagram for G4UA::G4AtlantisDumper:

Classes
struct	Config

Public Member Functions
	G4AtlantisDumper (const Config &config)
virtual void	UserSteppingAction (const G4Step *) override
virtual void	EndOfEventAction (const G4Event *) override
virtual void	BeginOfEventAction (const G4Event *) 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
void	initMessaging () const
	Initialize our message level and MessageSvc. More...

Private Attributes
Config	m_config
ServiceHandle< StoreGateSvc >	m_evtStore
	Pointer to StoreGate (event store by default) More...
ServiceHandle< StoreGateSvc >	m_detStore
	Pointer to StoreGate (detector store by default) More...
int	m_nsteps
std::ofstream *	m_file
std::string	m_filename
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 21 of file G4AtlantisDumper.h.

Constructor & Destructor Documentation

G4AtlantisDumper()

G4UA::G4AtlantisDumper::G4AtlantisDumper

(

const Config &

config

)

Definition at line 24 of file G4AtlantisDumper.cxx.

                                                        :
    AthMessaging(Gaudi::svcLocator()->service< IMessageSvc >( "MessageSvc" ),"G4AtlantisDumper"),
    m_config(config),
    m_evtStore("StoreGateSvc/StoreGateSvc","G4AtlantisDumper"),
    m_detStore("StoreGateSvc/DetectorStore","G4AtlantisDumper"),
    m_nsteps(0),
    m_file(NULL)
  {}

Member Function Documentation

BeginOfEventAction()

void G4UA::G4AtlantisDumper::BeginOfEventAction ( const G4Event *  )

overridevirtual

Definition at line 109 of file G4AtlantisDumper.cxx.

  {
    m_nsteps = 0;
 
    SG::ReadHandle<xAOD::EventInfo> evt("EventInfo");
    if(!evt.isValid()) {
      ATH_MSG_INFO("G4AtlantisDumper could not get event info!");
      return;
    }
 
    int athevent = evt->eventNumber();
    int athrun = evt->runNumber();
    ATH_MSG_INFO("G4AtlantisDumper: Athena run event is "<<athrun<<" "<<athevent);
    char buf[1000];
    int length=snprintf(buf, 1000, "G4Atlantis_%d_%d.txt", athrun, athevent);
    if(999<length) {
      ATH_MSG_WARNING( "BeginOfEvent:: m_filename string (" << buf << ") was truncated" );
    }
    ATH_MSG_INFO( "G4AtlantisDumper: Opening m_file "<<buf);
    m_file = new std::ofstream();
    m_file->open(buf);
    m_filename = buf;
  }

EndOfEventAction()

void G4UA::G4AtlantisDumper::EndOfEventAction ( const G4Event *  event )

overridevirtual

Definition at line 99 of file G4AtlantisDumper.cxx.

  {
    ATH_MSG_INFO( "Goodbye from G4AtlantisDumper, event "<<event->GetEventID());
    m_file->close();
    delete m_file;
 
    ATH_MSG_INFO("G4AtlantisDumper: zipping "<<m_filename);
    system(("zip -m "+m_filename+".zip "+m_filename).c_str());
  }

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

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

UserSteppingAction()

void G4UA::G4AtlantisDumper::UserSteppingAction ( const G4Step *  aStep )

overridevirtual

Definition at line 33 of file G4AtlantisDumper.cxx.

  {
    const G4Track *t=aStep->GetTrack();
    //if(m_nsteps>40) return;
    if (!t) return;
 
    int id=t->GetTrackID();
    int pid=t->GetParentID();
 
    int fs=aStep->IsFirstStepInVolume();
    int ls=aStep->IsLastStepInVolume();
    int stepn=t->GetCurrentStepNumber();
 
    const G4TrackVector *sec=aStep->GetSecondary();
    int nsec=0; if (sec) nsec=sec->size();
 
    double tedep=aStep->GetTotalEnergyDeposit();
    double tnedep=aStep->GetNonIonizingEnergyDeposit();
    double length=aStep->GetStepLength();
    int pdg=t->GetDefinition()->GetPDGEncoding();
 
    bool keep=false;
    if (fs||ls) keep=true;
    if (tedep>m_config.tedep_cut) keep=true;//default = 0.1 MeV
    if (length>m_config.length_cut) keep=true;//default = 1.0 mm
    if (aStep->GetPreStepPoint()->GetMomentum().mag()<m_config.p_cut) keep=false;//default = 10 MeV
    if (aStep->GetPreStepPoint()->GetKineticEnergy()<m_config.ke_cut) keep=false;//default = 10 MeV
    if (!keep) return;
 
    //const G4VPhysicalVolume *v =t->GetVolume();
    //const G4Material *m =t->GetMaterial();
 
    const G4VProcess *p =t->GetCreatorProcess();
    G4String pn="none"; if(p) pn=p->GetProcessName();
 
    G4ThreeVector pos1=aStep->GetPreStepPoint()->GetPosition();
    double x1=pos1.x();    double y1=pos1.y();    double z1=pos1.z();
    double t1=aStep->GetPreStepPoint()->GetGlobalTime();
    G4ThreeVector pos2=aStep->GetPostStepPoint()->GetPosition();
    double x2=pos2.x();    double y2=pos2.y();    double z2=pos2.z();
    double t2=aStep->GetPostStepPoint()->GetGlobalTime();
    G4ThreeVector p1=aStep->GetPreStepPoint()->GetMomentum();
    double px1=p1.x();    double py1=p1.y();    double pz1=p1.z();
    G4ThreeVector p2=aStep->GetPostStepPoint()->GetMomentum();
    double px2=p2.x();    double py2=p2.y();    double pz2=p2.z();
    double ke1=aStep->GetPreStepPoint()->GetKineticEnergy();
    double ke2=aStep->GetPostStepPoint()->GetKineticEnergy();
    //GetVelocity() //GetBeta() //GetGamma()
    //GetSensitiveDetector() //GetPolarization()
    double mass=aStep->GetPreStepPoint()->GetMass();
    double charge=aStep->GetPreStepPoint()->GetCharge();
    //GetMagneticMoment()
 
    (*m_file)<<"ACH_G4Step";
    (*m_file)<<" "<<m_nsteps;    (*m_file)<<" "<<id;    (*m_file)<<" "<<pid;    (*m_file)<<" "<<(int)fs;    (*m_file)<<" "<<(int)ls;    (*m_file)<<" "<<stepn;
    (*m_file)<<" "<<pn;    (*m_file)<<" "<<nsec;
    (*m_file)<<" "<<tedep;    (*m_file)<<" "<<tnedep;    (*m_file)<<" "<<length;    (*m_file)<<" "<<pdg;    (*m_file)<<" "<<mass;    (*m_file)<<" "<<charge;
    (*m_file)<<" "<<x1;    (*m_file)<<" "<<y1;    (*m_file)<<" "<<z1;    (*m_file)<<" "<<t1;
    (*m_file)<<" "<<px1;    (*m_file)<<" "<<py1;    (*m_file)<<" "<<pz1;    (*m_file)<<" "<<ke1;
    (*m_file)<<" "<<x2;    (*m_file)<<" "<<y2;    (*m_file)<<" "<<z2;    (*m_file)<<" "<<t2;
    (*m_file)<<" "<<px2;    (*m_file)<<" "<<py2;    (*m_file)<<" "<<pz2;    (*m_file)<<" "<<ke2;
    (*m_file)<<std::endl;
 
    ++m_nsteps;
  }

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_config

Config G4UA::G4AtlantisDumper::m_config

private

Definition at line 42 of file G4AtlantisDumper.h.

m_detStore

ServiceHandle<StoreGateSvc> G4UA::G4AtlantisDumper::m_detStore

private

Pointer to StoreGate (detector store by default)

Definition at line 47 of file G4AtlantisDumper.h.

m_evtStore

ServiceHandle<StoreGateSvc> G4UA::G4AtlantisDumper::m_evtStore

private

Pointer to StoreGate (event store by default)

Definition at line 45 of file G4AtlantisDumper.h.

m_file

std::ofstream* G4UA::G4AtlantisDumper::m_file

private

Definition at line 49 of file G4AtlantisDumper.h.

m_filename

std::string G4UA::G4AtlantisDumper::m_filename

private

Definition at line 50 of file G4AtlantisDumper.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_nm

std::string AthMessaging::m_nm

privateinherited

Message source name.

Definition at line 129 of file AthMessaging.h.

m_nsteps

int G4UA::G4AtlantisDumper::m_nsteps

private

Definition at line 48 of file G4AtlantisDumper.h.

---

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

• G4AtlantisDumper.h
• G4AtlantisDumper.cxx