Author: [log in to unmask]
Date: Mon Mar 21 17:40:38 2016
New Revision: 4316
Log:
Undo r4306 to fix conditions db issue.
Modified:
java/trunk/ecal-recon/src/main/java/org/hps/recon/ecal/EcalTimeWalk.java
Modified: java/trunk/ecal-recon/src/main/java/org/hps/recon/ecal/EcalTimeWalk.java
=============================================================================
--- java/trunk/ecal-recon/src/main/java/org/hps/recon/ecal/EcalTimeWalk.java (original)
+++ java/trunk/ecal-recon/src/main/java/org/hps/recon/ecal/EcalTimeWalk.java Mon Mar 21 17:40:38 2016
@@ -1,6 +1,4 @@
package org.hps.recon.ecal;
-
-import org.hps.conditions.database.DatabaseConditionsManager;
public final class EcalTimeWalk {
@@ -8,76 +6,62 @@
// Put parameters in the conditions database once functional form is fixed.
// From a fit of Pass0 Run 3261 (Mode-3 FADC):
- // (p0+p1*e+e^2*p2) * exp(-(p3*e+p4*e^2+p5*e^4))
- private static final double[] pars = { 3.64218e+01, -4.60756e+02,
- 9.18743e+03, 3.73873e+01, -6.57130e+01, 1.07182e+02 };
+ // (p0+p1*e+e^2*p2) * exp(-(p3*e+p4*e^2+p5*e^4))
+ private static final double[] pars = {
+ 3.64218e+01,
+ -4.60756e+02,
+ 9.18743e+03,
+ 3.73873e+01,
+ -6.57130e+01,
+ 1.07182e+02
+ };
private EcalTimeWalk() {
}
/**
* Perform Time Walk Correction
- *
- * @param time
- * - FADC Mode-3 Hit time (ns)
- * @param energy
- * - Pulse energy (GeV)
+ * @param time - FADC Mode-3 Hit time (ns)
+ * @param energy - Pulse energy (GeV)
* @return corrected time (ns)
*/
public static final double correctTimeWalk(double time, double energy) {
- final double poly1 = pars[0] + pars[1] * energy + pars[2] * energy
- * energy;
- final double poly2 = pars[3] * energy + pars[4] * energy * energy
- + pars[5] * Math.pow(energy, 4);
+ final double poly1 = pars[0] +
+ pars[1] * energy +
+ pars[2] * energy * energy;
+ final double poly2 = pars[3] * energy +
+ pars[4] * energy * energy +
+ pars[5] * Math.pow(energy, 4);
return time - poly1 * Math.exp(-poly2);
}
-
+
/*
* Time walk parameters for pulse fitting
*/
-
private static final double[] par = {
- DatabaseConditionsManager
- .getInstance()
- .getCachedConditions(
- org.hps.conditions.ecal.EcalTimeWalk.class,
- "ecal_time_walk").getCachedData().getP0(),
- DatabaseConditionsManager
- .getInstance()
- .getCachedConditions(
- org.hps.conditions.ecal.EcalTimeWalk.class,
- "ecal_time_walk").getCachedData().getP1(),
- DatabaseConditionsManager
- .getInstance()
- .getCachedConditions(
- org.hps.conditions.ecal.EcalTimeWalk.class,
- "ecal_time_walk").getCachedData().getP2(),
- DatabaseConditionsManager
- .getInstance()
- .getCachedConditions(
- org.hps.conditions.ecal.EcalTimeWalk.class,
- "ecal_time_walk").getCachedData().getP3(),
- DatabaseConditionsManager
- .getInstance()
- .getCachedConditions(
- org.hps.conditions.ecal.EcalTimeWalk.class,
- "ecal_time_walk").getCachedData().getP4() };
-
+ 0.9509,
+ -33.21,
+ 0.2614,
+ -0.9128,
+ 0.6251
+ };
+
/**
* Perform Time Walk Correction for Mode 1 hits using pulse fitting
- *
- * @param time
- * - FADC Mode 1 hit time from pulse fitting (ns)
- * @param energy
- * - Pulse energy from pulse fitting (GeV)
+ * @param time - FADC Mode 1 hit time from pulse fitting (ns)
+ * @param energy - Pulse energy from pulse fitting (GeV)
* @return corrected time (ns)
*/
- public static final double correctTimeWalkPulseFitting(double time,
- double energy) {
- final double polyA = par[0] + par[1] * energy;
- final double polyB = par[2] + par[3] * energy + par[4]
- * Math.pow(energy, 2);
+ public static final double correctTimeWalkPulseFitting(double time, double energy) {
+ final double polyA = par[0] + par[1]*energy;
+ final double polyB = par[2] + par[3] * energy + par[4] * Math.pow(energy, 2);
return time - (Math.exp(polyA) + polyB);
}
+
+
+
+
+
+
}
|