Hi Kerstin,
> We are trying to run over your signal MC and write out all events (also
> the events without Breco candidate with the aim to get unbiased
> information on mxhadgen for the unfolding of the m_X spectrum). Following
> the instructions of
> http://www.slac.stanford.edu/BFROOT/www/Physics/Analysis/AWG/InclusiveSL/islrecoil/islbrecoiluser.instructions.html
> but using analysis-13b instead of analysis-13 and replacing Linux2 by
> Linux24 in the following, the only other change we did was in IbuIsl.tcl:
> instead of
have you compared the log files?
Alessio, can you check if you used the same tags and options they are
using to run IslBrecoilUser? in particular for PID killing?
> To run anaQA on our file we needed to make two changes to it:
> as a default, fBRecoFlavor is =0 (now) and for events with
> dnB0Lund[indexbestB]==0 (all n for charged, n==1 for neutral B)
> fBrecoFlavor does not get changed (in util.icc) and
> tmppurB =
> brecosig[modeB-10000]/(brecosig[modeB-10000]+brecobkg[modeB-10000]);
> does not get set for modeB<10000 (in recoilNtp.cc).
> We checked that this and the original anaQA give the same results for the
> variables quoted below on your file.
I did not get this but if you have the same result it's fine?
> after running anaQA on these files:
> # events from mxhadgen 768 816
> # events w/ pcms>1. 23 18
> # events w/ pcms>0. 29 25
>
> For the # of events from mxhadgen we required brecoflav==0 here (without
> this requirement we get 9218 for our file).
which are the option you used to run anaQA? (-MC -d ? -pk?) are you using
the HEAD version of VubAnalysis?
> Since we did not understand that, we made a dummy check:
>
> We try to run and just write out only the events with Breco
> candidate and compare with your root-tuples, we already find differences
this is already the case... but I am probably missing something here
> For example:
> your file our file
> # events in file 849 896
> # B0 (from modeB0) 574 623
> # chB (from modeChB) 539 579
>
> after running anaQA on these files:
> # events from mxhadgen 768 816
> # events w/ pcms>1. 23 18
> # events w/ pcms>0. 29 25
I see the same numbers that you have shown before. Can you explain again
which is the difference between the two cases?
Daniele
|