Mc Truth Match [message #16669] |
Tue, 20 May 2014 13:57 |
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
|
|
|
|
|