OMSimSensitiveDetector.hh

#pragma once
 
#include "G4VSensitiveDetector.hh"
#include "G4ThreeVector.hh"
#include "OMSimPMTResponse.hh"
#include "OMSimOpBoundaryProcess.hh"
#include <vector>
 
class G4Step;
class G4TouchableHistory;
 
 
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; 
};
 
 
 
class OMSimSensitiveDetector : public G4VSensitiveDetector
{
public:
    OMSimSensitiveDetector(G4String pName, DetectorType pDetectorType);
    ~OMSimSensitiveDetector();
 
    G4bool ProcessHits(G4Step *pStep, G4TouchableHistory *pTouchableHistory) override;
    void setPMTResponse(OMSimPMTResponse *pResponse);
 
private:
    bool m_QEcut;
    OMSimPMTResponse *m_PMTResponse;
    DetectorType m_detectorType;
    thread_local static 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);
};