mcd-analysis/src/main/java/org/lcsim/mcd/analysis
diff -u -r1.1 -r1.2
--- HiggsAnalysisDriver.java 19 Mar 2013 19:58:49 -0000 1.1
+++ HiggsAnalysisDriver.java 19 Mar 2013 20:27:06 -0000 1.2
@@ -189,14 +189,18 @@
private void makeThrustPlots(List<MCParticle> mcpl, String prefix) {
EventShape es = getEventShape(mcpl);
Hep3Vector thrust = es.thrust();
- Hep3Vector majAxis = es.majorAxis();
- Hep3Vector minAxis = es.minorAxis();
+ Hep3Vector ThrustAxis = es.thrustAxis();
double ob = es.oblateness();
- if(thrust.x()>0) aida.cloud1D(prefix+"/Thrust").fill(thrust.x());
+ if(thrust.x()>0) aida.cloud1D(prefix+"/Thrst").fill(thrust.x());
aida.cloud1D(prefix+"/Major Axis").fill(thrust.y());
aida.cloud1D(prefix+"/Minor Axis").fill(thrust.z());
aida.cloud1D(prefix+"/Oblateness").fill(ob);
+ aida.cloud2D(prefix+"/Correl Thrust Major").fill(thrust.x(), thrust.y());
+ aida.cloud2D(prefix+"/Correl Thrust Oblateness").fill(thrust.x(), ob);
+ aida.cloud1D(prefix+"/Thrust Axis CosTheta").fill(VecOp.cosTheta(ThrustAxis));
+ aida.cloud2D(prefix+"/Correl Thrust vs Thrust Axis CosT"
+ ).fill(thrust.x(), VecOp.cosTheta(ThrustAxis));
}
private void makeJetPlots(EventHeader event, String prefix, String decay,String evt) {