SAXSSteppingAction Class Reference

Stepping Action. More...

#include <Doxymodules_exoticphysics.h>

Inheritance diagram for SAXSSteppingAction:

Public Member Functions
	SAXSSteppingAction (SAXSEventAction *)
virtual	~SAXSSteppingAction ()
virtual void	UserSteppingAction (const G4Step *)

Private Attributes
SAXSEventAction *	fEventAction
G4LogicalVolume *	fPhantom
G4int	fEventNumber
G4int	fNSe

Detailed Description

Stepping Action.

It is used to retrive information about every interaction of the photons occuring iside the phantom, and in particular, scattering events. The total number of scattering events (subdivided by tyoe) is passed to the EventAction class and stored at the end of the event.

Definition at line 75 of file Doxymodules_exoticphysics.h.

Constructor & Destructor Documentation

SAXSSteppingAction()

SAXSSteppingAction::SAXSSteppingAction

(

SAXSEventAction *

eventAction

)

Definition at line 62 of file SAXSSteppingAction.cc.

                                                                  :
  G4UserSteppingAction(),
  fEventAction(eventAction),
  fPhantom(0),
  fEventNumber(-1),
  fNSe(0)
{}

~SAXSSteppingAction()

SAXSSteppingAction::~SAXSSteppingAction ( )

virtual

Definition at line 72 of file SAXSSteppingAction.cc.

{}

Member Function Documentation

UserSteppingAction()

void SAXSSteppingAction::UserSteppingAction ( const G4Step *  step )

virtual

Definition at line 76 of file SAXSSteppingAction.cc.

{        
  //DetectorConstruction instance
  const SAXSDetectorConstruction* detectorConstruction = 
    static_cast<const SAXSDetectorConstruction*>
    (G4RunManager::GetRunManager()->GetUserDetectorConstruction());
        
  //get Phantom volume
  fPhantom = detectorConstruction->GetPhantom();
 
  //get pre and post step points
  G4StepPoint* preStepPoint = step->GetPreStepPoint();
  G4StepPoint* postStepPoint = step->GetPostStepPoint();
 
  //get the volume of the current step
  G4LogicalVolume* volume = preStepPoint->GetTouchableHandle()
    ->GetVolume()->GetLogicalVolume();
          
  //get track and particle name
  G4Track* track = step->GetTrack();
  G4int ID = track->GetTrackID();
  G4String partName = track->GetDefinition()->GetParticleName();
 
  //update the primary particle (event) weight
  G4double PrimaryWeight = 1.;
  if (partName == "gamma" && ID==1) {
    PrimaryWeight = track->GetWeight();
    fEventAction->UpdateEventWeight(PrimaryWeight);
  }
                                         
  //if the particle is not a primary photon stepping in the phantom, exit!
  if (volume != fPhantom || partName != "gamma" || ID!=1) return;
                
  //if it is a new event, reset variable values to 0
  G4int eventNumber = G4RunManager::GetRunManager()
    ->GetCurrentEvent()->GetEventID();
  if (eventNumber != fEventNumber) {
    fEventNumber = eventNumber;
    fNSe = 0;
  }
 
  //identify the process that occurred
  G4String Process; 
  if (postStepPoint->GetProcessDefinedStep() != NULL) {
    Process = postStepPoint->GetProcessDefinedStep()->GetProcessName();
  } else {
    Process = "UserLimit";
  }
  //G4cout << "Process: " << Process << G4endl;
  G4int ProcIndex = -1;
  if (Process == "Transportation" || Process == "CoupledTransportation" 
      || Process == "UserLimit") {
    ProcIndex = 0;
  }
  if (Process == "Rayl" || Process == "biasWrapper(Rayl)") {
    ProcIndex = 1;
    fNSe++;
    fEventAction->AddNRi();
  }
  if (Process == "compt" || Process == "biasWrapper(compt)") {
    ProcIndex = 2;
    fNSe++;
    fEventAction->AddNCi();
  }
  if (Process == "phot" || Process == "biasWrapper(phot)") ProcIndex = 3;
  if (Process == "conv" || Process == "biasWrapper(conv)") ProcIndex = 4;        
  if (Process == "photonNuclear" ||
      Process == "biasWrapper(photonNuclear)") {
    ProcIndex = 5;
  }
  if (Process == "PowderDiffraction" ||
      Process == "biasWrapper(PowderDiffraction)") {
    ProcIndex = 6;
    fNSe++;
    fEventAction->AddNDi();
  }
  //G4cout << "ProcIndex: " << ProcIndex << G4endl;
        
  //calculate the scattering angle
  G4ThreeVector mom1 = preStepPoint->GetMomentumDirection();
  G4ThreeVector mom2 = postStepPoint->GetMomentumDirection();        
  G4double theta = mom1.angle(mom2); 
 
  //fill the Scattering ntuple 
  G4AnalysisManager* analysisManager = G4AnalysisManager::Instance();
  if (theta > 1e-6) { //record only if the angle is not 0
    analysisManager->FillNtupleIColumn(1,0,ProcIndex);
    analysisManager->FillNtupleDColumn(1,1,
        preStepPoint->GetKineticEnergy()/CLHEP::keV);
    analysisManager->FillNtupleDColumn(1,2,theta/CLHEP::deg); 
    analysisManager->FillNtupleDColumn(1,3,track->GetWeight());    
    analysisManager->AddNtupleRow(1);
  }        
}

Member Data Documentation

fEventAction

SAXSEventAction* SAXSSteppingAction::fEventAction

private

Definition at line 58 of file SAXSSteppingAction.hh.

fPhantom

G4LogicalVolume* SAXSSteppingAction::fPhantom

private

Definition at line 60 of file SAXSSteppingAction.hh.

fEventNumber

G4int SAXSSteppingAction::fEventNumber

private

Definition at line 62 of file SAXSSteppingAction.hh.

fNSe

G4int SAXSSteppingAction::fNSe

private

Definition at line 63 of file SAXSSteppingAction.hh.

---

The documentation for this class was generated from the following files:

• .doxygen/Doxymodules_exoticphysics.h
• extended/exoticphysics/saxs/include/SAXSSteppingAction.hh
• extended/exoticphysics/saxs/src/SAXSSteppingAction.cc