LISTSERV mailing list manager LISTSERV 16.5

Help for HPS-SVN Archives


HPS-SVN Archives

HPS-SVN Archives


HPS-SVN@LISTSERV.SLAC.STANFORD.EDU


View:

Message:

[

First

|

Previous

|

Next

|

Last

]

By Topic:

[

First

|

Previous

|

Next

|

Last

]

By Author:

[

First

|

Previous

|

Next

|

Last

]

Font:

Proportional Font

LISTSERV Archives

LISTSERV Archives

HPS-SVN Home

HPS-SVN Home

HPS-SVN  March 2015

HPS-SVN March 2015

Subject:

r2645 - in /java/trunk: steering-files/src/main/resources/org/hps/steering/recon/ tracking/src/main/java/org/hps/recon/tracking/

From:

[log in to unmask]

Reply-To:

Notification of commits to the hps svn repository <[log in to unmask]>

Date:

Wed, 1 Apr 2015 02:29:10 -0000

Content-Type:

text/plain

Parts/Attachments:

Parts/Attachments

text/plain (180 lines)

Author: [log in to unmask]
Date: Tue Mar 31 19:29:03 2015
New Revision: 2645

Log:
define hit time cuts relative to other hits, not any absolute time

Modified:
    java/trunk/steering-files/src/main/resources/org/hps/steering/recon/HPS2014OfflineRecon.lcsim
    java/trunk/steering-files/src/main/resources/org/hps/steering/recon/HPS2014OfflineTruthRecon.lcsim
    java/trunk/tracking/src/main/java/org/hps/recon/tracking/HelicalTrackHitDriver.java
    java/trunk/tracking/src/main/java/org/hps/recon/tracking/HitTimeTrackCheck.java

Modified: java/trunk/steering-files/src/main/resources/org/hps/steering/recon/HPS2014OfflineRecon.lcsim
 =============================================================================
--- java/trunk/steering-files/src/main/resources/org/hps/steering/recon/HPS2014OfflineRecon.lcsim	(original)
+++ java/trunk/steering-files/src/main/resources/org/hps/steering/recon/HPS2014OfflineRecon.lcsim	Tue Mar 31 19:29:03 2015
@@ -37,13 +37,15 @@
         </driver>
         <driver name="HelicalTrackHitDriver" type="org.hps.recon.tracking.HelicalTrackHitDriver">
             <debug>false</debug>
-            <clusterTimeCut>12.0</clusterTimeCut>
+            <!--<clusterTimeCut>12.0</clusterTimeCut>-->
+            <maxDt>16.0</maxDt>
             <maxSeperation>20.0</maxSeperation>
             <tolerance>1.0</tolerance>
         </driver>
         <driver name="TrackerReconDriver" type="org.hps.recon.tracking.TrackerReconDriver">
             <strategyResource>HPS-Full.xml</strategyResource>
             <debug>false</debug>
+            <rmsTimeCut>8.0</rmsTimeCut>
         </driver>        
         <driver name="GBLDriver" type="org.hps.recon.tracking.gbl.GBLOutputDriver">
             <debug>0</debug>

Modified: java/trunk/steering-files/src/main/resources/org/hps/steering/recon/HPS2014OfflineTruthRecon.lcsim
 =============================================================================
--- java/trunk/steering-files/src/main/resources/org/hps/steering/recon/HPS2014OfflineTruthRecon.lcsim	(original)
+++ java/trunk/steering-files/src/main/resources/org/hps/steering/recon/HPS2014OfflineTruthRecon.lcsim	Tue Mar 31 19:29:03 2015
@@ -38,13 +38,15 @@
         </driver>
         <driver name="HelicalTrackHitDriver" type="org.hps.recon.tracking.HelicalTrackHitDriver">
             <debug>false</debug>
-            <clusterTimeCut>12.0</clusterTimeCut>
+            <!--<clusterTimeCut>12.0</clusterTimeCut>-->
+            <maxDt>16.0</maxDt>
             <maxSeperation>20.0</maxSeperation>
             <tolerance>1.0</tolerance>
         </driver>
         <driver name="TrackerReconDriver" type="org.hps.recon.tracking.TrackerReconDriver">
             <strategyResource>HPS-Full.xml</strategyResource>
             <debug>false</debug>
+            <rmsTimeCut>8.0</rmsTimeCut>
         </driver>        
         <driver name="GBLDriver" type="org.hps.recon.tracking.gbl.GBLOutputDriver">
             <debug>0</debug>

Modified: java/trunk/tracking/src/main/java/org/hps/recon/tracking/HelicalTrackHitDriver.java
 =============================================================================
--- java/trunk/tracking/src/main/java/org/hps/recon/tracking/HelicalTrackHitDriver.java	(original)
+++ java/trunk/tracking/src/main/java/org/hps/recon/tracking/HelicalTrackHitDriver.java	Tue Mar 31 19:29:03 2015
@@ -6,6 +6,7 @@
 import hep.physics.vec.VecOp;
 import java.util.ArrayList;
 import java.util.HashMap;
+import java.util.Iterator;
 import java.util.List;
 import java.util.Map;
 import org.hps.recon.tracking.axial.HelicalTrack2DHit;
@@ -49,6 +50,7 @@
     private boolean _debug = false;
     private double _clusterTimeCut = -99; // if negative, don't cut..otherwise,
     // dt cut time in ns
+    private double maxDt = -99; // max time difference between the two hits in a cross
     private String _subdetectorName = "Tracker";
     private final Map<String, String> _stereomap = new HashMap<String, String>();
     private List<SvtStereoLayer> stereoLayers = null;
@@ -87,6 +89,10 @@
      */
     public void setClusterTimeCut(double dtCut) {
         this._clusterTimeCut = dtCut;
+    }
+
+    public void setMaxDt(double maxDt) {
+        this.maxDt = maxDt;
     }
 
     /**
@@ -332,7 +338,12 @@
                 } // Loop over stereo pairs
             }
 
-            for (HelicalTrackCross cross : helicalTrackCrosses) {
+            for (Iterator<HelicalTrackCross> iter = helicalTrackCrosses.listIterator(); iter.hasNext();) {
+                HelicalTrackCross cross = iter.next();
+                if (maxDt > 0 && Math.abs(cross.getStrips().get(0).time() - cross.getStrips().get(1).time()) > maxDt) {
+                    iter.remove();
+                    continue;
+                }
                 if (cross.getMCParticles() != null)
                     for (MCParticle mcp : cross.getMCParticles())
                         mcrelations.add(new MyLCRelation((HelicalTrackHit) cross, mcp));

Modified: java/trunk/tracking/src/main/java/org/hps/recon/tracking/HitTimeTrackCheck.java
 =============================================================================
--- java/trunk/tracking/src/main/java/org/hps/recon/tracking/HitTimeTrackCheck.java	(original)
+++ java/trunk/tracking/src/main/java/org/hps/recon/tracking/HitTimeTrackCheck.java	Tue Mar 31 19:29:03 2015
@@ -16,7 +16,8 @@
 public class HitTimeTrackCheck implements TrackCheck {
 
     private final double rmsTimeCut;
-    private final int minTrackHits = 10;
+//    private final int minTrackHits = 10;
+    private final int minTrackHits = 6;
     private int seedsChecked = 0;
     private int seedsPassed = 0;
     private int tracksChecked = 0;
@@ -29,18 +30,29 @@
     @Override
     public boolean checkSeed(SeedCandidate candidate) {
 //        System.out.format("seed with %d hits\n", candidate.getHits().size());
-        seedsChecked++;
         int nStrips = 0;
+        double meanTime = 0;
+        for (HelicalTrackHit hth : candidate.getHits()) {
+            for (HelicalTrackStrip hts : ((HelicalTrackCross) hth).getStrips()) {
+                nStrips++;
+                meanTime += hts.time();
+            }
+        }
+        meanTime /= nStrips;
         double rmsTime = 0;
         for (HelicalTrackHit hth : candidate.getHits()) {
             for (HelicalTrackStrip hts : ((HelicalTrackCross) hth).getStrips()) {
-                nStrips++;
-                rmsTime += hts.time() * hts.time();
+                rmsTime += Math.pow(hts.time() - meanTime, 2);
+//                rmsTime += hts.time() * hts.time();
+//                rmsTime += Math.abs(hts.time());
             }
         }
+//        if (nStrips<6) return true;
+        seedsChecked++;
 //        rmsTime = Math.sqrt(rmsTime / nStrips);
 //        System.out.format("seed RMS %f on %d hits\n",rmsTime,nStrips);
         boolean passCheck = (rmsTime < minTrackHits * rmsTimeCut * rmsTimeCut);
+//        boolean passCheck = (rmsTime < minTrackHits * rmsTimeCut);
         if (passCheck) {
             seedsPassed++;
         }
@@ -56,15 +68,25 @@
 //        System.out.format("track with %d hits\n", track.getTrackerHits().size());
         tracksChecked++;
         int nStrips = 0;
+        double meanTime = 0;
+        for (TrackerHit hit : track.getTrackerHits()) {
+            for (HelicalTrackStrip hts : ((HelicalTrackCross) hit).getStrips()) {
+                nStrips++;
+                meanTime += hts.time();
+            }
+        }
+        meanTime /= nStrips;
         double rmsTime = 0;
         for (TrackerHit hit : track.getTrackerHits()) {
             for (HelicalTrackStrip hts : ((HelicalTrackCross) hit).getStrips()) {
-                nStrips++;
-                rmsTime += hts.time() * hts.time();
+                rmsTime += Math.pow(hts.time() - meanTime, 2);
+//                rmsTime += hts.time() * hts.time();
+//                rmsTime += Math.abs(hts.time());
             }
         }
         rmsTime = Math.sqrt(rmsTime / nStrips);
-        System.out.format("track RMS %f on %d hits\n", rmsTime, nStrips);
+//        rmsTime = rmsTime / nStrips;
+//        System.out.format("track RMS %f on %d hits\n", rmsTime, nStrips);
         boolean passCheck = (rmsTime < rmsTimeCut);
         if (passCheck) {
             tracksPassed++;

Top of Message | Previous Page | Permalink

Advanced Options


Options

Log In

Log In

Get Password

Get Password


Search Archives

Search Archives


Subscribe or Unsubscribe

Subscribe or Unsubscribe


Archives

November 2017
August 2017
July 2017
January 2017
December 2016
November 2016
October 2016
September 2016
August 2016
July 2016
June 2016
May 2016
April 2016
March 2016
February 2016
January 2016
December 2015
November 2015
October 2015
September 2015
August 2015
July 2015
June 2015
May 2015
April 2015
March 2015
February 2015
January 2015
December 2014
November 2014
October 2014
September 2014
August 2014
July 2014
June 2014
May 2014
April 2014
March 2014
February 2014
January 2014
December 2013
November 2013

ATOM RSS1 RSS2



LISTSERV.SLAC.STANFORD.EDU

Secured by F-Secure Anti-Virus CataList Email List Search Powered by the LISTSERV Email List Manager

Privacy Notice, Security Notice and Terms of Use