--- PamVMC/trk/include/PamVMCTrkSD.h	2009/02/18 17:41:08	1.1
+++ PamVMC/trk/include/PamVMCTrkSD.h	2009/06/12 18:39:49	1.5
@@ -66,10 +66,10 @@
  //which should present in the same directory then application library
  
   void Init(){ 
-    TString G4WORKDIR=gSystem->Getenv("G4WORKDIR");
+    TString PAM_VMC=gSystem->Getenv("PAM_VMC");
     TString PLATFORM=gSystem->Getenv("PLATFORM");
     TFile * resspe = 
-      new TFile(G4WORKDIR+"/lib/tgt_"+PLATFORM+"/resxy_new.root");
+      new TFile(PAM_VMC+"/lib/tgt_"+PLATFORM+"/resxy_new.root");
     if(resspe->IsOpen()){
       TH1F * hsx;
       TH1F * hsy;
@@ -98,7 +98,7 @@
        delete hsx;
        delete hsy;
     } else cout<<"!!!WARNING: SPE hit resolution data file resxy_new.root"
-	     <<" is not present in: " <<G4WORKDIR+"/lib/tgt_"+PLATFORM
+	     <<" is not present in: " <<PAM_VMC+"/lib/tgt_"+PLATFORM
 	     <<" or corrupted!!!"<<endl;
 
     delete resspe;
@@ -148,7 +148,7 @@
     return hit;
   }
 
-  void FillHit(fin f,TVirtualMC *g){
+ void FillHit(fin f,TVirtualMC *g, Bool_t is_prim){
     ffi=f;
     //g->SetMaxStep(0.001); // 10 um
     switch(f) {
@@ -158,15 +158,22 @@
       InitHit();
       //cout<<">>>ENTERING in PLANE: "<<GetPlaneID()<<endl;
     default: 
-      UpdateHit(g);
+      UpdateHit(g, is_prim);
       //TParticlePDG* particlePDG = TDatabasePDG::Instance()->GetParticle(gMC->TrackPid());
       //Double_t mass = particlePDG->Mass();
-      //Double_t x,y,z,e;
+      //Double_t x[3];
+      //Double_t B[3];
+      //Double_t * xp = &x[0];
+      //Double_t *Bp = &B[0];
+      //g->TrackPosition(x[0],x[1],x[2]);
       //g->TrackMomentum(x,y,z,e);
       //Double_t P0 = TMath::Sqrt(x*x+y*y+z*z);
       //Double_t energy =  TMath::Sqrt(P0*P0+mass*mass)-mass;
       //cout<<"stepping... Stepsize"<<gMC->TrackStep()*10000<<"um "<<" particle PDG is:"<<gMC->TrackPid()<<" Kin.energy:"<<energy*1000000<<" keV"<<endl;
-     
+      //PamVMCFieldMgr::Instance()->Field(xp,Bp);
+      //cout<<"X,Y,Z: "<<x[0]<<" "<<x[1]<<" "<<x[2]<<" "<<"Bx, By, Bz: "<<B[0]<<" "<<B[1]<<" "<<B[2]<<endl;
+      //fhit->Print();
+      //cout<<"Edep"<<g->Edep()*1e6<<" kev"<<endl;
       if (g->Edep()) PreDigit(); 
       
       break;      
@@ -176,7 +183,10 @@
     case EXITING:
       //cout<<"<<<EXITING out"<<gMC->TrackPid()<<endl;
     // Save hit if energy release is greater than zero
-      if(fhit->GetEREL()){
+      //cout<<"PDG:"<<fhit->GetPDG()<<"TOF="<<fhit->GetTOF()<<endl;
+      
+      if(fhit->GetEREL() ){
+	//cout<<"Saving.."<<endl;
 	SaveTrkHit(fdname.Data());
       }
       break;
@@ -286,6 +296,78 @@
 
 };
 
+#endif // PAMVMCTRKSD_H
 
 
-#endif // PAMVMCTRKSD_H
+
+#ifndef PAMVMCTPANSD_H
+#define PAMVMCTPANSD_H
+
+class PamVMCTPANSD: public PamVMCDetectorSD{
+
+public:
+ PamVMCTPANSD():PamVMCDetectorSD("PamVMCDetectorHit","TPAN",1000)
+    {
+    }; 
+  virtual void FillHit(fin f,TVirtualMC *g, Bool_t is_prim){
+    ffi=f;
+    switch(f) {
+    case ENTERING:
+      CleanHit();
+      FillVolID();
+      InitHit();         
+    default: 
+      UpdateHit(g, is_prim);
+      break;      
+    }
+    
+    switch(f){
+    case EXITING:
+      // Save hit if PATH, not energy release is greater than zero
+      if(fhit->GetPATH()){
+	SaveHit(fdname.Data());
+      }
+      break;
+    default:
+      break;
+    }
+  }      
+  ClassDef(PamVMCTPANSD,1)
+    
+    };
+
+#endif // PAMVMCTPAN_H
+
+
+#ifndef PAMVMCTRCNSD_H
+#define PAMVMCTRCNSD_H
+
+class PamVMCTRCNSD: public PamVMCDetectorSD{
+
+public:
+ PamVMCTRCNSD():PamVMCDetectorSD("PamVMCDetectorHit","TRCN",1000)
+    {
+      
+    };
+  ClassDef(PamVMCTRCNSD,1)
+    
+    };
+
+#endif // PAMVMCTRCN_H
+
+
+#ifndef PAMVMCTRSLSD_H
+#define PAMVMCTRSLSD_H
+
+class PamVMCTRSLSD: public PamVMCDetectorSD{
+
+public:
+ PamVMCTRSLSD():PamVMCDetectorSD("PamVMCDetectorHit","TRSL",1000)
+    {
+      
+    };
+  ClassDef(PamVMCTRSLSD,1)
+    
+    };
+
+#endif // PAMVMCTRSL_H