Skip to content
Snippets Groups Projects
Commit f642bc8d authored by Chao Peng's avatar Chao Peng
Browse files

move gas cherenkov plugins to this repo

parent 4afba872
No related branches found
No related tags found
No related merge requests found
Pipeline #47371 passed with warnings
Stage: build
Stage: test
Hide whitespace changes
Inline Side-by-side
CMakeLists.txt
+ 2
− 0
View file @ f642bc8d
......@@ -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}
......
compact/HGC.xml
+ 1
− 1
View file @ f642bc8d
<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">
......
compact/LGC.xml
+ 1
− 1
View file @ f642bc8d
<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">
......
src/HeavyGasCherenkov_geo.cpp 0 → 100644
+ 257
− 0
View file @ f642bc8d
#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)
src/ThresholdGasCherenkov_geo.cpp 0 → 100644
+ 259
− 0
View file @ f642bc8d
#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)
views/top_view/.DAWN_1.history
+ 4
− 4
View file @ f642bc8d
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
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment