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 /// \file electromagnetic/TestEm0/DirectAccess 26 /// \file electromagnetic/TestEm0/DirectAccess.cc 27 /// \brief Main program of the electromagnetic 27 /// \brief Main program of the electromagnetic/TestEm0 example 28 // 28 // 29 // 29 // 30 // << 30 // $Id: DirectAccess.cc 74920 2013-10-24 14:21:59Z gcosmo $ >> 31 // 31 // ------------------------------------------- 32 // ------------------------------------------------------------ 32 // 33 // 33 // To print cross sections per atom and mean 34 // To print cross sections per atom and mean free path for simple material 34 // 35 // 35 #include "G4BetheBlochModel.hh" << 36 #include "G4BetheHeitlerModel.hh" << 37 #include "G4BraggModel.hh" << 38 #include "G4DataVector.hh" << 39 #include "G4Electron.hh" << 40 #include "G4Gamma.hh" << 41 #include "G4KleinNishinaCompton.hh" << 42 #include "G4Material.hh" 36 #include "G4Material.hh" >> 37 >> 38 #include "G4PEEffectFluoModel.hh" >> 39 #include "G4KleinNishinaCompton.hh" >> 40 #include "G4BetheHeitlerModel.hh" >> 41 >> 42 #include "G4eeToTwoGammaModel.hh" >> 43 43 #include "G4MollerBhabhaModel.hh" 44 #include "G4MollerBhabhaModel.hh" >> 45 #include "G4SeltzerBergerModel.hh" >> 46 >> 47 #include "G4BetheBlochModel.hh" >> 48 #include "G4BraggModel.hh" >> 49 44 #include "G4MuBetheBlochModel.hh" 50 #include "G4MuBetheBlochModel.hh" 45 #include "G4MuBremsstrahlungModel.hh" 51 #include "G4MuBremsstrahlungModel.hh" 46 #include "G4MuPairProductionModel.hh" 52 #include "G4MuPairProductionModel.hh" 47 #include "G4MuonPlus.hh" << 53 48 #include "G4NistManager.hh" << 49 #include "G4PEEffectFluoModel.hh" << 50 #include "G4ParticleTable.hh" << 51 #include "G4Positron.hh" << 52 #include "G4Proton.hh" << 53 #include "G4SeltzerBergerModel.hh" << 54 #include "G4SystemOfUnits.hh" << 55 #include "G4UnitsTable.hh" << 56 #include "G4eeToTwoGammaModel.hh" << 57 #include "globals.hh" 54 #include "globals.hh" >> 55 #include "G4UnitsTable.hh" >> 56 #include "G4SystemOfUnits.hh" 58 57 59 int main() << 58 #include "G4Gamma.hh" 60 { << 59 #include "G4Positron.hh" 61 G4UnitDefinition::BuildUnitsTable(); << 60 #include "G4Electron.hh" >> 61 #include "G4Proton.hh" >> 62 #include "G4MuonPlus.hh" 62 63 63 G4ParticleDefinition* gamma = G4Gamma::Gamma << 64 int main() { 64 G4ParticleDefinition* posit = G4Positron::Po << 65 G4ParticleDefinition* elec = G4Electron::Ele << 66 G4ParticleDefinition* prot = G4Proton::Proto << 67 G4ParticleDefinition* muon = G4MuonPlus::Muo << 68 G4ParticleTable* partTable = G4ParticleTable << 69 partTable->SetReadiness(); << 70 65 71 G4DataVector cuts; << 66 G4UnitDefinition::BuildUnitsTable(); 72 cuts.push_back(1 * keV); << 73 67 74 // define materials 68 // define materials 75 // 69 // 76 G4Material* material = G4NistManager::Instan << 70 G4double Z, A; 77 71 78 G4cout << *(G4Material::GetMaterialTable()) << 72 G4Material* material = >> 73 new G4Material("Iodine", Z=53., A=126.90*g/mole, 4.93*g/cm3); 79 74 80 G4MaterialCutsCouple* couple = new G4Materia << 75 G4cout << *(G4Material::GetMaterialTable()) << G4endl; 81 couple->SetIndex(0); << 82 76 83 // work only for simple materials << 77 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 84 G4double Z = material->GetZ(); << 85 G4double A = material->GetA(); << 86 78 87 // initialise gamma processes (models) 79 // initialise gamma processes (models) 88 // 80 // >> 81 G4ParticleDefinition* gamma = G4Gamma::Gamma(); >> 82 89 G4VEmModel* phot = new G4PEEffectFluoModel() 83 G4VEmModel* phot = new G4PEEffectFluoModel(); 90 G4VEmModel* comp = new G4KleinNishinaCompton 84 G4VEmModel* comp = new G4KleinNishinaCompton(); 91 G4VEmModel* conv = new G4BetheHeitlerModel() << 85 G4VEmModel* conv = new G4BetheHeitlerModel(); 92 phot->Initialise(gamma, cuts); << 86 93 comp->Initialise(gamma, cuts); << 94 conv->Initialise(gamma, cuts); << 95 << 96 // valid pointer to a couple is needed for t << 97 phot->SetCurrentCouple(couple); << 98 << 99 // compute CrossSection per atom and MeanFre 87 // compute CrossSection per atom and MeanFreePath 100 // 88 // 101 G4double Emin = 1.01 * MeV, Emax = 2.01 * Me << 89 G4double Emin = 1.01*MeV, Emax = 2.01*MeV, dE = 100*keV; 102 90 103 G4cout << "\n #### Gamma : CrossSectionPerAt << 91 G4cout << "\n #### Gamma : CrossSectionPerAtom and MeanFreePath for " 104 << G4endl; << 92 << material->GetName() << G4endl; 105 G4cout << "\n Energy \t PhotoElec \t Compton 93 G4cout << "\n Energy \t PhotoElec \t Compton \t Conversion \t"; 106 G4cout << "\t PhotoElec \t Compton \t Conver << 94 G4cout << "\t PhotoElec \t Compton \t Conversion" << G4endl; 107 << 95 108 for (G4double Energy = Emin; Energy <= Emax; 96 for (G4double Energy = Emin; Energy <= Emax; Energy += dE) { 109 G4cout << "\n " << G4BestUnit(Energy, "Ene << 97 G4cout << "\n " << G4BestUnit (Energy, "Energy") 110 << G4BestUnit(phot->ComputeCrossSec << 98 << "\t" 111 << G4BestUnit(comp->ComputeCrossSec << 99 << G4BestUnit (phot->ComputeCrossSectionPerAtom(gamma,Energy,Z),"Surface") 112 << G4BestUnit(conv->ComputeCrossSec << 100 << "\t" 113 << G4BestUnit(phot->ComputeMeanFree << 101 << G4BestUnit (comp->ComputeCrossSectionPerAtom(gamma,Energy,Z),"Surface") 114 << G4BestUnit(comp->ComputeMeanFree << 102 << "\t" 115 << G4BestUnit(conv->ComputeMeanFree << 103 << G4BestUnit (conv->ComputeCrossSectionPerAtom(gamma,Energy,Z),"Surface") >> 104 << "\t \t" >> 105 << G4BestUnit (phot->ComputeMeanFreePath(gamma,Energy,material),"Length") >> 106 << "\t" >> 107 << G4BestUnit (comp->ComputeMeanFreePath(gamma,Energy,material),"Length") >> 108 << "\t" >> 109 << G4BestUnit (conv->ComputeMeanFreePath(gamma,Energy,material),"Length"); 116 } 110 } 117 111 118 G4cout << G4endl; 112 G4cout << G4endl; 119 113 >> 114 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... >> 115 120 // initialise positron annihilation (model) 116 // initialise positron annihilation (model) 121 // << 117 // >> 118 G4ParticleDefinition* posit = G4Positron::Positron(); >> 119 122 G4VEmModel* anni = new G4eeToTwoGammaModel() 120 G4VEmModel* anni = new G4eeToTwoGammaModel(); 123 anni->Initialise(posit, cuts); << 121 124 << 125 // compute CrossSection per atom and MeanFre 122 // compute CrossSection per atom and MeanFreePath 126 // 123 // 127 Emin = 1.01 * MeV; << 124 Emin = 1.01*MeV; Emax = 2.01*MeV; dE = 100*keV; 128 Emax = 2.01 * MeV; << 129 dE = 100 * keV; << 130 125 131 G4cout << "\n #### e+ annihilation : CrossSe 126 G4cout << "\n #### e+ annihilation : CrossSectionPerAtom and MeanFreePath" 132 << " for " << material->GetName() << 127 << " for " << material->GetName() << G4endl; 133 G4cout << "\n Energy \t e+ annihil \t"; 128 G4cout << "\n Energy \t e+ annihil \t"; 134 G4cout << "\t e+ annihil" << G4endl; << 129 G4cout << "\t e+ annihil" << G4endl; 135 << 130 136 for (G4double Energy = Emin; Energy <= Emax; 131 for (G4double Energy = Emin; Energy <= Emax; Energy += dE) { 137 G4cout << "\n " << G4BestUnit(Energy, "Ene << 132 G4cout << "\n " << G4BestUnit (Energy, "Energy") 138 << G4BestUnit(anni->ComputeCrossSec << 133 << "\t" 139 << G4BestUnit(anni->ComputeMeanFree << 134 << G4BestUnit (anni->ComputeCrossSectionPerAtom(posit,Energy,Z),"Surface") >> 135 << "\t \t" >> 136 << G4BestUnit (anni->ComputeMeanFreePath(posit,Energy,material),"Length"); 140 } 137 } 141 138 142 G4cout << G4endl; 139 G4cout << G4endl; 143 140 >> 141 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... >> 142 144 // initialise electron processes (models) 143 // initialise electron processes (models) 145 // << 144 // >> 145 G4ParticleDefinition* elec = G4Electron::Electron(); >> 146 146 G4VEmModel* ioni = new G4MollerBhabhaModel() 147 G4VEmModel* ioni = new G4MollerBhabhaModel(); 147 G4VEmModel* brem = new G4SeltzerBergerModel( 148 G4VEmModel* brem = new G4SeltzerBergerModel(); 148 ioni->Initialise(elec, cuts); << 149 149 brem->Initialise(elec, cuts); << 150 << 151 // compute CrossSection per atom and MeanFre 150 // compute CrossSection per atom and MeanFreePath 152 // 151 // 153 Emin = 1.01 * MeV; << 152 Emin = 1.01*MeV; Emax = 101.01*MeV; dE = 10*MeV; 154 Emax = 101.01 * MeV; << 153 G4double Ecut = 100*keV; 155 dE = 10 * MeV; << 156 G4double Ecut = 100 * keV; << 157 154 158 G4cout << "\n ####electron: CrossSection, Me 155 G4cout << "\n ####electron: CrossSection, MeanFreePath and StoppingPower" 159 << " for " << material->GetName() << << 156 << " for " << material->GetName() 160 << G4endl; << 157 << ";\tEnergy cut = " << G4BestUnit (Ecut, "Energy") << G4endl; 161 << 158 162 G4cout << "\n Energy \t ionization \t bremss 159 G4cout << "\n Energy \t ionization \t bremsstra \t"; 163 G4cout << "\t ionization \t bremsstra \t"; << 160 G4cout << "\t ionization \t bremsstra \t"; 164 G4cout << "\t ionization \t bremsstra" << G4 << 161 G4cout << "\t ionization \t bremsstra" << G4endl; 165 << 162 166 for (G4double Energy = Emin; Energy <= Emax; 163 for (G4double Energy = Emin; Energy <= Emax; Energy += dE) { 167 G4cout << "\n " << G4BestUnit(Energy, "Ene << 164 G4cout << "\n " << G4BestUnit (Energy, "Energy") 168 << G4BestUnit(ioni->ComputeCrossSec << 165 << "\t" 169 << "\t" << 166 << G4BestUnit (ioni->ComputeCrossSectionPerAtom(elec,Energy,Z,A,Ecut), 170 << G4BestUnit(brem->ComputeCrossSec << 167 "Surface") 171 << "\t \t" << 168 << "\t" 172 << G4BestUnit(ioni->ComputeMeanFree << 169 << G4BestUnit (brem->ComputeCrossSectionPerAtom(elec,Energy,Z,A,Ecut), 173 << G4BestUnit(brem->ComputeMeanFree << 170 "Surface") 174 << "\t \t" << 171 << "\t \t" 175 << G4BestUnit(ioni->ComputeDEDXPerV << 172 << G4BestUnit (ioni->ComputeMeanFreePath(elec,Energy,material,Ecut), 176 << "\t" << 173 "Length") 177 << G4BestUnit(brem->ComputeDEDXPerV << 174 << "\t" >> 175 << G4BestUnit (brem->ComputeMeanFreePath(elec,Energy,material,Ecut), >> 176 "Length") >> 177 << "\t \t" >> 178 << G4BestUnit (ioni->ComputeDEDXPerVolume(material,elec,Energy,Ecut), >> 179 "Energy/Length") >> 180 << "\t" >> 181 << G4BestUnit (brem->ComputeDEDXPerVolume(material,elec,Energy,Ecut), >> 182 "Energy/Length"); 178 } 183 } 179 << 184 180 G4cout << G4endl; 185 G4cout << G4endl; >> 186 >> 187 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 181 188 182 // initialise proton processes (models) << 183 // << 184 ioni = new G4BetheBlochModel(); << 185 ioni->Initialise(prot, cuts); << 186 189 >> 190 // initialise proton processes (models) >> 191 // >> 192 G4ParticleDefinition* prot = G4Proton::Proton(); >> 193 >> 194 ioni = new G4BetheBlochModel(prot); >> 195 187 // compute CrossSection per atom and MeanFre 196 // compute CrossSection per atom and MeanFreePath 188 // 197 // 189 Emin = 1.01 * MeV; << 198 Emin = 1.01*MeV; Emax = 102.01*MeV; dE = 10*MeV; 190 Emax = 102.01 * MeV; << 199 Ecut = 100*keV; 191 dE = 10 * MeV; << 192 Ecut = 100 * keV; << 193 200 194 G4cout << "\n #### proton : CrossSection, Me 201 G4cout << "\n #### proton : CrossSection, MeanFreePath and StoppingPower" 195 << " for " << material->GetName() << << 202 << " for " << material->GetName() 196 << G4endl; << 203 << ";\tEnergy cut = " << G4BestUnit (Ecut, "Energy") << G4endl; 197 << 204 198 G4cout << "\n Energy \t ionization \t"; 205 G4cout << "\n Energy \t ionization \t"; 199 G4cout << "\t ionization \t"; << 206 G4cout << "\t ionization \t"; 200 G4cout << "\t ionization" << G4endl; << 207 G4cout << "\t ionization" << G4endl; 201 << 208 202 for (G4double Energy = Emin; Energy <= Emax; 209 for (G4double Energy = Emin; Energy <= Emax; Energy += dE) { 203 G4cout << "\n " << G4BestUnit(Energy, "Ene << 210 G4cout << "\n " << G4BestUnit (Energy, "Energy") 204 << G4BestUnit(ioni->ComputeCrossSec << 211 << "\t" 205 << "\t \t" << 212 << G4BestUnit (ioni->ComputeCrossSectionPerAtom(prot,Energy,Z,A,Ecut), 206 << G4BestUnit(ioni->ComputeMeanFree << 213 "Surface") 207 << "\t \t" << 214 << "\t \t" 208 << G4BestUnit(ioni->ComputeDEDXPerV << 215 << G4BestUnit (ioni->ComputeMeanFreePath(prot,Energy,material,Ecut), >> 216 "Length") >> 217 << "\t \t" >> 218 << G4BestUnit (ioni->ComputeDEDXPerVolume(material,prot,Energy,Ecut), >> 219 "Energy/Length"); 209 } 220 } 210 << 221 211 G4cout << G4endl; 222 G4cout << G4endl; 212 << 223 213 // low energy : Bragg Model 224 // low energy : Bragg Model >> 225 214 ioni = new G4BraggModel(prot); 226 ioni = new G4BraggModel(prot); 215 ioni->Initialise(prot, cuts); << 227 216 << 217 // compute CrossSection per atom and MeanFre 228 // compute CrossSection per atom and MeanFreePath 218 // 229 // 219 Emin = 1.1 * keV; << 230 Emin = 1.1*keV; Emax = 2.01*MeV; dE = 300*keV; 220 Emax = 2.01 * MeV; << 231 Ecut = 10*keV; 221 dE = 300 * keV; << 232 222 Ecut = 10 * keV; << 223 << 224 G4cout << "\n #### proton : low energy model 233 G4cout << "\n #### proton : low energy model (Bragg) " 225 << ";\tEnergy cut = " << G4BestUnit(E << 234 << ";\tEnergy cut = " << G4BestUnit (Ecut, "Energy") << G4endl; 226 << 235 227 G4cout << "\n Energy \t ionization \t"; 236 G4cout << "\n Energy \t ionization \t"; 228 G4cout << "\t ionization \t"; << 237 G4cout << "\t ionization \t"; 229 G4cout << "\t ionization" << G4endl; << 238 G4cout << "\t ionization" << G4endl; 230 << 239 231 for (G4double Energy = Emin; Energy <= Emax; 240 for (G4double Energy = Emin; Energy <= Emax; Energy += dE) { 232 G4cout << "\n " << G4BestUnit(Energy, "Ene << 241 G4cout << "\n " << G4BestUnit (Energy, "Energy") 233 << G4BestUnit(ioni->ComputeCrossSec << 242 << "\t" 234 << "\t \t" << 243 << G4BestUnit (ioni->ComputeCrossSectionPerAtom(prot,Energy,Z,A,Ecut), 235 << G4BestUnit(ioni->ComputeMeanFree << 244 "Surface") 236 << "\t \t" << 245 << "\t \t" 237 << G4BestUnit(ioni->ComputeDEDXPerV << 246 << G4BestUnit (ioni->ComputeMeanFreePath(prot,Energy,material,Ecut), >> 247 "Length") >> 248 << "\t \t" >> 249 << G4BestUnit (ioni->ComputeDEDXPerVolume(material,prot,Energy,Ecut), >> 250 "Energy/Length"); 238 } 251 } 239 << 252 240 G4cout << G4endl; 253 G4cout << G4endl; 241 << 254 >> 255 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... >> 256 242 // initialise muon processes (models) 257 // initialise muon processes (models) 243 // << 258 // 244 ioni = new G4MuBetheBlochModel(); << 259 G4ParticleDefinition* muon = G4MuonPlus::MuonPlus(); 245 brem = new G4MuBremsstrahlungModel(); << 260 246 G4VEmModel* pair = new G4MuPairProductionMod << 261 ioni = new G4MuBetheBlochModel(muon); 247 ioni->Initialise(muon, cuts); << 262 brem = new G4MuBremsstrahlungModel(muon); 248 brem->Initialise(muon, cuts); << 263 G4VEmModel* pair = new G4MuPairProductionModel(muon); 249 pair->Initialise(muon, cuts); << 264 250 << 251 // compute CrossSection per atom and MeanFre 265 // compute CrossSection per atom and MeanFreePath 252 // 266 // 253 Emin = 1.01 * GeV; << 267 Emin = 1.01*GeV; Emax = 101.01*GeV; dE = 10*GeV; 254 Emax = 101.01 * GeV; << 268 Ecut = 10*MeV; 255 dE = 10 * GeV; << 256 Ecut = 10 * MeV; << 257 << 258 G4cout << "\n ####muon: CrossSection and Mea << 259 << ";\tEnergy cut = " << G4BestUnit(E << 260 269 >> 270 G4cout << "\n ####muon: CrossSection and MeanFreePath for " >> 271 << material->GetName() >> 272 << ";\tEnergy cut = " << G4BestUnit (Ecut, "Energy") << G4endl; >> 273 261 G4cout << "\n Energy \t ionization \t bremss 274 G4cout << "\n Energy \t ionization \t bremsstra \t pair_prod \t"; 262 G4cout << "\t ionization \t bremsstra \t pai << 275 G4cout << "\t ionization \t bremsstra \t pair_prod" << G4endl; 263 << 276 264 for (G4double Energy = Emin; Energy <= Emax; 277 for (G4double Energy = Emin; Energy <= Emax; Energy += dE) { 265 G4cout << "\n " << G4BestUnit(Energy, "Ene << 278 G4cout << "\n " << G4BestUnit (Energy, "Energy") 266 << G4BestUnit(ioni->ComputeCrossSec << 279 << "\t" 267 << "\t" << 280 << G4BestUnit (ioni->ComputeCrossSectionPerAtom(muon,Energy,Z,A,Ecut), 268 << G4BestUnit(brem->ComputeCrossSec << 281 "Surface") 269 << "\t" << 282 << "\t" 270 << G4BestUnit(pair->ComputeCrossSec << 283 << G4BestUnit (brem->ComputeCrossSectionPerAtom(muon,Energy,Z,A,Ecut), 271 << "\t \t" << 284 "Surface") 272 << G4BestUnit(ioni->ComputeMeanFree << 285 << "\t" 273 << G4BestUnit(brem->ComputeMeanFree << 286 << G4BestUnit (pair->ComputeCrossSectionPerAtom(muon,Energy,Z,A,Ecut), 274 << G4BestUnit(pair->ComputeMeanFree << 287 "Surface") >> 288 << "\t \t" >> 289 << G4BestUnit (ioni->ComputeMeanFreePath(muon,Energy,material,Ecut), >> 290 "Length") >> 291 << "\t" >> 292 << G4BestUnit (brem->ComputeMeanFreePath(muon,Energy,material,Ecut), >> 293 "Length") >> 294 << "\t" >> 295 << G4BestUnit (pair->ComputeMeanFreePath(muon,Energy,material,Ecut), >> 296 "Length"); 275 } 297 } 276 << 298 277 G4cout << G4endl; 299 G4cout << G4endl; 278 << 300 279 G4cout << "\n ####muon: StoppingPower for " << 301 G4cout << "\n ####muon: StoppingPower for " 280 << ";\tEnergy cut = " << G4BestUnit(E << 302 << material->GetName() 281 << 303 << ";\tEnergy cut = " << G4BestUnit (Ecut, "Energy") << G4endl; >> 304 282 G4cout << "\n Energy \t ionization \t bremss 305 G4cout << "\n Energy \t ionization \t bremsstra \t pair_prod \t" << G4endl; 283 << 306 284 for (G4double Energy = Emin; Energy <= Emax; 307 for (G4double Energy = Emin; Energy <= Emax; Energy += dE) { 285 G4cout << "\n " << G4BestUnit(Energy, "Ene << 308 G4cout << "\n " << G4BestUnit (Energy, "Energy") 286 << G4BestUnit(ioni->ComputeDEDXPerV << 309 << "\t" 287 << "\t" << 310 << G4BestUnit (ioni->ComputeDEDXPerVolume(material,muon,Energy,Ecut), 288 << G4BestUnit(brem->ComputeDEDXPerV << 311 "Energy/Length") 289 << "\t" << 312 << "\t" 290 << G4BestUnit(pair->ComputeDEDXPerV << 313 << G4BestUnit (brem->ComputeDEDXPerVolume(material,muon,Energy,Ecut), >> 314 "Energy/Length") >> 315 << "\t" >> 316 << G4BestUnit (pair->ComputeDEDXPerVolume(material,muon,Energy,Ecut), >> 317 "Energy/Length"); 291 } 318 } 292 << 319 293 G4cout << G4endl; 320 G4cout << G4endl; 294 return EXIT_SUCCESS; << 321 295 } << 296 << 297 //....oooOO0OOooo........oooOO0OOooo........oo 322 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... >> 323 >> 324 return EXIT_SUCCESS; >> 325 } 298 326