BesSteppingAction Class Reference

#include <BesSteppingAction.hh>

List of all members.

Public Member Functions
	BesSteppingAction ()
	~BesSteppingAction ()
void	UserSteppingAction (const G4Step *)

---

Detailed Description

Definition at line 7 of file BesSteppingAction.hh.

---

Constructor & Destructor Documentation

BesSteppingAction::BesSteppingAction

(

)

Definition at line 10 of file BesSteppingAction.cc.

{ }

BesSteppingAction::~BesSteppingAction

(

)

Definition at line 13 of file BesSteppingAction.cc.

{ }

---

Member Function Documentation

void BesSteppingAction::UserSteppingAction

(

const G4Step *

)

Definition at line 17 of file BesSteppingAction.cc.

References abs.

{

  G4StepPoint* postStep = currentStep->GetPostStepPoint();
  G4ThreeVector currentPosition = postStep->GetPosition();
  G4Track* currentTrack = currentStep->GetTrack();
  if(std::abs(currentPosition.x()) > 263.5*cm||
     std::abs(currentPosition.y()) > 263.5*cm||
     std::abs(currentPosition.z()) > 287.5*cm){
    G4cout<<"Out of World!!!"<<G4endl;
    currentTrack->SetTrackStatus(fKillTrackAndSecondaries);
  }else if(currentTrack->GetCurrentStepNumber()>=20000){ 
    G4cout<<"StepNumber>=20000 !!!"<<G4endl;
    currentTrack->SetTrackStatus(fKillTrackAndSecondaries);
  }
/*     G4cout.precision(15);
     G4cout<<"#Step# "<<currentTrack->GetCurrentStepNumber()<<" pName "<<currentTrack->GetDefinition()->GetParticleName()<<" prex prey prez "<<currentStep->GetPreStepPoint()->GetPosition().x()/mm<<" "<<currentStep->GetPreStepPoint()->GetPosition().y()/mm<<" "<<currentStep->GetPreStepPoint()->GetPosition().z()/mm<<G4endl;
     G4cout<<"prepx prepy prepz "<<currentStep->GetPreStepPoint()->GetMomentum().x()/GeV<<" "<<currentStep->GetPreStepPoint()->GetMomentum().y()/GeV<<" "<<currentStep->GetPreStepPoint()->GetMomentum().z()/GeV<<G4endl;
     G4cout<<"post step postx posty postz "<<postStep->GetPosition().x()/mm<<" "<<postStep->GetPosition().y()/mm<<" "<<postStep->GetPosition().z()/mm<<G4endl;
     G4cout<<"postpx postpy postpz "<<postStep->GetMomentum().x()/GeV<<" "<<postStep->GetMomentum().y()/GeV<<" "<<postStep->GetMomentum().z()/GeV<<G4endl;         
      G4cout << G4endl;
      G4cout << std::setw( 5) << "#Step#"     << " "
             << std::setw( 9) << "pName"      << " "
             << std::setw( 6) << "X"          << "    "
             << std::setw( 6) << "Y"          << "    "
             << std::setw( 6) << "Z"          << "    "
             << std::setw( 9) << "KineE"      << " "
             << std::setw( 9) << "dEStep"     << " "
             << std::setw(10) << "StepLeng"
             << std::setw(10) << "Volume"    << "  "
             << std::setw(10) << "Process"  
             << G4endl;
    

    G4cout << std::setw(5) << currentTrack->GetCurrentStepNumber() << " "
        << std::setw(9) << currentTrack->GetDefinition()->GetParticleName()
        << std::setw(6) << G4BestUnit(currentTrack->GetPosition().x(),"Length")
        << std::setw(6) << G4BestUnit(currentTrack->GetPosition().y(),"Length")
        << std::setw(6) << G4BestUnit(currentTrack->GetPosition().z(),"Length")
        << std::setw(6) << G4BestUnit(currentTrack->GetKineticEnergy(),"Energy")
        << std::setw(6) << G4BestUnit(currentStep->GetTotalEnergyDeposit(),"Energy")
        << std::setw(6) << G4BestUnit(currentStep->GetStepLength(),"Length");
    if( currentTrack->GetNextVolume() != 0 ) { 
        G4cout << std::setw(10) << currentTrack->GetVolume()->GetName();
     } else {
       G4cout << std::setw(10) << "OutOfWorld";
     }

    if(currentStep->GetPostStepPoint()->GetProcessDefinedStep() != NULL){
      G4cout << "  "
        << std::setw(10) << currentStep->GetPostStepPoint()->GetProcessDefinedStep()
                                        ->GetProcessName();
    } else {
      G4cout << "   UserLimit";
    }

    G4cout << G4endl;
   HepRandom::showEngineStatus();
*/
}

---