Modified CMakeLists.txt and Added ZeroDegreeCalorimeter detector.

src/GenericDetectors/CMakeLists.txt

@@ -39,6 +39,7 @@ dd4hep_add_plugin(GenDetectors
 	calorimeters/src/EndcapECAL_geo.cpp
         calorimeters/src/CrystalEndcapECAL_geo.cpp
         calorimeters/src/ScintillatingTileEndcapHCAL_geo.cpp
+	calorimeters/src/ZeroDegreeCAL_geo.cpp
 	beamline/src/Beampipe_geo.cpp
         beamline/src/IRChamber_geo.cpp
         beamline/src/Beampipe_geo.cpp

src/GenericDetectors/calorimeters/compact/ZeroDegreeCAL_example.xml

+<lccdd xmlns:compact="http://www.lcsim.org/schemas/compact/1.0" 
+       xmlns:xs="http://www.w3.org/2001/XMLSchema" 
+       xs:noNamespaceSchemaLocation="http://www.lcsim.org/schemas/compact/1.0/compact.xsd">
+
+  <!-- Some information about detector  -->
+  <info name="ZDC_example" title="Zero Degree Calorimeter detector example"
+        author="Jihee Kim"
+        url="https://eicweb.phy.anl.gov/EIC/NPDet"
+        status="development"
+        version="v1.1 2020-05-14">
+    <comment>Zero Degree Calorimeter detector</comment>        
+  </info>
+
+  <!-- Use DD4hep elements and materials definitions -->
+  <includes>
+    <gdmlFile  ref="elements.xml"/>
+    <gdmlFile  ref="materials.xml"/>
+  </includes>
+
+  <!-- Define the dimensions of the world volume -->
+  <define>
+    <constant name="world_side" value="10*m"/>
+    <constant name="world_x" value="world_side"/>
+    <constant name="world_y" value="world_side"/>
+    <constant name="world_z" value="world_side"/>
+    
+    <constant name="CrossingAngle" value="0.020*rad"/>
+    <constant name="tracker_region_zmax" value="5*m"/>
+    <constant name="tracker_region_rmax" value="5*m"/>
+  </define>
+
+  <limits>
+    <limitset name="cal_limits">
+      <limit name="step_length_max" particles="*" value="5.0" unit="mm" />
+    </limitset>
+    <limitset name="GenericRICHRegionLimitSet">
+      <limit name="step_length_max" particles="*" value="1.0" unit="mm" />
+      <limit name="track_length_max" particles="*" value="1.0" unit="mm" />
+      <limit name="time_max" particles="*" value="0.1" unit="ns" />
+      <limit name="ekin_min" particles="*" value="0.001" unit="MeV" />
+      <limit name="range_min" particles="*" value="0.1" unit="mm" />
+    </limitset>
+  </limits>
+
+  <regions>
+    <region name="GenericRICHRegion" eunit="MeV" lunit="mm" cut="0.0001" threshold="0.0001">
+      <limitsetref name="GenericRICHRegionLimitSet"/>
+    </region>
+  </regions>
+
+  <!-- Common Generic visualization attributes -->
+  <comment>Common Generic visualization attributes</comment>
+  <display>
+    <vis name="InvisibleNoDaughters"      showDaughters="false" visible="false"/>
+    <vis name="InvisibleWithDaughters"    showDaughters="true" visible="false"/>
+    <vis name="GreenVis"       alpha="0.5"  r= "0.0" g="1.0" b="0.0" showDaughters="true" visible="true"/>
+    <vis name="RedVis"         alpha="0.5"  r= "1.0" g="0.0" b="0.0" showDaughters="true" visible="true"/>
+    <vis name="BlueVis"        alpha="0.5"  r= "0.0" g="0.0" b="1.0" showDaughters="true" visible="true"/>
+    <vis name="OrangeVis"      alpha="0.5"  r= "1.0" g="0.45" b="0.0" showDaughters="true" visible="true"/>
+    <vis name="RedGreenVis"    alpha="0.5"  r= "1.0" g="1.0" b="0.0" showDaughters="true" visible="true"/>
+    <vis name="BlueGreenVis"   alpha="0.5"  r= "0.0" g="1.0" b="1.0" showDaughters="true" visible="true"/>
+    <vis name="PurpleVis"      alpha="0.5"  r= "1.0" g="0.0" b="1.0" showDaughters="true" visible="true"/>
+    <vis name="DoubleRedG"     alpha="0.5"  r= "2.0" g=".10" b="0.0" showDaughters="true" visible="true"/>
+    <vis name="RBG015"         alpha="0.5"  r= "0.0" g=".2"  b="1.0" showDaughters="true" visible="true"/>
+    <vis name="RBG510"         alpha="0.5"  r= "1.0" g=".2"  b="0.0" showDaughters="true" visible="true"/>
+    <vis name="RBG"            alpha="0.5"  r= "1.0" g="1.0" b="1.0" showDaughters="true" visible="true"/>
+    <vis name="GrayVis"        alpha="0.5"  r= "0.75" g="0.75" b="0.75" showDaughters="true" visible="true"/>
+  </display>
+
+  <!-- Define detector -->
+  <detectors>
+    <detector id="1" name="ZDC" type="ZeroDegreeCAL" readout="ZDCHits" vis="RedVis">
+      <module name="diamond_module" side="40*mm" layers="16"> 
+        <slice name="Tungsten_slice" thickness="7*mm" number="0" material="TungstenDens24"/>
+	<slice name="Scint_slice" thickness="3*mm" number="1" material="PlasticScint"/>
+	<slice name="SciFi_belt" thickness="1*mm" number="2" material="PlasticScint"/>
+	<dimensions rmin="0*cm" rmax="40*cm" zmin="-10*cm"  zmax="40*cm"/>
+      </module>
+    </detector>
+  </detectors>
+
+  <!--  Definition of the readout segmentation/definition  -->
+  <readouts>
+          <readout name="ZDCHits">
+		<segmentation type="CartesianGridXY" grid_size_x="3.5*cm" grid_size_y="3.5*cm" />
+		<id>system:8,module:14,layer:10,slice:2,x:48:-8,y:-8</id>  
+    </readout>
+  </readouts>
+
+  <plugins>
+    <!--
+    <plugin name="DD4hep_GenericSurfaceInstallerPlugin">
+      <argument value="ForwardRomanPot"/>
+      <argument value="dimension=2"/>
+      <argument value="u_x=-1."/>
+      <argument value="v_y=-1."/>
+      <argument value="n_z=1."/>
+    </plugin>	
+    <plugin name="InstallSurfaceManager"/>
+    -->
+  </plugins>
+
+  <fields>
+    <field name="GlobalSolenoid" type="solenoid" 
+      inner_field="4.0*tesla"
+      outer_field="-0.6*tesla" 
+      zmax="3*m"
+      outer_radius="2*m">
+    </field>
+  </fields>
+</lccdd>

src/GenericDetectors/calorimeters/compact/macro/gps_zdc.mac

+/run/beamOn 1
+/control/verbose 2
+/run/initialize
+
+/gps/verbose 2
+/gps/particle neutron 
+/gps/number 1
+
+/gps/ene/type Gauss
+/gps/ene/mono  2.5 GeV
+/gps/ene/sigma 2.0 GeV
+
+/gps/pos/type Volume
+/gps/pos/shape Cylinder
+/gps/pos/centre 0.0 0.0 0.0 cm
+/gps/pos/radius 0.01 cm
+/gps/pos/halfz  0.01 cm
+/gps/position 0 0 -1.0 cm
+
+/gps/direction 0.0 0.5 1.0 
+#/gps/ang/type iso
+
+
+/control/execute macro/vis_zdc.mac
+
+/run/beamOn 1

src/GenericDetectors/calorimeters/compact/macro/vis_zdc.mac

+/vis/open OGL 800x800-0+0
+
+/vis/drawVolume
+/vis/viewer/set/viewpointThetaPhi 270. 0.
+#/vis/viewer/zoom 30.
+#/vis/viewer/set/style wireframe
+#/vis/scene/add/axes 0 0 0 1 m
+/vis/scene/add/trajectories rich smooth
+/vis/modeling/trajectories/create/drawByCharge
+/vis/modeling/trajectories/drawByCharge-0/default/setDrawStepPts true
+/vis/modeling/trajectories/drawByCharge-0/default/setStepPtsSize 2
+/vis/scene/add/hits
+/vis/scene/endOfEventAction accumulate 10 
+
+/vis/viewer/flush
+
+
+

src/GenericDetectors/calorimeters/compact/scripts/run_example3

@@ -3,5 +3,5 @@
 ddsim  --runType vis  --enableG4GPS \
    --macroFile macro/gps3.mac \
    --compactFile $1 \
-   --outputFile test_calorimeter_sampling.root  || exit
+   --outputFile test_test.root  || exit
    #--outputFile test_tracker_disc.root  || exit

src/GenericDetectors/calorimeters/compact/scripts/run_example_zdc

+#!/bin/bash
+
+ddsim  --runType vis  --enableG4GPS \
+   --macroFile macro/gps_zdc.mac \
+   --compactFile $1 \
+   --outputFile test_zdc.root  || exit
+   #--outputFile test_tracker_disc.root  || exit