G4TankSimulatorOG/G4TankConstruction.h

Go to the documentation of this file.

#ifndef _G4TankSimulatorOG_G4TankConstruction_h_
#define _G4TankSimulatorOG_G4TankConstruction_h_

#include "G4TankSimulator.h"
#include "G4TankPMTAction.h"

#include <utl/Point.h>
#include <utl/TabulatedFunction.h>

// Geant 4 header files
#include <G4VUserDetectorConstruction.hh>
#include <G4Box.hh>
#include <G4Colour.hh>
#include <G4Element.hh>
#include <G4ElementTable.hh>
#include <G4LogicalBorderSurface.hh>
#include <G4LogicalVolume.hh>
#include <G4Material.hh>
#include <G4MaterialTable.hh>
#include <G4PVPlacement.hh>
#include <G4SDManager.hh>
#include <G4Sphere.hh>
#include <G4ThreeVector.hh>
#include <G4Tubs.hh>
#include <G4VisAttributes.hh>
#include <G4GeometryManager.hh>
#include <G4PhysicalVolumeStore.hh>
#include <G4LogicalVolumeStore.hh>
#include <G4SolidStore.hh>
#include <G4RunManager.hh>
#include <CLHEP/Units/SystemOfUnits.h>
#include <G4SubtractionSolid.hh>

#include "G4TankOpBoundaryProcess.h"
#include <G4Ellipsoid.hh>

namespace G4TankSimulatorOG {

  class G4TankConstruction : public G4VUserDetectorConstruction {

  public:

    G4TankConstruction();
    ~G4TankConstruction();

    G4VPhysicalVolume* Construct(void);
    void GetDataForThisTank(void);

  private:

    friend class G4TankFastCerenkov;
    void SetRequiredParameters(void);
    void SetDetectorParameters(void);
    void SetXMLParameters(void);
    void CreateElements(void);
    void CreateMaterials(void);
    void CleanupFreestore(void);
    G4VPhysicalVolume *CreateTank(void);
    
    void CreateAir(void);
    void CreateWater(void);
    void CreateVacuum(void);
    void CreatePyrex(void);
    void CreatePyrex1(void);
    void CreateLucite(void);
    void CreateInterface(void);
    void CreateHDPE(void);
    void CreateLiner(void);
    void CreateAluminium(void);
    void CreateDirt(void);
    void CreatePrimitives(void);
    void CreateHall(void);
    void AssembleTank(void);

    bool fElecBoxEnable, fGroundEnable, fSolarPanelEnable;
  
    G4Element *elN, *elO, *elAl, *elFe, *elH, *elSi, *elB, *elNa, *elC;
    G4Material *Air, *Water, *Vacuum, *SiO2, *B2O2, *Na2O, *Pyrex, *Pyrex1, *CH3;
    G4Material *CH2, *C, *CO2, *Lucite, *Interface, *HDPE, *Al, *Dirt;
    G4MaterialPropertiesTable *airMPT, *waterMPT, *pmtfaceMPT, *pmtfaceMPT1, *pmtdomeMPT;
    G4MaterialPropertiesTable *interfaceMPT, *linerMPT;
    G4MaterialPropertiesTable *linerOpticalMPT;
    G4OpticalSurface *OpLinerSurface;
    G4LogicalBorderSurface *topsurface, *bottomsurface, *sidesurface;
    
    G4Box *expHall_box;
    G4Box *solarPanel_solid, *elecBox_solid;
    G4Tubs *tank_solid, *top_solid, *side_solid, *ground_solid;
    G4Ellipsoid *inner_solid, *interface_in_aux, *interface_out_aux; 
    G4Ellipsoid *pmt_aux, *pmt_aux1, *dome_in_aux, *dome_out_aux; 

    G4SubtractionSolid *pmt_solid, *pmt_solid1, *interface_solid, *dome_solid;
    
    G4LogicalVolume *expHall_log, *tank_log, *top_log, *bottom_log, *side_log;
    G4LogicalVolume *inner_log, *inner1_log, *inner2_log, *inner3_log;
    G4LogicalVolume *pmt1_log, *pmt2_log, *pmt3_log, *interface1_log, *interface2_log;
    G4LogicalVolume *pmt1_log1, *pmt2_log1, *pmt3_log1;
    G4LogicalVolume *interface3_log, *dome1_log, *dome2_log, *dome3_log;
    G4LogicalVolume *ground_log;
    G4LogicalVolume *solarPanel_log, *elecBox_log;
    
    G4VPhysicalVolume *expHall_phys, *tank_phys, *top_phys, *bottom_phys;
    G4VPhysicalVolume *side_phys, *inner1_phys, *inner2_phys, *inner3_phys;
    G4VPhysicalVolume *pmt1_phys, *pmt2_phys, *pmt3_phys, *interface1_phys;
    G4VPhysicalVolume *pmt1_phys1, *pmt2_phys1, *pmt3_phys1;
    G4VPhysicalVolume *interface2_phys, *interface3_phys, *dome1_phys;
    G4VPhysicalVolume *dome2_phys, *dome3_phys, *ground_phys;
    G4VPhysicalVolume *solarPanel_phys, *elecBox_phys;
    
    static double fFaceRadius;
    static double fFaceRadiusz;
    static double fFaceActiveRadius;

    static double fExpHall_x, fExpHall_y, fExpHall_z;
    static double fGroundThickness;
    
    static double fTankRadius;
    static double fTankHalfHeight;
    static double fTankThickness;
    static double fTankPos_x, fTankPos_y, fTankPos_z;

    static utl::Point fPmt1, fPmt2, fPmt3;

    static double fPmtRmin, fPmtRmax;
    static double fInterfaceRmin, fInterfaceRmax, fInterfaceThickness;
    static double fGlassThickness, fDomeThickness;
    static double fDomeRmin, fDomeRmax;
    static double fMinPhi, fMaxPhi;
    static double fMinTheta, fMaxTheta;
    static double fPmtRzmin;
    static double fPmtRzmax;
    static double fInterfaceRzmin;
    static double fInterfaceRzmax;
    static double fDomeRzmin;
    static double fDomeRzmax;

    static double alpha;
    static double fHeightz;
    
    static utl::TabulatedFunction fPmtfaceRINDEX, fInterfaceRINDEX;
    static utl::TabulatedFunction fPmtdomeRINDEX;
    static utl::TabulatedFunction fPmtfaceABSORPTION, fInterfaceABSORPTION;
    static utl::TabulatedFunction fPmtdomeABSORPTION;
    static double fSIGMA_ALPHA;  
    static utl::TabulatedFunction fLinerREFLECTIVITY;
    static utl::TabulatedFunction fLinerSPECULARLOBECONSTANT;
    static utl::TabulatedFunction fLinerSPECULARSPIKECONSTANT;
    static utl::TabulatedFunction fLinerTYVEK_RINDEX;
    static utl::TabulatedFunction fLinerABSORPTION;
    static utl::TabulatedFunction fLinerBACKSCATTERCONSTANT;
    static utl::TabulatedFunction fWaterABSORPTION, fWaterRINDEX;
    
    // Solar panel setup
    //------------------
    
    static double fSolarPanelLength;
    static double fSolarPanelWidth;
    static double fSolarPanelThickness;
    
    static double fSolarPanelX, fSolarPanelY, fSolarPanelZ;
    static double fSolarPanelTiltAngle;    
    
    
    // electronics box setup
    // ---------------------
    
    static double fElecBoxLength;
    static double fElecBoxWidth;
    static double fElecBoxThickness;
    
    static double fElecBoxX, fElecBoxY, fElecBoxZ;
    static double fElecBoxTiltAngle;    
    
    
    // PMT setup
    //----------
    /*
      The "PMT setup" consists of the dome, 
      a thin interface layer, and the face.
      The pmt domes are hemispherical
      shells with the same surface area as the
      the actual domes.
      The pmt "faces" are hemispherical domes 
      that have the same active area as the 
      photocathodes.
      The origin of the following
      coordinates is the center of the 
      tank -> (0.0, 0.0, 60.0) cm 
      Why? Because in order to have the 
      optical properties of the dome, etc.
      work correctly, the volumes must 
      be defined and placed in the tank volume
      (as opposed to the experimental hall volume)      
      Since the dome, interface, and face are 
      defined as hemispherical shells placed in
       one another, they have they same
       x-y-z coordinates.
    */
        
    // PMT Optical properties

    /*
      Tyvek Refractive Index on liner surface vs. momentum
      (NOTE: refractive index is set to zero for the 
      following reason.  The only surfaces that work 
      with the unified model are ground and
      groundbackpainted (others are limited to either
      purely specular or Lambertian reflection.)  Both
      ground and groundbackpainted necessarily calculate
      transmissin coefficients based on index of refraction
      of the next material.  For groundbackpainted, you
      are allowd to define a 'thin' layer of dielectric
      between the ground surface and the painted one.   
      If this is defined to be zero, then all the photons
      undergo total internal reflection from the ground
      surface and the unfied model is used.  Note that
      this is a kludge, and the author of this code 
      needs to be contacted to discuss providing some
      additional flexibility in the optical model.
    */
  };

}
#endif