Geant4 Cross Reference |
1 // 1 // 2 // ******************************************* 2 // ******************************************************************** 3 // * License and Disclaimer << 3 // * DISCLAIMER * 4 // * 4 // * * 5 // * The Geant4 software is copyright of th << 5 // * The following disclaimer summarizes all the specific disclaimers * 6 // * the Geant4 Collaboration. It is provided << 6 // * of contributors to this software. The specific disclaimers,which * 7 // * conditions of the Geant4 Software License << 7 // * govern, are listed with their locations in: * 8 // * LICENSE and available at http://cern.ch/ << 8 // * http://cern.ch/geant4/license * 9 // * include a list of copyright holders. << 10 // * 9 // * * 11 // * Neither the authors of this software syst 10 // * Neither the authors of this software system, nor their employing * 12 // * institutes,nor the agencies providing fin 11 // * institutes,nor the agencies providing financial support for this * 13 // * work make any representation or warran 12 // * work make any representation or warranty, express or implied, * 14 // * regarding this software system or assum 13 // * regarding this software system or assume any liability for its * 15 // * use. Please see the license in the file << 14 // * use. * 16 // * for the full disclaimer and the limitatio << 17 // * 15 // * * 18 // * This code implementation is the result << 16 // * This code implementation is the intellectual property of the * 19 // * technical work of the GEANT4 collaboratio << 17 // * GEANT4 collaboration. * 20 // * By using, copying, modifying or distri << 18 // * By copying, distributing or modifying the Program (or any work * 21 // * any work based on the software) you ag << 19 // * based on the Program) you indicate your acceptance of this * 22 // * use in resulting scientific publicati << 20 // * statement, and all its terms. * 23 // * acceptance of all terms of the Geant4 Sof << 24 // ******************************************* 21 // ******************************************************************** 25 // 22 // >> 23 // >> 24 // $Id: F01RunAction.cc,v 1.9 2005/12/06 13:58:04 gcosmo Exp $ >> 25 // GEANT4 tag $Name: geant4-08-00-patch-01 $ >> 26 // >> 27 // >> 28 >> 29 26 #include "F01RunAction.hh" 30 #include "F01RunAction.hh" >> 31 #include "F01RunMessenger.hh" 27 32 28 #include "G4CoupledTransportation.hh" << 29 #include "G4Electron.hh" << 30 #include "G4ParticleDefinition.hh" << 31 #include "G4ProcessManager.hh" << 32 #include "G4Run.hh" 33 #include "G4Run.hh" 33 #include "G4SystemOfUnits.hh" << 34 #include "G4UImanager.hh" 34 #include "G4Transportation.hh" << 35 #include "G4VVisManager.hh" 35 #include "globals.hh" << 36 #include "G4ios.hh" >> 37 #include <iomanip> 36 38 37 //....oooOO0OOooo........oooOO0OOooo........oo << 39 #include "Randomize.hh" 38 40 39 void F01RunAction::BeginOfRunAction(const G4Ru << 41 ////////////////////////////////////////////////////////////////////////////// 40 { << 41 G4cout << "### Run " << aRun->GetRunID() << << 42 42 43 G4cout << " Calling F01RunAction::ChangeLoop << 43 F01RunAction::F01RunAction() 44 ChangeLooperParameters(G4Electron::Definitio << 44 : saveRndm(0) >> 45 { >> 46 runMessenger = new F01RunMessenger(this); 45 } 47 } 46 48 47 //....oooOO0OOooo........oooOO0OOooo........oo << 49 //////////////////////////////////////////////////////////////////////////// 48 50 49 void F01RunAction::ChangeLooperParameters(cons << 51 F01RunAction::~F01RunAction() 50 { 52 { 51 if (particleDef == nullptr) particleDef = G4 << 53 delete runMessenger; 52 auto transport = FindTransportation(particle << 54 } 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 55 83 if (fNumberOfTrials > 0) { << 56 ///////////////////////////////////////////////////////////////////////////// 84 transport->SetThresholdTrials(fNumberOfT << 85 57 86 G4cout << "-- Changed number of Trials ( << 58 void F01RunAction::BeginOfRunAction(const G4Run* aRun) 87 } << 59 { 88 } << 60 G4cout << "### Run " << aRun->GetRunID() << " start." << G4endl; >> 61 >> 62 // save Rndm status >> 63 if (saveRndm > 0) >> 64 { >> 65 CLHEP::HepRandom::showEngineStatus(); >> 66 CLHEP::HepRandom::saveEngineStatus("beginOfRun.rndm"); >> 67 } >> 68 G4UImanager* UI = G4UImanager::GetUIpointer(); >> 69 >> 70 G4VVisManager* pVVisManager = G4VVisManager::GetConcreteInstance(); 89 71 90 if (transport == nullptr) { << 72 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 } 73 } 99 74 100 //....oooOO0OOooo........oooOO0OOooo........oo << 75 ///////////////////////////////////////////////////////////////////////////// 101 76 102 void F01RunAction::EndOfRunAction(const G4Run* 77 void F01RunAction::EndOfRunAction(const G4Run*) 103 { 78 { 104 if (fVerboseLevel > 1) << 79 if (G4VVisManager::GetConcreteInstance()) 105 G4cout << G4endl << G4endl << " ########## << 80 { 106 << " ########### " << G4endl << " # << 81 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 } << 139 << 140 if (reportError && !transport) { << 141 G4cerr << "Unable to find Transportation p << 142 << particleDef->GetParticleName() < << 143 << " ) " << G4endl; << 144 } 82 } 145 83 146 return transport; << 84 // save Rndm status >> 85 >> 86 if (saveRndm == 1) >> 87 { >> 88 CLHEP::HepRandom::showEngineStatus(); >> 89 CLHEP::HepRandom::saveEngineStatus("endOfRun.rndm"); >> 90 } 147 } 91 } >> 92 >> 93 // >> 94 // >> 95 //////////////////////////////////////////////////////////////////////// 148 96