Geant4.10
 All Data Structures Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
Public Member Functions | Protected Member Functions | Protected Attributes
G4VAdjointReverseReaction Class Reference

#include <G4VAdjointReverseReaction.hh>

Inheritance diagram for G4VAdjointReverseReaction:
G4VDiscreteProcess G4VProcess G4eInverseBremsstrahlung G4eInverseCompton G4eInverseIonisation G4hInverseIonisation G4InversePEEffect G4IonInverseIonisation

Public Member Functions

 G4VAdjointReverseReaction (G4String process_name, G4bool whichScatCase)
 
virtual ~G4VAdjointReverseReaction ()
 
void PreparePhysicsTable (const G4ParticleDefinition &)
 
void BuildPhysicsTable (const G4ParticleDefinition &)
 
virtual G4VParticleChange * PostStepDoIt (const G4Track &, const G4Step &)
 
void SetIntegralMode (G4bool aBool)
 
- Public Member Functions inherited from G4VDiscreteProcess
 G4VDiscreteProcess (const G4String &, G4ProcessType aType=fNotDefined)
 
 G4VDiscreteProcess (G4VDiscreteProcess &)
 
virtual ~G4VDiscreteProcess ()
 
virtual G4double PostStepGetPhysicalInteractionLength (const G4Track &track, G4double previousStepSize, G4ForceCondition *condition)
 
virtual G4double AlongStepGetPhysicalInteractionLength (const G4Track &, G4double, G4double, G4double &, G4GPILSelection *)
 
virtual G4double AtRestGetPhysicalInteractionLength (const G4Track &, G4ForceCondition *)
 
virtual G4VParticleChange * AtRestDoIt (const G4Track &, const G4Step &)
 
virtual G4VParticleChange * AlongStepDoIt (const G4Track &, const G4Step &)
 
- Public Member Functions inherited from G4VProcess
 G4VProcess (const G4String &aName="NoName", G4ProcessType aType=fNotDefined)
 
 G4VProcess (const G4VProcess &right)
 
virtual ~G4VProcess ()
 
G4int operator== (const G4VProcess &right) const
 
G4int operator!= (const G4VProcess &right) const
 
G4double GetCurrentInteractionLength () const
 
void SetPILfactor (G4double value)
 
G4double GetPILfactor () const
 
G4double AlongStepGPIL (const G4Track &track, G4double previousStepSize, G4double currentMinimumStep, G4double &proposedSafety, G4GPILSelection *selection)
 
G4double AtRestGPIL (const G4Track &track, G4ForceCondition *condition)
 
G4double PostStepGPIL (const G4Track &track, G4double previousStepSize, G4ForceCondition *condition)
 
virtual G4bool IsApplicable (const G4ParticleDefinition &)
 
virtual G4bool StorePhysicsTable (const G4ParticleDefinition *, const G4String &, G4bool)
 
virtual G4bool RetrievePhysicsTable (const G4ParticleDefinition *, const G4String &, G4bool)
 
const G4String & GetPhysicsTableFileName (const G4ParticleDefinition *, const G4String &directory, const G4String &tableName, G4bool ascii=false)
 
const G4String & GetProcessName () const
 
G4ProcessType GetProcessType () const
 
void SetProcessType (G4ProcessType)
 
G4int GetProcessSubType () const
 
void SetProcessSubType (G4int)
 
virtual void StartTracking (G4Track *)
 
virtual void EndTracking ()
 
virtual void SetProcessManager (const G4ProcessManager *)
 
virtual const G4ProcessManager * GetProcessManager ()
 
virtual void ResetNumberOfInteractionLengthLeft ()
 
G4double GetNumberOfInteractionLengthLeft () const
 
G4double GetTotalNumberOfInteractionLengthTraversed () const
 
G4bool isAtRestDoItIsEnabled () const
 
G4bool isAlongStepDoItIsEnabled () const
 
G4bool isPostStepDoItIsEnabled () const
 
virtual void DumpInfo () const
 
void SetVerboseLevel (G4int value)
 
G4int GetVerboseLevel () const
 
virtual void SetMasterProcess (G4VProcess *masterP)
 
const G4VProcess * GetMasterProcess () const
 
virtual void BuildWorkerPhysicsTable (const G4ParticleDefinition &part)
 
virtual void PrepareWorkerPhysicsTable (const G4ParticleDefinition &)
 

Protected Member Functions

virtual G4double GetMeanFreePath (const G4Track &track, G4double previousStepSize, G4ForceCondition *condition)
 
- Protected Member Functions inherited from G4VProcess
void SubtractNumberOfInteractionLengthLeft (G4double previousStepSize)
 
void ClearNumberOfInteractionLengthLeft ()
 

Protected Attributes

G4VEmAdjointModel * theAdjointEMModel
 
G4ParticleChange * fParticleChange
 
G4AdjointCSManager * theAdjointCSManager
 
G4bool IsScatProjToProjCase
 
- Protected Attributes inherited from G4VProcess
const G4ProcessManager * aProcessManager
 
G4VParticleChange * pParticleChange
 
G4ParticleChange aParticleChange
 
G4double theNumberOfInteractionLengthLeft
 
G4double currentInteractionLength
 
G4double theInitialNumberOfInteractionLength
 
G4String theProcessName
 
G4String thePhysicsTableFileName
 
G4ProcessType theProcessType
 
G4int theProcessSubType
 
G4double thePILfactor
 
G4bool enableAtRestDoIt
 
G4bool enableAlongStepDoIt
 
G4bool enablePostStepDoIt
 
G4int verboseLevel
 

Additional Inherited Members

- Static Public Member Functions inherited from G4VProcess
static const G4String & GetProcessTypeName (G4ProcessType)
 

Detailed Description

Definition at line 74 of file G4VAdjointReverseReaction.hh.

Constructor & Destructor Documentation

G4VAdjointReverseReaction::G4VAdjointReverseReaction ( G4String  process_name,
G4bool  whichScatCase 
)

Definition at line 45 of file G4VAdjointReverseReaction.cc.

References fParticleChange, G4AdjointCSManager::GetAdjointCSManager(), IsScatProjToProjCase, and theAdjointCSManager.

45  :
46  G4VDiscreteProcess(process_name)
47 {theAdjointCSManager = G4AdjointCSManager::GetAdjointCSManager();
48  IsScatProjToProjCase=whichScatCase;
49  fParticleChange=new G4ParticleChange();
50  IsFwdCSUsed=false;
51  IsIntegralModeUsed=false;
52  lastCS=0.;
53 }
G4VAdjointReverseReaction::theAdjointCSManager
G4AdjointCSManager * theAdjointCSManager
Definition: G4VAdjointReverseReaction.hh:100
G4AdjointCSManager::GetAdjointCSManager
static G4AdjointCSManager * GetAdjointCSManager()
Definition: G4AdjointCSManager.cc:61
G4VAdjointReverseReaction::~G4VAdjointReverseReaction ( )
virtual

Definition at line 57 of file G4VAdjointReverseReaction.cc.

References fParticleChange.

58 { if (fParticleChange) delete fParticleChange;
59 }

Member Function Documentation

void G4VAdjointReverseReaction::BuildPhysicsTable ( const G4ParticleDefinition &  )
virtual

Reimplemented from G4VProcess.

Definition at line 67 of file G4VAdjointReverseReaction.cc.

References G4AdjointCSManager::BuildCrossSectionMatrices(), G4AdjointCSManager::BuildTotalSigmaTables(), and theAdjointCSManager.

68 {
69 
70  theAdjointCSManager->BuildCrossSectionMatrices(); //do not worry it will be done just once
71  theAdjointCSManager->BuildTotalSigmaTables();
72 
73 }
G4VAdjointReverseReaction::theAdjointCSManager
G4AdjointCSManager * theAdjointCSManager
Definition: G4VAdjointReverseReaction.hh:100
G4double G4VAdjointReverseReaction::GetMeanFreePath ( const G4Track &  track,
G4double  previousStepSize,
G4ForceCondition *  condition 
)
protectedvirtual

Implements G4VDiscreteProcess.

Definition at line 107 of file G4VAdjointReverseReaction.cc.

References G4VEmAdjointModel::GetAdjointCrossSection(), G4AdjointCSManager::GetCrossSectionCorrection(), G4Track::GetDefinition(), G4Track::GetKineticEnergy(), G4Track::GetMaterialCutsCouple(), IsScatProjToProjCase, python.hepunit::mm, NotForced, theAdjointCSManager, and theAdjointEMModel.

110 { *condition = NotForced;
111  G4double preStepKinEnergy = track.GetKineticEnergy();
112 
113  /*G4double Sigma =
114  theAdjointEMModel->AdjointCrossSection(track.GetMaterialCutsCouple(),preStepKinEnergy,IsScatProjToProjCase);*/
115 
116  G4double Sigma =
117  theAdjointEMModel->GetAdjointCrossSection(track.GetMaterialCutsCouple(),preStepKinEnergy,IsScatProjToProjCase);
118  G4double fwd_TotCS;
119  Sigma *= theAdjointCSManager->GetCrossSectionCorrection(track.GetDefinition(),preStepKinEnergy,track.GetMaterialCutsCouple(),IsFwdCSUsed, fwd_TotCS);
120  //G4cout<<fwd_TotCS<<G4endl;
121  /*if (IsFwdCSUsed && IsIntegralModeUsed){ //take the maximum cross section only for charged particle
122  G4double e_sigma_max, sigma_max;
123  theAdjointCSManager->GetMaxFwdTotalCS(track.GetDefinition(),
124  track.GetMaterialCutsCouple(), e_sigma_max, sigma_max);
125  if (e_sigma_max > preStepKinEnergy){
126  Sigma*=sigma_max/fwd_TotCS;
127  }
128  }
129  */
130 
131  G4double mean_free_path = 1.e60 *mm;
132  if (Sigma>0) mean_free_path = 1./Sigma;
133  lastCS=Sigma;
134 
135  /*G4cout<<"Sigma "<<Sigma<<G4endl;
136  G4cout<<"mean_free_path [mm] "<<mean_free_path/mm<<G4endl;
137  */
138 
139 
140  return mean_free_path;
141 }
condition
G4double condition(const G4ErrorSymMatrix &m)
G4Track::GetDefinition
G4ParticleDefinition * GetDefinition() const
G4AdjointCSManager::GetCrossSectionCorrection
G4double GetCrossSectionCorrection(G4ParticleDefinition *aPartDef, G4double PreStepEkin, const G4MaterialCutsCouple *aCouple, G4bool &fwd_is_used, G4double &fwd_TotCS)
Definition: G4AdjointCSManager.cc:465
G4VEmAdjointModel::GetAdjointCrossSection
virtual G4double GetAdjointCrossSection(const G4MaterialCutsCouple *aCouple, G4double primEnergy, G4bool IsScatProjToProjCase)
Definition: G4VEmAdjointModel.cc:83
G4Track::GetMaterialCutsCouple
const G4MaterialCutsCouple * GetMaterialCutsCouple() const
G4Track::GetKineticEnergy
G4double GetKineticEnergy() const
G4VAdjointReverseReaction::theAdjointCSManager
G4AdjointCSManager * theAdjointCSManager
Definition: G4VAdjointReverseReaction.hh:100
NotForced
Definition: G4ForceCondition.hh:56
G4double
double G4double
Definition: G4Types.hh:76
G4VParticleChange * G4VAdjointReverseReaction::PostStepDoIt ( const G4Track &  track,
const G4Step &   
)
virtual

Reimplemented from G4VDiscreteProcess.

Definition at line 76 of file G4VAdjointReverseReaction.cc.

References G4VProcess::ClearNumberOfInteractionLengthLeft(), fParticleChange, G4ParticleChange::Initialize(), IsScatProjToProjCase, G4VEmAdjointModel::SampleSecondaries(), and theAdjointEMModel.

77 {
78 
79 
80 
81  fParticleChange->Initialize(track);
82 
83  /* if (IsFwdCSUsed && IsIntegralModeUsed){ //INtegral mode still unstable
84  G4double Tkin = step.GetPostStepPoint()->GetKineticEnergy();
85  G4double fwdCS = theAdjointCSManager->GetTotalForwardCS(track.GetDefinition(), Tkin, track.GetMaterialCutsCouple());
86  //G4cout<<"lastCS "<<lastCS<<G4endl;
87  if (fwdCS<lastCS*G4UniformRand()) { // the reaction does not take place, same integral method as the one used for forward ionisation in G4
88  ClearNumberOfInteractionLengthLeft();
89  return fParticleChange;
90  }
91 
92  }
93  */
94 
95  theAdjointEMModel->SampleSecondaries(track,
96  IsScatProjToProjCase,
97  fParticleChange);
98 
99  ClearNumberOfInteractionLengthLeft();
100  return fParticleChange;
101 
102 
103 
104 }
G4VProcess::ClearNumberOfInteractionLengthLeft
void ClearNumberOfInteractionLengthLeft()
Definition: G4VProcess.hh:447
G4VEmAdjointModel::SampleSecondaries
virtual void SampleSecondaries(const G4Track &aTrack, G4bool IsScatProjToProjCase, G4ParticleChange *fParticleChange)=0
G4ParticleChange::Initialize
virtual void Initialize(const G4Track &)
Definition: G4ParticleChange.cc:228
void G4VAdjointReverseReaction::PreparePhysicsTable ( const G4ParticleDefinition &  )
virtual

Reimplemented from G4VProcess.

Definition at line 62 of file G4VAdjointReverseReaction.cc.

63 {;
64 }
void G4VAdjointReverseReaction::SetIntegralMode ( G4bool  aBool)
inline

Definition at line 89 of file G4VAdjointReverseReaction.hh.

Referenced by G4eInverseBremsstrahlung::G4eInverseBremsstrahlung(), G4eInverseCompton::G4eInverseCompton(), G4eInverseIonisation::G4eInverseIonisation(), G4hInverseIonisation::G4hInverseIonisation(), G4InversePEEffect::G4InversePEEffect(), and G4IonInverseIonisation::G4IonInverseIonisation().

89 {IsIntegralModeUsed = aBool;}

Field Documentation

G4ParticleChange* G4VAdjointReverseReaction::fParticleChange
protected

Definition at line 99 of file G4VAdjointReverseReaction.hh.

Referenced by G4VAdjointReverseReaction(), PostStepDoIt(), and ~G4VAdjointReverseReaction().

G4bool G4VAdjointReverseReaction::IsScatProjToProjCase
protected

Definition at line 101 of file G4VAdjointReverseReaction.hh.

Referenced by G4eInverseBremsstrahlung::G4eInverseBremsstrahlung(), G4VAdjointReverseReaction(), GetMeanFreePath(), and PostStepDoIt().

G4AdjointCSManager* G4VAdjointReverseReaction::theAdjointCSManager
protected

Definition at line 100 of file G4VAdjointReverseReaction.hh.

Referenced by BuildPhysicsTable(), G4VAdjointReverseReaction(), and GetMeanFreePath().

G4VEmAdjointModel* G4VAdjointReverseReaction::theAdjointEMModel
protected

Definition at line 98 of file G4VAdjointReverseReaction.hh.

Referenced by G4eInverseBremsstrahlung::G4eInverseBremsstrahlung(), G4eInverseCompton::G4eInverseCompton(), G4eInverseIonisation::G4eInverseIonisation(), G4hInverseIonisation::G4hInverseIonisation(), G4InversePEEffect::G4InversePEEffect(), G4IonInverseIonisation::G4IonInverseIonisation(), GetMeanFreePath(), and PostStepDoIt().

The documentation for this class was generated from the following files:
Generated on Wed Apr 30 2014 15:57:10 for Geant4.10 by   doxygen 1.8.7