Select Git revision
HeavyGasCherenkov_geo.cpp
Chao Peng authored
HeavyGasCherenkov_geo.cpp 11.01 KiB
#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)