Adding Cross section dummy function

diff --git a/TGeant3/TGeant3.cxx b/TGeant3/TGeant3.cxx
index 141e0d869e4af73d5b921ac845df24a73f95039a..199ed01b0ebb5581f5c1084b02858f32e94bae15 100644 (file)
--- a/TGeant3/TGeant3.cxx
+++ b/TGeant3/TGeant3.cxx
@@ -976,6 +976,19 @@ const char* TGeant3::VolName(Int_t id) const
   return name;
 }
 
+//_____________________________________________________________________________
+Float_t TGeant3::Xsec(char* reac, Float_t energy, Int_t part, Int_t mate)
+{
+  Int_t gpart = IdFromPDG(part);
+  if(!strcmp(reac,"PHOT"))
+  {
+    if(part!=22) {
+      Error("Xsec","Can calculate photoelectric only for photons\n");
+    }
+  }
+  return 0;
+}
+
 //_____________________________________________________________________________
 void TGeant3::TrackPosition(TLorentzVector &xyz) const
 {

diff --git a/TGeant3/TGeant3.h b/TGeant3/TGeant3.h
index a2e26129c5f2c024dcc6dd30806122aa994a394c..0ea4517e2737ffa4688a45aba2e62384136c7db0 100644 (file)
--- a/TGeant3/TGeant3.h
+++ b/TGeant3/TGeant3.h
@@ -555,6 +555,7 @@ public:
   Int_t PDGFromId(Int_t pdg) const;
   void  DefineParticles();
   const char* VolName(Int_t id) const;
+  Float_t Xsec(char*, Float_t, Int_t, Int_t);
   void  TrackPosition(TLorentzVector &xyz) const;
   void  TrackMomentum(TLorentzVector &xyz) const;  
   Int_t NofVolumes() const;