sync with Zhichen, fix W_lep order

Author: Moelf

src/ana.jl

@@ -15,5 +15,12 @@ function Find_Z_Pairs(v_l_pids, v_l_tlv)
             end
         end
     end
-    return Z_pair, setdiff(1:4, Z_pair), M
+    W_pair = setdiff(1:4, Z_pair)
+
+    # order W lep by pT
+    if pt(v_l_tlv[W_pair[1]]) < pt(v_l_tlv[W_pair[2]])
+        W_pair[1], W_pair[2] = W_pair[2], W_pair[1]
+    end
+
+    return Z_pair, W_pair, M
 end

src/mainlooper.jl

@@ -202,7 +202,6 @@ function main_looper(mytree, sumWeight, dict, pusher!, models,
         Wlep1_dphi = LorentzVectorHEP.phi_mpi_pi(Z_phi - Wlep1_phi)
         Wlep2_dphi = LorentzVectorHEP.phi_mpi_pi(Z_phi - Wlep2_phi)
         Wleps_deta = abs(Wlep1_eta - Wlep2_eta)
-        Wleps_dphi = LorentzVectorHEP.phi_mpi_pi(Wlep1_phi - Wlep2_phi)
         pt_4l = pt(sum(v_l_tlv))
 
         MET_dPhi = LorentzVectorHEP.phi_mpi_pi(Z_phi - METPhi)