G4PiMinusStopPb.cc

Go to the documentation of this file.
00001 //
00002 // ********************************************************************
00003 // * License and Disclaimer                                           *
00004 // *                                                                  *
00005 // * The  Geant4 software  is  copyright of the Copyright Holders  of *
00006 // * the Geant4 Collaboration.  It is provided  under  the terms  and *
00007 // * conditions of the Geant4 Software License,  included in the file *
00008 // * LICENSE and available at  http://cern.ch/geant4/license .  These *
00009 // * include a list of copyright holders.                             *
00010 // *                                                                  *
00011 // * Neither the authors of this software system, nor their employing *
00012 // * institutes,nor the agencies providing financial support for this *
00013 // * work  make  any representation or  warranty, express or implied, *
00014 // * regarding  this  software system or assume any liability for its *
00015 // * use.  Please see the license in the file  LICENSE  and URL above *
00016 // * for the full disclaimer and the limitation of liability.         *
00017 // *                                                                  *
00018 // * This  code  implementation is the result of  the  scientific and *
00019 // * technical work of the GEANT4 collaboration.                      *
00020 // * By using,  copying,  modifying or  distributing the software (or *
00021 // * any work based  on the software)  you  agree  to acknowledge its *
00022 // * use  in  resulting  scientific  publications,  and indicate your *
00023 // * acceptance of all terms of the Geant4 Software license.          *
00024 // ********************************************************************
00025 //
00026 //      File name:     G4PiMinusStopPb
00027 //
00028 //      Author:        Maria Grazia Pia (pia@genova.infn.it)
00029 // 
00030 //      Creation date: 8 May 1998
00031 //
00032 // -------------------------------------------------------------------
00033 
00034 #include "G4ios.hh"
00035 
00036 #include "G4PiMinusStopPb.hh"
00037 
00038 #include <vector>
00039 
00040 #include "globals.hh"
00041 #include "Randomize.hh"
00042 #include "G4Proton.hh"
00043 #include "G4Neutron.hh"
00044 #include "G4ParticleTypes.hh"
00045 #include "G4ReactionKinematics.hh"
00046 #include "G4DynamicParticleVector.hh"
00047 #include "G4LorentzVector.hh"
00048 #include "G4PiMinusStopMaterial.hh"
00049 #include "G4DistributionGenerator.hh"
00050 
00051 // np/pp production ratio
00052 // Experimental values: 
00053 // R(np/pp) (E. Gadioli et al., Phys Rev C36 (1987) 741) 
00054 G4double G4PiMinusStopPb::npRatio = 6;
00055  
00056 // Average numbers of final nucleons detected, for N-pair absorption
00057 // (R. Madey  et al., Phys Rev C25 (1982) 3050
00058 G4double G4PiMinusStopPb::nFinalNucleons = 1.67;
00059 
00060 // Kinetic energy (MeV) distributions measured for coincident nucleon 
00061 // emission
00062 // (R. Madey  et al., Phys Rev C25 (1982) 3050
00063 
00064 G4int G4PiMinusStopPb::eKinEntries = 10;
00065 
00066 G4double G4PiMinusStopPb::eKinData[10] = { 0.24, 0.5, 1.03, 1.28,
00067                                      1.44, 
00068                                      1.28, 1.03, 0.82, 0.50, 0.22};
00069 
00070 G4double G4PiMinusStopPb::eKin[11] = { 5.2, 15., 27., 41.5,
00071                                   49.6,
00072                                   57.7, 79.3, 94.4, 114., 125., 140.};
00073 
00074 
00075 // Opening angle distributions measured for coincident nucleon emission
00076 // (P.Heusi et al., Nucl. Phys. A407 (1983) 429
00077 
00078 G4int G4PiMinusStopPb::angleEntries = 7;
00079 
00080 G4double G4PiMinusStopPb::angleData[7] = 
00081 { 1.43, 1.67, 2.62, 4.29, 7.62, 11.90, 14.76 };
00082 
00083 G4double G4PiMinusStopPb::angle[8] = { 1.308997, 1.570796, 1.832596, 2.094395, 
00084                                   2.356194, 2.617994, 2.967060, 3.1415927 };
00085 
00086 
00087 
00088 // Constructor
00089 
00090 G4PiMinusStopPb::G4PiMinusStopPb()
00091   
00092 {
00093   // Cluster size: nucleon pair, alpha, triton etc.
00094   // First implementation: interaction with nucleon pair only
00095   _clusterSize = 2;
00096 
00097   // R ratio
00098   theR = 1. / (1. + npRatio);
00099 
00100   _definitions = new std::vector<G4ParticleDefinition*>();
00101   _momenta = new std::vector<G4LorentzVector*>();
00102 
00103   std::vector<double> eKinVector;
00104   std::vector<double> eKinDataVector;
00105   int i;
00106   for (i=0; i<eKinEntries; i++)
00107     {
00108       eKinVector.push_back(eKin[i]);
00109       eKinDataVector.push_back(eKinData[i]);
00110     }
00111   eKinVector.push_back(eKin[eKinEntries]);
00112   _distributionE = new G4DistributionGenerator(eKinVector,eKinDataVector);
00113 
00114   std::vector<double> angleVector;
00115   std::vector<double> angleDataVector;
00116   for (i=0; i<angleEntries; i++)
00117     {
00118       angleVector.push_back(angle[i]);
00119       angleDataVector.push_back(angleData[i]);
00120     }
00121   angleVector.push_back(angle[angleEntries]);
00122   _distributionAngle = new G4DistributionGenerator(angleVector,angleDataVector);
00123 }
00124 
00125 
00126 // Destructor
00127 
00128 G4PiMinusStopPb::~G4PiMinusStopPb()
00129 {}
00130 
00131 G4double G4PiMinusStopPb::FinalNucleons()
00132 {
00133   return nFinalNucleons;
00134 }
00135 

Generated on Mon May 27 17:49:21 2013 for Geant4 by  doxygen 1.4.7