Geant4  10.03.p03
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
Public Member Functions | List of all members
CexmcChargeExchangeProductionModel< OutputParticle > Class Template Reference

#include <CexmcChargeExchangeProductionModel.hh>

Inheritance diagram for CexmcChargeExchangeProductionModel< OutputParticle >:
Inheritance graph
[legend]
Collaboration diagram for CexmcChargeExchangeProductionModel< OutputParticle >:
Collaboration graph
[legend]

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 &aTrack, G4Nucleus &targetNucleus)
 
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)
 
G4int GetVerboseLevel () const
 
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
 
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
 
virtual void BuildPhysicsTable (const G4ParticleDefinition &)
 
virtual void InitialiseModel ()
 
- 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.

81  :
82  G4HadronicInteraction( CexmcChargeExchangeInteractionName ),
83  CexmcProductionModel( CexmcChargeExchangeProductionModelName ),
84  nucleusParticleMass( 0 ), phaseSpaceGenerator( NULL )
85 {
86  incidentParticle = G4PionMinus::Definition();
87  nucleusParticle = G4Proton::Definition();
88  outputParticle = OutputParticle::Definition();
89  nucleusOutputParticle = G4Neutron::Definition();
90 
91  nucleusParticleMass = nucleusParticle->GetPDGMass();
92 
93  productionModelData.incidentParticle = incidentParticle;
94  productionModelData.nucleusParticle = nucleusParticle;
95  productionModelData.outputParticle = outputParticle;
96  productionModelData.nucleusOutputParticle = nucleusOutputParticle;
97 
98  CexmcPhaseSpaceInVector inVec;
99 
100  inVec.push_back( &productionModelData.incidentParticleSCM );
101  inVec.push_back( &productionModelData.nucleusParticleSCM );
102 
103  CexmcPhaseSpaceOutVector outVec;
104 
105  outVec.push_back( CexmcPhaseSpaceOutVectorElement(
106  &productionModelData.outputParticleSCM,
107  outputParticle->GetPDGMass() ) );
108  outVec.push_back( CexmcPhaseSpaceOutVectorElement(
109  &productionModelData.nucleusOutputParticleSCM,
110  nucleusOutputParticle->GetPDGMass() ) );
111 
112 #ifdef CEXMC_USE_GENBOD
113  phaseSpaceGenerator = new CexmcGenbod;
114 #else
115  phaseSpaceGenerator = new CexmcReimplementedGenbod;
116 #endif
117 
118  phaseSpaceGenerator->SetParticles( inVec, outVec );
119 }
CexmcChargeExchangeInteractionName
const G4String CexmcChargeExchangeInteractionName("Cexmc"+CexmcChargeExchangeProductionModelName)
CexmcPhaseSpaceInVector
std::vector< const G4LorentzVector * > CexmcPhaseSpaceInVector
Definition: CexmcPhaseSpaceGenerator.hh:64
CexmcProductionModelData::nucleusParticle
const G4ParticleDefinition * nucleusParticle
Definition: CexmcProductionModelData.hh:97
CexmcProductionModel::nucleusOutputParticle
G4ParticleDefinition * nucleusOutputParticle
Definition: CexmcProductionModel.hh:135
CexmcPhaseSpaceOutVector
std::vector< CexmcPhaseSpaceOutVectorElement > CexmcPhaseSpaceOutVector
Definition: CexmcPhaseSpaceGenerator.hh:67
CexmcProductionModelData::incidentParticle
const G4ParticleDefinition * incidentParticle
Definition: CexmcProductionModelData.hh:95
CexmcProductionModel::productionModelData
CexmcProductionModelData productionModelData
Definition: CexmcProductionModel.hh:126
CexmcProductionModelData::outputParticle
const G4ParticleDefinition * outputParticle
Definition: CexmcProductionModelData.hh:99
G4Proton::Definition
static G4Proton * Definition()
Definition: G4Proton.cc:49
CexmcProductionModel::CexmcProductionModel
CexmcProductionModel(const G4String &name="unspecified", G4bool fermiMotionIsOn=false)
Definition: CexmcProductionModel.cc:49
CexmcProductionModel::nucleusParticle
G4ParticleDefinition * nucleusParticle
Definition: CexmcProductionModel.hh:131
G4HadronicInteraction::G4HadronicInteraction
G4HadronicInteraction(const G4String &modelName="HadronicModel")
Definition: G4HadronicInteraction.cc:42
G4ParticleDefinition::GetPDGMass
G4double GetPDGMass() const
Definition: G4ParticleDefinition.hh:122
CexmcProductionModel::incidentParticle
G4ParticleDefinition * incidentParticle
Definition: CexmcProductionModel.hh:129
G4PionMinus::Definition
static G4PionMinus * Definition()
Definition: G4PionMinus.cc:52
G4Neutron::Definition
static G4Neutron * Definition()
Definition: G4Neutron.cc:54
CexmcProductionModelData::nucleusOutputParticle
const G4ParticleDefinition * nucleusOutputParticle
Definition: CexmcProductionModelData.hh:101
CexmcPhaseSpaceGenerator::SetParticles
void SetParticles(const CexmcPhaseSpaceInVector &inVec_, const CexmcPhaseSpaceOutVector &outVec_)
Definition: CexmcPhaseSpaceGenerator.cc:58
CexmcChargeExchangeProductionModelName
const G4String CexmcChargeExchangeProductionModelName("ChargeExchange")
CexmcProductionModel::outputParticle
G4ParticleDefinition * outputParticle
Definition: CexmcProductionModel.hh:133

Here is the call graph for this function:

template<typename OutputParticle >
CexmcChargeExchangeProductionModel< OutputParticle >::~CexmcChargeExchangeProductionModel ( )

Definition at line 124 of file CexmcChargeExchangeProductionModel.hh.

125 {
126  delete phaseSpaceGenerator;
127 }

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.

134 {
135  theParticleChange.Clear();
136 
137  G4double kinEnergy( projectile.GetKineticEnergy() );
138  G4HadProjectile & theProjectile( const_cast< G4HadProjectile & >(
139  projectile ) );
140  const G4LorentzRotation & projToLab(
141  const_cast< const G4LorentzRotation & >(
142  theProjectile.GetTrafoToLab() ) );
143  productionModelData.incidentParticleLAB = projectile.Get4Momentum();
144  productionModelData.incidentParticleLAB.transform( projToLab );
145  productionModelData.nucleusParticleLAB.setPx( 0 );
146  productionModelData.nucleusParticleLAB.setPy( 0 );
147  productionModelData.nucleusParticleLAB.setPz( 0 );
148  productionModelData.nucleusParticleLAB.setE( nucleusParticleMass );
149 
150  if ( fermiMotionIsOn )
151  {
152  G4ThreeVector targetNucleusMomentum(
153  targetNucleus.GetFermiMomentum() );
154  G4double targetNucleusEnergy(
155  std::sqrt( targetNucleusMomentum.mag2() +
156  nucleusParticleMass * nucleusParticleMass ) );
157  productionModelData.nucleusParticleLAB = G4LorentzVector(
158  targetNucleusMomentum, targetNucleusEnergy );
159  }
160  productionModelData.nucleusParticleLAB.transform( projToLab );
161  G4LorentzVector lVecSum( productionModelData.incidentParticleLAB +
162  productionModelData.nucleusParticleLAB );
163  G4ThreeVector boostVec( lVecSum.boostVector() );
164 
165  productionModelData.incidentParticleSCM =
166  productionModelData.incidentParticleLAB;
167  productionModelData.nucleusParticleSCM =
168  productionModelData.nucleusParticleLAB;
169 
170  productionModelData.incidentParticleSCM.boost( -boostVec );
171  productionModelData.nucleusParticleSCM.boost( -boostVec );
172 
173  triggeredAngularRanges.clear();
174 
175  if ( ! phaseSpaceGenerator->CheckKinematics() )
176  {
177  theParticleChange.SetEnergyChange( kinEnergy );
178  theParticleChange.SetMomentumChange(
179  projectile.Get4Momentum().vect().unit());
180  return &theParticleChange;
181  }
182 
183  do
184  {
185  phaseSpaceGenerator->Generate();
186  G4double cosTheta( productionModelData.outputParticleSCM.cosTheta() );
187  for ( CexmcAngularRangeList::iterator k( angularRanges.begin() );
188  k != angularRanges.end(); ++k )
189  {
190  if ( cosTheta <= k->top && cosTheta > k->bottom )
191  triggeredAngularRanges.push_back( CexmcAngularRange(
192  k->top, k->bottom, k->index ) );
193  }
194  } while ( triggeredAngularRanges.empty() );
195 
196  productionModelData.outputParticleLAB =
197  productionModelData.outputParticleSCM;
198  productionModelData.nucleusOutputParticleLAB =
199  productionModelData.nucleusOutputParticleSCM;
200 
201  productionModelData.outputParticleLAB.boost( boostVec );
202  productionModelData.nucleusOutputParticleLAB.boost( boostVec );
203 
204  theParticleChange.SetStatusChange( stopAndKill );
205  theParticleChange.SetEnergyChange( 0.0 );
206 
207  G4DynamicParticle * secOutParticle( new G4DynamicParticle(
208  outputParticle,
209  productionModelData.outputParticleLAB ) );
210  theParticleChange.AddSecondary( secOutParticle );
211  G4DynamicParticle * secNeutron( new G4DynamicParticle(
212  nucleusOutputParticle,
213  productionModelData.nucleusOutputParticleLAB ) );
214  theParticleChange.AddSecondary( secNeutron );
215 
216  /* projectile->GetDefinition() shall always be identical to incidentParticle
217  * as far as CexmcHadronicProcess::IsApplicable() will check that only
218  * incidentParticle is allowed. Here is mostly unnecessary assignment. */
219  productionModelData.incidentParticle = projectile.GetDefinition();
220 
221  return &theParticleChange;
222 }
CexmcProductionModel::nucleusOutputParticle
G4ParticleDefinition * nucleusOutputParticle
Definition: CexmcProductionModel.hh:135
stopAndKill
Definition: G4HadFinalState.hh:42
CexmcProductionModelData::incidentParticle
const G4ParticleDefinition * incidentParticle
Definition: CexmcProductionModelData.hh:95
CexmcProductionModel::productionModelData
CexmcProductionModelData productionModelData
Definition: CexmcProductionModel.hh:126
G4HadFinalState::SetStatusChange
void SetStatusChange(G4HadFinalStateStatus aS)
Definition: G4HadFinalState.hh:67
CLHEP::HepLorentzVector::vect
Hep3Vector vect() const
G4HadProjectile::GetDefinition
const G4ParticleDefinition * GetDefinition() const
Definition: G4HadProjectile.hh:81
CLHEP::HepLorentzVector::boost
HepLorentzVector & boost(double, double, double)
Definition: LorentzVector.cc:59
G4HadProjectile::GetKineticEnergy
G4double GetKineticEnergy() const
Definition: G4HadProjectile.hh:106
G4HadProjectile::Get4Momentum
const G4LorentzVector & Get4Momentum() const
Definition: G4HadProjectile.hh:86
G4HadFinalState::SetEnergyChange
void SetEnergyChange(G4double anEnergy)
Definition: G4HadFinalState.cc:39
G4Nucleus::GetFermiMomentum
G4ThreeVector GetFermiMomentum()
Definition: G4Nucleus.cc:412
CLHEP::Hep3Vector::unit
Hep3Vector unit() const
CexmcProductionModel::angularRanges
CexmcAngularRangeList angularRanges
Definition: CexmcProductionModel.hh:120
CexmcProductionModel::triggeredAngularRanges
CexmcAngularRangeList triggeredAngularRanges
Definition: CexmcProductionModel.hh:124
CLHEP::HepLorentzVector::transform
HepLorentzVector & transform(const HepRotation &)
G4HadFinalState::AddSecondary
void AddSecondary(G4DynamicParticle *aP, G4int mod=-1)
Definition: G4HadFinalState.hh:61
CexmcPhaseSpaceGenerator::CheckKinematics
virtual G4bool CheckKinematics(void)
Definition: CexmcPhaseSpaceGenerator.cc:82
G4double
double G4double
Definition: G4Types.hh:76
G4HadFinalState::SetMomentumChange
void SetMomentumChange(const G4ThreeVector &aV)
Definition: G4HadFinalState.hh:56
CexmcPhaseSpaceGenerator::Generate
virtual G4double Generate(void)=0
CexmcProductionModel::outputParticle
G4ParticleDefinition * outputParticle
Definition: CexmcProductionModel.hh:133
CLHEP::HepLorentzVector::cosTheta
double cosTheta() const
G4LorentzVector
CLHEP::HepLorentzVector G4LorentzVector
Definition: G4LorentzVector.hh:35

Here is the call graph for this function:

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