GSI Forum - RDF feed
https://forum.gsi.de/index.php
negative chi2 in PndKinVtxFitter !
https://forum.gsi.de/index.phpindex.php?t=rview&goto=12751&th=3385#msg_12751
In the simulation (rev 13401) of pbar p --> Ds+ Ds-; Ds+ --> phi pi+, Ds- --> phi e- nuebar I found in some cases negative(!) chi2 values in the phi-->K+K- vertex fit, using PndKinVtxFitter. The occurrence is on the level of ~1% with respect to the number of all K+K- combinations, or even more than 5% with respect to the number of analyzed events. Sometimes extremely small positive chi2 values << 1 (e.g. 10**-5 or 10**-5) are seen which is also strange.
It also occurs if I fit the Ds+ and Ds- vertices with three charged tracks.
I then checked the macro for the CT TDR benchmark channel Ds+ Ds- in macro/run/tdrct/psi3770: it also occurs there. (I am now running the eta_c macro in macro/run/tdrct/eta_c, but I have no doubt it will also occur there.)
Usually people plot chi2 distributions only for positive values, so that this problem may have escaped our attention.
My question(s): has anybody else taken notice of this problem so far? Is someone willing to look into the PndKinVtxFitter code to find out where things go wrong? How much confidence can we have in a vertex fitter doing such things?
Albrecht
]]>Albrecht Gillitzer2011-10-26T08:21:16-00:00Re: negative chi2 in PndKinVtxFitter !
https://forum.gsi.de/index.phpindex.php?t=rview&goto=12752&th=3385#msg_12752
I can confirm negative chi2 values from the PndKinVtxFitter as you observed them also for PandaRoot revision 12725.
Marius]]>Marius Mertens2011-10-26T11:00:48-00:00Re: negative chi2 in PndKinVtxFitter !
https://forum.gsi.de/index.phpindex.php?t=rview&goto=12759&th=3385#msg_12759
Thanks for finding that. I want to make two remarks. First, the PndKinVtxfitter is not the "final product" as everybody wants it to be. As discussed with Jim, the developer will hopefully be invited for a period to finalize it. Second, the fitter is using always a list of true reconstruction tracks, not intermediate composite particles. Hence it tries to fit your K+K-pi+ to one vertex.
Greetings.
Ralf
]]>Ralf Kliemt2011-10-31T08:15:23-00:00Re: negative chi2 in PndKinVtxFitter !
https://forum.gsi.de/index.phpindex.php?t=rview&goto=12776&th=3385#msg_12776
do you happen to know if there is any fitter which is able to work with composite particles?
Marius]]>Marius Mertens2011-11-14T15:16:02-00:00Re: negative chi2 in PndKinVtxFitter !
https://forum.gsi.de/index.phpindex.php?t=rview&goto=12777&th=3385#msg_12777
Until now, I fear you have to artificially create a FairTrackPar for composite TCandidate's. Then the fitter will process it as a real measured track.
I think we might want to discuss how to handle such things on the Collaboration meeting. I think the handling is a bit messy because of the leftovers from Beta/Rho which are not used by Panda.
Kind regards.
Ralf]]>Ralf Kliemt2011-11-14T15:38:11-00:00Re: negative chi2 in PndKinVtxFitter !
https://forum.gsi.de/index.phpindex.php?t=rview&goto=12778&th=3385#msg_12778
thanks for the hint about manually creating the Track. However, it's not quite clear to me how this can then be converted to a TCandidate, or any other structure I can pass to the fitters. Is there any documentation on how to do that or something which can help to figure out the data flow?
Marius]]>Marius Mertens2011-11-14T18:11:06-00:00Re: negative chi2 in PndKinVtxFitter !
https://forum.gsi.de/index.phpindex.php?t=rview&goto=12779&th=3385#msg_12779
Skippong sercasm, I can only offer you to look it up in the next days.
Ralf]]>Ralf Kliemt2011-11-14T21:54:31-00:00Re: negative chi2 in PndKinVtxFitter !
https://forum.gsi.de/index.phpindex.php?t=rview&goto=12780&th=3385#msg_12780
Well, there was some nonzero chance for a tutorial or similarly helpful stuff to have escaped my attention. But I do get your point regarding documentation. I must admit that I got pretty lost when I tried to figure out how to transfer the information from object to object which convinces the fitter that the TCandidate it sees is a real (well, fake in this case) track which it can fit.