| 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 // G4AdjointSteppingAction class implementatio << 26 // $Id: G4AdjointSteppingAction.cc,v 1.3 2009/11/20 14:24:34 ldesorgh Exp $
>> 27 // GEANT4 tag $Name: geant4-09-03 $
27 // 28 //
28 // Author: L. Desorgher, SpaceIT GmbH << 29 /////////////////////////////////////////////////////////////////////////////
29 // Contract: ESA contract 21435/08/NL/AT << 30 // Class Name: G4AdjointSteppingAction
30 // Customer: ESA/ESTEC << 31 // Author: L. Desorgher
31 // ------------------------------------------- << 32 // Organisation: SpaceIT GmbH
>> 33 // Contract: ESA contract 21435/08/NL/AT
>> 34 // Customer: ESA/ESTEC
>> 35 /////////////////////////////////////////////////////////////////////////////
32 36
33 #include "G4AdjointSteppingAction.hh" 37 #include "G4AdjointSteppingAction.hh"
34 <<
35 #include "G4AdjointCrossSurfChecker.hh" <<
36 #include "G4AffineTransform.hh" <<
37 #include "G4PhysicalVolumeStore.hh" <<
38 #include "G4Track.hh" 38 #include "G4Track.hh"
39 << 39 #include "G4PhysicalVolumeStore.hh"
40 ////////////////////////////////////////////// << 40 #include "G4AffineTransform.hh"
>> 41 #include "G4AdjointCrossSurfChecker.hh"
>> 42 ////////////////////////////////////////////////////////////////////////////////////////////////
41 // 43 //
42 G4AdjointSteppingAction::G4AdjointSteppingActi 44 G4AdjointSteppingAction::G4AdjointSteppingAction()
43 { << 45 {
44 theG4AdjointCrossSurfChecker = G4AdjointCros 46 theG4AdjointCrossSurfChecker = G4AdjointCrossSurfChecker::GetInstance();
>> 47 did_adj_part_reach_ext_source =false;
>> 48 theUserAdjointSteppingAction =0;
45 } 49 }
>> 50 ////////////////////////////////////////////////////////////////////////////////////////////////
>> 51 //
>> 52 G4AdjointSteppingAction::~G4AdjointSteppingAction()
>> 53 {;}
46 54
47 ////////////////////////////////////////////// << 55 ////////////////////////////////////////////////////////////////////////////////////////////////
48 // 56 //
49 void G4AdjointSteppingAction::UserSteppingActi 57 void G4AdjointSteppingAction::UserSteppingAction(const G4Step* aStep)
50 { 58 {
51 G4Track* aTrack = aStep->GetTrack(); << 59 //Apply first the user adjoint stepping action
52 if (! is_adjoint_tracking_mode) // forward << 60 //---------------------------
53 { << 61 if (theUserAdjointSteppingAction) theUserAdjointSteppingAction->UserSteppingAction(aStep);
54 if (! did_one_adj_part_reach_ext_source_du <<
55 aTrack->SetTrackStatus(fStopAndKill); <<
56 return; <<
57 } <<
58 if (theUserFwdSteppingAction != nullptr) { <<
59 theUserFwdSteppingAction->UserSteppingAc <<
60 } <<
61 return; <<
62 } <<
63 62
64 // Apply first the user adjoint stepping act << 63 G4Track* aTrack =aStep->GetTrack();
65 // ----------------------------------------- << 64 G4double nb_nuc=1.;
66 <<
67 did_adj_part_reach_ext_source = false; <<
68 if (theUserAdjointSteppingAction != nullptr) <<
69 theUserAdjointSteppingAction->UserStepping <<
70 } <<
71 <<
72 G4double nb_nuc = 1.; <<
73 G4ParticleDefinition* thePartDef = aTrack->G 65 G4ParticleDefinition* thePartDef = aTrack->GetDefinition();
74 << 66
75 if (thePartDef->GetParticleType() == "adjoin << 67 if (thePartDef->GetParticleType() == "adjoint_nucleus"){
76 nb_nuc = G4double(thePartDef->GetBaryonNum << 68 nb_nuc=double(thePartDef->GetBaryonNumber());
77 } 69 }
78 << 70 //Kill conditions for adjoint particles reaching the maximum energy
79 // Kill conditions for adjoint particles rea << 71 //-----------------------------------------------------------------
80 // ----------------------------------------- << 72 if(aTrack->GetKineticEnergy() >= ext_sourceEMax*nb_nuc){
81 << 73 aTrack->SetTrackStatus(fStopAndKill);
82 if (aTrack->GetKineticEnergy() >= ext_source << 74 did_adj_part_reach_ext_source=false;
83 aTrack->SetTrackStatus(fStopAndKill); << 75 return;
84 did_adj_part_reach_ext_source = false; <<
85 return; <<
86 } 76 }
87 77
88 // Kill conditions for surface crossing << 78 G4double weight_factor = aTrack->GetWeight()/prim_weight;
89 // -------------------------------------- <<
90 79
>> 80 if ( (weight_factor>0 && weight_factor<=0) || weight_factor<= 1e-290 || weight_factor>1.e200)
>> 81 {
>> 82 //std::cout<<"Weight_factor problem! Value = "<<weight_factor<<std::endl;
>> 83 aTrack->SetTrackStatus(fStopAndKill);
>> 84 did_adj_part_reach_ext_source=false;
>> 85 return;
>> 86 }
>> 87
>> 88
>> 89 //Kill conditions for surface crossing
>> 90 //--------------------------------------
>> 91
91 G4String surface_name; 92 G4String surface_name;
92 G4double cos_to_surface; 93 G4double cos_to_surface;
93 G4bool GoingIn; 94 G4bool GoingIn;
94 G4ThreeVector crossing_pos; 95 G4ThreeVector crossing_pos;
95 if (theG4AdjointCrossSurfChecker->CrossingOn << 96 if (theG4AdjointCrossSurfChecker->CrossingOneOfTheRegisteredSurface(aStep, surface_name, crossing_pos, cos_to_surface, GoingIn) ){
96 aStep, surface_name, crossing_pos, cos << 97
97 { << 98 //G4cout<<"Test_step11"<<std::endl;
98 if (surface_name == "ExternalSource") { << 99 //G4cout<<surface_name<<std::endl;
99 // Registering still needed << 100 if (surface_name == "ExternalSource") {
100 did_adj_part_reach_ext_source = true; << 101 //Registering still needed
101 did_one_adj_part_reach_ext_source_during << 102 did_adj_part_reach_ext_source=true;
102 aTrack->SetTrackStatus(fStopAndKill); << 103 aTrack->SetTrackStatus(fStopAndKill);
103 << 104 //now register the adjoint particles reaching the external surface
104 // now register the adjoint particles re << 105 last_momentum =aTrack->GetMomentum();
105 last_momentum = aTrack->GetMomentum(); << 106 last_ekin=aTrack->GetKineticEnergy();
106 last_ekin = aTrack->GetKineticEnergy(); << 107 last_weight = aTrack->GetWeight();
107 last_weight = aTrack->GetWeight(); << 108 last_part_def = aTrack->GetDefinition();
108 last_part_def = aTrack->GetDefinition(); << 109 last_pos = crossing_pos;
109 last_pos = crossing_pos; << 110 return;
110 return; << 111 }
111 } << 112 else if (surface_name == "AdjointSource" && GoingIn) {
112 if (surface_name == "AdjointSource" && Goi << 113 did_adj_part_reach_ext_source=false;
113 did_adj_part_reach_ext_source = false; << 114 aTrack->SetTrackStatus(fStopAndKill);
114 aTrack->SetTrackStatus(fStopAndKill); << 115 return;
115 return; << 116 }
116 } <<
117 } <<
118 <<
119 // Check for reaching out of world <<
120 // <<
121 if (aStep->GetPostStepPoint()->GetStepStatus <<
122 did_adj_part_reach_ext_source = true; <<
123 did_one_adj_part_reach_ext_source_during_e <<
124 last_momentum = aTrack->GetMomentum(); <<
125 last_ekin = aTrack->GetKineticEnergy(); <<
126 last_weight = aTrack->GetWeight(); <<
127 last_part_def = aTrack->GetDefinition(); <<
128 last_pos = crossing_pos; <<
129 return; <<
130 } 117 }
131 } 118 }
>> 119
132 120