27 #include "EventAction.hh"
33 #include "G4HCofThisEvent.hh"
34 #include "G4VHitsCollection.hh"
43 #include "SensitiveDetectorHit.hh"
45 #include "Analysis.hh"
89 for(
int i2=0;i2<3;i2++){
90 for(
int i1=0;i1<n_hit_sd;i1++)
103 G4double angXin = (ssd[1].
x() - ssd[0].
x()) / (ssd[1].
z() - ssd[0].
z());
104 G4double angYin = (ssd[1].
y() - ssd[0].
y()) / (ssd[1].
z() - ssd[0].
z());
109 double posXin = ssd[1].
x() - angXin * ssd[1].
z();
110 double posYin = ssd[1].
y() - angYin * ssd[1].
z();
115 G4double angXout = (ssd[2].
x() - posXin) / (ssd[2].
z());
116 G4double angYout = (ssd[2].
y() - posYin) / (ssd[2].
z());