Commit in lcsim-contrib/src/main/java/org/lcsim/contrib/SteveMagill on MAIN | |||
CorrClusEsumDriver.java | +58 | added 1.1 |
diff -N CorrClusEsumDriver.java --- /dev/null 1 Jan 1970 00:00:00 -0000 +++ CorrClusEsumDriver.java 14 Feb 2012 19:40:47 -0000 1.1 @@ -0,0 +1,58 @@
+package org.lcsim.contrib.SteveMagill; + +import java.util.*; +import org.lcsim.util.aida.AIDA; +import org.lcsim.event.EventHeader; +import org.lcsim.util.Driver; +import org.lcsim.recon.cluster.util.BasicCluster; + +/** + * Diagnostic routine to evaluate performance of mip finder + */ + +public class CorrClusEsumDriver extends Driver +{ + public CorrClusEsumDriver() { + } + + private String _inclusname; + private AIDA aida = AIDA.defaultInstance(); + private boolean ccedb = true; + + public void process(EventHeader event) + { + int NClus = 0; + double ClEsum = 0; + int nhitsclus = 0; // number of hits in all mip clusters + + try + { + List<BasicCluster> inclusters = event.get(BasicCluster.class,_inclusname); + for (BasicCluster inclus : inclusters) + { + NClus++; + int numhits = inclus.getSize(); // number of hits in this cluster + if (ccedb) aida.cloud1D("Num Hits per Corrected Cluster").fill(numhits); + nhitsclus += inclus.getSize(); // sum hits from all clusters in event + ClEsum += inclus.getEnergy(); // energy of this cluster + if (ccedb) aida.cloud2D("E over Eclus vs Num hits").fill(numhits,10./inclus.getEnergy()); + if (ccedb && 10./inclus.getEnergy()<2.) aida.cloud2D("CF vs Num Hits").fill(numhits,10./inclus.getEnergy()); + } + + } + catch(java.lang.IllegalArgumentException ex) + { + System.out.println("requested object not found in event " + _inclusname); + } + + if (ccedb) aida.cloud1D("ESum of Corrected Clusters").fill(ClEsum); + if (ccedb) aida.cloud1D("Number of Corrected Clusters").fill(NClus); + + } + + public void setInputClusName(String inclusname) + { + _inclusname = inclusname; + } + +}
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