GSI Forum
GSI Helmholtzzentrum für Schwerionenforschung

Home » PANDA » PandaRoot » Analysis » McTruthMatch problems
McTruthMatch problems [message #16775] Thu, 05 June 2014 17:41 Go to previous message
Elisabetta Prencipe (2) is currently offline  Elisabetta Prencipe (2)
Messages: 214
Registered: February 2013
first-grade participant
From: *
Dear all,

I have a problem with the McTruth match. As preface, I have to say that it works pretty well in the standard trunk (25180) which I am currently using.
I substituted genfit with genfit2. If I try to run standard macros (sim-, digi-, reco-, pid-) evertyhing looks working smooth: variables are filled.

I run a standard analysis macro, making use of rho classes (e.g. ana_complete.C), and I modified it following the example /macro/scrut/ana_jpsi.C.
ntp1 (reco quantities) is filled properly.
nmc is filled properly.

The problems occur when I try to loop over the rho candidates, and ask the Mc true values in ntp1. For example:

RhoCandidate *truth = jpsi[j]->GetMcTruth();
TLorentzVector lv;
if (truth) lv = truth->P4();
qa.qaP4("true_jpsi_", lv, ntp1);


all variables labeled as true_jpsi_* look spikes centered in zero (but same number of entries as for reco variables).
Indeed, in the mc block:

theAnalysis->FillList(mclist, "McTruth");

nmc->Column("ev", (Int_t) i);
qa.qaMcList("", mclist, nmc);

variables are filled properly.
When running the standard test macros, I see often the message:

-W- PndMCTrackAssociator::Exec: Not identified Detector ID
-W- PndMCTrackAssociator::Exec: Not identified Detector ID
-W- PndMCTrackAssociator::Exec: Not identified Detector ID
-W- PndMCTrackAssociator::Exec: Not identified Detector ID

does it imply that the Mc truth association fails, here?

Even more tricky: when I try to run the fast sim, I don't get this problem at all. Of course, the fast sim is not using "kalman algorithms", except for one method in the PidCorrelator. So, I am pretty sure that I did not face properly something related to genfit2 inside pandaroot, for the MCtruth association. It is not clear to me how to debug this, and I need your help. Do you have any idea at which level the true association is going wrong? did anybody experienced before something like that?

thank you in advance for your help,

Read Message
Read Message
Read Message
Read Message
Previous Topic: hyperon fast simualtions
Next Topic: Overview of Different Methods of Running Over TTree / TNTuple
Goto Forum:

Current Time: Fri Jul 12 13:15:53 CEST 2024

Total time taken to generate the page: 0.00818 seconds