src/SensitiveDetector.cc

Go to the documentation of this file.

#include "SensitiveDetector.hh"
#include "DetectorConstruction.hh"
#include "RunAction.hh"

#include "G4Step.hh"
#include "Randomize.hh"
#include "G4RunManager.hh"

SensitiveDetector::SensitiveDetector(G4String SDname)
  : G4VSensitiveDetector(SDname)
{
  // -- obtain detector construction and run action pointers through the run manager
  m_detector = (DetectorConstruction*)(G4RunManager::GetRunManager()->GetUserDetectorConstruction());
  m_runAction = (RunAction*)(G4RunManager::GetRunManager()->GetUserRunAction());
}

SensitiveDetector::~SensitiveDetector()
{}

G4bool SensitiveDetector::ProcessHits(G4Step *step, G4TouchableHistory *)
{
  // ------------------------------- for first demo need (not to be kept !) ---------------------
  G4cout << std::setw( 5) << step->GetTrack()->GetCurrentStepNumber()
         << " ->> Sensitive detector `" << GetName() << "' ProcessHits(...) method called." << G4endl;
  // --------------------------------------------------------------------------------------------

  // -----------------
  // -- skeleton code:
  // -----------------

  // energy deposit in this step 
  G4double edep = step->GetTotalEnergyDeposit();

  if (edep <= 0.) return false;
  
  // sum edep in RunAction (in SD vector):
  // -->> this is a new method of RunAction for this exercise :  m_runAction->SumDepthDoseSD(zLocal,edep);

  return true;


  // ------------------------------------------------------------------------
  // --                    Recalling of SteppingAction code                --
  // ------------------------------------------------------------------------
  //  o Adapt this code to the example need.
  //  o Note : the pointers to the DetectorConstruction and RunAction
  //    objects, m_detector and m_runAction, respectively are provided to you
  //    (see this class constructor code).
  /* --->>

  // in which volume is the step ?
  G4TouchableHandle touchable = step->GetPreStepPoint()->GetTouchableHandle();
  G4VPhysicalVolume* volume   = touchable->GetVolume();
  
  // energy deposit in this step 
  G4double edep = step->GetTotalEnergyDeposit();
  
  // is astronaut ?
  if (volume == m_detector->GetAstronaut()) {
    ////G4cout << "--> edep in astronaut : " << edep/keV << " keV" << G4endl;
    
    // sum edep in EventAction
    m_eventAction->AddEdep(edep);
    
    //---------for depth dose distribution---------------------------------------
    
    // get step points in world coordinate system
    G4ThreeVector point1 = step->GetPreStepPoint()->GetPosition();
    G4ThreeVector point2 = step->GetPostStepPoint()->GetPosition();
    
    // randomize point of energy deposition
    G4ThreeVector pointE = point1 + G4UniformRand()*(point2 - point1);      
    // transform it in local coordinate system
    G4ThreeVector localPointE 
      = touchable->GetHistory()->GetTopTransform().TransformPoint(pointE);
    
    // extract z and offset
    G4double zLocal = 0.5*(m_detector->GetAstronautHeight()) - localPointE.z();
    
    // sum edep in RunAction
    m_runAction->SumDepthDose(zLocal,edep);
    
    //----------end of depth dose distribution-------------------------------------
  }

  <<--- */

}

---