diff --git a/interface/JMEJetAnalyzer.h b/interface/JMEJetAnalyzer.h
index a3aaddc..4c40654 100644
--- a/interface/JMEJetAnalyzer.h
+++ b/interface/JMEJetAnalyzer.h
@@ -51,6 +51,7 @@ class JMEJetAnalyzer : public JME::PhysicsObjectAnalyzer
         std::vector<float>& betaClassic = tree["betaClassic"].write<std::vector<float>>();
         std::vector<float>& betaStarClassic = tree["betaStarClassic"].write<std::vector<float>>();
         std::vector<float>& dZ = tree["dZ"].write<std::vector<float>>();
+        std::vector<float>& d0 = tree["d0"].write<std::vector<float>>();
         std::vector<float>& DRweighted = tree["DRweighted"].write<std::vector<float>>();
         std::vector<float>& fRing0 = tree["fRing0"].write<std::vector<float>>();
         std::vector<float>& fRing1 = tree["fRing1"].write<std::vector<float>>();
diff --git a/interface/PhysicsObjectAnalyzer.h b/interface/PhysicsObjectAnalyzer.h
index 5d64d2b..ab21517 100644
--- a/interface/PhysicsObjectAnalyzer.h
+++ b/interface/PhysicsObjectAnalyzer.h
@@ -58,11 +58,11 @@ namespace JME {
         protected:
             std::vector<PtEtaPhiEVector>& p4 = tree["p4"].write<std::vector<PtEtaPhiEVector>>();
             std::vector<float>& y = tree["y"].write<std::vector<float>>();
-            std::vector<int8_t>& charge = tree["charge"].write<std::vector<int8_t>>();
+            std::vector<float>& charge = tree["charge"].write<std::vector<float>>();
 
             std::vector<bool>& is_matched = tree["has_gen_particle"].write<std::vector<bool>>();
             std::vector<PtEtaPhiEVector>& gen_p4 = tree["gen_p4"].write<std::vector<PtEtaPhiEVector>>();
             std::vector<float>& gen_y = tree["gen_y"].write<std::vector<float>>();
-            std::vector<int8_t>& gen_charge = tree["gen_charge"].write<std::vector<int8_t>>();
+            std::vector<float>& gen_charge = tree["gen_charge"].write<std::vector<float>>();
     };
 }
diff --git a/src/JMEJetAnalyzer.cc b/src/JMEJetAnalyzer.cc
index 6fee80d..cf64ef9 100644
--- a/src/JMEJetAnalyzer.cc
+++ b/src/JMEJetAnalyzer.cc
@@ -294,7 +294,7 @@ void JMEJetAnalyzer::computeBetaStar(const pat::Jet& jet, const std::vector<reco
     int nCh_tmp(0), nNeutrals_tmp(0);
     float sumTkPt(0.0);
     float beta_tmp(0.0), betaStar_tmp(0.0), betaStarClassic_tmp(0.0), betaClassic_tmp(0.0);
-    float pTMax(0.0), pTMaxChg(0.0),pTMaxNeutral(0.0), dZ2(-999);
+    float pTMax(0.0), pTMaxChg(0.0),pTMaxNeutral(0.0), dZ2(-999),d02(-999);
     float sumW(0.0), sumW2(0.0), sumWdR2(0.0);
     float sum_deta(0.0),sum_dphi(0.0),sum_deta2(0.0),sum_dphi2(0.0),sum_detadphi(0.0),Teta(0.0),Tphi(0.0);
     float ave_deta(0.0), ave_dphi(0.0);
@@ -351,6 +351,7 @@ void JMEJetAnalyzer::computeBetaStar(const pat::Jet& jet, const std::vector<reco
                 if (lPack->pt() > pTMaxChg) {
                     pTMaxChg = lPack->pt();
                     dZ2 = lPack->dz();
+                    d02 = lPack->dxy();
                 }
 
                 float tkpt = lPack->pt();
@@ -454,6 +455,7 @@ void JMEJetAnalyzer::computeBetaStar(const pat::Jet& jet, const std::vector<reco
     }
 
     dZ.push_back(dZ2);
+    d0.push_back(d02);
     nCh.push_back(nCh_tmp);
     nNeutrals.push_back(nNeutrals_tmp);
     axisMajor.push_back(axis1);