Geant4 Cross Reference |
1 // 1 2 // ******************************************* 3 // * License and Disclaimer 4 // * 5 // * The Geant4 software is copyright of th 6 // * the Geant4 Collaboration. It is provided 7 // * conditions of the Geant4 Software License 8 // * LICENSE and available at http://cern.ch/ 9 // * include a list of copyright holders. 10 // * 11 // * Neither the authors of this software syst 12 // * institutes,nor the agencies providing fin 13 // * work make any representation or warran 14 // * regarding this software system or assum 15 // * use. Please see the license in the file 16 // * for the full disclaimer and the limitatio 17 // * 18 // * This code implementation is the result 19 // * technical work of the GEANT4 collaboratio 20 // * By using, copying, modifying or distri 21 // * any work based on the software) you ag 22 // * use in resulting scientific publicati 23 // * acceptance of all terms of the Geant4 Sof 24 // ******************************************* 25 // 26 /// \file OpNovice/src/OpNoviceStackingAction. 27 /// \brief Implementation of the OpNoviceStack 28 // 29 // 30 //....oooOO0OOooo........oooOO0OOooo........oo 31 //....oooOO0OOooo........oooOO0OOooo........oo 32 33 #include "OpNoviceStackingAction.hh" 34 35 #include "OpNoviceRun.hh" 36 37 #include "G4OpticalPhoton.hh" 38 #include "G4RunManager.hh" 39 #include "G4Track.hh" 40 #include "G4VProcess.hh" 41 #include "G4ios.hh" 42 43 //....oooOO0OOooo........oooOO0OOooo........oo 44 G4ClassificationOfNewTrack OpNoviceStackingAct 45 { 46 if (aTrack->GetDefinition() == G4OpticalPhot 47 { // particle is optical photon 48 if (aTrack->GetParentID() > 0) { // parti 49 if (aTrack->GetCreatorProcess()->GetProc 50 ++fScintillationCounter; 51 else if (aTrack->GetCreatorProcess()->Ge 52 ++fCerenkovCounter; 53 } 54 } 55 return fUrgent; 56 } 57 58 //....oooOO0OOooo........oooOO0OOooo........oo 59 void OpNoviceStackingAction::NewStage() 60 { 61 // G4cout << "Number of Scintillation photon 62 // << fScintillationCounter << G4endl; 63 // G4cout << "Number of Cerenkov photons pro 64 // << fCerenkovCounter << G4endl; 65 66 auto run = static_cast<OpNoviceRun*>(G4RunMa 67 run->AddScintillation((G4double)fScintillati 68 run->AddCerenkov((G4double)fCerenkovCounter) 69 } 70 71 //....oooOO0OOooo........oooOO0OOooo........oo 72 void OpNoviceStackingAction::PrepareNewEvent() 73 { 74 fScintillationCounter = 0; 75 fCerenkovCounter = 0; 76 } 77 //....oooOO0OOooo........oooOO0OOooo........oo 78