- Timestamp:
- Dec 22, 2010, 3:52:27 PM (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/source/processes/hadronic/models/abrasion/src/G4WilsonAbrasionModel.cc
r1340 r1347 34 34 // ******************************************************************** 35 35 // 36 // %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %36 // %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 37 37 // 38 38 // MODULE: G4WilsonAbrasionModel.cc … … 45 45 // Contract: 17191/03/NL/LvH 46 46 // 47 // %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %47 // %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 48 48 // 49 49 // CHANGE HISTORY … … 72 72 // Additional clauses have been included in 73 73 // G4WilsonAbrasionModel::GetNucleonInducedExcitation 74 // Previously it was possible to get sqrt of negative number as Wilson algorithm75 // not properly defined if either:74 // Previously it was possible to get sqrt of negative number as Wilson 75 // algorithm not properly defined if either: 76 76 // rT > rP && rsq < rTsq - rPsq) or (rP > rT && rsq < rPsq - rTsq) 77 77 // 78 // %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %79 /////////////////////////////////////////////////////////////////////////////// /80 // 78 // %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 79 /////////////////////////////////////////////////////////////////////////////// 80 81 81 #include "G4WilsonAbrasionModel.hh" 82 82 #include "G4WilsonRadius.hh" … … 99 99 #include "G4IonTable.hh" 100 100 #include "globals.hh" 101 //////////////////////////////////////////////////////////////////////////////// 102 // 101 102 103 103 G4WilsonAbrasionModel::G4WilsonAbrasionModel (G4bool useAblation1) 104 104 :G4HadronicInteraction("G4WilsonAbrasion") 105 105 { 106 // 107 // 108 // Send message to stdout to advise that the G4Abrasion model is being used. 109 // 106 // Send message to stdout to advise that the G4Abrasion model is being used. 110 107 PrintWelcomeMessage(); 111 // 112 // 113 // Set the default verbose level to 0 - no output. 114 // 108 109 // Set the default verbose level to 0 - no output. 115 110 verboseLevel = 0; 116 111 useAblation = useAblation1; 117 // 118 // 119 // No de-excitation handler has been supplied - define the default handler. 120 // 112 113 // No de-excitation handler has been supplied - define the default handler. 114 121 115 theExcitationHandler = new G4ExcitationHandler; 122 116 theExcitationHandlerx = new G4ExcitationHandler; … … 146 140 theExcitationHandlerx->SetMaxAandZForFermiBreakUp(12, 6); 147 141 } 148 // 149 // 150 // Set the minimum and maximum range for the model (despite nomanclature, this 151 // is in energy per nucleon number). 152 // 142 143 // Set the minimum and maximum range for the model (despite nomanclature, 144 // this is in energy per nucleon number). 145 153 146 SetMinEnergy(70.0*MeV); 154 147 SetMaxEnergy(10.1*GeV); 155 148 isBlocked = false; 156 // 157 // 158 // npK, when mutiplied by the nuclear Fermi momentum, determines the range of 159 // momentum over which the secondary nucleon momentum is sampled. 160 // 161 npK 162 B 163 third 164 fradius 165 conserveEnergy 149 150 // npK, when mutiplied by the nuclear Fermi momentum, determines the range of 151 // momentum over which the secondary nucleon momentum is sampled. 152 153 r0sq = 0.0; 154 npK = 5.0; 155 B = 10.0 * MeV; 156 third = 1.0 / 3.0; 157 fradius = 0.99; 158 conserveEnergy = false; 166 159 conserveMomentum = true; 167 160 } 168 //////////////////////////////////////////////////////////////////////////////// 169 // 161 162 170 163 G4WilsonAbrasionModel::G4WilsonAbrasionModel (G4ExcitationHandler *aExcitationHandler) 171 164 {
Note: See TracChangeset
for help on using the changeset viewer.