G4VHadDecayAlgorithm Class Reference

#include <G4VHadDecayAlgorithm.hh>

Inheritance diagram for G4VHadDecayAlgorithm:

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

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

Detailed Description

Definition at line 43 of file G4VHadDecayAlgorithm.hh.

Constructor & Destructor Documentation

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

inline

Definition at line 45 of file G4VHadDecayAlgorithm.hh.

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

virtual G4VHadDecayAlgorithm::~G4VHadDecayAlgorithm ( )

inlinevirtual

Definition at line 47 of file G4VHadDecayAlgorithm.hh.

{;}

Member Function Documentation

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

Definition at line 49 of file G4VHadDecayAlgorithm.cc.

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

Referenced by G4HadDecayGenerator::Generate().

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

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

protectedpure virtual

Implemented in G4CascadeFinalStateAlgorithm, G4HadPhaseSpaceKopylov, G4HadPhaseSpaceNBodyAsai, and G4HadPhaseSpaceGenbod.

Referenced by Generate().

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().

const G4String& G4VHadDecayAlgorithm::GetName ( void  ) const

inline

Definition at line 57 of file G4VHadDecayAlgorithm.hh.

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(), G4HadPhaseSpaceGenbod::GenerateMultiBody(), G4HadPhaseSpaceNBodyAsai::GenerateMultiBody(), G4HadPhaseSpaceKopylov::GenerateMultiBody(), G4CascadeFinalStateAlgorithm::GenerateMultiBody(), G4CascadeFinalStateAlgorithm::GenerateTwoBody(), G4HadDecayGenerator::GetAlgorithmName(), G4HadPhaseSpaceGenbod::Initialize(), IsDecayAllowed(), G4CascadeFinalStateAlgorithm::satisfyTriangle(), G4CascadeFinalStateAlgorithm::SaveKinematics(), and TwoBodyMomentum().

{ return name; }

G4int G4VHadDecayAlgorithm::GetVerboseLevel ( ) const

inline

Definition at line 56 of file G4VHadDecayAlgorithm.hh.

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(), G4HadPhaseSpaceGenbod::GenerateMultiBody(), G4HadPhaseSpaceKopylov::GenerateMultiBody(), G4HadPhaseSpaceNBodyAsai::GenerateMultiBody(), G4CascadeFinalStateAlgorithm::GenerateMultiBody(), G4VHadPhaseSpaceAlgorithm::GenerateTwoBody(), G4CascadeFinalStateAlgorithm::GenerateTwoBody(), G4HadPhaseSpaceGenbod::Initialize(), G4CascadeFinalStateAlgorithm::satisfyTriangle(), and G4CascadeFinalStateAlgorithm::SaveKinematics().

{ return verboseLevel; }

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

protectedvirtual

Definition at line 69 of file G4VHadDecayAlgorithm.cc.

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

Referenced by Generate().

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

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

protected

Definition at line 123 of file G4VHadDecayAlgorithm.cc.

References copy().

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

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

virtual void G4VHadDecayAlgorithm::SetVerboseLevel ( G4int  verbose )

inlinevirtual

Reimplemented in G4CascadeFinalStateAlgorithm.

Definition at line 55 of file G4VHadDecayAlgorithm.hh.

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

{ verboseLevel = verbose; }

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

protected

Definition at line 91 of file G4VHadDecayAlgorithm.cc.

References G4cout, G4endl, GetName(), python.hepunit::GeV, and python.hepunit::MeV.

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

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

G4double G4VHadDecayAlgorithm::UniformPhi ( ) const

protected

Definition at line 115 of file G4VHadDecayAlgorithm.cc.

References G4UniformRand, and python.hepunit::twopi.

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

                                                {
  return twopi*G4UniformRand();
}

G4double G4VHadDecayAlgorithm::UniformTheta ( ) const

protected

Definition at line 111 of file G4VHadDecayAlgorithm.cc.

References G4UniformRand.

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

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

---

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

• G4VHadDecayAlgorithm.hh
• G4VHadDecayAlgorithm.cc