Print

---

hier. wir koennen NICHTS reproduzieren von denen. NICHTS. nicht mal, wenn
es eigentlich identisch sein sollte.


On Fri, 7 Nov 2003, Kerstin Tackmann wrote:

>
> Hello experts,
>
> we need your help, please read this mail until the end.
>
> 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
> # -- Dump semileptonic events on generator level
> # mod talk CompBRecoNtpDump
> # mod talk CompIbuNtpDump
> #  verbose set f
> #  dumpAllSemiLep set 1
> #  show
> # exit
> we used
> # -- Dump semileptonic events on generator level
> # mod talk CompBRecoNtpDump
>  mod talk CompIbuNtpDump
> #  verbose set f
>   dumpAllSemiLep set 1
>   show
>  exit
>
>
> We found some inconsistencies when comparing the file for which the
> following tclfile was used:
> /afs/slac.stanford.edu/g/babar/www/Physics/Analysis/AWG/InclusiveSL/islrecoil/tcl/new2vubmix-2000-new2000-12.tcl
>
> with the following of your files:
> /nfs/farm/babar/AWG23/ISL/sx-080702/newsig/rootfiles/new2vubmix-2000-new2000-12.root
>
> For example:
>                                your file       our file
> # events in file                  849             896
> # B0 (from modeB0)                574             623
> # chB (from modeChB)              539             579
>
> 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.
>
> 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).
>
> So our file has more events and more B candidates but less leptons in
> these regions of pcms.
>
>
> 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
> in the number of events that were written when we run over the same
> collection.
>
> We used IbuIsl.tcl with
> # -- Dump semileptonic events on generator level
> # mod talk CompBRecoNtpDump
> # mod talk CompIbuNtpDump
> #  verbose set f
> #  dumpAllSemiLep set 1
> #  show
> # exit
>
> Here is what we did:
> We followed the instructions of
> http://www.slac.stanford.edu/BFROOT/www/Physics/Analysis/AWG/InclusiveSL/islrecoil/islbrecoiluser.instructions.html
> but took analysis-13b instead of analysis-13 (and replaced Linux2 by
> Linux24 in the following instructions as well as RH62 by RH72) and
> compiled locally. We used the run script for the submission of a job with
> the following tcl file:
> /afs/slac.stanford.edu/g/babar/www/Physics/Analysis/AWG/InclusiveSL/islrecoil/tcl/new2vubmix-2000-new2000-12.tcl
>
> with the following command:
> ../IslBrecoilUser/run -q bfobjy -l $SCRATCH/log -r
> $SCRATCH/output -w 20 ../bin/Linux24/IslBrecoilUserApp
> ../IslBrecoilUser/islNew.tcl
> tclpath/signconc/new2vubmix-2000-new2000-12.tcl
>
> after h2r we compared the resulting rootfile with the following of your
> rootfiles:
> /nfs/farm/babar/AWG23/ISL/sx-080702/newsig/rootfiles/new2vubmix-2000-new2000-12.root
>
> and we find some discrepancies.
> 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
>
> So our file, which should be identical to your original one, has more
> events and more B candidates but less leptons in these regions of pcms.
>
> Is there a way to understand these discrepancies? Is the file we use for
> comparison the correct one to use?  Are we doing anything wrongly?
>
> Any answer will be appreciated,
> thanks,
> Kerstin
>