FDsimG4PMT_SD.cc

Go to the documentation of this file.

#include "FDsimG4PMT_SD.hh"
#include "FDsimG4OpBoundaryProcess.hh"
#include "G4RunManager.hh"
#include "G4Material.hh"
#include "G4Step.hh"
#include "G4VTouchable.hh"
#include "G4TouchableHistory.hh"
#include "G4SDManager.hh"
#include "G4ios.hh"

#include <sstream>

using namespace TelescopeSimulatorLX;

FDsimG4PMT_SD::FDsimG4PMT_SD(G4String name):G4VSensitiveDetector(name)
{
  collectionName.insert("OpticalHits");
}


FDsimG4PMT_SD::~FDsimG4PMT_SD()
{
}


void FDsimG4PMT_SD::Initialize(G4HCofThisEvent* HCE)
{
  static int HCID1 = -1;
  fOpticalHitsCollection = 
    new FDsimG4OpticalHitsCollection(SensitiveDetectorName,collectionName[0]);

  if (HCID1 < 0) 
    HCID1 = GetCollectionID(0); 
  HCE->AddHitsCollection(HCID1,fOpticalHitsCollection);
}


G4bool FDsimG4PMT_SD::ProcessHits(G4Step* aStep,G4TouchableHistory*)
{ 
  G4bool status = false;
  G4String particleName = aStep->GetTrack()->GetDefinition()->GetParticleName();

  G4TouchableHistory* theTouchable
    = (G4TouchableHistory*)(aStep->GetPreStepPoint()->GetTouchable());

  FDsimG4OpBoundaryProcessStatus boundaryStatus = Undefined;
  static FDsimG4OpBoundaryProcess* boundary = NULL;

  if (particleName == "opticalphoton"){
    if (!boundary){
      G4ProcessManager* pm 
        = aStep->GetTrack()->GetDefinition()->GetProcessManager();
      G4int nprocesses = pm->GetProcessListLength();
      G4ProcessVector* pv = pm->GetProcessList();
      G4int i;
      for( i=0;i<nprocesses;i++){
        if ((*pv)[i]->GetProcessName() == "OpBoundary"){
          boundary = (FDsimG4OpBoundaryProcess*)(*pv)[i];
          break;
        }
      }
    }

    boundaryStatus=boundary->GetStatus();

    //Check to see if the partcile was actually at a boundary
    G4StepPoint* posStepPoint = aStep->GetPostStepPoint();
    
    if (posStepPoint->GetStepStatus() == fGeomBoundary) {
      switch(boundaryStatus){
      case Absorption:
        break;
      case Detection:
        status = ProcessOpticalHit(aStep,theTouchable);
        break;
      default:
        break;
      }
    }
  } 
  else
    status = true;

  return status;
}
 

G4bool FDsimG4PMT_SD::ProcessOpticalHit(G4Step* aStep, G4TouchableHistory* theTouchable)
{
  G4Track* track = aStep->GetTrack();
  G4String particleName = track->GetDefinition()->GetParticleName();

  G4int HitID;
  G4VPhysicalVolume* fPhysVolume = theTouchable->GetVolume(); 
  G4String VolumeName = fPhysVolume->GetName();
  G4int copyNo = fPhysVolume->GetCopyNo();

  if (particleName == "opticalphoton")
    HitID = 1;
  else
    HitID = 10;

  G4double fGlobalTOF = track->GetGlobalTime();
  G4double kineticEnergy = track->GetKineticEnergy();
  G4double Wavelength = h_Planck*c_light/kineticEnergy; 
  G4StepPoint* preStepPoint = aStep->GetPreStepPoint();
  G4ThreeVector worldPos = preStepPoint->GetPosition();
  G4ThreeVector HitPosition = 
    theTouchable->GetHistory()->GetTopTransform().TransformPoint(worldPos);
  G4ThreeVector HitDirection = aStep->GetPreStepPoint()->GetMomentumDirection();
  G4ThreeVector SDposition = fPhysVolume->GetObjectTranslation(); 
  G4double weight = track->GetWeight(); 

  FDsimG4OpticalHit* OpticalHit = new FDsimG4OpticalHit;

  OpticalHit->SetEnergy(kineticEnergy);
  OpticalHit->SetWavelength(Wavelength);
  OpticalHit->SetTime(fGlobalTOF);
  OpticalHit->SetPosition(HitPosition);
  OpticalHit->SetWorldPosition(worldPos);
  OpticalHit->SetDirection(HitDirection);
  OpticalHit->SetPixelPosition(SDposition);
  OpticalHit->SetHitID(HitID);
  OpticalHit->SetPMTid(copyNo);
  OpticalHit->SetWeight(weight);

  fOpticalHitsCollection->insert(OpticalHit);

  return true;
}


void FDsimG4PMT_SD::EndOfEvent(G4HCofThisEvent* HCE)
{
  static G4int HCID = -1;
  if (HCID < 0)
    HCID = G4SDManager::GetSDMpointer()->GetCollectionID(collectionName[0]);
    
  HCE->AddHitsCollection(HCID,fOpticalHitsCollection);
}


void FDsimG4PMT_SD::clear()
{
} 


void FDsimG4PMT_SD::DrawAll()
{
} 


void FDsimG4PMT_SD::PrintAll()
{
}