PMT encapsulation

ratdb/ENCAPSULATION.ratdb

@@ -0,0 +1,30 @@
+{
+name: "ENCAPSULATION", 
+index: "BUTTON_r7081pe",
+run_range: [0, 0],
+
+construction: "hemisphere",
+envelope_radius: 244.0,  // should be large enough to fit all objects (>flange_rmax)
+front_encapsulation_material: "acrylic_uvt_good",
+rear_encapsulation_material: "acrylic_black",
+flange_rmax: 243.0,  // applies to acrylic and metal flange
+encap_radius: 200.0,
+encap_thickness: 8.0,
+use_metal_flange: 1,
+metal_flange_material: "stainless_steel",
+metal_flange_dimensions: [ 220.0, 243.0, 4.0 ],  // G4Tubs dimensions rmin, rmax, z
+inside_encapsulation_material: "air",
+use_optical_gel: 1,
+optical_gel_material: "optical_grease",
+optical_gel_sub_height: 150.0,  // start from a dome full of gel of size encap_radius, e.g. height of ...
+//  gel dome left is (encap_radius - optical_gel_sub_height)
+pmtposoffset: [ 0.0, 0.0, 98.0 ],
+use_silica_bag: 1,
+silica_bag_material: "acrylic_white",
+silica_bag_dimensions: [ 18.0, 33.0, 3.0 ],  // G4Box dimensions x, y, z
+silica_bag_position: [ 0.0, 135.0, -74.0 ],  // x, y, z
+use_cable: 1,
+cable_material: "acrylic_black",
+cable_dimensions: [ 0.0, 6.5, 45.0 ], // G4Tubs dimensions rmin, rmax, z
+cable_position: [ 0.0, 0.0, -152.0 ],  // x, y, z
+}

src/geo/CMakeLists.txt

@@ -64,6 +64,8 @@ add_library(geo OBJECT
         src/pmt/RevolutionPMTConstruction.cc
         src/pmt/ToroidalPMTConstruction.cc
         src/pmt/CylindricalPMTConstruction.cc
+        src/pmt/PMTEncapsulation.cc 
+        src/pmt/HemisphereEncapsulation.cc
         )
 
 # Set our include directories

src/geo/include/RAT/HemisphereEncapsulation.hh

@@ -0,0 +1,106 @@
+// This is the encapsulation model that is used for BUTTON.
+// The 96 Hamamatsu r7081pe PMTs are encapsulated by two acryilic domes that are held together with metal flanges.
+// Created by Lewis Sexton (Sheffield) and Adam Tarrant (Liverpool)
+
+#ifndef __RAT_HemisphereEncapsulation__
+#define __RAT_HemisphereEncapsulation__
+
+#include <G4Box.hh>
+#include <G4GenericPolycone.hh>
+#include <G4LogicalVolume.hh>
+#include <G4Material.hh>
+#include <G4OpticalSurface.hh>
+#include <G4PVPlacement.hh>
+#include <G4Paraboloid.hh>
+#include <G4Sphere.hh>
+#include <G4SubtractionSolid.hh>
+#include <G4UnionSolid.hh>
+#include <G4VSensitiveDetector.hh>
+#include <G4VSolid.hh>
+#include <G4VisAttributes.hh>
+#include <RAT/DB.hh>
+#include <RAT/Factory.hh>
+#include <RAT/GLG4TorusStack.hh>
+#include <RAT/PMTEncapsulation.hh>
+#include <string>
+#include <vector>
+namespace RAT {
+
+struct HemisphereEncapsulationParams {
+  int useMetalFlange;
+  int useGel;
+  int useSilicaBag;
+  int useCable;
+
+  // encapsulation object dimensions
+  double envelope_radius;
+  double encap_radius;
+  double encap_thickness;
+  double flange_rmax;
+  double dome_flange_thickness;
+  double metal_flange_thickness;
+  double optical_gel_sub_height;
+  std::vector<double> metal_flange_dimensions;
+  std::vector<double> silica_bag_dimensions;
+  std::vector<double> silica_bag_position;
+  std::vector<double> cable_dimensions;
+  std::vector<double> cable_position;
+
+  // PMT Body
+  std::vector<double> zEdge;    // n+1
+  std::vector<double> rhoEdge;  // n+1
+  std::vector<double> zOrigin;  // n
+  G4ThreeVector pmtposvec;
+  std::vector<double> posvec;
+
+  G4Material *exterior_material;
+  G4Material *inner_encapsulation_material;
+  G4Material *front_encapsulation_material;
+  G4Material *rear_encapsulation_material;
+  G4Material *metal_flange_material;
+  G4Material *silica_bag_material;
+  G4Material *cable_material;
+  G4Material *optical_gel_material;
+
+  G4OpticalSurface *inner_encapsulation_surface;
+  G4OpticalSurface *front_encapsulation_surface;
+  G4OpticalSurface *rear_encapsulation_surface;
+  G4OpticalSurface *metal_flange_surface;
+  G4OpticalSurface *cable_surface;
+  G4OpticalSurface *silica_bag_surface;
+  G4OpticalSurface *optical_gel_surface;
+};
+
+class HemisphereEncapsulation : public PMTEncapsulation {
+ public:
+  HemisphereEncapsulation(DBLinkPtr encaptable, DBLinkPtr pmttable, G4LogicalVolume *mother);
+  virtual ~HemisphereEncapsulation() {}
+
+  virtual G4LogicalVolume *BuildVolume(const std::string &prefix);
+  virtual G4VSolid *BuildSolid(const std::string &prefix);
+  virtual G4PVPlacement *PlaceEncap(G4RotationMatrix *pmtrot, G4ThreeVector pmtpos, const std::string &name,
+                                    G4LogicalVolume *logi_pmt, G4VPhysicalVolume *mother_phys, bool booleanSolid,
+                                    int copyNo);
+
+ protected:
+  G4VSolid *NewEnvelopeSolid(const std::string &name);
+  G4VSolid *NewEncapsulationSolid(const std::string &name);
+  G4VSolid *optical_gel_height_subtraction(const std::string &_name);
+  G4VSolid *optical_gel_pmt_subtraction(const std::string &_name, GLG4TorusStack *body);
+
+  // phyiscal volumes
+  G4PVPlacement *inner_encapsulation_phys;
+  G4PVPlacement *optical_gel_encapsulation_phys;
+  G4PVPlacement *front_encapsulation_phys;
+  G4PVPlacement *rear_encapsulation_phys;
+  G4PVPlacement *front_metal_encapsulaion_flange_phys;
+  G4PVPlacement *rear_metal_encapsulation_flange_phys;
+  G4PVPlacement *silica_bag_encapsulation_phys;
+  G4PVPlacement *cable_encapsulation_phys;
+
+  HemisphereEncapsulationParams fParams;
+};
+
+}  // namespace RAT
+
+#endif

src/geo/include/RAT/PMTEncapsulation.hh

@@ -0,0 +1,35 @@
+#ifndef __RAT_PMTEncapsulation__
+#define __RAT_PMTEncapsulation__
+
+#include <G4LogicalVolume.hh>
+#include <G4PVPlacement.hh>
+#include <G4RotationMatrix.hh>
+#include <G4VSolid.hh>
+#include <RAT/DB.hh>
+#include <string>
+
+namespace RAT {
+
+class PMTEncapsulation {
+ public:
+  static PMTEncapsulation *NewConstruction(DBLinkPtr encaptable, DBLinkPtr pmttable, G4LogicalVolume *mother);
+
+  PMTEncapsulation(std::string _name) : name(_name) {}
+
+  virtual ~PMTEncapsulation() {}
+
+  virtual G4VSolid *BuildSolid(const std::string &prefix) = 0;
+
+  virtual G4LogicalVolume *BuildVolume(const std::string &prefix) = 0;
+
+  virtual G4PVPlacement *PlaceEncap(G4RotationMatrix *pmtrot, G4ThreeVector pmtpos, const std::string &name,
+                                    G4LogicalVolume *logi_pmt, G4VPhysicalVolume *mother_phys, bool booleanSolid,
+                                    int copyNo) = 0;
+
+ protected:
+  std::string name;
+};
+
+}  // namespace RAT
+
+#endif