GSI Forum
GSI Helmholtzzentrum für Schwerionenforschung

Home » PANDA » PandaRoot » Analysis » Try with no PID (Analysis with no PID problems/questions)
Re: Try with no PID [message #21548 is a reply to message #21547] Wed, 27 September 2017 11:45 Go to previous messageGo to previous message
Klaus Götzen is currently offline  Klaus Götzen
Messages: 293
Registered: June 2006
Location: GSI
first-grade participant
From: *
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: 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

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

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

Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Previous Topic: Fitters and hyperons
Next Topic: Negative lambda masses after VtxFit
Goto Forum:

Current Time: Wed Jan 16 20:10:54 CET 2019

Total time taken to generate the page: 0.02924 seconds