Ana1921 pass0

detector-data/detectors/HPS_Nominal_3pt7_iter0/SamplingFractions/Ecal.properties

@@ -0,0 +1 @@
+samplingFraction: 1.0

detector-data/detectors/HPS_Nominal_3pt7_iter0/SamplingFractions/Hodoscope.properties

@@ -0,0 +1 @@
+samplingFraction: 1.0

detector-data/detectors/HPS_Nominal_3pt7_iter0/detector.properties

@@ -0,0 +1 @@
+name: HPS_Nominal_3pt7_iter0

detector-data/detectors/HPS_Nominal_iter1/SamplingFractions/Ecal.properties

@@ -0,0 +1 @@
+samplingFraction: 1.0

detector-data/detectors/HPS_Nominal_iter1/SamplingFractions/Hodoscope.properties

@@ -0,0 +1 @@
+samplingFraction: 1.0

detector-data/detectors/HPS_Nominal_iter1/detector.properties

@@ -0,0 +1 @@
+name: HPS_Nominal_iter1

detector-data/detectors/HPS_Nominal_iter2/SamplingFractions/Ecal.properties

@@ -0,0 +1 @@
+samplingFraction: 1.0

detector-data/detectors/HPS_Nominal_iter2/SamplingFractions/Hodoscope.properties

@@ -0,0 +1 @@
+samplingFraction: 1.0

detector-data/detectors/HPS_Nominal_iter2/detector.properties

@@ -0,0 +1 @@
+name: HPS_Nominal_iter2

detector-data/detectors/HPS_PhysicsRun2019-v2-FEE-Pass2_iter17/SamplingFractions/Ecal.properties

@@ -0,0 +1 @@
+samplingFraction: 1.0

detector-data/detectors/HPS_PhysicsRun2019-v2-FEE-Pass2_iter17/SamplingFractions/Hodoscope.properties

@@ -0,0 +1 @@
+samplingFraction: 1.0

detector-data/detectors/HPS_PhysicsRun2019-v2-FEE-Pass2_iter17/detector.properties

@@ -0,0 +1 @@
+name: HPS_PhysicsRun2019-v2-FEE-Pass2_iter17

detector-data/detectors/HPS_PhysicsRun2019-v2-FEE-Pass2_iter17_fieldOff/SamplingFractions/Ecal.properties

@@ -0,0 +1 @@
+samplingFraction: 1.0

detector-data/detectors/HPS_PhysicsRun2019-v2-FEE-Pass2_iter17_fieldOff/SamplingFractions/Hodoscope.properties

@@ -0,0 +1 @@
+samplingFraction: 1.0

detector-data/detectors/HPS_PhysicsRun2019-v2-FEE-Pass2_iter17_fieldOff/detector.properties

@@ -0,0 +1,2 @@
+name: HPS_PhysicsRun2019-v2-FEE-Pass2_iter17_fieldOff
+

detector-data/detectors/HPS_PhysicsRun2019_Pass2/SamplingFractions/Ecal.properties

@@ -0,0 +1 @@
+samplingFraction: 1.0

detector-data/detectors/HPS_PhysicsRun2019_Pass2/SamplingFractions/Hodoscope.properties

@@ -0,0 +1 @@
+samplingFraction: 1.0

detector-data/detectors/HPS_PhysicsRun2019_Pass2/detector.properties

@@ -0,0 +1 @@
+name: HPS_PhysicsRun2019_Pass2

detector-data/detectors/HPS_Run2021Pass0_v1/SamplingFractions/Ecal.properties

@@ -0,0 +1 @@
+samplingFraction: 1.0

detector-data/detectors/HPS_Run2021Pass0_v1/SamplingFractions/Hodoscope.properties

@@ -0,0 +1 @@
+samplingFraction: 1.0

detector-data/detectors/HPS_Run2021Pass0_v1/detector.properties

@@ -0,0 +1 @@
+name: HPS_Run2021Pass0_v1

online-monitoring/src/main/java/org/hps/monitoring/drivers/trackrecon/SVTOpeningAlignment.java

@@ -13,6 +13,7 @@
 import hep.aida.IPlotterStyle;
 
 import java.io.IOException;
+import static java.lang.Math.abs;
 import java.util.ArrayList;
 import java.util.List;
 import java.util.logging.Level;
@@ -108,6 +109,8 @@ public class SVTOpeningAlignment extends Driver {
     String fullTrackCollectionName = "s234_c5_e167";
     double targetPosition = -5.0; //mm
 
+    boolean _writeRunAndEventNumbers = true;
+
     public SVTOpeningAlignment() {
     }
 
@@ -299,24 +302,30 @@ protected void detectorChanged(Detector detector) {
     @Override
     public void process(EventHeader event) {
         aida.tree().cd("/");
-        if (!event.hasCollection(HelicalTrackHit.class, helicalTrackHitCollectionName))
+        if (!event.hasCollection(HelicalTrackHit.class, helicalTrackHitCollectionName)) {
             return;
+        }
 
-        if (!event.hasCollection(Track.class, l0to3CollectionName))
+        if (!event.hasCollection(Track.class, l0to3CollectionName)) {
             return;
+        }
 
-        if (!event.hasCollection(Track.class, l4to6CollectionName))
+        if (!event.hasCollection(Track.class, l4to6CollectionName)) {
             return;
+        }
 
-        if (matchFullTracks)
-            if (!event.hasCollection(Track.class, fullTrackCollectionName))
+        if (matchFullTracks) {
+            if (!event.hasCollection(Track.class, fullTrackCollectionName)) {
                 return;
+            }
+        }
 
         List<Track> l0to3tracks = event.get(Track.class, l0to3CollectionName);
         List<Track> l4to6tracks = event.get(Track.class, l4to6CollectionName);
         List<Track> fulltracks = new ArrayList<Track>();
-        if (matchFullTracks)
+        if (matchFullTracks) {
             fulltracks = event.get(Track.class, fullTrackCollectionName);
+        }
 
         List<Track> l0to3tracksTop = splitTrackList(l0to3tracks, true);
         List<Track> l0to3tracksBot = splitTrackList(l0to3tracks, false);
@@ -365,8 +374,9 @@ public void process(EventHeader event) {
             for (Track trk03 : l0to3tracksTop) {
                 if (matchFullTracks) {
                     matchedTrack = checkFullTrack(fulltracks, trk03, trk46);
-                    if (matchedTrack == null)
+                    if (matchedTrack == null) {
                         continue;
+                    }
                 }
                 TrackState ts03 = trk03.getTrackStates().get(0);
                 double x0L03 = TrackUtils.getX0(ts03);
@@ -404,7 +414,10 @@ public void process(EventHeader event) {
                 yAtHingeL03VsL46Top.fill(yAtHingeL46, yAtHingeL03);
                 delYAtHingeVsL46SlopeTop.fill(deltaYAtHinge, slL46);
                 delYAtHingeVsL03SlopeTop.fill(deltaYAtHinge, slL03);
-
+                if (_writeRunAndEventNumbers) {
+                    System.out.println(event.getRunNumber() + " " + event.getEventNumber() + " t");
+                }
+                analyzeMatchedTrack("top", matchedTrack);
             }
         }
         fitAndPutParameters(deld0Top, fd0Top);
@@ -418,8 +431,9 @@ public void process(EventHeader event) {
             for (Track trk03 : l0to3tracksBot) {
                 if (matchFullTracks) {
                     matchedTrack = checkFullTrack(fulltracks, trk03, trk46);
-                    if (matchedTrack == null)
+                    if (matchedTrack == null) {
                         continue;
+                    }
                 }
                 TrackState ts03 = trk03.getTrackStates().get(0);
                 double x0L03 = TrackUtils.getX0(ts03);
@@ -457,6 +471,10 @@ public void process(EventHeader event) {
                 yAtHingeL03VsL46Bot.fill(yAtHingeL46, yAtHingeL03);
                 delYAtHingeVsL46SlopeBot.fill(deltaYAtHinge, slL46);
                 delYAtHingeVsL03SlopeBot.fill(deltaYAtHinge, slL03);
+                if (_writeRunAndEventNumbers) {
+                    System.out.println(event.getRunNumber() + " " + event.getEventNumber() + " b");
+                }
+                analyzeMatchedTrack("bottom", matchedTrack);
             }
         }
 
@@ -472,13 +490,14 @@ public void process(EventHeader event) {
 
     @Override
     public void endOfData() {
-        if (outputPlots != null)
+        if (outputPlots != null) {
             try {
                 plotterTop.writeToFile(outputPlots + "-deltasTop.gif");
                 plotterBot.writeToFile(outputPlots + "-deltasBottom.gif");
             } catch (IOException ex) {
                 Logger.getLogger(TrackingReconPlots.class.getName()).log(Level.SEVERE, null, ex);
             }
+        }
     }
 
     private List<Track> splitTrackList(List<Track> trks, boolean doTop) {
@@ -488,15 +507,22 @@ private List<Track> splitTrackList(List<Track> trks, boolean doTop) {
         for (Track trk : trks) {
             isTop = false;
             isBot = false;
-            for (TrackerHit hit : trk.getTrackerHits())
+            for (TrackerHit hit : trk.getTrackerHits()) {
                 if (hit.getPosition()[2] > 0)// remember, non-bend in tracking frame is z-direction
+                {
                     isTop = true;
-                else
+                } else {
                     isBot = true;
-            if (isTop == true && isBot != true && doTop == true)  // if we want top tracks and all hits are in top
+                }
+            }
+            if (isTop == true && isBot != true && doTop == true) // if we want top tracks and all hits are in top
+            {
                 tracksHalf.add(trk);
+            }
             if (isBot == true && isTop != true && doTop == false) // if we want bottom tracks and all hits are in bottom
+            {
                 tracksHalf.add(trk);
+            }
         }
         return tracksHalf;
     }
@@ -507,10 +533,30 @@ private Track checkFullTrack(List<Track> fullTracks, Track t03, Track t46) {
 
         for (Track fullTr : fullTracks) {
             List<TrackerHit> trkHitsFull = fullTr.getTrackerHits();
-            if (trkHitsFull.containsAll(trkHitsL03) && trkHitsFull.containsAll(trkHitsL46))
+            if (trkHitsFull.containsAll(trkHitsL03) && trkHitsFull.containsAll(trkHitsL46)) {
                 return fullTr;
+            }
         }
 
         return null;
     }
+
+    private void analyzeMatchedTrack(String torb, Track t) {
+        aida.tree().mkdirs("track quality");
+        aida.tree().cd("track quality");
+        String pdg = t.getCharge() < 0 ? " electron " : " positron ";
+        aida.histogram1D("chisq per ndf " + pdg + torb, 100, 0., 15.).fill(t.getChi2() / t.getNDF());
+        aida.histogram1D("number of hits " + pdg + torb, 20, -0.5, 19.5).fill(t.getTrackerHits().size());
+        aida.histogram1D("charge " + pdg + torb, 3, -1.5, 1.5).fill(t.getCharge());
+        TrackState ts = t.getTrackStates().get(0);
+        aida.histogram1D("tanLambda " + pdg + torb, 50, 0, 0.06).fill(abs(ts.getTanLambda()));
+        aida.histogram1D("z0 " + pdg + torb, 50, -1.0, 1.0).fill(ts.getZ0());
+        aida.histogram1D("d0 " + pdg + torb, 50, -3.0, 3.0).fill(ts.getD0());
+
+        aida.tree().cd("..");
+    }
+
+    public void setWriteRunAndEventNumbers(boolean b) {
+        _writeRunAndEventNumbers = b;
+    }
 }