add a couple correlation plots

eic · Oct 24, 2022 · a1c5846 · a1c5846
1 parent c16a4f4
commit a1c5846
Show file tree

Hide file tree

Showing 3 changed files with 14 additions and 2 deletions.
diff --git a/src/Analysis.cxx b/src/Analysis.cxx
@@ -527,6 +527,15 @@ void Analysis::Prepare() {
       HS->DefineHist1D("DihPhiH",   "#phi_{h}", "",    NBINS,   -TMath::Pi(), TMath::Pi());
       HS->DefineHist1D("DihPhiR",   "#phi_{R}", "",    NBINS,   -TMath::Pi(), TMath::Pi());
       HS->DefineHist1D("DihPhiS",   "#phi_{S}", "",    NBINS,   -TMath::Pi(), TMath::Pi());
+      HS->DefineHist2D("DihPhiHvsPhiR", "#phi_{R}", "#phi_{h}", "", "",
+          NBINS, -TMath::Pi(), TMath::Pi(),
+          NBINS, -TMath::Pi(), TMath::Pi()
+          );
+      HS->DefineHist2D("DihThetaVsPh", "P_{h}", "#theta", "GeV", "",
+          NBINS, 1e-2, 40,
+          NBINS, 0, TMath::Pi(),
+          true
+          );
     }
     // -- jet kinematics
 #ifndef EXCLUDE_DELPHES
@@ -857,6 +866,8 @@ void Analysis::FillHistosDihadrons() {
     H->FillHist1D("DihPhiH",   dih->PhiH,   wTrack);
     H->FillHist1D("DihPhiR",   dih->PhiR,   wTrack);
     H->FillHist1D("DihPhiS",   dih->PhiS,   wTrack);
+    H->FillHist2D("DihPhiHvsPhiR", dih->PhiR, dih->PhiH,  wTrack);
+    H->FillHist2D("DihThetaVsPh",  dih->Ph,   dih->Theta, wTrack);
     });
   HD->ExecuteOps(true);
 }

diff --git a/src/Dihadrons.cxx b/src/Dihadrons.cxx
@@ -29,8 +29,9 @@ Dihadron::Dihadron(TLorentzVector vecHad1, TLorentzVector vecHad2, Kinematics *K
   Z  = vecP.Dot(vecPh) / vecP.Dot(vecQ);
   for(int h=0; h<2; h++)
     hadZ[h] = vecP.Dot(vecH[h]) / vecP.Dot(vecQ);
+  Ph     = vecPh.Vect().Mag();
   PhPerp = Kinematics::Reject(IvecPh.Vect(),IvecQ.Vect()).Mag();
-  Zeta = 2 * vecR.Dot(vecP) / vecPh.Dot(vecP);
+  Zeta   = 2 * vecR.Dot(vecP) / vecPh.Dot(vecP);
 
   // PhiH
   PhiH = Kinematics::PlaneAngle(

diff --git a/src/Dihadrons.h b/src/Dihadrons.h
@@ -11,7 +11,7 @@ class Dihadron {
     Dihadron() {};
     Dihadron(TLorentzVector vecHad1, TLorentzVector vecHad2, Kinematics *K);
     ~Dihadron() {};
-    Double_t Mh, MX, Z, PhPerp, Theta, PhiH, PhiR, PhiS, Zeta;
+    Double_t Mh, MX, Z, Ph, PhPerp, Theta, PhiH, PhiR, PhiS, Zeta;
     Double_t hadZ[2];
   private:
   ClassDef(Dihadron,1);