Special particles added in AddParticlesToPdgDataBase()
authormorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 30 Oct 2007 15:58:48 +0000 (15:58 +0000)
committermorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 30 Oct 2007 15:58:48 +0000 (15:58 +0000)
TFluka/TFluka.cxx
TFluka/TFluka.h

index a9e6d12..d9492fb 100644 (file)
@@ -256,7 +256,7 @@ void TFluka::Init() {
     // Add ions to PDG Data base
     //
      AddParticlesToPdgDataBase();
-     //
+    //
 }
 
 
@@ -2542,6 +2542,15 @@ void TFluka::AddParticlesToPdgDataBase() const
                      khShGev/(12.33*kYear2Sec),6,"Ion",GetIonPdg(2,4));
   pdgDB->AddParticle("HE3","HE3",3*kAu2Gev+14.931e-3,kFALSE,
                      0,6,"Ion",GetIonPdg(2,3));
+//
+//
+//
+// Special particles
+//
+  pdgDB->AddParticle("Cherenkov","Cherenkov",0,kFALSE,
+                     0,0,"Special",GetSpecialPdg(50));
+  pdgDB->AddParticle("FeedbackPhoton","FeedbackPhoton",0,kFALSE,
+                     0,0,"Special",GetSpecialPdg(51));
 }
 
 //
@@ -2597,6 +2606,15 @@ Int_t TFluka::GetIonPdg(Int_t z, Int_t a, Int_t i) const
 
   return 1000000000 + 10*1000*z + 10*a + i;
 }  
+
+//__________________________________________________________________
+Int_t TFluka::GetSpecialPdg(Int_t number) const
+{
+// Numbering for special particles
+
+  return 50000000 + number;
+}                
+
      
 void  TFluka::PrimaryIonisationStepping(Int_t nprim)
 {
index b09857d..5e4e545 100644 (file)
@@ -409,6 +409,8 @@ class TFluka : public TVirtualMC {
   void  PrintHeader();
   void  AddParticlesToPdgDataBase() const;
   Int_t GetIonPdg(Int_t z, Int_t a, Int_t i = 0) const;
+  Int_t GetSpecialPdg(Int_t number) const;
+  
   Float_t* CreateFloatArray(Double_t* array, Int_t size) const;
   
   //