Compare revisions

1008a184 · 10e66bf1 · 10e66bf1 · 1008a184 · 1008a184 · 1008a184
--- a/calibrations/material-maps.json
+++ b/calibrations/material-maps.json
--- a/compact/cb_Solenoid.xml
+++ b/compact/cb_Solenoid.xml
-<lccdd>
-
-  <comment>Central Barrel Solenoid</comment>        
-
-  <define>
-  </define>
-
-  <limits>
-  </limits>
-
-  <regions>
-  </regions>
-
-  <display>
-  </display>
-
-  <!-- Define detector -->
-  <detectors>
-    <comment>
-      Central Barrel Solenoid Magnet
-    </comment>
-    <detector id="Solenoid_ID" name="cb_Solenoid" type="cb_Solenoid" insideTrackingVolume="false" vis="cb_SolenoidVis">
-	    <material name="Vacuum"/>  <!-- G4_Galactic -->	    
-      <dimensions rmin="cb_Solenoid_rmin" rmax="cb_Solenoid_rmax" z="cb_Solenoid_z" delta="cb_Solenoid_shift"/>
-    </detector>
-  </detectors>
-
-  <plugins>
-  </plugins>
-
-  <!-- TODO: Define field using field map-->
-  <fields>
-  </fields>
-</lccdd>
--- a/compact/ce_mrich.xml
+++ b/compact/ce_mrich.xml
-<?xml version="1.0" encoding="UTF-8"?>
-<lccdd>
-  <define>
-  </define>
-
-  <detectors>
-    <detector id="ce_MRICH_ID" name="ce_MRICH" type="refdet_ce_MRICH" readout="MRICHHits" vis="BlueVis" material="Air">
-      <dimensions rmin="ce_MRICHRMin" rmax="ce_MRICHRMax" length="ce_MRICHLength" z="ce_MRICHZMin"/>
-      <modules material="Quartz" thickness="10*cm" width="10*cm" gap="1*cm" />
-    </detector>
-  </detectors>
-
-  <readouts>
-    <readout name="MRICHHits">
-      <segmentation type="CartesianGridXY" grid_size_x="3*mm" grid_size_y="3*mm" />
-      <id>system:8,layer:4,piece:4,module:14,x:32:-16,y:-16</id>
-    </readout>
-  </readouts>
-
-</lccdd>
--- a/compact/central_tracker.xml
+++ b/compact/central_tracker.xml
--- a/compact/central_tracker_hybrid_v2.xml
+++ b/compact/central_tracker_hybrid_v2.xml
--- a/compact/colors.xml
+++ b/compact/colors.xml
+<display>
+
+  <comment> 
+    Colors tweaked to be displayable by geoViewer
+  </comment>
+  <vis name="AnlGreen"        alpha="1" r="100/256" g="200/256" b="0/256"/>
+  <vis name="AnlViolet"       alpha="1" r="100/256" g="0/256"   b="150/256"/>
+  <vis name="AnlProcess_Blue" alpha="1" r="0/256"   g="100/256" b="202/256"/>
+  <vis name="AnlOrange"       alpha="1" r="255/256" g="121/256" b="0"/>
+  <vis name="AnlRed"          alpha="1" r="205/256" g="45/256 " b="45/256"/>
+  <vis name="AnlGold"         alpha="1" r="248/256" g="188/256" b="0/256"/>
+  <vis name="AnlBlue"         alpha="1" r="0/256"   g="96/256"  b="156/256"/>
+  <vis name="AnlTeal"         alpha="1" r="0/256"   g="161/256" b="156/256"/>
+  <vis name="AnlGray"         alpha="1" r="102/256" g="102/256" b="102/256"/>
+  <vis name="AnlLight_Gray"   alpha="1" r="209/256" g="209/256" b="209/256"/>
+  <vis name="AnlLightGray"    ref="AnlLight_Gray" />
+  <!-- off-white (or white) don't work! -->
+  <!--vis name="AnlOff_White"    alpha="1" r="242/256" g="242/256" b="242/256"/-->
+  <vis name="AnlDarkRed"      alpha="1" r="161/256" g="45/256"  b="45/256"/>
+  <vis name="AnlDarkGreen"    alpha="1" r="0"       g="110/256" b="50/256"/>
+  <vis name="AnlDarkBlue"     alpha="1" r="45/256"   g="45/256"  b="110/256"/>
+  <vis name="AnlDelta_Red"    ref="AnlDarkRed"/>
+  <vis name="AnlDelta_Green"    ref="AnlDarkGreen"/>
+  <vis name="AnlDelta_Blue"    ref="AnlDarkBlue"/>
+
+</display>
--- a/compact/definitions.xml
+++ b/compact/definitions.xml
--- a/compact/dirc.xml
+++ b/compact/dirc.xml
+<lccdd>
+
+  <comment>
+    DIRC consists of 16 (default) identical modules making a barrel.
+    Each DIRC module consists of:
+    _______
+    [  |   \
+    [  |    \
+    [FD|Prizm|Lens|--------bars+glue-----|mirror]
+
+    Main constants to control DIRC global geometry:
+      - DIRC_length
+      - DIRC_offset
+      - DIRC_rmin
+      - DIRCBox_count  (16 full dirc, 1 - a single module)
+      - DIRC_rotation  (allows to flip the DIRC)
+
+    When DIRC_length is set, it affects bars length (all other lengths are fixed)
+    The width of each DIRC module is derived from DIRCPrism_width
+
+  </comment>
+
+  <define>
+    <constant name="DIRCFake_rmin"             value="BarrelTracking_rmax-11*cm"/>
+    <constant name="DIRC_rotation"         value="pi" comment="Allows DIRC flip"/>
+    <!-- <constant name="cb_DIRC_length"    value="DIRCLength"/> -->
+    <!-- CLEANUP THIS if not sure. This is from the initial implementation
+    <constant name="DIRC_length"            value="285.500*cm"/>
+    <constant name="DIRCFake_rmin"              value="82.00*cm"/>-->
+
+    <!-- Prism -->
+    <constant name="DIRCPrism_width"        value="360*mm"/>
+    <constant name="DIRCPrism_length"       value="300*mm"/>
+    <constant name="DIRCPrism_short_edge"   value="50*mm"/>
+    <constant name="DIRCPrism_angle"        value="32*deg"/>
+    <constant name="DIRCPrism_long_edge"    value="DIRCPrism_short_edge + DIRCPrism_length * tan(DIRCPrism_angle)"/>
+    <constant name="DIRCPrism_height"       value="DIRCPrism_long_edge"/>
+
+    <!-- DIRC length -->
+    <constant name="DIRCMain_length"           value="DIRC_length"/>
+    <constant name="DIRCFake_offset"           value="DIRC_offset+15*cm"/>
+
+    <!-- Box - main DIRC modules -->
+    <constant name="DIRCBox_count"          value="16"  comment="Number of DIRC boxes per... DIRC. 16 - default"/>
+    <constant name="DIRCBox_width"          value="DIRCPrism_width"/>
+
+    <!-- Mirror -->
+    <constant name="DIRCMirror_height"      value="20 * mm"/>
+    <constant name="DIRCMirror_width"       value="DIRCPrism_width"/>
+    <constant name="DIRCMirror_thickness"   value="1 * mm"/>
+
+    <!-- Lens -->
+    <constant name="DIRCLens_height"        value="DIRCPrism_long_edge"/>
+    <constant name="DIRCLens_width"         value="DIRCPrism_width"/>
+    <constant name="DIRCLens_thickness"     value="12 * mm"/>
+    <constant name="DIRCLens_r1"            value="33 * mm"/>
+    <constant name="DIRCLens_r2"            value="24 * mm"/>
+    <constant name="DIRCLens_shight"        value="25 * mm" comment="Is 'shight' from UK urban dictionary???"/>
+
+    <!-- FD (Foto detectors?) -->
+    <constant name="DIRCFd_height"          value="DIRCPrism_height"/>
+    <constant name="DIRCFd_width"           value="DIRCPrism_width"/>
+    <constant name="DIRCFd_thickness"       value="1*mm"/>
+
+    <!-- Bar - Each DIRC box consists of N "bars" -->
+    <!-- BarAssembly - Bars + Glue -->
+    <constant name="DIRCBarAssm_length"     value="DIRCMain_length - DIRCPrism_length - DIRCMirror_thickness - DIRCLens_thickness - DIRCFd_thickness" comment="Length of bars+glue assembly"/>
+    <constant name="DIRCBar_count"          value="11" comment="Number of bars per box"/>
+    <constant name="DIRCBar_gap"            value="0.15 * mm"/>
+    <constant name="DIRCBar_height"         value="17 * mm"/>
+    <constant name="DIRCBar_width"          value="(DIRCPrism_width - (DIRCBar_count - 1) * DIRCBar_gap) / DIRCBar_count"/>
+    <constant name="DIRCGlue_thickness"     value="0.05 * mm"/>
+    <constant name="DIRCBar_length"         value="(DIRCBarAssm_length-4*DIRCGlue_thickness)/4"/>
+
+    <!-- dirc outer volume -->
+    <comment>
+      RMax is calculated according to "rectangle inside ring" problem solution
+      https://math.stackexchange.com/questions/4222684/calculate-rectangle-inside-ring-parameters/4222691#4222691
+    </comment>
+    <constant name="DIRCFake_rmax"              value="sqrt( (DIRCFake_rmin+DIRCPrism_height)^2 + 0.25*(DIRCPrism_height)^2 )"/>
+  </define>
+
+  <materials>
+    <comment>
+      Only non optical material variants are here.
+      See optical_materials.xml for optical ones
+      -
+    </comment>
+
+    <material name="Epotek">
+      <D type="density" value="1.2" unit="g/cm3"/>
+      <composite n="5" ref="H"/>
+      <composite n="3" ref="C"/>
+      <composite n="2" ref="O"/>
+    </material>
+    <material name="Nlak33a">
+      <D type="density" value="4.220" unit="g/cm3"/>
+      <composite n="1" ref="Si"/>
+      <composite n="2" ref="O"/>
+    </material>
+  </materials>
+
+  <limits>
+  </limits>
+
+  <regions>
+  </regions>
+
+  <display>
+    <vis name="DIRCTube"   ref="AnlTeal" alpha="0.1" visible="true"  showDaughters="true"   lineStyle="solid" drawingStyle="solid" />
+    <vis name="DIRCBox"    ref="AnlLight_Gray" alpha="0.1" visible="true"  showDaughters="true"   lineStyle="solid" drawingStyle="solid" />
+    <vis name="DIRCLens1"  ref="AnlProcess_Blue" visible="true"  showDaughters="false"  lineStyle="solid" drawingStyle="solid" />
+    <vis name="DIRCLens2"  ref="AnlProcess_Blue" visible="true"  showDaughters="false"  lineStyle="solid" drawingStyle="solid" />
+    <vis name="DIRCLens3"  ref="AnlTeal" visible="true"  showDaughters="false"  lineStyle="solid" drawingStyle="solid" />
+    <vis name="DIRCBar"    ref="AnlTeal" visible="true"  showDaughters="false"  lineStyle="solid" drawingStyle="solid" />
+    <vis name="DIRCGlue"   ref="AnlViolet" visible="true"  showDaughters="false"  lineStyle="solid" drawingStyle="solid" />
+    <vis name="DIRCMirror" ref="AnlGray" visible="true"  showDaughters="false"  lineStyle="solid" drawingStyle="solid" />
+    <vis name="DIRCPrism"  ref="AnlTeal" visible="true"  showDaughters="false"  lineStyle="solid" drawingStyle="solid" />
+    <vis name="DIRCFd"     ref="AnlRed" visible="true"  showDaughters="false" />
+  </display>
+
+  <detectors>
+    <detector id="BarrelDIRC_ID" name="cb_DIRC" type="cb_DIRC" readout="DIRCBarHits" vis="DIRCTube">
+      <dimensions rmin="DIRCFake_rmin" rmax="DIRCFake_rmax" length="DIRCMain_length" radius="DIRCRadius" dx="DIRCBars_DX" dy="DIRCBars_DY" number="DIRCNum" deltaphi="DIRCDPhi"/>
+      <position x="0" y="0" z="DIRCFake_offset"/>
+      <rotation theta="DIRC_rotation"/>
+      <module name="DircBox" R="DIRCRadius" repeat="DIRCBox_count" width="DIRCPrism_width + 1*mm" height="DIRCPrism_height*2" length="DIRCBarAssm_length + 550*mm" vis="DIRCBox">
+        <!-- Mirror (at the end of the module) -->
+        <mirror
+          height="DIRCMirror_height"
+          width="DIRCMirror_width"
+          thickness="DIRCMirror_thickness"
+          material="Aluminum"
+          vis="DIRCMirror"
+        />
+
+        <!-- DIRC Bar -->
+        <!-- TODO parametrize as bar+glue assembly with layers -->
+        <bar
+          height="DIRCBar_height"
+          width="DIRCBar_width"
+          length="DIRCBar_length"
+          repeat="DIRCBar_count"
+          gap="DIRCBar_gap"
+          material="Quartz"
+          vis="DIRCBar"
+        />
+
+        <glue thickness="DIRCGlue_thickness" material="Epotek" vis="DIRCGlue"/>
+
+        <!-- lens -->
+        <lens
+          height="DIRCLens_height"
+          width="DIRCLens_width"
+          thickness="DIRCLens_thickness"
+          r1="DIRCLens_r1"
+          r2="DIRCLens_r2"
+        />
+
+        <prism/> <!-- TODO parametrize here -->
+
+        <fd
+          height="DIRCFd_height"
+          width="DIRCFd_width"
+          thickness="DIRCFd_thickness"
+          vis="DIRCFd"
+          material="Quartz"
+        />
+      </module>
+    </detector>
+  </detectors>
+
+  <readouts>
+    <readout name="DIRCBarHits">
+      <segmentation type="CartesianGridXY" grid_size_x="3.0*mm" grid_size_y="3.0*mm" />
+      <id>system:8,module:8,section:8,bar:8,x:32:-16,y:-16</id>
+    </readout>
+  </readouts>
+
+  <plugins>
+  </plugins>
+
+  <fields>
+  </fields>
+</lccdd>
--- a/compact/display.xml
+++ b/compact/display.xml
--- a/compact/display_detailed.xml
+++ b/compact/display_detailed.xml
--- a/compact/display_geoviewer.xml
+++ b/compact/display_geoviewer.xml
--- a/compact/drich.xml
+++ b/compact/drich.xml
--- a/compact/ecal.xml
+++ b/compact/ecal.xml
--- a/compact/ecal_backward_hybrid.xml
+++ b/compact/ecal_backward_hybrid.xml
+<lccdd>
+  <define>
+    <comment>
+      Transition area.
+      The idea behind this parametrization is that:
+      one glass module with its wrap is always
+      a size of 4 crystal modules with its wraps.
+      Then the transition area (where glass meets crystals) has no gaps
+
+      +----------------+----------------+
+      | +----+  +----+ |  +----------+  |
+      | |    |  |    | |  |          |  |
+      | +----+  +----+ |  |          |  |
+      | +----+  +----+ |  |          |  |
+      | |    |  |    | |  |          |  |
+      | +----+  +----+ |  +----------+  |
+      +----------------+----------------+
+           crystal           glass
+
+      This implies that:
+        GlassModule_wrap = 2*CrystalModule_wrap
+        GlassModule_sx = 2*CrystalModule_sx
+        GlassModule_sy = 2*CrystalModule_sy
+
+    </comment>
+
+    <constant name="CrystalModule_width" value="20.00*mm"/>
+    <constant name="CrystalModule_length" value="200.00*mm"/>
+    <constant name="CrystalModule_wrap" value="0.50*mm"/>
+    <constant name="CrystalModule_z0" value="0.0*cm"/>
+
+    <constant name="GlassModule_width" value="2*CrystalModule_width"/>
+    <constant name="GlassModule_length" value="55.00*cm"/>
+    <constant name="GlassModule_wrap" value="2*CrystalModule_wrap"/>
+    <constant name="GlassModule_z0" value="0.0*cm"/>
+
+    <constant name="EcalEndcapNIonCutout_dphi" value="30*degree"/>
+
+    <constant name="EcalEndcapN_thickness" value="GlassModule_length"/>
+    <constant name="EcalEndcapN_z0" value="-EcalEndcapN_zmin - EcalEndcapN_thickness/2"/>
+    <constant name="EcalEndcapNCrystal_rmax" value="40*cm"/>
+
+    <constant name="CrystalModule_distance" value="CrystalModule_width + CrystalModule_wrap"/>
+    <constant name="GlassModule_distance" value="GlassModule_width + GlassModule_wrap"/>
+  </define>
+
+  <display>
+    <vis name="HybridEcalOuterVis" alpha="0.5"  r= "0.3"  g="0.3"  b="0.3"  showDaughters="true" visible="true"/>
+  </display>
+
+  <detectors>
+
+    <documentation level="10">
+      #### Backwards Endcap EM Calorimeter
+
+      Backwards Endcap EM Calorimeter, placements generated by script
+    </documentation>
+    <detector id="ECalEndcapN_ID" name="EcalEndcapN" type="HybridCalorimeter" readout="EcalEndcapNHits">
+      <position x="0" y="0" z="EcalEndcapN_z0"/>
+      <rotation x="0" y="0" z="0"/>
+    </detector>
+  </detectors>
+  <readouts>
+    <comment>Effectively no segmentation, the segmentation is used to provide cell dimension info</comment>
+    <readout name="EcalEndcapNHits">
+      <segmentation type="MultiSegmentation" key="sector">
+        <segmentation name="CrystalSeg" key_value="1" type="CartesianGridXY" grid_size_x="CrystalModule_distance" grid_size_y="CrystalModule_distance"/>
+        <segmentation name="GlassSeg" key_value="2" type="CartesianGridXY" grid_size_x="GlassModule_distance" grid_size_y="GlassModule_distance"/>
+      </segmentation>
+      <id>system:8,sector:4,module:20,x:32:-16,y:-16</id>
+    </readout>
+  </readouts>
+</lccdd>
--- a/compact/ecal_barrel.xml
+++ b/compact/ecal_barrel.xml
--- a/compact/ecal_barrel_interlayers.xml
+++ b/compact/ecal_barrel_interlayers.xml
--- a/compact/ecal_forward_homogenous.xml
+++ b/compact/ecal_forward_homogenous.xml
--- a/compact/ecal_forward_scfi.xml
+++ b/compact/ecal_forward_scfi.xml
--- a/compact/fake_dirc.xml
+++ b/compact/fake_dirc.xml
--- a/compact/ffi_ZDC.xml
+++ b/compact/ffi_ZDC.xml
-<lccdd>
-
-  <comment>Far Forward Ion Zero Degree Calorimeter</comment>        
-
-  <define>
-  </define>
-
-  <limits>
-  </limits>
-
-  <regions>
-  </regions>
-
-  <display>
-  </display>
-
-  <detectors>
-    <detector id="ffi_ZDC_ID" name="ffi_ZDC" type="ffi_ZDC" readout="ffi_ZDC_Hits" vis="ffi_ZDCVis">
-      <position x="ffi_ZDC_x_pos" y="ffi_ZDC_y_pos" z="ffi_ZDC_z_pos"/>
-      <rotation x="ffi_ZDC_rotateX_angle" y="ffi_ZDC_rotateY_angle" z="ffi_ZDC_rotateZ_angle"/>
-      <dimensions x="ffi_ZDC_width" z="ffi_ZDC_thickness"/>
-      <module name="ffi_ZDC_ECAL_module" vis="ffi_ZDCmoduleVis" material="PbWO4" thickness="ffi_ZDC_ECAL_thickness" rmin="ffi_ZDC_ECAL_inner_radius" rmax="ffi_ZDC_ECAL_outer_radius" width="ffi_ZDC_ECAL_width" gap="ffi_ZDC_ECAL_gap" ntower="ffi_ZDC_ECAL_ntower"/>
-    </detector>
-  </detectors>
-
-  <readouts>
-    <readout name="ffi_ZDC_Hits">
-      <segmentation type="CartesianGridXY" grid_size_x="1.0*mm" grid_size_y="1.0*mm" />
-      <id>system:8,sector:6,module:6,x:32:-16,y:-16</id>  
-    </readout>
-  </readouts>
-
-  <plugins>
-  </plugins>
-
-  <fields>
-  </fields>
-</lccdd>
No results found