45 #include "G4HCofThisEvent.hh"
46 #include "G4VHitsCollection.hh"
61 #include <AIDA/AIDA.h>
67 :calorimeterCollID(-1),drawFlag(
"all"),printModulo(10), StepAction(SA), eventMessenger(0)
77 delete eventMessenger;
88 if (evtNb%printModulo == 0)
90 G4cout <<
"\n---> Begin of event: " << evtNb+1 <<
G4endl;
94 NTracksOutOfWorld = 0;
104 #ifdef G4ANALYSIS_USE
110 G4cout <<
"N Tracks out of world " << NTracksOutOfWorld <<
G4endl;
114 G4String FileName1 =
"OutTracks_802_1mm.dat";
115 std::ios::openmode iostemp1;
117 iostemp1 = std::ios::out;
122 std::ofstream OutTracks(FileName1, iostemp1);
124 OutTracks << NTracksOutOfWorld <<
G4endl;
127 for(i=1; i<= NTracksOutOfWorld ; i++){
128 for(j=1; j<11 ; j++) {
131 OutTracks << OutOfWorld <<
" " ;
133 OutTracks << std::endl;
138 #ifdef G4ANALYSIS_USE
139 analysis->getfhisto_1()->fill(NTracksOutOfWorld);
147 G4String FileName2 =
"SecndTracks_802_1mm.dat";
148 std::ios::openmode iostemp2;
150 iostemp2 = std::ios::out;
156 std::ofstream SecndTracks(FileName2, iostemp2);
158 SecndTracks << NSecondaries << std::endl;
161 for(i=1; i<= NSecondaries ; i++){
162 for(j=1; j<11 ; j++) {
164 SecndTracks << Secondary <<
" " ;
166 SecndTracks << std::endl;
171 #ifdef G4ANALYSIS_USE
172 analysis->getfhisto_2()->fill(NSecondaries);
176 G4String FileName3 =
"EdepFCAL_802_1mm.dat";
177 std::ios::openmode iostemp3;
179 iostemp3 = std::ios::out;
185 std::ofstream EdepFCAL(FileName3, iostemp3);
190 EdepFCAL << EmEdep <<
" ";
192 EdepFCAL << std::endl;
195 #ifdef G4ANALYSIS_USE
196 analysis->getfhisto_3()->fill(EmEdep);
197 analysis->getfhisto_4()->fill(HadEdep);