Changeset 1196 for trunk/source/processes/electromagnetic/adjoint/include/G4AdjointAlongStepWeightCorrection.hh
- Timestamp:
- Nov 25, 2009, 5:13:58 PM (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/source/processes/electromagnetic/adjoint/include/G4AdjointAlongStepWeightCorrection.hh
r966 r1196 24 24 // ******************************************************************** 25 25 // 26 // $Id: G4AdjointAlongStepWeightCorrection.hh,v 1.4 2009/11/20 10:31:20 ldesorgh Exp $ 27 // GEANT4 tag $Name: geant4-09-03-cand-01 $ 28 // 26 29 ///////////////////////////////////////////////////////////////////////////////// 27 // Module: G4AdjointAlongStepWeightCorrection.hh30 // Class: G4AdjointAlongStepWeightCorrection 28 31 // Author: L. Desorgher 29 // Date: 10 May 200730 32 // Organisation: SpaceIT GmbH 33 // Contract: ESA contract 21435/08/NL/AT 31 34 // Customer: ESA/ESTEC 32 35 ///////////////////////////////////////////////////////////////////////////////// … … 35 38 // -------------- 36 39 // ChangeHistory: 37 // 10 May 2007 creation by L. Desorgher 40 // 10 May 2007 creation by L. Desorgher 41 // October 2009 implementation of the mode where the total adjoint and forward cross sections are equivalent. L. Desorgher 38 42 // 39 43 //------------------------------------------------------------- 40 44 // Documentation: 41 45 // Continuous processes acting on adjoint particles to correct continuously their weight during the adjoint reverse tracking. 46 // Thi process is needed whene the adjoint cross section are not scaled such that the total adjoint cross section match the total forward cross section. 47 // By default the mode where the total adjoint cross section is equal to the total forward cross section is used an therefore this along step weight 48 // correction factor is 1. 49 // However in some cases (some energy ranges) the total forward cross section or the total adjoint cross section can be null, in this case the along step 50 // weight correction is neede and is given by exp(-(Sigma_tot_adj-Sigma_tot_fwd).dx) 51 // 52 // 42 53 // 43 54 … … 120 131 currentMaterial = couple->GetMaterial(); 121 132 currentMaterialIndex = couple->GetIndex(); 122 //G4cout<<"Define Material"<<std::endl; 123 //if(!meanFreePath) ResetNumberOfInteractionLengthLeft(); 133 124 134 } 125 135 } 126 136 127 137 128 /////////////////////////////////////////////////////// 129 // 130 inline G4double G4AdjointAlongStepWeightCorrection::GetContinuousStepLimit(const G4Track& track, 131 G4double , G4double , G4double& ) 132 { 133 G4double x = DBL_MAX; 134 DefineMaterial(track.GetMaterialCutsCouple()); 135 preStepKinEnergy = track.GetKineticEnergy(); 136 return x; 137 } 138 138 139 #endif
Note: See TracChangeset
for help on using the changeset viewer.