Skip to content
Snippets Groups Projects
Commit 7aacc2f2 authored by Ziyue Zhang's avatar Ziyue Zhang
Browse files

WIP: Replace nu plot with y plot

parent cacdcff8
No related branches found
No related tags found
1 merge request!39Replacing nu with y
This commit is part of merge request !39. Comments created here will be created in the context of that merge request.
...@@ -209,6 +209,7 @@ namespace util { ...@@ -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 //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 q(parts[0] - parts[2]);
ROOT::Math::PxPyPzMVector k(parts[0]);
ROOT::Math::PxPyPzMVector P(parts[1]); ROOT::Math::PxPyPzMVector P(parts[1]);
ROOT::Math::PxPyPzMVector Delta(parts[3] - parts[1]);//exact 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 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 { ...@@ -221,7 +222,7 @@ namespace util {
//}else{ //}else{
t = Delta.Dot(Delta); 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}; inv_quant quantities = {nu, Q2, Q2/2./P.mass()/nu, y, t};
return quantities; return quantities;
} }
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment