| Geant4 Cross Reference |
1 // 1 //
2 // ******************************************* 2 // ********************************************************************
3 // * License and Disclaimer 3 // * License and Disclaimer *
4 // * 4 // * *
5 // * The Geant4 software is copyright of th 5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided 6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License 7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/ 8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. 9 // * include a list of copyright holders. *
10 // * 10 // * *
11 // * Neither the authors of this software syst 11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing fin 12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warran 13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assum 14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file 15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitatio 16 // * for the full disclaimer and the limitation of liability. *
17 // * 17 // * *
18 // * This code implementation is the result 18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboratio 19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distri 20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you ag 21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publicati 22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Sof 23 // * acceptance of all terms of the Geant4 Software license. *
24 // ******************************************* 24 // ********************************************************************
25 // 25 //
26 // Geant4 class file << 26 // $Id: G4AdjointhMultipleScattering.cc 73245 2013-08-22 13:13:20Z gcosmo $
>> 27 //
>> 28
>> 29 //
>> 30 // GEANT4 Class file
27 // 31 //
28 // File name: G4AdjointhMultipleScattering 32 // File name: G4AdjointhMultipleScattering
29 // 33 //
30 // Author: Desorgher Laurent 34 // Author: Desorgher Laurent
31 // 35 //
32 // Creation date: 03.06.2009 cloned from G4hMu << 36 // Creation date: 03.06.2009 cloned from G4hMultipleScattering by U.Laszlo with slight modification for adjoint_ion.
33 // slight modification for adjoint_ion. <<
34 // 37 //
35 // ------------------------------------------- 38 // -----------------------------------------------------------------------------
>> 39 //
>> 40 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
>> 41 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
36 42
37 #include "G4AdjointhMultipleScattering.hh" 43 #include "G4AdjointhMultipleScattering.hh"
38 <<
39 #include "G4MscStepLimitType.hh" <<
40 #include "G4SystemOfUnits.hh" 44 #include "G4SystemOfUnits.hh"
41 #include "G4UrbanMscModel.hh" 45 #include "G4UrbanMscModel.hh"
>> 46 #include "G4MscStepLimitType.hh"
42 47
43 G4AdjointhMultipleScattering::G4AdjointhMultip << 48 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
44 const G4String& processName) << 49
>> 50 using namespace std;
>> 51
>> 52 G4AdjointhMultipleScattering::G4AdjointhMultipleScattering(const G4String& processName)
45 : G4VMultipleScattering(processName) 53 : G4VMultipleScattering(processName)
46 {} << 54 {
>> 55 isInitialized = false;
>> 56 isIon = false;
>> 57 SetStepLimitType(fMinimal);
>> 58
>> 59 dtrl=0.;
>> 60 lambdalimit=0.;
>> 61 samplez=0.;
>> 62 }
47 63
48 //....oooOO0OOooo........oooOO0OOooo........oo 64 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
49 G4AdjointhMultipleScattering::~G4AdjointhMulti << 65
>> 66 G4AdjointhMultipleScattering::~G4AdjointhMultipleScattering()
>> 67 {}
50 68
51 //....oooOO0OOooo........oooOO0OOooo........oo 69 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
52 void G4AdjointhMultipleScattering::ProcessDesc << 70
>> 71 G4bool G4AdjointhMultipleScattering::IsApplicable (const G4ParticleDefinition& p)
53 { 72 {
54 out << "Inverse multiple scattering process << 73 return (p.GetPDGCharge() != 0.0 && !p.IsShortLived());
55 StreamProcessInfo(out); <<
56 } 74 }
57 75
58 //....oooOO0OOooo........oooOO0OOooo........oo 76 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
59 void G4AdjointhMultipleScattering::StreamProce << 77
>> 78 void G4AdjointhMultipleScattering::InitialiseProcess(const G4ParticleDefinition* p)
60 { 79 {
61 out << " RangeFactor= " << RangeFactor( << 80 // Modification of parameters between runs
62 << ", step limit type: " << StepLimitTyp << 81 if(isInitialized) {
63 << ", lateralDisplacement: " << LateralD << 82 if (p->GetParticleType() != "adjoint_nucleus" && p->GetPDGMass() < GeV) {
64 << ", skin= " << Skin() << G4endl; << 83 mscUrban->SetStepLimitType(StepLimitType());
>> 84 mscUrban->SetLateralDisplasmentFlag(LateralDisplasmentFlag());
>> 85 mscUrban->SetSkin(Skin());
>> 86 mscUrban->SetRangeFactor(RangeFactor());
>> 87 mscUrban->SetGeomFactor(GeomFactor());
>> 88 }
>> 89 return;
>> 90 }
>> 91
>> 92 // defaults for ions, which cannot be overwritten
>> 93 if (p->GetParticleType() == "adjoint_nucleus" || p->GetPDGMass() > GeV) {
>> 94 SetStepLimitType(fMinimal);
>> 95 SetLateralDisplasmentFlag(false);
>> 96 //SetBuildLambdaTable(false);
>> 97 if(p->GetParticleType() == "adjoint_nucleus") isIon = true;
>> 98 }
>> 99
>> 100 // initialisation of parameters
>> 101 G4String part_name = p->GetParticleName();
>> 102 mscUrban = new G4UrbanMscModel();
>> 103
>> 104 mscUrban->SetStepLimitType(StepLimitType());
>> 105 mscUrban->SetLateralDisplasmentFlag(LateralDisplasmentFlag());
>> 106 mscUrban->SetSkin(Skin());
>> 107 mscUrban->SetRangeFactor(RangeFactor());
>> 108 mscUrban->SetGeomFactor(GeomFactor());
>> 109
>> 110 AddEmModel(1,mscUrban);
>> 111 isInitialized = true;
65 } 112 }
66 113
67 //....oooOO0OOooo........oooOO0OOooo........oo 114 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
68 G4bool G4AdjointhMultipleScattering::IsApplica << 115
>> 116 void G4AdjointhMultipleScattering::PrintInfo()
69 { 117 {
70 return (p.GetPDGCharge() != 0.0 && !p.IsShor << 118 G4cout << " RangeFactor= " << RangeFactor()
>> 119 << ", step limit type: " << StepLimitType()
>> 120 << ", lateralDisplacement: " << LateralDisplasmentFlag()
>> 121 << ", skin= " << Skin()
>> 122 // << ", geomFactor= " << GeomFactor()
>> 123 << G4endl;
71 } 124 }
72 125
73 //....oooOO0OOooo........oooOO0OOooo........oo 126 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
74 void G4AdjointhMultipleScattering::InitialiseP << 127
75 const G4ParticleDefinition*) << 128 /*G4double G4AdjointhMultipleScattering::AlongStepGetPhysicalInteractionLength(
>> 129 const G4Track& track,
>> 130 double,
>> 131 G4double currentMinimalStep,
>> 132 G4double& currentSafety,
>> 133 G4GPILSelection* selection)
76 { 134 {
77 if(fIsInitialized) << 135 // get Step limit proposed by the process
78 { << 136 valueGPILSelectionMSC = NotCandidateForSelection;
79 return; << 137
80 } << 138 G4double escaled = track.GetKineticEnergy();
81 AddEmModel(1, new G4UrbanMscModel()); << 139 if(isIon) escaled *= track.GetDynamicParticle()->GetMass()/proton_mass_c2;
82 fIsInitialized = true; << 140
>> 141 G4double steplength = GetMscContinuousStepLimit(track,
>> 142 escaled,
>> 143 currentMinimalStep,
>> 144 currentSafety);
>> 145 // G4cout << "StepLimit= " << steplength << G4endl;
>> 146 // set return value for G4GPILSelection
>> 147 *selection = valueGPILSelectionMSC;
>> 148 return steplength;
83 } 149 }
>> 150 */
>> 151 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
84 152