Day5, 10/04/2011

go back

/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();

 }