lcsim/src/org/lcsim/recon/postrecon/leptonID/electron
diff -u -r1.3 -r1.4
--- IsolatedHighPElectronIdentifier.java 17 Oct 2012 06:36:23 -0000 1.3
+++ IsolatedHighPElectronIdentifier.java 23 Apr 2013 22:30:00 -0000 1.4
@@ -53,6 +53,9 @@
//
// Match Tracks with Recon Particles
//
+ System.out.println("This is the one!");
+
+ // if (1==1) return;
List<Track> allTracks = event.get(Track.class,"Tracks");
List<Track> unusedTracks = new ArrayList<Track>(allTracks);
@@ -101,8 +104,13 @@
{
Hep3Vector cp = new BasicHep3Vector(c.getPosition());
double ct = (tp.x()*cp.x()+tp.y()*cp.y()+tp.z()*cp.z())/tp.magnitude()/cp.magnitude();
+ // double hec = 0.0;
+ // double hec = c.getSubdetectorEnergies()[3] + c.getSubdetectorEnergies()[7];
+ double hec = c.getSubdetectorEnergies()[1];
+ // System.out.println(" ....... hec from c.getSubdetectorEnergies()[1] = " + hec);
- double hec = 0.0;
+
+ double hithec = 0.0;
List<CalorimeterHit> calhitsList = c.getCalorimeterHits();
for (CalorimeterHit pfochit : calhitsList) {
@@ -110,10 +118,14 @@
String dname = _decoder.getSubdetector().getName();
// System.out.println(" name of subdetector for this hit " + dname);
if (dname.equals("HcalBarrel") || dname.equals("HcalEndcap")) {
- hec += pfochit.getCorrectedEnergy();
+ hithec += pfochit.getCorrectedEnergy();
// System.out.println(" HCAL energy contribution = "+pfochit.getCorrectedEnergy());
}
}
+ // System.out.println(" hithec = "+hithec);
+
+ // System.out.println(" . hec from c.getSubdetectorEnergies()[1] = " + hec + " ... hithec from CalorimeterHit getCorrectedEnergy() for HcalBarrel+HcalEndcap = " + hithec );
+
if(ct > conect)
{
@@ -124,6 +136,7 @@
double f = HcalE/tp.magnitude();
if(!eid)
{
+ System.out.println(" f = "+f);
if(f < 0.04)xtrae.add(trck); // .04
}
ttocll.put(trck,tcl);