Implementierung der CLs-Methode in ROOT zur statistischen ...

RooDataHist impData_2 = measure_2 . binnedClone ( ) ;
Der Import der Daten geschieht aus dem Tree des Rootfiles über ein RooDataSet in
ein RooDataHist, sodass ein RooPlot des Histogramms impData_2 identisches ergibt
wie in Abb. 4. Um dies zu erreichen sind 25 Bins nötig. Ein Blick in die Statistik des
Histogramms enthüllt 926 Einträge, so dass der Initialisierungswert von 700 Ereignissen
für die Normierung bg_yield in einem möglichen Bereich von 0 bis 1000 Ereignissen
in Ordnung ist. Da das Signal als klein angenommen wird ist hier die Anzahl der
Gesamtereignisse auf 7 mit einem möglichen Bereich von 0 bis 50 Ereignissen gesetzt.
H→
Events / ( 2 )
γ γ
80
70
60
50
40
30
20
10
0
100 105 110 115 120 125 130 135 140 145 150
in GeV
Abbildung 4: Daten nach den Schnitten aus [ATLAS11] für den Übergang eines theoretischen
Higgs-Bosons H → γγ
HybridCalculatorOriginal myHybridCalc ( impData_2 , model , bg_ext_pdf , ←↪
&nuisance_parameters , &prior_pdf ) ;
//for 2*lnQ
myHybridCalc . SetTestStatistic ( 1 ) ;
myHybridCalc . SetNumberOfToys ( ntoys ) ;
myHybridCalc . UseNuisance ( true ) ;
myHybridCalc . SetGenerateBinned ( true ) ;
std : : cout