2 #ifndef CBM_ANAJPSI_KINEMATIC_PARAMS_H
3 #define CBM_ANAJPSI_KINEMATIC_PARAMS_H
7 #include "TLorentzVector.h"
10 #define M2E 2.6112004954086e-7
30 Double_t energyP = TMath::Sqrt(momP.Mag2() +
M2E);
31 TLorentzVector lorVecP(momP, energyP);
35 Double_t energyM = TMath::Sqrt(momM.Mag2() +
M2E);
36 TLorentzVector lorVecM(momM, energyM);
38 TVector3 momPair = momP + momM;
39 Double_t energyPair = energyP + energyM;
40 Double_t ptPair = momPair.Perp();
41 Double_t pzPair = momPair.Pz();
43 0.5 * TMath::Log((energyPair + pzPair) / (energyPair - pzPair));
44 Double_t anglePair = lorVecM.Angle(lorVecP.Vect());
45 Double_t theta = 180. * anglePair /
TMath::Pi();
47 2. * TMath::Sin(anglePair / 2.) * TMath::Sqrt(momM.Mag() * momP.Mag());
70 Double_t ptPair = momPair.Perp();
71 Double_t pzPair = momPair.Pz();
73 0.5 * TMath::Log((energyPair + pzPair) / (energyPair - pzPair));
74 Double_t anglePair = lorVecM.Angle(lorVecP.Vect());
75 Double_t theta = 180. * anglePair /
TMath::Pi();
77 2. * TMath::Sin(anglePair / 2.)