Skip to content
Snippets Groups Projects
Commit c315293d authored by Jihee Kim's avatar Jihee Kim
Browse files

Resolve "Implement cb_VTX"

parent 0fcb2b44
No related branches found
No related tags found
1 merge request!36Resolve "Implement cb_VTX"
Pipeline #8214 passed with stages
Stage: config
    Stage: build
      Stage: docs
        Stage: test
          Stage: finalize
            in 2 minutes and 44 seconds
            Hide whitespace changes
            Inline Side-by-side
            compact/cb_VTX_Barrel.xml 0 → 100644
            + 39
            − 0
            View file @ c315293d
            <lccdd>
            <comment> Central Barrel Vertex Detector </comment>
            <define>
            </define>
            <limits>
            </limits>
            <regions>
            </regions>
            <display>
            </display>
            <detectors>
            <detector id="cb_VTX_Barrel_ID" name="cb_VTX_Barrel" type="cb_VTX_Barrel" readout="cb_VTX_Barrel_Hits" vis="cb_VTX_BarrelVis" insideTrackingVloume="true">
            <dimensions rmin="cb_VTX_Barrel_rmin" rmax="cb_VTX_Barrel_rmax" length="cb_VTX_Barrel_length"/>
            <position x="cb_VTX_Barrel_x_pos" y="cb_VTX_Barrel_y_pos" z="cb_VTX_Barrel_z_pos"/>
            <layer repeat="cb_VTX_Barrel_layer" vis="cb_VTX_Barrel_layerVis">
            <slice name="Silicon_slice" material="Si" sensitive="true"/>
            </layer>
            </detector>
            </detectors>
            <readouts>
            <readout name="cb_VTX_Barrel_Hits">
            <segmentation type="CartesianGridXY" grid_size_x="1.0*mm" grid_size_y="1.0*mm" />
            <id>system:8,barrel:3,layer:4,module:8,x:32:-16,y:-16</id>
            </readout>
            </readouts>
            <plugins>
            </plugins>
            <fields>
            </fields>
            </lccdd>
            compact/definitions.xml
            + 11
            − 2
            View file @ c315293d
            ...@@ -68,7 +68,7 @@ ...@@ -68,7 +68,7 @@
            Layer 4 subassembly ID: 29 Layer 4 subassembly ID: 29
            Layer 5 subassembly ID: 30 Layer 5 subassembly ID: 30
            </comment> </comment>
            <constant name="SiVertexSubAssembly_ID" value="25"/> <constant name="cb_VTX_Barrel_ID" value="25"/>
            <comment> <comment>
            Unused values reserved for vertex: Unused values reserved for vertex:
            constant name="SiVertexSubAssemblyLayer1_ID" value="121" constant name="SiVertexSubAssemblyLayer1_ID" value="121"
            ...@@ -337,7 +337,16 @@ ...@@ -337,7 +337,16 @@
            <constant name="VertexTrackerOuterRadius" value="90.0*mm"/> <constant name="VertexTrackerOuterRadius" value="90.0*mm"/>
            <constant name="VertexTrackerOuterBarrelLength" value="400.0*mm"/> <constant name="VertexTrackerOuterBarrelLength" value="400.0*mm"/>
            <constant name="cb_VTX_Barrel_rmin" value="33.0*mm"/>
            <constant name="cb_VTX_Barrel_rmax" value="165.0*mm"/>
            <constant name="cb_VTX_Barrel_length" value="740.0*mm"/>
            <constant name="cb_VTX_Barrel_x_pos" value="0.0*mm"/>
            <constant name="cb_VTX_Barrel_y_pos" value="0.0*mm"/>
            <constant name="cb_VTX_Barrel_z_pos" value="0.0*mm"/>
            <constant name="cb_VTX_Barrel_layer" value="6"/>
            <comment> <comment>
            -------------------------- --------------------------
            Silicon Tracker Parameters Silicon Tracker Parameters
            ...@@ -370,7 +379,7 @@ ...@@ -370,7 +379,7 @@
            The cb_CTD_* parameters can probably be moved to the detector specific xml file. The cb_CTD_* parameters can probably be moved to the detector specific xml file.
            </comment> </comment>
            <constant name="cb_CTD_rmin" value="SiliconTrackerInnerRadius"/> <constant name="cb_CTD_rmin" value="SiliconTrackerInnerRadius + 70.0 * mm"/>
            <constant name="cb_CTD_rmax" value="SiliconTrackerOuterRadius"/> <constant name="cb_CTD_rmax" value="SiliconTrackerOuterRadius"/>
            <constant name="cb_CTD_zmax" value="SiliconTrackerOuterBarrel_zmax"/> <constant name="cb_CTD_zmax" value="SiliconTrackerOuterBarrel_zmax"/>
            <constant name="cb_CTD_length" value="SiliconTrackerOuterBarrelLength"/> <constant name="cb_CTD_length" value="SiliconTrackerOuterBarrelLength"/>
            ......
            compact/display.xml
            + 2
            − 0
            View file @ c315293d
            ...@@ -55,6 +55,8 @@ ...@@ -55,6 +55,8 @@
            <vis name="cb_CTD_Si_layerVis" alpha="1.0" r= "0.9" g="1.0" b="0.1" showDaughters="true" visible="true"/> <vis name="cb_CTD_Si_layerVis" alpha="1.0" r= "0.9" g="1.0" b="0.1" showDaughters="true" visible="true"/>
            <vis name="ce_GEMVis" alpha="0.1" r= "0.1" g="0.0" b="1.0" showDaughters="true" visible="true"/> <vis name="ce_GEMVis" alpha="0.1" r= "0.1" g="0.0" b="1.0" showDaughters="true" visible="true"/>
            <vis name="cb_GEM_layerVis" alpha="0.8" r= "0.8" g="0.4" b="0.3" showDaughters="true" visible="true"/> <vis name="cb_GEM_layerVis" alpha="0.8" r= "0.8" g="0.4" b="0.3" showDaughters="true" visible="true"/>
            <vis name="cb_VTX_BarrelVis" alpha="0.1" r= "0.1" g="0.0" b="1.0" showDaughters="true" visible="true"/>
            <vis name="cb_VTX_Barrel_layerVis" alpha="2.0" r= "0.0" g="0.2" b="0.8" showDaughters="true" visible="true"/>
            <vis name="ci_GEMVis" r= "0.8" g="0.4" b="0.3" alpha="0.8" showDaughters="true" visible="true"/> <vis name="ci_GEMVis" r= "0.8" g="0.4" b="0.3" alpha="0.8" showDaughters="true" visible="true"/>
            <vis name="ci_HCALVis" r= "0.6" g="0" b="0.6" alpha="1.0" showDaughters="true" visible="true"/> <vis name="ci_HCALVis" r= "0.6" g="0" b="0.6" alpha="1.0" showDaughters="true" visible="true"/>
            ......
            reference_detector.xml
            + 2
            − 1
            View file @ c315293d
            ...@@ -8,7 +8,7 @@ ...@@ -8,7 +8,7 @@
            url="https://eicweb.phy.anl.gov/EIC/detectors/reference_detector.git" url="https://eicweb.phy.anl.gov/EIC/detectors/reference_detector.git"
            status="development" status="development"
            version="v1 2021-03-16"> version="v1 2021-03-16">
            <comment>Test Detector</comment> <comment>Reference Detector</comment>
            </info> </info>
            <define> <define>
            ...@@ -125,6 +125,7 @@ ...@@ -125,6 +125,7 @@
            <include ref="compact/hcal.xml"/> <include ref="compact/hcal.xml"/>
            <include ref="compact/forward_trd.xml"/> <include ref="compact/forward_trd.xml"/>
            <include ref="compact/B0_tracker.xml"/> <include ref="compact/B0_tracker.xml"/>
            <include ref="compact/cb_VTX_Barrel.xml"/>
            <!-- <!--
            <include ref="compact/ci_HCAL.xml"/> <include ref="compact/ci_HCAL.xml"/>
            <include ref="compact/forward_rich.xml"/> <include ref="compact/forward_rich.xml"/>
            ......
            src/cb_VTX_Barrel.cpp 0 → 100644
            + 117
            − 0
            View file @ c315293d
            #include <XML/Helper.h>
            #include "DDRec/Surface.h"
            #include "DDRec/DetectorData.h"
            #include "DD4hep/OpticalSurfaces.h"
            #include "DD4hep/DetFactoryHelper.h"
            #include "DD4hep/Printout.h"
            //////////////////////////////////
            // Central Barrel Vertex Detector
            //////////////////////////////////
            using namespace std;
            using namespace dd4hep;
            static Ref_t createDetector(Detector& desc, xml_h e, SensitiveDetector sens)
            {
            xml_det_t x_det = e;
            string detName = x_det.nameStr();
            int detID = x_det.id();
            xml_dim_t dim = x_det.dimensions();
            double RIn = dim.rmin();
            double ROut = dim.rmax();
            double SizeZ = dim.length();
            xml_dim_t pos = x_det.position();
            Material Vacuum = desc.material("Vacuum");
            // Create Global Volume
            Tube cb_VTX_Barrel_GVol_Solid(RIn, ROut, SizeZ / 2.0, 0., 360.0 * deg);
            Volume detVol("cb_VTX_Barrel_GVol_Logic", cb_VTX_Barrel_GVol_Solid, Vacuum);
            detVol.setVisAttributes(desc.visAttributes(x_det.visStr()));
            //////////////////
            // Barrel Ladder
            //////////////////
            xml_comp_t x_layer = x_det.child(_U(layer));
            const int repeat = x_layer.repeat();
            xml_comp_t x_slice = x_layer.child(_U(slice));
            Material slice_mat = desc.material(x_slice.materialStr());
            double x = 0.0 * cm;
            double y = 0.0 * cm;
            double z = 0.0 * cm;
            int FDIV = 0;
            double dR = 0.0;
            double length = 0.0;
            double phi = 0.0;
            // Ladder Layer Parameters
            double lay_Dx[6];
            double lay_Dy[6];
            double lay_Dz[6];
            double lay_Rin[6];
            lay_Dx[0] = 0.050 * mm; lay_Dy[0] = 2.0 * cm; lay_Dz[0] = 10.0 * cm; lay_Rin[0] = 3.5 * cm;
            lay_Dx[1] = 0.050 * mm; lay_Dy[1] = 2.0 * cm; lay_Dz[1] = 11.0 * cm; lay_Rin[1] = 4.5 * cm;
            lay_Dx[2] = 0.150 * mm; lay_Dy[2] = 4.0 * cm; lay_Dz[2] = 18.0 * cm; lay_Rin[2] = 6.5 * cm;
            lay_Dx[3] = 0.150 * mm; lay_Dy[3] = 4.0 * cm; lay_Dz[3] = 24.0 * cm; lay_Rin[3] = 10.5 * cm;
            lay_Dx[4] = 0.150 * mm; lay_Dy[4] = 4.0 * cm; lay_Dz[4] = 36.0 * cm; lay_Rin[4] = 13.5 * cm;
            lay_Dx[5] = 0.150 * mm; lay_Dy[5] = 4.0 * cm; lay_Dz[5] = 48.0 * cm; lay_Rin[5] = 15.5 * cm;
            int i_layer = 0;
            int i_module = 0;
            // Loop over layers
            for(int i = 0; i < repeat; i++) {
            double cb_VTX_Barrel_ladder_DZ = lay_Dz[i];
            double cb_VTX_Barrel_ladder_DY = lay_Dy[i];
            double cb_VTX_Barrel_ladder_Thickness = lay_Dx[i];
            dR = lay_Rin[i];
            length = 2.0 * 3.1415 * dR;
            int laddersCount = length / cb_VTX_Barrel_ladder_DY;
            for (int i = 0; i < 2; i++) {
            double LN = cb_VTX_Barrel_ladder_DY * laddersCount;
            double LN1 = cb_VTX_Barrel_ladder_DY * (laddersCount + 1.0 + i);
            if (LN/LN1 > 0.8)
            laddersCount = laddersCount + 1;
            }
            double cb_VTX_Barrel_ladder_deltaphi = 2.0 * 3.1415926 / laddersCount;
            string ladderBoxName = detName + _toString(i, "_ladder_Solid_%d");
            string ladderName = detName + _toString(i, "_ladder_Logic_%d");
            Volume ladderVol(ladderName, Box(cb_VTX_Barrel_ladder_Thickness * 0.5, cb_VTX_Barrel_ladder_DY * 0.5, cb_VTX_Barrel_ladder_DZ * 0.5), slice_mat);
            ladderVol.setVisAttributes(desc,x_layer.visStr());
            sens.setType("tracker");
            ladderVol.setSensitiveDetector(sens);
            i_layer++;
            for (int ia = 0; ia < laddersCount; ia++) {
            phi = (ia * (cb_VTX_Barrel_ladder_deltaphi));
            x = - dR * cos(phi);
            y = - dR * sin(phi);
            RotationZYX ladder_rot = RotationZYX(cb_VTX_Barrel_ladder_deltaphi * ia, 0.0, 0.0);
            Position ladder_pos = Position(x, y, z);
            string ladderName = detName + _toString(i, "_ladder_Phys_%d") + _toString(ia, "_%d");
            PlacedVolume ladderPV = detVol.placeVolume(ladderVol, Transform3D(ladder_rot, ladder_pos));
            i_module++;
            ladderPV.addPhysVolID("layer", i_layer).addPhysVolID("module", i_module);
            }
            }
            // TODO: Pixels
            DetElement det(detName, detID);
            Volume motherVol = desc.pickMotherVolume(det);
            Transform3D tr(RotationZYX(0.0, 0.0, 0.0), Position(0.0, 0.0, pos.z()));
            PlacedVolume detPV = motherVol.placeVolume(detVol, tr);
            detPV.addPhysVolID("system", detID);
            detPV.addPhysVolID("barrel", 1);
            det.setPlacement(detPV);
            return det;
            }
            DECLARE_DETELEMENT(cb_VTX_Barrel, createDetector)
            0% or .
            You are about to add 0 people to the discussion. Proceed with caution.
            Finish editing this message first!
            Please register or to comment