Difference: MassFitting (10 vs. 11)

Revision 112010-09-24 - LarsEklund

Line: 1 to 1

META TOPICPARENT	name="WebHome"

G-Fact Mass Fitting

Validation of signal fraction fitter (Lars)

Line: 27 to 27

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.

Changed:

<
<

Paul's study with the signal fraction fitting

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

>
>

Sensitivity to initial values of the fit

Added:

>
>

To study how sensitive the mass fitter is to the initial values given to the fit, a series of 15 fits were made to all 940 Mix_Bd2pipiKpiNBs2KK toy data sets. Four sets of completely random intitial values were used and one set of 'sensible' values for reference. And the fits were perfomed with 10000, 1000 and 100 events (5x3=15 sets of fits). The table below shos the five sets of initial values and the true values of the signal fractions.

Signal class	random 1	random 2	random 3	random 4	sensible	true value
Bd -> pi+ pi1	0.020	0.245	0.279	0.157	0.200	0.240
Bd -> K+ pi	0.428	0.231	0.125	0.254	0.150	0.110
Bd -> K- pi+	0.253	0.354	0.121	0.099	0.150	0.090
Bs -> K+ K-	0.194	0.061	0.250	0.240	0.250	0.210
Combinatoric	0.195	0.110	0.225	0.250	0.250	0.350

Added:

>
>

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. The statistical error is as expected completely independent of the initial values and depend only on the number of events.
The mean of the pull distribution is shown in this graph. The mean of the pull scatters between +/-0.1 for all signal classes and initial values, with a slightly larger spread when only using 100 events.
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 (of %). As seen from the graph, the bias in signal fractions are less than 0.1%-units if 1000 or more events are used and less than 0.5%-units if 100 events are used.
The sigma of the pull distribution is seen in this graph. Seems to be pretty independent of the number of events and initial values.

-- AlisonBates - 2010-08-06

META FILEATTACHMENT	attachment="mcEventClasses.png" attr="" comment="" date="1281428071" name="mcEventClasses.png" path="mcEventClasses.png" size="20615" stream="mcEventClasses.png" tmpFilename="/usr/tmp/CGItemp29263" user="LarsEklund" version="1"

Line: 58 to 64

META FILEATTACHMENT	attachment="Mix4_Bs2KminusPiplus_pull.png" attr="" comment="" date="1285081659" name="Mix4_Bs2KminusPiplus_pull.png" path="Mix4_Bs2KminusPiplus_pull.png" size="23586" stream="Mix4_Bs2KminusPiplus_pull.png" tmpFilename="/usr/tmp/CGItemp29127" user="LarsEklund" version="1"
META FILEATTACHMENT	attachment="Mix4_Bs2KplusPiminus_pull.png" attr="" comment="" date="1285081666" name="Mix4_Bs2KplusPiminus_pull.png" path="Mix4_Bs2KplusPiminus_pull.png" size="23580" stream="Mix4_Bs2KplusPiminus_pull.png" tmpFilename="/usr/tmp/CGItemp29140" user="LarsEklund" version="1"
META FILEATTACHMENT	attachment="Mix4_Bs2KK_pull.png" attr="" comment="" date="1285081885" name="Mix4_Bs2KK_pull.png" path="Mix4_Bs2KK_pull.png" size="23849" stream="Mix4_Bs2KK_pull.png" tmpFilename="/usr/tmp/CGItemp29124" user="LarsEklund" version="1"

Changed:

<
<

META FILEATTACHMENT	attachment="FitErrorVsNEvents.pdf" attr="" comment="" date="1285240808" name="FitErrorVsNEvents.pdf" path="FitErrorVsNEvents.pdf" size="9257" stream="FitErrorVsNEvents.pdf" tmpFilename="/usr/tmp/CGItemp27693" user="LarsEklund" version="1"
META FILEATTACHMENT	attachment="MeanOfPullVsNEvents.pdf" attr="" comment="" date="1285319161" name="MeanOfPullVsNEvents.pdf" path="MeanOfPullVsNEvents.pdf" size="8824" stream="MeanOfPullVsNEvents.pdf" tmpFilename="/usr/tmp/CGItemp27401" user="LarsEklund" version="1"
META FILEATTACHMENT	attachment="SigmaOfPullVsNEvents.pdf" attr="" comment="" date="1285319168" name="SigmaOfPullVsNEvents.pdf" path="SigmaOfPullVsNEvents.pdf" size="8956" stream="SigmaOfPullVsNEvents.pdf" tmpFilename="/usr/tmp/CGItemp27469" user="LarsEklund" version="1"
META FILEATTACHMENT	attachment="AbsBiasVsNEvents.pdf" attr="" comment="" date="1285319173" name="AbsBiasVsNEvents.pdf" path="AbsBiasVsNEvents.pdf" size="9420" stream="AbsBiasVsNEvents.pdf" tmpFilename="/usr/tmp/CGItemp27422" user="LarsEklund" version="1"

>
>

META FILEATTACHMENT	attachment="FitErrorVsNEvents.pdf" attr="" comment="" date="1285349079" name="FitErrorVsNEvents.pdf" path="AbsBiasVsNEvents.pdf" size="8992" stream="AbsBiasVsNEvents.pdf" tmpFilename="/usr/tmp/CGItemp24058" user="LarsEklund" version="2"
META FILEATTACHMENT	attachment="MeanOfPullVsNEvents.pdf" attr="" comment="" date="1285349096" name="MeanOfPullVsNEvents.pdf" path="MeanOfPullVsNEvents.pdf" size="8283" stream="MeanOfPullVsNEvents.pdf" tmpFilename="/usr/tmp/CGItemp24243" user="LarsEklund" version="2"
META FILEATTACHMENT	attachment="SigmaOfPullVsNEvents.pdf" attr="" comment="" date="1285349123" name="SigmaOfPullVsNEvents.pdf" path="SigmaOfPullVsNEvents.pdf" size="8373" stream="SigmaOfPullVsNEvents.pdf" tmpFilename="/usr/tmp/CGItemp24389" user="LarsEklund" version="2"
META FILEATTACHMENT	attachment="AbsBiasVsNEvents.pdf" attr="" comment="" date="1285349061" name="AbsBiasVsNEvents.pdf" path="AbsBiasVsNEvents.pdf" size="8992" stream="AbsBiasVsNEvents.pdf" tmpFilename="/usr/tmp/CGItemp23890" user="LarsEklund" version="2"
META FILEATTACHMENT	attachment="AbsBiasVsInitValues.pdf" attr="" comment="" date="1285350317" name="AbsBiasVsInitValues.pdf" path="AbsBiasVsInitValues.pdf" size="10698" stream="AbsBiasVsInitValues.pdf" tmpFilename="/usr/tmp/CGItemp27751" user="LarsEklund" version="1"
META FILEATTACHMENT	attachment="FitErrorVsInitValues.pdf" attr="" comment="" date="1285350324" name="FitErrorVsInitValues.pdf" path="FitErrorVsInitValues.pdf" size="10162" stream="FitErrorVsInitValues.pdf" tmpFilename="/usr/tmp/CGItemp27758" user="LarsEklund" version="1"
META FILEATTACHMENT	attachment="MeanOfPullVsInitValues.pdf" attr="" comment="" date="1285350336" name="MeanOfPullVsInitValues.pdf" path="MeanOfPullVsInitValues.pdf" size="9709" stream="MeanOfPullVsInitValues.pdf" tmpFilename="/usr/tmp/CGItemp27717" user="LarsEklund" version="1"
META FILEATTACHMENT	attachment="SigmaOfPullVsInitValues.pdf" attr="" comment="" date="1285350343" name="SigmaOfPullVsInitValues.pdf" path="SigmaOfPullVsInitValues.pdf" size="9976" stream="SigmaOfPullVsInitValues.pdf" tmpFilename="/usr/tmp/CGItemp27775" user="LarsEklund" version="1"

View topic | History: r14 < r13 < r12 < r11 | More topic actions...