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 #include "F01RunAction.hh" 26 #include "F01RunAction.hh" >> 27 #include "globals.hh" >> 28 #include "G4SystemOfUnits.hh" >> 29 #include "G4Run.hh" 27 30 28 #include "G4CoupledTransportation.hh" << 29 #include "G4Electron.hh" << 30 #include "G4ParticleDefinition.hh" 31 #include "G4ParticleDefinition.hh" >> 32 #include "G4Electron.hh" 31 #include "G4ProcessManager.hh" 33 #include "G4ProcessManager.hh" 32 #include "G4Run.hh" << 33 #include "G4SystemOfUnits.hh" << 34 #include "G4Transportation.hh" 34 #include "G4Transportation.hh" 35 #include "globals.hh" << 35 #include "G4CoupledTransportation.hh" 36 36 37 //....oooOO0OOooo........oooOO0OOooo........oo << 37 F01RunAction::F01RunAction() { >> 38 theWarningEnergy = 1.0 * CLHEP::kiloelectronvolt; // Arbitrary >> 39 theImportantEnergy = 10.0 * CLHEP::kiloelectronvolt; // Arbitrary >> 40 theNumberOfTrials = 15; // Arbitrary >> 41 // Applications should determine these thresholds according to >> 42 // - physics requirements, and >> 43 // - the computing cost of continuing integration for looping tracks >> 44 } 38 45 39 void F01RunAction::BeginOfRunAction(const G4Ru << 46 F01RunAction::~F01RunAction() {} 40 { << 47 >> 48 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo..... >> 49 >> 50 void F01RunAction::BeginOfRunAction( const G4Run* aRun ) { 41 G4cout << "### Run " << aRun->GetRunID() << 51 G4cout << "### Run " << aRun->GetRunID() << " start." << G4endl; 42 52 43 G4cout << " Calling F01RunAction::ChangeLoop 53 G4cout << " Calling F01RunAction::ChangeLooperParameters() " << G4endl; 44 ChangeLooperParameters(G4Electron::Definitio << 54 ChangeLooperParameters( G4Electron::Definition() ); 45 } 55 } 46 56 47 //....oooOO0OOooo........oooOO0OOooo........oo << 57 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo..... 48 58 49 void F01RunAction::ChangeLooperParameters(cons << 59 void F01RunAction:: >> 60 ChangeLooperParameters(const G4ParticleDefinition* particleDef ) 50 { 61 { 51 if (particleDef == nullptr) particleDef = G4 << 62 if( particleDef == nullptr ) 52 auto transport = FindTransportation(particle << 63 particleDef = G4Electron::Definition(); 53 << 64 auto transportPair= findTransportation( particleDef ); 54 // Note that all particles that share the sa << 65 auto transport = transportPair.first; 55 // (or Coupled Transportation) process (wit << 66 auto coupledTransport = transportPair.second; 56 // will be affected by this change << 67 57 // -- to change these parameters by partic << 68 if( transport != nullptr ) 58 // must instead be changed/reset in a ( << 69 { 59 << 70 // Change the values of the looping particle parameters of Transportation 60 G4cout << " ChangeLooperParameters called wi << 71 if( theWarningEnergy >= 0.0 ) 61 << " transport process: "; << 72 transport->SetThresholdWarningEnergy( theWarningEnergy ); 62 if (transport != nullptr) { << 73 if( theImportantEnergy >= 0.0 ) 63 G4cout << transport->GetProcessName(); << 74 transport->SetThresholdImportantEnergy( theImportantEnergy ); 64 } << 75 if( theNumberOfTrials > 0 ) 65 else { << 76 transport->SetThresholdTrials( theNumberOfTrials ); 66 G4cout << " UNKNOWN -- it is neither G4Tra << 77 } 67 } << 78 else if( coupledTransport != nullptr ) 68 G4cout << G4endl; << 79 { 69 << 80 // Change the values for Coupled Transport 70 if (transport != nullptr) { << 81 if( theWarningEnergy >= 0.0 ) 71 if (fWarningEnergy >= 0.0) { << 82 coupledTransport->SetThresholdWarningEnergy( theWarningEnergy ); 72 transport->SetThresholdWarningEnergy(fWa << 83 if( theImportantEnergy >= 0.0 ) 73 G4cout << "-- Changed Threshold Warning << 84 coupledTransport->SetThresholdImportantEnergy( theImportantEnergy ); 74 << fWarningEnergy / CLHEP::MeV << << 85 if( theNumberOfTrials > 0 ) 75 } << 86 coupledTransport->SetThresholdTrials( theNumberOfTrials ); 76 if (fImportantEnergy >= 0.0) { << 77 transport->SetThresholdImportantEnergy(f << 78 << 79 G4cout << "-- Changed Threshold Importan << 80 << fImportantEnergy / CLHEP::MeV << 81 } << 82 << 83 if (fNumberOfTrials > 0) { << 84 transport->SetThresholdTrials(fNumberOfT << 85 << 86 G4cout << "-- Changed number of Trials ( << 87 } << 88 } << 89 << 90 if (transport == nullptr) { << 91 if (fWarningEnergy >= 0.0) << 92 G4cerr << " Unknown transport process> C << 93 if (fImportantEnergy >= 0.0) << 94 G4cerr << " Unknown transport process> C << 95 if (fNumberOfTrials > 0) << 96 G4cerr << " Unknown transport process> C << 97 } 87 } 98 } 88 } 99 89 100 //....oooOO0OOooo........oooOO0OOooo........oo << 90 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo..... 101 91 102 void F01RunAction::EndOfRunAction(const G4Run* << 92 void F01RunAction::EndOfRunAction( const G4Run* ) { 103 { << 93 if( theVerboseLevel > 1 ) 104 if (fVerboseLevel > 1) << 94 G4cout << G4endl << G4endl 105 G4cout << G4endl << G4endl << " ########## << 95 << " ########### Track Statistics for Transportation process(es) " 106 << " ########### " << G4endl << " # << 96 << " ########### " << G4endl 107 << " ####################### " << G << 97 << " ############################################## " 108 << 98 << " ####################### " << G4endl << G4endl; 109 auto transport = FindTransportation(G4Electr << 99 110 if (transport) { << 100 auto transportPair= findTransportation( G4Electron::Definition() ); 111 transport->PrintStatistics(G4cout); << 101 auto transport = transportPair.first; 112 } << 102 auto coupledTransport = transportPair.second; >> 103 if( transport) { transport->PrintStatistics(G4cout); } >> 104 else if( coupledTransport) { coupledTransport->PrintStatistics(G4cout); } 113 } 105 } 114 106 115 //....oooOO0OOooo........oooOO0OOooo........oo 107 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 116 108 117 G4Transportation* F01RunAction::FindTransporta << 109 std::pair<G4Transportation*, G4CoupledTransportation*> 118 << 110 F01RunAction::findTransportation( const G4ParticleDefinition* particleDef, >> 111 bool reportError ) 119 { 112 { 120 const auto* partPM = particleDef->GetProcess << 113 const auto *partPM= particleDef->GetProcessManager(); 121 << 114 122 G4VProcess* partTransport = partPM->GetProce 115 G4VProcess* partTransport = partPM->GetProcess("Transportation"); 123 auto transport = dynamic_cast<G4Transportati << 116 auto transport= dynamic_cast<G4Transportation*>(partTransport); 124 << 125 if (!transport) { << 126 partTransport = partPM->GetProcess("Couple << 127 auto coupledTransport = dynamic_cast<G4Cou << 128 if (coupledTransport) { << 129 transport = coupledTransport; << 130 } << 131 else { << 132 partTransport = partPM->GetProcess("Tran << 133 auto transportWithMsc = dynamic_cast<G4C << 134 if (transportWithMsc) { << 135 transport = transportWithMsc; << 136 } << 137 } << 138 } << 139 << 140 if (reportError && !transport) { << 141 G4cerr << "Unable to find Transportation p << 142 << particleDef->GetParticleName() < << 143 << " ) " << G4endl; << 144 } << 145 117 146 return transport; << 118 partTransport = partPM->GetProcess("CoupledTransportation"); >> 119 auto coupledTransport= >> 120 dynamic_cast<G4CoupledTransportation*>(partTransport); >> 121 >> 122 if( reportError && !transport && !coupledTransport ) >> 123 { >> 124 G4cerr << "Unable to find Transportation process for particle type " >> 125 << particleDef->GetParticleName() >> 126 << " ( PDG code = " << particleDef->GetPDGEncoding() << " ) " >> 127 << G4endl; >> 128 } >> 129 >> 130 return >> 131 std::make_pair( transport, coupledTransport ); >> 132 // <G4Transportation*, G4CoupledTransportation*> 147 } 133 } 148 134