G4RPGReaction.hh

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 // $Id$
00027 //
00028 // Author: D. H. Wright
00029 // Date:   26 May 2007
00030 //
00031 
00032 #ifndef G4RPGReaction_h
00033 #define G4RPGReaction_h 1
00034 
00035 // Class Description:
00036 //
00037 // Base class providing methods for various stages of the re-parameterized 
00038 // Gheisha model calculation of primary and secondary final state momenta
00039 
00040 #include "G4DynamicParticle.hh"
00041 #include "G4ReactionProduct.hh"
00042 #include "G4Nucleus.hh"
00043 #include "G4FastVector.hh"
00044 #include "G4HadProjectile.hh"
00045 
00046 
00047 class G4RPGReaction 
00048 {
00049 public:  // with description
00050     
00051   G4RPGReaction() {}
00052     
00053   virtual ~G4RPGReaction() {}
00054 
00055   G4bool ReactionStage(const G4HadProjectile* /*originalIncident*/,
00056                        G4ReactionProduct& /*modifiedOriginal*/,
00057                        G4bool& /*incidentHasChanged*/,
00058                        const G4DynamicParticle* /*originalTarget*/,
00059                        G4ReactionProduct& /*targetParticle*/,
00060                        G4bool& /*targetHasChanged*/,
00061                        const G4Nucleus& /*targetNucleus*/,
00062                        G4ReactionProduct& /*currentParticle*/,
00063                        G4FastVector<G4ReactionProduct,256>& /*vec*/,
00064                        G4int& /*vecLen*/,
00065                        G4bool /*leadFlag*/,
00066                        G4ReactionProduct& /*leadingStrangeParticle*/);
00067 
00068 
00069   void AddBlackTrackParticles(const G4double /*epnb*/,
00070                               const G4int /*npnb*/,
00071                               const G4double /*edta*/,
00072                               const G4int /*ndta*/,
00073                               const G4ReactionProduct& /*modifiedOriginal*/,
00074                               G4int /*PinNucleus*/,
00075                               G4int /*NinNucleus*/,
00076                               const G4Nucleus& /*aNucleus*/,
00077                               G4FastVector<G4ReactionProduct,256>& /*vec*/,
00078                               G4int& /*vecLen*/ );
00079 
00080 
00081   G4double GenerateNBodyEvent(const G4double totalEnergy,
00082                               const G4bool constantCrossSection,
00083                               G4FastVector<G4ReactionProduct,256> &vec,
00084                               G4int& vecLen);
00085     
00086   G4double GenerateNBodyEventT(const G4double totalEnergy,
00087                                const G4bool constantCrossSection,
00088                                std::vector<G4ReactionProduct*>& list);
00089     
00090   void NuclearReaction(G4FastVector<G4ReactionProduct,4> &vec,
00091                        G4int& vecLen,
00092                        const G4HadProjectile* originalIncident,
00093                        const G4Nucleus& aNucleus,
00094                        const G4double theAtomicMass,
00095                        const G4double* massVec);
00096     
00097 protected:  // with description
00098     
00099   void Rotate(const G4double numberofFinalStateNucleons,
00100               const G4ThreeVector& temp,
00101               const G4ReactionProduct& modifiedOriginal,
00102               const G4HadProjectile* originalIncident,
00103               const G4Nucleus& targetNucleus,
00104               G4ReactionProduct &currentParticle,
00105               G4ReactionProduct &targetParticle,
00106               G4FastVector<G4ReactionProduct,256> &vec,
00107               G4int& vecLen );
00108     
00109   void Defs1(const G4ReactionProduct& modifiedOriginal,
00110              G4ReactionProduct& currentParticle,
00111              G4ReactionProduct& targetParticle,
00112              G4FastVector<G4ReactionProduct,256> &vec,
00113              G4int& vecLen);
00114     
00115   std::pair<G4int, G4int> GetFinalStateNucleons(
00116              const G4DynamicParticle* originalTarget,
00117              const G4FastVector<G4ReactionProduct,256>& vec,
00118              const G4int& vecLen );
00119 
00120   void MomentumCheck(const G4ReactionProduct &modifiedOriginal,
00121                      G4ReactionProduct &currentParticle,
00122                      G4ReactionProduct &targetParticle,
00123                      G4FastVector<G4ReactionProduct,256> &vec,
00124                      G4int& vecLen);
00125     
00126   G4double normal();
00127     
00128   G4ThreeVector Isotropic(const G4double&);
00129  
00130 };
00131  
00132 #endif
00133  

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