Author: [log in to unmask]
Date: Tue Sep 8 16:45:15 2015
New Revision: 3557
Log:
now everything is fine
Modified:
java/trunk/tracking/src/main/java/org/hps/recon/tracking/gbl/FittedGblTrajectory.java
java/trunk/tracking/src/main/java/org/hps/recon/tracking/gbl/GBLOutput.java
java/trunk/tracking/src/main/java/org/hps/recon/tracking/gbl/MakeGblTracks.java
Modified: java/trunk/tracking/src/main/java/org/hps/recon/tracking/gbl/FittedGblTrajectory.java
=============================================================================
--- java/trunk/tracking/src/main/java/org/hps/recon/tracking/gbl/FittedGblTrajectory.java (original)
+++ java/trunk/tracking/src/main/java/org/hps/recon/tracking/gbl/FittedGblTrajectory.java Tue Sep 8 16:45:15 2015
@@ -10,7 +10,7 @@
*/
public class FittedGblTrajectory {
public static enum GBLPARIDX {
- QOVERP(0),XTPRIME(1),YTPRIME(2),XT(3),YT(4);
+ QOVERP(0),YTPRIME(1),XTPRIME(2),XT(3),YT(4);
private int _value;
private GBLPARIDX(int value) {
_value = value;
Modified: java/trunk/tracking/src/main/java/org/hps/recon/tracking/gbl/GBLOutput.java
=============================================================================
--- java/trunk/tracking/src/main/java/org/hps/recon/tracking/gbl/GBLOutput.java (original)
+++ java/trunk/tracking/src/main/java/org/hps/recon/tracking/gbl/GBLOutput.java Tue Sep 8 16:45:15 2015
@@ -908,35 +908,35 @@
double qOverP = q / htf.p(Math.abs(B));
double phi = htf.phi0();
- _params.setElement(0, 0, qOverP);
- _params.setElement(0, 1, lambda);
- _params.setElement(0, 2, phi);
- _params.setElement(0, 3, xT);
- _params.setElement(0, 4, yT);
+ _params.setElement(0, FittedGblTrajectory.GBLPARIDX.QOVERP.getValue(), qOverP);
+ _params.setElement(0, FittedGblTrajectory.GBLPARIDX.YTPRIME.getValue(), lambda);
+ _params.setElement(0, FittedGblTrajectory.GBLPARIDX.XTPRIME.getValue(), phi);
+ _params.setElement(0, FittedGblTrajectory.GBLPARIDX.XT.getValue(), xT);
+ _params.setElement(0, FittedGblTrajectory.GBLPARIDX.YT.getValue(), yT);
}
public BasicMatrix getParams() {
return _params;
}
+ double getQoverP() {
+ return _params.e(0, FittedGblTrajectory.GBLPARIDX.QOVERP.getValue());
+ }
+
double getLambda() {
- return _params.e(0, 1);
- }
-
- double getQoverP() {
- return _params.e(0, 0);
+ return _params.e(0, FittedGblTrajectory.GBLPARIDX.YTPRIME.getValue());
}
double getPhi() {
- return _params.e(0, 2);
+ return _params.e(0, FittedGblTrajectory.GBLPARIDX.XTPRIME.getValue());
}
double getXt() {
- return _params.e(0, 3);
+ return _params.e(0, FittedGblTrajectory.GBLPARIDX.XT.getValue());
}
double getYt() {
- return _params.e(0, 4);
+ return _params.e(0, FittedGblTrajectory.GBLPARIDX.YT.getValue());
}
}
Modified: java/trunk/tracking/src/main/java/org/hps/recon/tracking/gbl/MakeGblTracks.java
=============================================================================
--- java/trunk/tracking/src/main/java/org/hps/recon/tracking/gbl/MakeGblTracks.java (original)
+++ java/trunk/tracking/src/main/java/org/hps/recon/tracking/gbl/MakeGblTracks.java Tue Sep 8 16:45:15 2015
@@ -5,12 +5,10 @@
import hep.physics.vec.Hep3Matrix;
import hep.physics.vec.Hep3Vector;
import hep.physics.vec.VecOp;
-
import java.util.ArrayList;
import java.util.List;
import java.util.logging.Level;
import java.util.logging.Logger;
-
import org.apache.commons.math3.util.Pair;
import org.hps.recon.tracking.StrategyType;
import org.hps.recon.tracking.TrackType;
@@ -22,6 +20,7 @@
import org.lcsim.event.EventHeader;
import org.lcsim.event.Track;
import org.lcsim.event.TrackerHit;
+import org.lcsim.event.base.BaseTrack;
import org.lcsim.fit.helicaltrack.HelicalTrackFit;
import org.lcsim.fit.helicaltrack.HelicalTrackHit;
import org.lcsim.fit.helicaltrack.HelixUtils;
@@ -88,8 +87,8 @@
SeedTrack seedTrack = (SeedTrack) fittedTraj.get_seed();
SeedCandidate trackseed = seedTrack.getSeedCandidate();
- // Create a new SeedTrack (SeedTrack extends BaseTrack)
- SeedTrack trk = new SeedTrack();
+ // Create a new BaeTrack
+ BaseTrack trk = new BaseTrack();
// Add the hits to the track
for (HelicalTrackHit hit : trackseed.getHits()) {
@@ -110,10 +109,10 @@
trk.setRefPointIsDCA(true);
// Set the strategy used to find this track
- trk.setStratetgy(trackseed.getSeedStrategy());
+// trk.setStratetgy(trackseed.getSeedStrategy());
// Set the SeedCandidate this track is based on
- trk.setSeedCandidate(trackseed);
+// trk.setSeedCandidate(trackseed);
// Add the GBL flag to the track type.
trk.setTrackType(TrackType.setGBL(seedTrack.getType(), true));
@@ -121,7 +120,7 @@
// Check the track - hook for plugging in external constraint
//if ((_trackCheck != null) && (! _trackCheck.checkTrack(trk))) continue;
// Add the track to the list of tracks
- tracks.add((Track) trk);
+ tracks.add(trk);
logger.info(String.format("helix chi2 %f ndf %d gbl chi2 %f ndf %d\n", helix.chisqtot(), helix.ndf()[0] + helix.ndf()[1], trk.getChi2(), trk.getNDF()));
if (logger.getLevel().intValue() <= Level.INFO.intValue()) {
for (int i = 0; i < 5; ++i) {
|