hps-java/src/main/java/org/lcsim/hps/users/phansson
diff -u -r1.2 -r1.3
--- ParticleHelixProducer.java 10 Oct 2012 18:53:10 -0000 1.2
+++ ParticleHelixProducer.java 10 Oct 2012 19:57:34 -0000 1.3
@@ -58,6 +58,14 @@
// Name of StripHit1D output collection.
private String trackOutputCollectionName = "MCParticle_HelicalTrackFit";
+ private AIDA aida = AIDA.defaultInstance();
+ IPlotter plotter_trkparams;
+ ICloud1D h_pt;
+ ICloud1D h_phi0;
+ ICloud1D h_d0;
+ ICloud1D h_R;
+ ICloud1D h_slope;
+ ICloud1D h_z0;
public void setDebug(boolean debug) {
this.debug = debug;
@@ -88,6 +96,23 @@
detToTrk = trackerhitutils.detToTrackRotationMatrix();
+ plotter_trkparams = aida.analysisFactory().createPlotterFactory().create();
+ plotter_trkparams.createRegions(3,2);
+ h_pt = aida.cloud1D("Track pT");
+ h_R = aida.cloud1D("Track R");
+ h_phi0 = aida.cloud1D("Track phi0");
+ h_d0 = aida.cloud1D("Track d0");
+ h_slope = aida.cloud1D("Track slope");
+ h_z0 = aida.cloud1D("Track z0");
+ plotter_trkparams.region(0).plot(h_pt);
+ plotter_trkparams.region(1).plot(h_R);
+ plotter_trkparams.region(2).plot(h_phi0);
+ plotter_trkparams.region(3).plot(h_d0);
+ plotter_trkparams.region(4).plot(h_slope);
+ plotter_trkparams.region(5).plot(h_z0);
+ plotter_trkparams.show();
+
+
}
/**
@@ -163,6 +188,12 @@
}
+ h_pt.fill(htf.pT(bfield));
+ h_R.fill(htf.R());
+ h_d0.fill(htf.dca());
+ h_phi0.fill(htf.phi0());
+ h_slope.fill(htf.slope());
+ h_z0.fill(htf.z0());
}
}