lcsim/src/org/lcsim/recon/muon

diff -u -r1.5 -r1.6
--- MuonReco.java	6 Dec 2005 16:24:08 -0000	1.5
+++ MuonReco.java	12 Dec 2005 22:36:46 -0000	1.6
@@ -89,7 +89,7 @@

             
             if(!stpr.getStopTkELow())
                 muSegmentFinder.findSegment(event,stpr);

-            

+            rpVMR = stpr.getNewRp();

             // If segments pass cuts, then create a MuonCandidate object
             // and add it to the list
             System.out.println("MR:MuonCandidate");

@@ -105,10 +105,18 @@

                 muSegmentFinder.getNHits(), muSegmentFinder.getNLayersHit(),
                 muSegmentFinder.getHits());
                 muonList.add(muon);

+		double x = rpVMR[0];
+		double y = rpVMR[1];
+		double rho = Math.sqrt(x*x+y*y);
+		double theta = Math.atan2(rho,rpVMR[2])*180/Math.PI;

 		System.out.println("MuonReco: #muons="+muonList.getNMuons()

-		   +", #hits: Ecal="+ecalSegmentFinder.getNHits()
-		   +", Hcal="+hcalSegmentFinder.getNHits()
-		   +", muon="+muSegmentFinder.getNHits()  );

+ 	          +", #hits: Ecal="+ecalSegmentFinder.getNHits()
+		  +", Hcal="+hcalSegmentFinder.getNHits()
+		   +", muon="+muSegmentFinder.getNHits()
+				   +", E="+stpr.partPabs()
+				   +", z="+rpVMR[2]
+				   +", rho="+rho
+				   +", theta="+theta  );

             }
             // Reset for the next track
             ecalSegmentFinder.reset();		// C.M.- Include the Ecal inthe code