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 // G4DNADoubleIonisation.cc 27 // 28 // Created at 2024/04/03 (Thu.) 29 // Author: Shogo OKADA @KEK-CRC (shogo.okada@kek.jp) 30 // 31 32 #include "G4DNADoubleIonisation.hh" 33 #include "G4DNADoubleIonisationModel.hh" 34 #include "G4LowEnergyEmProcessSubType.hh" 35 #include "G4SystemOfUnits.hh" 36 #include "G4GenericIon.hh" 37 #include "G4Positron.hh" 38 39 //------------------------------------------------------------------------------ 40 G4DNADoubleIonisation::G4DNADoubleIonisation( 41 const G4String& pname, G4ProcessType type) 42 : G4VEmProcess(pname, type), 43 is_initialized_(false) 44 { 45 SetProcessSubType(fLowEnergyDoubleIonisation); 46 } 47 48 //------------------------------------------------------------------------------ 49 G4bool G4DNADoubleIonisation::IsApplicable(const G4ParticleDefinition& p) 50 { 51 return ( 52 &p == G4Proton::Proton() || 53 &p == G4DNAGenericIonsManager::Instance()->GetIon("alpha++") || 54 &p == G4GenericIon::GenericIonDefinition() 55 ); 56 } 57 58 //------------------------------------------------------------------------------ 59 void G4DNADoubleIonisation::InitialiseProcess(const G4ParticleDefinition* p) 60 { 61 if (is_initialized_) { return; } 62 63 is_initialized_ = true; 64 SetBuildTableFlag(false); 65 66 const auto& name = p->GetParticleName(); 67 68 if (name == "proton") { 69 70 if (!EmModel()) { 71 auto ptr = new G4DNADoubleIonisationModel(); 72 SetEmModel(ptr); 73 ptr->SetLowEnergyLimit(0.0 * keV); 74 ptr->SetHighEnergyLimit(3.0 * MeV); 75 } 76 77 AddEmModel(1, EmModel()); 78 79 } else if (name == "alpha") { 80 81 if (!EmModel()) { 82 auto ptr = new G4DNADoubleIonisationModel(); 83 SetEmModel(ptr); 84 ptr->SetLowEnergyLimit(0.0 * keV); 85 ptr->SetHighEnergyLimit(23.0 * MeV); 86 } 87 88 AddEmModel(1, EmModel()); 89 90 } else if (name == "GenericIon") { 91 92 // for carbon ions (12C6+) 93 if (!EmModel()) { 94 auto ptr = new G4DNADoubleIonisationModel(); 95 SetEmModel(ptr); 96 ptr->SetLowEnergyLimit(0.0 * keV); 97 ptr->SetHighEnergyLimit(120 * MeV); 98 } 99 100 AddEmModel(1, EmModel()); 101 } 102 103 } 104 105 //------------------------------------------------------------------------------ 106 void G4DNADoubleIonisation::PrintInfo() 107 { 108 if (EmModel(1)) { 109 G4cout << " Total cross sections computed from " << EmModel(0)->GetName() 110 << " and " << EmModel(1)->GetName() << " models" << G4endl; 111 } else { 112 G4cout << " Total cross sections computed from " << EmModel()->GetName() 113 << G4endl; 114 } 115 } 116