OMSimMDOM.hh

Go to the documentation of this file.

 
#pragma once
#include "OMSimMDOMFlasher.hh"
#include "OMSimPMTConstruction.hh"
#include "OMSimOpticalModule.hh"
 
class mDOMHarness;
 
class mDOM : public OMSimOpticalModule
{
 
public:
    mDOM(G4bool p_placeHarness = true);
    ~mDOM();
    void construction();
    double getPressureVesselWeight() { return 13.0 * kg; };
    int getNumberOfPMTs() { return m_totalNumberPMTs; };
    G4String getName()
    {
        std::stringstream ss;
        ss << "/mDOM/" << m_index;
        return ss.str();
    }
 
    G4int m_NrTotalLED;
    std::vector<G4Transform3D> m_LEDTransformers;           
    std::vector<std::vector<G4double>> m_LEDAngFromSphere; 
 
    void runBeamOnFlasher(G4int pModuleIndex, G4int pLEDIndex) { m_flashers->runBeamOnFlasher(this, pModuleIndex, pLEDIndex); }
 
private:
    mDOMFlasher *m_flashers;
    mDOMHarness *m_harness;
    G4SubtractionSolid *equatorialReflector(G4VSolid *pSupportStructure, G4Cons *pReflCone, G4double pAngle, G4String pSuffix);
    void setPMTPositions();
    G4UnionSolid *pressureVessel(const G4double pOutRad, G4String pSuffix);
    G4SubtractionSolid *substractHarnessPlug(G4VSolid *pSolid);
    std::tuple<G4SubtractionSolid *, G4UnionSolid *> supportStructure();
    G4SubtractionSolid *substractFlashers(G4VSolid *lSupStructureSolid);
    void setLEDPositions();
 
    G4bool m_placeHarness = true;
    G4bool m_harnessUnion = true; // it should be true for the first module that you build, and then false
    std::vector<G4ThreeVector> m_positionsPMT;
    std::vector<G4RotationMatrix> m_PMTRotations;
    std::vector<G4RotationMatrix> m_PMTRotPhi;
    std::vector<G4ThreeVector> m_reflectorPositions;
 
    G4double m_PMToffset;
    G4double m_refConeIdealInRad;
    const G4double m_glassThick = 13.5 * mm; 
 
public:
    const G4double m_cylinderAngle = 2.8 * deg; 
    const G4double m_glassOutRad = 176.5 * mm;  
    const G4double m_cylinderHeight = 27.5 * mm;       
    const G4double m_glassInRad = m_glassOutRad - m_glassThick;
 
private:
    const G4double m_gelThicknessFrontPMT = 3.6 * mm;  
    const G4double m_gelThickness = 4.5 * mm;          
    const G4double m_EqPMTrOffset = 2.6 * mm;          
    const G4double m_EqPMTzOffset = 10.0 * mm;         
    const G4double m_reflectorHalfZ = 15 * mm;           
    const G4double m_reflectorConeSheetThickness = 0.5 * mm; 
    const G4double m_reflectorConeToHolder = 1.55 * mm;      
    const G4double m_thetaPolar = 33.0 * deg;
    const G4double m_thetaEquatorial = 72.0 * deg;
    const G4int m_numberPolarPMTs = 4;
    const G4int m_numberEqPMTs = 8;
    const G4int mRefConeAngle = 51;
    const G4double m_polarEquatorialPMTphiPhase = 0 * deg;
    const G4int m_totalNumberPMTs = (m_numberPolarPMTs + m_numberEqPMTs) * 2;
    const G4double m_supStructureRad = m_glassOutRad - m_glassThick - m_gelThickness;
 
    const G4double m_thetaEqLED = 61 * deg;   // 61 upper sphere, 180-61 lower sphere
    const G4double m_thetaPolLED = 8.2 * deg; // 8.2 upper sphere, 180-8.2 lower sphere
};