Geant4 Cross Reference |
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 // ------------------------------------------------------------------- 28 // 29 // GEANT4 Class file 30 // 31 // 32 // File name: G4PreCompoundAlpha 33 // 34 // Author: V.Lara 35 // 36 // Modified: 37 // 21.08.2008 J. M. Quesada add choice of options 38 // 20.08.2010 V.Ivanchenko added G4Pow and G4PreCompoundParameters pointers 39 // use int Z and A and cleanup 40 // 05.07.2013 J.M. Quesada FactorialFactor fixed 41 // 42 43 #include "G4PreCompoundAlpha.hh" 44 #include "G4CoulombBarrier.hh" 45 #include "G4Alpha.hh" 46 47 G4PreCompoundAlpha::G4PreCompoundAlpha() 48 : G4PreCompoundIon(G4Alpha::Alpha(), new G4CoulombBarrier(4, 2)) 49 {} 50 51 G4double G4PreCompoundAlpha::FactorialFactor(G4int N, G4int P) const 52 { 53 return static_cast<G4double>(((N-4)*(P-3)*(N-3)*(P-2))*((N-2)*(P-1)*(N-1)*P)) 54 /144.0; 55 } 56 57 G4double G4PreCompoundAlpha::CoalescenceFactor(G4int A) const 58 { 59 return 4096.0/static_cast<G4double>(A*A*A); 60 } 61 62 G4double G4PreCompoundAlpha::GetAlpha() const 63 { 64 G4double C = 0.0; 65 if (theFragZ <= 30) 66 { 67 C = 0.10; 68 } 69 else if (theFragZ <= 50) 70 { 71 C = 0.1 - (theFragZ-30)*0.001; 72 } 73 else if (theFragZ < 70) 74 { 75 C = 0.08 - (theFragZ-50)*0.001; 76 } 77 else 78 { 79 C = 0.06; 80 } 81 return 1.0+C; 82 } 83 84 G4double G4PreCompoundAlpha::GetRj(G4int nParticles, G4int nCharged) const 85 { 86 G4double rj = 0.0; 87 if(nCharged >=2 && (nParticles-nCharged) >=2 ) { 88 G4double denominator = (G4double) 89 (nParticles*(nParticles-1))*((nParticles-2)*(nParticles-3)); 90 rj = (6.0*nCharged*(nCharged-1))* 91 ((nParticles-nCharged)*(nParticles-nCharged-1))/denominator; 92 } 93 return rj; 94 } 95 96