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 // 26 // 27 // ******************************************* 27 // ********************************************************************** 28 // * 28 // * * 29 // * GEANT 4 xray_telescope 29 // * GEANT 4 xray_telescope advanced example * 30 // * 30 // * * 31 // * MODULE: XrayTelSteppingAction. 31 // * MODULE: XrayTelSteppingAction.cc * 32 // * ------- 32 // * ------- * 33 // * 33 // * * 34 // * Version: 0.4 34 // * Version: 0.4 * 35 // * Date: 06/11/00 35 // * Date: 06/11/00 * 36 // * Author: R Nartallo 36 // * Author: R Nartallo * 37 // * Organisation: ESA/ESTEC, Noordwijk, 37 // * Organisation: ESA/ESTEC, Noordwijk, THe Netherlands * 38 // * 38 // * * 39 // ******************************************* 39 // ********************************************************************** 40 // 40 // 41 // CHANGE HISTORY 41 // CHANGE HISTORY 42 // -------------- 42 // -------------- 43 // 43 // 44 // 05.12.2001 R. Nartallo 44 // 05.12.2001 R. Nartallo 45 // - Return condition for entering detector (c 45 // - Return condition for entering detector (cf LowEnTest) 46 // - Remove line to kill track if too many ste 46 // - Remove line to kill track if too many steps 47 // 47 // 48 // 07.11.2001 M.G. Pia 48 // 07.11.2001 M.G. Pia 49 // - Modified the analysis management 49 // - Modified the analysis management 50 // - Small design iteration 50 // - Small design iteration 51 // 51 // 52 // 30.11.2000 R. Nartallo 52 // 30.11.2000 R. Nartallo 53 // - Add pre-processor directives to compile w 53 // - Add pre-processor directives to compile without analysis option 54 // 54 // 55 // 16.11.2000 A. Pfeiffer 55 // 16.11.2000 A. Pfeiffer 56 // - Implementation of analysis manager call 56 // - Implementation of analysis manager call 57 // 57 // 58 // 06.11.2000 R.Nartallo 58 // 06.11.2000 R.Nartallo 59 // - First implementation of xray_telescope Ph 59 // - First implementation of xray_telescope Physics list 60 // - Based on Chandra and XMM models 60 // - Based on Chandra and XMM models 61 // 61 // 62 // ******************************************* 62 // ********************************************************************** 63 63 64 #include "XrayTelSteppingAction.hh" 64 #include "XrayTelSteppingAction.hh" 65 #include "globals.hh" 65 #include "globals.hh" 66 #include "G4Track.hh" 66 #include "G4Track.hh" 67 #include "G4Step.hh" 67 #include "G4Step.hh" 68 #include "G4Threading.hh" 68 #include "G4Threading.hh" 69 69 70 #include "XrayTelAnalysis.hh" 70 #include "XrayTelAnalysis.hh" 71 71 72 XrayTelSteppingAction::XrayTelSteppingAction() 72 XrayTelSteppingAction::XrayTelSteppingAction() 73 {;} 73 {;} 74 74 75 75 76 XrayTelSteppingAction::~XrayTelSteppingAction( 76 XrayTelSteppingAction::~XrayTelSteppingAction() 77 { } 77 { } 78 78 79 79 80 void XrayTelSteppingAction::UserSteppingAction 80 void XrayTelSteppingAction::UserSteppingAction(const G4Step* step) 81 { 81 { 82 G4bool entering = false; 82 G4bool entering = false; 83 G4Track* track = step->GetTrack(); 83 G4Track* track = step->GetTrack(); 84 84 85 G4String volName; 85 G4String volName; 86 if (track->GetVolume()) volName = track->Ge 86 if (track->GetVolume()) volName = track->GetVolume()->GetName(); 87 G4String nextVolName; 87 G4String nextVolName; 88 if (track->GetNextVolume()) nextVolName = t 88 if (track->GetNextVolume()) nextVolName = track->GetNextVolume()->GetName(); 89 89 90 XrayTelAnalysis* analysis = XrayTelAnalysis: 90 XrayTelAnalysis* analysis = XrayTelAnalysis::getInstance(); 91 91 92 // Entering Detector 92 // Entering Detector 93 if (volName != "Detector_P" && nextVolName = 93 if (volName != "Detector_P" && nextVolName == "Detector_P") 94 { 94 { 95 entering = true; 95 entering = true; 96 analysis->Update(track->GetKineticEnergy 96 analysis->Update(track->GetKineticEnergy(),G4Threading::G4GetThreadId()); 97 } 97 } 98 98 99 // Do the analysis related to this step 99 // Do the analysis related to this step 100 100 101 analysis->analyseStepping(*track,entering); 101 analysis->analyseStepping(*track,entering); 102 102 103 } 103 } 104 104 105 105 106 106 107 107 108 108