for (HPSFittedRawTrackerHit hrth : fittedrawHits) { SiSensor sensor = (SiSensor) hrth.getRawTrackerHit().getDetectorElement(); int layer = hrth.getRawTrackerHit().getLayerNumber();

// pedestalShifts[module][layer] = aida.histogram2D(HPSSVTDAQMaps.sensorArray[module][layer].getName() + " Pedestal Shifts", 640, -0.5, 639.5, 100, -500, 500); plotter.region(region).plot(pedestalShifts[module][layer]);

plotter4.region(region).plot(t0a[module][layer]); // plotter3.region(region).style().zAxisStyle().setParameter("labelVertical", "false"); // for (String str : plotter3.region(region).style().zAxisStyle().availableParameterOptions("labelVertical"))

double x = htc.getPosition()[0]; double y = htc.getPosition()[1]; SiSensor sensor = ((SiSensor) ((RawTrackerHit) htc.getRawHits().get(0)).getDetectorElement());

layersTop[htc.Layer() - 1]++; Hep3Vector sensorPos = ((SiSensor) ((RawTrackerHit) htc.getRawHits().get(0)).getDetectorElement()).getGeometry().getPosition(); if (htc.Layer() == 1) {

ChargeCarrier carrier = ChargeCarrier.getCarrier(_sid_helper.getSideValue(id)); SiSensorElectrodes electrodes = ((SiSensor) hit.getDetectorElement()).getReadoutElectrodes(carrier);

hps-java/src/main/java/org/lcsim/hps/recon/tracking

SvtUtils.java added at 1.1

diff -N SvtUtils.java
--- /dev/null	1 Jan 1970 00:00:00 -0000
+++ SvtUtils.java	24 May 2012 17:08:18 -0000	1.1
@@ -0,0 +1,220 @@

+package org.lcsim.hps.recon.tracking;
+
+//--- java ---//
+import java.util.HashMap;
+import java.util.HashSet;
+import java.util.Map;
+import java.util.Set;
+
+//--- org.lcsim ---//
+import org.lcsim.conditions.ConditionsManager;
+import org.lcsim.conditions.ConditionsSet;
+import org.lcsim.detector.IDetectorElement;
+import org.lcsim.detector.identifier.IIdentifier;
+import org.lcsim.detector.tracker.silicon.SiSensor;
+import org.lcsim.detector.tracker.silicon.SiTrackerIdentifierHelper;
+import org.lcsim.geometry.Detector;
+
+//--- hps-java ---//
+import org.lcsim.hps.util.Pair;
+
+/**
+ * 
+ * @author Omar Moreno
+ * @version $Id: SvtUtils.java,v 1.1 2012/05/24 17:08:18 omoreno Exp $
+ */
+public class SvtUtils {
+		
+	// Use eager Singleton initialization since an instance will always be needed
+	private static final SvtUtils INSTANCE = new SvtUtils();
+	
+    // Set of sensors
+    private Set<SiSensor> sensors = new HashSet<SiSensor>();
+    
+    // Map from Sensor to Hybrid/FPGA pair
+    private Map<SiSensor /* Sensor */, Pair<Integer /* Hybrid */, Integer /* FPGA */>> sensorToDAQPair
+    	= new HashMap<SiSensor, Pair<Integer, Integer>>();
+    
+    // Map from Hybrid/FPGA pair
+    private Map<Pair<Integer /* Hybrid */, Integer /* FPGA */>, SiSensor /* Sensor*/> daqPairToSensor
+   		= new HashMap<Pair<Integer, Integer>, SiSensor>();
+    
+    // Map from SVT top layer to Hybrid/FPGA pair
+    private Map<Integer /* Layer # */, Pair<Integer /* Hybrid */, Integer /* FPGA */>> topLayerToDAQ
+    	= new HashMap<Integer, Pair<Integer, Integer>>();
+    
+    // Map from SVT bottom layer to Hybrid/FPGA pair
+    private Map<Integer /* Layer # */, Pair<Integer /* Hybrid */, Integer /* FPGA */>> bottomLayerToDAQ
+    	= new HashMap<Integer, Pair<Integer, Integer>>();
+    
+    // Map sensor to an SVT top layer
+    private Map<SiSensor /* Sensor */, Integer /* Layer */> sensorToTopLayer
+    	= new HashMap<SiSensor, Integer>();
+    
+    // Map sensor to an SVT bottom layer
+    private Map<SiSensor /* Sensor */, Integer /* Layer */> sensorToBottomLayer
+    	= new HashMap<SiSensor, Integer>();
+
+    private SiSensor[][] sensorArray = new SiSensor[2][10];
+    
+	boolean debug = true;
+	String subdetectorName = "Tracker";
+	ConditionsManager mgr = ConditionsManager.defaultInstance();
+    
+    private SvtUtils(){}
+    
+    /**
+     * 
+     */
+    public static SvtUtils getInstance(){
+    	return INSTANCE;
+    }
+    
+    /**
+     * 
+     * @param sensor
+     * @return
+     */
+    public int getLayerNumber(SiSensor sensor){
+    	if(sensorToTopLayer.containsKey(sensor)) return sensorToTopLayer.get(sensor);
+		else if(sensorToBottomLayer.containsKey(sensor)) return sensorToBottomLayer.get(sensor);
+		else throw new RuntimeException("There is no layer associated with sensor " + sensor.getName());
+    }
+    
+    /**
+     * 
+     */
+    public boolean isTopLayer(SiSensor sensor){
+    	if(sensorToTopLayer.containsKey(sensor)) return true;
+		else if(sensorToBottomLayer.containsKey(sensor)) return false;
+		else throw new RuntimeException("There is no layer associated with sensor " + sensor.getName());
+    }
+    
+    /**
+     * 
+     */
+    public boolean isAxial(SiSensor sensor){
+    	if(this.isTopLayer(sensor) && this.getLayerNumber(sensor)%2 == 1) return true;
+    	else if(!this.isTopLayer(sensor) && this.getLayerNumber(sensor)%2 == 0) return true;
+    	return false;
+    }
+    
+    /**
+     * 
+     */
+    public SiSensor getSensor(int module, int layer){
+    	if(module < 0 || module > 1)
+    		throw new RuntimeException("Module number " + module + " is out of range!");
+    	else if(layer < 0 || layer > 10)
+    		throw new RuntimeException("Layer number " + layer + " is out of range!");
+    	return sensorArray[module][layer];
+    }
+    
+    /**
+     * 
+     */
+    public SiSensor getSensor(Pair<Integer, Integer> daqPair){
+    	return this.daqPairToSensor.get(daqPair);
+    }
+    
+    /**
+     * 
+     * @return
+     */
+    public Set<SiSensor> getSensors(){
+    	return sensors;
+    }
+    
+    /**
+     * 
+     */
+    public int getFPGA(SiSensor sensor){
+    	return this.sensorToDAQPair.get(sensor).getSecondElement();
+    }
+    
+    /**
+     * 
+     */
+    public int getHybrid(SiSensor sensor){
+    	return this.sensorToDAQPair.get(sensor).getFirstElement();
+    }
+    
+    /**
+     * 
+     */
+    public Pair<Integer, Integer> getDaqPair(SiSensor sensor){
+    	return this.sensorToDAQPair.get(sensor);
+    }
+    
+    /**
+     * 
+     */
+    protected void setupDaqMaps(Detector detector){
+    	
+        if (subdetectorName == null)
+            throw new RuntimeException("The subdetector name was not set!");
+        
+        // Get the SVT DAQ Maps
+        ConditionsSet topSvtDaqMap = mgr.getConditions("daqmap/SVTDAQMapTop");
+        ConditionsSet bottomSvtDaqMap = mgr.getConditions("daqmap/SVTDAQMapBottom");
+        String[] daqPair;
+        
+        // Loop over the top layers and fill the DAQ map
+        for (Object topLayer : topSvtDaqMap.keySet()) {
+            daqPair = topSvtDaqMap.getString(topLayer.toString()).split(" ");
+            if (debug)
+                System.out.println("Layer: " + topLayer.toString() + " FPGA: " + daqPair[1] +  " Hybrid: " + daqPair[0] );
+            topLayerToDAQ.put(Integer.parseInt(topLayer.toString()), new Pair<Integer, Integer>(Integer.parseInt(daqPair[0]), Integer.parseInt(daqPair[1])));
+        }
+                
+    	// Loop over the bottom layers and fill the DAQ Map
+        for (Object bottomLayer : bottomSvtDaqMap.keySet()) {
+            daqPair = bottomSvtDaqMap.getString(bottomLayer.toString()).split(" ");
+            if (debug)
+                System.out.println("Layer: " + bottomLayer.toString() + " FPGA: " + daqPair[1] + " Hybrid: " + daqPair[0]);
+            bottomLayerToDAQ.put(Integer.parseInt(bottomLayer.toString()), new Pair<Integer, Integer>(Integer.parseInt(daqPair[0]), Integer.parseInt(daqPair[1])));
+        }
+        
+        // Get the SVT
+        IDetectorElement detectorElement = detector.getDetectorElement().findDetectorElement(subdetectorName);
+
+        // Get all SVT sensors
+        sensors.addAll(detectorElement.findDescendants(SiSensor.class));    
+    
+        IIdentifier sensorIdent;
+        SiTrackerIdentifierHelper sensorHelper;
+        
+        // Create the maps
+        for(SiSensor sensor : sensors){
+        	
+        	// Get the sensor identifier
+        	sensorIdent = sensor.getIdentifier();
+        	
+            // Get the sensor identifier helper
+            sensorHelper = (SiTrackerIdentifierHelper) sensor.getIdentifierHelper();
+            
+            int layerNumber = sensorHelper.getLayerValue(sensorIdent);
+            int moduleNumber = sensorHelper.getModuleValue(sensorIdent);
+        	sensorArray[moduleNumber][layerNumber - 1] = sensor;
+            
+            if (moduleNumber == 0) {
+            	sensorToTopLayer.put(sensor, layerNumber);
+            	sensorToDAQPair.put(sensor, this.topLayerToDAQ.get(layerNumber));
+            	daqPairToSensor.put(topLayerToDAQ.get(layerNumber), sensor);
+            } else if(moduleNumber == 1){
+                sensorToBottomLayer.put(sensor, layerNumber);
+            	sensorToDAQPair.put(sensor, this.bottomLayerToDAQ.get(layerNumber));
+            	daqPairToSensor.put(bottomLayerToDAQ.get(layerNumber), sensor);
+            } else throw new RuntimeException("Invalid Module Number: " + sensorHelper.getModuleValue(sensorIdent));
+        }
+    
+        for (int module = 0; module < 2; module++) {
+            for (int layer = 0; layer < 10; layer++) {
+                if (sensorArray[module][layer] == null) {
+                    throw new RuntimeException("Couldn't find sensor with module " + module + ", layer " + layer + 1);
+                }
+            }
+        }
+    }
+    
+}

+ + for(Map.Entry<String, List<HelicalTrackStrip>> entry : striplistmap.entrySet()){ + String ident = entry.getKey(); + List<HelicalTrackStrip> hts = entry.getValue(); + for(HelicalTrackStrip hstrip : hts ){ + } + System.out.println(this.getClass().getSimpleName() + ": Total number of HelicalTrackStrips: " + ident + ": " + hts.size()); + }

- Pair<Integer, Integer> daqPair = new Pair<Integer, Integer>(svtData.getHybridNumber(), svtData.getFPGAAddress()); - SiSensor sensor = HPSSVTDAQMaps.daqPairToSensor.get(daqPair); - if (sensor == null) {

+ if(debug) + System.out.println(this.getClass().getSimpleName() + ": FPGA: " + svtData.getFPGAAddress() + " Hybrid: " + svtData.getHybridNumber()); + + Pair<Integer, Integer> daqPair = new Pair<Integer, Integer>(svtData.getHybridNumber(), svtData.getFPGAAddress()); + SiSensor sensor = SvtUtils.getInstance().getSensor(daqPair); + + if (sensor == null) {

+ if(debug){ + if(SvtUtils.getInstance().isTopLayer(sensor)){ + System.out.println(this.getClass().getSimpleName() + ": Top Layer Number: " + SvtUtils.getInstance().getLayerNumber(sensor)); + } else { + System.out.println(this.getClass().getSimpleName() + ": Bottom Layer Number: " + SvtUtils.getInstance().getLayerNumber(sensor)); + } + } +

int sensorChannel = svtData.getSensorChannel(); if (sensorChannel >= HPSSVTConstants.TOTAL_STRIPS_PER_SENSOR) { continue; //strip 639 exists as an APV channel but not as a sensor strip; drop it

+ + if(debug) + System.out.println(this.getClass().getSimpleName() + "Cell ID Module " + myHit.getIdentifierFieldValue("module") + "Cell ID Layer: " + myHit.getIdentifierFieldValue("layer") + "Strip Number: " + myHit.getIdentifierFieldValue("strip"));

aida.histogram2D("Track Position at Ecal: Two Tracks").fill(positionEcal.y(), positionEcal.z()); aida.histogram2D("Track Position at Harp: Two Tracks").fill(positionConverter.y(), positionConverter.z()); aida.histogram1D("Px: Two Tracks").fill(track.getPX());

+ + double tracksRatio = ntracks/nevents; + double tracksTopRatio = ntracksTop/nevents; + double tracksBottomRatio = ntracksBottom/nevents; + double twoTrackRatio = nTwoTracks/nevents; + System.out.println("Number of tracks per event: " + tracksRatio); + System.out.println("Number of top tracks per event: " + tracksTopRatio); + System.out.println("Number of bottom tracks per event: " + tracksBottomRatio); + System.out.println("Number of two track events: " + twoTrackRatio);

Commit in `hps-java/src/main/java/org/lcsim/hps` on MAIN
`monitoring/HPSCalibrationDriver.java`	+2	-2	1.1 -> 1.2
`monitoring/svt/SVTSimpleEventDisplay.java`	+4	-3	1.6 -> 1.7
`/SVTHitReconstructionPlots.java`	+2	-1	1.15 -> 1.16
`/SVTMonitoringPlots.java`	+5	-4	1.4 -> 1.5
`/TrackingReconstructionPlots.java`	+1	-1	1.8 -> 1.9
`/SVTEventInfo.java`	+2	-1	1.6 -> 1.7
`recon/tracking/SvtUtils.java`	+220		added 1.1
`/HPSMakeSVTCellID.java`	+2	-2	1.4 -> 1.5
`/HPSDataProcessingModule.java`	+4	-4	1.13 -> 1.14
`/HPSHelicalTrackHitDriver.java`	+10	-1	1.6 -> 1.7
`/HPSSVTSensorSetup.java`	+6	-1	1.7 -> 1.8
`/HPSSVTCalibrationConstants.java`	+3	-3	1.11 -> 1.12
`/SVTDataToRawTrackerHitDriver.java`	+19	-4	1.9 -> 1.10
`users/omoreno/SvtTrackAnalysis.java`	+22	-3	1.2 -> 1.3
`recon/tracking/apv25/HPSSiSensorReadout.java`	+4	-3	1.8 -> 1.9
	+306	-33