musrSteppingAction.cc

Go to the documentation of this file.

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

#include "musrSteppingAction.hh"
#include "G4SteppingManager.hh"
#include "G4UnitsTable.hh"

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

musrSteppingAction::musrSteppingAction()
{ 
 pointer=this;
}

musrSteppingAction::~musrSteppingAction()
{ 
 ;
}


 musrSteppingAction* musrSteppingAction::pointer=0;
 musrSteppingAction* musrSteppingAction::GetInstance()
{
  return pointer;
  }


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

void musrSteppingAction::UserSteppingAction(const G4Step* aStep)
{ 
  G4String p_name;
  if (aStep->GetTrack()->GetDefinition())
    {
    p_name=aStep->GetTrack()->GetDynamicParticle()->GetDefinition()->GetParticleName();

    if (p_name == "mu+")
      {
    G4ThreeVector position = aStep->GetPostStepPoint()->GetPosition();
    //    G4double tof=aStep->GetTrack()->GetLocalTime();
    //    G4double energy=aStep->GetTrack()->GetDynamicParticle()->GetKineticEnergy();
    G4ThreeVector polarization=aStep->GetTrack()->GetDynamicParticle()->GetPolarization();

    //      G4cout << "Stepping action: mu+ properties \n"
    //       << "position "      << G4BestUnit(position,"Length") <<"; \n"
    //        << "time of flight " << G4BestUnit(tof,"Time")        <<"; \n"
    //        << "kinetic energy  "<< G4BestUnit(energy,"Energy")   <<"; \n"
    //        <<  "polarization   "<< polarization <<";\n"
    //       <<G4endl;
      }

    }

}

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

---