Author: [log in to unmask] Date: Thu Sep 10 11:00:28 2015 New Revision: 3578 Log: undo clusterer changes for Holly Modified: java/trunk/ecal-recon/src/main/java/org/hps/recon/ecal/cluster/ReconClusterer.java java/trunk/integration-tests/src/test/java/org/hps/test/it/EcalSimReconTest.java Modified: java/trunk/ecal-recon/src/main/java/org/hps/recon/ecal/cluster/ReconClusterer.java ============================================================================= --- java/trunk/ecal-recon/src/main/java/org/hps/recon/ecal/cluster/ReconClusterer.java (original) +++ java/trunk/ecal-recon/src/main/java/org/hps/recon/ecal/cluster/ReconClusterer.java Thu Sep 10 11:00:28 2015 @@ -64,7 +64,7 @@ double minTime = 0.0; // Maximum time cut window range. Units in ns. - double timeWindow = 6.0; + double timeWindow = 8.0; // Make a map for quick calculation of the x-y position of crystal face Map<Point, double[]> correctedPositionMap = new HashMap<Point, double[]>(); @@ -73,7 +73,7 @@ ReconClusterer() { super(new String[] { "hitEnergyThreshold", "seedEnergyThreshold", "clusterEnergyThreshold", "minTime", "timeWindow" }, - new double[] { 0.0075, 0.05, 0.1, 0.0, 6.0 }); + new double[] { 0.0075, 0.05, 0.1, 0.0, 8.0 }); } void setUseTimeCut(boolean useTimeCut) { @@ -383,16 +383,10 @@ // Consider time cut-is this hit in same time window as seed? if (useTimeCut){ - if(ihit.getCorrectedEnergy() < 0.1 && (Math.abs(ihit.getTime() - iseed.getTime()) < timeWindow)) + if(Math.abs(ihit.getTime() - iseed.getTime()) < timeWindow) { icluster.addHit(ihit); } - - else if (ihit.getCorrectedEnergy() > 0.1 && (Math.abs(ihit.getTime() - iseed.getTime()) < 2.0)) - { - icluster.addHit(ihit); - } - } // end of using time cut else {icluster.addHit(ihit);} } @@ -406,29 +400,15 @@ boolean inTimeWithA = false; boolean inTimeWithB = false; // In time window with seedA? - if (commHit.getKey().getCorrectedEnergy() < 0.1){ - if (Math.abs(commHit.getKey().getTime() - seedA.getTime()) < timeWindow){ - inTimeWithA = true; - } + if (Math.abs(commHit.getKey().getTime() - seedA.getTime()) < timeWindow){ + inTimeWithA = true; + } - // In time window with seedB? - if (Math.abs(commHit.getKey().getTime() - seedB.getTime()) < timeWindow){ - inTimeWithB = true; - } - } - - else if (commHit.getKey().getCorrectedEnergy() > 0.1) - { - if (Math.abs(commHit.getKey().getTime() - seedA.getTime()) < 2.0){ - inTimeWithA = true; - } - - // In time window with seedB? - if (Math.abs(commHit.getKey().getTime() - seedB.getTime()) < 2.0){ - inTimeWithB = true; - } - } - + // In time window with seedB? + if (Math.abs(commHit.getKey().getTime() - seedB.getTime()) < timeWindow){ + inTimeWithB = true; + } + double eclusterA = seedToCluster.get(seedA).getEnergy(); double eclusterB = seedToCluster.get(seedB).getEnergy(); double fractionA = eclusterA / (eclusterA + eclusterB); Modified: java/trunk/integration-tests/src/test/java/org/hps/test/it/EcalSimReconTest.java ============================================================================= --- java/trunk/integration-tests/src/test/java/org/hps/test/it/EcalSimReconTest.java (original) +++ java/trunk/integration-tests/src/test/java/org/hps/test/it/EcalSimReconTest.java Thu Sep 10 11:00:28 2015 @@ -104,7 +104,7 @@ // Check high cluster time mean. TestCase.assertEquals("High cluster mean time does not match.", CLUS_MEAN_T, - clusTimeH1D.mean(), 0.3); + clusTimeH1D.mean(), 0.03); // Check mean number of clusters per event. TestCase.assertEquals("Mean number of clusters per event does not match.", CLUS_COUNT_MEAN,