NeutronFastSim Class Reference

#include <NeutronFastSim.h>

Inheritance diagram for NeutronFastSim:

Collaboration diagram for NeutronFastSim:

Public Member Functions
	NeutronFastSim (const std::string &name, G4Region *region, const std::string &fsSDname, const double etaCut, const double timeCut)
	~NeutronFastSim ()
G4bool	IsApplicable (const G4ParticleDefinition &) override final
virtual G4bool	ModelTrigger (const G4FastTrack &) override final
void	DoIt (const G4FastTrack &, G4FastStep &) override final

Protected Attributes
G4double	m_Energy
TrackFastSimSD *	m_fsSD
bool	m_init
std::string	m_fsSDname
double	m_etaCut
double	m_timeCut

Detailed Description

Definition at line 17 of file NeutronFastSim.h.

Constructor & Destructor Documentation

NeutronFastSim()

NeutronFastSim::NeutronFastSim	(	const std::string &	name,
		G4Region *	region,
		const std::string &	fsSDname,
		const double	etaCut,
		const double	timeCut )

Definition at line 16 of file NeutronFastSim.cxx.

  : G4VFastSimulationModel(name, region)
  , m_Energy(5)
  , m_fsSD(0)
  , m_init(false)
  , m_fsSDname(fsSDname)
  , m_etaCut(etaCut)
  , m_timeCut(timeCut)
{
}

~NeutronFastSim()

NeutronFastSim::~NeutronFastSim ( )

inline

Definition at line 22 of file NeutronFastSim.h.

{}

Member Function Documentation

DoIt()

void NeutronFastSim::DoIt ( const G4FastTrack & fastTrack, G4FastStep & fastStep )

finaloverride

Definition at line 66 of file NeutronFastSim.cxx.

{
  if (m_fsSD) m_fsSD->WriteTrack( fastTrack.GetPrimaryTrack() , false , false );
  fastStep.KillPrimaryTrack();
}

IsApplicable()

G4bool NeutronFastSim::IsApplicable ( const G4ParticleDefinition & )

finaloverride

Definition at line 27 of file NeutronFastSim.cxx.

{
  if (!m_init){
    m_init = true;
 
    G4SDManager *sdm = G4SDManager::GetSDMpointer();
    G4VSensitiveDetector * vsd = sdm->FindSensitiveDetector( m_fsSDname );
    if (!vsd) {
      G4cout << "NeutronFastSim::IsApplicable WARNING Could not get TrackFastSimSD sensitive detector.  If you are not writing track records this is expected." << G4endl;
    } else {
      m_fsSD = dynamic_cast<TrackFastSimSD*>(vsd);
      if (!m_fsSD) {
        G4cout << "NeutronFastSim::IsApplicable WARNING Could not cast the SD.  If you are not writing track records this is expected." << G4endl;
      }
    } // found the SD
  } // End of lazy init
  return true;
}

ModelTrigger()

G4bool NeutronFastSim::ModelTrigger ( const G4FastTrack & fastTrack )

finaloverridevirtual

Definition at line 46 of file NeutronFastSim.cxx.

{
  // Trigger if the neutron energy is below our threshold or if the time is over 150 ns
  if (fastTrack.GetPrimaryTrack()->GetDefinition() == G4Neutron::NeutronDefinition() ){
    return (m_Energy<0?true:fastTrack.GetPrimaryTrack()->GetKineticEnergy()<m_Energy) || fastTrack.GetPrimaryTrack()->GetGlobalTime()>m_timeCut;
  }
 
  // Not a neutron... Pick it up if the primary had eta>6.0
  TrackHelper trackHelper(fastTrack.GetPrimaryTrack());
  HepMC::ConstGenParticlePtr primaryGenParticle = trackHelper.GetPrimaryGenParticle();
  if (!primaryGenParticle) {
    return false;
  }
  if (std::abs(primaryGenParticle->momentum().eta())>m_etaCut && !HepMC::is_simulation_particle(primaryGenParticle)){
    return true;
  } else {
    return false;
  }
}

Member Data Documentation

m_Energy

G4double NeutronFastSim::m_Energy

protected

Definition at line 31 of file NeutronFastSim.h.

m_etaCut

double NeutronFastSim::m_etaCut

protected

Definition at line 35 of file NeutronFastSim.h.

m_fsSD

TrackFastSimSD* NeutronFastSim::m_fsSD

protected

Definition at line 32 of file NeutronFastSim.h.

m_fsSDname

std::string NeutronFastSim::m_fsSDname

protected

Definition at line 34 of file NeutronFastSim.h.

m_init

bool NeutronFastSim::m_init

protected

Definition at line 33 of file NeutronFastSim.h.

m_timeCut

double NeutronFastSim::m_timeCut

protected

Definition at line 36 of file NeutronFastSim.h.

---

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

• NeutronFastSim.h
• NeutronFastSim.cxx