Geant4  10.03
G4PolarizedAnnihilationModel.hh
Go to the documentation of this file.
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 // $Id: G4PolarizedAnnihilationModel.hh 96114 2016-03-16 18:51:33Z gcosmo $
27 //
28 // -------------------------------------------------------------------
29 //
30 // GEANT4 Class header file
31 //
32 //
33 // File name: G4PolarizedAnnihilationModel
34 //
35 // Author: Andreas Schaelicke and Pavel Starovoitov
36 //
37 // Creation date: 01.05.2005
38 //
39 // Modifications:
40 // 18-07-06 use newly calculated cross sections (P. Starovoitov)
41 // 21-08-06 update interface to geant4.8.1 (A. Schaelicke)
42 // 10-07-07 copied Initialise() method from G4eeToTwoGammaModel to provide a
43 // ParticleChangeForGamma object (A. Schaelicke)
44 //
45 //
46 // Class Description:
47 //
48 // Implementation of polarized gamma Annihilation scattering on free electron
49 //
50 
51 // -------------------------------------------------------------------
52 //
53 
54 #ifndef G4PolarizedAnnihilationModel_h
55 #define G4PolarizedAnnihilationModel_h 1
56 
57 #include "G4eeToTwoGammaModel.hh"
58 #include "G4StokesVector.hh"
59 
62 
64 {
65 
66 public:
67 
68  explicit G4PolarizedAnnihilationModel(const G4ParticleDefinition* p = nullptr,
69  const G4String& nam = "Polarized-Annihilation");
70 
72 
73  virtual void Initialise(const G4ParticleDefinition*,
74  const G4DataVector&) override;
76  const G4ParticleDefinition*,
77  G4double kinEnergy,
78  G4double cut,
79  G4double emax) override;
80  void ComputeAsymmetriesPerElectron(G4double gammaEnergy,
81  G4double & valueX,
82  G4double & valueA,
83  G4double & valueT);
84 
85  virtual void SampleSecondaries(std::vector<G4DynamicParticle*>*,
86  const G4MaterialCutsCouple*,
87  const G4DynamicParticle*,
88  G4double tmin,
89  G4double maxEnergy) override;
90 
91  // polarized routines
92  inline void SetTargetPolarization(const G4ThreeVector & pTarget);
93  inline void SetBeamPolarization(const G4ThreeVector & pBeam);
94  inline const G4ThreeVector & GetTargetPolarization() const;
95  inline const G4ThreeVector & GetBeamPolarization() const;
96  inline const G4ThreeVector & GetFinalGamma1Polarization() const;
97  inline const G4ThreeVector & GetFinalGamma2Polarization() const;
98 private:
99 
100  // hide assignment operator
104 
106  // incomming
109  // outgoing
112 
114 
117 };
118 
119 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
120 
122 {
123  theTargetPolarization = pTarget;
124 }
126 {
127  theBeamPolarization = pBeam;
128 }
130 {
131  return theTargetPolarization;
132 }
134 {
135  return theBeamPolarization;
136 }
138 {
140 }
142 {
144 }
145 
146 #endif
void SetBeamPolarization(const G4ThreeVector &pBeam)
const G4ThreeVector & GetFinalGamma2Polarization() const
const G4ThreeVector & GetBeamPolarization() const
G4PolarizedAnnihilationModel & operator=(const G4PolarizedAnnihilationModel &right)=delete
CLHEP::Hep3Vector G4ThreeVector
virtual void SampleSecondaries(std::vector< G4DynamicParticle * > *, const G4MaterialCutsCouple *, const G4DynamicParticle *, G4double tmin, G4double maxEnergy) override
virtual G4double ComputeCrossSectionPerElectron(const G4ParticleDefinition *, G4double kinEnergy, G4double cut, G4double emax) override
int G4int
Definition: G4Types.hh:78
G4ParticleChangeForGamma * gParticleChange
void SetTargetPolarization(const G4ThreeVector &pTarget)
bool G4bool
Definition: G4Types.hh:79
const G4ThreeVector & GetFinalGamma1Polarization() const
G4PolarizedAnnihilationCrossSection * crossSectionCalculator
void ComputeAsymmetriesPerElectron(G4double gammaEnergy, G4double &valueX, G4double &valueA, G4double &valueT)
virtual void Initialise(const G4ParticleDefinition *, const G4DataVector &) override
static const G4double emax
G4PolarizedAnnihilationModel(const G4ParticleDefinition *p=nullptr, const G4String &nam="Polarized-Annihilation")
double G4double
Definition: G4Types.hh:76
const G4ThreeVector & GetTargetPolarization() const