fix kappa extraction for dyGroomer

mverwe · Mar 6, 2021 · 6ac884d · 6ac884d
1 parent eb4d7ab
commit 6ac884d
Show file tree

Hide file tree

Showing 2 changed files with 6 additions and 4 deletions.
diff --git a/include/dyGroomer.hh b/include/dyGroomer.hh
@@ -215,7 +215,8 @@ std::vector<fastjet::PseudoJet> dyGroomer::doGrooming()
      drBranches_.push_back(ndrop);
      zg_.push_back(zg);
      dr12_.push_back(deltaR);
-     kappa_.push_back(kappa);
+     //kappa_.push_back(kappa);
+     kappa_.push_back(1./getKappa(pt,deltaR,zg));
    }
 
     // Compute the n-subjettiness ratio

diff --git a/runMiguelsVariablesBkg.cc b/runMiguelsVariablesBkg.cc
@@ -31,6 +31,7 @@ using namespace fastjet;
 // This class runs time reclustering with background
 
 // ./runMiguelsVariablesBkg -hard samples/PythiaEventsTune14PtHat120_10k.pu14 -pileup samples/ThermalEventsMult7000PtAv1.20_0.pu14 -nev 10
+// ./runMiguelsVariablesBkg -hard samples/nikhef/PythiaEventsTune14PtHat120_0.pu14 -pileup samples/nikhef/ThermalEventsMult7000PtAv1.20_0.pu14 -nev 10
 //https://indico.cern.ch/event/974749/contributions/4104667/attachments/2147059/3619205/Strong2020-JetObservables.pdf
 //(SD_tau2, kappaTD), (tau2,SD_pt), (SD_rz, kappaTD), (SD_ptD, kappaktD)
 
@@ -180,7 +181,7 @@ int main (int argc, char ** argv) {
     dyGroomer dygTDSig(2);
     jetCollection jetCollectionSigDYTD(dygTDSig.doGrooming(jetCollectionSig));
     jetCollectionSigDYTD.addVector("kappaSigDYTD",    dygTDSig.getKappas());
-
+    
     dyGroomer dygKTDSig(1);
     jetCollection jetCollectionSigDYKTD(dygKTDSig.doGrooming(jetCollectionSig));
     jetCollectionSigDYKTD.addVector("kappaSigDYKTD",    dygKTDSig.getKappas());
@@ -268,8 +269,8 @@ int main (int argc, char ** argv) {
     //   }
     //   ++ijc;
     // }
-    // jetCollectionCSSDBeta00Z01.addVector("widthCSSD", widthSig);
-    // jetCollectionCSSDBeta00Z01.addVector("pTDCSSD", pTDSig);
+    // jetCollectionCSSDBeta00Z01.addVector("widthCSSD", widthSigCSSD);
+    // jetCollectionCSSDBeta00Z01.addVector("pTDCSSD", pTDSigCSSD);
     // jetCollectionCSSDBeta00Z01.addVector("tau2CSSD", tau2CSSD);
 
     // //---------------------------------------------------------------------------