GSI Forum
GSI Helmholtzzentrum für Schwerionenforschung

Home » PANDA » PandaRoot » Analysis » negative chi2 in PndKinVtxFitter !
negative chi2 in PndKinVtxFitter ! [message #12751] Wed, 26 October 2011 10:21 Go to next message
Albrecht Gillitzer is currently offline  Albrecht Gillitzer
Messages: 49
Registered: March 2004
continuous participant
From: *ikp.kfa-juelich.de
Hi

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
Re: negative chi2 in PndKinVtxFitter ! [message #12752 is a reply to message #12751] Wed, 26 October 2011 13:00 Go to previous messageGo to next message
Marius Mertens is currently offline  Marius Mertens
Messages: 55
Registered: January 2009
continuous participant
From: *ikp.kfa-juelich.de
Hi Albrecht,

I can confirm negative chi2 values from the PndKinVtxFitter as you observed them also for PandaRoot revision 12725.

Marius
Re: negative chi2 in PndKinVtxFitter ! [message #12759 is a reply to message #12751] Mon, 31 October 2011 09:15 Go to previous messageGo to next message
Ralf Kliemt is currently offline  Ralf Kliemt
Messages: 507
Registered: May 2007
Location: GSI, Darmstadt
first-grade participant

From: *pool.mediaWays.net
Hello Albrecht,

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
Re: negative chi2 in PndKinVtxFitter ! [message #12776 is a reply to message #12759] Mon, 14 November 2011 16:16 Go to previous messageGo to next message
Marius Mertens is currently offline  Marius Mertens
Messages: 55
Registered: January 2009
continuous participant
From: *ikp.kfa-juelich.de
Hi Ralf,

do you happen to know if there is any fitter which is able to work with composite particles?

Marius
Re: negative chi2 in PndKinVtxFitter ! [message #12777 is a reply to message #12751] Mon, 14 November 2011 16:38 Go to previous messageGo to next message
Ralf Kliemt is currently offline  Ralf Kliemt
Messages: 507
Registered: May 2007
Location: GSI, Darmstadt
first-grade participant

From: *pool.mediaWays.net
Hi Marius,

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
Re: negative chi2 in PndKinVtxFitter ! [message #12778 is a reply to message #12777] Mon, 14 November 2011 19:11 Go to previous messageGo to next message
Marius Mertens is currently offline  Marius Mertens
Messages: 55
Registered: January 2009
continuous participant
From: *ext.kfa-juelich.de
Hi Ralf,

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
Re: negative chi2 in PndKinVtxFitter ! [message #12779 is a reply to message #12751] Mon, 14 November 2011 22:54 Go to previous messageGo to next message
Ralf Kliemt is currently offline  Ralf Kliemt
Messages: 507
Registered: May 2007
Location: GSI, Darmstadt
first-grade participant

From: *pool.mediaWays.net
Hi Marius.

Skippong sercasm, I can only offer you to look it up in the next days.

Ralf
Re: negative chi2 in PndKinVtxFitter ! [message #12780 is a reply to message #12779] Tue, 15 November 2011 00:05 Go to previous message
Marius Mertens is currently offline  Marius Mertens
Messages: 55
Registered: January 2009
continuous participant
From: *ikp.kfa-juelich.de
Hi Ralf,

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.

Thanks a lot!

Marius
Previous Topic: Beta analysis package
Next Topic: Tutorial feb2012
Goto Forum:
  


Current Time: Wed Nov 27 19:51:48 CET 2024

Total time taken to generate the page: 0.00809 seconds