GSI Forum
GSI Helmholtzzentrum für Schwerionenforschung

Home » PANDA » PandaRoot » Tracking » RhoError class in pandaroot
RhoError class in pandaroot [message #16766] Tue, 03 June 2014 11:52 Go to previous message
Elisabetta Prencipe (2) is currently offline  Elisabetta Prencipe (2)
Messages: 214
Registered: February 2013
first-grade participant
From: *ikp.kfa-juelich.de
Dear Ralf, dear rho-developers,

I have a question how to evaluate the error of the momentum, position, theta, phi, distributions using rho candidates and rho lists, and I am wondering if below is the best way how to proceed. For example, if I want to evaluate the pull of quantities like momentum, or theta and phi, I need information about the reco variables (easy), the true value (easy) and the error distributions. I see that in RhoMath/ the information about the covariance error matrix is accessible. So I could write:


while (theAnalysis->GetEvent() && i++<nevts)
{

theAnalysis->FillList(muminus, "MuonAllMinus");

for (j=0;j<muminus.GetLength();++j)
{

// reco variables
hmomtrk->Fill(muminus[j]->P());
hthttrk->Fill(muminus[j]->P4().Theta());
hphitrk->Fill(muminus[j]->P4().Phi());

// error matrix : variance
RhoError tempvar = muminus[j]->Cov7(); //variance of 7 parameters: x,y,z,px,py,pz,E
double var_px = tempvar(3,3);
double var_py = tempvar(4,4);
double var_pz = tempvar(5,5);
double var_energy = tempvar(6,6);

}
}


This shoud give my access to the variance of px, py, pz, energy. Am I correct?
Sqrt() of what I get from here should deliver the error distribution of px, py, pz, E. Now my question is: is this the way to proceed to get the error distirbutions? and what about Theta() and Phi()? Should I combine, then, the information which I obtain, and get the error distribution for the angular variables? or is there a smarter way/funcion implemented to obtain the error distribution of kinematic variables, in pandaroot?

Thank you for your help,

Elisabetta
 
Read Message
Read Message
Read Message
Previous Topic: genfit in pandaroot - trunk rev 25180: test macros
Next Topic: [FIXED] PndTrkTracking2 statusflag array is not initialised
Goto Forum:
  


Current Time: Sat Dec 14 07:39:25 CET 2024

Total time taken to generate the page: 0.00764 seconds