1 package org.lcsim.mc.fast.cluster.ronan;
2
3 import org.lcsim.util.aida.AIDA;
4 import hep.physics.particle.Particle;
5 import java.util.Random;
6
7 public class ReconHADCluster extends ReconCluster {
8 ReconHADCluster(ClusterResolutionTables parm, Random rand, Particle mcp, boolean hist) {
9 super(parm, rand, mcp, hist);
10
11 a = parm.getHADResolution();
12 b = parm.getHADConstantTerm();
13 c = parm.getHADPositionError();
14 d = parm.getHADAlignmentError();
15
16 smear(rand, hist);
17 }
18
19 protected void smearPosition(Random rand, double E, boolean hist) {
20 double transSigma = c / Math.sqrt(E) + d;
21 transDist = transSigma * rand.nextGaussian();
22 if (hist)
23 AIDA.defaultInstance().cloud1D("HAD: transDist").fill(transDist);
24 smearPosition(rand);
25 }
26 }