PamCAL/include/PamVMCCaloSD.h

#ifndef PAMVMCCALOSD_H
#define PAMVMCCALOSD_H
#include <iostream>

#include "PamVMCDetectorSD.h"

#define NOCALHIT 4224

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

class PamVMCCaloSD: public PamVMCDetectorSD{
public:
  PamVMCCaloSD(): PamVMCDetectorSD("pHitData","CAST",NOCALHIT){
    //  for(Int_t i=0; i< NOCALHIT; ++i) CreateHit(i);
  }


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

    switch(f) {
    case ENTERING:
      CleanHit();
      FillVolID();
      InitHit();
      
     
      //FillVolID();
      // cout<<"["<<GetViewID()<<"] ["<<GetPlaneID()<<"] ["<<GetStripID()<<"] "<<fhit.GetPOS()<<"\n";
       UpdateHit(g); 
    default:
       UpdateHit(g);
      break;      
    }

    switch(f){
    case EXITING:
       UpdateHit(g);
       // AddGlobalTime(g);
      // Save hit if energy release is greater than zero 
      if(fhit.GetEREL()){
        if(fdetID->FillVolID()){
           SaveHit(fdetID->GetID( ));
           //           cout<<fhit.GetPOS()<<"\n";
        } else {
          cout << "Calorimeter FillVolID false" <<endl; 
        }
      }
      break;
    default:
      break;
    }
  }

  //virtual void UpdateHit(TVirtualMC *g){


  //}
  
  //   virtual void ClearHitColl(){
     //do nothing... for calorimerer only
  //    } 

  Int_t GetPlaneID(){
    return ((int)fhit.GetPOS()/192);
  }

  Int_t GetViewID(){
    return ((int)fhit.GetPOS()/96%2==0?0:1);
  }
  Int_t GetStripID(){
    return (fhit.GetPOS()%96);
  }

 
  ClassDef(PamVMCCaloSD,1)

};

#endif // PAMVMCCALOSD_H
1	nikolas	1.1	#ifndef PAMVMCCALOSD_H
2			#define PAMVMCCALOSD_H
3			#include <iostream>
4
5			#include "PamVMCDetectorSD.h"
6
7			#define NOCALHIT 4224
8
9			using std::cout;
10			using std::endl;
11
12			class PamVMCCaloSD: public PamVMCDetectorSD{
13			public:
14			PamVMCCaloSD(): PamVMCDetectorSD("pHitData","CAST",NOCALHIT){
15			// for(Int_t i=0; i< NOCALHIT; ++i) CreateHit(i);
16			}
17
18
19			virtual void FillHit(fin f,TVirtualMC *g){
20
21			switch(f) {
22			case ENTERING:
23			CleanHit();
24			FillVolID();
25			InitHit();
26
27
28			//FillVolID();
29			// cout<<"["<<GetViewID()<<"] ["<<GetPlaneID()<<"] ["<<GetStripID()<<"] "<<fhit.GetPOS()<<"\n";
30			UpdateHit(g);
31			default:
32			UpdateHit(g);
33			break;
34			}
35
36			switch(f){
37			case EXITING:
38			UpdateHit(g);
39			// AddGlobalTime(g);
40			// Save hit if energy release is greater than zero
41			if(fhit.GetEREL()){
42			if(fdetID->FillVolID()){
43			SaveHit(fdetID->GetID( ));
44			// cout<<fhit.GetPOS()<<"\n";
45			} else {
46			cout << "Calorimeter FillVolID false" <<endl;
47			}
48			}
49			break;
50			default:
51			break;
52			}
53			}
54
55			//virtual void UpdateHit(TVirtualMC *g){
56
57
58			//}
59
60			// virtual void ClearHitColl(){
61			//do nothing... for calorimerer only
62			// }
63
64			Int_t GetPlaneID(){
65			return ((int)fhit.GetPOS()/192);
66			}
67
68			Int_t GetViewID(){
69			return ((int)fhit.GetPOS()/96%2==0?0:1);
70			}
71			Int_t GetStripID(){
72			return (fhit.GetPOS()%96);
73			}
74
75
76
77			ClassDef(PamVMCCaloSD,1)
78
79			};
80
81			#endif // PAMVMCCALOSD_H