| 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 // S. Incerti (incerti@cenbg.in2p3.fr)
27 //
28
29 #include "G4EmDNAPhysics.hh"
30 #include "G4EmDNABuilder.hh"
31 #include "G4SystemOfUnits.hh"
32
33 // ions
34 #include "G4Alpha.hh"
35 #include "G4DNAGenericIonsManager.hh"
36
37 // utilities
38 #include "G4EmParameters.hh"
39 #include "G4PhysicsListHelper.hh"
40 #include "G4BuilderType.hh"
41 #include "G4EmBuilder.hh"
42 #include "G4PhysListUtil.hh"
43
44 // factory
45 #include "G4PhysicsConstructorFactory.hh"
46 //
47 G4_DECLARE_PHYSCONSTR_FACTORY(G4EmDNAPhysics);
48
49 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
50
51 G4EmDNAPhysics::G4EmDNAPhysics(G4int ver, const G4String& name)
52 : G4VPhysicsConstructor(name)
53 {
54 // instantiate singletones for physics
55 G4PhysListUtil::InitialiseParameters();
56
57 // parameters for DNA and for option3 EM physics
58 SetVerboseLevel(ver);
59 G4EmParameters* param = G4EmParameters::Instance();
60 param->SetDefaults();
61 param->SetMinEnergy(10*CLHEP::eV);
62 param->SetMaxEnergy(600*CLHEP::MeV);
63 param->SetLowestElectronEnergy(0*CLHEP::eV);
64 param->SetNumberOfBinsPerDecade(20);
65 param->ActivateAngularGeneratorForIonisation(true);
66 param->SetStepFunction(0.2, 10*CLHEP::um);
67 param->SetStepFunctionMuHad(0.1, 50*CLHEP::um);
68 param->SetStepFunctionLightIons(0.1, 20*CLHEP::um);
69 param->SetStepFunctionIons(0.1, 1*CLHEP::um);
70 param->SetUseICRU90Data(true);
71 param->SetUseMottCorrection(true);
72 param->SetMscStepLimitType(fUseSafetyPlus);
73 param->SetMscSkin(3);
74 param->SetMscRangeFactor(0.08);
75 param->SetFluo(true);
76 param->SetAuger(true);
77 param->SetDeexcitationIgnoreCut(true);
78 param->SetPositronAtRestModelType(fAllisonPositronium);
79 param->ActivateDNA();
80 param->SetDNAFast(false);
81 param->SetDNAStationary(false);
82
83 SetPhysicsType(bElectromagnetic);
84 }
85
86 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
87
88 void G4EmDNAPhysics::ConstructParticle()
89 {
90 G4EmDNABuilder::ConstructDNAParticles();
91 }
92
93 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
94
95 void G4EmDNAPhysics::ConstructProcess()
96 {
97 // parameters
98 G4EmParameters* param = G4EmParameters::Instance();
99 const G4double emaxDNA = 1.*CLHEP::MeV;
100 const G4double emaxIonDNA = 300.*CLHEP::MeV;
101 const G4double emaxLightIonDNA = 400.*CLHEP::MeV;
102 const G4double eminBorn = 500.*CLHEP::keV;
103 const G4bool fast = param->DNAFast();
104 const G4bool st = param->DNAStationary();
105 if(verboseLevel > 1) {
106 G4cout << "### " << GetPhysicsName()
107 << " Construct Processes EmaxDNA(MeV)= "
108 << emaxDNA/CLHEP::MeV << "; useMSC: " << fast
109 << "; stationary: " << st << G4endl;
110 }
111 G4DNAGenericIonsManager* genericIonsManager
112 = G4DNAGenericIonsManager::Instance();
113
114 // standard physics
115 G4EmDNABuilder::ConstructStandardEmPhysics(emaxDNA, emaxIonDNA,
116 emaxIonDNA, emaxIonDNA,
117 dnaGS, fast);
118
119 // DNA physics
120 G4EmDNABuilder::ConstructDNAElectronPhysics(emaxDNA, 0, fast, st);
121 G4EmDNABuilder::ConstructDNAProtonPhysics(eminBorn, emaxIonDNA, 0, fast, st);
122 G4EmDNABuilder::ConstructDNAIonPhysics(emaxIonDNA, st);
123
124 G4ParticleDefinition* part = genericIonsManager->GetIon("hydrogen");
125 G4EmDNABuilder::ConstructDNALightIonPhysics(part, 0, 0, emaxIonDNA, fast, st);
126
127 part = G4Alpha::Alpha();
128 G4EmDNABuilder::ConstructDNALightIonPhysics(part, 2, 0, emaxLightIonDNA, fast, st);
129
130 part = genericIonsManager->GetIon("alpha+");
131 G4EmDNABuilder::ConstructDNALightIonPhysics(part, 1, 0, emaxLightIonDNA, fast, st);
132
133 part = genericIonsManager->GetIon("helium");
134 G4EmDNABuilder::ConstructDNALightIonPhysics(part, 0, 0, emaxLightIonDNA, fast, st);
135 }
136
137 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
138