Geant4 Cross Reference |
1 // 2 // ******************************************************************** 3 // * License and Disclaimer * 4 // * * 5 // * The Geant4 software is copyright of the Copyright Holders of * 6 // * the Geant4 Collaboration. It is provided under the terms and * 7 // * conditions of the Geant4 Software License, included in the file * 8 // * LICENSE and available at http://cern.ch/geant4/license . These * 9 // * include a list of copyright holders. * 10 // * * 11 // * Neither the authors of this software system, nor their employing * 12 // * institutes,nor the agencies providing financial support for this * 13 // * work make any representation or warranty, express or implied, * 14 // * regarding this software system or assume any liability for its * 15 // * use. Please see the license in the file LICENSE and URL above * 16 // * for the full disclaimer and the limitation of liability. * 17 // * * 18 // * This code implementation is the result of the scientific and * 19 // * technical work of the GEANT4 collaboration. * 20 // * By using, copying, modifying or distributing the software (or * 21 // * any work based on the software) you agree to acknowledge its * 22 // * use in resulting scientific publications, and indicate your * 23 // * acceptance of all terms of the Geant4 Software license. * 24 // ******************************************************************** 25 // 26 27 #include "G4VAdjointReverseReaction.hh" 28 29 #include "G4AdjointCSManager.hh" 30 #include "G4ParticleChange.hh" 31 #include "G4VEmAdjointModel.hh" 32 33 G4VAdjointReverseReaction::G4VAdjointReverseReaction(const G4String& process_name, 34 G4bool whichScatCase) 35 : G4VDiscreteProcess(process_name) 36 { 37 fCSManager = G4AdjointCSManager::GetAdjointCSManager(); 38 fIsScatProjToProj = whichScatCase; 39 fParticleChange = new G4ParticleChange(); 40 } 41 42 ////////////////////////////////////////////////////////////////////////////// 43 G4VAdjointReverseReaction::~G4VAdjointReverseReaction() 44 { 45 if(fParticleChange) 46 delete fParticleChange; 47 } 48 49 ////////////////////////////////////////////////////////////////////////////// 50 void G4VAdjointReverseReaction::BuildPhysicsTable(const G4ParticleDefinition&) 51 { 52 fCSManager->BuildCrossSectionMatrices(); // it will be done just once 53 fCSManager->BuildTotalSigmaTables(); 54 } 55 56 ////////////////////////////////////////////////////////////////////////////// 57 G4VParticleChange* G4VAdjointReverseReaction::PostStepDoIt(const G4Track& track, 58 const G4Step&) 59 { 60 fParticleChange->Initialize(track); 61 fAdjointModel->SampleSecondaries(track, fIsScatProjToProj, fParticleChange); 62 63 ClearNumberOfInteractionLengthLeft(); 64 return fParticleChange; 65 } 66 67 ////////////////////////////////////////////////////////////////////////////// 68 G4double G4VAdjointReverseReaction::GetMeanFreePath(const G4Track& track, 69 G4double, 70 G4ForceCondition* condition) 71 { 72 *condition = NotForced; 73 G4double preStepKinEnergy = track.GetKineticEnergy(); 74 75 if(track.GetTrackID() != fTrackId) 76 { 77 fTrackId = track.GetTrackID(); 78 } 79 G4double sigma = fAdjointModel->AdjointCrossSection( 80 track.GetMaterialCutsCouple(), preStepKinEnergy, fIsScatProjToProj); 81 82 G4double corr = fCSManager->GetCrossSectionCorrection( 83 track.GetDefinition(), preStepKinEnergy, track.GetMaterialCutsCouple(), 84 fIsFwdCSUsed); 85 86 if(std::fabs(corr) > 100.) 87 { 88 sigma = 0.0; 89 } 90 else 91 { 92 sigma *= corr; 93 } 94 95 G4double mean_free_path = 1.e60; 96 if(sigma > 0.) 97 mean_free_path = 1. / sigma; 98 99 return mean_free_path; 100 } 101