From a59ef7798010080fad57a00c863c27b461d9b42a Mon Sep 17 00:00:00 2001
From: Kate Evans <ktevans@ifarm1901.jlab.org>
Date: Tue, 13 Feb 2024 12:40:18 -0500
Subject: [PATCH] I added an "Optics3" target that is a single carbon foil in
 the center of the z region of the target ladder. It it placed under the
 "Optics2" target in y on the ladder. I also edited Optics1 and Optics2 to be
 single foil targets that sample from the DSC target. Also, I grabbed the most
 up-to-date sieve design from Vassu Doomra.

---
 geometry/target/targetLadder.gdml |  50 +++-
 geometry/upstream/sieve.gdml      | 365 ++++++++++++++++++++----------
 macros/target/Optics2.mac         |   2 +-
 macros/target/Optics3.mac         |   3 +
 4 files changed, 297 insertions(+), 123 deletions(-)
 create mode 100644 macros/target/Optics3.mac

diff --git a/geometry/target/targetLadder.gdml b/geometry/target/targetLadder.gdml
index 485ae4357..ca057b68e 100644
--- a/geometry/target/targetLadder.gdml
+++ b/geometry/target/targetLadder.gdml
@@ -121,7 +121,7 @@
   <volume name="TargetAlDummy4pctUS_logical">
     <materialref ref="G4_Al"/>
     <solidref ref="TargetAlDummy4pct_solid"/>
-    <auxiliary auxtype="TargetSamplingVolume" auxvalue="USAl"/>
+    <auxiliary auxtype="TargetSamplingVolume" auxvalue="DSAl"/>
   </volume>
   <volume name="TargetAlDummy2pctDS_logical">
     <materialref ref="G4_Al"/>
@@ -175,38 +175,71 @@
     <solidref ref="TargetCFoil_solid"/>
     <auxiliary auxtype="TargetSamplingVolume" auxvalue="USC"/>
   </volume>
+
+  <volume name="TargetCFoilMS_logical">
+    <materialref ref="G4_C"/>
+    <solidref ref="TargetCFoil_solid"/>
+    <auxiliary auxtype="TargetSamplingVolume" auxvalue="MSC"/>
+  </volume>
+
   <volume name="TargetCFoilDS_logical">
     <materialref ref="G4_C"/>
     <solidref ref="TargetCFoil_solid"/>
     <auxiliary auxtype="TargetSamplingVolume" auxvalue="DSC"/>
   </volume>
+
+
   <volume name="TargetPositionOptics1_logical">
     <materialref ref="G4_Galactic"/>
     <solidref ref="TargetPosition_solid"/>
-    <physvol>
+    <!--<physvol>
       <volumeref ref="TargetCFoilUS_logical"/>
       <position z="-30*cm"/>
-    </physvol>
+    </physvol>-->
+
     <physvol>
       <volumeref ref="TargetCFoilDS_logical"/>
       <position z="+60*cm"/>
     </physvol>
+
     <auxiliary auxtype="TargetSystem" auxvalue="Optics1"/>
   </volume>
+
+
   <volume name="TargetPositionOptics2_logical">
     <materialref ref="G4_Galactic"/>
     <solidref ref="TargetPosition_solid"/>
-    <physvol>
+
+    <!--<physvol>
       <volumeref ref="TargetCFoilUS_logical"/>
       <position z="-60*cm"/>
-    </physvol>
+    </physvol>-->
+  
     <physvol>
       <volumeref ref="TargetCFoilDS_logical"/>
-      <position z="+30*cm"/>
+      <position z="-60*cm"/>
     </physvol>
+
     <auxiliary auxtype="TargetSystem" auxvalue="Optics2"/>
   </volume>
 
+  <volume name="TargetPositionOptics3_logical">
+    <materialref ref="G4_Galactic"/>
+    <solidref ref="TargetPosition_solid"/>
+
+    <!--<physvol>
+      <volumeref ref="TargetCFoilUS_logical"/>
+      <position z="-60*cm"/>
+    </physvol>-->
+  
+    <physvol>
+      <volumeref ref="TargetCFoilDS_logical"/>
+      <position z="0*cm"/>
+    </physvol>
+
+    <auxiliary auxtype="TargetSystem" auxvalue="Optics3"/>
+  </volume>
+
   <volume name="TargetLadder_logical">
     <materialref ref="G4_Galactic"/>
     <solidref ref="TargetLadder_solid"/>
@@ -255,6 +288,11 @@
       <position y="-63*cm"/>
     </physvol>
 
+    <physvol>
+      <volumeref ref="TargetPositionOptics3_logical"/>
+      <position y="-70*cm"/>
+    </physvol>
+
     <auxiliary auxtype="TargetLadder" auxvalue=""/>
   </volume>
 
diff --git a/geometry/upstream/sieve.gdml b/geometry/upstream/sieve.gdml
index 80775f63a..59ba8e713 100644
--- a/geometry/upstream/sieve.gdml
+++ b/geometry/upstream/sieve.gdml
@@ -5,35 +5,86 @@
 
 <define>
   <!-- Full sieve system dimensions -->
-  <quantity name="sieve_system_inner_radius" type="length" value="1.390" unit="in"/>
-  <quantity name="sieve_system_outer_radius" type="length" value="3.858" unit="in"/>
-  <quantity name="sieve_system_thickness" type="length" value="4.724" unit="in"/>
-
-  <!-- Sieve sector dimensions -->
-  <quantity name="sieve_sector_inner_radius" type="length" value="1.390" unit="in"/>
-  <quantity name="sieve_sector_outer_radius" type="length" value="3.858" unit="in"/>
-  <!-- Radiation length in W is 0.3504 cm [1], so 40 X_0 is about 15 cm.
-   [1] http://pdg.lbl.gov/2019/AtomicNuclearProperties/HTML/tungsten_W.html
-   -->
-  <quantity name="sieve_sector_thickness" type="length" value="4.724" unit="in"/>
-  <quantity name="sieve_sector_angle" type="angle" value="25.7" unit="deg"/>
+  <quantity name="sieve_system_inner_radius" type="length" value="26.5" unit="mm"/>
+  <quantity name="sieve_system_outer_radius" type="length" value="98" unit="mm"/>
+  <quantity name="sieve_system_thickness" type="length" value="100" unit="mm"/>
+
+  <quantity name="sieve_hole_diameter" type="length" value="5" unit="mm"/>
+
 
   <!-- Sieve hole dimensions -->
-  <quantity name="sieve_hole1_diameter" type="length" value="0.5" unit="in"/>
-  <quantity name="sieve_hole1_radial_position" type="length" value="1.8" unit="in"/>
-  <quantity name="sieve_hole1_angular_offset_phi" type="angle" value="0" unit="deg"/>
 
-  <quantity name="sieve_hole2_diameter" type="length" value="0.5" unit="in"/>
-  <quantity name="sieve_hole2_radial_position" type="length" value="3.5" unit="in"/>
-  <quantity name="sieve_hole2_angular_offset_phi" type="angle" value="-5" unit="deg"/>
+  <!-- sector 7 -->
+  <quantity name="sieve_hole71_radial_position" type="length" value="44" unit="mm"/>
+  <quantity name="sieve_hole71_angular_offset_phi" type="angle" value="0 + 360/14 - 360/7" unit="deg"/>
+
+  <quantity name="sieve_hole72_radial_position" type="length" value="68" unit="mm"/>
+  <quantity name="sieve_hole72_angular_offset_phi" type="angle" value="8 + 360/14 - 360/7" unit="deg"/>
+
+  <quantity name="sieve_hole73_radial_position" type="length" value="84.5" unit="mm"/>
+  <quantity name="sieve_hole73_angular_offset_phi" type="angle" value="-8 + 360/14 - 360/7" unit="deg"/>
+
+  <!-- sector 6 -->
+  <quantity name="sieve_hole61_radial_position" type="length" value="39" unit="mm"/>
+  <quantity name="sieve_hole61_angular_offset_phi" type="angle" value="0 - (360/7)*2 + 360/14" unit="deg"/>
+
+  <quantity name="sieve_hole62_radial_position" type="length" value="60" unit="mm"/>
+  <quantity name="sieve_hole62_angular_offset_phi" type="angle" value="8 - (360/7)*2 + 360/14" unit="deg"/>
+
+  <quantity name="sieve_hole63_radial_position" type="length" value="75" unit="mm"/>
+  <quantity name="sieve_hole63_angular_offset_phi" type="angle" value="-8 - (360/7)*2 + 360/14" unit="deg"/>
+
+  <!-- sector 5 -->
+  <quantity name="sieve_hole51_radial_position" type="length" value="56" unit="mm"/>
+  <quantity name="sieve_hole51_angular_offset_phi" type="angle" value="-8 + 360/14 -  (360/7)*3" unit="deg"/>
+
+  <quantity name="sieve_hole52_radial_position" type="length" value="56" unit="mm"/>
+  <quantity name="sieve_hole52_angular_offset_phi" type="angle" value="8 + 360/14 -  (360/7)*3" unit="deg"/>
+
+  <quantity name="sieve_hole53_radial_position" type="length" value="80" unit="mm"/>
+  <quantity name="sieve_hole53_angular_offset_phi" type="angle" value="0 + 360/14 -  (360/7)*3" unit="deg"/>
+
+  <!-- sector 4 -->
+  <quantity name="sieve_hole41_radial_position" type="length" value="50" unit="mm"/>
+  <quantity name="sieve_hole41_angular_offset_phi" type="angle" value="-8 + 360/14 - (360/7)*4" unit="deg"/>
+
+  <quantity name="sieve_hole42_radial_position" type="length" value="60" unit="mm"/>
+  <quantity name="sieve_hole42_angular_offset_phi" type="angle" value="0 + 360/14 - (360/7)*4" unit="deg"/>
+
+  <quantity name="sieve_hole43_radial_position" type="length" value="70" unit="mm"/>
+  <quantity name="sieve_hole43_angular_offset_phi" type="angle" value="8 + 360/14 - (360/7)*4" unit="deg"/>
+
+  <!-- sector 3 -->
+  <quantity name="sieve_hole31_radial_position" type="length" value="39" unit="mm"/>
+  <quantity name="sieve_hole31_angular_offset_phi" type="angle" value="0 + 360/14 - (360/7)*5" unit="deg"/>
+
+  <quantity name="sieve_hole32_radial_position" type="length" value="63" unit="mm"/>
+  <quantity name="sieve_hole32_angular_offset_phi" type="angle" value="8 + 360/14 - (360/7)*5" unit="deg"/>
+
+  <quantity name="sieve_hole33_radial_position" type="length" value="80" unit="mm"/>
+  <quantity name="sieve_hole33_angular_offset_phi" type="angle" value="-8 + 360/14 - (360/7)*5" unit="deg"/>
+
+  <!-- sector 2 -->
+  <quantity name="sieve_hole21_radial_position" type="length" value="50" unit="mm"/>
+  <quantity name="sieve_hole21_angular_offset_phi" type="angle" value="-8 + 360/14 - (360/7)*6" unit="deg"/>
+
+  <quantity name="sieve_hole22_radial_position" type="length" value="56" unit="mm"/>
+  <quantity name="sieve_hole22_angular_offset_phi" type="angle" value="8 + 360/14 - (360/7)*6" unit="deg"/>
+
+ <quantity name="sieve_hole23_radial_position" type="length" value="75" unit="mm"/>
+  <quantity name="sieve_hole23_angular_offset_phi" type="angle" value="0 + 360/14 - (360/7)*6" unit="deg"/>
+
+  <!-- sector 1 -->
+  <quantity name="sieve_hole11_radial_position" type="length" value="35" unit="mm"/>
+  <quantity name="sieve_hole11_angular_offset_phi" type="angle" value="0 + 360/14 - (360/7)*7" unit="deg"/>
+
+  <quantity name="sieve_hole13_radial_position" type="length" value="84.5" unit="mm"/>
+  <quantity name="sieve_hole13_angular_offset_phi" type="angle" value="8 + 360/14 - (360/7)*7" unit="deg"/>
+  
+  <quantity name="sieve_hole12_radial_position" type="length" value="58" unit="mm"/>
+  <quantity name="sieve_hole12_angular_offset_phi" type="angle" value="-8 + 360/14 - (360/7)*7" unit="deg"/>
 
-  <quantity name="sieve_hole3_diameter" type="length" value="0.5" unit="in"/>
-  <quantity name="sieve_hole3_radial_position" type="length" value="2.6" unit="in"/>
-  <quantity name="sieve_hole3_angular_offset_phi" type="angle" value="2.5" unit="deg"/>
-  <quantity name="sieve_hole3_axis_theta" type="angle" value="2.5" unit="deg"/>
-  <quantity name="sieve_hole3_axis_phi" type="angle" value="0" unit="deg"/>
 
-  <variable name="i" value="1"/>
 </define>
 
 <materials>
@@ -41,7 +92,7 @@
 
 
 <solids>
-  <!-- Sieve system solid -->
+ 
   <tube name="sieve_system_solid"
     rmin="sieve_system_inner_radius"
     rmax="sieve_system_outer_radius"
@@ -49,128 +100,210 @@
     startphi="0" deltaphi="360.0" aunit="deg"
     />
 
-  <!-- Sieve sector solid -->
-  <tube name="sieve_sector_solid"
-    rmin="sieve_sector_inner_radius"
-    rmax="sieve_sector_outer_radius"
-    z="sieve_sector_thickness"
-    startphi="-sieve_sector_angle/2" deltaphi="sieve_sector_angle"
-    />
-
   <!-- Sieve hole solids -->
-  <tube name="sieve_hole1_solid"
-    rmin="0."
-    rmax="sieve_hole1_diameter/2"
-    z="sieve_sector_thickness"
-    startphi="0" deltaphi="360.0" aunit="deg"
-    />
-  <tube name="sieve_hole2_solid"
-    rmin="0."
-    rmax="sieve_hole2_diameter/2"
-    z="sieve_sector_thickness"
-    startphi="0" deltaphi="360.0" aunit="deg"
-    />
-  <!-- Support for non-zero sieve_hole_axis_theta and sieve_hole_axis_phi
-    requires implementation of solid with cutTube, but this is slower so avoid
-    if not necessary -->
-  <cutTube name="sieve_hole3_solid"
+  <tube name="sieve_hole_solid"
     rmin="0."
-    rmax="sieve_hole3_diameter/2"
-    z="sieve_sector_thickness/cos(sieve_hole3_axis_theta)"
+    rmax="sieve_hole_diameter/2"
+    z="sieve_system_thickness+2"
     startphi="0" deltaphi="360.0" aunit="deg"
-    lowX="-sin(-sieve_hole3_axis_theta)*cos(-sieve_hole3_axis_phi)"
-    lowY="-sin(-sieve_hole3_axis_theta)*sin(-sieve_hole3_axis_phi)"
-    lowZ="-cos(-sieve_hole3_axis_theta)"
-    highX="sin(-sieve_hole3_axis_theta)*cos(-sieve_hole3_axis_phi)"
-    highY="sin(-sieve_hole3_axis_theta)*sin(-sieve_hole3_axis_phi)"
-    highZ="cos(-sieve_hole3_axis_theta)"
     />
+
 </solids>
 
 <structure>
-  <!-- Define the holes as detectors to enable mc-truth analysis where
-    one requires that the track went through a specific sieve hole -->
-  <volume name="sieve_hole1_logic">
-    <materialref ref="G4_Galactic"/>
-    <solidref ref="sieve_hole1_solid"/>
-    <auxiliary auxtype="Color" auxvalue="magenta"/>
-    <auxiliary auxtype="Alpha" auxvalue="0.1"/>
-    <auxiliary auxtype="SensDet" auxvalue="sieveDet"/>
-    <auxiliary auxtype="DetType" auxvalue="secondaries"/>
-    <auxiliary auxtype="DetType" auxvalue="boundaryhits"/>
-    <auxiliary auxtype="DetType" auxvalue="lowenergyneutral"/>
-    <auxiliary auxtype="DetNo" auxvalue="1001"/>
-  </volume>
-  <volume name="sieve_hole2_logic">
+  
+  <volume name="sieve_hole_logic">
     <materialref ref="G4_Galactic"/>
-    <solidref ref="sieve_hole2_solid"/>
+    <solidref ref="sieve_hole_solid"/>
     <auxiliary auxtype="Color" auxvalue="magenta"/>
     <auxiliary auxtype="Alpha" auxvalue="0.1"/>
     <auxiliary auxtype="SensDet" auxvalue="sieveDet"/>
     <auxiliary auxtype="DetType" auxvalue="secondaries"/>
     <auxiliary auxtype="DetType" auxvalue="boundaryhits"/>
     <auxiliary auxtype="DetType" auxvalue="lowenergyneutral"/>
-    <auxiliary auxtype="DetNo" auxvalue="1002"/>
   </volume>
-  <volume name="sieve_hole3_logic">
-    <materialref ref="G4_Galactic"/>
-    <solidref ref="sieve_hole3_solid"/>
-    <auxiliary auxtype="Color" auxvalue="magenta"/>
-    <auxiliary auxtype="Alpha" auxvalue="0.1"/>
+  
+  <volume name="sieve_system_logic">
+    <materialref ref="CW90"/>
+    <solidref ref="sieve_system_solid"/>
+    <auxiliary auxtype="Alpha" auxvalue="0.8"/>
     <auxiliary auxtype="SensDet" auxvalue="sieveDet"/>
     <auxiliary auxtype="DetType" auxvalue="secondaries"/>
     <auxiliary auxtype="DetType" auxvalue="boundaryhits"/>
     <auxiliary auxtype="DetType" auxvalue="lowenergyneutral"/>
-    <auxiliary auxtype="DetNo" auxvalue="1003"/>
-  </volume>
 
-  <!-- Place the holes inside the sector volume -->
-  <volume name="sieve_sector_logic">
-    <materialref ref="G4_W"/>
-    <solidref ref="sieve_sector_solid"/>
-    <physvol name="sieve_hole1_physical">
-      <volumeref ref="sieve_hole1_logic"/>
-      <position x="sieve_hole1_radial_position"/>
+    <physvol name="sieve_hole71_physical">
+      <volumeref ref="sieve_hole_logic"/>
+      <position
+        x="sieve_hole71_radial_position*cos(sieve_hole71_angular_offset_phi)"
+        y="sieve_hole71_radial_position*sin(sieve_hole71_angular_offset_phi)"
+        />
+    </physvol>
+
+    <physvol name="sieve_hole72_physical">
+      <volumeref ref="sieve_hole_logic"/>
+      <position
+        x="sieve_hole72_radial_position*cos(sieve_hole72_angular_offset_phi)"
+        y="sieve_hole72_radial_position*sin(sieve_hole72_angular_offset_phi)"
+        />
+
     </physvol>
 
-    <physvol name="sieve_hole2_physical">
-      <volumeref ref="sieve_hole2_logic"/>
+ <physvol name="sieve_hole73_physical">
+      <volumeref ref="sieve_hole_logic"/>
       <position
-        x="sieve_hole2_radial_position*cos(sieve_hole2_angular_offset_phi)"
-        y="sieve_hole2_radial_position*sin(sieve_hole2_angular_offset_phi)"
+        x="sieve_hole73_radial_position*cos(sieve_hole73_angular_offset_phi)"
+        y="sieve_hole73_radial_position*sin(sieve_hole73_angular_offset_phi)"
         />
     </physvol>
 
-    <physvol name="sieve_hole3_physical">
-      <volumeref ref="sieve_hole3_logic"/>
+   <physvol name="sieve_hole61_physical">
+      <volumeref ref="sieve_hole_logic"/>
       <position
-        x="sieve_hole3_radial_position*cos(sieve_hole3_angular_offset_phi)"
-        y="sieve_hole3_radial_position*sin(sieve_hole3_angular_offset_phi)"
+        x="sieve_hole61_radial_position*cos(sieve_hole61_angular_offset_phi)"
+        y="sieve_hole61_radial_position*sin(sieve_hole61_angular_offset_phi)"
         />
-      <rotation
-        y="-(sieve_hole3_axis_theta)"
+    </physvol>
+
+    <physvol name="sieve_hole62_physical">
+      <volumeref ref="sieve_hole_logic"/>
+      <position
+        x="sieve_hole62_radial_position*cos(sieve_hole62_angular_offset_phi)"
+        y="sieve_hole62_radial_position*sin(sieve_hole62_angular_offset_phi)"
         />
     </physvol>
-  </volume>
 
-  <!-- Place the sector volumes around the ring -->
-  <volume name="sieve_system_logic">
-    <materialref ref="G4_Galactic"/>
-    <solidref ref="sieve_system_solid"/>
-    <auxiliary auxtype="Alpha" auxvalue="0.1"/>
-    <auxiliary auxtype="SensDet" auxvalue="sieveDet"/>
-    <auxiliary auxtype="DetType" auxvalue="secondaries"/>
-    <auxiliary auxtype="DetType" auxvalue="boundaryhits"/>
-    <auxiliary auxtype="DetType" auxvalue="lowenergyneutral"/>
-    <auxiliary auxtype="DetNo" auxvalue="1000"/>
-    <loop for="i" from="1" to="7" step="1">
-      <physvol name="sieve_sector_physical[i]">
-        <volumeref ref="sieve_sector_logic"/>
-        <rotation unit="deg"
-          z="360.0/14.0+360.0/7.0*(i-1)"
-          />
-      </physvol>
-    </loop>
+   <physvol name="sieve_hole63_physical">
+      <volumeref ref="sieve_hole_logic"/>
+      <position
+        x="sieve_hole63_radial_position*cos(sieve_hole63_angular_offset_phi)"
+        y="sieve_hole63_radial_position*sin(sieve_hole63_angular_offset_phi)"
+        />
+    </physvol>
+
+   <physvol name="sieve_hole51_physical">
+      <volumeref ref="sieve_hole_logic"/>
+      <position
+        x="sieve_hole51_radial_position*cos(sieve_hole51_angular_offset_phi)"
+        y="sieve_hole51_radial_position*sin(sieve_hole51_angular_offset_phi)"
+        />
+    </physvol>
+
+    <physvol name="sieve_hole52_physical">
+      <volumeref ref="sieve_hole_logic"/>
+      <position
+        x="sieve_hole52_radial_position*cos(sieve_hole52_angular_offset_phi)"
+        y="sieve_hole52_radial_position*sin(sieve_hole52_angular_offset_phi)"
+        />
+    </physvol>
+    
+    <physvol name="sieve_hole53_physical">
+      <volumeref ref="sieve_hole_logic"/>
+      <position
+        x="sieve_hole53_radial_position*cos(sieve_hole53_angular_offset_phi)"
+        y="sieve_hole53_radial_position*sin(sieve_hole53_angular_offset_phi)"
+        />
+    </physvol>
+
+  <physvol name="sieve_hole41_physical">
+      <volumeref ref="sieve_hole_logic"/>
+      <position
+        x="sieve_hole41_radial_position*cos(sieve_hole41_angular_offset_phi)"
+        y="sieve_hole41_radial_position*sin(sieve_hole41_angular_offset_phi)"
+        />
+    </physvol>
+
+    <physvol name="sieve_hole42_physical">
+      <volumeref ref="sieve_hole_logic"/>
+      <position
+        x="sieve_hole42_radial_position*cos(sieve_hole42_angular_offset_phi)"
+        y="sieve_hole42_radial_position*sin(sieve_hole42_angular_offset_phi)"
+        />
+    </physvol>
+
+    <physvol name="sieve_hole43_physical">
+      <volumeref ref="sieve_hole_logic"/>
+      <position
+        x="sieve_hole43_radial_position*cos(sieve_hole43_angular_offset_phi)"
+        y="sieve_hole43_radial_position*sin(sieve_hole43_angular_offset_phi)"
+        />
+    </physvol>
+
+  <physvol name="sieve_hole31_physical">
+      <volumeref ref="sieve_hole_logic"/>
+      <position
+        x="sieve_hole31_radial_position*cos(sieve_hole31_angular_offset_phi)"
+        y="sieve_hole31_radial_position*sin(sieve_hole31_angular_offset_phi)"
+        />
+    </physvol>
+
+    <physvol name="sieve_hole32_physical">
+      <volumeref ref="sieve_hole_logic"/>
+      <position
+        x="sieve_hole32_radial_position*cos(sieve_hole32_angular_offset_phi)"
+        y="sieve_hole32_radial_position*sin(sieve_hole32_angular_offset_phi)"
+        />
+    </physvol>
+
+  <physvol name="sieve_hole33_physical">
+      <volumeref ref="sieve_hole_logic"/>
+      <position
+        x="sieve_hole33_radial_position*cos(sieve_hole33_angular_offset_phi)"
+        y="sieve_hole33_radial_position*sin(sieve_hole33_angular_offset_phi)"
+        />
+    </physvol>
+
+  <physvol name="sieve_hole21_physical">
+      <volumeref ref="sieve_hole_logic"/>
+      <position
+        x="sieve_hole21_radial_position*cos(sieve_hole21_angular_offset_phi)"
+        y="sieve_hole21_radial_position*sin(sieve_hole21_angular_offset_phi)"
+        />
+    </physvol>
+
+
+  <physvol name="sieve_hole22_physical">
+      <volumeref ref="sieve_hole_logic"/>
+      <position
+        x="sieve_hole22_radial_position*cos(sieve_hole22_angular_offset_phi)"
+        y="sieve_hole22_radial_position*sin(sieve_hole22_angular_offset_phi)"
+        />
+    </physvol>
+
+    <physvol name="sieve_hole23_physical">
+      <volumeref ref="sieve_hole_logic"/>
+      <position
+        x="sieve_hole23_radial_position*cos(sieve_hole23_angular_offset_phi)"
+        y="sieve_hole23_radial_position*sin(sieve_hole23_angular_offset_phi)"
+        />
+    </physvol>
+
+  <physvol name="sieve_hole11_physical">
+      <volumeref ref="sieve_hole_logic"/>
+      <position
+        x="sieve_hole11_radial_position*cos(sieve_hole11_angular_offset_phi)"
+        y="sieve_hole11_radial_position*sin(sieve_hole11_angular_offset_phi)"
+        />
+    </physvol>
+
+    <physvol name="sieve_hole12_physical">
+      <volumeref ref="sieve_hole_logic"/>
+      <position
+        x="sieve_hole12_radial_position*cos(sieve_hole12_angular_offset_phi)"
+        y="sieve_hole12_radial_position*sin(sieve_hole12_angular_offset_phi)"
+        />
+    </physvol>
+
+    <physvol name="sieve_hole13_physical">
+      <volumeref ref="sieve_hole_logic"/>
+      <position
+        x="sieve_hole13_radial_position*cos(sieve_hole13_angular_offset_phi)"
+        y="sieve_hole13_radial_position*sin(sieve_hole13_angular_offset_phi)"
+        />
+    </physvol>
+
+
+   
   </volume>
 </structure>
 
diff --git a/macros/target/Optics2.mac b/macros/target/Optics2.mac
index 2f0b72c07..ca3bc05de 100644
--- a/macros/target/Optics2.mac
+++ b/macros/target/Optics2.mac
@@ -1,3 +1,3 @@
 /remoll/geometry/absolute_position targetLadder (0,630,0)
 /remoll/target/mother Optics2
-/remoll/target/volume USC
+/remoll/target/volume DSC
diff --git a/macros/target/Optics3.mac b/macros/target/Optics3.mac
new file mode 100644
index 000000000..f7bd47078
--- /dev/null
+++ b/macros/target/Optics3.mac
@@ -0,0 +1,3 @@
+/remoll/geometry/absolute_position targetLadder (0,700,0)
+/remoll/target/mother Optics3
+/remoll/target/volume DSC