Geant4  10.02
G4UnstableFermiFragment.cc
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: G4UnstableFermiFragment.cc 85607 2014-10-31 11:21:24Z gcosmo $
27 //
28 // Hadronic Process: Nuclear De-excitations
29 // by V. Lara (Nov 1998)
30 //
31 // Modifications:
32 // 01.04.2011 General cleanup by V.Ivanchenko: integer Z and A, constructor
33 
35 #include "G4FermiFragmentsPool.hh"
36 
38  G4int Pol, G4double ExE)
39  : G4VFermiFragment(anA,aZ,Pol,ExE)
40 {
41  isStable = false;
42 }
43 
45 {}
46 
48  const G4LorentzVector& aMomentum) const
49 {
50  const G4FermiPhaseSpaceDecay* thePhaseSpace
52 
53  std::vector<G4LorentzVector*> * P = thePhaseSpace->Decay(aMomentum.m(), Masses);
54 
55  G4ThreeVector Beta = aMomentum.boostVector();
56 
57  size_t N = P->size();
58 
59  for (size_t i=0; i<N; ++i)
60  {
61  G4LorentzVector* v = (*P)[i];
62  v->boost(Beta);
63  theResult->push_back(new G4Fragment(AtomNum[i],Charges[i],*v));
64 
65  delete v;
66  }
67  delete P;
68 }
std::vector< G4double > Masses
CLHEP::Hep3Vector G4ThreeVector
int G4int
Definition: G4Types.hh:78
static double P[]
G4UnstableFermiFragment(G4int anA, G4int aZ, G4int Pol, G4double ExE)
std::vector< G4Fragment * > G4FragmentVector
Definition: G4Fragment.hh:63
void FillFragment(G4FragmentVector *, const G4LorentzVector &aMomentum) const
std::vector< G4LorentzVector * > * Decay(G4double parent_mass, const std::vector< G4double > &fragment_masses) const
double G4double
Definition: G4Types.hh:76
static G4FermiFragmentsPool * Instance()
const G4FermiPhaseSpaceDecay * GetFermiPhaseSpaceDecay() const
CLHEP::HepLorentzVector G4LorentzVector