iFatras::G4ParticleDecayHelper Class Reference

#include <G4ParticleDecayHelper.h>

Inheritance diagram for iFatras::G4ParticleDecayHelper:

Collaboration diagram for iFatras::G4ParticleDecayHelper:

Public Member Functions
	G4ParticleDecayHelper (const std::string &, const std::string &, const IInterface *)
	AlgTool constructor for ParticleDecayHelper. More...
	~G4ParticleDecayHelper ()
	Destructor. More...
StatusCode	initialize ()
	AlgTool initailize method. More...
StatusCode	finalize ()
	AlgTool finalize method. More...
double	freePath (const ISF::ISFParticle &isp) const
	free path estimator (-1 for stable particle) More...
void	decay (const ISF::ISFParticle &isp, const Amg::Vector3D &vertex, const Amg::Vector3D &mom, double timeStamp=0) const
	decay handling secondaries More...
std::vector< ISF::ISFParticle * >	decayParticle (const ISF::ISFParticle &parent, const Amg::Vector3D &vertex, const Amg::Vector3D &mom, double timeStamp=0) const
	decay More...

Private Member Functions
bool	initG4RunManager () const
	initialize G4RunManager on first call if not done by then More...
void	handleDecayParticles (const ISF::ISFParticle &isp, const ISF::ISFParticleVector &children) const
	fill decay products: into broker svc, truth svc More...

Private Attributes
ServiceHandle< ISF::IParticleBroker >	m_particleBroker
	Broker Svc for ISF particles. More...
ServiceHandle< ISF::ITruthSvc >	m_truthRecordSvc
	Truth Svc for truth tree. More...
ServiceHandle< IAtRndmGenSvc >	m_rndmSvc
	Random Svc. More...
CLHEP::HepRandomEngine *	m_randomEngine
	Random engine (updated to streams) More...
std::string	m_randomEngineName
	Name of the random number stream. More...
std::string	m_G4RandomEngineName
	Name of the random number stream for G4 tools. More...
ToolHandle< ISF::IG4RunManagerHelper >	m_g4RunManagerHelper
	G4RunManager needs to be initialized before G4 tables are accessed. More...
ToolHandle< PDGToG4Particle >	m_pdgToG4Conv
	Handle for the PDGToG4Particle converter tool. More...
bool	m_validationMode
	Validation output with histogram service. More...
ToolHandle< IPhysicsValidationTool >	m_validationTool
	the ntuple More...

Detailed Description

This Tool is a helper tool wrapping the IParticleDecayCreator, and organised the flight distance to

Author

Andre.nosp@m.as.S.nosp@m.alzbu.nosp@m.rger.nosp@m.@cern.nosp@m..ch, Joerg.nosp@m..Mec.nosp@m.hnich.nosp@m.@cer.nosp@m.n.ch

Definition at line 45 of file G4ParticleDecayHelper.h.

Constructor & Destructor Documentation

G4ParticleDecayHelper()

iFatras::G4ParticleDecayHelper::G4ParticleDecayHelper	(	const std::string &	t,
		const std::string &	n,
		const IInterface *	p
	)

AlgTool constructor for ParticleDecayHelper.

Definition at line 45 of file G4ParticleDecayHelper.cxx.

                                                                                                              :    
   base_class(t,n,p),
   m_particleBroker("ISF_ParticleBroker", n),
   m_truthRecordSvc("ISF_TruthRecordSvc", n),
   m_rndmSvc("AtDSFMTGenSvc", n),
   m_randomEngine(0),
   m_randomEngineName("FatrasRnd"),
   m_G4RandomEngineName("FatrasG4"),
   m_g4RunManagerHelper("iGeant4::G4RunManagerHelper/G4RunManagerHelper"),
   m_pdgToG4Conv("iFatras::PDGToG4Particle/PDGToG4ParticleConverter"),
   m_validationMode(false),
   m_validationTool("")
{
  // ISF: truth and broker service
  declareProperty("ParticleBroker",           m_particleBroker,       "ISF Particle Broker Svc");
  declareProperty("ParticleTruthSvc",         m_truthRecordSvc,       "ISF Particle Truth Svc");
  // random number initializations 
  declareProperty("RandomNumberService",      m_rndmSvc,              "Random number generator");
  declareProperty("RandomStreamName",         m_randomEngineName,     "Name of the random number stream");
  declareProperty("G4RandomStreamName",       m_G4RandomEngineName,   "Name of the random number stream for G4 tools");
  
  // tool declarations -------------------------------------------------
  declareProperty("G4RunManagerHelper"      , m_g4RunManagerHelper);
  declareProperty("PDGToG4ParticleConverter", m_pdgToG4Conv);
  // validation output section
  declareProperty( "ValidationMode",  m_validationMode );
  declareProperty( "PhysicsValidationTool",     m_validationTool  );
 
}

~G4ParticleDecayHelper()

iFatras::G4ParticleDecayHelper::~G4ParticleDecayHelper

(

)

Destructor.

Definition at line 79 of file G4ParticleDecayHelper.cxx.

{}

Member Function Documentation

decay()

void iFatras::G4ParticleDecayHelper::decay	(	const ISF::ISFParticle &	isp,
		const Amg::Vector3D &	vertex,
		const Amg::Vector3D &	mom,
		double	timeStamp = 0
	)		const

decay handling secondaries

Definition at line 195 of file G4ParticleDecayHelper.cxx.

{
 
  /*-------------------------------------------------------------------
   *  The actual decay
   *-------------------------------------------------------------------*/
  ATH_MSG_DEBUG( "[ decay ] calling G4ParticleDecayCreator with " << particleToDecay );
 
  // perform the decay 
  const ISF::ISFParticleVector decayProducts = decayParticle(particleToDecay,vertex,momentum,timeStamp);
 
  // fill them into broker & truth svc
  handleDecayParticles(particleToDecay,decayProducts); // Registers TruthIncident internally
}

decayParticle()

std::vector< ISF::ISFParticle * > iFatras::G4ParticleDecayHelper::decayParticle	(	const ISF::ISFParticle &	parent,
		const Amg::Vector3D &	vertex,
		const Amg::Vector3D &	mom,
		double	timeStamp = 0
	)		const

decay

Definition at line 272 of file G4ParticleDecayHelper.cxx.

{
  // Called from McMaterialEffectsUpdator::interact, McMaterialEffectsUpdator::interactLay and G4ParticleDecayHelper::decay
  // return vector for children
  std::vector<ISF::ISFParticle*> children;
 
  // initialize G4RunManager if not done already
  bool g4mgr = initG4RunManager();
  if (!g4mgr) ATH_MSG_WARNING(  "initialization of G4RunManager failed in G4ParticleDecayHelper" );
 
  if (!m_pdgToG4Conv) {
 
    if( m_pdgToG4Conv.retrieve().isFailure()) return children;
 
    if( msgLvl(MSG::VERBOSE)) {
      ATH_MSG_VERBOSE(  "List of particles with decay info:" );
      m_pdgToG4Conv->printListOfParticles( true);
    }
  }
 
  int pdgCode = parent.pdgCode();
    
  G4ParticleDefinition* pDef = m_pdgToG4Conv->getParticleDefinition( pdgCode);
  if( !pDef)
  {
    ATH_MSG_WARNING( "[ decay ] could not find geant4 equivalent"
          << " for particle with pdg id " << pdgCode );
    return children;
  }
  
  ATH_MSG_DEBUG(  "[ decay ] Decaying " << pDef->GetParticleName() );
  
  G4DecayTable* dt = pDef->GetDecayTable();
  if( !dt)
  {
    ATH_MSG_WARNING( "[ decay ] empty decay table for"
          << " particle with pdg id " << pdgCode );
    return children;
  }
    
  if( msgLvl(MSG::VERBOSE))
  {
    ATH_MSG_VERBOSE(  "[ decay ] Decay table:" );
    dt->DumpInfo();
  }
  
  G4VDecayChannel* channel = dt->SelectADecayChannel();
  if( !channel)
  {
    ATH_MSG_ERROR( "[ decay ] error selecting decay channel for"
          << " particle with pdg id " << pdgCode );
    return children;
  }
    
  if( msgLvl(MSG::DEBUG))
  {
    ATH_MSG_DEBUG(  "[ decay ] Decay channel:" );
    channel->DumpInfo();
  }
  
  G4DecayProducts* products = channel->DecayIt();
  if( !products)
  {
    ATH_MSG_ERROR( "[ decay ] error in decay product generation"
          << " for particle with pdg id " << pdgCode );
    return children;
  }
  
  if( msgLvl(MSG::VERBOSE))
  {
    ATH_MSG_VERBOSE(  "[ decay ] Decay products:" );
    products->DumpInfo();
  }
 
  // the parent energy
  double parentE = std::sqrt( std::pow( pDef->GetPDGMass(), 2) +
                              momentum.mag2());
  
  products->Boost( momentum.x()/parentE,
                   momentum.y()/parentE,
                   momentum.z()/parentE);
  if( msgLvl(MSG::VERBOSE))
  {
    ATH_MSG_VERBOSE( "[ decay ] Decay products after boost:" );
    products->DumpInfo();
  }
 
  G4int nProducts = products->entries();
  for( G4int i=0; i < nProducts; ++i)
  {
    G4DynamicParticle* prod = products->PopProducts();
    if( !prod)
    {
      ATH_MSG_WARNING( "[ decay ] Could not retrieve product " << i);
      continue;
    }
    
    // the decay product
    const G4ThreeVector &mom= prod->GetMomentum();
    Amg::Vector3D amgMom( mom.x(), mom.y(), mom.z() );
 
    const int status = 1 + HepMC::SIM_STATUS_THRESHOLD;
    const int id = HepMC::UNDEFINED_ID;
    ISF::ISFParticle* childParticle = new ISF::ISFParticle( vertex,
                                                            amgMom,
                                                            prod->GetMass(),
                                                            prod->GetCharge(),
                                                            prod->GetPDGcode(),
                                                            status,
                                                            timeStamp,
                                                            parent,
                                                            id );
 
    children.push_back( childParticle);
  }
  return children;
}

finalize()

StatusCode iFatras::G4ParticleDecayHelper::finalize

(

)

AlgTool finalize method.

Definition at line 146 of file G4ParticleDecayHelper.cxx.

{
  ATH_MSG_DEBUG( "finalize() successful" );
  return StatusCode::SUCCESS;
}

freePath()

double iFatras::G4ParticleDecayHelper::freePath

(

const ISF::ISFParticle &

isp

)

const

free path estimator (-1 for stable particle)

Definition at line 152 of file G4ParticleDecayHelper.cxx.

                                                                             {
  // get the particle properties    
  int pdgCode =  isp.pdgCode();
 
  // initialize G4RunManager if not done already
  bool g4mgr = initG4RunManager();
  if (!g4mgr) ATH_MSG_WARNING(  "initialization of G4RunManager failed in G4ParticleDecayHelper" );
 
  if (!m_pdgToG4Conv && m_pdgToG4Conv.retrieve().isFailure()) return -1.;
 
  G4ParticleDefinition* pDef = m_pdgToG4Conv->getParticleDefinition(pdgCode);
  
  // for the processing    
  bool isStable = true;
 
  if( !pDef )      {
    ATH_MSG_WARNING( "[ decay ] could not find particle properties" << " for PDG code " << pdgCode );
    return -1.;
  }
  // set is stable from Geant4 particle definition
  isStable = pDef->GetPDGStable();
  // set muons to stable
  if ( std::abs(pdgCode) == 13) isStable = true;
  // we have a stable particle --- return
  if (isStable) return -1.;
  
  double mass = pDef->GetPDGMass();
  // take momentum from ParticleState rather than associated truth
  CLHEP::HepLorentzVector particleMom( isp.momentum().x(),
                                       isp.momentum().y(),
                                       isp.momentum().z(),
                                       sqrt(isp.momentum().mag2()+mass*mass) );
 
  // get average lifetime
  double tau = pDef->GetPDGLifeTime()*1e-9;
  // sample the lifetime
  double lifeTime = -tau*log(CLHEP::RandFlat::shoot(m_randomEngine));
  // distance in mm
  double pathlength =
    lifeTime*s_speedOfLightSI*particleMom.gamma()*particleMom.beta();
  return pathlength;
}

handleDecayParticles()

void iFatras::G4ParticleDecayHelper::handleDecayParticles ( const ISF::ISFParticle &  isp, const ISF::ISFParticleVector &  children  ) const

private

fill decay products: into broker svc, truth svc

Definition at line 214 of file G4ParticleDecayHelper.cxx.

                                                                                                           {
   // process the decay products ---------------------------------------
  int                 process = 201;
 
  // (i) none
   if (!decayProducts.size()) {
        ATH_MSG_WARNING("[ decay ] Particle Decay Creator did not return any"
            << " decay products for particle with PDG code "
            << particle.pdgCode() );
   } else {
   // (ii) many
      std::ostringstream productSummaryString;
      productSummaryString << "[ decay ] products:";
 
      for (ISF::ISFParticle *decayProduct : decayProducts) {
        productSummaryString << "     - "  << (*decayProduct) << '\n';
        // in the validation mode, add process info
        if (m_validationMode) {
          ISF::ParticleUserInformation* validInfo = new ISF::ParticleUserInformation();
          validInfo->setProcess(process);
          if (particle.getUserInformation()) validInfo->setGeneration(particle.getUserInformation()->generation()+1);
          else validInfo->setGeneration(1);     // assume parent is a primary track
          decayProduct->setUserInformation(validInfo);
        }
        // register next geo (is current), next flavor can be defined by filter
        decayProduct->setNextGeoID( particle.nextGeoID() );
      }//loop over all decay products
 
      // register TruthIncident
      ISF::ISFTruthIncident truth( const_cast<ISF::ISFParticle&>(particle),
                                   decayProducts,
                                   process,
                                   particle.nextGeoID(),  // inherits from the parent
                                   ISF::fKillsPrimary );
      m_truthRecordSvc->registerTruthIncident( truth);
      // At this point we need to update the properties of the
      // ISFParticles produced in the interaction
      truth.updateChildParticleProperties();
 
      // simulate the tracks of the daughter particles ------- run over decay products
      for (ISF::ISFParticle *decayProduct : decayProducts) {
        // feed it the particle broker with parent information
        m_particleBroker->push(decayProduct, &particle);
      }//loop over all decay products
      ATH_MSG_VERBOSE(  productSummaryString.str() );
 
      // save info for validation
      if (m_validationMode && m_validationTool.isEnabled()) {
        Amg::Vector3D* nMom = 0;
        m_validationTool->saveISFVertexInfo(process,particle.position(),particle,particle.momentum(),nMom,decayProducts);
        delete nMom;
      }
 
    }
}

initG4RunManager()

bool iFatras::G4ParticleDecayHelper::initG4RunManager ( ) const

private

initialize G4RunManager on first call if not done by then

Definition at line 393 of file G4ParticleDecayHelper.cxx.

                                                          {
 
  // safe because the static initialization ensures that this is only called once
  static const G4RunManager* const g4runManager ATLAS_THREAD_SAFE = [&] ATLAS_NOT_THREAD_SAFE () {
    auto helper_nc = const_cast<ISF::IG4RunManagerHelper*>(m_g4RunManagerHelper.get());
    return helper_nc->fastG4RunManager();
  }();
 
  return g4runManager ? true : false;
}

initialize()

StatusCode iFatras::G4ParticleDecayHelper::initialize

(

)

AlgTool initailize method.

Definition at line 87 of file G4ParticleDecayHelper.cxx.

{
 
  ATH_MSG_DEBUG( "initialize()" );
 
  ATH_CHECK( m_validationTool.retrieve( DisableTool{ m_validationTool.empty() || !m_validationMode } ) );
  // ISF Services
  if ( m_particleBroker.retrieve().isFailure()){
      ATH_MSG_FATAL( "Could not retrieve " <<  m_particleBroker );
      return StatusCode::FAILURE;
  }
  if (m_truthRecordSvc.retrieve().isFailure()){
      ATH_MSG_FATAL( "Could not retrieve " << m_truthRecordSvc );
      return StatusCode::FAILURE;
  }
 
  // Random number service
  if ( m_rndmSvc.retrieve().isFailure() ) {
    ATH_MSG_FATAL(  "Could not retrieve " << m_rndmSvc );
    return StatusCode::FAILURE;
  }
  //Get own engine with own seeds:
  m_randomEngine = m_rndmSvc->GetEngine(m_randomEngineName);
  if (!m_randomEngine) {
    ATH_MSG_FATAL(  "Could not get random engine '" << m_randomEngineName << "'" );
    return StatusCode::FAILURE;
  }
  // G4ParticleDecayCreator uses functions from G4 (G4VDecayChannel::DecayIt() and G4DecayTable::SelectADecayChannel())
  // which use CLHEP::HepRandom directly with the default engine. Therefore we set the default engine to a special stream,
  // analog to G4AtlasAlg.
  // Set the random number generator to AtRndmGen
  CLHEP::HepRandomEngine* engine = m_rndmSvc->GetEngine(m_G4RandomEngineName);
  if (!engine) {
    ATH_MSG_FATAL(  "Could not get random engine '" << m_G4RandomEngineName << "'" );
    return StatusCode::FAILURE;
  }
  CLHEP::HepRandom::setTheEngine(engine);
 
  // get G4RunManagerHelper ( no G4RunManager initialization done )
  if (m_g4RunManagerHelper.retrieve().isFailure()) {
    ATH_MSG_FATAL(  "Could not get g4RunManagerHelper '" << m_g4RunManagerHelper << "'" );
    return StatusCode::FAILURE;
  } else {
    std::cout <<"f4dec: g4RunManagerHelper retrieved:"<< m_g4RunManagerHelper<< std::endl;
  }
 
  // Disable auto-retrieve of this tool: Geant will raise an exception
  // if this happens before the run manager is created.
  m_pdgToG4Conv.disable();
 
  ATH_MSG_DEBUG("initialize() successful");
  return StatusCode::SUCCESS;
}

Member Data Documentation

m_G4RandomEngineName

std::string iFatras::G4ParticleDecayHelper::m_G4RandomEngineName

private

Name of the random number stream for G4 tools.

Definition at line 91 of file G4ParticleDecayHelper.h.

m_g4RunManagerHelper

ToolHandle<ISF::IG4RunManagerHelper> iFatras::G4ParticleDecayHelper::m_g4RunManagerHelper

private

G4RunManager needs to be initialized before G4 tables are accessed.

Definition at line 96 of file G4ParticleDecayHelper.h.

m_particleBroker

ServiceHandle<ISF::IParticleBroker> iFatras::G4ParticleDecayHelper::m_particleBroker

private

Broker Svc for ISF particles.

Definition at line 84 of file G4ParticleDecayHelper.h.

m_pdgToG4Conv

ToolHandle<PDGToG4Particle> iFatras::G4ParticleDecayHelper::m_pdgToG4Conv

private

Handle for the PDGToG4Particle converter tool.

Definition at line 97 of file G4ParticleDecayHelper.h.

m_randomEngine

CLHEP::HepRandomEngine* iFatras::G4ParticleDecayHelper::m_randomEngine

private

Random engine (updated to streams)

Definition at line 89 of file G4ParticleDecayHelper.h.

m_randomEngineName

std::string iFatras::G4ParticleDecayHelper::m_randomEngineName

private

Name of the random number stream.

Definition at line 90 of file G4ParticleDecayHelper.h.

m_rndmSvc

ServiceHandle<IAtRndmGenSvc> iFatras::G4ParticleDecayHelper::m_rndmSvc

private

Random Svc.

Definition at line 88 of file G4ParticleDecayHelper.h.

m_truthRecordSvc

ServiceHandle<ISF::ITruthSvc> iFatras::G4ParticleDecayHelper::m_truthRecordSvc

private

Truth Svc for truth tree.

Definition at line 85 of file G4ParticleDecayHelper.h.

m_validationMode

bool iFatras::G4ParticleDecayHelper::m_validationMode

private

Validation output with histogram service.

turn validation mode on/off

Definition at line 100 of file G4ParticleDecayHelper.h.

m_validationTool

ToolHandle<IPhysicsValidationTool> iFatras::G4ParticleDecayHelper::m_validationTool

private

the ntuple

Definition at line 101 of file G4ParticleDecayHelper.h.

---

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

• G4ParticleDecayHelper.h
• G4ParticleDecayHelper.cxx