Geant4  10.00.p01
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 67983 2013-03-13 10:42:03Z 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 
37  G4int Pol, G4double ExE)
38  : G4VFermiFragment(anA,aZ,Pol,ExE)
39 {
40  isStable = false;
41 }
42 
44 {}
45 
48 {
49  std::vector<G4LorentzVector*> * P = thePhaseSpace.Decay(aMomentum.m(), Masses);
50 
51  G4ThreeVector Beta = aMomentum.boostVector();
52 
53  G4FragmentVector * theResult = new G4FragmentVector;
54  size_t N = P->size();
55 
56  for (size_t i=0; i<N; ++i)
57  {
58  G4LorentzVector* v = (*P)[i];
59  v->boost(Beta);
60  theResult->push_back(new G4Fragment(AtomNum[i],Charges[i],*v));
61 
62  delete v;
63  }
64  delete P;
65 
66  return theResult;
67 }
std::vector< G4double > Masses
CLHEP::Hep3Vector G4ThreeVector
int G4int
Definition: G4Types.hh:78
std::vector< G4LorentzVector * > * Decay(const G4double, const std::vector< G4double > &) const
G4UnstableFermiFragment(G4int anA, G4int aZ, G4int Pol, G4double ExE)
std::vector< G4Fragment * > G4FragmentVector
Definition: G4Fragment.hh:65
double G4double
Definition: G4Types.hh:76
virtual G4FragmentVector * GetFragment(const G4LorentzVector &) const
G4FermiPhaseSpaceDecay thePhaseSpace
CLHEP::HepLorentzVector G4LorentzVector