G4PiMinusStopCo.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:     G4PiMinusStopCo
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 "G4PiMinusStopCo.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 G4PiMinusStopCo::npRatio = 2.5;
00055 
00056 
00057  
00058 // Average numbers of final nucleons detected, for N-pair absorption
00059 // (Hartmann et al., Nucl. Phys. A300 (1978) 345
00060 G4double G4PiMinusStopCo::nFinalNucleons = 1.38;
00061 
00062 // Kinetic energy (MeV) distributions measured for coincident nucleon 
00063 // emission
00064 // P. Heusi et al., Nucl. Phys. A407(1983) 429
00065 
00066 G4int G4PiMinusStopCo::eKinEntries = 11;
00067 
00068 G4double G4PiMinusStopCo::eKinData[11] = {0.085, 0.09, 0.09, 0.15, 
00069                                       0.1, 0.09, 0.08, 
00070                                       0.04,  0.03,  0.02, 0.01};
00071 
00072 G4double G4PiMinusStopCo::eKin[12] = { 15., 17.5, 25.,  33.,  
00073                                    42.,  52.,  62., 
00074                                    75., 85., 95., 105. };
00075 
00076 
00077 // Opening angle distributions measured for coincident nucleon emission
00078 // (P.Heusi et al., Nucl. Phys. A407 (1983) 429
00079 
00080 G4int G4PiMinusStopCo::angleEntries = 7;
00081 
00082 G4double G4PiMinusStopCo::angleData[7] = 
00083 {6., 8., 9., 10., 25., 40., 45. };
00084 
00085 G4double G4PiMinusStopCo::angle[8] = { 0.5, 0.7, 0.87, 1.4, 2.1, 
00086                                    2.44, 2.8, 3.1415927 };
00087 
00088 
00089 
00090 // Constructor
00091 
00092 G4PiMinusStopCo::G4PiMinusStopCo()
00093   
00094 {
00095   // Cluster size: nucleon pair, alpha, triton etc.
00096   // First implementation: interaction with nucleon pair only
00097   _clusterSize = 2;
00098 
00099   // R ratio
00100   theR = 1. / (1. + npRatio);
00101 
00102   _definitions = new std::vector<G4ParticleDefinition*>();
00103   _momenta = new std::vector<G4LorentzVector*>();
00104 
00105   std::vector<double> eKinVector;
00106   std::vector<double> eKinDataVector;
00107   int i;
00108   for (i=0; i<eKinEntries; i++)
00109     {
00110       eKinVector.push_back(eKin[i]);
00111       eKinDataVector.push_back(eKinData[i]);
00112     }
00113   eKinVector.push_back(eKin[eKinEntries]);
00114   _distributionE = new G4DistributionGenerator(eKinVector,eKinDataVector);
00115 
00116   std::vector<double> angleVector;
00117   std::vector<double> angleDataVector;
00118   for (i=0; i<angleEntries; i++)
00119     {
00120       angleVector.push_back(angle[i]);
00121       angleDataVector.push_back(angleData[i]);
00122     }
00123   angleVector.push_back(angle[angleEntries]);
00124   _distributionAngle = new G4DistributionGenerator(angleVector,angleDataVector);
00125 }
00126 
00127 
00128 // Destructor
00129 
00130 G4PiMinusStopCo::~G4PiMinusStopCo()
00131 {}
00132 
00133 G4double G4PiMinusStopCo::FinalNucleons()
00134 {
00135   return nFinalNucleons;
00136 }
00137 

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