Changeset 1337 for trunk/source/processes/hadronic/models/chiral_inv_phase_space/processes/src/G4QLowEnergy.cc

Timestamp:

Sep 30, 2010, 2:47:17 PM (14 years ago)

Author:

garnier

Message:

tag geant4.9.4 beta 1 + modifs locales

File:

• trunk/source/processes/hadronic/models/chiral_inv_phase_space/processes/src/G4QLowEnergy.cc (modified) (4 diffs)

trunk/source/processes/hadronic/models/chiral_inv_phase_space/processes/src/G4QLowEnergy.cc

@@ -24 +24 @@
 // ********************************************************************
 //
-// $Id: G4QLowEnergy.cc,v 1.4 2010/06/07 15:19:55 mkossov Exp $
-// GEANT4 tag $Name: geant4-09-04-beta-cand-01 $
+// $Id: G4QLowEnergy.cc,v 1.5 2010/06/14 16:11:27 mkossov Exp $
+// GEANT4 tag $Name: geant4-09-04-beta-01 $
 //
 //      ---------------- G4QLowEnergy class -----------------
@@ -254 +254 @@
   static const G4ThreeVector zeroMom(0.,0.,0.);
   static G4ParticleDefinition* aGamma    = G4Gamma::Gamma();
+  static G4ParticleDefinition* aPiZero   = G4PionZero::PionZero();
+  static G4ParticleDefinition* aPiPlus   = G4PionPlus::PionPlus();
+  static G4ParticleDefinition* aPiMinus  = G4PionMinus::PionMinus();
   static G4ParticleDefinition* aProton   = G4Proton::Proton();
   static G4ParticleDefinition* aNeutron  = G4Neutron::Neutron();
@@ -1853 +1856 @@
     G4QHadron* rHadron = new G4QHadron(90000000+999*rZ+rA,res4Mom); // Input hadron-nucleus
     G4QHadronVector* evaHV = new G4QHadronVector; // Output vector of hadrons (delete!)
-    Nuc.EvaporateNucleus(rHadron, evaHV);
+    Nuc.EvaporateNucleus(rHadron, evaHV); // here a pion can appear !
     G4int nOut=evaHV->size();
     for(G4int i=0; i<nOut; i++)
@@ -1868 +1871 @@
       else if(hPDG==91000000 || hPDG==3122) theDefinition = aLambda;
       else if(hPDG==     22 )               theDefinition = aGamma;
+      else if(hPDG==     111)               theDefinition = aPiZero;
+      else if(hPDG==     211)               theDefinition = aPiPlus;
+      else if(hPDG==    -211)               theDefinition = aPiMinus;
       else
       {