From 5f0879c6e4aff60ecc8a87c4795de23d9cc6df38 Mon Sep 17 00:00:00 2001
From: myliu <201916234@mail.sdu.edu.cn>
Date: Fri, 15 Mar 2024 13:12:36 +0800
Subject: [PATCH] Add some description
---
.../src/driftchamber/DriftChamber.cpp | 9 +
.../src/driftchamber/DriftChamber_stero.cpp | 307 ------------------
Digitisers/DCHDigi/src/DCHDigiAlg.h | 3 +
3 files changed, 12 insertions(+), 307 deletions(-)
delete mode 100644 Detector/DetDriftChamber/src/driftchamber/DriftChamber_stero.cpp
diff --git a/Detector/DetDriftChamber/src/driftchamber/DriftChamber.cpp b/Detector/DetDriftChamber/src/driftchamber/DriftChamber.cpp
index ab80570f..1ffdd929 100644
--- a/Detector/DetDriftChamber/src/driftchamber/DriftChamber.cpp
+++ b/Detector/DetDriftChamber/src/driftchamber/DriftChamber.cpp
@@ -1,3 +1,12 @@
+//====================================================================
+// Detector description implementation of the Drift Chamber
+//--------------------------------------------------------------------
+//
+// Author: Tao Lin
+// Mengyao Liu
+//
+//====================================================================
+
#include "DD4hep/DetFactoryHelper.h"
#include "DD4hep/Printout.h"
#include "XML/Layering.h"
diff --git a/Detector/DetDriftChamber/src/driftchamber/DriftChamber_stero.cpp b/Detector/DetDriftChamber/src/driftchamber/DriftChamber_stero.cpp
deleted file mode 100644
index 0b9061a9..00000000
--- a/Detector/DetDriftChamber/src/driftchamber/DriftChamber_stero.cpp
+++ /dev/null
@@ -1,307 +0,0 @@
-//====================================================================
-// Detector description implementation of the Drift Chamber
-//--------------------------------------------------------------------
-//
-// Author: Tao Lin
-// Mengyao Liu
-//
-//====================================================================
-
-#include "DD4hep/DetFactoryHelper.h"
-#include "DD4hep/Printout.h"
-#include "XML/Layering.h"
-#include "XML/Utilities.h"
-#include "XML/XMLElements.h"
-#include "DDRec/DetectorData.h"
-#include "DDSegmentation/Segmentation.h"
-#include "DetSegmentation/GridDriftChamber.h"
-
-using namespace dd4hep;
-using namespace dd4hep::detail;
-using namespace dd4hep::rec ;
-
-#define MYDEBUG(x) std::cout << __FILE__ << ":" << __LINE__ << ": " << x << std::endl;
-#define MYDEBUGVAL(x) std::cout << __FILE__ << ":" << __LINE__ << ": " << #x << ": " << x << std::endl;
-
-static dd4hep::Ref_t create_detector(dd4hep::Detector& theDetector,
- xml_h e,
- dd4hep::SensitiveDetector sens) {
- // ------- Lambda functions ---- //
- auto delta_a_func = [](auto x, auto y) { return 0.5 * ( x + y ); };
- auto delta_b_func = [](auto x, auto y) { return 2 * std::sqrt((x + y) * (x - y)); };
- auto epsilon_func = [](auto delta_a, auto L) { return std::atan(delta_a / L); };
-
- // =======================================================================
- // Parameter Definition
- // =======================================================================
-
- xml_det_t x_det = e;
-
- xml_coll_t c(x_det,_U(chamber));
- xml_comp_t x_chamber = c;
-
- xml_coll_t cc(x_det,_U(side));
- xml_comp_t x_side = cc;
-
- std::string det_name = x_det.nameStr();
- std::string det_type = x_det.typeStr();
-
- dd4hep::SensitiveDetector sd = sens;
-
- // - global
- double chamber_half_length = theDetector.constant<double>("DC_half_length");
- double chamber_length = theDetector.constant<double>("DC_length");
-
- // - chamber
- double chamber_radius_min = theDetector.constant<double>("SDT_chamber_radius_min");
- double chamber_radius_max = theDetector.constant<double>("SDT_chamber_radius_max");
- double SDT_half_length = theDetector.constant<double>("SDT_chamber_half_length");
- int chamberID = x_chamber.id();
-
- // - layer
- double chamber_layer_width = theDetector.constant<double>("SDT_chamber_layer_width");
- double chamber_cell_width = theDetector.constant<double>("SDT_chamber_cell_width");
- double chamber_layer_rbegin = theDetector.constant<double>("DC_chamber_layer_rbegin");
- double chamber_layer_rend = theDetector.constant<double>("DC_chamber_layer_rend");
- int chamber_layer_number = floor((chamber_layer_rend-chamber_layer_rbegin)/chamber_layer_width);
-
- double cellsize = theDetector.constant<double>("SDT_chamber_cell_width");
-
- double safe_distance = theDetector.constant<double>("DC_safe_distance");
- double alpha = theDetector.constant<double>("Alpha");
-
- // =======================================================================
- // Detector Construction
- // =======================================================================
-
- dd4hep::DetElement sdet(det_name, x_det.id());
-
- Volume envelope = dd4hep::xml::createPlacedEnvelope( theDetector, e , sdet ) ;
-
- dd4hep::xml::setDetectorTypeFlag( e, sdet ) ;
-
- if( theDetector.buildType() == BUILD_ENVELOPE ) return sdet ;
-
-
-// dd4hep::Material det_mat(theDetector.material("Air"));
- dd4hep::Material det_mat(theDetector.material(x_det.materialStr()));
-// dd4hep::Material chamber_mat(theDetector.material("GasHe_90Isob_10"));
- dd4hep::Material chamber_mat = theDetector.material(x_chamber.materialStr());
-
- // - global
- Assembly det_vol( det_name+"_assembly" ) ;
-
- // - chamber volume
- dd4hep::Tube det_chamber_solid(chamber_radius_min, chamber_radius_max, chamber_half_length);
- dd4hep::Volume det_chamber_vol(det_name+"_chamber_vol", det_chamber_solid, chamber_mat);
-
- // - wall
- double chamber_inner_wall_rmin = theDetector.constant<double>("SDT_chamber_inner_wall_radius_min");
- double chamber_inner_wall_rmax = theDetector.constant<double>("SDT_chamber_inner_wall_radius_max");
- double chamber_outer_wall_rmin = theDetector.constant<double>("SDT_chamber_outer_wall_radius_min");
- double chamber_outer_wall_rmax = theDetector.constant<double>("SDT_chamber_outer_wall_radius_max");
-
-// dd4hep::Material wall_mat(theDetector.material("CarbonFiber"));
- dd4hep::Material wall_mat(theDetector.material(x_side.materialStr()));
-
- double wall_rmin[2] = {chamber_inner_wall_rmin, chamber_outer_wall_rmin};
- double wall_rmax[2] = {chamber_inner_wall_rmax, chamber_outer_wall_rmax};
-
- // - wire
- dd4hep::Volume module_vol;
- dd4hep::Volume Module_vol;
- for(xml_coll_t c(x_det,_U(module)); c; ++c) {
- xml_comp_t x_module = c;
- double module_rmin = x_module.rmin();
- double module_rmax = x_module.rmax();
- std::string module_name = x_module.nameStr();
- dd4hep::Tube module_solid(module_rmin,module_rmax,chamber_half_length);
- if(x_module.id()==0) {
- module_vol = dd4hep::Volume(module_name,module_solid,det_mat);
- module_vol.setVisAttributes(theDetector.visAttributes(x_module.visStr()));
- } else {
- Module_vol = dd4hep::Volume(module_name,module_solid,det_mat);
- Module_vol.setVisAttributes(theDetector.visAttributes(x_module.visStr()));
- }
-
- for(xml_coll_t l(x_module,_U(tubs)); l; ++l) {
- xml_comp_t x_tube =l;
- double tube_rmin = x_tube.rmin();
- double tube_rmax = x_tube.rmax();
- std::string tube_name = x_tube.nameStr();
- std::string wire_name= module_name + tube_name;
- dd4hep::Material tube_mat = theDetector.material(x_tube.materialStr());
- dd4hep::Tube wire_solid(tube_rmin,tube_rmax,chamber_half_length);
- dd4hep::Volume wire_vol(wire_name,wire_solid,tube_mat);
- dd4hep::Transform3D transform_wire(dd4hep::Rotation3D(),dd4hep::Position(0.,0.,0.));
- dd4hep::PlacedVolume wire_phy;
- if(x_module.id()==0) {
- wire_phy = module_vol.placeVolume(wire_vol,transform_wire);
- } else {
- wire_phy = Module_vol.placeVolume(wire_vol,transform_wire);
- }
- }
- }
-
- // End cap
- double Endcap_rmin = theDetector.constant<double>("DC_Endcap_rmin");
- double Endcap_rmax = theDetector.constant<double>("DC_Endcap_rmax");
- double Endcap_z = theDetector.constant<double>("DC_Endcap_dz");
- dd4hep::Tube det_Endcap_solid(Endcap_rmin,Endcap_rmax,0.5*Endcap_z);
- dd4hep::Volume det_Endcap_vol(det_name+"Endcap",det_Endcap_solid,det_mat);
- det_Endcap_vol.setVisAttributes(theDetector,"YellowVis");
-
- //Initialize the segmentation
- dd4hep::Readout readout = sd.readout();
- dd4hep::Segmentation geomseg = readout.segmentation();
- dd4hep::Segmentation* _geoSeg = &geomseg;
-
- auto DCHseg = dynamic_cast<dd4hep::DDSegmentation::GridDriftChamber*>(_geoSeg->segmentation());
-
- // - layer
- int chamber_id = 0;
- int layerIndex = -1;
- for(int layer_id = 0; layer_id < chamber_layer_number; layer_id++) {
-
- double rmin,rmax,offset=0;
-
- double sign_eps = 1;
-
- dd4hep::Volume* current_vol_ptr = nullptr;
-// current_vol_ptr = &det_chamber_vol;
- rmin = chamber_layer_rbegin+(layer_id*chamber_layer_width);
- rmax = rmin+chamber_layer_width;
- layerIndex = layer_id;
-
- //Construction of drift chamber layers
- double rmid = delta_a_func(rmin,rmax);
- double Rmid = rmid/std::cos(alpha/2);
- double ilayer_cir = 2 * M_PI * rmid;
- double ncell = ilayer_cir / chamber_layer_width;
- int ncell_layer = floor(ncell);
- int numWire = ncell_layer;
- double layer_Phi = 2*M_PI / ncell_layer;
-
- if(layer_id %2 ==0)
- {
- offset = 0.;
- sign_eps = -1;
- }
- else { offset = 0.5 * layer_Phi; }
-
-
- double epsilon = sign_eps*std::atan(2 * Rmid * std::tan(alpha / 2.0) / chamber_length);
- double alpha0 = 2*std::asin(chamber_length * std::tan(epsilon)/(2*Rmid));
-
- DCHseg->setGeomParams(chamber_id, layerIndex, layer_Phi, rmid, epsilon, offset);
- DCHseg->setWiresInLayer(chamber_id, layerIndex, numWire);
-
- double r_in_test = rmid*std::cos(alpha / 2.0);
-
- double r_in0 = rmid - cellsize / 2.0;
- double r_in = r_in0 / std::cos(alpha / 2.0);
-
- double r_out0 = rmid + cellsize / 2.0;
- double r_out = r_out0 / std::cos(alpha / 2.0);
-
- double delta_a_in = delta_b_func(r_in, r_in0);
- double delta_a_out = delta_b_func(r_out, r_out0);
-
- double eps_in = epsilon_func(delta_a_in, chamber_length );
- double eps_out = epsilon_func(delta_a_out, chamber_length );
-
-// dd4hep::Tube layer_vol_solid(rmin,rmax,chamber_half_length);
- dd4hep::Hyperboloid layer_vol_solid(r_in0, eps_in, r_out0, eps_out, chamber_half_length);
- dd4hep::Volume layer_vol(det_name+"_layer_vol",layer_vol_solid,det_mat);
- current_vol_ptr = &layer_vol;
-
- if ( x_det.isSensitive() ) {
- layer_vol.setRegion(theDetector,x_det.regionStr());
- layer_vol.setLimitSet(theDetector,x_det.limitsStr());
- layer_vol.setSensitiveDetector(sens);
- sd.setType("tracker");
- }
-
- // - wire vol
- //phi <-------------------> -phi
- // | F8 F7 F6 F5| Only on the outermost layer.
- // | |
- // | S F4|
- // | |
- // | F0 F1 F2 F3|
- // -----------------------
- for(int icell=0; icell< numWire; icell++) {
- double wire_phi = (icell+0.5)*layer_Phi + offset;
-
- // - signal wire
- dd4hep::RotationZ rz(wire_phi);
- dd4hep::RotationY ry(epsilon);
- dd4hep::Position tr3D = Position(rmid*std::cos(0.5*M_PI+wire_phi),rmid*std::sin(0.5*M_PI+wire_phi),0.);
- dd4hep::Transform3D transform_signal_wire(rz*ry,tr3D);
-
- dd4hep::PlacedVolume module_phy = (*current_vol_ptr).placeVolume(module_vol,transform_signal_wire);
- // - Field wire
- dd4hep::PlacedVolume Module_phy;
- double radius[5] = {rmid-chamber_layer_width*0.5+safe_distance,rmid-chamber_layer_width*0.5+safe_distance,rmid,rmid+chamber_layer_width*0.5-safe_distance,rmid+chamber_layer_width*0.5-safe_distance};
- double phi[5] = {wire_phi,wire_phi-layer_Phi*0.5,wire_phi-layer_Phi*0.5,wire_phi-layer_Phi*0.5,wire_phi};
- int num = 3;
- if(layer_id==(chamber_layer_number-1)) {
- num = 5;
- }
- for(int i=0; i<num ; i++) {
- dd4hep::RotationZ rz_field(phi[i]);
- dd4hep::RotationY ry_field(epsilon);
- dd4hep::Position tr3D_field = Position(radius[i]*std::cos(0.5*M_PI+phi[i]),radius[i]*std::sin(0.5*M_PI+phi[i]),0.);
- dd4hep::Transform3D transform_field_wire(rz_field*ry_field,tr3D_field);
-
- Module_phy = (*current_vol_ptr).placeVolume(Module_vol,transform_field_wire);
- }
- }
- dd4hep::Transform3D transform_layer(dd4hep::Rotation3D(),
- dd4hep::Position(0,0,0));
- dd4hep::PlacedVolume layer_phy = det_chamber_vol.placeVolume(layer_vol,transform_layer);
- layer_phy.addPhysVolID("layer", layer_id);
- }
-
- // - place in det
- // - chamber
- dd4hep::Transform3D transform_chamber(dd4hep::Rotation3D(),
- dd4hep::Position(0,0,0));
- dd4hep::PlacedVolume det_chamber_phy = det_vol.placeVolume(det_chamber_vol,
- transform_chamber);
-
- det_chamber_phy.addPhysVolID("chamber", chamberID);
-
- // - place in world
- dd4hep::Transform3D transform(dd4hep::Rotation3D(),
- dd4hep::Position(0,0,0));
- dd4hep::PlacedVolume phv = envelope.placeVolume(det_vol,transform);
- // - place wall
- dd4hep::PlacedVolume wall_phy;
- for(int i=0; i<2; i++) {
- dd4hep::Tube wall_solid(wall_rmin[i],wall_rmax[i],chamber_half_length);
- dd4hep::Volume wall_vol(det_name+"_wall_vol",wall_solid,wall_mat);
- wall_vol.setVisAttributes(theDetector,"VisibleGreen");
- wall_phy = envelope.placeVolume(wall_vol,transform);
- }
-
- // - place Endcap
- double endcap_pos[2] = {chamber_half_length+Endcap_z*0.5,-chamber_half_length-Endcap_z*0.5};
- dd4hep::PlacedVolume endcap_phy;
- for(int i=0; i<2; i++) {
- dd4hep::Transform3D Endcap_transform(dd4hep::Rotation3D(),dd4hep::Position(0,0,endcap_pos[i]));
- endcap_phy = envelope.placeVolume(det_Endcap_vol,Endcap_transform);
- }
-
- if ( x_det.hasAttr(_U(id)) ) {
- phv.addPhysVolID("system",x_det.id());
- }
- DetElement assDE( sdet , det_name+"_assembly" , x_det.id() ) ;
- assDE.setPlacement(phv);
-
- MYDEBUG("Build Detector Drift Chamber successfully.");
- return sdet;
-
-}
-
-DECLARE_DETELEMENT(DriftChamber_Stero, create_detector);
diff --git a/Digitisers/DCHDigi/src/DCHDigiAlg.h b/Digitisers/DCHDigi/src/DCHDigiAlg.h
index cc6a5312..946363c1 100644
--- a/Digitisers/DCHDigi/src/DCHDigiAlg.h
+++ b/Digitisers/DCHDigi/src/DCHDigiAlg.h
@@ -54,10 +54,13 @@ protected:
NTuple::Array<int > m_chamber ;
NTuple::Array<int > m_layer ;
NTuple::Array<int > m_cell ;
+ //The position of wire at -z
NTuple::Array<float> m_cell_x ;
NTuple::Array<float> m_cell_y ;
+ //The position of wire at +z
NTuple::Array<float> m_cell1_x ;
NTuple::Array<float> m_cell1_y ;
+
NTuple::Array<float> m_simhit_x ;
NTuple::Array<float> m_simhit_y ;
NTuple::Array<float> m_simhit_z ;
--
GitLab