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