FDsimG4MirrorSD.cc

Go to the documentation of this file.

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

//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
FDsimG4MirrorSD::FDsimG4MirrorSD(G4String name):G4VSensitiveDetector(name)
{
  collectionName.insert("MirrorHits");

  G4cerr << "Created " << name << " Sensitive Detector." << G4endl ;
}


//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....

FDsimG4MirrorSD::~FDsimG4MirrorSD()
{;}

//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....

void FDsimG4MirrorSD::Initialize(G4HCofThisEvent* HCE)
{

  static int HCID1 = -1;

  // SensitiveDetectorName e collectionName sao data members de G4VSensitiveDetector

  MirrorHitsCollection =
    new FDsimG4OpticalHitsCollection(SensitiveDetectorName,collectionName[0]);

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

}

//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....

G4bool FDsimG4MirrorSD::ProcessHits(G4Step* aStep,G4TouchableHistory*){

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


  // Get Material

  G4String thisVolume = theTouchable->GetVolume()->GetName() ;
  G4String particleName = aStep->GetTrack()->GetDefinition()->GetParticleName();


#ifdef VERBOSE_USE
  G4cerr << "Volume name: " << thisVolume << G4endl ;;
#endif


  G4int HitID ;

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

  // G4cerr << "Volume name: " << VolumeName << " " << copyNo << G4endl ;

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

  G4double      fGlobalTOF    = aStep->GetTrack()->GetGlobalTime()       ;
  G4double      kineticEnergy = aStep->GetTrack()->GetKineticEnergy()    ;
  G4double      Wavelength    =  h_Planck*c_light/kineticEnergy ;

  //  G4StepPoint* preStepPoint ;  // unused
  G4StepPoint* posStepPoint ;
  //  preStepPoint = aStep->GetPreStepPoint(); // unused
  posStepPoint = aStep->GetPostStepPoint();

  G4ThreeVector worldPos    = posStepPoint->GetPosition();
  G4ThreeVector localPos    = theTouchable->GetHistory()->GetTopTransform().TransformPoint(worldPos);

  G4ThreeVector HitPosition   = localPos;


  G4ThreeVector HitDirection  = posStepPoint->GetMomentumDirection() ;
  G4ThreeVector SDposition    = fPhysVolume->GetObjectTranslation()      ;


  FDsimG4OpticalHit* Hit = new FDsimG4OpticalHit ;

  Hit->SetEnergy(kineticEnergy);
  Hit->SetWavelength(Wavelength);
  Hit->SetTime(fGlobalTOF);
  Hit->SetPosition(HitPosition);
  Hit->SetWorldPosition(worldPos);
  Hit->SetDirection(HitDirection);
  Hit->SetPixelPosition(SDposition);
  Hit->SetHitID(HitID);

  Hit->SetPMTid(copyNo) ;


  MirrorHitsCollection->insert(Hit);


#ifdef VERBOSE_USE

  G4cerr << fPhysVolume->GetObjectTranslation() << " "
         << fPhysVolume->GetFrameTranslation()  << " "
         << fPhysVolume->GetTranslation()       << G4endl ;

  G4cerr << fPhysVolume->GetMultiplicity() << " "
         << fPhysVolume->GetCopyNo()       << G4endl ;

#endif
  return true ;

}



//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....

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

}

//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....

void FDsimG4MirrorSD::clear()
{;}

//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....

void FDsimG4MirrorSD::DrawAll()
{;}

//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....

void FDsimG4MirrorSD::PrintAll()
{;}

//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....