GSI Forum
GSI Helmholtzzentrum für Schwerionenforschung

Home » PANDA » PandaRoot » Analysis » vertex fitting
vertex fitting [message #17255] Mon, 08 September 2014 13:26 Go to previous message
Alexandros is currently offline  Alexandros
Messages: 51
Registered: January 2014
Location: KVI
continuous participant
From: *kvi-cart.rug.nl
Hi all,

Is there something wrong with the code below??Because I get a segmentation break with these lines and without everything works perfectly..
The problem seems to be in the red line...
Any ideas??

// store the 4-vector of the truth matched candidate (or a dummy, if not matched to keep ntuple consistent)
RhoCandidate *truth = d0[j]->GetMcTruth();
TLorentzVector lv;
if (truth) lv = truth->P4();
qa.qaP4("trd0", lv, nd0);

PndKinVtxFitter *vtxfitter=new PndKinVtxFitter(truth); // instantiate a vertex fitter
vtxfitter->Fit();


RhoCandidate *d0fitvtx = truth->GetFit(); // access the fitted cand
Float_t chi2_vtx = vtxfitter->GetChi2(); // access chi2 of fit
Float_t prob_vtx = vtxfitter->GetProb(); // access probability of fit

TVector3 d0Vtx = d0fitvtx->Pos(); // and the decay vertex position
nd0->Column("vtxprob", (Float_t) prob_vtx);
nd0->Column("vtxchi2", (Float_t) chi2_vtx);
nd0->Column("vtxposx", (Float_t) d0Vtx.Dead));
nd0->Column("vtxposy", (Float_t) d0Vtx.Y());
nd0->Column("vtxposz", (Float_t) d0Vtx.Z());
 
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Previous Topic: Is there a fast way of histogramming RhoCandLists?
Next Topic: Ds reconstruction - test with trunk rev-25813
Goto Forum:
  


Current Time: Thu Dec 13 16:58:32 CET 2018

Total time taken to generate the page: 0.02528 seconds