PamCAL/include/PamVMCTofSD.h

#ifndef PAMVMCTOFSD_H
#define PAMVMCTOFSD_H
#include <iostream>

#include "PamVMCDetectorSD.h"

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

class PamVMCTofSD: public PamVMCDetectorSD{

    const char* cname;

public:
  PamVMCTofSD(char *dname):PamVMCDetectorSD("pHitData",dname,1000)
{
     cname=dname;
};


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

  
     switch(f) {
    case ENTERING:
     
       CleanHit();
      FillVolID();
      InitHit();
    
    default:
     UpdateHit(g);
    
   cout<<"NEW HIT ENTERING OR  INSIDE "<<fnohit<<" in "<<cname<<" PDG="<<fhit.GetPDG()<<" GT "<<g->GetStack()->GetCurrentTrack()->T()<<" 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()){
    // AddGlobalTime(g);
    //cout<<"saving"<<fdetID->GetID()<<endl;
    cout<<"NEW HIT EXITING "<<fnohit<<" in "<<cname<<" PDG="<<fhit.GetPDG()<<" at"<<fhit.GetTOF()<<" GT  "<<g->GetStack()->GetCurrentTrack()->T() <<endl;
          SaveHit(cname);
  } else {
    cout << "TOF FillVolID false" <<endl; 
  }
    }
    break;
  default:
    break;
  }
  }
  
  ClassDef(PamVMCTofSD,1)

    };

#endif // PAMVMCTOFSD_H
1	#ifndef PAMVMCTOFSD_H
2	#define PAMVMCTOFSD_H
3	#include <iostream>
4
5	#include "PamVMCDetectorSD.h"
6
7	using std::cout;
8	using std::endl;
9
10	class PamVMCTofSD: public PamVMCDetectorSD{
11
12	const char* cname;
13
14	public:
15	PamVMCTofSD(char *dname):PamVMCDetectorSD("pHitData",dname,1000)
16	{
17	cname=dname;
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
34	cout<<"NEW HIT ENTERING OR INSIDE "<<fnohit<<" in "<<cname<<" PDG="<<fhit.GetPDG()<<" GT "<<g->GetStack()->GetCurrentTrack()->T()<<" 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	// AddGlobalTime(g);
45	//cout<<"saving"<<fdetID->GetID()<<endl;
46	cout<<"NEW HIT EXITING "<<fnohit<<" in "<<cname<<" PDG="<<fhit.GetPDG()<<" at"<<fhit.GetTOF()<<" GT "<<g->GetStack()->GetCurrentTrack()->T() <<endl;
47	SaveHit(cname);
48	} else {
49	cout << "TOF FillVolID false" <<endl;
50	}
51	}
52	break;
53	default:
54	break;
55	}
56	}
57
58	ClassDef(PamVMCTofSD,1)
59
60	};
61
62	#endif // PAMVMCTOFSD_H