TWiki> LHCb Web>MassFitting (revision 10) Edit Raw edit Attach Print version

G-Fact Mass Fitting

Validation of signal fraction fitter (Lars)

To validate the signal fraction fitter (aka mass fitter) a series of data sets were generated with different mixes of signals and combinatoric background were generated. Each data set has approximately 1000 toy data sets with 100k events each. All data sets are located in /data/lhcb01/eklund/B2hhToyData/, in sub-directories as given below.

Event mix: Bd -> π π, Bd -> K π and combinatoric background

Directory: BdMix _561pipi_103Kpi/, 1000 data sets, 100k events each


True S.F. Initial value Mean [%]
Sigma [%]
Pull mean Pull sigma
Bd->π+π- 56.1 55 56.12 0.19 0.11+/-0.03 1.00+/-0.02
Bd->K+π- 5.665 5 5.657 0.084 -0.12+/-0.03 1.03+/-0.02
Bd->K-π+ 4.635 5 4.619 0.078 -0.21+/-0.03 0.99+/-0.02
Combinatoric 33.6 35 33.60

Event mix: Bd -> π π, Bd -> K π, Bs -> K K and combinatoric background

Directory: Mix_Bd2pipiKpiNBs2KK/ , 940 data sets, 100k events each


True S.F. Initial value Mean [%]
Sigma [%]
Pull mean Pull sigma
Bd->π+π- 24
20 23.99 0.18 -0.04+/-0.03 1.01+/-0.03
Bd->K+π- 11
10 10.98 0.11 -0.16+/-0.03 1.01+/-0.02
Bd->K-π+ 9 10 8.98 0.11 -0.21+/-0.04 1.05+/-0.03
Bs->K+K- 21
25 21.04 0.14 0.27+/-0.03 1.00+/-0.03
Combinatoric 35 35 35.02

Event mix: Bd -> π π, Bd -> K π, Bs -> K K, Bs -> K π and combinatoric background

Directory: Mix_BdNBs2pipiKpiKK/ , 860 data sets, 100k events each


True S.F. Initial value Mean [%]
Sigma [%]
Pull mean Pull sigma
Bd->π+π- 20
23 19.98 0.17 -0.12+/-0.03 0.98+/-0.03
Bd->K+π- 13.75
12 13.75 0.14 -0.01+/-0.03 0.99+/-0.02
Bd->K-π+ 11.25 12 11.24 0.13 -0.08+/-0.03 0.97+/-0.03
Bs->K+K- 15
12 15.03 0.14 0.21+/-0.03 0.99+/-0.03
Bs->K+π- 5.56 5 5.56 0.10 0.01+/-0.04 1.02+/-0.03
Bs->K-π+ 2.44 5 2.43 0.07 -0.11+/-0.03 0.98+/-0.03
Combinatoric 32 35 35.02


Event mix: Bd -> π π, Bd -> K π, Bs -> K K, Bs -> K π, Bd -> π π π and combinatoric background

Directory: fullMFMixAug2010/ , 1000 data sets, 100k events each


True S.F. Initial value Mean [%]
Sigma [%]
Pull mean Pull sigma
Bd->π+π- 8.47
10 8.45 0.13 -0.13+/-0.03 1.02+/-0.02
Bd->K+π- 17.82
15 17.84 0.14 0.12+/-0.03 0.99+/-0.03
Bd->K-π+ 14.58 15 14.58 0.13 -0.02+/-0.03 1.01+/-0.03
Bs->K+K- 8.47
10 8.47 0.10 0.01+/-0.03 1.02+/-0.03
Bs->K+π- 1.62 1 1.62 0.07 -0.03+/-0.04 0.96+/-0.03
Bs->K-π+ 0.72 1 0.71 0.05 -0.32+/-0.03 1.01+/-0.02
Bd->π+π-π0 15.0 10 14.98 0.16 -0.11+/-0.03 1.02+/-0.02
Combinatoric 33.32 38 35.02


Sensitivity to number of events

A series of fits to the data set Mix_Bd2pipiKpiNBs2KK (above) was done using a decreasing number of events, from 100 000 to 100 events. The initial values and fit configuration was identical to the fit described above. The results are summarised in the following graphs.

  • The error from the fit as reported by MINUIT versus number of events is show in this graph. As expected, the statistical error from the fit follows well the 1/sqrt(N) law. The unit here is % signal fraction, hence 1% error means that the signal fraction Bd -> pi pi is estimated to be 24 +/- 1 %.
  • The mean of the pull distribution is shown in this graph. The mean is less that 0.1 apart from the fits with the largest number of events. Since the statistical error is very small (<2*10^-3 for 100k evt) even the smallest bias is seen. At 100k event the observed absolute bias is < 4*10^-4.
  • The bias in absolute numbers is shown in this graph. This is calculated as the statistical error of the fit times the mean of the pull distribution. Hence this gives the bias in absoulte units. The bias is less than 0.1 % if more than 1000 events are used, below that number a measurable bias is seen.
  • The sigma of the pull distribution is seen in this graph. Seems to be pretty independent of the number of events.

Paul's study with the signal fraction fitting

Just start typing here paul....This is a link.

-- AlisonBates - 2010-08-06

Topic attachments
I Attachment History Action Size Date Who Comment
PDFpdf AbsBiasVsNEvents.pdf r1 manage 9.2 K 2010-09-24 - 09:06 LarsEklund  
PDFpdf FitErrorVsNEvents.pdf r1 manage 9.0 K 2010-09-23 - 11:20 LarsEklund  
PDFpdf MeanOfPullVsNEvents.pdf r1 manage 8.6 K 2010-09-24 - 09:06 LarsEklund  
PNGpng Mix1_Bd2KminusPiplus_pull.png r1 manage 21.9 K 2010-09-21 - 10:16 LarsEklund  
PNGpng Mix1_Bd2KplusPiminus_pull.png r1 manage 23.7 K 2010-09-21 - 10:16 LarsEklund  
PNGpng Mix1_Bd2PiPi_pull.png r1 manage 20.2 K 2010-09-21 - 10:16 LarsEklund  
PNGpng Mix2_Bd2KminusPiplus_pull.png r1 manage 21.5 K 2010-09-21 - 11:23 LarsEklund  
PNGpng Mix2_Bd2KplusPiminus_pull.png r1 manage 22.1 K 2010-09-21 - 11:23 LarsEklund  
PNGpng Mix2_Bd2PiPi_pull.png r1 manage 20.3 K 2010-09-21 - 11:23 LarsEklund  
PNGpng Mix2_Bs2KK_pull.png r1 manage 20.9 K 2010-09-21 - 12:41 LarsEklund  
PNGpng Mix3_Bd2KminusPiplus_pull.png r1 manage 21.9 K 2010-09-21 - 13:00 LarsEklund  
PNGpng Mix3_Bd2KplusPiminus_pull.png r1 manage 21.5 K 2010-09-21 - 13:00 LarsEklund  
PNGpng Mix3_Bd2PiPi_pull.png r1 manage 21.1 K 2010-09-21 - 13:00 LarsEklund  
PNGpng Mix3_Bs2KK_pull.png r1 manage 22.2 K 2010-09-21 - 13:01 LarsEklund  
PNGpng Mix3_Bs2KminusPiplus_pull.png r1 manage 22.7 K 2010-09-21 - 13:01 LarsEklund  
PNGpng Mix3_Bs2KplusPiminus_pull.png r1 manage 21.9 K 2010-09-21 - 13:01 LarsEklund  
PNGpng Mix4_Bd2KminusPiplus_pull.png r1 manage 22.8 K 2010-09-21 - 15:07 LarsEklund  
PNGpng Mix4_Bd2KplusPiminus_pull.png r1 manage 23.4 K 2010-09-21 - 15:07 LarsEklund  
PNGpng Mix4_Bd2PiPiPi_pull.png r1 manage 23.1 K 2010-09-21 - 15:07 LarsEklund  
PNGpng Mix4_Bd2PiPi_pull.png r1 manage 23.8 K 2010-09-21 - 15:07 LarsEklund  
PNGpng Mix4_Bs2KK_pull.png r1 manage 23.3 K 2010-09-21 - 15:11 LarsEklund  
PNGpng Mix4_Bs2KminusPiplus_pull.png r1 manage 23.0 K 2010-09-21 - 15:07 LarsEklund  
PNGpng Mix4_Bs2KplusPiminus_pull.png r1 manage 23.0 K 2010-09-21 - 15:07 LarsEklund  
PDFpdf SigmaOfPullVsNEvents.pdf r1 manage 8.7 K 2010-09-24 - 09:06 LarsEklund  
GIFgif c1.gif r1 manage 7.7 K 2010-08-16 - 14:27 AlisonBates  
PNGpng mcEventClasses.png r1 manage 20.1 K 2010-08-10 - 08:14 LarsEklund  
PNGpng mcEventRatios.png r1 manage 14.4 K 2010-08-10 - 08:14 LarsEklund  
Edit | Attach | Print version | History: r14 | r12 < r11 < r10 < r9 | Backlinks | Raw View | Raw edit | More topic actions...
Topic revision: r10 - 2010-09-24 - LarsEklund
 
This site is powered by the TWiki collaboration platform Powered by PerlCopyright © 2008-2020 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding TWiki? Send feedback