GSI Forum
GSI Helmholtzzentrum für Schwerionenforschung

Home » PANDA » PandaRoot » Fast Simulations » Mc Truth Match
Mc Truth Match [message #16669] Tue, 20 May 2014 13:57 Go to next message
Elisa Fioravanti is currently offline  Elisa Fioravanti
Messages: 84
Registered: January 2008
continuous participant
From: *fe.infn.it
Dear all,
I'm simulating the channel
J/psi pi0 pi0 at the energy in the center of mass of the X(3872).
I select j/psi->e+e- in this way:

TString pidalg = "PidChargedProbability";
theAnalysis->FillList(eminus, "ElectronLooseMinus", pidalg);
theAnalysis->FillList(eplus, "ElectronTightPlus", pidalg);
SelectTruePid(theAnalysis, eplus);
SelectTruePid(theAnalysis, eminus);
jpsiee.Combine(eminus,eplus);
for(Int_t j=0;j<jpsiee.GetLength();++j){
if (theAnalysis->McTruthMatch(jpsiee[j])){
ntp1->Column("jpsi_before_fit",(Float_t) jpsiee[j]->M());
ntp1->DumpData();
}
}
in this way only 20% of J/psi are reconstructed.
If I did not (theAnalysis->McTruthMatch(jpsiee[j]))
I have 66% of J/psi reconstructed.

Maybe something is going wrong with the McTruthMatch?
Do you have an idea?

Thanks a lot
Elisa
Re: Mc Truth Match [message #16674 is a reply to message #16669] Tue, 20 May 2014 16:21 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 Elisa,


if you have enabled 'Photos' in EvtGen, the truth match fails if you don't add the bremsstrahlung gamma to the J/psi reco candidate. I.e. use 'noPhotos' in the decay file and try again, whether the number of truth matches increases.


Best,
Klaus
Re: Mc Truth Match [message #17282 is a reply to message #16674] Fri, 12 September 2014 10:31 Go to previous message
StefanoSpataro is currently offline  StefanoSpataro
Messages: 2736
Registered: June 2005
Location: Torino
first-grade participant

From: *lnf.infn.it
Did it help?
Previous Topic: Running analysis task in simfast macro
Next Topic: Switch off MCTrack persistency
Goto Forum:
  


Current Time: Sat Oct 12 12:31:19 CEST 2024

Total time taken to generate the page: 0.00622 seconds