| 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 /// \file GB05BOptrSplitAndKillByCrossSection.hh
27 /// \brief Definition of the GB05BOptrSplitAndKillByCrossSection class
28 //
29 //---------------------------------------------------------------
30 //
31 // GB05BOptrSplitAndKillByCrossSection
32 //
33 // Class Description:
34 // A G4VBiasingOperator concrete implementation example to
35 // illustrate how to bias physics processes cross-section for
36 // one particle type.
37 // The G4VBiasingOperation G4BOptnChangeCrossSection is
38 // selected by this operator, and is sent to each process
39 // calling the operator.
40 // A simple constant bias to the cross-section is applied,
41 // but more sophisticated changes can be applied.
42 //
43 //---------------------------------------------------------------
44 //
45
46 #ifndef GB05BOptrSplitAndKillByCrossSection_hh
47 #define GB05BOptrSplitAndKillByCrossSection_hh 1
48
49 #include "G4VBiasingOperator.hh"
50 class G4BOptnChangeCrossSection;
51 class G4ParticleDefinition;
52 class G4VProcess;
53 class GB05BOptnSplitAndKillByCrossSection;
54 #include <map>
55
56 class GB05BOptrSplitAndKillByCrossSection : public G4VBiasingOperator
57 {
58 public:
59 // ------------------------------------------------------------
60 // -- Constructor: takes the name of the particle type to bias:
61 // ------------------------------------------------------------
62 GB05BOptrSplitAndKillByCrossSection(G4String particleToBias,
63 G4String name = "SplitAndKillByXS");
64 virtual ~GB05BOptrSplitAndKillByCrossSection();
65
66 // -- method called at beginning of run:
67 virtual void StartRun();
68
69 private:
70 // -----------------------------
71 // -- Mandatory from base class:
72 // -----------------------------
73 // -- Not used:
74 virtual G4VBiasingOperation*
75 ProposeOccurenceBiasingOperation(const G4Track*, const G4BiasingProcessInterface*) final
76 {
77 return 0;
78 }
79
80 // -- Not used:
81 virtual G4VBiasingOperation*
82 ProposeFinalStateBiasingOperation(const G4Track*, const G4BiasingProcessInterface*) final
83 {
84 return 0;
85 }
86
87 // -- Used method : it will return the biasing operation that will split particles
88 // -- with a probabilty depending on the total absorption cross-section.
89 virtual G4VBiasingOperation*
90 ProposeNonPhysicsBiasingOperation(const G4Track*, const G4BiasingProcessInterface*) final;
91
92 // ---------------------------------------
93 // -- Method specific to this application:
94 // ---------------------------------------
95 // -- Each "absorbing" process that the biasing has to counterbalance
96 // -- its action for has to be passed to the biasing operator
97 public:
98 void AddProcessToEquipoise(G4String processName);
99
100 private:
101 GB05BOptnSplitAndKillByCrossSection* fSplitAndKillByCrossSection;
102 const G4ParticleDefinition* fParticleToBias;
103 std::vector<G4String> fProcessesToEquipoise;
104 G4bool fSetup;
105 std::vector<const G4VProcess*> fProcesses;
106 };
107
108 #endif
109