70 for (j =0 ; j<8 ; j++) {
ETailVis[j] = 0.;};
71 for (j =0 ; j<7 ; j++) {
ETailDep[j] = 0.;};
81 if (edep==0.)
return false;
88 if(strcmp(physVol->
GetName(),
"ScintS1Physical")==0) {
91 if(strcmp(physVol->
GetName(),
"ScintS2Physical")==0) {
94 if(strcmp(physVol->
GetName(),
"ScintS3Physical")==0) {
98 if(strcmp(physVol->
GetName(),
"HoleCntrScintPhysical")==0){
107 if(strcmp(physVol->
GetName(),
"BigScintPhysical")==0) {
111 if(strcmp(physVol->
GetName(),
"SmallScintPhysical")==0) {
116 if(strcmp(physVol->
GetName(),
"BigIronPhysical")==0) {
120 if(strcmp(physVol->
GetName(),
"SmallIronPhysical")==0) {
132 G4cout <<
" Visisble Energy in S1 , S2 , S3 in (MeV)" <<
G4endl;
135 G4cout <<
" Visible Energy in Hole Counter (MeV) " <<
G4endl;
138 G4cout <<
" Visible Energy in Upstream Dead Materials " <<
G4endl;
141 G4cout <<
" Visible Energy in Tail Catcher Scintillator" <<
G4endl;
145 G4cout <<
" Visible Energy in Tail Catcher Absorber" <<
G4endl;
150 G4String FileName =
"Beam_802_1mm.dat";
151 std::ios::openmode iostemp;
153 iostemp = std::ios::out;
159 std::ofstream BeamDatafile(FileName, iostemp);
165 for (j=1; j<8 ; j++) { BeamDatafile <<
ETailVis[j]/
MeV <<
" " ;} ; BeamDatafile << std::endl;
166 for (j=1; j<7 ; j++) { BeamDatafile <<
ETailDep[j]/
MeV <<
" " ;} ; BeamDatafile << std::endl;
168 BeamDatafile.close();