G4PiMinusStopTa.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:     G4PiMinusStopTa
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 "G4PiMinusStopTa.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) = 6.3 +- 1.4 (R. Madey  et al., Phys Rev C25 (1982) 3050
00054 G4double G4PiMinusStopTa::npRatio = 17.;
00055  
00056 // Average numbers of final nucleons detected, for N-pair absorption
00057 // (R. Madey  et al., Phys Rev C25 (1982) 3050
00058 G4double G4PiMinusStopTa::nFinalNucleons = 1.8;
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 G4PiMinusStopTa::eKinEntries = 10;
00065 
00066 G4double G4PiMinusStopTa::eKinData[10] = { 0.24, 0.59, 1.13, 1.38,
00067                                      1.50, 
00068                                      1.38, 1.13, 0.83, 0.49, 0.22};
00069 
00070 G4double G4PiMinusStopTa::eKin[11] = { 5.2, 12., 25., 41.5,
00071                                   49.6,
00072                                   57.7, 79.3, 94.4, 104., 120., 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 G4PiMinusStopTa::angleEntries = 7;
00079 
00080 G4double G4PiMinusStopTa::angleData[7] = 
00081 { 1.43, 1.67, 2.62, 4.29, 7.62, 11.90, 14.76 };
00082 
00083 G4double G4PiMinusStopTa::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 G4PiMinusStopTa::G4PiMinusStopTa()
00091 {
00092   // Cluster size: nucleon pair, alpha, triton etc.
00093   // First implementation: interaction with nucleon pair only
00094   _clusterSize = 2;
00095 
00096   // R ratio
00097   theR = 1. / (1. + npRatio);
00098 
00099   _definitions = new std::vector<G4ParticleDefinition*>();
00100   _momenta = new std::vector<G4LorentzVector*>();
00101 
00102   std::vector<double> eKinVector;
00103   std::vector<double> eKinDataVector;
00104   int i;
00105   for (i=0; i<eKinEntries; i++)
00106     {
00107       eKinVector.push_back(eKin[i]);
00108       eKinDataVector.push_back(eKinData[i]);
00109     }
00110   eKinVector.push_back(eKin[eKinEntries]);
00111   _distributionE = new G4DistributionGenerator(eKinVector,eKinDataVector);
00112 
00113   std::vector<double> angleVector;
00114   std::vector<double> angleDataVector;
00115   for (i=0; i<angleEntries; i++)
00116     {
00117       angleVector.push_back(angle[i]);
00118       angleDataVector.push_back(angleData[i]);
00119     }
00120   angleVector.push_back(angle[angleEntries]);
00121   _distributionAngle = new G4DistributionGenerator(angleVector,angleDataVector);
00122 }
00123 
00124 
00125 // Destructor
00126 
00127 G4PiMinusStopTa::~G4PiMinusStopTa()
00128 {}
00129 
00130 G4double G4PiMinusStopTa::FinalNucleons()
00131 {
00132   return nFinalNucleons;
00133 }
00134 

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