simulations/wavepid/include/OMSimSensitiveDetector.hh

Go to the documentation of this file.

 
#pragma once
 
#include "G4VSensitiveDetector.hh"
#include "G4ThreeVector.hh"
#include "OMSimPMTResponse.hh"
 
#include <vector>
#include <string>
 
class G4Step;
class G4Track;
class G4TouchableHistory;
class G4OpBoundaryProcess;
 
enum class DetectorType {
    PMT,
    VolumePhotonDetector,
    BoundaryPhotonDetector,
    BoundaryShellDetector,
    PerfectPMT
};
 
struct PhotonInfo {
    G4int eventID;
    G4double globalTime;
    G4double localTime;
    G4double trackLength;
    G4double kineticEnergy;
    G4double wavelength;
    G4ThreeVector globalPosition;
    G4ThreeVector localPosition;
    G4ThreeVector momentumDirection;
    G4ThreeVector deltaPosition;
    G4int pmtNumber;
    G4int detectorID;
    OMSimPMTResponse::PMTPulse PMTResponse;
    // WavePID extensions for photon origin tracking
    G4String photonOrigin;      
    G4int parentID;             
    std::string parentType;     
    G4String parentProcess;     
    G4double entryTime;         
};
 
class OMSimSensitiveDetector : public G4VSensitiveDetector
{
public:
    OMSimSensitiveDetector(G4String pName, DetectorType pDetectorType);
    ~OMSimSensitiveDetector() override;
 
    G4bool ProcessHits(G4Step* pStep, G4TouchableHistory* pTouchableHistory) override;
    void setPMTResponse(OMSimPMTResponse* pResponse);
 
private:
    bool m_QEcut;
    OMSimPMTResponse* m_PMTResponse;
    DetectorType m_detectorType;
 
    static thread_local G4OpBoundaryProcess* m_boundaryProcess;
 
    G4bool checkVolumeAbsorption(G4Step* pStep);
    G4bool checkBoundaryAbsorption(G4Step* pStep);
    PhotonInfo getPhotonInfo(G4Step* pStep);
    G4bool handlePMT(G4Step* pStep, G4TouchableHistory* pTouchableHistory);
    G4bool handleGeneralPhotonDetector(G4Step* pStep, G4TouchableHistory* pTouchableHistory);
    G4bool handleShellDetector(G4Step* pStep, G4TouchableHistory* pTouchableHistory);
    bool isPhotonDetected(double p_efficiency);
    void storePhotonHit(PhotonInfo& pInfo);
    void fetchBoundaryProcess();
    void killParticle(G4Track* pTrack);
};