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 /// \file field/field04/src/F04SteppingAction. 27 /// \file field/field04/src/F04SteppingAction.cc 28 /// \brief Implementation of the F04SteppingAc 28 /// \brief Implementation of the F04SteppingAction class 29 // 29 // 30 30 31 #include "F04SteppingAction.hh" << 31 #include "G4Track.hh" 32 32 33 #include "F04UserTrackInformation.hh" << 33 #include "F04SteppingAction.hh" >> 34 #include "G4SteppingManager.hh" 34 35 35 #include "G4LogicalVolumeStore.hh" << 36 #include "G4ParticleTypes.hh" 36 #include "G4ParticleTypes.hh" 37 #include "G4SteppingManager.hh" << 37 #include "G4LogicalVolumeStore.hh" 38 #include "G4Track.hh" << 38 >> 39 #include "F04UserTrackInformation.hh" 39 40 40 //....oooOO0OOooo........oooOO0OOooo........oo 41 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 41 42 42 void F04SteppingAction::UserSteppingAction(con 43 void F04SteppingAction::UserSteppingAction(const G4Step* theStep) 43 { 44 { 44 G4Track* theTrack = theStep->GetTrack(); 45 G4Track* theTrack = theStep->GetTrack(); 45 46 46 // Get pointers to test volumes (only once) 47 // Get pointers to test volumes (only once) 47 if (!fTargetVolume) { << 48 if ( ! fTargetVolume ) { 48 fTargetVolume = G4LogicalVolumeStore::GetI << 49 fTargetVolume 49 fTestPlaneVolume = G4LogicalVolumeStore::G << 50 = G4LogicalVolumeStore::GetInstance()->GetVolume("Target"); >> 51 fTestPlaneVolume >> 52 = G4LogicalVolumeStore::GetInstance()->GetVolume("TestPlane"); 50 } 53 } 51 54 52 if (theTrack->GetParentID() == 0) { << 55 if (theTrack->GetParentID()==0) { 53 // This is a primary track << 56 //This is a primary track 54 G4LogicalVolume* theVolume = << 57 G4LogicalVolume* theVolume 55 theStep->GetPreStepPoint()->GetPhysicalV << 58 = theStep->GetPreStepPoint()->GetPhysicalVolume()->GetLogicalVolume(); 56 if (theVolume != fTargetVolume) { 59 if (theVolume != fTargetVolume) { 57 theTrack->SetTrackStatus(fStopAndKill); << 60 theTrack->SetTrackStatus(fStopAndKill); 58 return; << 61 return; 59 } 62 } 60 } 63 } 61 64 62 // G4ParticleDefinition* particleType = the << 65 // G4ParticleDefinition* particleType = theTrack->GetDefinition(); 63 66 64 // check if it is entering the test volume 67 // check if it is entering the test volume 65 G4StepPoint* thePrePoint = theStep->GetPreSt 68 G4StepPoint* thePrePoint = theStep->GetPreStepPoint(); 66 G4VPhysicalVolume* thePrePV = thePrePoint->G 69 G4VPhysicalVolume* thePrePV = thePrePoint->GetPhysicalVolume(); 67 G4LogicalVolume* thePreLV = thePrePV->GetLog 70 G4LogicalVolume* thePreLV = thePrePV->GetLogicalVolume(); 68 71 69 G4LogicalVolume* thePostLV = nullptr; 72 G4LogicalVolume* thePostLV = nullptr; 70 G4StepPoint* thePostPoint = theStep->GetPost 73 G4StepPoint* thePostPoint = theStep->GetPostStepPoint(); 71 74 72 if (thePostPoint) { 75 if (thePostPoint) { 73 G4VPhysicalVolume* thePostPV = thePostPoin << 76 G4VPhysicalVolume* thePostPV = thePostPoint->GetPhysicalVolume(); 74 if (thePostPV) thePostLV = thePostPV->GetL << 77 if (thePostPV) thePostLV = thePostPV->GetLogicalVolume(); 75 } 78 } 76 79 77 if (thePostLV == fTestPlaneVolume && thePreL << 80 if (thePostLV == fTestPlaneVolume && 78 // G4double x = theTrack->GetPosition( << 81 thePreLV != fTestPlaneVolume) { 79 // G4double y = theTrack->GetPosition( << 82 >> 83 // G4double x = theTrack->GetPosition().x(); >> 84 // G4double y = theTrack->GetPosition().y(); 80 85 81 G4ThreeVector theMomentumDirection = theTr << 86 G4ThreeVector theMomentumDirection = theTrack-> 82 // G4double theTotalMomentum = theTrac << 87 GetDynamicParticle()->GetMomentumDirection(); 83 // << 88 // G4double theTotalMomentum = theTrack->GetDynamicParticle()-> >> 89 // GetTotalMomentum(); >> 90 >> 91 // then kill the track >> 92 theTrack->SetTrackStatus(fStopAndKill); >> 93 return; 84 94 85 // then kill the track << 86 theTrack->SetTrackStatus(fStopAndKill); << 87 return; << 88 } 95 } 89 96 90 // G4double z = theTrack->GetPosition().z() << 97 // G4double z = theTrack->GetPosition().z(); 91 98 92 auto trackInformation = (F04UserTrackInforma << 99 auto trackInformation = >> 100 (F04UserTrackInformation*)theTrack->GetUserInformation(); 93 101 94 if (trackInformation->GetTrackStatusFlag() ! 102 if (trackInformation->GetTrackStatusFlag() != reverse) { 95 if (thePreLV != fTargetVolume) { << 103 if (thePreLV != fTargetVolume ) { 96 if (theTrack->GetMomentumDirection().z() << 104 if ( theTrack-> GetMomentumDirection().z()>0.0 && 97 && theTrack->GetVertexMomentumDirect << 105 theTrack->GetVertexMomentumDirection().z()<0.0 ) 98 { << 106 { 99 trackInformation->SetTrackStatusFlag(r << 107 trackInformation->SetTrackStatusFlag(reverse); 100 } << 108 } 101 } << 109 } 102 } 110 } 103 111 104 // check if it is alive 112 // check if it is alive 105 if (theTrack->GetTrackStatus() == fAlive) { << 113 if (theTrack->GetTrackStatus() == fAlive) { return; } 106 return; << 107 } << 108 114 109 if (thePostPoint->GetProcessDefinedStep() != 115 if (thePostPoint->GetProcessDefinedStep() != nullptr) { 110 if (thePostPoint->GetProcessDefinedStep()- << 116 if (thePostPoint->GetProcessDefinedStep()->GetProcessName() != "Decay") >> 117 return; 111 } 118 } >> 119 112 } 120 } 113 121