CexmcChargeExchangeProductionModel< OutputParticle > Class Template Reference

#include <CexmcChargeExchangeProductionModel.hh>

Inheritance diagram for CexmcChargeExchangeProductionModel< OutputParticle >:

Public Member Functions
	CexmcChargeExchangeProductionModel ()
	~CexmcChargeExchangeProductionModel ()
G4HadFinalState *	ApplyYourself (const G4HadProjectile &projectile, G4Nucleus &targetNucleus)
Public Member Functions inherited from G4HadronicInteraction
	G4HadronicInteraction (const G4String &modelName="HadronicModel")
virtual	~G4HadronicInteraction ()
virtual G4double	SampleInvariantT (const G4ParticleDefinition *p, G4double plab, G4int Z, G4int A)
virtual G4bool	IsApplicable (const G4HadProjectile &, G4Nucleus &)
G4double	GetMinEnergy () const
G4double	GetMinEnergy (const G4Material *aMaterial, const G4Element *anElement) const
void	SetMinEnergy (G4double anEnergy)
void	SetMinEnergy (G4double anEnergy, const G4Element *anElement)
void	SetMinEnergy (G4double anEnergy, const G4Material *aMaterial)
G4double	GetMaxEnergy () const
G4double	GetMaxEnergy (const G4Material *aMaterial, const G4Element *anElement) const
void	SetMaxEnergy (const G4double anEnergy)
void	SetMaxEnergy (G4double anEnergy, const G4Element *anElement)
void	SetMaxEnergy (G4double anEnergy, const G4Material *aMaterial)
const G4HadronicInteraction *	GetMyPointer () const
virtual G4int	GetVerboseLevel () const
virtual void	SetVerboseLevel (G4int value)
const G4String &	GetModelName () const
void	DeActivateFor (const G4Material *aMaterial)
void	ActivateFor (const G4Material *aMaterial)
void	DeActivateFor (const G4Element *anElement)
void	ActivateFor (const G4Element *anElement)
G4bool	IsBlocked (const G4Material *aMaterial) const
G4bool	IsBlocked (const G4Element *anElement) const
void	SetRecoilEnergyThreshold (G4double val)
G4double	GetRecoilEnergyThreshold () const
G4bool	operator== (const G4HadronicInteraction &right) const
G4bool	operator!= (const G4HadronicInteraction &right) const
virtual const std::pair < G4double, G4double >	GetFatalEnergyCheckLevels () const
virtual std::pair< G4double, G4double >	GetEnergyMomentumCheckLevels () const
void	SetEnergyMomentumCheckLevels (G4double relativeLevel, G4double absoluteLevel)
virtual void	ModelDescription (std::ostream &outFile) const
Public Member Functions inherited from CexmcProductionModel
	CexmcProductionModel (const G4String &name="unspecified", G4bool fermiMotionIsOn=false)
virtual	~CexmcProductionModel ()
void	ApplyFermiMotion (G4bool on, G4bool fromMessenger=true)
void	SetAngularRange (G4double top, G4double bottom, G4int nmbOfDivs)
void	SetAngularRanges (const CexmcAngularRangeList &angularRanges_)
void	AddAngularRange (G4double top, G4double bottom, G4int nmbOfDivs)
void	SetProductionModelData (const CexmcProductionModelData &productionModelData_)
void	PrintInitialData (void) const
const CexmcAngularRangeList &	GetAngularRanges (void) const
const CexmcAngularRangeList &	GetTriggeredAngularRanges (void) const
const CexmcProductionModelData &	GetProductionModelData (void) const
G4bool	IsFermiMotionOn (void) const
void	SetTriggeredAngularRanges (G4double opCosThetaSCM)
const G4String &	GetName (void) const
G4ParticleDefinition *	GetIncidentParticle (void) const
G4ParticleDefinition *	GetNucleusParticle (void) const
G4ParticleDefinition *	GetOutputParticle (void) const
G4ParticleDefinition *	GetNucleusOutputParticle (void) const

Additional Inherited Members
Protected Member Functions inherited from G4HadronicInteraction
void	SetModelName (const G4String &nam)
G4bool	IsBlocked () const
void	Block ()
Protected Member Functions inherited from CexmcProductionModel
virtual void	FermiMotionStatusChangeHook (void)
Protected Attributes inherited from G4HadronicInteraction
G4HadFinalState	theParticleChange
G4int	verboseLevel
G4double	theMinEnergy
G4double	theMaxEnergy
G4bool	isBlocked
Protected Attributes inherited from CexmcProductionModel
G4String	name
G4bool	fermiMotionIsOn
CexmcAngularRangeList	angularRanges
CexmcAngularRangeList	angularRangesRef
CexmcAngularRangeList	triggeredAngularRanges
CexmcProductionModelData	productionModelData
G4ParticleDefinition *	incidentParticle
G4ParticleDefinition *	nucleusParticle
G4ParticleDefinition *	outputParticle
G4ParticleDefinition *	nucleusOutputParticle

Detailed Description

template<typename OutputParticle>
class CexmcChargeExchangeProductionModel< OutputParticle >

Definition at line 60 of file CexmcChargeExchangeProductionModel.hh.

Constructor & Destructor Documentation

template<typename OutputParticle >

CexmcChargeExchangeProductionModel< OutputParticle >::CexmcChargeExchangeProductionModel

(

)

Definition at line 81 of file CexmcChargeExchangeProductionModel.hh.

References G4PionMinus::Definition(), G4Proton::Definition(), G4Neutron::Definition(), G4ParticleDefinition::GetPDGMass(), CexmcProductionModelData::incidentParticle, CexmcProductionModel::incidentParticle, CexmcProductionModelData::incidentParticleSCM, CexmcProductionModelData::nucleusOutputParticle, CexmcProductionModel::nucleusOutputParticle, CexmcProductionModelData::nucleusOutputParticleSCM, CexmcProductionModelData::nucleusParticle, CexmcProductionModel::nucleusParticle, CexmcProductionModelData::nucleusParticleSCM, CexmcProductionModelData::outputParticle, CexmcProductionModel::outputParticle, CexmcProductionModelData::outputParticleSCM, CexmcProductionModel::productionModelData, and CexmcPhaseSpaceGenerator::SetParticles().

                                                                             :
    G4HadronicInteraction( CexmcChargeExchangeInteractionName ),
    CexmcProductionModel( CexmcChargeExchangeProductionModelName ),
    nucleusParticleMass( 0 ), phaseSpaceGenerator( NULL )
{
    incidentParticle = G4PionMinus::Definition();
    nucleusParticle = G4Proton::Definition();
    outputParticle = OutputParticle::Definition();
    nucleusOutputParticle = G4Neutron::Definition();

    nucleusParticleMass = nucleusParticle->GetPDGMass();

    productionModelData.incidentParticle = incidentParticle;
    productionModelData.nucleusParticle = nucleusParticle;
    productionModelData.outputParticle = outputParticle;
    productionModelData.nucleusOutputParticle = nucleusOutputParticle;

    CexmcPhaseSpaceInVector   inVec;

    inVec.push_back( &productionModelData.incidentParticleSCM );
    inVec.push_back( &productionModelData.nucleusParticleSCM );

    CexmcPhaseSpaceOutVector  outVec;

    outVec.push_back( CexmcPhaseSpaceOutVectorElement(
                            &productionModelData.outputParticleSCM,
                            outputParticle->GetPDGMass() ) );
    outVec.push_back( CexmcPhaseSpaceOutVectorElement(
                            &productionModelData.nucleusOutputParticleSCM,
                            nucleusOutputParticle->GetPDGMass() ) );

#ifdef CEXMC_USE_GENBOD
    phaseSpaceGenerator = new CexmcGenbod;
#else
    phaseSpaceGenerator = new CexmcReimplementedGenbod;
#endif

    phaseSpaceGenerator->SetParticles( inVec, outVec );
}

template<typename OutputParticle >

CexmcChargeExchangeProductionModel< OutputParticle >::~CexmcChargeExchangeProductionModel

(

)

Definition at line 124 of file CexmcChargeExchangeProductionModel.hh.

{
    delete phaseSpaceGenerator;
}

Member Function Documentation

template<typename OutputParticle >

G4HadFinalState * CexmcChargeExchangeProductionModel< OutputParticle >::ApplyYourself ( const G4HadProjectile &  projectile, G4Nucleus &  targetNucleus  )

virtual

Implements G4HadronicInteraction.

Definition at line 132 of file CexmcChargeExchangeProductionModel.hh.

References CLHEP::HepLorentzVector::boostVector(), G4HadFinalState::Clear(), G4HadProjectile::Get4Momentum(), G4HadProjectile::GetDefinition(), G4Nucleus::GetFermiMomentum(), G4HadProjectile::GetKineticEnergy(), CLHEP::HepLorentzVector::setPx(), stopAndKill, CLHEP::HepLorentzVector::transform(), CLHEP::Hep3Vector::unit(), and CLHEP::HepLorentzVector::vect().

{
    theParticleChange.Clear();

    G4double           kinEnergy( projectile.GetKineticEnergy() );
    G4HadProjectile &  theProjectile( const_cast< G4HadProjectile & >(
                                                                projectile ) );
    const G4LorentzRotation &  projToLab(
                                    const_cast< const G4LorentzRotation & >(
                                            theProjectile.GetTrafoToLab() ) );
    productionModelData.incidentParticleLAB = projectile.Get4Momentum();
    productionModelData.incidentParticleLAB.transform( projToLab );
    productionModelData.nucleusParticleLAB.setPx( 0 );
    productionModelData.nucleusParticleLAB.setPy( 0 );
    productionModelData.nucleusParticleLAB.setPz( 0 );
    productionModelData.nucleusParticleLAB.setE( nucleusParticleMass );

    if ( fermiMotionIsOn )
    {
        G4ThreeVector  targetNucleusMomentum(
                                        targetNucleus.GetFermiMomentum() );
        G4double       targetNucleusEnergy(
                            std::sqrt( targetNucleusMomentum.mag2() +
                                nucleusParticleMass * nucleusParticleMass ) );
        productionModelData.nucleusParticleLAB = G4LorentzVector(
                                targetNucleusMomentum, targetNucleusEnergy );
    }
    productionModelData.nucleusParticleLAB.transform( projToLab );
    G4LorentzVector  lVecSum( productionModelData.incidentParticleLAB +
                              productionModelData.nucleusParticleLAB );
    G4ThreeVector    boostVec( lVecSum.boostVector() );

    productionModelData.incidentParticleSCM =
            productionModelData.incidentParticleLAB;
    productionModelData.nucleusParticleSCM =
            productionModelData.nucleusParticleLAB;

    productionModelData.incidentParticleSCM.boost( -boostVec );
    productionModelData.nucleusParticleSCM.boost( -boostVec );

    triggeredAngularRanges.clear();

    if ( ! phaseSpaceGenerator->CheckKinematics() )
    {
        theParticleChange.SetEnergyChange( kinEnergy );
        theParticleChange.SetMomentumChange(
                                    projectile.Get4Momentum().vect().unit());
        return &theParticleChange;
    }

    do
    {
        phaseSpaceGenerator->Generate();
        G4double  cosTheta( productionModelData.outputParticleSCM.cosTheta() );
        for ( CexmcAngularRangeList::iterator  k( angularRanges.begin() );
                                                k != angularRanges.end(); ++k )
        {
            if ( cosTheta <= k->top && cosTheta > k->bottom )
                triggeredAngularRanges.push_back( CexmcAngularRange(
                                                k->top, k->bottom, k->index ) );
        }
    } while ( triggeredAngularRanges.empty() );

    productionModelData.outputParticleLAB =
            productionModelData.outputParticleSCM;
    productionModelData.nucleusOutputParticleLAB =
            productionModelData.nucleusOutputParticleSCM;

    productionModelData.outputParticleLAB.boost( boostVec );
    productionModelData.nucleusOutputParticleLAB.boost( boostVec );

    theParticleChange.SetStatusChange( stopAndKill );
    theParticleChange.SetEnergyChange( 0.0 );

    G4DynamicParticle *  secOutParticle( new G4DynamicParticle(
                            outputParticle,
                            productionModelData.outputParticleLAB ) );
    theParticleChange.AddSecondary( secOutParticle );
    G4DynamicParticle *  secNeutron( new G4DynamicParticle(
                            nucleusOutputParticle,
                            productionModelData.nucleusOutputParticleLAB ) );
    theParticleChange.AddSecondary( secNeutron );

    /* projectile->GetDefinition() shall always be identical to incidentParticle
     * as far as CexmcHadronicProcess::IsApplicable() will check that only
     * incidentParticle is allowed. Here is mostly unnecessary assignment. */
    productionModelData.incidentParticle = projectile.GetDefinition();

    return &theParticleChange;
}

---

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

• CexmcChargeExchangeProductionModel.hh