Skip to content

Commit cab8668

Browse files
author
owen
committed
 Owen : Files needed to run the RA2b likelihood analysis. Initial version is to run the analysis for AN-11-257, v1. Will update to AN-11-257, v2 soon.
1 parent 17f0bcb commit cab8668

33 files changed

+6421
-0
lines changed

‎an-11-257-v2-files/README.txt

+39
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,39 @@
1+
2+
3+
August 6, 2011
4+
5+
Owen Long, UCR
6+
7+
This directory and its subdirectories contains all of the input files and
8+
analysis macros needed to run the RA2b likelihood analysis and produce
9+
susy exclusion plots (though currently only for the old susy scan MC).
10+
11+
In everything below, I'll assume you are on a unix/linux machine, that
12+
you have the path to root set up, and that you are using tcsh (not bash).
13+
If your shell is not tcsh, just type "tcsh" to start one.
14+
15+
First, set up the output directories by doing this command from this
16+
directory
17+
18+
source sourceme-setup-directories
19+
20+
21+
All other things should be run from the directory below this one,
22+
which is UserCode/SusyAnalysis/RA2b/Statistics.
23+
24+
25+
To run the likelihood analysis using the data inputs
26+
27+
source an-11-257-v2-files/analysis-macros/sourceme-data-results
28+
29+
30+
To run the likelihood analysis using the MC inputs
31+
32+
source an-11-257-v2-files/analysis-macros/sourceme-mc-results
33+
34+
35+
To run the susy scans
36+
37+
source an-11-257-v2-files/analysis-macros/sourceme-scans
38+
39+

‎an-11-257-v2-files/analysis-macros/MCSMonly-lm9-loose-sel.c

+110
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,110 @@
1+
2+
3+
void MCSMonly_lm9_loose_sel() {
4+
5+
gROOT->LoadMacro("ra2bRoostatsClass3.c+") ;
6+
7+
float low,high ;
8+
float susylow, susyhigh ;
9+
10+
// ra2bRoostatsClass3( bool ArgUseSigTtwjVar=false, bool ArgUseLdpVars=true ) ;
11+
12+
ra2bRoostatsClass3 rfitNoSusy(1,0,1) ;
13+
rfitNoSusy.initialize("an-11-257-v2-files/input-files/mc-inputs-SMonly-loose-sel.txt") ;
14+
rfitNoSusy.setAndFixSusySig(0.) ;
15+
16+
rfitNoSusy.doFit() ;
17+
printf("\n\n ======== best fit, no susy included ==============\n\n") ;
18+
rfitNoSusy.parameterSnapshot() ;
19+
rfitNoSusy.fitQualityPlot(0,"an-11-257-v2-files/output-files/fitQual-MCSMonly-BestFit-NoSusy-loose-sel.png") ;
20+
rfitNoSusy.fitQualityPlot(1,"an-11-257-v2-files/output-files/fitQualNorm-MCSMonly-BestFit-NoSusy-loose-sel.png") ;
21+
rfitNoSusy.reinitialize() ;
22+
rfitNoSusy.setAndFixSusySig(0.) ;
23+
rfitNoSusy.doFit() ;
24+
rfitNoSusy.profilettwjSig(low,high,1,"an-11-257-v2-files/output-files/ttwjSigProf-MCSMonly-NoSusy-loose-sel.png", 120. ) ;
25+
rfitNoSusy.reinitialize() ;
26+
rfitNoSusy.setAndFixSusySig(0.) ;
27+
rfitNoSusy.doFit() ;
28+
rfitNoSusy.profileqcdSig(low,high,1,"an-11-257-v2-files/output-files/qcdSigProf-MCSMonly-NoSusy-loose-sel.png", 20. ) ;
29+
rfitNoSusy.reinitialize() ;
30+
rfitNoSusy.setAndFixSusySig(0.) ;
31+
rfitNoSusy.doFit() ;
32+
rfitNoSusy.profileZnnSig(low,high,1,"an-11-257-v2-files/output-files/znnSigProf-MCSMonly-NoSusy-loose-sel.png", 40. ) ;
33+
34+
35+
36+
37+
38+
39+
40+
41+
42+
//-------------------------
43+
44+
ra2bRoostatsClass3 rfitSusyFloat(1,0,1) ;
45+
rfitSusyFloat.initialize("an-11-257-v2-files/input-files/mc-inputs-SMonly-loose-sel.txt") ;
46+
rfitSusyFloat.setSusyScanPoint("an-11-257-v2-files/input-files/susy-lm9-nb1-ht350-sb150_200-sig200.txt",0,0) ;
47+
48+
rfitSusyFloat.doFit() ;
49+
printf("\n\n ======== best fit, susy floating ==============\n\n") ;
50+
rfitSusyFloat.parameterSnapshot() ;
51+
rfitSusyFloat.fitQualityPlot(0,"an-11-257-v2-files/output-files/fitQual-MCSMonly-BestFit-susyFloat-lm9-loose-sel.png") ;
52+
rfitSusyFloat.fitQualityPlot(1,"an-11-257-v2-files/output-files/fitQualNorm-MCSMonly-BestFit-susyFloat-lm9-loose-sel.png") ;
53+
rfitSusyFloat.profileSusySig(susylow,susyhigh,1,"an-11-257-v2-files/output-files/susySigProf-MCSMonly-susyFloat-lm9-loose-sel.png", 120. ) ;
54+
rfitSusyFloat.reinitialize() ;
55+
rfitSusyFloat.doFit() ;
56+
rfitSusyFloat.profilettwjSig(low,high,1,"an-11-257-v2-files/output-files/ttwjSigProf-MCSMonly-susyFloat-lm9-loose-sel.png", 120. ) ;
57+
rfitSusyFloat.reinitialize() ;
58+
rfitSusyFloat.doFit() ;
59+
rfitSusyFloat.profileqcdSig(low,high,1,"an-11-257-v2-files/output-files/qcdSigProf-MCSMonly-susyFloat-lm9-loose-sel.png", 20. ) ;
60+
rfitSusyFloat.reinitialize() ;
61+
rfitSusyFloat.doFit() ;
62+
rfitSusyFloat.profileZnnSig(low,high,1,"an-11-257-v2-files/output-files/znnSigProf-MCSMonly-susyFloat-lm9-loose-sel.png", 40. ) ;
63+
64+
65+
rfitSusyFloat.reinitialize() ;
66+
rfitSusyFloat.doFit() ;
67+
rfitSusyFloat.setAndFixSusySig(susylow) ;
68+
rfitSusyFloat.doFit() ;
69+
printf("\n\n ======== lower limit , susy floating==============\n\n") ;
70+
rfitSusyFloat.parameterSnapshot() ;
71+
rfitSusyFloat.fitQualityPlot(0,"an-11-257-v2-files/output-files/fitQual-MCSMonly-susyFloat-LowerLimit-lm9-loose-sel.png") ;
72+
rfitSusyFloat.fitQualityPlot(1,"an-11-257-v2-files/output-files/fitQualNorm-MCSMonly-susyFloat-LowerLimit-lm9-loose-sel.png") ;
73+
74+
rfitSusyFloat.setAndFixSusySig(susyhigh) ;
75+
rfitSusyFloat.doFit() ;
76+
printf("\n\n ======== upper limit , susy floating==============\n\n") ;
77+
rfitSusyFloat.parameterSnapshot() ;
78+
rfitSusyFloat.fitQualityPlot(0,"an-11-257-v2-files/output-files/fitQual-MCSMonly-susyFloat-UpperLimit-lm9-loose-sel.png") ;
79+
rfitSusyFloat.fitQualityPlot(1,"an-11-257-v2-files/output-files/fitQualNorm-MCSMonly-susyFloat-UpperLimit-lm9-loose-sel.png") ;
80+
81+
82+
//-------------------------
83+
84+
ra2bRoostatsClass3 rfitSusyFix(1,0,1) ;
85+
rfitSusyFix.initialize("an-11-257-v2-files/input-files/mc-inputs-SMonly-loose-sel.txt") ;
86+
rfitSusyFix.setSusyScanPoint("an-11-257-v2-files/input-files/susy-lm9-nb1-ht350-sb150_200-sig200.txt",0,0) ;
87+
rfitSusyFix.setAndFixSusySigToPredictedValue() ;
88+
89+
rfitSusyFix.doFit() ;
90+
printf("\n\n ======== best fit, susy fixed ==============\n\n") ;
91+
rfitSusyFix.parameterSnapshot() ;
92+
rfitSusyFix.fitQualityPlot(0,"an-11-257-v2-files/output-files/fitQual-MCSMonly-BestFit-susyFixed-lm9-loose-sel.png") ;
93+
rfitSusyFix.fitQualityPlot(1,"an-11-257-v2-files/output-files/fitQualNorm-MCSMonly-BestFit-susyFixed-lm9-loose-sel.png") ;
94+
rfitSusyFix.reinitialize() ;
95+
rfitSusyFix.setAndFixSusySigToPredictedValue() ;
96+
rfitSusyFix.doFit() ;
97+
rfitSusyFix.profilettwjSig(low,high,1,"an-11-257-v2-files/output-files/ttwjSigProf-MCSMonly-susyFixed-lm9-loose-sel.png", 120. ) ;
98+
rfitSusyFix.reinitialize() ;
99+
rfitSusyFix.setAndFixSusySigToPredictedValue() ;
100+
rfitSusyFix.doFit() ;
101+
rfitSusyFix.profileqcdSig(low,high,1,"an-11-257-v2-files/output-files/qcdSigProf-MCSMonly-susyFixed-lm9-loose-sel.png", 20. ) ;
102+
rfitSusyFix.reinitialize() ;
103+
rfitSusyFix.setAndFixSusySigToPredictedValue() ;
104+
rfitSusyFix.doFit() ;
105+
rfitSusyFix.profileZnnSig(low,high,1,"an-11-257-v2-files/output-files/znnSigProf-MCSMonly-susyFixed-lm9-loose-sel.png", 40. ) ;
106+
107+
108+
109+
110+
}

‎an-11-257-v2-files/analysis-macros/MCSMonly-lm9-tight-sel.c

+110
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,110 @@
1+
2+
3+
void MCSMonly_lm9_tight_sel() {
4+
5+
gROOT->LoadMacro("ra2bRoostatsClass3.c+") ;
6+
7+
float low,high ;
8+
float susylow, susyhigh ;
9+
10+
// ra2bRoostatsClass3( bool ArgUseSigTtwjVar=false, bool ArgUseLdpVars=true ) ;
11+
12+
ra2bRoostatsClass3 rfitNoSusy(1,0,2) ;
13+
rfitNoSusy.initialize("an-11-257-v2-files/input-files/mc-inputs-SMonly-tight-sel.txt") ;
14+
rfitNoSusy.setAndFixSusySig(0.) ;
15+
16+
rfitNoSusy.doFit() ;
17+
printf("\n\n ======== best fit, no susy included ==============\n\n") ;
18+
rfitNoSusy.parameterSnapshot() ;
19+
rfitNoSusy.fitQualityPlot(0,"an-11-257-v2-files/output-files/fitQual-MCSMonly-BestFit-NoSusy-tight-sel.png") ;
20+
rfitNoSusy.fitQualityPlot(1,"an-11-257-v2-files/output-files/fitQualNorm-MCSMonly-BestFit-NoSusy-tight-sel.png") ;
21+
rfitNoSusy.reinitialize() ;
22+
rfitNoSusy.setAndFixSusySig(0.) ;
23+
rfitNoSusy.doFit() ;
24+
rfitNoSusy.profilettwjSig(low,high,1,"an-11-257-v2-files/output-files/ttwjSigProf-MCSMonly-NoSusy-tight-sel.png", 25. ) ;
25+
rfitNoSusy.reinitialize() ;
26+
rfitNoSusy.setAndFixSusySig(0.) ;
27+
rfitNoSusy.doFit() ;
28+
rfitNoSusy.profileqcdSig(low,high,1,"an-11-257-v2-files/output-files/qcdSigProf-MCSMonly-NoSusy-tight-sel.png", 5. ) ;
29+
rfitNoSusy.reinitialize() ;
30+
rfitNoSusy.setAndFixSusySig(0.) ;
31+
rfitNoSusy.doFit() ;
32+
rfitNoSusy.profileZnnSig(low,high,1,"an-11-257-v2-files/output-files/znnSigProf-MCSMonly-NoSusy-tight-sel.png", 10. ) ;
33+
34+
35+
36+
37+
38+
39+
40+
41+
42+
//-------------------------
43+
44+
ra2bRoostatsClass3 rfitSusyFloat(1,0,2) ;
45+
rfitSusyFloat.initialize("an-11-257-v2-files/input-files/mc-inputs-SMonly-tight-sel.txt") ;
46+
rfitSusyFloat.setSusyScanPoint("an-11-257-v2-files/input-files/susy-lm9-nb1-ht500-sb150_200-sig300.txt",0,0) ;
47+
48+
rfitSusyFloat.doFit() ;
49+
printf("\n\n ======== best fit, susy floating ==============\n\n") ;
50+
rfitSusyFloat.parameterSnapshot() ;
51+
rfitSusyFloat.fitQualityPlot(0,"an-11-257-v2-files/output-files/fitQual-MCSMonly-BestFit-susyFloat-lm9-tight-sel.png") ;
52+
rfitSusyFloat.fitQualityPlot(1,"an-11-257-v2-files/output-files/fitQualNorm-MCSMonly-BestFit-susyFloat-lm9-tight-sel.png") ;
53+
rfitSusyFloat.profileSusySig(susylow,susyhigh,1,"an-11-257-v2-files/output-files/susySigProf-MCSMonly-susyFloat-lm9-tight-sel.png", 25. ) ;
54+
rfitSusyFloat.reinitialize() ;
55+
rfitSusyFloat.doFit() ;
56+
rfitSusyFloat.profilettwjSig(low,high,1,"an-11-257-v2-files/output-files/ttwjSigProf-MCSMonly-susyFloat-lm9-tight-sel.png", 25. ) ;
57+
rfitSusyFloat.reinitialize() ;
58+
rfitSusyFloat.doFit() ;
59+
rfitSusyFloat.profileqcdSig(low,high,1,"an-11-257-v2-files/output-files/qcdSigProf-MCSMonly-susyFloat-lm9-tight-sel.png", 5. ) ;
60+
rfitSusyFloat.reinitialize() ;
61+
rfitSusyFloat.doFit() ;
62+
rfitSusyFloat.profileZnnSig(low,high,1,"an-11-257-v2-files/output-files/znnSigProf-MCSMonly-susyFloat-lm9-tight-sel.png", 10. ) ;
63+
64+
65+
rfitSusyFloat.reinitialize() ;
66+
rfitSusyFloat.doFit() ;
67+
rfitSusyFloat.setAndFixSusySig(susylow) ;
68+
rfitSusyFloat.doFit() ;
69+
printf("\n\n ======== lower limit , susy floating==============\n\n") ;
70+
rfitSusyFloat.parameterSnapshot() ;
71+
rfitSusyFloat.fitQualityPlot(0,"an-11-257-v2-files/output-files/fitQual-MCSMonly-susyFloat-LowerLimit-lm9-tight-sel.png") ;
72+
rfitSusyFloat.fitQualityPlot(1,"an-11-257-v2-files/output-files/fitQualNorm-MCSMonly-susyFloat-LowerLimit-lm9-tight-sel.png") ;
73+
74+
rfitSusyFloat.setAndFixSusySig(susyhigh) ;
75+
rfitSusyFloat.doFit() ;
76+
printf("\n\n ======== upper limit , susy floating==============\n\n") ;
77+
rfitSusyFloat.parameterSnapshot() ;
78+
rfitSusyFloat.fitQualityPlot(0,"an-11-257-v2-files/output-files/fitQual-MCSMonly-susyFloat-UpperLimit-lm9-tight-sel.png") ;
79+
rfitSusyFloat.fitQualityPlot(1,"an-11-257-v2-files/output-files/fitQualNorm-MCSMonly-susyFloat-UpperLimit-lm9-tight-sel.png") ;
80+
81+
82+
//-------------------------
83+
84+
ra2bRoostatsClass3 rfitSusyFix(1,0,2) ;
85+
rfitSusyFix.initialize("an-11-257-v2-files/input-files/mc-inputs-SMonly-tight-sel.txt") ;
86+
rfitSusyFix.setSusyScanPoint("an-11-257-v2-files/input-files/susy-lm9-nb1-ht500-sb150_200-sig300.txt",0,0) ;
87+
rfitSusyFix.setAndFixSusySigToPredictedValue() ;
88+
89+
rfitSusyFix.doFit() ;
90+
printf("\n\n ======== best fit, susy fixed ==============\n\n") ;
91+
rfitSusyFix.parameterSnapshot() ;
92+
rfitSusyFix.fitQualityPlot(0,"an-11-257-v2-files/output-files/fitQual-MCSMonly-BestFit-susyFixed-lm9-tight-sel.png") ;
93+
rfitSusyFix.fitQualityPlot(1,"an-11-257-v2-files/output-files/fitQualNorm-MCSMonly-BestFit-susyFixed-lm9-tight-sel.png") ;
94+
rfitSusyFix.reinitialize() ;
95+
rfitSusyFix.setAndFixSusySigToPredictedValue() ;
96+
rfitSusyFix.doFit() ;
97+
rfitSusyFix.profilettwjSig(low,high,1,"an-11-257-v2-files/output-files/ttwjSigProf-MCSMonly-susyFixed-lm9-tight-sel.png", 25. ) ;
98+
rfitSusyFix.reinitialize() ;
99+
rfitSusyFix.setAndFixSusySigToPredictedValue() ;
100+
rfitSusyFix.doFit() ;
101+
rfitSusyFix.profileqcdSig(low,high,1,"an-11-257-v2-files/output-files/qcdSigProf-MCSMonly-susyFixed-lm9-tight-sel.png", 5. ) ;
102+
rfitSusyFix.reinitialize() ;
103+
rfitSusyFix.setAndFixSusySigToPredictedValue() ;
104+
rfitSusyFix.doFit() ;
105+
rfitSusyFix.profileZnnSig(low,high,1,"an-11-257-v2-files/output-files/znnSigProf-MCSMonly-susyFixed-lm9-tight-sel.png", 10. ) ;
106+
107+
108+
109+
110+
}

‎an-11-257-v2-files/analysis-macros/MCSMsusy-lm9-loose-sel.c

+110
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,110 @@
1+
2+
3+
void MCSMsusy_lm9_loose_sel() {
4+
5+
gROOT->LoadMacro("ra2bRoostatsClass3.c+") ;
6+
7+
float low,high ;
8+
float susylow, susyhigh ;
9+
10+
// ra2bRoostatsClass3( bool ArgUseSigTtwjVar=false, bool ArgUseLdpVars=true ) ;
11+
12+
ra2bRoostatsClass3 rfitNoSusy(1,0,1) ;
13+
rfitNoSusy.initialize("an-11-257-v2-files/input-files/mc-inputs-SM+lm9-loose-sel.txt") ;
14+
rfitNoSusy.setAndFixSusySig(0.) ;
15+
16+
rfitNoSusy.doFit() ;
17+
printf("\n\n ======== best fit, no susy included ==============\n\n") ;
18+
rfitNoSusy.parameterSnapshot() ;
19+
rfitNoSusy.fitQualityPlot(0,"an-11-257-v2-files/output-files/fitQual-MCSMsusy-BestFit-NoSusy-loose-sel.png") ;
20+
rfitNoSusy.fitQualityPlot(1,"an-11-257-v2-files/output-files/fitQualNorm-MCSMsusy-BestFit-NoSusy-loose-sel.png") ;
21+
rfitNoSusy.reinitialize() ;
22+
rfitNoSusy.setAndFixSusySig(0.) ;
23+
rfitNoSusy.doFit() ;
24+
rfitNoSusy.profilettwjSig(low,high,1,"an-11-257-v2-files/output-files/ttwjSigProf-MCSMsusy-NoSusy-loose-sel.png", 220. ) ;
25+
rfitNoSusy.reinitialize() ;
26+
rfitNoSusy.setAndFixSusySig(0.) ;
27+
rfitNoSusy.doFit() ;
28+
rfitNoSusy.profileqcdSig(low,high,1,"an-11-257-v2-files/output-files/qcdSigProf-MCSMsusy-NoSusy-loose-sel.png", 20. ) ;
29+
rfitNoSusy.reinitialize() ;
30+
rfitNoSusy.setAndFixSusySig(0.) ;
31+
rfitNoSusy.doFit() ;
32+
rfitNoSusy.profileZnnSig(low,high,1,"an-11-257-v2-files/output-files/znnSigProf-MCSMsusy-NoSusy-loose-sel.png", 40. ) ;
33+
34+
35+
36+
37+
38+
39+
40+
41+
42+
//-------------------------
43+
44+
ra2bRoostatsClass3 rfitSusyFloat(1,0,1) ;
45+
rfitSusyFloat.initialize("an-11-257-v2-files/input-files/mc-inputs-SM+lm9-loose-sel.txt") ;
46+
rfitSusyFloat.setSusyScanPoint("an-11-257-v2-files/input-files/susy-lm9-nb1-ht350-sb150_200-sig200.txt",0,0) ;
47+
48+
rfitSusyFloat.doFit() ;
49+
printf("\n\n ======== best fit, susy floating ==============\n\n") ;
50+
rfitSusyFloat.parameterSnapshot() ;
51+
rfitSusyFloat.fitQualityPlot(0,"an-11-257-v2-files/output-files/fitQual-MCSMsusy-BestFit-susyFloat-lm9-loose-sel.png") ;
52+
rfitSusyFloat.fitQualityPlot(1,"an-11-257-v2-files/output-files/fitQualNorm-MCSMsusy-BestFit-susyFloat-lm9-loose-sel.png") ;
53+
rfitSusyFloat.profileSusySig(susylow,susyhigh,1,"an-11-257-v2-files/output-files/susySigProf-MCSMsusy-susyFloat-lm9-loose-sel.png", 220. ) ;
54+
rfitSusyFloat.reinitialize() ;
55+
rfitSusyFloat.doFit() ;
56+
rfitSusyFloat.profilettwjSig(low,high,1,"an-11-257-v2-files/output-files/ttwjSigProf-MCSMsusy-susyFloat-lm9-loose-sel.png", 220. ) ;
57+
rfitSusyFloat.reinitialize() ;
58+
rfitSusyFloat.doFit() ;
59+
rfitSusyFloat.profileqcdSig(low,high,1,"an-11-257-v2-files/output-files/qcdSigProf-MCSMsusy-susyFloat-lm9-loose-sel.png", 20. ) ;
60+
rfitSusyFloat.reinitialize() ;
61+
rfitSusyFloat.doFit() ;
62+
rfitSusyFloat.profileZnnSig(low,high,1,"an-11-257-v2-files/output-files/znnSigProf-MCSMsusy-susyFloat-lm9-loose-sel.png", 40. ) ;
63+
64+
65+
rfitSusyFloat.reinitialize() ;
66+
rfitSusyFloat.doFit() ;
67+
rfitSusyFloat.setAndFixSusySig(susylow) ;
68+
rfitSusyFloat.doFit() ;
69+
printf("\n\n ======== lower limit , susy floating==============\n\n") ;
70+
rfitSusyFloat.parameterSnapshot() ;
71+
rfitSusyFloat.fitQualityPlot(0,"an-11-257-v2-files/output-files/fitQual-MCSMsusy-susyFloat-LowerLimit-lm9-loose-sel.png") ;
72+
rfitSusyFloat.fitQualityPlot(1,"an-11-257-v2-files/output-files/fitQualNorm-MCSMsusy-susyFloat-LowerLimit-lm9-loose-sel.png") ;
73+
74+
rfitSusyFloat.setAndFixSusySig(susyhigh) ;
75+
rfitSusyFloat.doFit() ;
76+
printf("\n\n ======== upper limit , susy floating==============\n\n") ;
77+
rfitSusyFloat.parameterSnapshot() ;
78+
rfitSusyFloat.fitQualityPlot(0,"an-11-257-v2-files/output-files/fitQual-MCSMsusy-susyFloat-UpperLimit-lm9-loose-sel.png") ;
79+
rfitSusyFloat.fitQualityPlot(1,"an-11-257-v2-files/output-files/fitQualNorm-MCSMsusy-susyFloat-UpperLimit-lm9-loose-sel.png") ;
80+
81+
82+
//-------------------------
83+
84+
ra2bRoostatsClass3 rfitSusyFix(1,0,1) ;
85+
rfitSusyFix.initialize("an-11-257-v2-files/input-files/mc-inputs-SM+lm9-loose-sel.txt") ;
86+
rfitSusyFix.setSusyScanPoint("an-11-257-v2-files/input-files/susy-lm9-nb1-ht350-sb150_200-sig200.txt",0,0) ;
87+
rfitSusyFix.setAndFixSusySigToPredictedValue() ;
88+
89+
rfitSusyFix.doFit() ;
90+
printf("\n\n ======== best fit, susy fixed ==============\n\n") ;
91+
rfitSusyFix.parameterSnapshot() ;
92+
rfitSusyFix.fitQualityPlot(0,"an-11-257-v2-files/output-files/fitQual-MCSMsusy-BestFit-susyFixed-lm9-loose-sel.png") ;
93+
rfitSusyFix.fitQualityPlot(1,"an-11-257-v2-files/output-files/fitQualNorm-MCSMsusy-BestFit-susyFixed-lm9-loose-sel.png") ;
94+
rfitSusyFix.reinitialize() ;
95+
rfitSusyFix.setAndFixSusySigToPredictedValue() ;
96+
rfitSusyFix.doFit() ;
97+
rfitSusyFix.profilettwjSig(low,high,1,"an-11-257-v2-files/output-files/ttwjSigProf-MCSMsusy-susyFixed-lm9-loose-sel.png", 220. ) ;
98+
rfitSusyFix.reinitialize() ;
99+
rfitSusyFix.setAndFixSusySigToPredictedValue() ;
100+
rfitSusyFix.doFit() ;
101+
rfitSusyFix.profileqcdSig(low,high,1,"an-11-257-v2-files/output-files/qcdSigProf-MCSMsusy-susyFixed-lm9-loose-sel.png", 20. ) ;
102+
rfitSusyFix.reinitialize() ;
103+
rfitSusyFix.setAndFixSusySigToPredictedValue() ;
104+
rfitSusyFix.doFit() ;
105+
rfitSusyFix.profileZnnSig(low,high,1,"an-11-257-v2-files/output-files/znnSigProf-MCSMsusy-susyFixed-lm9-loose-sel.png", 40. ) ;
106+
107+
108+
109+
110+
}

0 commit comments

Comments
 (0)