ECCE @ EIC Software
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
G4ITReactionChange.hh
Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file G4ITReactionChange.hh
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 // Author: Mathieu Karamitros
28 
29 // The code is developed in the framework of the ESA AO7146
30 //
31 // We would be very happy hearing from you, send us your feedback! :)
32 //
33 // In order for Geant4-DNA to be maintained and still open-source,
34 // article citations are crucial.
35 // If you use Geant4-DNA chemistry and you publish papers about your software,
36 // in addition to the general paper on Geant4-DNA:
37 //
38 // Int. J. Model. Simul. Sci. Comput. 1 (2010) 157–178
39 //
40 // we would be very happy if you could please also cite the following
41 // reference papers on chemistry:
42 //
43 // J. Comput. Phys. 274 (2014) 841-882
44 // Prog. Nucl. Sci. Tec. 2 (2011) 503-508
45 
46 #ifndef G4ITReactionChange_H
47 #define G4ITReactionChange_H
48 
49 #include "globals.hh"
50 #include "G4ParticleChange.hh"
51 #include <vector>
52 
59 {
60 public:
64  virtual ~G4ITReactionChange();
65 
66  //_____________________________________________________
67  // To be used by reaction processes
68  void Initialize(const G4Track&,
69  const G4Track&,
70  G4VParticleChange* particleChangeA = 0,
71  G4VParticleChange* particleChangeB = 0);
72  void AddSecondary(G4Track* aSecondary);
73  inline void KillParents(G4bool);
74 
75  // If both parents are not killed therefore
76  // we can used the particle change
77  // you will have to give the pointers
78  // of your particleChange while intializing
79  // but it might tell you that energy is not
80  // conserved for instance
82 
83  //_____________________________________________________
84  // Not to be used in reaction processes
85  void UpdateStepInfo(G4Step*, G4Step*);
86  G4Track* GetSecondary(G4int) const;
87  std::vector<G4Track*>* GetfSecondary();
88 
90  G4bool WereParentsKilled() const;
91 
92  const G4Track* GetTrackA();
93  const G4Track* GetTrackB();
94 
95 protected:
105  // equal/unequal operator
107  G4bool operator!=(const G4ITReactionChange &right) const;
108  // "equal" means that the objects have the same pointer.
109 protected:
110  std::map<const G4Track*, G4VParticleChange*> fParticleChange;
111  std::vector<G4Track*>* fSecondaries;
115 };
116 
118 {
119  if (fSecondaries) return (*fSecondaries)[anIndex];
120  else return 0;
121 }
122 
124 {
125  return fNumberOfSecondaries;
126 }
127 
129 {
130  fKillParents = kill;
131 }
132 
134 {
135  return fKillParents;
136 }
137 
138 inline std::vector<G4Track*>* G4ITReactionChange::GetfSecondary()
139 {
140  return fSecondaries;
141 }
142 
143 #endif // G4ITReactionChange_H