diff --git a/src/MVAvarsJABDTthw.cpp b/src/MVAvarsJABDTthw.cpp
index f9867abc1b2f66e41553782bc6ebcff0baee6590..58f4e92c1f19369f338284d129a623b7cccbe2e8 100644
--- a/src/MVAvarsJABDTthw.cpp
+++ b/src/MVAvarsJABDTthw.cpp
@@ -85,6 +85,8 @@ MVAvarsJABDTthw::MVAvarsJABDTthw()
                     ,"Reco_JABDT_tHW_log_h_pt"
                     ,"Reco_JABDT_tHW_log_wb_m"
                     ,"Reco_JABDT_tHW_log_wb_pt"
+                    ,"Reco_JABDT_tHW_log_whad_m"
+                    ,"Reco_JABDT_tHW_log_whad_pt"
 
                     ,"Reco_JABDT_tHW_abs_top_eta__M__wb_eta"
                     ,"Reco_JABDT_tHW_abs_btop_eta"
@@ -200,6 +202,8 @@ void MVAvarsJABDTthw::FillMVAvarMap(const std::vector<TLorentzVector> &selectedL
     variableMap["Reco_JABDT_tHW_log_h_pt"]        = log(vectors["higg"].Pt());
     variableMap["Reco_JABDT_tHW_log_wb_m"]        = log(vectors["wb"].M());
     variableMap["Reco_JABDT_tHW_log_wb_pt"]       = log(vectors["wb"].Pt());
+    variableMap["Reco_JABDT_tHW_log_whad_m"]      = log(vectors["whad"].M());
+    variableMap["Reco_JABDT_tHW_log_whad_pt"]     = log(vectors["whad"].Pt());
 
     variableMap["Reco_JABDT_tHW_abs_top_eta__M__wb_eta"]      = fabs(vectors["top"].Eta() - vectors["wb"].Eta());
     variableMap["Reco_JABDT_tHW_abs_btop_eta"]                = fabs(vectors["btop"].Eta());