From f642bc8d34dc97dc1eb0e45e52c5aeb7d8ff547b Mon Sep 17 00:00:00 2001
From: Chao Peng <cpeng@anl.gov>
Date: Mon, 23 Jan 2023 22:37:50 -0600
Subject: [PATCH] move gas cherenkov plugins to this repo
---
CMakeLists.txt | 2 +
compact/HGC.xml | 2 +-
compact/LGC.xml | 2 +-
src/HeavyGasCherenkov_geo.cpp | 257 +++++++++++++++++++++++++++++
src/ThresholdGasCherenkov_geo.cpp | 259 ++++++++++++++++++++++++++++++
views/top_view/.DAWN_1.history | 8 +-
6 files changed, 524 insertions(+), 6 deletions(-)
create mode 100644 src/HeavyGasCherenkov_geo.cpp
create mode 100644 src/ThresholdGasCherenkov_geo.cpp
diff --git a/CMakeLists.txt b/CMakeLists.txt
index 5f10882..dff6aa5 100644
--- a/CMakeLists.txt
+++ b/CMakeLists.txt
@@ -29,6 +29,8 @@ dd4hep_add_plugin(${a_lib_name} SOURCES
src/GEMTrackerDiscSOLID_geo.cpp
src/NH3Target_geo.cpp
src/StandInGeometrySOLID_geo.cpp
+ src/HeavyGasCherenkov_geo.cpp
+ src/ThresholdGasCherenkov_geo.cpp
USES ActsCore ActsPluginDD4hep
)
target_link_libraries(${a_lib_name}
diff --git a/compact/HGC.xml b/compact/HGC.xml
index db7e604..8928b4d 100644
--- a/compact/HGC.xml
+++ b/compact/HGC.xml
@@ -1,6 +1,6 @@
<lccdd>
<detectors>
- <detector id="8" name="HeavyGasCherenkov" type="HeavyGasCherenkov" readout="HeavyGasCherenkovHits" vis="BlueVis" material="N2cherenkov">
+ <detector id="8" name="HeavyGasCherenkov" type="SoLID_HGC" readout="HeavyGasCherenkovHits" vis="BlueVis" material="N2cherenkov">
<placement x="0*cm" y="0*cm" z="370*cm" />
<dimensions rmin1="80*cm" rmin2="94*cm" rmax1="265*cm" rmax2="265*cm" zmin="306*cm" zmax="406*cm"/>
<mirror id="1">
diff --git a/compact/LGC.xml b/compact/LGC.xml
index 9e36717..f65fb4d 100644
--- a/compact/LGC.xml
+++ b/compact/LGC.xml
@@ -1,6 +1,6 @@
<lccdd>
<detectors>
- <detector id="7" name="LightGasCherenkov" type="ThresholdGasCherenkov" readout="LightGasCherenkovHits" vis="RedVis" material="N2cherenkov">
+ <detector id="7" name="LightGasCherenkov" type="SoLID_LGC" readout="LightGasCherenkovHits" vis="RedVis" material="N2cherenkov">
<placement x="0*cm" y="0*cm" z="260*cm" />
<dimensions rmin="10*cm" rmax1="80*cm" rmax2="120*cm" zmin="20*cm" zmax="120*cm"/>
<mirror id="1">
diff --git a/src/HeavyGasCherenkov_geo.cpp b/src/HeavyGasCherenkov_geo.cpp
new file mode 100644
index 0000000..5679761
--- /dev/null
+++ b/src/HeavyGasCherenkov_geo.cpp
@@ -0,0 +1,257 @@
+#include "DD4hep/DetFactoryHelper.h"
+#include "DD4hep/Printout.h"
+#include <XML/Helper.h>
+#include "TMath.h"
+#include "DDRec/Surface.h"
+
+#include "DD4hep/OpticalSurfaces.h"
+#include "DDRec/DetectorData.h"
+
+using namespace std;
+using namespace dd4hep;
+using namespace dd4hep::rec;
+
+using namespace dd4hep;
+
+/** \addtogroup PID Particle ID Detectors
+ */
+/** \addtogroup HGC Heavy Gas threshold Cherenkov ring .
+ * \brief Type: **HeavyGasCherenkov**.
+ * \ingroup PID
+ *
+ * \code
+ * <detector>
+ * </detector>
+ * \endcode
+ * @{
+ */
+static Ref_t createDetector(Detector& desc, xml::Handle_t handle, SensitiveDetector sens) {
+ xml::DetElement detElem = handle;
+ xml_det_t x_det = handle;
+
+ std::string detName = detElem.nameStr();
+ int detID = detElem.id();
+
+ DetElement det(detName, detID);
+ xml::Component dims = detElem.dimensions();
+ double rInner1 = dims.rmin1();
+ double rInner2 = dims.rmin2();
+ double rOuter1 = dims.rmax1();
+ double rOuter2 = dims.rmax2();
+ double zMin = dims.zmin();
+ double zMax = dims.zmax();
+ double tank_length = zMax - zMin;
+
+ std::map<int,Position> mirror_positions;
+ std::map<int,std::array<double,3>> mirror_rotations;
+
+ for(xml_coll_t i(x_det,Unicode("mirror")); i; ++i){
+ xml_comp_t x_mir = i;
+ //std::cout << "mirror " << x_mir.id() << "\n";
+ xml_dim_t mir_pos = x_mir.child(_U(placement));
+ xml_dim_t mir_rot = x_mir.child(_U(rotation));
+ mirror_positions[x_mir.id()] = Position(mir_pos.x(), mir_pos.y(), mir_pos.z());
+ mirror_rotations[x_mir.id()] = {mir_rot.x(), mir_rot.y(), mir_rot.z()};
+ }
+
+ xml_dim_t pos = x_det.child(_U(placement));
+ double pos_x = pos.x();
+ double pos_y = pos.y();
+ double pos_z = pos.z();
+
+ Material air = desc.air();
+ Material PyrexGlass = desc.material("PyrexGlass");
+ Material N2cherenkov = desc.material("N2cherenkov");
+ Material Copper = desc.material("Copper");
+
+ double LGC_inner_radius1 = 71.0*cm;
+ double LGC_inner_radius2 = 85.0*cm;
+ double LGC_outer_radius1 = 265.0*cm;
+ double LGC_main_length = 105.0*cm;
+ double LGC_snout_length = 107.0*cm;
+ double LGC_snout_inner_radius1 = 58.0*cm;
+ double LGC_snout_inner_radius2 = LGC_inner_radius1;
+ double LGC_snout_outer_radius1 = 127.0*cm;
+ double LGC_snout_outer_radius2 = 144.0*cm;
+ double LGC_entrance_window_thickness = 0.05*mm; // something tells this might be 5 mil, not mm
+ double LGC_exit_window_thickness = 0.1*mm; // same here
+
+ double LGC_mirror1_radius = 210.0*cm;
+ double LGC_mirror2_radius = 157.99*cm;
+
+ double LGC_mirror1_length = 114.53*cm;
+ double LGC_mirror2_length = 59.260*cm;
+ double LGC_mirror1_width1 = 16.26*cm;
+ double LGC_mirror1_width2 = 36.03*cm;
+ double LGC_mirror2_width1 = 37.06*cm;
+ double LGC_mirror2_width2 = 45.95*cm;
+ double LGC_mirror1_thickness = 2.0*mm;
+ double LGC_mirror2_thickness = 2.0*mm;
+
+ double LGC_sector_angle = M_PI * 15.0 / 180.0;
+ //double LGC_scattering_angle = 11.0 * M_PI / 180.0;
+
+ double LGC_mirror1_tilt_angle = mirror_rotations[1][0];//25.0 * M_PI / 180.0;
+ double LGC_mirror2_tilt_angle = mirror_rotations[2][0];//2.0 * M_PI / 180.0;
+ double LGC_pmt_tilt_angle = mirror_rotations[3][0];//45.0 * M_PI / 180.0;
+
+ double LGC_pmt_z_pos = mirror_positions[3].z();//-30.0*cm;
+ double LGC_pmt_y_pos = mirror_positions[3].y();//LGC_outer_radius1 - 20.0*cm;
+ double LGC_pmt_array_size = 20.0*cm;
+
+ // the gas tank
+ ConeSegment tank_main(0.5 * tank_length, rInner1, rOuter1,
+ rInner2, rOuter2);
+ Volume v_lgc_tank("v_lgc_tank_gas", tank_main, N2cherenkov);
+
+ v_lgc_tank.setVisAttributes(desc.visAttributes(detElem.visStr()));
+ // Everything that goes in the tank will be copies of the sector assembly volume
+ Assembly v_sector("cherenkov_sector_1");
+ DetElement de_sector("de_sector"+std::to_string(1),1);
+
+
+ // mirrors
+ Sphere mirror1_shell(LGC_mirror1_radius, LGC_mirror1_radius + LGC_mirror1_thickness,
+ 0.0, M_PI / 2);
+ Trd1 mirror1_cutout(LGC_mirror1_width1 / 2.0, LGC_mirror1_width2 / 2.0,
+ LGC_mirror1_length / 2.0, LGC_mirror1_length / 2.0);
+ IntersectionSolid mirror1_shape(mirror1_cutout, mirror1_shell,
+ RotationX(M_PI/2.0)*Transform3D(Position(0, 0, -LGC_mirror1_radius)));
+ //Sphere mirror2_shell(LGC_mirror2_radius, LGC_mirror2_radius + LGC_mirror2_thickness,
+ // 0.0, M_PI / 2);
+ //Trd1 mirror2_cutout(LGC_mirror2_width1 / 2.0, LGC_mirror2_width2 / 2.0,
+ // LGC_mirror2_length / 2.0, LGC_mirror2_length / 2.0);
+ //IntersectionSolid mirror2_shape(mirror2_cutout, mirror2_shell,
+ // RotationX(M_PI/2.0)*Transform3D(Position(0, 0, -LGC_mirror2_radius)));
+
+ double z_mirror1 = mirror_positions[1].z();
+ //double z_mirror2 = mirror_positions[2].z();
+ double y_mirror1 = mirror_positions[1].y();
+ //double y_mirror2 = mirror_positions[2].y();
+
+ Volume v_mirror1_shape("v_mirror1_shape", mirror1_shape, PyrexGlass);
+ PlacedVolume pv_mirror1_shape = v_sector.placeVolume(
+ v_mirror1_shape, Transform3D(Position(0, y_mirror1, z_mirror1)) *
+ RotationX(-M_PI / 2.0 + LGC_mirror1_tilt_angle));
+
+ DetElement de_mirror1_shape(det,"de_mirror1_shape"+std::to_string(1),1);
+ pv_mirror1_shape.addPhysVolID("mirror", 1);
+ de_mirror1_shape.setPlacement(pv_mirror1_shape);
+ sens.setType("photoncounter");
+ v_mirror1_shape.setSensitiveDetector(sens);
+
+ //Volume v_mirror2_shape("v_mirror2_shape", mirror2_shape, PyrexGlass);
+ //PlacedVolume pv_mirror2_shape = v_sector.placeVolume(
+ // v_mirror2_shape, Transform3D(Position(0, y_mirror2, z_mirror2)) *
+ // RotationX(-M_PI / 2.0 + LGC_mirror2_tilt_angle));
+
+ //DetElement de_mirror2_shape(det,"de_mirror2_shape"+std::to_string(2),2);
+ //pv_mirror2_shape.addPhysVolID("mirror", 2);
+ //de_mirror2_shape.setPlacement(pv_mirror2_shape);
+ //sens.setType("photoncounter");
+ //v_mirror2_shape.setSensitiveDetector(sens);
+
+ // ---------------
+ // Winston Cone
+ double LGC_winston_cone_thickness = 4*mm;
+ double LGC_winston_tube_inner_radius = 11.28*cm;
+ double LGC_winston_tube_length = 30.0*cm;
+ double LGC_winston_cone_length = 30.0*cm;
+ double LGC_winston_cone_inner_radius1 = 7.8*cm;
+ double LGC_winston_cone_inner_radius2 = 21.0*cm;
+ double LGC_winston_cone_inset_length = 7.90909*cm;
+ DetElement de_winston_cone(det,"de_winston_cone1",1);
+ Tube winston_tube(LGC_winston_tube_inner_radius,
+ LGC_winston_tube_inner_radius + LGC_winston_cone_thickness,
+ LGC_winston_tube_length / 2.0);
+ //Cone winston_cone(LGC_winston_cone_length / 2.0, LGC_winston_cone_inner_radius1,
+ // LGC_winston_cone_inner_radius1 + LGC_winston_cone_thickness,
+ // LGC_winston_cone_inner_radius2,
+ // LGC_winston_cone_inner_radius2 + LGC_winston_cone_thickness );
+ //UnionSolid winston_cone_solid(winston_tube,winston_cone,Position(0,0,LGC_winston_tube_length / 2.0 - LGC_winston_cone_inset_length));
+ Paraboloid winston_cone1(LGC_winston_cone_inner_radius1 + LGC_winston_cone_thickness,
+ LGC_winston_cone_inner_radius2 + LGC_winston_cone_thickness,
+ LGC_winston_cone_length / 2.0 );
+ Paraboloid winston_cone2(LGC_winston_cone_inner_radius1,
+ LGC_winston_cone_inner_radius2,
+ LGC_winston_cone_length / 2.0 );
+
+ SubtractionSolid winston_cone(winston_cone1, winston_cone2);
+
+ Volume v_winston_cone_solid("v_winston_cone_solid", winston_cone, PyrexGlass);
+ PlacedVolume pv_winston_cone_solid = v_sector.placeVolume(
+ v_winston_cone_solid, Transform3D(Position(0, LGC_pmt_y_pos, LGC_pmt_z_pos)) *
+ RotationX(LGC_pmt_tilt_angle) *
+ Transform3D(Position(0, 0, LGC_winston_tube_length / 2.0 + 5.0 * mm)));
+ de_winston_cone.setPlacement(pv_winston_cone_solid);
+
+ //std::cout << " LGC_pmt_y_pos/cm " << LGC_pmt_y_pos/cm << "\n";
+ //std::cout << " LGC_pmt_z_pos/cm " << LGC_pmt_z_pos/cm << "\n";
+ //mirrorPV.addPhysVolID("layer", 2).addPhysVolID("module", 1);
+ //mirror_DE.setPlacement(mirrorPV);
+ //sens.setType("photoncounter");
+ //mirrorVol.setSensitiveDetector(sens);
+
+ //// ---------------
+
+ //// ---------------
+ //// Dummy PMT surface
+ DetElement de_pmt_array(det, "PMT_DE", 1);
+ Box pmt_array(LGC_pmt_array_size / 2.0, LGC_pmt_array_size / 2.0, 5 * mm / 2.0);
+ Volume v_pmt_array("v_pmt_array", pmt_array, N2cherenkov);
+ PlacedVolume pv_pmt_array =
+ v_sector.placeVolume(v_pmt_array, Transform3D(Position(0, LGC_pmt_y_pos, LGC_pmt_z_pos)) *
+ RotationX(LGC_pmt_tilt_angle));
+
+ pv_pmt_array.addPhysVolID("mirror", 3);
+ de_pmt_array.setPlacement(pv_pmt_array);
+ sens.setType("photoncounter");
+ v_pmt_array.setSensitiveDetector(sens);
+
+ // copper layer inside to stop photons
+ Box pmt_array_backing(LGC_pmt_array_size/2.0, LGC_pmt_array_size/2.0, 1*mm/2.0);
+ Volume v_pmt_array_backing("v_pmt_array_backing", pmt_array_backing, Copper);
+ PlacedVolume pv_pmt_array_backing = v_pmt_array.placeVolume(v_pmt_array_backing, Position(0,0,0));
+
+ // Optical Surfaces
+
+ OpticalSurfaceManager surfMgr = desc.surfaceManager();
+ OpticalSurface mirrorSurf = surfMgr.opticalSurface("MirrorOpticalSurface");
+ OpticalSurface pmtSurf = surfMgr.opticalSurface("PMTOpticalSurface");
+ //BorderSurface mirrorBorder_Surf = BorderSurface(desc, det, "RICHmirror", mirrorSurf, mirrorPV, envPV);
+ SkinSurface mirrorBorder_Surf(desc,de_mirror1_shape,"RICHmirror", mirrorSurf, v_mirror1_shape);
+ SkinSurface winstonBorder_Surf(desc,de_winston_cone,"WinstonCone", mirrorSurf, v_winston_cone_solid);
+ SkinSurface pmtBorder_Surf(desc,de_pmt_array,"PMTsurface", pmtSurf, v_pmt_array);
+ //BorderSurface bubbleSurf = BorderSurface(description, sdet, "TankBubble", airSurf, bubblePlace, tankPlace);
+ mirrorBorder_Surf.isValid();
+ winstonBorder_Surf.isValid();
+ pmtBorder_Surf.isValid();
+ //tankSurf.isValid();
+
+
+ // all sectors
+ for (int i_sector = 1; i_sector <= 30; i_sector++) {
+ //std::cout << i_sector << " sector\n";
+ PlacedVolume pv =
+ v_lgc_tank.placeVolume(v_sector, Transform3D(RotationZ((i_sector - 1) * LGC_sector_angle)));
+ pv.addPhysVolID("sector", i_sector);
+ auto amod = (i_sector == 1 ? de_sector : de_sector.clone("de_sector" + std::to_string(i_sector), i_sector));
+ amod.setPlacement(pv);
+ det.add(amod);
+ }
+
+
+ //// ---------------
+
+
+ Volume motherVol = desc.pickMotherVolume(det);
+ PlacedVolume envPV = motherVol.placeVolume(v_lgc_tank, Position(pos_x, pos_y, pos_z));
+ envPV.addPhysVolID("system", detID);
+ det.setPlacement(envPV);
+
+ return det;
+}
+//@}
+// clang-format off
+DECLARE_DETELEMENT(SoLID_HGC, createDetector)
+
diff --git a/src/ThresholdGasCherenkov_geo.cpp b/src/ThresholdGasCherenkov_geo.cpp
new file mode 100644
index 0000000..7ddf202
--- /dev/null
+++ b/src/ThresholdGasCherenkov_geo.cpp
@@ -0,0 +1,259 @@
+#include "DD4hep/DetFactoryHelper.h"
+#include "DD4hep/Printout.h"
+#include <XML/Helper.h>
+#include "TMath.h"
+#include "DDRec/Surface.h"
+
+#include "DD4hep/OpticalSurfaces.h"
+#include "DDRec/DetectorData.h"
+
+using namespace std;
+using namespace dd4hep;
+using namespace dd4hep::rec;
+
+using namespace dd4hep;
+
+/** \addtogroup PID Particle ID Detectors
+ */
+/** \addtogroup ThresholdGasCherenkov Light Gas (threshold) Cherenkov detector.
+ * \brief Type: **ThresholdGasCherenkov**.
+ * \ingroup PID
+ *
+ * \code
+ * <detector>
+ * </detector>
+ * \endcode
+ *
+ * @{
+ */
+static Ref_t createDetector(Detector& desc, xml::Handle_t handle, SensitiveDetector sens) {
+ xml::DetElement detElem = handle;
+ xml_det_t x_det = handle;
+
+ std::string detName = detElem.nameStr();
+ int detID = detElem.id();
+
+ DetElement det(detName, detID);
+ xml::Component dims = detElem.dimensions();
+ double rInner = dims.rmin();
+ double rOuter1 = dims.rmax1();
+ double rOuter2 = dims.rmax2();
+ double zMin = dims.zmin();
+ double zMax = dims.zmax();
+
+ std::map<int,Position> mirror_positions;
+ std::map<int,std::array<double,3>> mirror_rotations;
+
+ for(xml_coll_t i(x_det,Unicode("mirror")); i; ++i){
+ xml_comp_t x_mir = i;
+ //std::cout << "mirror " << x_mir.id() << "\n";
+ xml_dim_t mir_pos = x_mir.child(_U(placement));
+ xml_dim_t mir_rot = x_mir.child(_U(rotation));
+ mirror_positions[x_mir.id()] = Position(mir_pos.x(), mir_pos.y(), mir_pos.z());
+ mirror_rotations[x_mir.id()] = {mir_rot.x(), mir_rot.y(), mir_rot.z()};
+ }
+
+ xml_dim_t pos = x_det.child(_U(placement));
+ double pos_x = pos.x();
+ double pos_y = pos.y();
+ double pos_z = pos.z();
+
+ Material air = desc.air();
+ Material PyrexGlass = desc.material("PyrexGlass");
+ Material N2cherenkov = desc.material("N2cherenkov");
+ Material Copper = desc.material("Copper");
+
+ double LGC_inner_radius1 = 71.0*cm;
+ double LGC_inner_radius2 = 85.0*cm;
+ double LGC_outer_radius1 = 265.0*cm;
+ double LGC_main_length = 105.0*cm;
+ double LGC_snout_length = 107.0*cm;
+ double LGC_snout_inner_radius1 = 58.0*cm;
+ double LGC_snout_inner_radius2 = LGC_inner_radius1;
+ double LGC_snout_outer_radius1 = 127.0*cm;
+ double LGC_snout_outer_radius2 = 144.0*cm;
+ double LGC_entrance_window_thickness = 0.05*mm; // something tells this might be 5 mil, not mm
+ double LGC_exit_window_thickness = 0.1*mm; // same here
+ double LGC_mirror1_radius = 277.51*cm;
+ double LGC_mirror2_radius = 157.99*cm;
+ double LGC_mirror1_length = 114.53*cm;
+ double LGC_mirror2_length = 59.260*cm;
+ double LGC_mirror1_width1 = 16.26*cm;
+ double LGC_mirror1_width2 = 36.03*cm;
+ double LGC_mirror2_width1 = 37.06*cm;
+ double LGC_mirror2_width2 = 45.95*cm;
+ double LGC_mirror1_thickness = 2.0*mm;
+ double LGC_mirror2_thickness = 2.0*mm;
+
+ double LGC_sector_angle = M_PI * 15.0 / 180.0;
+ //double LGC_scattering_angle = 11.0 * M_PI / 180.0;
+
+ double LGC_mirror1_tilt_angle = mirror_rotations[1][0];//25.0 * M_PI / 180.0;
+ double LGC_mirror2_tilt_angle = mirror_rotations[2][0];//2.0 * M_PI / 180.0;
+ double LGC_pmt_tilt_angle = mirror_rotations[3][0];//45.0 * M_PI / 180.0;
+
+ double LGC_pmt_z_pos = mirror_positions[3].z();//-30.0*cm;
+ double LGC_pmt_y_pos = mirror_positions[3].y();//LGC_outer_radius1 - 20.0*cm;
+ double LGC_pmt_array_size = 20.0*cm;
+
+ // the gas tank
+ ConeSegment tank_main(0.5 * LGC_main_length, LGC_inner_radius1, LGC_outer_radius1,
+ LGC_inner_radius2, LGC_outer_radius1);
+ // M_PI / 2.0 - LGC_sector_angle / 2.0,
+ //M_PI / 2.0 + LGC_sector_angle / 2.0);
+ ConeSegment tank_snout(0.5 * LGC_snout_length, LGC_snout_inner_radius1, LGC_snout_outer_radius1,
+ LGC_snout_inner_radius2, LGC_snout_outer_radius2);
+ //M_PI / 2.0 - LGC_sector_angle / 2.0,
+ //M_PI / 2.0 + LGC_sector_angle / 2.0);
+ UnionSolid sidis_tank(tank_main,tank_snout,Position(0,0,-0.5 * LGC_main_length -0.5 * LGC_snout_length));
+ Volume v_lgc_tank("v_lgc_tank_gas", sidis_tank, N2cherenkov);
+ v_lgc_tank.setVisAttributes(desc,dd4hep::getAttrOrDefault(x_det, _Unicode(vis), "BlueVis"));
+
+ // Everything that goes in the tank will be copies of the sector assembly volume
+ Assembly v_sector("cherenkov_sector_1");
+ DetElement de_sector("de_sector"+std::to_string(1),1);
+
+ // mirrors
+ Sphere mirror1_shell(LGC_mirror1_radius, LGC_mirror1_radius + LGC_mirror1_thickness,
+ 0.0, M_PI / 2);
+ Trd1 mirror1_cutout(LGC_mirror1_width1 / 2.0, LGC_mirror1_width2 / 2.0,
+ LGC_mirror1_length / 2.0, LGC_mirror1_length / 2.0);
+ IntersectionSolid mirror1_shape(mirror1_cutout, mirror1_shell,
+ RotationX(M_PI/2.0)*Transform3D(Position(0, 0, -LGC_mirror1_radius)));
+ Sphere mirror2_shell(LGC_mirror2_radius, LGC_mirror2_radius + LGC_mirror2_thickness,
+ 0.0, M_PI / 2);
+ Trd1 mirror2_cutout(LGC_mirror2_width1 / 2.0, LGC_mirror2_width2 / 2.0,
+ LGC_mirror2_length / 2.0, LGC_mirror2_length / 2.0);
+ IntersectionSolid mirror2_shape(mirror2_cutout, mirror2_shell,
+ RotationX(M_PI/2.0)*Transform3D(Position(0, 0, -LGC_mirror2_radius)));
+
+ double z_mirror1 = mirror_positions[1].z();
+ double z_mirror2 = mirror_positions[2].z();
+ double y_mirror1 = mirror_positions[1].y();
+ double y_mirror2 = mirror_positions[2].y();
+
+ Volume v_mirror1_shape("v_mirror1_shape", mirror1_shape, PyrexGlass);
+ PlacedVolume pv_mirror1_shape = v_sector.placeVolume(
+ v_mirror1_shape, Transform3D(Position(0, y_mirror1, z_mirror1)) *
+ RotationX(-M_PI / 2.0 + LGC_mirror1_tilt_angle));
+
+ DetElement de_mirror1_shape(det,"de_mirror1_shape"+std::to_string(1),1);
+ pv_mirror1_shape.addPhysVolID("mirror", 1);
+ de_mirror1_shape.setPlacement(pv_mirror1_shape);
+ sens.setType("photoncounter");
+ v_mirror1_shape.setSensitiveDetector(sens);
+
+ Volume v_mirror2_shape("v_mirror2_shape", mirror2_shape, PyrexGlass);
+ PlacedVolume pv_mirror2_shape = v_sector.placeVolume(
+ v_mirror2_shape, Transform3D(Position(0, y_mirror2, z_mirror2)) *
+ RotationX(-M_PI / 2.0 + LGC_mirror2_tilt_angle));
+
+ DetElement de_mirror2_shape(det,"de_mirror2_shape"+std::to_string(2),2);
+ pv_mirror2_shape.addPhysVolID("mirror", 2);
+ de_mirror2_shape.setPlacement(pv_mirror2_shape);
+ sens.setType("photoncounter");
+ v_mirror2_shape.setSensitiveDetector(sens);
+
+ // ---------------
+ // Winston Cone
+ double LGC_winston_cone_thickness = 4*mm;
+ double LGC_winston_tube_inner_radius = 11.28*cm;
+ double LGC_winston_tube_length = 30.0*cm;
+ double LGC_winston_cone_length = 30.0*cm;
+ double LGC_winston_cone_inner_radius1 = 7.8*cm;
+ double LGC_winston_cone_inner_radius2 = 21.0*cm;
+ double LGC_winston_cone_inset_length = 7.90909*cm;
+ DetElement de_winston_cone(det,"de_winston_cone1",1);
+ Tube winston_tube(LGC_winston_tube_inner_radius,
+ LGC_winston_tube_inner_radius + LGC_winston_cone_thickness,
+ LGC_winston_tube_length / 2.0);
+ //Cone winston_cone(LGC_winston_cone_length / 2.0, LGC_winston_cone_inner_radius1,
+ // LGC_winston_cone_inner_radius1 + LGC_winston_cone_thickness,
+ // LGC_winston_cone_inner_radius2,
+ // LGC_winston_cone_inner_radius2 + LGC_winston_cone_thickness );
+ //UnionSolid winston_cone_solid(winston_tube,winston_cone,Position(0,0,LGC_winston_tube_length / 2.0 - LGC_winston_cone_inset_length));
+ Paraboloid winston_cone1(LGC_winston_cone_inner_radius1 + LGC_winston_cone_thickness,
+ LGC_winston_cone_inner_radius2 + LGC_winston_cone_thickness,
+ LGC_winston_cone_length / 2.0 );
+ Paraboloid winston_cone2(LGC_winston_cone_inner_radius1,
+ LGC_winston_cone_inner_radius2,
+ LGC_winston_cone_length / 2.0 );
+
+ SubtractionSolid winston_cone(winston_cone1, winston_cone2);
+
+ Volume v_winston_cone_solid("v_winston_cone_solid", winston_cone, PyrexGlass);
+ PlacedVolume pv_winston_cone_solid = v_sector.placeVolume(
+ v_winston_cone_solid, Transform3D(Position(0, LGC_pmt_y_pos, LGC_pmt_z_pos)) *
+ RotationX(LGC_pmt_tilt_angle) *
+ Transform3D(Position(0, 0, LGC_winston_tube_length / 2.0 + 5.0 * mm)));
+
+ //std::cout << " LGC_pmt_y_pos/cm " << LGC_pmt_y_pos/cm << "\n";
+ //std::cout << " LGC_pmt_z_pos/cm " << LGC_pmt_z_pos/cm << "\n";
+ //mirrorPV.addPhysVolID("layer", 2).addPhysVolID("module", 1);
+ //mirror_DE.setPlacement(mirrorPV);
+ //sens.setType("photoncounter");
+ //mirrorVol.setSensitiveDetector(sens);
+
+ //// ---------------
+
+ //// ---------------
+ //// Dummy PMT surface
+ DetElement de_pmt_array(det, "PMT_DE", 1);
+ Box pmt_array(LGC_pmt_array_size / 2.0, LGC_pmt_array_size / 2.0, 5 * mm / 2.0);
+ Volume v_pmt_array("v_pmt_array", pmt_array, N2cherenkov);
+ PlacedVolume pv_pmt_array =
+ v_sector.placeVolume(v_pmt_array, Transform3D(Position(0, LGC_pmt_y_pos, LGC_pmt_z_pos)) *
+ RotationX(LGC_pmt_tilt_angle));
+
+ pv_pmt_array.addPhysVolID("mirror", 3);
+ de_pmt_array.setPlacement(pv_pmt_array);
+ sens.setType("photoncounter");
+ v_pmt_array.setSensitiveDetector(sens);
+
+ // copper layer inside to stop photons
+ Box pmt_array_backing(LGC_pmt_array_size/2.0, LGC_pmt_array_size/2.0, 1*mm/2.0);
+ Volume v_pmt_array_backing("v_pmt_array_backing", pmt_array_backing, Copper);
+ PlacedVolume pv_pmt_array_backing = v_pmt_array.placeVolume(v_pmt_array_backing, Position(0,0,0));
+
+ // Optical Surfaces
+
+ OpticalSurfaceManager surfMgr = desc.surfaceManager();
+ OpticalSurface mirrorSurf = surfMgr.opticalSurface("MirrorOpticalSurface");
+ OpticalSurface pmtSurf = surfMgr.opticalSurface("PMTOpticalSurface");
+ //BorderSurface mirrorBorder_Surf = BorderSurface(desc, det, "RICHmirror", mirrorSurf, mirrorPV, envPV);
+ SkinSurface mirrorBorder_Surf(desc,de_mirror1_shape,"LGCmirror", mirrorSurf, v_mirror1_shape);
+ SkinSurface winstonBorder_Surf(desc,de_winston_cone,"LGCWinstonCone", mirrorSurf, v_winston_cone_solid);
+ SkinSurface pmtBorder_Surf(desc,de_pmt_array,"LGCPMTsurface", pmtSurf, v_pmt_array);
+ //BorderSurface bubbleSurf = BorderSurface(description, sdet, "TankBubble", airSurf, bubblePlace, tankPlace);
+ mirrorBorder_Surf.isValid();
+ winstonBorder_Surf.isValid();
+ pmtBorder_Surf.isValid();
+ //tankSurf.isValid();
+
+
+ // all sectors
+ for (int i_sector = 1; i_sector <= 30; i_sector++) {
+ //std::cout << i_sector << " sector\n";
+ PlacedVolume pv =
+ v_lgc_tank.placeVolume(v_sector, Transform3D(RotationZ((i_sector - 1) * LGC_sector_angle)));
+ pv.addPhysVolID("sector", i_sector);
+ auto amod = (i_sector == 1 ? de_sector : de_sector.clone("de_sector" + std::to_string(i_sector), i_sector));
+ amod.setPlacement(pv);
+ det.add(amod);
+ }
+
+
+ //// ---------------
+
+
+ Volume motherVol = desc.pickMotherVolume(det);
+ PlacedVolume envPV = motherVol.placeVolume(v_lgc_tank, Position(pos_x, pos_y, pos_z));
+ envPV.addPhysVolID("system", detID);
+ det.setPlacement(envPV);
+
+ return det;
+}
+//@}
+// clang-format off
+DECLARE_DETELEMENT(SoLID_LGC, createDetector)
+
diff --git a/views/top_view/.DAWN_1.history b/views/top_view/.DAWN_1.history
index f4980ff..35e45a6 100644
--- a/views/top_view/.DAWN_1.history
+++ b/views/top_view/.DAWN_1.history
@@ -1,11 +1,11 @@
1.34392e+07
90
-180
+90
0
0
0
0
-1
+10
1
0.001
0
@@ -15,8 +15,8 @@
0.5
0.5
0.5
-19
-71
+90
+180
0.001
0.001
0.001
--
GitLab