Geant4  10.02.p03
ExExChEventAction Class Reference

#include <ExExChEventAction.hh>

Inheritance diagram for ExExChEventAction:
Collaboration diagram for ExExChEventAction:

Public Member Functions

 ExExChEventAction ()
 
virtual ~ExExChEventAction ()
 
virtual void BeginOfEventAction (const G4Event *)
 
virtual void EndOfEventAction (const G4Event *)
 
void SetVerbose (G4int val)
 
G4int GetVerbose () const
 
- Public Member Functions inherited from G4UserEventAction
 G4UserEventAction ()
 
virtual ~G4UserEventAction ()
 
virtual void SetEventManager (G4EventManager *value)
 

Private Attributes

G4int fSD_ID
 
G4int fVerboseLevel
 

Additional Inherited Members

- Protected Attributes inherited from G4UserEventAction
G4EventManagerfpEventManager
 

Detailed Description

Definition at line 38 of file ExExChEventAction.hh.

Constructor & Destructor Documentation

◆ ExExChEventAction()

ExExChEventAction::ExExChEventAction ( )

Definition at line 48 of file ExExChEventAction.cc.

49 {
50  fSD_ID = -1;
51  fVerboseLevel = 0;
52 }

◆ ~ExExChEventAction()

ExExChEventAction::~ExExChEventAction ( )
virtual

Definition at line 56 of file ExExChEventAction.cc.

56  {
57 }

Member Function Documentation

◆ BeginOfEventAction()

void ExExChEventAction::BeginOfEventAction ( const G4Event )
virtual

Reimplemented from G4UserEventAction.

Definition at line 61 of file ExExChEventAction.cc.

61  {
62 }

◆ EndOfEventAction()

void ExExChEventAction::EndOfEventAction ( const G4Event evt)
virtual

Reimplemented from G4UserEventAction.

Definition at line 66 of file ExExChEventAction.cc.

67 {
69 
70  if(fSD_ID == -1) {
71  G4String sdName;
72  if(SDman->FindSensitiveDetector(sdName="telescope",0)){
73  fSD_ID = SDman->GetCollectionID(sdName="telescope/collection");
74  }
75  }
76 
78 
79  G4HCofThisEvent * HCE = evt->GetHCofThisEvent();
80 
81  if(HCE)
82  {
83  if(fSD_ID != -1){
84  G4VHitsCollection* aHCSD = HCE->GetHC(fSD_ID);
86  }
87  }
88 
89  G4ThreeVector SSDposition[3];
90  SSDposition[0]= G4ThreeVector(0.,0.,0.);
91  SSDposition[1]= G4ThreeVector(0.,0.,0.);
92  SSDposition[2]= G4ThreeVector(0.,0.,0.);
93 
94  if(fSD)
95  {
96  int bTotalHits = 0;
97 
98  int n_hit_sd = fSD->entries();
99  for(int i2=0;i2<3;i2++){
100  for(int i1=0;i1<n_hit_sd;i1++)
101  {
102  ExExChSensitiveDetectorHit* aHit = (*fSD)[i1];
103  if(aHit->GetLayerID()==i2) {
104  SSDposition[i2] = aHit->GetWorldPos();
105  bTotalHits++;
106  }
107  }
108  }
109 
110  if(bTotalHits == 3){
111  double fAngXin = (SSDposition[1].x() - SSDposition[0].x());
112  fAngXin /= (SSDposition[1].z() - SSDposition[0].z());
113  double fAngYin = (SSDposition[1].y() - SSDposition[0].y());
114  fAngYin /= (SSDposition[1].z() - SSDposition[0].z());
115  double fPosXin = SSDposition[1].x();
116  double fPosYin = SSDposition[1].y();
117  double fAngXout = (SSDposition[2].x() - SSDposition[1].x());
118  fAngXout /= (SSDposition[2].z() - SSDposition[1].z());
119  double fAngYout = (SSDposition[2].y() - SSDposition[1].y());
120  fAngYout /= (SSDposition[2].z() - SSDposition[1].z());
121 
122  G4AnalysisManager* analysisManager = G4AnalysisManager::Instance();
123 
124  analysisManager->FillNtupleDColumn(0, fAngXin * 1.E6 * CLHEP::rad);
125  analysisManager->FillNtupleDColumn(1, fAngYin * 1.E6 * CLHEP::rad);
126  analysisManager->FillNtupleDColumn(2, fPosXin / CLHEP::mm);
127  analysisManager->FillNtupleDColumn(3, fPosYin / CLHEP::mm);
128  analysisManager->FillNtupleDColumn(4, fAngXout * 1.E6 * CLHEP::rad);
129  analysisManager->FillNtupleDColumn(5, fAngYout * 1.E6 * CLHEP::rad);
130  analysisManager->AddNtupleRow();
131  }
132  }
133 }
G4int GetCollectionID(G4String colName)
Definition: G4SDManager.cc:135
CLHEP::Hep3Vector G4ThreeVector
static const double rad
G4bool FillNtupleDColumn(G4int id, G4double value)
G4VSensitiveDetector * FindSensitiveDetector(G4String dName, G4bool warning=true)
Definition: G4SDManager.cc:128
double x() const
static const double mm
Definition: SystemOfUnits.h:94
double y() const
static G4SDManager * GetSDMpointer()
Definition: G4SDManager.cc:40
G4HCofThisEvent * GetHCofThisEvent() const
Definition: G4Event.hh:185
double z() const
Here is the call graph for this function:

◆ GetVerbose()

G4int ExExChEventAction::GetVerbose ( ) const
inline

Definition at line 53 of file ExExChEventAction.hh.

53 { return fVerboseLevel; }

◆ SetVerbose()

void ExExChEventAction::SetVerbose ( G4int  val)
inline

Definition at line 52 of file ExExChEventAction.hh.

52 { fVerboseLevel = val; }

Member Data Documentation

◆ fSD_ID

G4int ExExChEventAction::fSD_ID
private

Definition at line 48 of file ExExChEventAction.hh.

◆ fVerboseLevel

G4int ExExChEventAction::fVerboseLevel
private

Definition at line 49 of file ExExChEventAction.hh.


The documentation for this class was generated from the following files: