PamCAL/include/PamVMCTrkSD.h

#ifndef PAMVMCTRKSD_H
#define PAMVMCTRKSD_H
#include <iostream>

#include "PamVMCDetectorSD.h"

using std::cout;
using std::endl;

class PamVMCTrkSD: public PamVMCDetectorSD{

 
public:
  PamVMCTrkSD():PamVMCDetectorSD("pHitData","TSPA",1000)
{
  
};


  /* virtual void FillHit(fin f,TVirtualMC *g){ 

  
     switch(f) {
    case ENTERING:
     
       CleanHit();
      FillVolID();
      InitHit(); 
    
    default:
     UpdateHit(g);
     cout<<"plane "<<GetPlaneID()<<" pad "<<GetPadID();
     cout<<" NEW HIT ENTERING OR  INSIDE "<<fnohit<<" PDG="<<fhit.GetPDG()<<" at "<<g->TrackTime()<<" fTOF is: "<<fhit.GetTOF()<<endl;
     break;      
   }

  switch(f){
  case EXITING:
  // Save hit if energy release is greater than zero
    UpdateHit(g); 
   if(fhit.GetEREL()){
  if(fdetID->FillVolID()){
      cout<<"plane "<<GetPlaneID()<<" pad "<<GetPadID();
      cout<<" NEW HIT EXITING "<<fnohit<<" PDG="<<fhit.GetPDG()<<" at "<<g->TrackTime()<<" fTOF is: "<<fhit.GetTOF()<<endl;
          SaveHit("TSPA");
  } else {
    cout << "TOF FillVolID false" <<endl; 
  }
    }
    break;
  default:
    break;
  }
  }

  */
Int_t GetPlaneID(){
    return ((int)fhit.GetPOS()/6+1);
}

Int_t GetPadID(){
    return (fhit.GetPOS()-(GetPlaneID()-1)*6);
}

  ClassDef(PamVMCTrkSD,1)

    };

#endif // PAMVMCTRKSD_H
1	#ifndef PAMVMCTRKSD_H
2	#define PAMVMCTRKSD_H
3	#include <iostream>
4
5	#include "PamVMCDetectorSD.h"
6
7	using std::cout;
8	using std::endl;
9
10	class PamVMCTrkSD: public PamVMCDetectorSD{
11
12
13
14	public:
15	PamVMCTrkSD():PamVMCDetectorSD("pHitData","TSPA",1000)
16	{
17
18	};
19
20
21	/* virtual void FillHit(fin f,TVirtualMC *g){
22
23
24	switch(f) {
25	case ENTERING:
26
27	CleanHit();
28	FillVolID();
29	InitHit();
30
31	default:
32	UpdateHit(g);
33	cout<<"plane "<<GetPlaneID()<<" pad "<<GetPadID();
34	cout<<" NEW HIT ENTERING OR INSIDE "<<fnohit<<" PDG="<<fhit.GetPDG()<<" at "<<g->TrackTime()<<" fTOF is: "<<fhit.GetTOF()<<endl;
35	break;
36	}
37
38	switch(f){
39	case EXITING:
40	// Save hit if energy release is greater than zero
41	UpdateHit(g);
42	if(fhit.GetEREL()){
43	if(fdetID->FillVolID()){
44	cout<<"plane "<<GetPlaneID()<<" pad "<<GetPadID();
45	cout<<" NEW HIT EXITING "<<fnohit<<" PDG="<<fhit.GetPDG()<<" at "<<g->TrackTime()<<" fTOF is: "<<fhit.GetTOF()<<endl;
46	SaveHit("TSPA");
47	} else {
48	cout << "TOF FillVolID false" <<endl;
49	}
50	}
51	break;
52	default:
53	break;
54	}
55	}
56
57	*/
58	Int_t GetPlaneID(){
59	return ((int)fhit.GetPOS()/6+1);
60	}
61
62	Int_t GetPadID(){
63	return (fhit.GetPOS()-(GetPlaneID()-1)*6);
64	}
65
66	ClassDef(PamVMCTrkSD,1)
67
68	};
69
70	#endif // PAMVMCTRKSD_H