6 #include "G4VPhysicalVolume.hh"
7 #include "G4VTouchable.hh"
13 G4VPhysicalVolume* physVol;
16 unsigned int roman_pot = 0;
19 const G4VTouchable* touch = aStep->GetPreStepPoint()->GetTouchable();
21 for (
int ii = 0;
ii < touch->GetHistoryDepth();
ii++) {
22 physVol = touch->GetVolume(
ii);
23 if (physVol->GetName() ==
"RP_Silicon_Detector") {
25 }
else if (physVol->GetName() ==
"RP_box_primary_vacuum") {
26 int cpy_no = physVol->GetCopyNo();
27 arm = (cpy_no / 100) % 10;
29 roman_pot = cpy_no % 10;