| 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 //
28
29 #ifndef G4FTFParticipants_h
30 #define G4FTFParticipants_h 1
31
32 // ------------------------------------------------------------
33 // GEANT 4 class header file
34 //
35 // ---------------- G4FTFParticipants----------------
36 // by Gunter Folger, June 1998.
37 // class finding colliding particles in FTFPartonStringModel
38 // ------------------------------------------------------------
39
40 #include "G4VParticipants.hh"
41 #include "G4FTFParameters.hh"
42 #include <vector>
43 #include "G4Nucleon.hh"
44 #include "G4V3DNucleus.hh"
45 #include "G4Fancy3DNucleus.hh"
46 #include "G4ReactionProduct.hh"
47 #include "G4InteractionContent.hh"
48
49
50 class G4FTFParticipants : public G4VParticipants {
51 public:
52 G4FTFParticipants();
53 ~G4FTFParticipants();
54
55 const G4FTFParticipants& operator=( const G4FTFParticipants& right ) = delete;
56 G4bool operator==( const G4FTFParticipants& right ) const = delete;
57 G4bool operator!=( const G4FTFParticipants& right ) const = delete;
58 G4FTFParticipants( const G4FTFParticipants& right ) = delete;
59
60 // New methods to get/set the impact parameter.
61 // (Note: to get the impact parameter in fermi units, do:
62 // GetImpactParameter() / fermi )
63 // By default, SampleBinInterval() is false and therefore the sampling of
64 // the impact parameter is done as usual by FTF taking into account the
65 // radius of target nucleus, and, in the case of nucleus-nucleus collisions,
66 // also of radius of the projectile nucleus; if, instead, SampleBinInterval()
67 // is true, then the square of the impact parameter is drawn from a flat
68 // distribution in the specified interval [Bmin2, Bmax2].
69 void SetImpactParameter( const G4double b_value );
70 G4double GetImpactParameter() const;
71 void SetBminBmax( const G4double bmin_value, const G4double bmax_value );
72 G4bool SampleBinInterval() const;
73 G4double GetBmin2() const; // Minimum value of the square of the impact parameter
74 G4double GetBmax2() const; // Maximum value of the square of the impact parameter
75
76 void GetList( const G4ReactionProduct& thePrimary, G4FTFParameters* theParameters );
77 void StartLoop();
78 G4bool Next();
79 void SortInteractionsIncT();
80 void ShiftInteractionTime();
81 G4InteractionContent& GetInteraction();
82 void Clean();
83
84 private:
85 G4double Bimpact;
86 G4bool BinInterval;
87 G4double Bmin2;
88 G4double Bmax2;
89
90 std::vector< G4InteractionContent* > theInteractions;
91 G4int currentInteraction;
92 };
93
94
95 inline void G4FTFParticipants::SetImpactParameter( const G4double b_value ) {
96 Bimpact = b_value;
97 }
98
99 inline G4double G4FTFParticipants::GetImpactParameter() const {
100 return Bimpact;
101 }
102
103 inline void G4FTFParticipants::SetBminBmax( const G4double bmin_value, const G4double bmax_value ) {
104 BinInterval = false;
105 if ( bmin_value < 0.0 || bmax_value < 0.0 || bmax_value < bmin_value ) return;
106 BinInterval = true;
107 Bmin2 = bmin_value * bmin_value;
108 Bmax2 = bmax_value * bmax_value;
109 }
110
111 inline G4bool G4FTFParticipants::SampleBinInterval() const {
112 return BinInterval;
113 }
114
115 inline G4double G4FTFParticipants::GetBmin2() const {
116 return Bmin2;
117 }
118
119 inline G4double G4FTFParticipants::GetBmax2() const {
120 return Bmax2;
121 }
122
123 inline void G4FTFParticipants::StartLoop() {
124 currentInteraction = -1;
125 }
126
127 inline G4bool G4FTFParticipants::Next() {
128 return ++currentInteraction < static_cast< G4int >( theInteractions.size() );
129 }
130
131 inline G4InteractionContent& G4FTFParticipants::GetInteraction() {
132 return *theInteractions[ currentInteraction ];
133 }
134
135 #endif
136
137