OpDetSensitiveDetector.cxx

Go to the documentation of this file.

//
// Implementation of the OpDetSensitiveDetector
//
// See comments in OpDetSensitiveDetector.h
//
// Ben Jones, MIT, 06/04/2010
//



#include "larsim/LArG4/OpDetSensitiveDetector.h"
#include "larsim/LArG4/OpDetPhotonTable.h"
#include "larsim/LArG4/OpDetLookup.h"
#include "lardataobj/Simulation/SimPhotons.h"
#include "Geant4/G4SDManager.hh"

namespace larg4{


  OpDetSensitiveDetector::OpDetSensitiveDetector(G4String DetectorUniqueName)
    : G4VSensitiveDetector(DetectorUniqueName)
  {
    // Register self with sensitive detector manager
    G4SDManager::GetSDMpointer()->AddNewDetector(this);

    // Get instances of singleton classes
    fTheOpDetLookup        = OpDetLookup::Instance();
    fThePhotonTable        = OpDetPhotonTable::Instance();
    
  }

  
  //--------------------------------------------------------

  G4bool OpDetSensitiveDetector::ProcessHits(G4Step * aStep, G4TouchableHistory *)
  {
    sim::OnePhoton ThePhoton;
    
    
    // Get photon data to store in the hit

    ThePhoton.SetInSD      = true;

    ThePhoton.InitialPosition     = TVector3(
                                              aStep->GetTrack()->GetVertexPosition().x(),
                                              aStep->GetTrack()->GetVertexPosition().y(),
                                              aStep->GetTrack()->GetVertexPosition().z()        
                                              );
    
    //ThePhoton.Time                = aStep->GetTrack()->GetGlobalTime() - fGlobalTimeOffset;
    ThePhoton.Time                = aStep->GetTrack()->GetGlobalTime();
    

    ThePhoton.Energy              =  aStep->GetTrack()->GetVertexKineticEnergy();

    // Lookup which OpDet we are in
    G4StepPoint *preStepPoint = aStep->GetPreStepPoint();

    int OpDet = fTheOpDetLookup->GetOpDet(preStepPoint->GetPhysicalVolume());

    // Store relative position on the photon detector
    G4ThreeVector worldPosition  = preStepPoint->GetPosition();
    G4ThreeVector localPosition  = preStepPoint->GetTouchableHandle()->GetHistory()->GetTopTransform().TransformPoint(worldPosition);
    ThePhoton.FinalLocalPosition = TVector3(localPosition.x()/CLHEP::cm, localPosition.y()/CLHEP::cm, localPosition.z()/CLHEP::cm);
    

    // Add this photon to the detected photons table
    fThePhotonTable->AddPhoton(OpDet, std::move(ThePhoton));
    
    // Kill this photon track
    aStep->GetTrack()->SetTrackStatus(fStopAndKill);

    return true;
    

  }

  //--------------------------------------------------------
  
  void OpDetSensitiveDetector::Initialize(G4HCofThisEvent *)
  {

  }

}