diff --git a/compact/cb_CTD_Si.xml b/compact/cb_CTD_Si.xml index 3505aa3fad614d89d53223692b46031734b668f9..06b4f825556ef7b826860bef188cdf11521ec55d 100644 --- a/compact/cb_CTD_Si.xml +++ b/compact/cb_CTD_Si.xml @@ -15,8 +15,8 @@ </display> <detectors> - <detector id="cb_CTD_Si_ID" name="cb_CTD_Si" type="cb_CTD_Si" readout="cb_CTD_Si_Hits" vis="cb_CTD_SiVis" insideTrackingVloume="true"> - <dimensions rmin="cb_CTD_Si_rmin" rmax="cb_CTD_Si_rmax" length="cb_CTD_Si_length" zmax="cb_CTD_Si_zmax" gap="cb_CTD_Si_gap"/> + <detector id="cb_CTD_Si_ID" name="cb_CTD_Si" type="cb_CTD_Si" readout="cb_CTD_Si_Hits" vis="cb_CTDVis" insideTrackingVloume="true"> + <dimensions rmin="cb_CTD_rmin" rmax="cb_CTD_rmax" length="cb_CTD_length" zmax="cb_CTD_zmax" gap="cb_CTD_Si_gap"/> <layer repeat="cb_CTD_Si_layer" vis="cb_CTD_Si_layerVis"> <slice name="Silicon_slice" material="Si" sensitive="true"/> </layer> diff --git a/compact/definitions.xml b/compact/definitions.xml index df5f8113072b810cb1179d34cc9c462869f8609e..1f454dc65c626784ec861a8f1a4054df08c30999 100644 --- a/compact/definitions.xml +++ b/compact/definitions.xml @@ -321,22 +321,23 @@ <constant name="SiliconTrackerOuterBarrelLength" value="SolenoidLength/2.0"/> <comment> - ---------------------------- - Central Barrel Tracker - Si - ---------------------------- + ------------------------------- + Central Barrel Tracker - Global + ------------------------------- </comment> - <constant name="cb_CTD_Si_rmin" value="21.0 * cm"/> - <constant name="cb_CTD_Si_rmax" value="80.0 * cm"/> - <constant name="cb_CTD_Si_zmax" value="200.0 * cm"/> + <constant name="cb_CTD_rmin" value="21.0 * cm"/> + <constant name="cb_CTD_rmax" value="80.0 * cm"/> + <constant name="cb_CTD_zmax" value="200.0 * cm"/> <comment> Original Global parameter </comment> <!-- - <constant name="cb_CTD_Si_zmax" value="120.0 * cm"/> + <constant name="cb_CTD_zmax" value="120.0 * cm"/> --> - <constant name="cb_CTD_Si_length" value="SolenoidLength - cb_CTD_Si_zmax"/> + <constant name="cb_CTD_length" value="SolenoidLength - cb_CTD_zmax"/> <comment> After fix cb_Solenoid parameters </comment> <!-- - <constant name="cb_CTD_Si_length" value="cb_Solenoid_z - cb_CTD_Si_zmax"/> + <constant name="cb_CTD_length" value="cb_Solenoid_z - cb_CTD_zmax"/> --> + <comment> Si layer opiton </comment> <constant name="cb_CTD_Si_gap" value="5.0 * cm"/> <constant name="cb_CTD_Si_layer" value="15"/> diff --git a/compact/display.xml b/compact/display.xml index 2a7fa01def07962300863523683fb5ea8956be09..81130895a31b7a1ae52aa2daa5af2ce38d196641 100644 --- a/compact/display.xml +++ b/compact/display.xml @@ -51,7 +51,7 @@ <vis name="cb_SolenoidVis" alpha="0.4" r= "0.1" g="0.0" b="0.1" showDaughters="true" visible="true"/> <vis name="ffi_ZDCVis" alpha="0.1" r= "0.1" g="0.0" b="1.0" showDaughters="true" visible="true"/> <vis name="ffi_ZDCmoduleVis" alpha="1.0" r= "0.1" g="1.0" b="0.9" showDaughters="true" visible="true"/> - <vis name="cb_CTD_SiVis" alpha="0.1" r= "0.1" g="0.0" b="1.0" showDaughters="true" visible="true"/> + <vis name="cb_CTDVis" alpha="0.1" r= "0.1" g="0.0" b="1.0" 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"/> <comment> diff --git a/src/cb_CTD_Si.cpp b/src/cb_CTD_Si.cpp index d0f7025f9c3246f31a578252bfe8ceffe8947e50..4aed32b4ac10437149b7555c171b3f2253c9fa31 100644 --- a/src/cb_CTD_Si.cpp +++ b/src/cb_CTD_Si.cpp @@ -26,11 +26,41 @@ static Ref_t createDetector(Detector& desc, xml_h e, SensitiveDetector sens) Volume detVol("cb_CTD_GVol_Logic", cb_CTD_GVol_Solid, Vacuum); detVol.setVisAttributes(desc.visAttributes(x_det.visStr())); + // Construct Silicon Layers + 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 layerRIn[100]; + double layerROut[100]; + + // Loop over layers + for(int i = 0; i < repeat; i++) { + layerRIn[i] = RIn + (SiLayerGap * i); + layerROut[i] = RIn + (0.01 + SiLayerGap * i); + + if (layerROut[i] > ROut) + continue; + + string logic_layer_name = detName + _toString(i, "_Logic_lay_%d"); + Volume layerVol(logic_layer_name,Tube(layerRIn[i], layerROut[i], SizeZ / 2.0, 0.0, 360.0 * deg), slice_mat); + layerVol.setVisAttributes(desc,x_layer.visStr()); + sens.setType("tracker"); + layerVol.setSensitiveDetector(sens); + + Position layer_pos = Position(0.0, 0.0, 0.0); + PlacedVolume layerPV = detVol.placeVolume(layerVol, layer_pos); + layerPV.addPhysVolID("layer", i+1); + } + DetElement det(detName, detID); Volume motherVol = desc.pickMotherVolume(det); Transform3D tr(RotationZYX(0.0, 0.0, 0.0), Position(0.0, 0.0, 0.0)); PlacedVolume detPV = motherVol.placeVolume(detVol, tr); detPV.addPhysVolID("system", detID); + detPV.addPhysVolID("barrel", 1); det.setPlacement(detPV); return det; }