Author: [log in to unmask]
Date: Mon May 16 12:44:19 2016
New Revision: 4352
Log: (empty)
Modified:
java/trunk/users/src/main/java/org/hps/users/spaul/moller/MollerBeamtiltAnalysis.java
Modified: java/trunk/users/src/main/java/org/hps/users/spaul/moller/MollerBeamtiltAnalysis.java
=============================================================================
--- java/trunk/users/src/main/java/org/hps/users/spaul/moller/MollerBeamtiltAnalysis.java (original)
+++ java/trunk/users/src/main/java/org/hps/users/spaul/moller/MollerBeamtiltAnalysis.java Mon May 16 12:44:19 2016
@@ -86,6 +86,23 @@
pypz_vs_mass[i].fill(mass, pypz );
timediff[i].fill(top.getClusters().get(0).getCalorimeterHits().get(0).getTime()
-bottom.getClusters().get(0).getCalorimeterHits().get(0).getTime());
+
+
+ double theta, phi;
+ double pxpz_top = top.getMomentum().x()/top.getMomentum().z();
+ double pypz_top = top.getMomentum().y()/top.getMomentum().z();
+
+ theta = Math.hypot(pypz_top-(-.0008), pxpz_top-.0295);
+ phi = Math.atan2(pypz_top-(-.0008), pxpz_top-.0295);
+ thetaPhi[i].fill(theta, phi);
+
+ double pxpz_bot = bottom.getMomentum().x()/bottom.getMomentum().z();
+ double pypz_bot = bottom.getMomentum().y()/bottom.getMomentum().z();
+
+ theta = Math.hypot(pypz_bot-(-.0008), pxpz_bot-.0295);
+ phi = Math.atan2(pypz_bot-(-.0008), pxpz_bot-.0295);
+ thetaPhi[i].fill(theta, phi);
+
/*if(moreEnergetic.getMomentum().y() > 0)
{
pypz_tophighE.fill(pypz);
@@ -140,6 +157,8 @@
hpypz_botHighE[], hpxpz_botHighE[],
hpypz_mid[], hpxpz_mid[];
+ private IHistogram2D[] thetaPhi;
+
public double getMaxVtxChi2() {
return _maxVtxChi2;
}
@@ -211,6 +230,11 @@
@Override
public void startOfData(){
AIDA aida = AIDA.defaultInstance();
+
+
+
+ thetaPhi = new IHistogram2D[3];
+
hpypz = new IHistogram1D[3];
hpxpz = new IHistogram1D[3];
hpypz_mid = new IHistogram1D[3];
@@ -240,6 +264,8 @@
timediff= new IHistogram1D[3];
for(int i = 0; i< 3; i++){
+
+ thetaPhi[i] = aida.histogram2D(mollerCollections[i]+"/"+"theta vs phi", 100, 0, .2, 628, -3.14, 3.14);
hpypz[i] = aida.histogram1D(mollerCollections[i]+"/"+"pypz", 60, -.005,.005);
hpxpz[i] = aida.histogram1D(mollerCollections[i]+"/"+"pxpz", 60, .025,.035);
|