| Geant4 Cross Reference |
1 //
2 // ********************************************************************
3 // * License and Disclaimer *
4 // * *
5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. *
10 // * *
11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitation of liability. *
17 // * *
18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Software license. *
24 // ********************************************************************
25 //
26 //
27 /// \file GB06BOptnSplitAndKillByImportance.hh
28 /// \brief Definition of the GB06BOptnSplitAndKillByImportance class
29
30 #ifndef GB06BOptnSplitAndKillByImportance_hh
31 #define GB06BOptnSplitAndKillByImportance_hh 1
32
33 #include "G4ParticleChange.hh"
34 #include "G4ParticleChangeForNothing.hh"
35 #include "G4TouchableHandle.hh"
36 #include "G4VBiasingOperation.hh"
37 class G4BiasingProcessSharedData;
38 #include <map>
39
40 class GB06BOptnSplitAndKillByImportance : public G4VBiasingOperation
41 {
42 public:
43 // -- Constructor :
44 GB06BOptnSplitAndKillByImportance(G4String name);
45 // -- destructor:
46 virtual ~GB06BOptnSplitAndKillByImportance();
47
48 public:
49 // ----------------------------------------------
50 // -- Methods from G4VBiasingOperation interface:
51 // ----------------------------------------------
52 // -- Unused:
53 virtual const G4VBiasingInteractionLaw*
54 ProvideOccurenceBiasingInteractionLaw(const G4BiasingProcessInterface*, G4ForceCondition&) final
55 {
56 return 0;
57 }
58 virtual G4VParticleChange* ApplyFinalStateBiasing(const G4BiasingProcessInterface*,
59 const G4Track*, const G4Step*, G4bool&) final
60 {
61 return 0;
62 }
63
64 // -- Used methods ("non-physics biasing methods"):
65 // ------------------------------------------------
66 // -- Method to return the distance or the condition under which
67 // -- requesting the biasing.
68 // -- Here we use the condition "forced" and the distance returned
69 // -- is "dummy" (DBL_MAX).
70 virtual G4double DistanceToApplyOperation(const G4Track*, G4double,
71 G4ForceCondition* condition) final;
72 // -- Method the generate the final state, which is:
73 // -- - made of the primary with half of its original weight, and a clone of it in
74 // -- case of splitting
75 // -- - the primary with increased weight or the primary killed, in case of killing
76 virtual G4VParticleChange* GenerateBiasingFinalState(const G4Track*, const G4Step*) final;
77
78 // -- Specific to this example:
79 // ----------------------------
80 // -- The parallel world index refers to the index in the list of parallel worlds
81 // -- handled by the G4ParallelGeometriesLimiterProcess intance:
82 void SetParallelWorldIndex(G4int parallelWorldIndex)
83 {
84 fParallelWorldIndex = parallelWorldIndex;
85 }
86 G4int GetParallelWorldIndex() const { return fParallelWorldIndex; }
87 // -- Method to set the biasing shared data, they contain many
88 // -- of the information we need related to parallel geometry:
89 void SetBiasingSharedData(const G4BiasingProcessSharedData* sharedData)
90 {
91 fBiasingSharedData = sharedData;
92 }
93 // -- Set the importance map to be used:
94 void SetImportanceMap(std::map<G4int, G4int>* importanceMap) { fImportanceMap = importanceMap; }
95
96 private:
97 G4int fParallelWorldIndex;
98 const G4BiasingProcessSharedData* fBiasingSharedData;
99 G4TouchableHandle fPreStepTouchableHistory;
100 G4TouchableHandle fPostStepTouchableHistory;
101 G4ParticleChange fParticleChange;
102 G4ParticleChangeForNothing fDummyParticleChange;
103 std::map<G4int, G4int>* fImportanceMap;
104 };
105
106 #endif
107