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 biasing/ReverseMC01/src/RMC01Adjoint 26 /// \file biasing/ReverseMC01/src/RMC01AdjointEventAction.cc 27 /// \brief Implementation of the RMC01AdjointE 27 /// \brief Implementation of the RMC01AdjointEventAction class 28 // 28 // 29 // 29 // 30 ////////////////////////////////////////////// 30 ////////////////////////////////////////////////////////////// 31 // Class Name: RMC01AdjointEventAc 31 // Class Name: RMC01AdjointEventAction 32 // Author: L. Desorgher 32 // Author: L. Desorgher 33 // Organisation: SpaceIT GmbH 33 // Organisation: SpaceIT GmbH 34 // Contract: ESA contract 21435/ 34 // Contract: ESA contract 21435/08/NL/AT 35 // Customer: ESA/ESTEC 35 // Customer: ESA/ESTEC 36 ////////////////////////////////////////////// 36 ////////////////////////////////////////////////////////////// 37 37 38 //....oooOO0OOooo........oooOO0OOooo........oo 38 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 39 //....oooOO0OOooo........oooOO0OOooo........oo 39 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 40 40 41 #include "RMC01AdjointEventAction.hh" 41 #include "RMC01AdjointEventAction.hh" 42 << 43 #include "RMC01AnalysisManager.hh" << 44 #include "RMC01SD.hh" << 45 << 46 #include "G4Colour.hh" << 47 #include "G4Event.hh" 42 #include "G4Event.hh" 48 #include "G4Polyline.hh" << 49 #include "G4Polymarker.hh" << 50 #include "G4PrimaryParticle.hh" << 51 #include "G4PrimaryVertex.hh" 43 #include "G4PrimaryVertex.hh" 52 #include "G4SDManager.hh" << 44 #include "G4PrimaryParticle.hh" 53 #include "G4Trajectory.hh" << 54 #include "G4TrajectoryContainer.hh" 45 #include "G4TrajectoryContainer.hh" 55 #include "G4TrajectoryPoint.hh" << 56 #include "G4UnitsTable.hh" << 57 #include "G4VTrajectory.hh" 46 #include "G4VTrajectory.hh" 58 #include "G4VVisManager.hh" 47 #include "G4VVisManager.hh" >> 48 #include "G4Polyline.hh" >> 49 #include "G4Polymarker.hh" >> 50 #include "G4Colour.hh" 59 #include "G4VisAttributes.hh" 51 #include "G4VisAttributes.hh" >> 52 #include "G4Trajectory.hh" >> 53 #include "G4TrajectoryPoint.hh" >> 54 #include "G4UnitsTable.hh" >> 55 #include "G4SDManager.hh" 60 #include "Randomize.hh" 56 #include "Randomize.hh" 61 << 62 #include <iomanip> 57 #include <iomanip> >> 58 #include "RMC01SD.hh" >> 59 #include "RMC01AnalysisManager.hh" 63 60 64 //....oooOO0OOooo........oooOO0OOooo........oo 61 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 65 62 66 RMC01AdjointEventAction::RMC01AdjointEventActi << 63 RMC01AdjointEventAction::RMC01AdjointEventAction() 67 { << 64 : G4UserEventAction() 68 ; << 65 {; 69 } 66 } 70 << 67 >> 68 71 //....oooOO0OOooo........oooOO0OOooo........oo 69 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 72 70 73 RMC01AdjointEventAction::~RMC01AdjointEventAct 71 RMC01AdjointEventAction::~RMC01AdjointEventAction() 74 { << 72 {; 75 ; << 76 } 73 } 77 74 78 //....oooOO0OOooo........oooOO0OOooo........oo 75 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 79 76 80 void RMC01AdjointEventAction::BeginOfEventActi << 77 void RMC01AdjointEventAction::BeginOfEventAction(const G4Event* anEvent ) 81 { << 78 { G4int i_event=anEvent->GetEventID(); 82 G4int i_event = anEvent->GetEventID(); << 79 G4bool print_nb=false; 83 G4bool print_nb = false; << 80 if (i_event <100) print_nb =true; 84 if (i_event < 100) << 81 else if (i_event<500 && (i_event/100)*100 == i_event) print_nb = true; 85 print_nb = true; << 82 else if (i_event<5000 && (i_event/500)*500 == i_event) print_nb = true; 86 else if (i_event < 500 && (i_event / 100) * << 83 else if ((i_event/5000)*5000 == i_event) print_nb = true; 87 print_nb = true; << 84 if (print_nb) G4cout<<"nb event "<<i_event<<std::endl; 88 else if (i_event < 5000 && (i_event / 500) * << 89 print_nb = true; << 90 else if ((i_event / 5000) * 5000 == i_event) << 91 print_nb = true; << 92 if (print_nb) G4cout << "nb event " << i_eve << 93 } 85 } 94 86 95 //....oooOO0OOooo........oooOO0OOooo........oo 87 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 96 88 97 void RMC01AdjointEventAction::EndOfEventAction << 89 void RMC01AdjointEventAction::EndOfEventAction(const G4Event* ) 98 { << 90 { ; 99 ; << 91 } 100 } << 101 92 102 //....oooOO0OOooo........oooOO0OOooo........oo 93 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 103 94