FDsimG4MercedesSD.cc

Go to the documentation of this file.

#include "FDsimG4MercedesSD.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;

FDsimG4MercedesSD::FDsimG4MercedesSD(G4String name):G4VSensitiveDetector(name)
{
  collectionName.insert("MercedesHits");
}


FDsimG4MercedesSD::~FDsimG4MercedesSD()
{
}


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

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


G4bool FDsimG4MercedesSD::ProcessHits(G4Step* aStep,G4TouchableHistory*)
{ 
  G4TouchableHistory* theTouchable
    = (G4TouchableHistory*)(aStep->GetPreStepPoint()->GetTouchable());
  
  G4String thisVolume = theTouchable->GetVolume()->GetName();  
  G4String particleName = aStep->GetTrack()->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 globalTOF = aStep->GetTrack()->GetGlobalTime();
  G4double kineticEnergy = aStep->GetTrack()->GetKineticEnergy();
  G4double wavelength =  h_Planck*c_light/kineticEnergy; 
  G4StepPoint* preStepPoint = aStep->GetPreStepPoint();
  G4ThreeVector worldPos = preStepPoint->GetPosition();
  G4ThreeVector localPos = 
    theTouchable->GetHistory()->GetTopTransform().TransformPoint(worldPos);
  G4ThreeVector HitPosition = localPos;
  G4ThreeVector HitDirection = aStep->GetPreStepPoint()->GetMomentumDirection();
  G4ThreeVector SDposition = fPhysVolume->GetObjectTranslation(); 

  FDsimG4OpticalHit* Hit = new FDsimG4OpticalHit;

  Hit->SetEnergy(kineticEnergy);
  Hit->SetWavelength(wavelength);
  Hit->SetTime(globalTOF);
  Hit->SetPosition(HitPosition);
  Hit->SetWorldPosition(worldPos);
  Hit->SetDirection(HitDirection);
  Hit->SetPixelPosition(SDposition);
  Hit->SetHitID(HitID);
  Hit->SetPMTid(copyNo);

  MercedesHitsCollection->insert(Hit);

  return true;
 }


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


void FDsimG4MercedesSD::clear()
{
} 


void FDsimG4MercedesSD::DrawAll()
{
} 


void FDsimG4MercedesSD::PrintAll()
{
}