Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/processes/hadronic/models/parton_string/qgsm/src/G4BaryonSplitter.cc

Version: [ ReleaseNotes ] [ 1.0 ] [ 1.1 ] [ 2.0 ] [ 3.0 ] [ 3.1 ] [ 3.2 ] [ 4.0 ] [ 4.0.p1 ] [ 4.0.p2 ] [ 4.1 ] [ 4.1.p1 ] [ 5.0 ] [ 5.0.p1 ] [ 5.1 ] [ 5.1.p1 ] [ 5.2 ] [ 5.2.p1 ] [ 5.2.p2 ] [ 6.0 ] [ 6.0.p1 ] [ 6.1 ] [ 6.2 ] [ 6.2.p1 ] [ 6.2.p2 ] [ 7.0 ] [ 7.0.p1 ] [ 7.1 ] [ 7.1.p1 ] [ 8.0 ] [ 8.0.p1 ] [ 8.1 ] [ 8.1.p1 ] [ 8.1.p2 ] [ 8.2 ] [ 8.2.p1 ] [ 8.3 ] [ 8.3.p1 ] [ 8.3.p2 ] [ 9.0 ] [ 9.0.p1 ] [ 9.0.p2 ] [ 9.1 ] [ 9.1.p1 ] [ 9.1.p2 ] [ 9.1.p3 ] [ 9.2 ] [ 9.2.p1 ] [ 9.2.p2 ] [ 9.2.p3 ] [ 9.2.p4 ] [ 9.3 ] [ 9.3.p1 ] [ 9.3.p2 ] [ 9.4 ] [ 9.4.p1 ] [ 9.4.p2 ] [ 9.4.p3 ] [ 9.4.p4 ] [ 9.5 ] [ 9.5.p1 ] [ 9.5.p2 ] [ 9.6 ] [ 9.6.p1 ] [ 9.6.p2 ] [ 9.6.p3 ] [ 9.6.p4 ] [ 10.0 ] [ 10.0.p1 ] [ 10.0.p2 ] [ 10.0.p3 ] [ 10.0.p4 ] [ 10.1 ] [ 10.1.p1 ] [ 10.1.p2 ] [ 10.1.p3 ] [ 10.2 ] [ 10.2.p1 ] [ 10.2.p2 ] [ 10.2.p3 ] [ 10.3 ] [ 10.3.p1 ] [ 10.3.p2 ] [ 10.3.p3 ] [ 10.4 ] [ 10.4.p1 ] [ 10.4.p2 ] [ 10.4.p3 ] [ 10.5 ] [ 10.5.p1 ] [ 10.6 ] [ 10.6.p1 ] [ 10.6.p2 ] [ 10.6.p3 ] [ 10.7 ] [ 10.7.p1 ] [ 10.7.p2 ] [ 10.7.p3 ] [ 10.7.p4 ] [ 11.0 ] [ 11.0.p1 ] [ 11.0.p2 ] [ 11.0.p3, ] [ 11.0.p4 ] [ 11.1 ] [ 11.1.1 ] [ 11.1.2 ] [ 11.1.3 ] [ 11.2 ] [ 11.2.1 ] [ 11.2.2 ] [ 11.3.0 ]

Diff markup

Differences between /processes/hadronic/models/parton_string/qgsm/src/G4BaryonSplitter.cc (Version 11.3.0) and /processes/hadronic/models/parton_string/qgsm/src/G4BaryonSplitter.cc (Version 6.0.p1)


  1 //                                                  1 //
  2 // *******************************************      2 // ********************************************************************
  3 // * License and Disclaimer                    <<   3 // * DISCLAIMER                                                       *
  4 // *                                                4 // *                                                                  *
  5 // * The  Geant4 software  is  copyright of th <<   5 // * The following disclaimer summarizes all the specific disclaimers *
  6 // * the Geant4 Collaboration.  It is provided <<   6 // * of contributors to this software. The specific disclaimers,which *
  7 // * conditions of the Geant4 Software License <<   7 // * govern, are listed with their locations in:                      *
  8 // * LICENSE and available at  http://cern.ch/ <<   8 // *   http://cern.ch/geant4/license                                  *
  9 // * include a list of copyright holders.      << 
 10 // *                                                9 // *                                                                  *
 11 // * Neither the authors of this software syst     10 // * Neither the authors of this software system, nor their employing *
 12 // * institutes,nor the agencies providing fin     11 // * institutes,nor the agencies providing financial support for this *
 13 // * work  make  any representation or  warran     12 // * work  make  any representation or  warranty, express or implied, *
 14 // * regarding  this  software system or assum     13 // * regarding  this  software system or assume any liability for its *
 15 // * use.  Please see the license in the file  <<  14 // * use.                                                             *
 16 // * for the full disclaimer and the limitatio << 
 17 // *                                               15 // *                                                                  *
 18 // * This  code  implementation is the result  <<  16 // * This  code  implementation is the  intellectual property  of the *
 19 // * technical work of the GEANT4 collaboratio <<  17 // * GEANT4 collaboration.                                            *
 20 // * By using,  copying,  modifying or  distri <<  18 // * By copying,  distributing  or modifying the Program (or any work *
 21 // * any work based  on the software)  you  ag <<  19 // * based  on  the Program)  you indicate  your  acceptance of  this *
 22 // * use  in  resulting  scientific  publicati <<  20 // * statement, and all its terms.                                    *
 23 // * acceptance of all terms of the Geant4 Sof << 
 24 // *******************************************     21 // ********************************************************************
 25 //                                                 22 //
 26 // Split baryon (antibaryon) into quark and di <<  23 // Split barion (antibarion) into quark and diquark (antidiquark and antiqaurk ) 
 27 // based on prototype, needs clean up of inter     24 // based on prototype, needs clean up of interfaces HPW Feb 1999  
 28 // Numbers verified and errors corrected, HPW      25 // Numbers verified and errors corrected, HPW Dec 1999
 29                                                    26 
 30 #include "G4BaryonSplitter.hh"                     27 #include "G4BaryonSplitter.hh"
 31 #include "G4ParticleTable.hh"                      28 #include "G4ParticleTable.hh"
 32                                                    29 
 33 G4BaryonSplitter::                                 30 G4BaryonSplitter::
 34 G4BaryonSplitter()                                 31 G4BaryonSplitter()
 35 {                                                  32 {
 36   theBaryons.insert(new G4SPBaryon(G4Proton::P     33   theBaryons.insert(new G4SPBaryon(G4Proton::Proton()));
 37   theBaryons.insert(new G4SPBaryon(G4Neutron::     34   theBaryons.insert(new G4SPBaryon(G4Neutron::Neutron()));
 38   theBaryons.insert(new G4SPBaryon(G4AntiProto     35   theBaryons.insert(new G4SPBaryon(G4AntiProton::AntiProton()));
 39   theBaryons.insert(new G4SPBaryon(G4AntiNeutr     36   theBaryons.insert(new G4SPBaryon(G4AntiNeutron::AntiNeutron()));
 40   theBaryons.insert(new G4SPBaryon(G4ParticleT <<  37   theBaryons.insert(new G4SPBaryon(G4ParticleTable::GetParticleTable()->FindParticle(2224))); // D++
 41   theBaryons.insert(new G4SPBaryon(G4ParticleT <<  38   theBaryons.insert(new G4SPBaryon(G4ParticleTable::GetParticleTable()->FindParticle(2214))); // D+
 42   theBaryons.insert(new G4SPBaryon(G4ParticleT <<  39   theBaryons.insert(new G4SPBaryon(G4ParticleTable::GetParticleTable()->FindParticle(2114))); // D0
 43   theBaryons.insert(new G4SPBaryon(G4ParticleT <<  40   theBaryons.insert(new G4SPBaryon(G4ParticleTable::GetParticleTable()->FindParticle(1114))); // D-
 44   theBaryons.insert(new G4SPBaryon(G4ParticleT <<  41   theBaryons.insert(new G4SPBaryon(G4ParticleTable::GetParticleTable()->FindParticle(-2224))); // anti D++
 45   theBaryons.insert(new G4SPBaryon(G4ParticleT <<  42   theBaryons.insert(new G4SPBaryon(G4ParticleTable::GetParticleTable()->FindParticle(-2214))); // anti D+
 46   theBaryons.insert(new G4SPBaryon(G4ParticleT <<  43   theBaryons.insert(new G4SPBaryon(G4ParticleTable::GetParticleTable()->FindParticle(-2114))); // anti D0
 47   theBaryons.insert(new G4SPBaryon(G4ParticleT <<  44   theBaryons.insert(new G4SPBaryon(G4ParticleTable::GetParticleTable()->FindParticle(-1114))); // anti D-
 48   theBaryons.insert(new G4SPBaryon(G4Lambda::L     45   theBaryons.insert(new G4SPBaryon(G4Lambda::Lambda()));
 49   theBaryons.insert(new G4SPBaryon(G4AntiLambd     46   theBaryons.insert(new G4SPBaryon(G4AntiLambda::AntiLambda()));
 50   theBaryons.insert(new G4SPBaryon(G4SigmaPlus     47   theBaryons.insert(new G4SPBaryon(G4SigmaPlus::SigmaPlus()));
 51   theBaryons.insert(new G4SPBaryon(G4SigmaZero     48   theBaryons.insert(new G4SPBaryon(G4SigmaZero::SigmaZero()));
 52   theBaryons.insert(new G4SPBaryon(G4SigmaMinu     49   theBaryons.insert(new G4SPBaryon(G4SigmaMinus::SigmaMinus()));
 53   theBaryons.insert(new G4SPBaryon(G4AntiSigma     50   theBaryons.insert(new G4SPBaryon(G4AntiSigmaPlus::AntiSigmaPlus()));
 54   theBaryons.insert(new G4SPBaryon(G4AntiSigma     51   theBaryons.insert(new G4SPBaryon(G4AntiSigmaZero::AntiSigmaZero()));
 55   theBaryons.insert(new G4SPBaryon(G4AntiSigma     52   theBaryons.insert(new G4SPBaryon(G4AntiSigmaMinus::AntiSigmaMinus()));
                                                   >>  53   theBaryons.insert(new G4SPBaryon(G4ParticleTable::GetParticleTable()->FindParticle(3224))); // S+*
                                                   >>  54   theBaryons.insert(new G4SPBaryon(G4ParticleTable::GetParticleTable()->FindParticle(3214))); // S0*
                                                   >>  55   theBaryons.insert(new G4SPBaryon(G4ParticleTable::GetParticleTable()->FindParticle(3114))); // S-*
                                                   >>  56   theBaryons.insert(new G4SPBaryon(G4ParticleTable::GetParticleTable()->FindParticle(-3224))); // anti S+*
                                                   >>  57   theBaryons.insert(new G4SPBaryon(G4ParticleTable::GetParticleTable()->FindParticle(-3214))); // anti S0*
                                                   >>  58   theBaryons.insert(new G4SPBaryon(G4ParticleTable::GetParticleTable()->FindParticle(-3114))); // anti S-*
 56   theBaryons.insert(new G4SPBaryon(G4XiMinus::     59   theBaryons.insert(new G4SPBaryon(G4XiMinus::XiMinus()));
 57   theBaryons.insert(new G4SPBaryon(G4XiZero::X     60   theBaryons.insert(new G4SPBaryon(G4XiZero::XiZero()));
 58   theBaryons.insert(new G4SPBaryon(G4AntiXiMin     61   theBaryons.insert(new G4SPBaryon(G4AntiXiMinus::AntiXiMinus()));
 59   theBaryons.insert(new G4SPBaryon(G4AntiXiZer     62   theBaryons.insert(new G4SPBaryon(G4AntiXiZero::AntiXiZero()));
 60   theBaryons.insert(new G4SPBaryon(G4OmegaMinu     63   theBaryons.insert(new G4SPBaryon(G4OmegaMinus::OmegaMinus()));
 61   theBaryons.insert(new G4SPBaryon(G4AntiOmega     64   theBaryons.insert(new G4SPBaryon(G4AntiOmegaMinus::AntiOmegaMinus()));
 62   theBaryons.insert(new G4SPBaryon(G4LambdacPl <<  65   theBaryons.insert(new G4SPBaryon(G4ParticleTable::GetParticleTable()->FindParticle(3324))); // X0*
 63   theBaryons.insert(new G4SPBaryon(G4AntiLambd <<  66   theBaryons.insert(new G4SPBaryon(G4ParticleTable::GetParticleTable()->FindParticle(3314))); // X-*
 64   theBaryons.insert(new G4SPBaryon(G4SigmacPlu <<  67   theBaryons.insert(new G4SPBaryon(G4ParticleTable::GetParticleTable()->FindParticle(-3324))); // anti X0*
 65   theBaryons.insert(new G4SPBaryon(G4AntiSigma <<  68   theBaryons.insert(new G4SPBaryon(G4ParticleTable::GetParticleTable()->FindParticle(-3314))); // anti X-*
 66   theBaryons.insert(new G4SPBaryon(G4SigmacPlu << 
 67   theBaryons.insert(new G4SPBaryon(G4AntiSigma << 
 68   theBaryons.insert(new G4SPBaryon(G4SigmacZer << 
 69   theBaryons.insert(new G4SPBaryon(G4AntiSigma << 
 70   theBaryons.insert(new G4SPBaryon(G4XicPlus:: << 
 71   theBaryons.insert(new G4SPBaryon(G4AntiXicPl << 
 72   theBaryons.insert(new G4SPBaryon(G4XicZero:: << 
 73   theBaryons.insert(new G4SPBaryon(G4AntiXicZe << 
 74   theBaryons.insert(new G4SPBaryon(G4OmegacZer << 
 75   theBaryons.insert(new G4SPBaryon(G4AntiOmega << 
 76   theBaryons.insert(new G4SPBaryon(G4Lambdab:: << 
 77   theBaryons.insert(new G4SPBaryon(G4AntiLambd << 
 78   theBaryons.insert(new G4SPBaryon(G4SigmabPlu << 
 79   theBaryons.insert(new G4SPBaryon(G4AntiSigma << 
 80   theBaryons.insert(new G4SPBaryon(G4SigmabZer << 
 81   theBaryons.insert(new G4SPBaryon(G4AntiSigma << 
 82   theBaryons.insert(new G4SPBaryon(G4SigmabMin << 
 83   theBaryons.insert(new G4SPBaryon(G4AntiSigma << 
 84   theBaryons.insert(new G4SPBaryon(G4XibZero:: << 
 85   theBaryons.insert(new G4SPBaryon(G4AntiXibZe << 
 86   theBaryons.insert(new G4SPBaryon(G4XibMinus: << 
 87   theBaryons.insert(new G4SPBaryon(G4AntiXibMi << 
 88   theBaryons.insert(new G4SPBaryon(G4OmegabMin << 
 89   theBaryons.insert(new G4SPBaryon(G4AntiOmega << 
 90 }                                                  69 }
 91                                                    70 
 92 G4bool G4BaryonSplitter::                          71 G4bool G4BaryonSplitter::
 93 SplitBarion(G4int PDGCode, G4int& q_or_qqbar,  <<  72 SplitBarion(G4int PDGCode, G4int* q_or_qqbar, G4int* qbar_or_qq)
 94 {                                                  73 {
 95   const G4SPBaryon * aBaryon = theBaryons.GetB     74   const G4SPBaryon * aBaryon = theBaryons.GetBaryon(G4ParticleTable::GetParticleTable()->FindParticle(PDGCode));
 96                                                <<  75   if(aBaryon==NULL)
 97   if(aBaryon == nullptr)                       <<  76   {
                                                   >>  77     return FALSE;
                                                   >>  78   }
                                                   >>  79   else
 98   {                                                80   {
 99     return false;                              <<  81     aBaryon->SampleQuarkAndDiquark(*q_or_qqbar, *qbar_or_qq);
100   } else {                                     <<  82     return TRUE;
101     aBaryon->SampleQuarkAndDiquark(q_or_qqbar, << 
102     return true;                               << 
103   }                                                83   }
104 }                                                  84 }
105                                                    85 
106                                                    86 
107 // Get the splittable baryon for a PDG code.       87 // Get the splittable baryon for a PDG code.
108 const G4SPBaryon & G4BaryonSplitter::              88 const G4SPBaryon & G4BaryonSplitter::
109 GetSPBaryon(G4int PDGCode)                         89 GetSPBaryon(G4int PDGCode)
110 {                                                  90 {
111   return *theBaryons.GetBaryon(G4ParticleTable <<  91 return *theBaryons.GetBaryon(G4ParticleTable::GetParticleTable()->FindParticle(PDGCode));
112 }                                                  92 }
113                                                    93 
114                                                    94 
115 // Find rest diquark in given barion after qua     95 // Find rest diquark in given barion after quark - antiquark annihilation  
116 G4bool G4BaryonSplitter::FindDiquark(G4int PDG <<  96 G4bool G4BaryonSplitter::
                                                   >>  97 FindDiquark(G4int PDGCode, G4int Quark, G4int* Diquark)
117 {                                                  98 {
118   const G4SPBaryon * aBaryon = theBaryons.GetB     99   const G4SPBaryon * aBaryon = theBaryons.GetBaryon(G4ParticleTable::GetParticleTable()->FindParticle(PDGCode));
119   if(aBaryon)                                  << 100   if(aBaryon==NULL)
                                                   >> 101   {
                                                   >> 102     return FALSE;
                                                   >> 103   }
                                                   >> 104   else
120   {                                               105   {
121     aBaryon->FindDiquark(Quark, Diquark);      << 106     aBaryon->FindDiquark(Quark, *Diquark);
122     return true;                               << 107     return TRUE;
123   }                                               108   }
124   return false;                                << 
125 }                                                 109 }
126                                                   110 
127                                                   111