Geant4  10.03.p03
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
G4eeTo3PiModel Class Reference

#include <G4eeTo3PiModel.hh>

Inheritance diagram for G4eeTo3PiModel:
Collaboration diagram for G4eeTo3PiModel:

Public Member Functions

 G4eeTo3PiModel (G4eeCrossSections *, G4double, G4double)
 
virtual ~G4eeTo3PiModel ()
 
virtual G4double PeakEnergy () const override
 
virtual G4double ComputeCrossSection (G4double) const override
 
virtual void SampleSecondaries (std::vector< G4DynamicParticle * > *, G4double, const G4ThreeVector &) override
 
- Public Member Functions inherited from G4Vee2hadrons
 G4Vee2hadrons (G4eeCrossSections *cr, G4double vlowEnergy, G4double vhighEnergy, G4double vdelta)
 
virtual ~G4Vee2hadrons ()
 
G4PhysicsVectorPhysicsVector () const
 
G4double LowEnergy () const
 
G4double HighEnergy () const
 

Additional Inherited Members

- Protected Attributes inherited from G4Vee2hadrons
G4eeCrossSectionscross
 

Detailed Description

Definition at line 59 of file G4eeTo3PiModel.hh.

Constructor & Destructor Documentation

G4eeTo3PiModel::G4eeTo3PiModel ( G4eeCrossSections cr,
G4double  maxkinEnergy,
G4double  binWidth 
)
explicit

Definition at line 66 of file G4eeTo3PiModel.cc.

69 : G4Vee2hadrons(cr,
70  0.41612*GeV, //threshold
71  maxkinEnergy,
72  binWidth)
73 {
74  G4cout << "####G4eeTo3PiModel####" << G4endl;
75 
76  massPi = G4PionPlus::PionPlus()->GetPDGMass();
77  massPi0 = G4PionZero::PionZero()->GetPDGMass();
78  massOm = 782.62*MeV;
79  massPhi = 1019.46*MeV;
80  gmax = 3.0e-8;
81 }
G4GLOB_DLL std::ostream G4cout
static G4PionPlus * PionPlus()
Definition: G4PionPlus.cc:98
static G4PionZero * PionZero()
Definition: G4PionZero.cc:108
G4Vee2hadrons(G4eeCrossSections *cr, G4double vlowEnergy, G4double vhighEnergy, G4double vdelta)
G4double GetPDGMass() const
static constexpr double GeV
Definition: G4SIunits.hh:217
#define G4endl
Definition: G4ios.hh:61
static constexpr double MeV
Definition: G4SIunits.hh:214

Here is the call graph for this function:

G4eeTo3PiModel::~G4eeTo3PiModel ( )
virtual

Definition at line 85 of file G4eeTo3PiModel.cc.

86 {}

Member Function Documentation

G4double G4eeTo3PiModel::ComputeCrossSection ( G4double  e) const
overridevirtual

Implements G4Vee2hadrons.

Definition at line 99 of file G4eeTo3PiModel.cc.

100 {
101  return cross->CrossSection3pi(e);
102 }
G4double CrossSection3pi(G4double)
G4eeCrossSections * cross

Here is the call graph for this function:

G4double G4eeTo3PiModel::PeakEnergy ( ) const
overridevirtual

Implements G4Vee2hadrons.

Definition at line 90 of file G4eeTo3PiModel.cc.

91 {
92  G4double e = massOm;
93  if(HighEnergy() > massPhi) { e = massPhi; }
94  return e;
95 }
G4double HighEnergy() const
double G4double
Definition: G4Types.hh:76

Here is the call graph for this function:

void G4eeTo3PiModel::SampleSecondaries ( std::vector< G4DynamicParticle * > *  newp,
G4double  e,
const G4ThreeVector direction 
)
overridevirtual

Implements G4Vee2hadrons.

Definition at line 106 of file G4eeTo3PiModel.cc.

108 {
109 
110  G4double x0 = massPi0/e;
111  G4double x1 = massPi/e;
112 
113  G4LorentzVector w0(0.,0.,0.,0.), w1(0.,0.,0.,0.), w2(0.,0.,0.,0.);
114  G4ThreeVector dir0, dir1;
115  G4double e0, p0, e2, p, gg, m01, m02, m12;
116 
117  // max pi0 energy
118  G4double edel = 0.5*e*(1.0 + x0*x0 - 4.0*x1*x1) - massPi0;
119 
120  const G4int nmax = 200;
121  G4int nn = 0;
122  do {
123  ++nn;
124  // pi0 sample
125  e0 = edel*G4UniformRand() + massPi0;
126  p0 = sqrt(e0*e0 - massPi0*massPi0);
127  dir0 = G4RandomDirection();
128  w0 = G4LorentzVector(p0*dir0.x(),p0*dir0.y(),p0*dir0.z(),e0);
129 
130  // pi+pi- pair
131  w1 = G4LorentzVector(-p0*dir0.x(),-p0*dir0.y(),-p0*dir0.z(),e-e0);
132  G4ThreeVector bst = w1.boostVector();
133  e2 = 0.25*w1.m2();
134 
135  // pi+
136  p = sqrt(e2 - massPi*massPi);
137  dir1 = G4RandomDirection();
138  w2 = G4LorentzVector(p*dir1.x(),p*dir1.y(),p*dir1.z(),sqrt(e2));
139  // pi-
140  w1.set(-w2.px(), -w2.py(), -w2.pz(), w2.e());
141 
142  w1.boost(bst);
143  w2.boost(bst);
144 
145  G4double px2 = w2.x();
146  G4double py2 = w2.y();
147  G4double pz2 = w2.z();
148 
149  G4double px1 = w1.x();
150  G4double py1 = w1.y();
151  G4double pz1 = w1.z();
152 
153  m01 = w0*w1;
154  m02 = w0*w2;
155  m12 = w1*w2;
156 
157  G4double px = py1*pz2 - py2*pz1;
158  G4double py = pz1*px2 - pz2*px1;
159  G4double pz = px1*py2 - px2*py1;
160 
161  gg = (px*px + py*py + pz*pz)*
162  norm( 1.0/cross->DpRho(m01) + 1.0/cross->DpRho(m02)
163  + 1.0/cross->DpRho(m12) );
164 
165  if(gg > gmax) {
166  G4cout << "G4eeTo3PiModel::SampleSecondaries WARNING matrix element g= "
167  << gg << " > " << gmax << " (majoranta)" << G4endl;
168  gmax = gg;
169  }
170  // Loop checking, 07-Aug-2015, Vladimir Ivanchenko
171  } while( gmax*G4UniformRand() > gg || nn < nmax);
172 
173  w0.rotateUz(direction);
174  w1.rotateUz(direction);
175  w2.rotateUz(direction);
176 
177  // create G4DynamicParticle objects
178  G4DynamicParticle* dp0 =
180  G4DynamicParticle* dp1 =
182  G4DynamicParticle* dp2 =
184  newp->push_back(dp0);
185  newp->push_back(dp1);
186  newp->push_back(dp2);
187 }
Hep3Vector boostVector() const
double x() const
const char * p
Definition: xmltok.h:285
G4ThreeVector G4RandomDirection()
int G4int
Definition: G4Types.hh:78
double z() const
std::complex< G4double > DpRho(G4double e)
#define G4UniformRand()
Definition: Randomize.hh:97
G4GLOB_DLL std::ostream G4cout
HepLorentzVector & boost(double, double, double)
static G4PionPlus * PionPlus()
Definition: G4PionPlus.cc:98
const G4int nmax
static G4PionZero * PionZero()
Definition: G4PionZero.cc:108
void set(double x, double y, double z, double t)
double m2() const
static G4PionMinus * PionMinus()
Definition: G4PionMinus.cc:98
double y() const
G4eeCrossSections * cross
#define G4endl
Definition: G4ios.hh:61
double G4double
Definition: G4Types.hh:76
CLHEP::HepLorentzVector G4LorentzVector

Here is the call graph for this function:


The documentation for this class was generated from the following files: