Difference between revisions of "Day5, 10/04/2011"

From New IAC Wiki
Jump to navigation Jump to search
(Created page with '[https://wiki.iac.isu.edu/index.php/Simulations_of_Particle_Interactions_with_Matter._Dr._Forest._Fall_2011 go back]')
 
 
(13 intermediate revisions by the same user not shown)
Line 1: Line 1:
 
[https://wiki.iac.isu.edu/index.php/Simulations_of_Particle_Interactions_with_Matter._Dr._Forest._Fall_2011 go back]
 
[https://wiki.iac.isu.edu/index.php/Simulations_of_Particle_Interactions_with_Matter._Dr._Forest._Fall_2011 go back]
 +
 +
 +
/home/foretony/geant4/NucSim/Homewrk5
 +
 +
 +
* gedit src/ExN02SteppingVerbose.cc
 +
//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 +
 +
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::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();
 +
 +
  }
 +
 +
 +
** 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;
 +
 +
};

Latest revision as of 22:07, 4 October 2011

go back


/home/foretony/geant4/NucSim/Homewrk5


  • gedit src/ExN02SteppingVerbose.cc

//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

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


    • 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;

};