Geant4  10.03.p03
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
Public Member Functions | Protected Member Functions | Protected Attributes | List of all members
G4VAdjointReverseReaction Class Reference

#include <G4VAdjointReverseReaction.hh>

Inheritance diagram for G4VAdjointReverseReaction:
Inheritance graph
[legend]
Collaboration diagram for G4VAdjointReverseReaction:
Collaboration graph
[legend]

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.

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  trackid = nstep = 0;
54 }
G4VAdjointReverseReaction::theAdjointCSManager
G4AdjointCSManager * theAdjointCSManager
Definition: G4VAdjointReverseReaction.hh:100
G4AdjointCSManager::GetAdjointCSManager
static G4AdjointCSManager * GetAdjointCSManager()
Definition: G4AdjointCSManager.cc:61

Here is the call graph for this function:

G4VAdjointReverseReaction::~G4VAdjointReverseReaction ( )
virtual

Definition at line 58 of file G4VAdjointReverseReaction.cc.

59 { if (fParticleChange) delete fParticleChange;
60 }

Member Function Documentation

void G4VAdjointReverseReaction::BuildPhysicsTable ( const G4ParticleDefinition &  )
virtual

Reimplemented from G4VProcess.

Definition at line 68 of file G4VAdjointReverseReaction.cc.

69 {
70 
71  theAdjointCSManager->BuildCrossSectionMatrices(); //do not worry it will be done just once
72  theAdjointCSManager->BuildTotalSigmaTables();
73 
74 }
G4VAdjointReverseReaction::theAdjointCSManager
G4AdjointCSManager * theAdjointCSManager
Definition: G4VAdjointReverseReaction.hh:100

Here is the call graph for this function:

G4double G4VAdjointReverseReaction::GetMeanFreePath ( const G4Track &  track,
G4double  previousStepSize,
G4ForceCondition *  condition 
)
protectedvirtual

Implements G4VDiscreteProcess.

Definition at line 106 of file G4VAdjointReverseReaction.cc.

109 { *condition = NotForced;
110  G4double preStepKinEnergy = track.GetKineticEnergy();
111 
112  if(track.GetTrackID() != trackid) {
113  trackid = track.GetTrackID();
114  nstep = 0;
115  }
116  ++nstep;
117 
118 
119 
120  /*G4double Sigma =
121  theAdjointEMModel->AdjointCrossSection(track.GetMaterialCutsCouple(),preStepKinEnergy,IsScatProjToProjCase);*/
122 
123  G4double Sigma =
124  theAdjointEMModel->GetAdjointCrossSection(track.GetMaterialCutsCouple(),preStepKinEnergy,IsScatProjToProjCase);
125 
126  //G4double sig = Sigma;
127 
128  G4double fwd_TotCS;
129  G4double corr = theAdjointCSManager->GetCrossSectionCorrection(track.GetDefinition(),preStepKinEnergy,track.GetMaterialCutsCouple(),IsFwdCSUsed, fwd_TotCS);
130 
131  if(std::fabs(corr) > 100.) { Sigma = 0.0; }
132  else { Sigma *= corr; }
133 
134  //G4cout<<fwd_TotCS<<G4endl;
135  /*if (IsFwdCSUsed && IsIntegralModeUsed){ //take the maximum cross section only for charged particle
136  G4double e_sigma_max, sigma_max;
137  theAdjointCSManager->GetMaxFwdTotalCS(track.GetDefinition(),
138  track.GetMaterialCutsCouple(), e_sigma_max, sigma_max);
139  if (e_sigma_max > preStepKinEnergy){
140  Sigma*=sigma_max/fwd_TotCS;
141  }
142  }
143  */
144 
145  G4double mean_free_path = 1.e60 *mm;
146  if (Sigma>0) mean_free_path = 1./Sigma;
147  lastCS=Sigma;
148  /*
149  if(nstep > 100) {
150 
151  G4cout << "#* " << track.GetDefinition()->GetParticleName()
152  << " " << GetProcessName()
153  << " Nstep " << nstep
154  << " E(MeV)= " << preStepKinEnergy << " Sig0= " << sig
155  << " sig1= " << Sigma << " mfp= " << mean_free_path << G4endl;
156 
157  }
158  if (nstep > 20000) {
159  exit(1);
160  }
161  */
162  /*G4cout<<"Sigma "<<Sigma<<G4endl;
163  G4cout<<"mean_free_path [mm] "<<mean_free_path/mm<<G4endl;
164  */
165 
166 
167  return mean_free_path;
168 }
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:467
mm
static constexpr double mm
Definition: G4SIunits.hh:115
G4VEmAdjointModel::GetAdjointCrossSection
virtual G4double GetAdjointCrossSection(const G4MaterialCutsCouple *aCouple, G4double primEnergy, G4bool IsScatProjToProjCase)
Definition: G4VEmAdjointModel.cc:84
G4Track::GetMaterialCutsCouple
const G4MaterialCutsCouple * GetMaterialCutsCouple() const
G4Track::GetKineticEnergy
G4double GetKineticEnergy() const
G4VAdjointReverseReaction::theAdjointCSManager
G4AdjointCSManager * theAdjointCSManager
Definition: G4VAdjointReverseReaction.hh:100
G4Track::GetTrackID
G4int GetTrackID() const
NotForced
Definition: G4ForceCondition.hh:56
G4double
double G4double
Definition: G4Types.hh:76

Here is the call graph for this function:

G4VParticleChange * G4VAdjointReverseReaction::PostStepDoIt ( const G4Track &  track,
const G4Step &   
)
virtual

Reimplemented from G4VDiscreteProcess.

Definition at line 77 of file G4VAdjointReverseReaction.cc.

78 {
79 
80  fParticleChange->Initialize(track);
81 
82  /* if (IsFwdCSUsed && IsIntegralModeUsed){ //INtegral mode still unstable
83  G4double Tkin = step.GetPostStepPoint()->GetKineticEnergy();
84  G4double fwdCS = theAdjointCSManager->GetTotalForwardCS(track.GetDefinition(), Tkin, track.GetMaterialCutsCouple());
85  //G4cout<<"lastCS "<<lastCS<<G4endl;
86  if (fwdCS<lastCS*G4UniformRand()) { // the reaction does not take place, same integral method as the one used for forward ionisation in G4
87  ClearNumberOfInteractionLengthLeft();
88  return fParticleChange;
89  }
90 
91  }
92  */
93 
94  theAdjointEMModel->SampleSecondaries(track,
95  IsScatProjToProjCase,
96  fParticleChange);
97 
98  ClearNumberOfInteractionLengthLeft();
99  return fParticleChange;
100 
101 
102 
103 }
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

Here is the call graph for this function:

void G4VAdjointReverseReaction::PreparePhysicsTable ( const G4ParticleDefinition &  )
virtual

Reimplemented from G4VProcess.

Definition at line 63 of file G4VAdjointReverseReaction.cc.

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

Definition at line 89 of file G4VAdjointReverseReaction.hh.

89 {IsIntegralModeUsed = aBool;}

Here is the caller graph for this function:

Member Data Documentation

G4ParticleChange* G4VAdjointReverseReaction::fParticleChange
protected

Definition at line 99 of file G4VAdjointReverseReaction.hh.

G4bool G4VAdjointReverseReaction::IsScatProjToProjCase
protected

Definition at line 101 of file G4VAdjointReverseReaction.hh.

G4AdjointCSManager* G4VAdjointReverseReaction::theAdjointCSManager
protected

Definition at line 100 of file G4VAdjointReverseReaction.hh.

G4VEmAdjointModel* G4VAdjointReverseReaction::theAdjointEMModel
protected

Definition at line 98 of file G4VAdjointReverseReaction.hh.

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