(r4) AsymmetryFitting < LHCb

LHCb Web>AsymmetryFitting (revision 4) (raw view)~~EditAttach~~

---+ B2hh GFact Asymmetry Fitting page

The aims for the asymmetry fitter are the following
   * Perform basic checks of current fitting package with toy data for Bd and Bs asymmetries (see below).
   * Add proper time pdf to the fitting package (started looking at this below).
   * Add Bs&rarr;KK pdf to fit package
   * Think about adding the other background pdfs (3 body, &Lambda;'s)
   * Alter the parameters that are being varied by the minimiser to d, &theta; and &gamma;
   * Move to MC data/ real data (Bd) for first looks. So far 0.3fb-1 data is enough for asymmetry fitter

Most of the following is on my web page http://ppewww.physics.gla.ac.uk/~alison/B2hh.html where all the attached files are there and work. I have now stopped using my web page and will now just post all my results here.
---++ Basic Checks with asymmetry fitter
---+++ 1. B<sub>d</sub> &rarr; &pi; &pi;

50 datasets were generated with 6099 events per data set using 0.68 b2pipi (sf1=1 and ssb = 0.68) and 0.32 combinatoric background (MassDistribution) mass distribution. The asymmetry fitter was used to fit for Bd asymmetries only with sf1'=0.68 and obtained the following:<br /> Adir Bd - input was *0.38* and fitted was *0.3953 ± 0.0476* ( [[%ATTACHURL%/AdirBdforjustBdandCB.gif][A<sup>dir</sup><sub>Bd</sub>]])<br /> Amix Bd - input was *0.61* and fitted was <strong>0.6503 ± 0.0311 </strong>( [[%ATTACHURL%/AmixBdforjustBdandCB.gif][A<sup>mix</sup><sub>Bd</sub>]])<br /> Then re-ran same job with 150 data sets with 2000 events per data set and fitted for the following:<br /> Adir Bd - input was *0.38* and fitted was <strong> 0.4338 ± 0.0196 </strong>( [[%ATTACHURL%/ABddirjustBdppipiCB.gif][A<sup>dir</sup><sub>Bd</sub>]])<br /> Amix Bd - input was *0.61* and fitted was <strong> 0.6619 ± 0.0197 </strong>( [[%ATTACHURL%/ABdmixjustBdppipiCB.gif][A<sup>mix</sup><sub>Bd</sub>]]) <br /> The asymmetry plot is shown [[%ATTACHURL%/asymplotBdCB.gif][ here ]].

******100 data sets with 10k events per data set using 0.68 b2pipi (sf1=1 and ssb = 0.68) and 0.32 combinatoric background and mass window of 5229 - 5329 MeV, I fit the following asymmetries:

Adir Bd - input was *0.38* and fitted was <strong> 0.4078 ± 0.0077 </strong>( [[%ATTACHURL%/justbdandcbnarrowmassAdir.gif][A<sup>dir</sup><sub>Bd</sub>]])<br /> Amix Bd - input was *0.61* and fitted was *0.599 ± 0.006* ( [[%ATTACHURL%/justbdandcbnarrowmassAmix.gif][A<sup>mix</sup><sub>Bd</sub>]])

---+++ 2. B<sub>s</sub> &rarr; K K

50 datasets were generated with 6099 events per data set using 0.68 bs2KK (sf3=1 and ssb = 0.68) and 0.32 combinatoric background ( [[%ATTACHURL%/bs2KKdatagenwithcombkgnd.gif][ mass distribution ]]). The asymmetry fitter was used to fit for Bs asymmetries only with sf3'=0.68 and obtained the following [[%ATTACHURL%/AdirforBsandCB.gif][Adir for Bd]] and [[%ATTACHURL%/AmixforBsandCB.gif][Amix for Bs]] which in numbers are:<br /> Adir Bs - input was *0.1* and fitted was <strong> 0.0811 ± 0.0275 </strong><br /> Amix Bs - input was *0.25* and fitted was *0.1799 ± 0.0346*

Ran a job with 100 data sets with 6000 events per data set using 0.664 Bs2KK and 0.336 combinatoric background.

Adir Bs - input was *0.1* and fitted was <strong> 0.06436 ± 0.00630 </strong>( [[%ATTACHURL%/AdirBs2KKCBwithsmearing.gif][A<sup>dir</sup><sub>Bs</sub>]])<br /> Amix Bs - input was *0.25* and fitted was <strong> 0.1721 ± 0.0076 </strong>( [[%ATTACHURL%/AmixBs2KKCBwithsmearing.gif][A<sup>mix</sup><sub>Bs</sub>]])

Now ran the same job as above but removed the Gaussian smearing of the proper time in the toy data generation and re-ran the job and got the exact same asymmetries and uncertainties. This should make a difference surely? Still looking in to it.
---+++ 3. B<sub>d</sub> &rarr; &pi; &pi; and B<sub>d</sub> &rarr; K &pi;

50 data sets were generated with 6099 events per data set using 0.68 bd2pipi and bd2Kpi (sf1=sf2=0.5 and ssb = 0.68) and 0.32 combinatoric background ( [[%ATTACHURL%/bothBddatagenwithCBforMassdistn.gif][ mass distribution ]]).The asymmetry fitter was used to fit for Bd asymmetries only with sf1'=sf2'=0.34 and obtained the following [[%ATTACHURL%/bothBdwithCBAdir.gif][Adir for Bd]] and [[%ATTACHURL%/bothBdwithCBAmix.gif][Amix for Bd]] which in numbers are:<br /> Adir Bd - input was *0.38* and fitted was <strong> 0.298 ± 0.206 </strong><br /> Amix Bd - input was *0.61* and fitted was *0.660 ± 0.062*

Now using 150 data sets with 100k events per data set using bd2pipi and bd2Kpi (sf1= 0.845 and sf2=0.155 and ssb = 0.664, see Laurence's thesis, page 183). The asymmetry fitter was used to fit for Bd asymmetries with sf1'=0.561 and sf2'=0.103 and obtained the following:<br /> Adir Bd - input was *0.38* and fitted was *0.4241 ± 0.0017* ( [[%ATTACHURL%/AdirB2pipiBd2pipiBd2kpiCBalaL.gif][A<sup>dir</sup><sub>Bd</sub>]]) and Laurence got <strong> 0.407 ± 0.033 </strong><br /> Amix Bd - input was *0.61* and fitted was <strong> 0.6818 ± 0.0018 </strong>( [[%ATTACHURL%/AmixB2pipiBd2pipiBd2kpiCBalaL.gif][A<sup>mix</sup><sub>Bd</sub>]]) and Laurence got <strong> 0.655 ± 0.024 </strong><br /> The asymmetry plot is shown [[%ATTACHURL%/AsymplotLaurencecfjobforBd.gif][ here ]]. The pull plots show a large bias. The only difference was the large mass window I used (4800 - 5800 MeV) compared with Laurence (5229 - 5329 MeV). Hence I have re-run the same signal fractions in a 5fb-1 job (288 data sets with 100k events per data set) with the much tighter mass window and found the following:

Adir Bd - input was *0.38* and fitted was *0.4245 ± 0.0022* ( [[%ATTACHURL%/AdirBigJob2BddecaysCB.gif][A<sup>dir</sup><sub>Bd</sub>]]) and Laurence got *0.407 ± 0.033*

Amix Bd - input was *0.61* and fitted was <strong> 0.6811 ± 0.0019 </strong>( [[%ATTACHURL%/AmixBigJob2BddecaysCB.gif][A<sup>mix</sup><sub>Bd</sub>]]) and Laurence got *0.655 ± 0.024*

The asymmetry plot is shown [[%ATTACHURL%/AsymPLotBigJob2BddecaysCB.gif][ here ]]. However the pull plots are still showing a large bias in the mean of the pull plot (the widths are approximately ok). What's causing the bias in the pull plots? It is not present in the first case when only B2pipi and combinatoric background are used....looking in to this now. Interesteingly, in an earlier version of Laurence's thesis he has numbers like mine then they suddenly got really good. Waiting to hear back from him to see if he can remember what he changed.
---+++ 4. B<sub>d</sub> &rarr; &pi; &pi;, B<sub>d</sub> &rarr; K &pi; and B<sub>s</sub> &rarr; K K

50 datasets were generated with 6099 events per data set using all 3 channels with sf1=sf2=sf3=0.333 and ssb = 0.68) and 0.32 combinatoric background ( [[%ATTACHURL%/alldatagenandCBmass.gif][ mass distribution ]]). The asymmetry fitter was used to fit for Bd asymmetries only with sf1'=sf2'=sf3'=0.226 and obtained the following [[%ATTACHURL%/AdirBdforall3ofthem.gif][Adir for Bd]] and [[%ATTACHURL%/AmixBdforall3ofthem.gif][Amix for Bd]] which in numbers are:<br /> Adir Bd - input was *0.38* and fitted was <strong>0.4917 ± 0.0638 </strong><br /> Amix Bd - input was *0.61* and fitted was <strong>0.6852 ± 0.0456 </strong><br /> Then re-ran the asymmetry fitter for the B<sub>s</sub> asymmetries and found;<br /> Adir Bs - input was *0.1* and fitted was <strong> 0.0416 ± 0.0606 </strong>( [[%ATTACHURL%/AdirBsforall3ofthem.gif][A<sup>dir</sup><sub>Bs</sub>]])<br /> Amix Bs - input was *0.25* and fitted was <strong> 0.1505 ± 0.0598</strong>( [[%ATTACHURL%/AmixBsforall3ofthem.gif][A<sup>mix</sup><sub>Bs</sub>]])
---+++ 5. Full toy data sample

-- Main.AlisonBates - 2010-08-06

Attachments

Topic attachments
I	Attachment	History	Action	Size	Date	Who
gif	ABddirjustBdppipiCB.gif	r1	manage	14.9 K	2010-08-06 - 15:00	AlisonBates
gif	ABdmixjustBdppipiCB.gif	r1	manage	14.7 K	2010-08-06 - 15:00	AlisonBates
gif	AdirB2pipiBd2pipiBd2kpiCBalaL.gif	r1	manage	14.2 K	2010-08-06 - 15:00	AlisonBates
gif	AdirBd.gif	r1	manage	14.8 K	2010-08-06 - 15:03	AlisonBates
gif	AdirBdforall3ofthem.gif	r1	manage	15.1 K	2010-08-06 - 15:01	AlisonBates
gif	AdirBdforjustBdandCB.gif	r1	manage	15.7 K	2010-08-06 - 15:03	AlisonBates
gif	AdirBigJob2BddecaysCB.gif	r1	manage	14.2 K	2010-08-09 - 08:42	AlisonBates
gif	AdirBs.gif	r1	manage	13.6 K	2010-08-06 - 15:04	AlisonBates
gif	AdirBs2KKCBwithsmearing.gif	r1	manage	14.0 K	2010-08-09 - 11:32	AlisonBates
gif	AdirBsforall3ofthem.gif	r1	manage	14.3 K	2010-08-06 - 15:12	AlisonBates
gif	AdirforBsandCB.gif	r1	manage	15.6 K	2010-08-06 - 15:13	AlisonBates
gif	AmixB2pipiBd2pipiBd2kpiCBalaL.gif	r1	manage	14.1 K	2010-08-06 - 15:02	AlisonBates
gif	AmixBdforall3ofthem.gif	r1	manage	15.1 K	2010-08-06 - 15:03	AlisonBates
gif	AmixBdforjustBdandCB.gif	r1	manage	15.1 K	2010-08-06 - 15:11	AlisonBates
gif	AmixBigJob2BddecaysCB.gif	r1	manage	14.1 K	2010-08-09 - 08:43	AlisonBates
gif	AmixBs2KKCBwithsmearing.gif	r1	manage	14.1 K	2010-08-09 - 11:32	AlisonBates
gif	AmixBsforall3ofthem.gif	r1	manage	15.5 K	2010-08-06 - 15:04	AlisonBates
gif	AmixforBsandCB.gif	r1	manage	14.6 K	2010-08-06 - 15:05	AlisonBates
gif	AsymPLotBigJob2BddecaysCB.gif	r1	manage	10.0 K	2010-08-09 - 08:43	AlisonBates
gif	AsymplotLaurencecfjobforBd.gif	r1	manage	17.1 K	2010-08-06 - 15:06	AlisonBates
gif	alldatagenandCBmass.gif	r1	manage	7.7 K	2010-08-06 - 15:02	AlisonBates
gif	asymplotBdCB.gif	r1	manage	19.6 K	2010-08-06 - 15:05	AlisonBates
gif	b2pipidatagenwithcombkgnd.gif	r1	manage	8.0 K	2010-08-06 - 15:06	AlisonBates
gif	bothBddatagenwithCBforMassdistn.gif	r1	manage	7.7 K	2010-08-06 - 15:07	AlisonBates
gif	bothBdwithCBAdir.gif	r1	manage	14.9 K	2010-08-06 - 15:08	AlisonBates
gif	bothBdwithCBAmix.gif	r1	manage	14.7 K	2010-08-06 - 15:08	AlisonBates
gif	bs2KKdatagenwithcombkgnd.gif	r1	manage	7.8 K	2010-08-06 - 15:08	AlisonBates
gif	justbdandcbnarrowmassAdir.gif	r1	manage	14.8 K	2010-08-16 - 14:59	AlisonBates
gif	justbdandcbnarrowmassAmix.gif	r1	manage	15.1 K	2010-08-16 - 14:59	AlisonBates

Topic revision: r4 - 2010-08-16 - AlisonBates

LHCb