GSI Forum
GSI Helmholtzzentrum für Schwerionenforschung

Home » PANDA » PandaRoot » Analysis » Try with no PID (Analysis with no PID problems/questions)
Try with no PID [message #21542] Wed, 27 September 2017 09:48 Go to next message
Alexandros is currently offline  Alexandros
Messages: 51
Registered: January 2014
Location: KVI
continuous participant
From: *kvi-cart.rug.nl
Hi,

This will be a long message.
I am trying to run simulation and analysis for the channel ppbar->D0D0bar->K-pi+K+pi-
I want to run this with no PID at all.
I am running this, and without using PID, i fill my list with 4 candidates(minus1, plus1, minus2, plus2) since I have 4 charged tracks.
I use the "KaonAllMinus", "PionAllPlus", "KaonAllPlus", "PionAllMinus" for my 4 candidates.
Now the problem I have is the following: I am running 35000 events with EvtGen. The possible combinations for the D0 without PID should be: K-pi+(the correct one) and K+pi-(the wrong one). There must not be other combinations since in the plot attached (Firstplot) you will see that there is a cut on the reconstructed mass at the limit of one pion and one kaon, thus no pi-pi+ and subsequently K-K+ are tried. In this way, I should get around 70000 (2 times the Kpi combinations) D0 mesons. What I get is 140000 events. So there is a double counting somewhere.
Now, I have also tried the same with "AllMinus", "AllPlus", "AllPlus", "AllMinus" for the minus1, plus1, plus2, minus2 candidates. There I would expect to have 4 times the combinations since besides the K-pi+ and the K+pi- now I also have pi-pi+ and k-k+..I can see that since the lower limit now is at the two pion masses and the number or reconstructed D0 is 4x35000 events.(Secondplot)
I have attached both plots and my analysis macro too.
Do you have any idea whether there is something wrong or whether I am doing something wrong?
Thank you in advance.
Alex
Re: Try with no PID [message #21544 is a reply to message #21542] Wed, 27 September 2017 10:14 Go to previous messageGo to next message
Klaus Götzen is currently offline  Klaus Götzen
Messages: 293
Registered: June 2006
Location: GSI
first-grade participant
From: *gsi.de
Hi Alex,


maybe there is some misunderstanding in the keywords 'KaonAllPlus', 'PionAllPlus', etc. What happens e.g. for 'KaonAllPlus' is, that the energy components of the 4-vectors of all positively charged tracks (pi+ and K+ in your case) are set with the kaon mass hypothesis to sqrt(m_K^2 + p^2). So your lists 'plus1' and 'plus2' contain exactly the same (namely all positive) tracks, in the former with pion and in the latter with kaon mass set.

For the two positive and two negative tracks per event, you get then 4 t+ t- combinations (as you already wrote). The kinematic limit of m_K + m_pi in Firstplot.png is just the result of the fact, that one track in each combination has kaon mass, the other one pion mass, so the smallest possible invariant mass is m_K + m_pi by construction.

Concerning your approach with 'AllPlus' and 'AllMinus': Here, all tracks are assigned the default mass hypo, which is pion mass. Therefore the phase space limit is 2m_pi.

I hope this explanation clarifies things a bit.


Best regards,
Klaus
Re: Try with no PID [message #21547 is a reply to message #21544] Wed, 27 September 2017 10:24 Go to previous messageGo to next message
Alexandros is currently offline  Alexandros
Messages: 51
Registered: January 2014
Location: KVI
continuous participant
From: *kvi-cart.rug.nl
Hmmmmmm.
I see what you mean. So there is a double counting effect since for each list I have the k and pi particles.
Now I have a second question.
Is there a way to remove from my ntuple this double counting effect so that i get only the K-pi+ and the K+pi- combinations for e.g the D0 reconstruction?
Thanks again.
Re: Try with no PID [message #21548 is a reply to message #21547] Wed, 27 September 2017 11:45 Go to previous messageGo to next message
Klaus Götzen is currently offline  Klaus Götzen
Messages: 293
Registered: June 2006
Location: GSI
first-grade participant
From: *gsi.de
Hi Alex,

actually it is not double counting. Double counting means (for my understanding), that you enter twice the same number (i.e. the same invariant mass) in your spectrum. This does not happen in your case. Considering all particles under all possible hypothesis is completely legitim and the only option, if an experiment does not have any means of PID.

Is it possible that you mean ideal PID instead of no PID? I.e. you only want true kaons in your plus2 list? To test that you could try

theAnalysis->FillList(plus1,  "PionLoosePlus",  "PidAlgoIdealCharged");
theAnalysis->FillList(minus1, "KaonLooseMinus", "PidAlgoIdealCharged");
theAnalysis->FillList(plus2,  "KaonLoosePlus",  "PidAlgoIdealCharged");
theAnalysis->FillList(minus2, "PionLooseMinus", "PidAlgoIdealCharged");


I guess this should work. Ideal PID assigns P=1 to the correct PID and P=0 to the others, and the Loose criterion cuts on P>0.25, thus selecting only the correct tracks. However, there might still be mis-ID even in truth matched PID (the track truth match is not perfect). Alternatively, you can apply a MC-truth match like explained here: https://panda-wiki.gsi.de/foswiki/bin/view/Computing/PandaRootRhoTutoria l#A_2.4._Monte_Carlo_Truth_Match. This takes into account the decay pattern D0 -> K- pi+ as well.

Essentially it works by e.g. adding a flag in your macro inside your D0/D0bar loop

theAnalysis->McTruthMatch(d0[j]);
nd0->Column("d0mct", (int) (d0[j]->GetMcTruth()!=0) );

and plot with
root [0] nd0->Draw("d0m","d0mct")

Best,
Klaus
Re: Try with no PID [message #21550 is a reply to message #21548] Wed, 27 September 2017 13:03 Go to previous messageGo to next message
Alexandros is currently offline  Alexandros
Messages: 51
Registered: January 2014
Location: KVI
continuous participant
From: *kvi-cart.rug.nl
The point is that on purpose I do not want to use PID at all and be able to reconstruct my mesons with the correct combination (removing all the other 3 combinations present).
I will try to figure something out.
Thanks anyway
Re: Try with no PID [message #21556 is a reply to message #21550] Wed, 27 September 2017 17:55 Go to previous message
Klaus Götzen is currently offline  Klaus Götzen
Messages: 293
Registered: June 2006
Location: GSI
first-grade participant
From: *dynamic.qsc.de
But without PID there is basically no other option than trying all combinations. What you could do as well is to reconstruct the full pbarp-> D0 D0bar system and apply a 4C fit and select the best candidate. That is something like PID on a kinematic way.
Previous Topic: Fitters and hyperons
Next Topic: Negative lambda masses after VtxFit
Goto Forum:
  


Current Time: Thu Dec 02 10:41:24 CET 2021

Total time taken to generate the page: 0.02299 seconds