Print

Print


Commit in lcsim/src/org/lcsim/recon/postrecon/leptonID/electron on MAIN
IsolatedHighPElectronIdentifier.java+15-21.3 -> 1.4
Updated to be compatible with current SID production output.

lcsim/src/org/lcsim/recon/postrecon/leptonID/electron
IsolatedHighPElectronIdentifier.java 1.3 -> 1.4
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);
CVSspam 0.2.12


Use REPLY-ALL to reply to list

To unsubscribe from the LCD-CVS list, click the following link:
https://listserv.slac.stanford.edu/cgi-bin/wa?SUBED1=LCD-CVS&A=1