def v15_corrections(self, vertex_z):
if abs(self.etas1) < 1.5:
R = self.v15_RZ_1stSampling_cscopt2
Z = R * sinh(self.etas1)
else:
Z = self.v15_RZ_1stSampling_cscopt2
R = Z / sinh(self.etas1)
eta_corrected = asinh((Z - vertex_z) / R)
E_corrected = self.v15_E_corrected
pt_corrected = E_corrected / cosh(eta_corrected)
v = Fourvec_PtEtaPhiE(pt_corrected, eta_corrected, self.phi, E_corrected)
v.isConv = self.isConv
return v
def v16_corrections(self):
E_corrected = self.v16_E_corrected()
v = Fourvec_PtEtaPhiE(self.cl.pt, self.etas1, self.phi, E_corrected)
v.isConv = self.isConv
return v