Day5, 10/04/2011
Jump to navigation
Jump to search
/home/foretony/geant4/NucSim/Homewrk5
- gedit src/ExN02SteppingVerbose.cc
ExN02SteppingVerbose::ExN02SteppingVerbose() {
outfile.open("sigmaR.txt"); //added 10/4/11
}
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
ExN02SteppingVerbose::~ExN02SteppingVerbose() {
outfile.close(); //added 10/4/11
}
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
- void ExN02SteppingVerbose::TrackingStarted() print first step
- gedit include/ExN02SteppingVerbose.hh
class ExN02SteppingVerbose : public G4SteppingVerbose {
public: ExN02SteppingVerbose(); ~ExN02SteppingVerbose();
void StepInfo(); void TrackingStarted();
Float SigmaRoverR; //added 10/4/11 Float FinalEnergy; //added 10/4/11 Float PreviousEnergy; //added 10/4/11
private: //added 10/4/11 std::ofstream outfile;
};
- gedit src/ExN02SteppingVerbose.cc
void ExN02SteppingVerbose::StepInfo() {
CopyState(); G4int prec = G4cout.precision(3);
//outfile << "Roman Olga: " << fTrack->GetVolume()->GetName() << G4endl; //added 10/4/11 try for fun //G4cout << "Roman Olga: " << fTrack->GetVolume()->GetName() << G4endl; //added 10/4/11 try for fun
// if( verboseLevel >= 1 ){ if( fTrack->GetVolume()->GetName() =="World") { if(SigmaRoverR) { outfile << SigmaRoverR << G4endl; G4cout << " R= " << SigmaRoverR << G4endl; } SigmaRoverR=0; } // if( fTrack->GetVolume()->GetName() =="Target" && fStep->GetPostStepPoint()->GetProcessDefinedStep()->GetProcessName() == "hIoni"){ if( fTrack->GetVolume()->GetName() =="Target" ){ SigmaRoverR+=fStep->GetStepLength();
}