i have created 1000 events for each event there are produced particles and i have set Px , Py, Pz for all those with the help of poisson distribution. and mass with then help of grandom now i want to reconstruct mass of dimuon using lorentz vector, how i can do?
Hi, welcome to the ROOT forum!
You can take a look at this RDataFrame tutorial for an example analysis with invariant masses:
https://root.cern.ch/doc/master/df106__HiggsToFourLeptons_8py.html
The key part is this one:
ROOT::Math::PtEtaPhiEVector p1(pt[0], eta[0], phi[0], e[0]);
ROOT::Math::PtEtaPhiEVector p2(pt[1], eta[1], phi[1], e[1]);
ROOT::Math::PtEtaPhiEVector p3(pt[2], eta[2], phi[2], e[2]);
ROOT::Math::PtEtaPhiEVector p4(pt[3], eta[3], phi[3], e[3]);
return (p1 + p2 + p3 + p4).M() / 1000;
So you use your kinematic variables to build Lorentz vectors, which you sum and then you get the invariant mass with .M().
Note that since, you have Px, Py, Pz, you need to use the ROOT::Math::PxPyPzMVector instead (I presume that you either have the mass available too, or you’re going to assume that the particles are massless).
How i can choose these p1, p2 vector among all of available vector that i have constructed?
That depends, how did you construct the vectors, and how do you want to select them? Are you using RDataFrame?
Your question is a bit too vague to be answered, maybe it helps if you explain again more concisely and attach your code so we understand better what you mean?