Public Member Functions
	BertiniEventTopologyProcess (bool count_light_ions=true)

virtual bool	acceptEvent () const =0

virtual bool	acceptProjectile (const G4HadProjectile &projectile) const =0

virtual bool	acceptTarget (const G4Nucleus &targetNucleus) const =0

G4HadFinalState *	ApplyYourself (const G4HadProjectile &projectile, G4Nucleus &targetNucleus) override

void	cleanupSecondaries ()

constexpr bool	isLightIon (const int pdgCode) const
	Check if the PDG code corresponds to a light ion nucleus.

constexpr bool	skipCountingParticle (const int pdgcode) const
	Whether or not to include a particular particle type in any counting.

virtual void	incrementEventWeight (int N)

Private Attributes
bool	count_light_ions_

Detailed Description

Definition at line 34 of file BertiniEventTopologyProcess.h.

Constructor & Destructor Documentation

◆ BertiniEventTopologyProcess()

simcore::BertiniEventTopologyProcess::BertiniEventTopologyProcess ( bool count_light_ions = true )

inline

Definition at line 36 of file BertiniEventTopologyProcess.h.

37 : G4CascadeInterface{}, count_light_ions_{count_light_ions} {}

Member Function Documentation

◆ ApplyYourself()

G4HadFinalState * simcore::BertiniEventTopologyProcess::ApplyYourself	(	const G4HadProjectile &	projectile,
		G4Nucleus &	targetNucleus )

override

Definition at line 19 of file BertiniEventTopologyProcess.cxx.

                                                                 {
  int attempts{1};
  if (!acceptProjectile(projectile) || !acceptTarget(targetNucleus)) {
    // Bertini will handle the particle change on its own here
    return G4CascadeInterface::ApplyYourself(projectile, targetNucleus);
  }
 
  while (true) {
    theParticleChange.Clear();
    theParticleChange.SetStatusChange(stopAndKill);
    G4CascadeInterface::ApplyYourself(projectile, targetNucleus);
    if (acceptEvent()) {
      incrementEventWeight(attempts);
      return &theParticleChange;
    }
    attempts++;
    cleanupSecondaries();
  }
}

◆ cleanupSecondaries()

void simcore::BertiniEventTopologyProcess::cleanupSecondaries ( )

Definition at line 4 of file BertiniEventTopologyProcess.cxx.

                                                     {
  int secondaries{theParticleChange.GetNumberOfSecondaries()};
  // Geant4 won't clean up this memory for us by default
  for (int i{0}; i < secondaries; ++i) {
    // disable clang's use-after-free warning for this delete
    // as far as we (@tomeichlersmith and @EinarElen) can tell,
    // this is correct and clang is just getting confused since
    // Geant4's inlined delete method is written in a wacky way
#ifndef __clang_analyzer__
    auto secondary{theParticleChange.GetSecondary(i)->GetParticle()};
    delete secondary;
#endif
  }
}

◆ incrementEventWeight()

virtual void simcore::BertiniEventTopologyProcess::incrementEventWeight ( int N )

inlinevirtual

Definition at line 126 of file BertiniEventTopologyProcess.h.

                                           {
    auto event_info{static_cast<UserEventInformation*>(
        G4EventManager::GetEventManager()->GetUserInformation())};
    event_info->incWeight(1. / N);
  }

◆ isLightIon()

bool simcore::BertiniEventTopologyProcess::isLightIon ( const int pdgCode ) const

inlineconstexpr

Check if the PDG code corresponds to a light ion nucleus.

Nuclear PDG codes are given by ±10LZZZAAAI So to find the atomic number, we first divide by 10 (to lose the I-component) and then take the modulo with 1000.

Definition at line 93 of file BertiniEventTopologyProcess.h.

                                                     {
    //
    if (pdgCode > 1000000000) {
      // Check if the atomic number is less than or equal to 4
      return ((pdgCode / 10) % 1000) <= 4;
    }
    return false;
  }

Referenced by skipCountingParticle().

◆ skipCountingParticle()

bool simcore::BertiniEventTopologyProcess::skipCountingParticle ( const int pdgcode ) const

inlineconstexpr

Whether or not to include a particular particle type in any counting.

Unless count_light_ions_ is set, we don't count anything with a nuclear PDG code. This is consistent with the counting behaviour used in the PhotoNuclearDQM.

If count_light_ions_ is set, we also match PDG codes for nuclei with atomic number < 4.

See also: isLightIon

Definition at line 114 of file BertiniEventTopologyProcess.h.

                                                               {
    return !(pdgcode < 10000 || (count_light_ions_ && isLightIon(pdgcode)));
  }

References isLightIon().

Member Data Documentation

◆ count_light_ions_

bool simcore::BertiniEventTopologyProcess::count_light_ions_

private

Definition at line 133 of file BertiniEventTopologyProcess.h.

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

SimCore/include/SimCore/PhotoNuclearModels/BertiniEventTopologyProcess.h
SimCore/src/SimCore/PhotoNuclearModels/BertiniEventTopologyProcess.cxx

Public Member Functions

Private Attributes