--- PamVMC/include/PamVMCDetectorHit.h	2009/02/19 16:50:36	1.1
+++ PamVMC/include/PamVMCDetectorHit.h	2009/06/12 18:39:17	1.5
@@ -28,10 +28,11 @@
   Double_t fXIN;    
   Double_t fYIN;
   Double_t fZIN;
+  Bool_t   fPRIM; //1 if it was a hit from primary particle
 
   pHitData():  fPOS(0), fPDG(0), fPATH(0.), fTOF(0.), fPIN(0.), fXOUT (0.), 
-    fYOUT (0.),  fZOUT (0.), fEREL (0.), 
-      fXIN(0.),  fYIN (0.), fZIN (0.) { };
+    fYOUT (0.),  fZOUT (0.), fEREL (0.),
+    fXIN(0.),  fYIN (0.), fZIN (0.), fPRIM(kFALSE) { };
 
 
   // Kept for back-compatibility
@@ -40,6 +41,7 @@
     fPATH= fTOF= fPIN= fXOUT = 
       fYOUT =  fZOUT = fEREL = 
       fXIN=  fYIN = fZIN = 0.;
+    fPRIM = kFALSE;
   }
 
   // Inherithed from TObject
@@ -60,7 +62,8 @@
          <<"ZOUT: " << fZOUT << " (cm)"<<endl
          <<"PATH: " << fPATH << " (cm)"<<endl
          <<"EREL: " << fEREL*1E6 << " (keV)"<<endl
-         <<"TOF: " << fTOF << " (s)"<<endl; } 
+         <<"TOF: " << fTOF << " (s)"<<endl
+	 <<"PRIMARY: "<<fPRIM<<endl; } 
   
   ClassDef(pHitData,1)
 
@@ -96,7 +99,7 @@
 
     virtual void Print(Option_t* option = "") const;
 
-    void FillHit(TVirtualMC *g) {
+    void FillHit(TVirtualMC *g, Bool_t is_prim) {
       fnocall++;
       double x,y,z,e;
       g->TrackPosition(x,y,z);
@@ -112,9 +115,10 @@
       AddTOF((g->GetStack()->GetCurrentTrack()->T())/1.0E9 /*TG4G3Units::Time()*/+(g->TrackTime()));
       AddPATH(g->TrackStep());
       AddEREL(g->Edep());
-
+      SetPRIM(is_prim);
     }    
   
+
   void SetHit( const pHitData &phd){ 
       fnocall++; 
       SetPOS(phd.fPOS);
@@ -129,6 +133,7 @@
       SetXOUT( phd.fXOUT );
       SetYOUT( phd.fYOUT );
       SetZOUT( phd.fZOUT );
+      SetPRIM( phd.fPRIM );
     }
 
     void UpdXIN(Double_t x){ 
@@ -153,6 +158,7 @@
     void SetXOUT   (Double_t xout)  { fhitdata.fXOUT = xout; }
     void SetYOUT   (Double_t yout)  { fhitdata.fYOUT = yout; }
     void SetZOUT   (Double_t zout)  { fhitdata.fZOUT = zout; }
+    void SetPRIM   (Bool_t prim)    { fhitdata.fPRIM = prim; }
     void Reset(); 
                  
     // Get methods
@@ -170,6 +176,7 @@
     Double_t GetXOUT()     { return fhitdata.fXOUT; };
     Double_t GetYOUT()     { return fhitdata.fYOUT; };
     Double_t GetZOUT()     { return fhitdata.fZOUT; };
+    Bool_t GetPRIM()       { return fhitdata.fPRIM; };
     
  protected: