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 // >> 27 // $Id: F01RunAction.cc,v 1.10 2006/06/29 17:16:49 gunter Exp $ >> 28 // GEANT4 tag $Name: geant4-08-01-patch-01 $ >> 29 // >> 30 // >> 31 >> 32 26 #include "F01RunAction.hh" 33 #include "F01RunAction.hh" >> 34 #include "F01RunMessenger.hh" 27 35 28 #include "G4CoupledTransportation.hh" << 29 #include "G4Electron.hh" << 30 #include "G4ParticleDefinition.hh" << 31 #include "G4ProcessManager.hh" << 32 #include "G4Run.hh" 36 #include "G4Run.hh" 33 #include "G4SystemOfUnits.hh" << 37 #include "G4UImanager.hh" 34 #include "G4Transportation.hh" << 38 #include "G4VVisManager.hh" 35 #include "globals.hh" << 39 #include "G4ios.hh" >> 40 #include <iomanip> 36 41 37 //....oooOO0OOooo........oooOO0OOooo........oo << 42 #include "Randomize.hh" 38 43 39 void F01RunAction::BeginOfRunAction(const G4Ru << 44 ////////////////////////////////////////////////////////////////////////////// 40 { << 41 G4cout << "### Run " << aRun->GetRunID() << << 42 45 43 G4cout << " Calling F01RunAction::ChangeLoop << 46 F01RunAction::F01RunAction() 44 ChangeLooperParameters(G4Electron::Definitio << 47 : saveRndm(0) >> 48 { >> 49 runMessenger = new F01RunMessenger(this); 45 } 50 } 46 51 47 //....oooOO0OOooo........oooOO0OOooo........oo << 52 //////////////////////////////////////////////////////////////////////////// 48 53 49 void F01RunAction::ChangeLooperParameters(cons << 54 F01RunAction::~F01RunAction() 50 { 55 { 51 if (particleDef == nullptr) particleDef = G4 << 56 delete runMessenger; 52 auto transport = FindTransportation(particle << 57 } 53 << 54 // Note that all particles that share the sa << 55 // (or Coupled Transportation) process (wit << 56 // will be affected by this change << 57 // -- to change these parameters by partic << 58 // must instead be changed/reset in a ( << 59 << 60 G4cout << " ChangeLooperParameters called wi << 61 << " transport process: "; << 62 if (transport != nullptr) { << 63 G4cout << transport->GetProcessName(); << 64 } << 65 else { << 66 G4cout << " UNKNOWN -- it is neither G4Tra << 67 } << 68 G4cout << G4endl; << 69 << 70 if (transport != nullptr) { << 71 if (fWarningEnergy >= 0.0) { << 72 transport->SetThresholdWarningEnergy(fWa << 73 G4cout << "-- Changed Threshold Warning << 74 << fWarningEnergy / CLHEP::MeV << << 75 } << 76 if (fImportantEnergy >= 0.0) { << 77 transport->SetThresholdImportantEnergy(f << 78 << 79 G4cout << "-- Changed Threshold Importan << 80 << fImportantEnergy / CLHEP::MeV << 81 } << 82 58 83 if (fNumberOfTrials > 0) { << 59 ///////////////////////////////////////////////////////////////////////////// 84 transport->SetThresholdTrials(fNumberOfT << 85 60 86 G4cout << "-- Changed number of Trials ( << 61 void F01RunAction::BeginOfRunAction(const G4Run* aRun) 87 } << 62 { 88 } << 63 G4cout << "### Run " << aRun->GetRunID() << " start." << G4endl; >> 64 >> 65 // save Rndm status >> 66 if (saveRndm > 0) >> 67 { >> 68 CLHEP::HepRandom::showEngineStatus(); >> 69 CLHEP::HepRandom::saveEngineStatus("beginOfRun.rndm"); >> 70 } >> 71 G4UImanager* UI = G4UImanager::GetUIpointer(); >> 72 >> 73 G4VVisManager* pVVisManager = G4VVisManager::GetConcreteInstance(); 89 74 90 if (transport == nullptr) { << 75 if(pVVisManager) UI->ApplyCommand("/vis/scene/notifyHandlers"); 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 } << 98 } 76 } 99 77 100 //....oooOO0OOooo........oooOO0OOooo........oo << 78 ///////////////////////////////////////////////////////////////////////////// 101 79 102 void F01RunAction::EndOfRunAction(const G4Run* 80 void F01RunAction::EndOfRunAction(const G4Run*) 103 { 81 { 104 if (fVerboseLevel > 1) << 82 if (G4VVisManager::GetConcreteInstance()) 105 G4cout << G4endl << G4endl << " ########## << 83 { 106 << " ########### " << G4endl << " # << 84 G4UImanager::GetUIpointer()->ApplyCommand("/vis/viewer/update"); 107 << " ####################### " << G << 108 << 109 auto transport = FindTransportation(G4Electr << 110 if (transport) { << 111 transport->PrintStatistics(G4cout); << 112 } << 113 } << 114 << 115 //....oooOO0OOooo........oooOO0OOooo........oo << 116 << 117 G4Transportation* F01RunAction::FindTransporta << 118 << 119 { << 120 const auto* partPM = particleDef->GetProcess << 121 << 122 G4VProcess* partTransport = partPM->GetProce << 123 auto transport = dynamic_cast<G4Transportati << 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 } 85 } 139 86 140 if (reportError && !transport) { << 87 // save Rndm status 141 G4cerr << "Unable to find Transportation p << 88 142 << particleDef->GetParticleName() < << 89 if (saveRndm == 1) 143 << " ) " << G4endl; << 90 { 144 } << 91 CLHEP::HepRandom::showEngineStatus(); 145 << 92 CLHEP::HepRandom::saveEngineStatus("endOfRun.rndm"); 146 return transport; << 93 } 147 } 94 } >> 95 >> 96 // >> 97 // >> 98 //////////////////////////////////////////////////////////////////////// 148 99