G4UA::iGeant4::PhysicsValidationUserAction Class Reference

#include <PhysicsValidationUserAction.h>

Inheritance diagram for G4UA::iGeant4::PhysicsValidationUserAction:

Collaboration diagram for G4UA::iGeant4::PhysicsValidationUserAction:

Classes
struct	Config

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

Private Attributes
Config	m_config
ISF::IGeoIDSvc *	m_geoIDSvcQuick
	access to the central ISF GeoID serice
TTree *	m_particles
	ROOT tree containing track info.
int	m_pdg
int	m_scIn
int	m_scEnd
int	m_gen
int	m_geoID
float	m_pth
float	m_pph
float	m_p
float	m_eloss
float	m_radloss
float	m_ionloss
float	m_wzOaTr
float	m_thIn
float	m_phIn
float	m_dIn
float	m_thEnd
float	m_phEnd
float	m_dEnd
float	m_X0
float	m_L0
float	m_wZ
float	m_dt
TTree *	m_interactions
	ROOT tree containing vertex info.
int	m_process
int	m_pdg_mother
int	m_gen_mother
int	m_nChild
float	m_vtx_dist
float	m_vtx_theta
float	m_vtx_phi
float	m_vtx_e_diff
float	m_vtx_p_diff
float	m_vtx_plong_diff
float	m_vtx_pperp_diff
float	m_p_mother
float	m_radLength
int	m_pdg_child [MAXCHILDREN] ={}
float	m_fp_child [MAXCHILDREN] ={}
float	m_oa_child [MAXCHILDREN] ={}
int	m_volumeOffset
int	m_minHistoryDepth
int	m_currentTrack
std::map< int, int >	m_trackGenMap
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 33 of file PhysicsValidationUserAction.h.

Constructor & Destructor Documentation

PhysicsValidationUserAction()

G4UA::iGeant4::PhysicsValidationUserAction::PhysicsValidationUserAction

(

const Config &

config

)

Definition at line 43 of file PhysicsValidationUserAction.cxx.

      : AthMessaging("PhysicsValidationUserAction")
      , m_config(config)
      , m_geoIDSvcQuick(nullptr)
      // branches
      , m_particles(nullptr)
      , m_pdg(0)
      , m_scIn(0)
      , m_scEnd(0)
      , m_gen(0)
      , m_geoID(0)
      , m_pth(0.)
      , m_pph(0.)
      , m_p(0.)
      , m_eloss(0.)
      , m_radloss(0.)
      , m_ionloss(0.)
      , m_wzOaTr(0), m_thIn(0), m_phIn(0), m_dIn(0)
      , m_thEnd(0), m_phEnd(0), m_dEnd(0)
      , m_X0(0), m_L0(0), m_wZ(0), m_dt(0)
      // more branches
      , m_interactions(nullptr)
      , m_process(0), m_pdg_mother(0), m_gen_mother(0), m_nChild(0)
      , m_vtx_dist(0), m_vtx_theta(0), m_vtx_phi(0), m_vtx_e_diff(0)
      , m_vtx_p_diff(0), m_vtx_plong_diff(0), m_vtx_pperp_diff(0)
      , m_p_mother(0), m_radLength(0)
      , m_volumeOffset(1)
      , m_minHistoryDepth(0)
      , m_currentTrack(0)
    {
      setLevel(m_config.verboseLevel);
    }

Member Function Documentation

BeginOfEventAction()

void G4UA::iGeant4::PhysicsValidationUserAction::BeginOfEventAction ( const G4Event * )

finaloverridevirtual

Definition at line 76 of file PhysicsValidationUserAction.cxx.

    {
 
 
      ATH_MSG_VERBOSE("idRmax: "<<m_config.idR      <<", idZmax: "<<m_config.idZ);
      ATH_MSG_VERBOSE("caloR : "<<m_config.caloRmean<<", caloZ : "<<m_config.caloZmean);
      ATH_MSG_VERBOSE("muonR : "<<m_config.muonRmean<<", muonZ : "<<m_config.muonZmean);
      ATH_MSG_VERBOSE("muonR : "<<m_config.cavernRmean<<", muonZ : "<<m_config.cavernZmean);
      ATH_MSG_VERBOSE("cavernR : "<<m_config.cavernRmean<<", cavernZ : "<<m_config.cavernZmean);
      
      m_currentTrack = -1;
      m_trackGenMap.clear();
      //m_idToStackParticleMap.clear();
      
      
    }

BeginOfRunAction()

void G4UA::iGeant4::PhysicsValidationUserAction::BeginOfRunAction ( const G4Run * )

finaloverridevirtual

Definition at line 106 of file PhysicsValidationUserAction.cxx.

    {
 
      if (m_config.geoIDSvc.retrieve().isFailure()) {
    ATH_MSG_FATAL("Could not retrieve ISF GeoID Svc: " << m_config.geoIDSvc);
    return;
      }
 
      m_geoIDSvcQuick = &(*m_config.geoIDSvc);
 
      // setup for validation mode
      if ( m_config.validationOutput) {
    
    // retrieve the histogram service
    if ( m_config.thistSvc.retrieve().isSuccess() ) {
      // Create the prefix of histogram names for the THistSvc
      const char *treeName="particles";
      const std::string prefix = "/" + m_config.validationStream + "/"+ treeName;
      m_particles = new TTree( treeName, treeName );
      m_particles->Branch("pdg"       , &m_pdg    , "pdg/I"         );         // pdg id
      m_particles->Branch("scIn"      , &m_scIn   , "scIn/I"         );        // input process
      m_particles->Branch("scOut"     , &m_scEnd  , "scOut/I"         );       // endpoint process
      m_particles->Branch("gen"       , &m_gen    , "gen/I"         );         // generation (0 for primary)
      m_particles->Branch("geoID"     , &m_geoID  , "geoID/I"         );       // subdetector id
      m_particles->Branch("pth"       , &m_pth     , "pth/F"            );     // input momentum polar angle
      m_particles->Branch("pph"       , &m_pph     , "pph/F"            );     // input momemtum azimuthal angle
      m_particles->Branch("p"         , &m_p      , "p/F"            );        // input momentum
      m_particles->Branch("eloss"     , &m_eloss  , "eloss/F"            );    // energy loss
      m_particles->Branch("radloss"     , &m_radloss  , "radloss/F"            );   // radiative eloss
      m_particles->Branch("ionloss"     , &m_ionloss  , "ionloss/F"            );   // ionization eloss
      m_particles->Branch("wzOaTr"     , &m_wzOaTr  , "wzOaTr/F"            );      // zOverZtimesRho times dInX0
      m_particles->Branch("X0"        , &m_X0     , "X0/F"            );            // dInX0 (material thickness)
      m_particles->Branch("L0"        , &m_L0     , "L0/F"            );            // dInL0
      m_particles->Branch("wZ"        , &m_wZ     , "wZ/F"            );            // averageZ time dInX0
      m_particles->Branch("dt"        , &m_dt     , "dt/F"            );            // time interval
      m_particles->Branch("thIn"      , &m_thIn   , "thIn/F"            );          // polar angle input position
      m_particles->Branch("phIn"      , &m_phIn   , "phIn/F"            );          // azimuthal angle input position
      m_particles->Branch("dIn"       , &m_dIn    , "dIn/F"            );           // distance input position
      m_particles->Branch("thEnd"     , &m_thEnd  , "thEnd/F"            );         // polar angle exit position
      m_particles->Branch("phEnd"     , &m_phEnd  , "phEnd/F"            );         // azimuthal angle exit position
      m_particles->Branch("dEnd"      , &m_dEnd   , "dEnd/F"            );          // distance exit position
      
      // register the Tree to the THistSvc and return it's StatusCode
      if(m_config.thistSvc->regTree( prefix, m_particles)!=StatusCode::SUCCESS){
        ATH_MSG_ERROR("Cannot register validation tree");
      }
      
      m_X0=0.;
      m_L0=0.;
      m_wZ=0.;
      m_wzOaTr=0.;
      
      const char *treeNameInt="interactions";
      const std::string prefixInt = "/" + m_config.validationStream + "/"+ treeNameInt;
      m_interactions = new TTree( treeNameInt, treeNameInt );
      m_interactions->Branch("process"     , &m_process  , "process/I"         );
      m_interactions->Branch("pdg_mother"  , &m_pdg_mother  , "pdg_mother/I"         );
      m_interactions->Branch("gen_mother"  , &m_gen_mother  , "gen_mother/I"         );
      m_interactions->Branch("nChild"      , &m_nChild      , "nch/I"         );
      m_interactions->Branch("vtx_dist"    , &m_vtx_dist    , "vtx_dist/F"         );
      m_interactions->Branch("vtx_theta"   , &m_vtx_theta   , "vtx_theta/F"         );
      m_interactions->Branch("vtx_phi"     , &m_vtx_phi     , "vtx_phi/F"         );
      m_interactions->Branch("vtx_e_diff"  , &m_vtx_e_diff  , "vtx_e_diff/F"         );
      m_interactions->Branch("vtx_p_diff"  , &m_vtx_p_diff  , "vtx_p_diff/F"         );
      m_interactions->Branch("vtx_plong_diff"  , &m_vtx_plong_diff  , "vtx_plong_diff/F"         );
      m_interactions->Branch("vtx_pperp_diff"  , &m_vtx_pperp_diff  , "vtx_pperp_diff/F"         );
      m_interactions->Branch("radLength"    , &m_radLength     , "radLength/F"         );
      m_interactions->Branch("p_mother"    , &m_p_mother    , "p_mother/F"         );
      m_interactions->Branch("pdg_child"   , m_pdg_child    , "pdg_child[nch]/I"         );
      m_interactions->Branch("fp_child"    , m_fp_child     , "fp_child[nch]/F"         );
      m_interactions->Branch("oa_child"    , m_oa_child     , "oa_child[nch]/F"         );
      // register the Tree to the THistSvc and return it's StatusCode
      if(m_config.thistSvc->regTree( prefixInt, m_interactions)!=StatusCode::SUCCESS){
        ATH_MSG_ERROR("Cannot register validation tree");
      }
      
      m_radLength = 0.;
      
    }
      }
      
      
      // get the geometry manager and check how many layers are present.
      G4TransportationManager *transportationManager(G4TransportationManager::GetTransportationManager());
      G4LogicalVolume *world((*(transportationManager->GetWorldsIterator()))->GetLogicalVolume());
      ATH_MSG_VERBOSE("World G4LogicalVolume Name: " << world->GetName() << " has " << world->GetNoDaughters() << " daughters.");
      if ("World::World"==world->GetName())
    {
      ATH_MSG_INFO("Atlas::Atlas is not the world volume, so assume we are in a cosmics job.");
      //Cosmics-specific configuration.
      m_volumeOffset=2;
      m_minHistoryDepth=1;
    }
      
    }

EndOfEventAction()

void G4UA::iGeant4::PhysicsValidationUserAction::EndOfEventAction ( const G4Event * )

finaloverridevirtual

Definition at line 93 of file PhysicsValidationUserAction.cxx.

    {
      
      m_X0=0.;
      m_L0=0.;
      m_wZ=0.;
      m_wzOaTr=0.;
      
      m_radLength=0;
      
      return;
    }

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

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 163 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 178 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_lvl <= lvl) {
    msg() << lvl;
    return true;
  } else {
    return false;
  }
}

PreUserTrackingAction()

void G4UA::iGeant4::PhysicsValidationUserAction::PreUserTrackingAction ( const G4Track * )

finaloverridevirtual

Definition at line 504 of file PhysicsValidationUserAction.cxx.

    {
      return;
    }

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::iGeant4::PhysicsValidationUserAction::UserSteppingAction ( const G4Step * aStep )

finaloverridevirtual

Definition at line 202 of file PhysicsValidationUserAction.cxx.

    {
      //std::cout<<"PhysicsValidationUserAction::SteppingAction"<<std::endl;
      
      // process identification
      const G4VProcess * process = aStep->GetPostStepPoint()->GetProcessDefinedStep();
      
      // material info
      G4StepPoint *preStep=aStep->GetPreStepPoint();
      G4StepPoint *postStep=aStep->GetPostStepPoint();
      
      G4ThreeVector mom = preStep->GetMomentum();
      const G4ThreeVector& pos = preStep->GetPosition();
      
      // info about generating particle
      G4Track * track = aStep->GetTrack();
      
      int trackID=track->GetTrackID();
      const bool trackIsAlive = (track->GetTrackStatus() == fAlive || track->GetTrackStatus() == fStopButAlive);
      
      if (trackID != m_currentTrack) {     // for new G4Track only
    
    m_pdg = track->GetDefinition()->GetPDGEncoding();
    const G4VProcess* creation = track->GetCreatorProcess();
    m_scIn = creation? creation->GetProcessSubType() : -1;
    
    VTrackInformation * trackInfo= static_cast<VTrackInformation*>(track->GetUserInformation());
    HepMC::GenParticlePtr currentGenParticle= trackInfo ? trackInfo->GetCurrentGenParticle() : nullptr;
    HepMC::GenVertexPtr vtx = currentGenParticle ? currentGenParticle->production_vertex() : nullptr;
    m_gen = currentGenParticle? 0 : -1;
    
    if (currentGenParticle)  { // mc truth known
      while (currentGenParticle && vtx ) {
        int pdgID=currentGenParticle->pdg_id();
#ifdef HEPMC3
        const HepMC::GenParticlePtr  genmom = vtx->particles_in().size()>0 ? vtx->particles_in().front() : nullptr;
        if ( genmom && pdgID!=genmom->pdg_id() ) m_gen++;
        else if (vtx->particles_out().size()>0 && currentGenParticle!=vtx->particles_out().front()) m_gen++;
 
#else
        HepMC::GenParticlePtr genmom = vtx->particles_in_size()>0 ? *(vtx->particles_in_const_begin()) : nullptr;
        if ( genmom && pdgID!=genmom->pdg_id() ) m_gen++;
        else if (vtx->particles_out_size()>0 && currentGenParticle!=*(vtx->particles_out_const_begin())) m_gen++;
#endif
        vtx = genmom ? genmom->production_vertex() : nullptr;
        currentGenParticle = genmom;
      }
    } else {
      // retrieve info from parent track
      int parentID=track->GetParentID();
      std::map<int, int>::iterator genIt = m_trackGenMap.find(parentID);
      if ( genIt != m_trackGenMap.end()) m_gen = (genIt->second >= 0) ? genIt->second+1 : genIt->second-1;
    }
    
    m_trackGenMap.try_emplace(trackID, m_gen);
    
    m_currentTrack=trackID;
    
    m_radLength = 0.;
      }
      
      AtlasDetDescr::AtlasRegion geoID = m_geoIDSvcQuick->identifyNextGeoID( pos.x(),
                                         pos.y(),
                                         pos.z(),
                                         mom.x(),
                                         mom.y(),
                                         mom.z() );
      
      
      
      double stepLength = aStep->GetStepLength();
      double radLengthInX0  = preStep->GetMaterial()->GetRadlen();
      double l0  = preStep->GetMaterial()->GetNuclearInterLength();
      float stepInX0   = stepLength/radLengthInX0;
      
      if (stepInX0>1.e-06) {
    
    m_X0 += stepInX0;
    m_radLength += stepInX0;
    
    if (l0>0.) {
      m_L0 += stepLength/l0;
      // average Z/A over fraction of atoms rather than weight fraction
      // const G4double* fVec = preStep->GetMaterial()->GetFractionVector();      // mass fraction
      double totNA =  preStep->GetMaterial()->GetTotNbOfAtomsPerVolume();
      const G4ElementVector* eVec = preStep->GetMaterial()->GetElementVector();
      const G4double* atVector = preStep->GetMaterial() ->GetVecNbOfAtomsPerVolume();
      float mFactor =stepInX0* preStep->GetMaterial()->GetDensity();
      
      float zOverA = 0.;    float frSum = 0.;
      for (unsigned int i=0; i<eVec->size(); i++) {
        float fEl = atVector ?  atVector[i]/totNA : 0.;
        m_wZ += stepInX0*fEl*((*eVec)[i]->GetZ());
        //std::cout <<"elements:"<<i<<","<<fVec[i]<<":"<<(*eVec)[i]->GetZ()<< ","<<m_wZ<<","<<m_wZ/m_X0<<std::endl;
        //m_wA += stepInX0*fVec[i]*((*eVec)[i]->GetA());
        zOverA += fEl*((*eVec)[i]->GetZ())/((*eVec)[i]->GetA());
        frSum += fEl;
      }
      if (fabs(frSum-1.)>0.01)  ATH_MSG_DEBUG("G4 material description inconsistent, sum of element fractions:"<< frSum);
      m_wzOaTr += mFactor*zOverA;
    }
    
      }
      
      
      // save interaction info (if any)
      if ( process && process->GetProcessSubType()>0 && process->GetProcessSubType()!=91) {
    
    float eloss = postStep->GetMomentum().mag()-preStep->GetMomentum().mag();
    
    if (process->GetProcessSubType()==2 ) m_ionloss+=eloss;
    if (process->GetProcessSubType()==3 ) m_radloss+=eloss;
    
    AtlasG4EventUserInfo* atlasG4EvtUserInfo = static_cast<AtlasG4EventUserInfo*> (G4EventManager::GetEventManager()->GetConstCurrentEvent()->GetUserInformation());
    VTrackInformation * trackInfo = static_cast<VTrackInformation*>(track->GetUserInformation());
    ::iGeant4::Geant4TruthIncident truth( aStep, *trackInfo->GetBaseISFParticle(), geoID, atlasG4EvtUserInfo);
    unsigned int nSec = truth.numberOfChildren();
    if (nSec>0 || !trackIsAlive ) {      // save interaction info
      //std::cout <<"interaction:"<< process->GetProcessSubType() <<":"<<nSec<< std::endl;
      m_process=process->GetProcessSubType();
      m_pdg_mother = track->GetDefinition()->GetPDGEncoding();
      m_gen_mother = m_gen;
      G4ThreeVector mom = preStep->GetMomentum();
      m_p_mother   = mom.mag();
      
      m_vtx_dist    = postStep->GetPosition().mag();
      m_vtx_theta   = postStep->GetPosition().theta();
      m_vtx_phi     = postStep->GetPosition().phi();
      
      int iPrimSurv = !trackIsAlive ? 0 : 1;
      m_nChild     = nSec+iPrimSurv;
      
      G4ThreeVector pbal(mom);
      
      if (iPrimSurv>0) {
        m_pdg_child[0] = m_pdg_mother;
        m_fp_child[0] = postStep->GetMomentum().mag()/m_p_mother;
        m_oa_child[0] = mom*postStep->GetMomentum()/m_p_mother/postStep->GetMomentum().mag();
        pbal -= postStep->GetMomentum();
      }
      
      unsigned int nSecMax = nSec+iPrimSurv> MAXCHILDREN ? MAXCHILDREN-iPrimSurv : nSec;
      for (unsigned int isec=0; isec< nSec; isec++) {
        G4ThreeVector secMom = truth.childP(isec);
        if (isec<nSecMax) {
          m_pdg_child[isec+iPrimSurv] = truth.childPdgCode(isec);
          m_fp_child[isec+iPrimSurv] = secMom.mag()/m_p_mother;
          m_oa_child[isec+iPrimSurv] = secMom*mom/m_p_mother/secMom.mag();
        }
        pbal -= secMom;
      }
      
      m_vtx_p_diff = pbal.mag();
      m_vtx_plong_diff = pbal*mom/m_p_mother;
      m_vtx_pperp_diff = std::sqrt(m_vtx_p_diff*m_vtx_p_diff-m_vtx_plong_diff*m_vtx_plong_diff);
      
      m_interactions->Fill();
      
      // reset the radiation length
      if (m_process==3) m_radLength = 0.;
    }
      }
      
      // crossing subdetector boundary ?
      G4VPhysicalVolume *preVol=preStep->GetPhysicalVolume();
      G4VPhysicalVolume *postVol=postStep->GetPhysicalVolume();
      
      if (postVol==0) {
    // left world -save info
    m_scEnd = 0;
    m_geoID = geoID;
    m_dt = track->GetLocalTime();
    
    m_pth = track->GetVertexMomentumDirection().theta();
    m_pph = track->GetVertexMomentumDirection().phi();
    m_p   = track->GetVertexKineticEnergy();
    m_eloss = track->GetKineticEnergy()-m_p;
    
    m_thIn= track->GetVertexPosition().theta();
    m_phIn= track->GetVertexPosition().phi();
    m_dIn=  track->GetVertexPosition().mag();
    
    m_thEnd=postStep->GetPosition().theta();
    m_phEnd=postStep->GetPosition().phi();
    m_dEnd=postStep->GetPosition().mag();
    
    m_particles->Fill();
    m_X0 = 0.;
    m_L0 = 0.;
    m_wZ = 0.;
    m_wzOaTr = 0.;
    
    m_radloss = 0.;
    m_ionloss = 0.;
    
    return;
      }
      
      // if particle killed, save the info
      if ( !trackIsAlive ) {
    m_scEnd = process? process->GetProcessSubType() : -1;
    m_geoID = geoID;
    m_dt = track->GetLocalTime();
    
    m_pth = track->GetVertexMomentumDirection().theta();
    m_pph = track->GetVertexMomentumDirection().phi();
    m_p   = track->GetVertexKineticEnergy();
    m_eloss = track->GetKineticEnergy()-m_p;
    
    m_thIn= track->GetVertexPosition().theta();
    m_phIn= track->GetVertexPosition().phi();
    m_dIn=  track->GetVertexPosition().mag();
    
    m_thEnd=postStep->GetPosition().theta();
    m_phEnd=postStep->GetPosition().phi();
    m_dEnd=postStep->GetPosition().mag();
    
    m_particles->Fill();
    m_X0 = 0.;
    m_L0 = 0.;
    m_wZ = 0.;
    m_radloss = 0.;
    m_ionloss = 0.;
    m_wzOaTr = 0.;
    m_radLength = 0.;
      }
      
      if ( preVol==postVol ) return;      // assume boundaries do not cross G4Volumes
      
      // Detector boundaries defined by central GeoIDSvc
      const G4ThreeVector &postPos  = postStep->GetPosition();
      const G4ThreeVector &postMom  = postStep->GetMomentum();
      AtlasDetDescr::AtlasRegion  nextGeoID = m_geoIDSvcQuick->identifyNextGeoID( postPos.x(),
                                          postPos.y(),
                                          postPos.z(),
                                          postMom.x(),
                                          postMom.y(),
                                          postMom.z() );
      
      // save info if leaving the subdetector
      if ( nextGeoID == geoID) {
    ATH_MSG_DEBUG("track stays inside "<<geoID);
      } else {
        ATH_MSG_DEBUG("track moves from "<<geoID<<" to "<<nextGeoID);
 
        // Don't save if doing a zero step ?
        //if (aStep->GetTrack()->GetTrackLength() == 0.) {
        if (aStep->GetStepLength() == 0.) {
          return;
        }
 
        // don't change geometry assignment for validation ntuple
        m_geoID = geoID;
 
        //4ParticleDefinition* particleDefinition = track->GetDefinition();
 
        const G4ThreeVector& g4pos = track->GetPosition();
        //const double gTime = track->GetGlobalTime();
        const HepGeom::Point3D<double> position(g4pos.x(),g4pos.y(),g4pos.z());
 
        G4ThreeVector g4mom = track->GetMomentum();
        const HepGeom::Vector3D<double> momentum(g4mom.x(),g4mom.y(),g4mom.z());
 
        bool dead=false;
        if (track->GetTrackStatus()==fStopAndKill) {
          dead=true;
        }
 
        if (!dead) {
 
          // track info
          //VTrackInformation * trackInfo= static_cast<VTrackInformation*>(track->GetUserInformation());
          m_scEnd = 0;
          m_dt = track->GetLocalTime();
 
          m_pth = track->GetVertexMomentumDirection().theta();
          m_pph = track->GetVertexMomentumDirection().phi();
          m_p   = track->GetVertexKineticEnergy();
          m_eloss = track->GetKineticEnergy()-m_p;
 
          m_thIn= track->GetVertexPosition().theta();
          m_phIn= track->GetVertexPosition().phi();
          m_dIn=  track->GetVertexPosition().mag();
 
          m_thEnd=postStep->GetPosition().theta();
          m_phEnd=postStep->GetPosition().phi();
          m_dEnd=postStep->GetPosition().mag();
 
          m_particles->Fill();
          m_X0 = 0.;
          m_L0 = 0.;
          m_wZ = 0.;
          m_radloss = 0.;
          m_ionloss = 0.;
          m_wzOaTr= 0.;
        }
        m_X0=0.;
 
      }
      
    }

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::iGeant4::PhysicsValidationUserAction::m_config

private

Definition at line 66 of file PhysicsValidationUserAction.h.

m_currentTrack

int G4UA::iGeant4::PhysicsValidationUserAction::m_currentTrack

private

Definition at line 117 of file PhysicsValidationUserAction.h.

m_dEnd

float G4UA::iGeant4::PhysicsValidationUserAction::m_dEnd

private

Definition at line 89 of file PhysicsValidationUserAction.h.

m_dIn

float G4UA::iGeant4::PhysicsValidationUserAction::m_dIn

private

Definition at line 86 of file PhysicsValidationUserAction.h.

m_dt

float G4UA::iGeant4::PhysicsValidationUserAction::m_dt

private

Definition at line 93 of file PhysicsValidationUserAction.h.

m_eloss

float G4UA::iGeant4::PhysicsValidationUserAction::m_eloss

private

Definition at line 80 of file PhysicsValidationUserAction.h.

m_fp_child

float G4UA::iGeant4::PhysicsValidationUserAction::m_fp_child[MAXCHILDREN] ={}

private

Definition at line 110 of file PhysicsValidationUserAction.h.

{};      // fraction of incoming momentum

m_gen

int G4UA::iGeant4::PhysicsValidationUserAction::m_gen

private

Definition at line 75 of file PhysicsValidationUserAction.h.

m_gen_mother

int G4UA::iGeant4::PhysicsValidationUserAction::m_gen_mother

private

Definition at line 98 of file PhysicsValidationUserAction.h.

m_geoID

int G4UA::iGeant4::PhysicsValidationUserAction::m_geoID

private

Definition at line 76 of file PhysicsValidationUserAction.h.

m_geoIDSvcQuick

ISF::IGeoIDSvc* G4UA::iGeant4::PhysicsValidationUserAction::m_geoIDSvcQuick

private

access to the central ISF GeoID serice

quickaccess avoiding gaudi ovehead

Definition at line 69 of file PhysicsValidationUserAction.h.

m_imsg

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

mutableprivateinherited

MessageSvc pointer.

Definition at line 135 of file AthMessaging.h.

{ nullptr };

m_interactions

TTree* G4UA::iGeant4::PhysicsValidationUserAction::m_interactions

private

ROOT tree containing vertex info.

Definition at line 95 of file PhysicsValidationUserAction.h.

m_ionloss

float G4UA::iGeant4::PhysicsValidationUserAction::m_ionloss

private

Definition at line 82 of file PhysicsValidationUserAction.h.

m_L0

float G4UA::iGeant4::PhysicsValidationUserAction::m_L0

private

Definition at line 91 of file PhysicsValidationUserAction.h.

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_minHistoryDepth

int G4UA::iGeant4::PhysicsValidationUserAction::m_minHistoryDepth

private

Definition at line 115 of file PhysicsValidationUserAction.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_nChild

int G4UA::iGeant4::PhysicsValidationUserAction::m_nChild

private

Definition at line 99 of file PhysicsValidationUserAction.h.

m_nm

std::string AthMessaging::m_nm

privateinherited

Message source name.

Definition at line 129 of file AthMessaging.h.

m_oa_child

float G4UA::iGeant4::PhysicsValidationUserAction::m_oa_child[MAXCHILDREN] ={}

private

Definition at line 111 of file PhysicsValidationUserAction.h.

{};      // opening angle wrt the mother

m_p

float G4UA::iGeant4::PhysicsValidationUserAction::m_p

private

Definition at line 79 of file PhysicsValidationUserAction.h.

m_p_mother

float G4UA::iGeant4::PhysicsValidationUserAction::m_p_mother

private

Definition at line 107 of file PhysicsValidationUserAction.h.

m_particles

TTree* G4UA::iGeant4::PhysicsValidationUserAction::m_particles

private

ROOT tree containing track info.

Definition at line 71 of file PhysicsValidationUserAction.h.

m_pdg

int G4UA::iGeant4::PhysicsValidationUserAction::m_pdg

private

Definition at line 72 of file PhysicsValidationUserAction.h.

m_pdg_child

int G4UA::iGeant4::PhysicsValidationUserAction::m_pdg_child[MAXCHILDREN] ={}

private

Definition at line 109 of file PhysicsValidationUserAction.h.

{};     // decay product pdg code

m_pdg_mother

int G4UA::iGeant4::PhysicsValidationUserAction::m_pdg_mother

private

Definition at line 97 of file PhysicsValidationUserAction.h.

m_phEnd

float G4UA::iGeant4::PhysicsValidationUserAction::m_phEnd

private

Definition at line 88 of file PhysicsValidationUserAction.h.

m_phIn

float G4UA::iGeant4::PhysicsValidationUserAction::m_phIn

private

Definition at line 85 of file PhysicsValidationUserAction.h.

m_pph

float G4UA::iGeant4::PhysicsValidationUserAction::m_pph

private

Definition at line 78 of file PhysicsValidationUserAction.h.

m_process

int G4UA::iGeant4::PhysicsValidationUserAction::m_process

private

Definition at line 96 of file PhysicsValidationUserAction.h.

m_pth

float G4UA::iGeant4::PhysicsValidationUserAction::m_pth

private

Definition at line 77 of file PhysicsValidationUserAction.h.

m_radLength

float G4UA::iGeant4::PhysicsValidationUserAction::m_radLength

private

Definition at line 108 of file PhysicsValidationUserAction.h.

m_radloss

float G4UA::iGeant4::PhysicsValidationUserAction::m_radloss

private

Definition at line 81 of file PhysicsValidationUserAction.h.

m_scEnd

int G4UA::iGeant4::PhysicsValidationUserAction::m_scEnd

private

Definition at line 74 of file PhysicsValidationUserAction.h.

m_scIn

int G4UA::iGeant4::PhysicsValidationUserAction::m_scIn

private

Definition at line 73 of file PhysicsValidationUserAction.h.

m_thEnd

float G4UA::iGeant4::PhysicsValidationUserAction::m_thEnd

private

Definition at line 87 of file PhysicsValidationUserAction.h.

m_thIn

float G4UA::iGeant4::PhysicsValidationUserAction::m_thIn

private

Definition at line 84 of file PhysicsValidationUserAction.h.

m_trackGenMap

std::map<int, int> G4UA::iGeant4::PhysicsValidationUserAction::m_trackGenMap

private

Definition at line 118 of file PhysicsValidationUserAction.h.

m_volumeOffset

int G4UA::iGeant4::PhysicsValidationUserAction::m_volumeOffset

private

Definition at line 114 of file PhysicsValidationUserAction.h.

m_vtx_dist

float G4UA::iGeant4::PhysicsValidationUserAction::m_vtx_dist

private

Definition at line 100 of file PhysicsValidationUserAction.h.

m_vtx_e_diff

float G4UA::iGeant4::PhysicsValidationUserAction::m_vtx_e_diff

private

Definition at line 103 of file PhysicsValidationUserAction.h.

m_vtx_p_diff

float G4UA::iGeant4::PhysicsValidationUserAction::m_vtx_p_diff

private

Definition at line 104 of file PhysicsValidationUserAction.h.

m_vtx_phi

float G4UA::iGeant4::PhysicsValidationUserAction::m_vtx_phi

private

Definition at line 102 of file PhysicsValidationUserAction.h.

m_vtx_plong_diff

float G4UA::iGeant4::PhysicsValidationUserAction::m_vtx_plong_diff

private

Definition at line 105 of file PhysicsValidationUserAction.h.

m_vtx_pperp_diff

float G4UA::iGeant4::PhysicsValidationUserAction::m_vtx_pperp_diff

private

Definition at line 106 of file PhysicsValidationUserAction.h.

m_vtx_theta

float G4UA::iGeant4::PhysicsValidationUserAction::m_vtx_theta

private

Definition at line 101 of file PhysicsValidationUserAction.h.

m_wZ

float G4UA::iGeant4::PhysicsValidationUserAction::m_wZ

private

Definition at line 92 of file PhysicsValidationUserAction.h.

m_wzOaTr

float G4UA::iGeant4::PhysicsValidationUserAction::m_wzOaTr

private

Definition at line 83 of file PhysicsValidationUserAction.h.

m_X0

float G4UA::iGeant4::PhysicsValidationUserAction::m_X0

private

Definition at line 90 of file PhysicsValidationUserAction.h.

---

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

• PhysicsValidationUserAction.h
• PhysicsValidationUserAction.cxx