G4VHadDecayAlgorithm Class Reference

#include <G4VHadDecayAlgorithm.hh>

Inheritance diagram for G4VHadDecayAlgorithm:

Public Member Functions
	G4VHadDecayAlgorithm (const G4String &algName, G4int verbose=0)
void	Generate (G4double initialMass, const std::vector< G4double > &masses, std::vector< G4LorentzVector > &finalState)
const G4String &	GetName () const
G4int	GetVerboseLevel () const
virtual void	SetVerboseLevel (G4int verbose)
virtual	~G4VHadDecayAlgorithm ()

Protected Member Functions
virtual void	GenerateMultiBody (G4double initialMass, const std::vector< G4double > &masses, std::vector< G4LorentzVector > &finalState)=0
virtual void	GenerateTwoBody (G4double initialMass, const std::vector< G4double > &masses, std::vector< G4LorentzVector > &finalState)=0
virtual G4bool	IsDecayAllowed (G4double initialMass, const std::vector< G4double > &masses) const
void	PrintVector (const std::vector< G4double > &v, const G4String &name, std::ostream &os) const
G4double	TwoBodyMomentum (G4double M0, G4double M1, G4double M2) const
G4double	UniformPhi () const
G4double	UniformTheta () const

Private Attributes
G4String	name
G4int	verboseLevel

Detailed Description

Definition at line 42 of file G4VHadDecayAlgorithm.hh.

Constructor & Destructor Documentation

G4VHadDecayAlgorithm()

G4VHadDecayAlgorithm::G4VHadDecayAlgorithm ( const G4String &  algName, G4int  verbose = 0  )

inline

Definition at line 44 of file G4VHadDecayAlgorithm.hh.

    : name(algName), verboseLevel(verbose) {;}

~G4VHadDecayAlgorithm()

virtual G4VHadDecayAlgorithm::~G4VHadDecayAlgorithm ( )

inlinevirtual

Definition at line 46 of file G4VHadDecayAlgorithm.hh.

{;}

Member Function Documentation

Generate()

void G4VHadDecayAlgorithm::Generate	(	G4double	initialMass,
		const std::vector< G4double > &	masses,
		std::vector< G4LorentzVector > &	finalState
	)

Definition at line 48 of file G4VHadDecayAlgorithm.cc.

                                                            {
  if (verboseLevel) G4cout << GetName() << "::Generate" << G4endl;
 
  // Initialization and sanity check
  finalState.clear();
  if (!IsDecayAllowed(initialMass, masses)) return;
 
  // Allow different procedures for two-body or N-body distributions
  if (masses.size() == 2U) 
    GenerateTwoBody(initialMass, masses, finalState);
  else 
    GenerateMultiBody(initialMass, masses, finalState);
}

References G4cout, G4endl, GenerateMultiBody(), GenerateTwoBody(), GetName(), IsDecayAllowed(), and verboseLevel.

Referenced by G4HadDecayGenerator::Generate().

GenerateMultiBody()

virtual void G4VHadDecayAlgorithm::GenerateMultiBody ( G4double  initialMass, const std::vector< G4double > &  masses, std::vector< G4LorentzVector > &  finalState  )

protectedpure virtual

Implemented in G4CascadeFinalStateAlgorithm, G4HadPhaseSpaceGenbod, G4HadPhaseSpaceKopylov, and G4HadPhaseSpaceNBodyAsai.

Referenced by Generate().

GenerateTwoBody()

virtual void G4VHadDecayAlgorithm::GenerateTwoBody ( G4double  initialMass, const std::vector< G4double > &  masses, std::vector< G4LorentzVector > &  finalState  )

protectedpure virtual

Implemented in G4CascadeFinalStateAlgorithm, and G4VHadPhaseSpaceAlgorithm.

Referenced by Generate().

GetName()

const G4String & G4VHadDecayAlgorithm::GetName ( ) const

inline

Definition at line 56 of file G4VHadDecayAlgorithm.hh.

{ return name; }

References name.

Referenced by G4HadPhaseSpaceGenbod::AcceptEvent(), G4HadPhaseSpaceGenbod::AccumulateFinalState(), G4CascadeFinalStateAlgorithm::ChooseGenerators(), G4HadPhaseSpaceGenbod::ComputeWeight(), G4HadPhaseSpaceGenbod::ComputeWeightScale(), G4CascadeFinalStateAlgorithm::Configure(), G4CascadeFinalStateAlgorithm::FillDirections(), G4CascadeFinalStateAlgorithm::FillDirManyBody(), G4CascadeFinalStateAlgorithm::FillDirThreeBody(), G4HadPhaseSpaceGenbod::FillEnergySteps(), G4CascadeFinalStateAlgorithm::FillMagnitudes(), G4HadPhaseSpaceGenbod::FillRandomBuffer(), G4CascadeFinalStateAlgorithm::FillUsingKopylov(), G4HadDecayGenerator::G4HadDecayGenerator(), Generate(), G4CascadeFinalStateAlgorithm::GenerateCosTheta(), G4HadPhaseSpaceGenbod::GenerateMomenta(), G4CascadeFinalStateAlgorithm::GenerateMultiBody(), G4HadPhaseSpaceGenbod::GenerateMultiBody(), G4HadPhaseSpaceKopylov::GenerateMultiBody(), G4HadPhaseSpaceNBodyAsai::GenerateMultiBody(), G4CascadeFinalStateAlgorithm::GenerateTwoBody(), G4HadDecayGenerator::GetAlgorithmName(), G4HadPhaseSpaceGenbod::Initialize(), IsDecayAllowed(), G4CascadeFinalStateAlgorithm::satisfyTriangle(), G4CascadeFinalStateAlgorithm::SaveKinematics(), and TwoBodyMomentum().

GetVerboseLevel()

G4int G4VHadDecayAlgorithm::GetVerboseLevel ( ) const

inline

Definition at line 55 of file G4VHadDecayAlgorithm.hh.

{ return verboseLevel; }

References verboseLevel.

Referenced by G4HadPhaseSpaceGenbod::AcceptEvent(), G4HadPhaseSpaceGenbod::AccumulateFinalState(), G4CascadeFinalStateAlgorithm::ChooseGenerators(), G4HadPhaseSpaceGenbod::ComputeWeight(), G4HadPhaseSpaceGenbod::ComputeWeightScale(), G4CascadeFinalStateAlgorithm::Configure(), G4CascadeFinalStateAlgorithm::FillDirections(), G4CascadeFinalStateAlgorithm::FillDirManyBody(), G4CascadeFinalStateAlgorithm::FillDirThreeBody(), G4HadPhaseSpaceGenbod::FillEnergySteps(), G4CascadeFinalStateAlgorithm::FillMagnitudes(), G4HadPhaseSpaceGenbod::FillRandomBuffer(), G4CascadeFinalStateAlgorithm::FillUsingKopylov(), G4CascadeFinalStateAlgorithm::GenerateCosTheta(), G4HadPhaseSpaceGenbod::GenerateMomenta(), G4CascadeFinalStateAlgorithm::GenerateMultiBody(), G4HadPhaseSpaceGenbod::GenerateMultiBody(), G4HadPhaseSpaceKopylov::GenerateMultiBody(), G4HadPhaseSpaceNBodyAsai::GenerateMultiBody(), G4CascadeFinalStateAlgorithm::GenerateTwoBody(), G4VHadPhaseSpaceAlgorithm::GenerateTwoBody(), G4HadPhaseSpaceGenbod::Initialize(), G4CascadeFinalStateAlgorithm::satisfyTriangle(), and G4CascadeFinalStateAlgorithm::SaveKinematics().

IsDecayAllowed()

G4bool G4VHadDecayAlgorithm::IsDecayAllowed ( G4double  initialMass, const std::vector< G4double > &  masses  ) const

protectedvirtual

Definition at line 67 of file G4VHadDecayAlgorithm.cc.

                                                    {
  G4bool okay =
    (initialMass > 0. && masses.size() >= 2 &&
     initialMass >= std::accumulate(masses.begin(),masses.end(),0.));
 
  if (verboseLevel) {
    G4cout << GetName() << "::IsDecayAllowed? initialMass " << initialMass
       << " " << masses.size() << " masses sum "
       << std::accumulate(masses.begin(),masses.end(),0.) << G4endl;
 
    if (verboseLevel>1) PrintVector(masses," ",G4cout);
 
    G4cout << " Returning " << okay << G4endl;
  }
 
  return okay;
}

References G4cout, G4endl, GetName(), PrintVector(), and verboseLevel.

Referenced by Generate().

PrintVector()

void G4VHadDecayAlgorithm::PrintVector ( const std::vector< G4double > &  v, const G4String &  name, std::ostream &  os  ) const

protected

Definition at line 121 of file G4VHadDecayAlgorithm.cc.

                                                     {
  os << " " << vname << "(" << v.size() << ") ";
  std::copy(v.begin(), v.end(), std::ostream_iterator<G4double>(os, " "));
  os << std::endl;
}

References field_utils::copy().

Referenced by G4HadPhaseSpaceGenbod::FillEnergySteps(), G4HadPhaseSpaceGenbod::FillRandomBuffer(), G4HadPhaseSpaceNBodyAsai::GenerateMultiBody(), G4HadPhaseSpaceGenbod::Initialize(), and IsDecayAllowed().

SetVerboseLevel()

virtual void G4VHadDecayAlgorithm::SetVerboseLevel ( G4int  verbose )

inlinevirtual

Reimplemented in G4CascadeFinalStateAlgorithm.

Definition at line 54 of file G4VHadDecayAlgorithm.hh.

{ verboseLevel = verbose; }

References verboseLevel.

Referenced by G4CascadeFinalStateAlgorithm::SetVerboseLevel(), and G4HadDecayGenerator::SetVerboseLevel().

TwoBodyMomentum()

G4double G4VHadDecayAlgorithm::TwoBodyMomentum ( G4double  M0, G4double  M1, G4double  M2  ) const

protected

Definition at line 90 of file G4VHadDecayAlgorithm.cc.

                                              {
  G4double PSQ = (M0+M1+M2)*(M0+M1-M2)*(M0-M1+M2)*(M0-M1-M2);
  if (PSQ < 0.) {
    G4cout << GetName() << ":  problem of decay of M(GeV) " << M0/GeV 
       << " to M1(GeV) " << M1/GeV << " and M2(GeV) " << M2/GeV
       << " PSQ(MeV) " << PSQ/MeV << " < 0" << G4endl;
    // exception only if the problem is numerically significant
    if (PSQ < -CLHEP::eV) {
      throw G4HadronicException(__FILE__, __LINE__,"Error in decay kinematics");
    }
 
    PSQ = 0.;
  }
 
  return std::sqrt(PSQ)/(2.*M0);
}

References CLHEP::eV, G4cout, G4endl, GetName(), GeV, and MeV.

Referenced by G4HadPhaseSpaceGenbod::ComputeWeightScale(), G4HadPhaseSpaceGenbod::FillEnergySteps(), G4CascadeFinalStateAlgorithm::FillUsingKopylov(), G4HadPhaseSpaceKopylov::GenerateMultiBody(), G4HadPhaseSpaceNBodyAsai::GenerateMultiBody(), G4CascadeFinalStateAlgorithm::GenerateTwoBody(), and G4VHadPhaseSpaceAlgorithm::GenerateTwoBody().

UniformPhi()

G4double G4VHadDecayAlgorithm::UniformPhi ( ) const

protected

Definition at line 114 of file G4VHadDecayAlgorithm.cc.

                                                {
  return twopi*G4UniformRand();
}

References G4UniformRand, and twopi.

Referenced by G4CascadeFinalStateAlgorithm::FillUsingKopylov(), G4CascadeFinalStateAlgorithm::GenerateTwoBody(), and G4VHadPhaseSpaceAlgorithm::UniformVector().

UniformTheta()

G4double G4VHadDecayAlgorithm::UniformTheta ( ) const

protected

Definition at line 110 of file G4VHadDecayAlgorithm.cc.

                                                  {
  return std::acos(2.0*G4UniformRand() - 1.0);
}

References G4UniformRand.

Referenced by G4CascadeFinalStateAlgorithm::FillUsingKopylov(), and G4VHadPhaseSpaceAlgorithm::UniformVector().

Field Documentation

name

G4String G4VHadDecayAlgorithm::name

private

Definition at line 85 of file G4VHadDecayAlgorithm.hh.

Referenced by source.g4viscp.G4Scene::create_scene(), GetName(), mcscore.MCParticle::printout(), and source.g4viscp.G4Scene::update_scene().

verboseLevel

G4int G4VHadDecayAlgorithm::verboseLevel

private

Definition at line 86 of file G4VHadDecayAlgorithm.hh.

Referenced by Generate(), GetVerboseLevel(), IsDecayAllowed(), and SetVerboseLevel().

---

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

• source/processes/hadronic/util/include/G4VHadDecayAlgorithm.hh
• source/processes/hadronic/util/src/G4VHadDecayAlgorithm.cc