LISTSERV 16.5 - LCD-CVS Archives

+ * <p> + * Within that part, each sensor is identified by a unique non-negative integer + * <tt>postfix</tt>. Subclusses should override three methods to define a particular + * segmentation: {@link #getSensorPostfix(SimTrackerHit hit)} that returns <tt>postfix</tt> + * given the simulated hit position, {@link #getMaxPostfix()} that returns maximum + * possible postfix value, and {@link #makeSensor(int postfix)} that creates a new + * {@link Sensor} objects given <tt>postfix</tt>. + * <p> + * {@link SegmentationManager} will assign <tt>prefix</tt> value and overall ID length + * to this <tt>Segmenter</tt> through calls to {@link #setPrefix(int prefix)} and + * {@link #setPostfixLength(int length)}. + * After that, {@link #getSensorID(SimTrackerHit hit)} and {@link #getSensor(int sensorID)} + * can be used (again, by {@link SegmentationManager}) to calculate <tt>SensorID</tt> + * (unique within the detector) and to create a new {@link Sensor}. + * <p> + * If a subclass requires any detector dependent initialization, it should override + * {@link #detectorChanged(Detector)} but include a call to <tt>super.detectorChanged(Detector)</tt>.

- * Returns sensor ID postfix corresponding to the given position. This postfix must - * be positive integer, unique within the part of the detector handled by this - * <tt>Segmenter</tt> object. The final Sensor ID will be constructed taking into - * account the prefix set through a call to {@link #setPrefix()} method.

+ * Returns <tt>postfix</tt> corresponding to the position of the given simulated + * hit. The <tt>postfix</tt> must be non-negative integer, + * unique within the part of the detector handled by this <tt>Segmenter</tt> + * object. The final <tt>SensorID</tt> will be constructed taking into account + * the prefix set through a call to {@link #setPrefix(int prefix)} method.

+ System.out.println("Hit layer: " + hit.getLayer() +" Sensor ID: " + ((_prefix << _postfixLength) | getSensorPostfix(hit)) + + " Prefix: " + _prefix + " Postfix " + getSensorPostfix(hit));

* For any ID returned by {@link #getSensorID(SimTrackerHit)} method, * the number of lower bits set throught a call to {@link #setPostfixLength(int)} * will contain the postfix returned by {@link #getSensorPostfix(SimTrackerHit)},

* For any ID returned by {@link #getSensorID(SimTrackerHit)} method, * the number of lower bits set throught a call to this method will contain the * postfix returned by {@link #getSensorPostfix(SimTrackerHit)}, preceded by

+ * Simplistic {@link Segmenter} that tiles barrel cylinders with Z-parallel strips or pixels. + * <p> + * Each barrel layer will correspond to a single {@link Sensor} object, with the + * <tt>postfix</tt> equal to layer number.

+ + /** + * Creates a new instance of CylindricalBarrelSegmenter. + * Subdetector name supplied to the constructor is used to provide reasonable + * defaults for strip width and length. + */

public void detectorChanged(Detector detector) { super.detectorChanged(detector); List<IDetectorElement> sensors = SegmentationManager.getLeaves(detector.getSubdetector(_sdName).getDetectorElement());

+ Collections.sort(sensors, new Comparator<IDetectorElement>() { + public int compare(IDetectorElement s1, IDetectorElement s2) { + return (int)Math.signum(((Tube)(s1.getGeometry().getLogicalVolume().getSolid())).getInnerRadius() + - ((Tube)(s2.getGeometry().getLogicalVolume().getSolid())).getInnerRadius()); + } + });

for (IDetectorElement del : sensors) { IGeometryInfo gInfo = del.getGeometry(); Tube solid = (Tube) gInfo.getLogicalVolume().getSolid(); double rInner = solid.getInnerRadius(); double rOuter = solid.getOuterRadius();

lcsim/src/org/lcsim/contrib/onoprien/tracking/geom/segmentation

RingSegmenter.java 1.4 -> 1.5

diff -u -r1.4 -r1.5
--- RingSegmenter.java	31 May 2007 09:38:54 -0000	1.4
+++ RingSegmenter.java	12 Jul 2007 05:00:22 -0000	1.5
@@ -1,6 +1,6 @@

 package org.lcsim.contrib.onoprien.tracking.geom.segmentation;

-import java.util.List;

+import java.util.*;

 
 import hep.physics.vec.BasicHep3Vector;
 import hep.physics.vec.Hep3Vector;

@@ -24,64 +24,90 @@

 
 /**
  *

- * Fix relation between layer returned by hits and copy number.

+ * Simplistic {@link Segmenter} that tiles endcap disks with strips or pixels.
+ * <p>
+ * Each disk will correspond to a {@link Sensor} object, with the <tt>postfix</tt>
+ * equal to layer number for disks in Z-positive endcap, and to layer number plus
+ * number of layers in Z-negative endcap. Disks corresponding to odd-numbered layers
+ * are rotated by an angle set through a call to {@link #setStereoAngle(double angle)}.

  *
  * @author D.Onoprienko

- * @version $Id: RingSegmenter.java,v 1.4 2007/05/31 09:38:54 onoprien Exp $

+ * @version $Id: RingSegmenter.java,v 1.5 2007/07/12 05:00:22 onoprien Exp $

  */
 public class RingSegmenter extends Segmenter {
   
 // -- Constructors and initialization :  ---------------------------------------

-  /** Creates a new instance of CylindricalBarrelSegmenter */

+  /** 
+   * Creates a new instance of RingSegmenter.
+   * Subdetector name supplied to the constructor is used to provide reasonable 
+   * defaults for strip width, length, and stereo angle.
+   */

   public RingSegmenter(String subdetectorName) {
     _sdName = subdetectorName;
     if (_sdName == "TrackerEndcap") {
       setStripWidth(25.*Const.micrometer);
       setStripLength(10.*Const.cm);

+      setStereoAngle(Math.PI / 2.);

     } else if (_sdName == "VertexEndcap" || _sdName == "TrackerForward") {
       setStripWidth(25.*Const.micrometer);
       setStripLength(25.*Const.micrometer);

+      setStereoAngle(0.);

-    setStereoAngle(Math.PI / 2.);

   /** 
    * Detector-dependent initialization.
    */
   public void detectorChanged(Detector detector) {
     super.detectorChanged(detector);
     List<IDetectorElement> sensors = SegmentationManager.getLeaves(detector.getSubdetector(_sdName).getDetectorElement());

-    _nDisks = sensors.size();

+    Collections.sort(sensors, new Comparator<IDetectorElement>() {
+      public int compare(IDetectorElement s1, IDetectorElement s2) {
+        return (int)Math.signum(s1.getGeometry().getPosition().z()-s2.getGeometry().getPosition().z());
+      }
+    });
+    int _nDisks = sensors.size();

     _radiusInner = new double[_nDisks];
     _radiusOuter = new double[_nDisks];
     _thickness = new double[_nDisks];
     _z = new double[_nDisks];

-    int id = 0;

+    _nLayers = _nDisks/2;
+    int iDisk = 0;

     for (IDetectorElement del : sensors) {

+      int postfix = (iDisk < _nLayers) ? _nDisks - iDisk - 1 : iDisk - _nLayers;

       IGeometryInfo gInfo = del.getGeometry();
       Tube solid = (Tube) gInfo.getLogicalVolume().getSolid();

-      _radiusInner[id] = solid.getInnerRadius();
-      _radiusOuter[id] = solid.getOuterRadius();
-      _thickness[id] = solid.getZHalfLength()*2.;
-      _z[id] = gInfo.getPosition().z();
-      id++;

+      _radiusInner[postfix] = solid.getInnerRadius();
+      _radiusOuter[postfix] = solid.getOuterRadius();
+      _thickness[postfix] = solid.getZHalfLength()*2.;
+      _z[postfix] = gInfo.getPosition().z();
+      iDisk++;

     }
   }
   
 // -- Setters :  ---------------------------------------------------------------

-  /** Set strip width. Default is 5 micron. */

+  /**
+   * Set strip width.
+   * Default is 25 micron.
+   */

   public void setStripWidth(double pitch) {
     _stripWidth = pitch;
   }

-  /** Set strip length. Default is 10 cm. */

+  /**
+   * Set strip length. 
+   * Default is 10 cm for "TrackerEndcap", 25 micron for "VertexEndcap" and "TrackerForward".
+   */

   public void setStripLength(double length) {
     _stripLength = length;
   }

-  /** Set stereo angle. Default is 90 degrees */

+  /**
+   * Set stereo angle.
+   * Default is 90 degrees for "TrackerEndcap", 0 for "VertexEndcap" and "TrackerForward".
+   */

   public void setStereoAngle(double angle) {
     _rot1 = new Rotation3D(Axis.Z, angle);
   }

@@ -89,14 +115,10 @@

 // -----------------------------------------------------------------------------
 
   /**

-   * Returns sensor ID postfix corresponding to the given position. This postfix must 
-   * be positive  integer, unique within the part of the detector handled by this 
-   * <tt>Segmenter</tt>  object. The final Sensor ID will be constructed taking into
-   * account the prefix set through a call to {@link #setPrefix()} method.

+   * Returns sensor ID postfix corresponding to the given position.

    */
   public int getSensorPostfix(SimTrackerHit hit) {

-    int offSet = hit.getPoint()[2] > 0. ? 0 : _nDisks;
-    return hit.getLayer() + 1 + offSet;

+    return hit.getPoint()[2] > 0. ? hit.getLayer() : _nLayers + hit.getLayer();

   }
   
   /**

@@ -104,15 +126,14 @@

    * {@link #getSensorPostfix(SimTrackerHit)} method of this <tt>Segmenter</tt>object.
    */
   public int getMaxPostfix() {

-    return _nDisks * 2;

+    return _nLayers * 2;

   }
 
   /** Creates a {@link Sensor} object given the ID. */
   public Sensor makeSensor(int postfix) {

-    int id = postfix - 1;
-    SensorType type = new Ring(_radiusInner[id], _radiusOuter[id], _thickness[id], _stripLength, _stripWidth);
-    Hep3Vector trans = new BasicHep3Vector(0.,0.,_z[id]);
-    Rotation3D rot = ((postfix % _nDisks) % 2 == 1) ? _rot0 : _rot1 ;

+    SensorType type = new Ring(_radiusInner[postfix], _radiusOuter[postfix], _thickness[postfix], _stripLength, _stripWidth);
+    Hep3Vector trans = new BasicHep3Vector(0.,0.,_z[postfix]);
+    Rotation3D rot = ((postfix % _nLayers) % 2 == 0) ? _rot0 : _rot1 ;

     return new Sensor(type, trans, rot);
   }

@@ -120,7 +141,7 @@

   
   String _sdName;

-  private int _nDisks;

+  private int _nLayers;

   private double[] _radiusInner;
   private double[] _radiusOuter;
   private double[] _thickness;

Commit in `lcsim/src/org/lcsim/contrib/onoprien/tracking/geom` on MAIN
`SegmentationManager.java`	+2	-1	1.3 -> 1.4
`Segmenter.java`	+34	-10	1.2 -> 1.3
`segmentation/CylindricalBarrelSegmenter.java`	+25	-12	1.4 -> 1.5
`/RingSegmenter.java`	+49	-28	1.4 -> 1.5
	+110	-51