Hi all,
I think I have found an inconstiency in the way we get and use the
weight for B decays in VVF.
In recoilDsys ctor we initialize the B->D BF in a vector using this
order (as by the comment)
0 = B-> X_c lnu
1 = B-> D l nu
2 = B-> D* l nu
3 = Not used
4 = D1 l nu
5 = D2* l nu
we do this for the existing measurements and for generation values.
Then we compute the weights as ratio of measurement/generation value.
The way we use these weights in VVF is via recoilDsys::weight(int bmode)
function, where bmode is actually the variable Gvxbtyp in our MC ntuples.
I've had a look in BRecoilTools/RecoilMCUtil.cc and from there it seems
to me that when Gvxbtyp=4 we have actually a B->D2* l nu decay (not a
D1 l nu decay, that is the weight picked up by recoilDsys::weight(int
bmode) function) and when Gvxbtyp=5 we have a B->D1 l nu decay.
Are there intermediate steps that I'm missing somewhere or are the two
weights for D1 and D*2 switched?
Thank you,
Antonio
|