The Clusterer currently just adds the energy in the included cells,
which is in raw energy. You need to correct for the sampling fraction,
which in turn depends on which detector you are using. I think there are
some appropriate numbers in one of Richard Dubois' old talks which are
linked to the LCD web site.
Tony
> -----Original Message-----
> From: Jeff B. Gronberg [mailto:[log in to unmask]]
> Sent: Wednesday, April 12, 2000 4:05 PM
> To: [log in to unmask]
> Subject: difficulty with SimpleClusterBuilder
>
>
> Hello all,
>
> I'm having difficulty interpreting the
> output of SimpleClusterBuilder. When I add
> the energy in all EM and HAD clusters in an event
> I get a number that varies between 0. and 13.
> Since I am running on pybms-ZZ-500-990115-Small-sim
> I would naively expect something more in the region
> of 500.
>
> Thanks, Jeff
>
> My code is:
>
> final public class LCDTrackDriver extends Driver
> {
> public LCDTrackDriver()
> {
> setHistogramLevel(ProcessorContext.LEVEL_HIGH);
> add(new MCFast()); add(new TrackAnalyzer());
> add(new SimpleClusterBuilder()); add(new ClusterAnalyzer());
> add(new JBGAnalyzer());
> }
> }
> final class JBGAnalyzer extends AbstractProcessor
> {
> public final void process(LCDEvent event)
> {
> ClusterList EMcl = event.getEMClusters();
> ClusterList HADcl = event.getHADClusters();
> double totE = 0.0;
> {
> Enumeration e = EMcl.getClusters();
> while (e.hasMoreElements())
> {
> Cluster c = (Cluster) e.nextElement();
> totE += c.getEnergy();
> }
> }
> {
> Enumeration e = HADcl.getClusters();
> while (e.hasMoreElements())
> {
> Cluster c = (Cluster) e.nextElement();
> totE += c.getEnergy();
> }
> }
> histogram("Tot shw E").fill(totE);
> }
> }
>
