lcsim/src/org/lcsim/contrib/uiowa/MuonFinder
diff -u -r1.2 -r1.3
--- MipTrackMap.java 30 Sep 2008 01:11:18 -0000 1.2
+++ MipTrackMap.java 30 Sep 2008 01:25:23 -0000 1.3
@@ -384,4 +384,23 @@
return layer;
}
+ protected Hep3Vector getClusterDirection(Cluster clus){
+
+ Hep3Vector trackMuonUnit = new BasicHep3Vector();
+ if( clus.getCalorimeterHits().size() > 1){
+ CalorimeterHit muonhit1 = clus.getCalorimeterHits().get(clus.getCalorimeterHits().size()-1);
+ CalorimeterHit muonhit0 = clus.getCalorimeterHits().get(0);
+ Hep3Vector muonpos1 = new BasicHep3Vector(muonhit1.getPosition());
+ Hep3Vector muonpos0 = new BasicHep3Vector(muonhit0.getPosition());
+ Hep3Vector trackMuon = new BasicHep3Vector();
+ if(muonpos0.magnitude() > muonpos1.magnitude()){
+ trackMuon = VecOp.sub(muonpos0, muonpos1);
+ }else{
+ trackMuon = VecOp.sub(muonpos1, muonpos0);
+ }
+ trackMuonUnit = VecOp.unit(trackMuon);
+ } else { trackMuonUnit = null;}
+
+ return trackMuonUnit;
+ }
}