From 7aacc2f29ff54229eb6809ca5d5d13bf9f739a8c Mon Sep 17 00:00:00 2001
From: Ziyue Zhang <Ziyue_Zhang@localhost.localdomain>
Date: Thu, 11 Mar 2021 14:25:56 -0600
Subject: [PATCH] WIP: Replace nu plot with y plot
---
benchmarks/dvmp/analysis/dvmp.h | 3 ++-
1 file changed, 2 insertions(+), 1 deletion(-)
diff --git a/benchmarks/dvmp/analysis/dvmp.h b/benchmarks/dvmp/analysis/dvmp.h
index 1784c289..5c743eae 100644
--- a/benchmarks/dvmp/analysis/dvmp.h
+++ b/benchmarks/dvmp/analysis/dvmp.h
@@ -209,6 +209,7 @@ namespace util {
{
//0:e0 1:p0 2:e1 3:p1 4:recoil system (without p1) 5:l1 from 4 6:l2 from 4
ROOT::Math::PxPyPzMVector q(parts[0] - parts[2]);
+ ROOT::Math::PxPyPzMVector k(parts[0]);
ROOT::Math::PxPyPzMVector P(parts[1]);
ROOT::Math::PxPyPzMVector Delta(parts[3] - parts[1]);//exact
ROOT::Math::PxPyPzMVector Delta_prime(parts[0] - parts[2] - parts[5] - parts[6]);//exclude gamma radiation in jpsi decay
@@ -221,7 +222,7 @@ namespace util {
//}else{
t = Delta.Dot(Delta);
//}
- double y = q.Dot(P) / parts[0].Dot(P);
+ double y = q.Dot(P) / k.Dot(P);
inv_quant quantities = {nu, Q2, Q2/2./P.mass()/nu, y, t};
return quantities;
}
--
GitLab