lcsim/src/org/lcsim/contrib/uiowa/template
diff -u -r1.4 -r1.5
--- NonTrivialPFA.java 4 Feb 2006 01:25:31 -0000 1.4
+++ NonTrivialPFA.java 8 Feb 2006 19:06:25 -0000 1.5
@@ -4,6 +4,7 @@
import org.lcsim.util.*;
import org.lcsim.event.*;
+import org.lcsim.event.util.*;
import org.lcsim.recon.cluster.mst.*;
import org.lcsim.recon.cluster.mipfinder.*;
import org.lcsim.recon.cluster.clumpfinder.*;
@@ -260,6 +261,26 @@
*/
add(new DebugInfoParticleList("hadron particles"));
add(new EnergySumPlotter("hadron particles", "test.aida"));
+
+ // Go to work finding the correction for missing energy:
+ hitmap.MapToHitMapDriver convertECAL = new hitmap.MapToHitMapDriver();
+ hitmap.MapToHitMapDriver convertHCAL = new hitmap.MapToHitMapDriver();
+ convertECAL.setInputMap("input hit map ecal");
+ convertHCAL.setInputMap("input hit map hcal");
+ convertECAL.setOutputHitMap("converted ecal");
+ convertHCAL.setOutputHitMap("converted hcal");
+ add(convertECAL);
+ add(convertHCAL);
+ hitmap.HitMapAddDriver adder = new hitmap.HitMapAddDriver();
+ adder.addInputHitMap("converted ecal");
+ adder.addInputHitMap("converted hcal");
+ adder.setOutputHitMap("converted all");
+ add(adder);
+ // Set up the MC list
+ CreateFinalStateMCParticleList mcListMaker = new CreateFinalStateMCParticleList("Gen");
+ add(mcListMaker);
+ // Make plots:
+ add(new CorrectedEnergySumPlotter("converted all", "hadron particles", "GenFinalStateParticles", "reco-corrected.aida"));
}
public void process(EventHeader event) {