GSI Forum
GSI Helmholtzzentrum für Schwerionenforschung

Home » PANDA » PandaRoot » Analysis » Decay analysis (question about number of particles obtained as a result, and also MCTruth stuff)
Decay analysis [message #15959] Thu, 13 March 2014 18:38 Go to next message
Alexandros is currently offline  Alexandros
Messages: 51
Registered: January 2014
Location: KVI
continuous participant
From: *
Hi all,

So I am new to panda and I am trying to perform some simulations..
I am interested in the decay psi(4160)->D0_antiD0 and D0 goes to K- π+ and antiD0 goes to K+ π-...
I have tried to follow the example rho analysis to have my first results..
Some of the results are ok but I have several problems..
First: why is the efficiency for psi(4160) so low???less than 20%...Is it normal or do i have mistakes in my code??
Second: for each event I don't get one D0 and one antid0...there are some(not so few) events that give 2 or 3 D0's and some other events 2 or 3 antiD0's...I suppose this is not normal right???I should get one D0 per event..
Third: How does this McTruth work exactly???I have tried some things for D0's, antiD0's and psi(4160) but nothing seems to work...I also tried to follow the way that the mctruth is shown in the slides from Klaus in PANDA Computing Workshop 2012 Torino but it didnt work either..
Can you tell me what is going wrong???
I will also attach my code to have a look...
For the simulation I use the tut_sim.C where I put in the evtgendirect my decay file as input...
The decay file is called psi4160.dec...
The analysis file is called psi4160......Analysis.C and the drawin part is called psi4160.........DrawPlots.C

I am looking forward for your answers....
Thank you guys and girls for your time
Re: Decay analysis [message #15960 is a reply to message #15959] Fri, 14 March 2014 08:00 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: *
Hi Alexandros,

from a quick view I would say your analysis in general looks fine. Concerning the issues you observe I have some ideas:

- Low psi efficiency: The reason for this could be your PID selection (KaonLoose, PionLoose) requiring P>0.25. There is currently a bug (being fixed at the moment) in the assignment of DIRC PID info to tracks leading to a drop in efficiency for selection of pions, kaons and protons. Try with KaonAll and PionAll and look, whether your efficiency increases. You can tighten your PID requirements in the analysis of your ntuple to find out whether this is the problem.

- Multiple MCT matches: The full mct tree match requires, that the assignement of reco particles to truth objects is correct. This is not always the case as far as I've seen, so that sometimes multiple recos are assigned to the same McTruth object. This of course also leads to multiple tree matches. I think the truth match itself works fine.

- Concerning your analysis script: I saw tiny problems in lines 149ff and 240ff, where you store info about the truth candidates, if the match was successful. When using the RhoTuple it is more safe to always fill all columns, and not some just under certain conditions, since in the latter case for some ntuple entries the branches might not be filled properly. So I usually do some 'fake' fill in case the mct match failed...

Hope that helps a bit.

Re: Decay analysis [message #15997 is a reply to message #15960] Mon, 17 March 2014 11:41 Go to previous message
Alexandros is currently offline  Alexandros
Messages: 51
Registered: January 2014
Location: KVI
continuous participant
From: *
thanks a lot Klaus...
everything seems to work now...
Previous Topic: Talk at Hirschegg 2014 - Perspectives of Open Charm Physics at PANDA
Next Topic: invariant mass
Goto Forum:

Current Time: Tue Jan 18 02:47:11 CET 2022

Total time taken to generate the page: 0.02189 seconds