Day5, 10/04/2011

From New IAC Wiki
Revision as of 22:06, 4 October 2011 by Shaproma (talk | contribs)
Jump to navigation Jump to search

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