Changeset 1315 for trunk/source/physics_lists/lists/src/HadronPhysicsFTFP_BERT.cc

Timestamp:

Jun 18, 2010, 11:42:07 AM (14 years ago)

Author:

garnier

Message:

update geant4-09-04-beta-cand-01 interfaces-V09-03-09 vis-V09-03-08

File:

• trunk/source/physics_lists/lists/src/HadronPhysicsFTFP_BERT.cc (modified) (5 diffs)

trunk/source/physics_lists/lists/src/HadronPhysicsFTFP_BERT.cc

@@ -24 +24 @@
 // ********************************************************************
 //
-// $Id: HadronPhysicsFTFP_BERT.cc,v 1.1 2007/10/19 15:35:08 gunter Exp $
-// GEANT4 tag $Name: geant4-09-03 $
+// $Id: HadronPhysicsFTFP_BERT.cc,v 1.3 2010/06/03 10:42:44 gunter Exp $
+// GEANT4 tag $Name: geant4-09-04-beta-cand-01 $
 //
 //---------------------------------------------------------------------------
@@ -49 +49 @@
 #include "G4BaryonConstructor.hh"
 #include "G4ShortLivedConstructor.hh"
+
+#include "G4QHadronInelasticDataSet.hh"
+
+HadronPhysicsFTFP_BERT::HadronPhysicsFTFP_BERT(G4int)
+                    :  G4VPhysicsConstructor("hInelastic FTFP_BERT")
+                     , QuasiElastic(false)
+{}
 
 HadronPhysicsFTFP_BERT::HadronPhysicsFTFP_BERT(const G4String& name, G4bool quasiElastic)
@@ -79 +86 @@
   theBertiniPiK->SetMaxEnergy(5*GeV);
   
-  theMiscLHEP=new G4MiscLHEPBuilder;
+  theMiscCHIPS=new G4MiscCHIPSBuilder;
 }
 
@@ -96 +103 @@
   delete theFTFPPro;    
     
-  delete theMiscLHEP;
+  delete theMiscCHIPS;
+  delete theCHIPSInelastic;
 }
 
@@ -118 +126 @@
   thePro->Build();
   thePiK->Build();
-  theMiscLHEP->Build();
+  // use CHIPS cross sections also for Kaons
+  theCHIPSInelastic = new G4QHadronInelasticDataSet();
+  
+  FindInelasticProcess(G4KaonMinus::KaonMinus())->AddDataSet(theCHIPSInelastic);
+  FindInelasticProcess(G4KaonPlus::KaonPlus())->AddDataSet(theCHIPSInelastic);
+  FindInelasticProcess(G4KaonZeroShort::KaonZeroShort())->AddDataSet(theCHIPSInelastic);
+  FindInelasticProcess(G4KaonZeroLong::KaonZeroLong())->AddDataSet(theCHIPSInelastic);
+
+  theMiscCHIPS->Build();
+}
+G4HadronicProcess* 
+HadronPhysicsFTFP_BERT::FindInelasticProcess(const G4ParticleDefinition* p)
+{
+  G4HadronicProcess* had = 0;
+  if(p) {
+     G4ProcessVector*  pvec = p->GetProcessManager()->GetProcessList();
+     size_t n = pvec->size();
+     if(0 < n) {
+       for(size_t i=0; i<n; ++i) {
+         if(fHadronInelastic == ((*pvec)[i])->GetProcessSubType()) {
+           had = static_cast<G4HadronicProcess*>((*pvec)[i]);
+           break;
+         }
+       }
+     }
+  }
+  return had;
 }