GSI Forum - RDF feed
https://forum.gsi.de/index.php
[FIXED] 4Cfit
https://forum.gsi.de/index.phpindex.php?t=rview&goto=16590&th=4277#msg_16590
I am looking in the channel psi(4160)->D0-antiD0 etc.
I am trying to apply 4C fit in psi(4160).
So what I add in my macro is:
PndKinFitter fitter(psi4160[j]);
fitter.Add4MomConstraint(ini); // set 4 constraint
fitter.Fit();
RhoCandidate *d0fit = psi4160[j]->Daughter(0)->GetFit(); // get fitted D0
RhoCandidate *antid0fit = psi4160[j]->Daughter(1)->GetFit(); // get fitted antiD0
Error in <TDecompLU::DecomposeLUCrout>: matrix is singular
Error in <TDecompLU::InvertLU>: matrix is singular, 0 diag elements < tolerance of 2.2204e-16
It is the same going over and over again.
Is there something wrong???
]]>Alexandros2014-05-13T12:20:12-00:00Re: 4Cfit
https://forum.gsi.de/index.phpindex.php?t=rview&goto=16591&th=4277#msg_16591
for (j=0;j<psi4160.GetLength();++j)
{
SORRY!!!]]>Alexandros2014-05-13T12:25:44-00:00Re: 4Cfit
https://forum.gsi.de/index.phpindex.php?t=rview&goto=16592&th=4277#msg_16592
StefanoSpataro2014-05-13T12:42:47-00:00Re: 4Cfit
https://forum.gsi.de/index.phpindex.php?t=rview&goto=16593&th=4277#msg_16593
Alexandros2014-05-13T12:46:04-00:00Re: 4Cfit
https://forum.gsi.de/index.phpindex.php?t=rview&goto=16594&th=4277#msg_16594
StefanoSpataro2014-05-13T12:47:41-00:00Re: 4Cfit
https://forum.gsi.de/index.phpindex.php?t=rview&goto=16595&th=4277#msg_16595
double m0_p = TDatabasePDG::Instance()->GetParticle("proton")->Mass(); // Get nominal PDG mass of the proton
TLorentzVector ini(0, 0, pbarmom, sqrt(m0_p*m0_p + pbarmom*pbarmom) + m0_p);]]>Alexandros2014-05-13T12:50:11-00:00Re: 4Cfit
https://forum.gsi.de/index.phpindex.php?t=rview&goto=16596&th=4277#msg_16596
StefanoSpataro2014-05-13T12:53:31-00:00Re: 4Cfit
https://forum.gsi.de/index.phpindex.php?t=rview&goto=16597&th=4277#msg_16597
if this appears only rarely, this is not a problem - for some combinations the fit fails.
Best,
Klaus]]>Klaus Götzen2014-05-13T13:13:16-00:00Re: 4Cfit
https://forum.gsi.de/index.phpindex.php?t=rview&goto=16598&th=4277#msg_16598
But from the EvtGen I get p=8.352 and energy=9.342...
Do I have to put these numbers in my initial 4vector??]]>Alexandros2014-05-13T13:23:58-00:00Re: 4Cfit
https://forum.gsi.de/index.phpindex.php?t=rview&goto=16599&th=4277#msg_16599
Maybe it is better that you upload your sim + dec files.]]>StefanoSpataro2014-05-13T13:28:51-00:00Re: 4Cfit
https://forum.gsi.de/index.phpindex.php?t=rview&goto=16600&th=4277#msg_16600
simfast("psi4160","psi4160.dec",8.00000,10000,"psi(4160)")
I have ran everything before without the 4-5 new lines for the 4Cfit and I get nice results and everything works fine.
Now that I added in the for loop for psi4160 the 4c fit:
it gives me always the error I mentioned before but without crashing!!!!
So I get a column called "psi4160fit" but the result is exactly the same like the normal psi4160m!!!]]>Alexandros2014-05-13T13:57:17-00:00Re: 4Cfit
https://forum.gsi.de/index.phpindex.php?t=rview&goto=16601&th=4277#msg_16601
( p=8.352 and energy=9.342). Then most probably 8 GeV/c are not enough to form such a state.
The alternative is to start from a pbarpSystem (in both sim and dec), and in such a case your 8 GeV/c will be set. But i think that with this low momentum evtgen will now have the psi(4160) mass.]]>StefanoSpataro2014-05-13T14:01:55-00:00Re: 4Cfit
https://forum.gsi.de/index.phpindex.php?t=rview&goto=16603&th=4277#msg_16603
But still I have the same kind of error for the 4C fitter... ]]>Alexandros2014-05-13T15:12:30-00:00Re: 4Cfit
https://forum.gsi.de/index.phpindex.php?t=rview&goto=16604&th=4277#msg_16604
StefanoSpataro2014-05-13T15:14:53-00:00Re: 4Cfit
https://forum.gsi.de/index.phpindex.php?t=rview&goto=16605&th=4277#msg_16605
Alexandros2014-05-13T15:23:46-00:00Re: 4Cfit
https://forum.gsi.de/index.phpindex.php?t=rview&goto=16606&th=4277#msg_16606
StefanoSpataro2014-05-13T15:26:13-00:00Re: 4Cfit
https://forum.gsi.de/index.phpindex.php?t=rview&goto=16607&th=4277#msg_16607
Alexandros2014-05-13T15:44:00-00:00Re: 4Cfit
https://forum.gsi.de/index.phpindex.php?t=rview&goto=16608&th=4277#msg_16608
StefanoSpataro2014-05-13T15:45:34-00:00Re: 4Cfit
https://forum.gsi.de/index.phpindex.php?t=rview&goto=16609&th=4277#msg_16609
Alexandros2014-05-13T15:54:14-00:00Re: 4Cfit
https://forum.gsi.de/index.phpindex.php?t=rview&goto=16611&th=4277#msg_16611
Are you sure that you don't have some broken file? Or maybe a screwed installation?]]>StefanoSpataro2014-05-13T19:10:39-00:00Re: 4Cfit
https://forum.gsi.de/index.phpindex.php?t=rview&goto=16613&th=4277#msg_16613
I will try tomorrow and let you know so that you can close the topic.
Thanks for your help.]]>Alexandros2014-05-13T21:37:48-00:00Re: 4Cfit
https://forum.gsi.de/index.phpindex.php?t=rview&goto=16624&th=4277#msg_16624
i was using march14 version and not the scrut14...
with scrut 14 4Cfit works...with march14 not..
you can close the topic!!
thanks again...]]>Alexandros2014-05-14T12:25:00-00:00