| 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 // Code developed by 26 // Code developed by
27 // Silvia Pozzi (1), silvia.pozzi@iss.it 27 // Silvia Pozzi (1), silvia.pozzi@iss.it
28 // Barbara Caccia (1), barbara.caccia@iss.it 28 // Barbara Caccia (1), barbara.caccia@iss.it
29 // Carlo Mancini Terracciano (2), carlo.mancin 29 // Carlo Mancini Terracciano (2), carlo.mancini.terracciano@roma1.infn.it
30 // (1) Istituto Superiore di Sanita' and INFN 30 // (1) Istituto Superiore di Sanita' and INFN Roma, Italy
31 // (2) Univ. La Sapienza and INFN Roma, Italy 31 // (2) Univ. La Sapienza and INFN Roma, Italy
32 32
33 #include "DetectorConstruction.hh" 33 #include "DetectorConstruction.hh"
34 #include "SteppingAction.hh" 34 #include "SteppingAction.hh"
35 35
36 #include <G4RunManager.hh> 36 #include <G4RunManager.hh>
37 #include <G4Step.hh> 37 #include <G4Step.hh>
38 #include <G4SystemOfUnits.hh> 38 #include <G4SystemOfUnits.hh>
39 39
40 void SteppingAction::UserSteppingAction(const 40 void SteppingAction::UserSteppingAction(const G4Step* aStep)
41 { 41 {
42 DetectorConstruction* detector = (DetectorCo 42 DetectorConstruction* detector = (DetectorConstruction*)
43 G4RunManager::GetRunManager() -> GetUserDe 43 G4RunManager::GetRunManager() -> GetUserDetectorConstruction();
44 44
45 G4double fFilterZ = detector -> GetFFilterZ( 45 G4double fFilterZ = detector -> GetFFilterZ();
46 G4double fFilterRadius = detector -> GetFFil 46 G4double fFilterRadius = detector -> GetFFilterRadius();
47 47
48 G4String particleName; 48 G4String particleName;
49 49
50 G4VPhysicalVolume* volume = aStep -> GetPost 50 G4VPhysicalVolume* volume = aStep -> GetPostStepPoint() -> GetTouchable() -> GetVolume();
51 51
52 const G4ThreeVector momentumDir = aStep -> G 52 const G4ThreeVector momentumDir = aStep -> GetTrack()-> GetMomentumDirection();
53 if (momentumDir.z() < 0) 53 if (momentumDir.z() < 0)
54 { 54 {
55 aStep->GetTrack()->SetTrackStatus(fStopAnd 55 aStep->GetTrack()->SetTrackStatus(fStopAndKill);
56 } 56 }
57 57
58 if ( volume != nullptr ) 58 if ( volume != nullptr )
59 { 59 {
60 G4String volumeName = volume -> GetName( 60 G4String volumeName = volume -> GetName();
61 61
62 if (volumeName == "Jaw1XPV" || volumeNam 62 if (volumeName == "Jaw1XPV" || volumeName == "Jaw2XPV"
63 || volumeName == "Jaw1YPV" || volumeName = 63 || volumeName == "Jaw1YPV" || volumeName == "Jaw2YPV" )
64 { 64 {
65 aStep->GetTrack()->SetTrackStatus(fStopAnd 65 aStep->GetTrack()->SetTrackStatus(fStopAndKill);
66 } 66 }
67 } 67 }
68 68
69 if (aStep -> GetPreStepPoint() -> GetTouchab 69 if (aStep -> GetPreStepPoint() -> GetTouchable() -> GetVolume() != nullptr )
70 { 70 {
71 if (aStep -> GetPreStepPoint() -> GetTou 71 if (aStep -> GetPreStepPoint() -> GetTouchable() -> GetVolume() -> GetName() == "TargetPV" ||
72 aStep -> GetPreStepPoint() -> GetTouchable 72 aStep -> GetPreStepPoint() -> GetTouchable() -> GetVolume() -> GetName() == "PriCollUpperPV" ||
73 aStep -> GetPreStepPoint() -> GetTouchable 73 aStep -> GetPreStepPoint() -> GetTouchable() -> GetVolume() -> GetName() == "PriCollMiddlePV" ||
74 aStep -> GetPreStepPoint() -> GetTouchable 74 aStep -> GetPreStepPoint() -> GetTouchable() -> GetVolume() -> GetName() == "PriCollLowererPV" ||
75 aStep -> GetPreStepPoint() -> GetTouchable 75 aStep -> GetPreStepPoint() -> GetTouchable() -> GetVolume() -> GetName() == "ffPV")
76 { 76 {
77 if (std::abs(momentumDir.getX()/momentumDi 77 if (std::abs(momentumDir.getX()/momentumDir.getZ()) > fFilterRadius/fFilterZ ||
78 std::abs(momentumDir.getY()/momentumDi 78 std::abs(momentumDir.getY()/momentumDir.getZ()) > fFilterRadius/fFilterZ )
79 { 79 {
80 aStep->GetTrack()->SetTrackStatus(fSto 80 aStep->GetTrack()->SetTrackStatus(fStopAndKill);
81 } 81 }
82 } 82 }
83 } 83 }
84 } 84 }
85 85