ViewVC Help

View File | Revision Log | Show Annotations | Root Listing

root/cvsroot/UserCode/cbrown/AnalysisFramework/Plotting/various_studies.C

Comparing UserCode/cbrown/AnalysisFramework/Plotting/various_studies.C (file contents):
Revision 1.5 by buchmann, Tue Jul 19 07:06:51 2011 UTC vs.
Revision 1.15 by buchmann, Thu Jul 28 19:26:12 2011 UTC

#	Line 9 | Line 9 | to compile: just run make ... that shoul
9		#include "Modules/GeneralToolBox.C"
10		#include "Modules/SampleClass.C"
11		#include "Modules/PeakFinder.C"
12	–	#include "Modules/HistoDB.C"
12		#include "Modules/Poisson_Calculator.C"
13		#include "Modules/setTDRStyle.C"
14		#include "Modules/ActiveSamples.C"
15		#include "Modules/UpperLimitsWithShape.C"
16		#include "Modules/Setup.C"
17	+	#include "Modules/external/cl95cms.c"
18		#include "Modules/Plotting_Functions.C" //also included for peak finding etc.
19		#include "Modules/StudyModule.C" //also included for peak finding etc.
20	+	#include "Modules/ExperimentalModule.C"
21	+	#include "Modules/ResultModule.C"
22	+	#include "Modules/LimitCalculation.C"
23	+	#include "Modules/Systematics.C"
24	+	#include "Modules/SUSYScan.C"
25
26		#include <TCut.h>
27		#include <TROOT.h>
#	Line 48 | Line 53 | int main()
53		gStyle->SetTextFont(42);
54		bool showList=false;
55		set_treename("events");//you can set the treename here to be used; options are "events" (for reco) for "PFevents" (for particle flow)
56	<	define_samples(showList,allsamples,signalsamples);
56	>	define_samples(showList,allsamples,signalsamples,scansample);
57		setlumi(luminosity);
58		do_png(true);
59		do_pdf(false);
#	Line 66 | Line 71 | int main()
71		bool do_test=false;/// DONE    (just any test you wish)
72		bool do_jzb_correction=false; //use this to find out the correction factor; originally this was done automatically but now you should do this separately and update Modules/Setup.C
73		bool do_find_sideband_definition=false;
74	<	bool do_pick_up_events=true; ///DONE
75	<
74	>	bool do_pick_up_events=false; ///DONE
75	>	bool do_ttbar_with_shapes=false;
76	>	bool do_upper_limits=false;
77	>	bool do_run_check=false;
78	>	bool do_plot_list=true;
79	>
80	>	bool do_show_dibosons=false;
81	>
82	>	bool do_efficiency_scan_in_susy_space=false; /// NOW OUTSOURCED TO T3
83	>	bool requireZ=true;
84	>
85		//**** part 1 : peak finding
86		float MCPeak=0,MCPeakError=0,DataPeak=0,DataPeakError=0,MCSigma=10,DataSigma=10;
87		method=Kostasmethod;//Kostasmethod;//dogaus3sigma;// options: dogaus,doKM,dogaus2sigma,dogaus3sigma
88	<	if(do_peak_finding||calculate_pred_and_observed||do_all) find_peaks(MCPeak,MCPeakError, DataPeak, DataPeakError,MCSigma,DataSigma,resultsummary);
88	>	if(do_peak_finding||do_show_dibosons||calculate_pred_and_observed||do_all) find_peaks(MCPeak,MCPeakError, DataPeak, DataPeakError,MCSigma,DataSigma,resultsummary);
89
90		stringstream datajzb;
91		if(DataPeak>0) datajzb<<"("<<jzbvariabledata<<"-"<<TMath::Abs(DataPeak)<<")";
#	Line 80 | Line 94 | int main()
94		if(MCPeak>0) mcjzb<<"("<<jzbvariablemc<<"-"<<TMath::Abs(MCPeak)<<")";
95		else mcjzb<<"("<<jzbvariablemc<<"+"<<TMath::Abs(MCPeak)<<")";
96
97	<	cout << "With peak correction, we get : " << endl;
98	<	cout << "    Data : " << datajzb.str() << endl;
99	<	cout << "    MC : " << datajzb.str() << endl;
97	>	dout << "With peak correction, we get : " << endl;
98	>	dout << "    Data : " << datajzb.str() << endl;
99	>	dout << "    MC : " << datajzb.str() << endl;
100
101		if(do_find_sideband_definition) find_sideband_definition();
102
103		if(do_jzb_correction) {
104		find_correction_factors(jzbvariabledata,jzbvariablemc);
105	<	cout << "Please update Modules/Setup.C to reflect the following values (round them first ... )" <<endl;
106	<	cout << "Corrected JZB variable definition " << endl << "  Data: " << jzbvariabledata << "  and MC: " << jzbvariablemc << endl;
107	<	cout << endl;
108	<	cout << "If you're feeling lazy, copy & paste this : " << endl;
109	<	cout << "  string jzbvariabledata=\""<<jzbvariabledata<<"\";"<<endl;
110	<	cout << "  string jzbvariablemc=\""  <<jzbvariablemc<<"\";"<<endl;
105	>	dout << "Please update Modules/Setup.C to reflect the following values (round them first ... )" <<endl;
106	>	dout << "Corrected JZB variable definition " << endl << "  Data: " << jzbvariabledata << "  and MC: " << jzbvariablemc << endl;
107	>	dout << endl;
108	>	dout << "If you're feeling lazy, copy & paste this : " << endl;
109	>	dout << "  string jzbvariabledata=\""<<jzbvariabledata<<"\";"<<endl;
110	>	dout << "  string jzbvariablemc=\""  <<jzbvariablemc<<"\";"<<endl;
111		}
112
113		if(study_sidebands) look_at_sidebands(mcjzb.str(),datajzb.str());
114
115	+	vector<float>jzb_limit_bins;
116	+	jzb_limit_bins.push_back(50);
117	+	//jzb_limit_bins.push_back(75);jzb_limit_bins.push_back(100);
118	+	//jzb_limit_bins.push_back(150);jzb_limit_bins.push_back(200);jzb_limit_bins.push_back(500);
119	+
120	+
121	+	if(do_ttbar_with_shapes) prepare_ttbar_limits(mcjzb.str(),datajzb.str(),DataPeakError,MCPeakError,jzb_limit_bins);
122	+
123	+	if(do_upper_limits) calculate_upper_limits(mcjzb.str(),datajzb.str());
124	+
125	+	if(do_efficiency_scan_in_susy_space) efficiency_scan_in_susy_space(mcjzb.str(),datajzb.str(),requireZ,MCPeakError);
126	+
127		if(do_pick_up_events) {
128
129	<	cout << "Observed: " << endl;
129	>	dout << "Observed: " << endl;
130		pick_up_events((const char*)(cutmass&&cutOSSF&&cutnJets&&basiccut&&"((jzb[1]+0.06*pt-2.84727)>100)"));
131	<	cout << "Predicted (JZB<-100) OSSF" << endl;
131	>	dout << "Predicted (JZB<-100) OSSF" << endl;
132		pick_up_events((const char*)(cutmass&&cutOSSF&&cutnJets&&basiccut&&"((jzb[1]+0.06*pt-2.84727)<-100)"));
133	<	cout << "Predicted (emu, JZB>100) OSOF" << endl;
133	>	dout << "Predicted (emu, JZB>100) OSOF" << endl;
134		pick_up_events((const char*)(cutmass&&cutOSOF&&cutnJets&&basiccut&&"((jzb[1]+0.06*pt-2.84727)>100)"));
135	<	cout << "Predicted (emu, JZB<-100) OSOF" << endl;
135	>	dout << "Predicted (emu, JZB<-100) OSOF" << endl;
136		pick_up_events((const char*)(cutmass&&cutOSOF&&cutnJets&&basiccut&&"((jzb[1]+0.06*pt-2.84727)<-100)"));
137	<	cout << "Predicted (SB emu, JZB>100) OSOF" << endl;
137	>	dout << "Predicted (SB emu, JZB>100) OSOF" << endl;
138		pick_up_events((const char*)(sidebandcut&&cutOSOF&&cutnJets&&basiccut&&"((jzb[1]+0.06*pt-2.84727)>100)"));
139	<	cout << "Predicted (SB emu, JZB<-100) OSOF" << endl;
139	>	dout << "Predicted (SB emu, JZB<-100) OSOF" << endl;
140		pick_up_events((const char*)(cutOSOF&&cutnJets&&basiccut&&sidebandcut&&"((jzb[1]+0.06*pt-2.84727)<-100)"));
141	<	cout << "Predicted (SB SF, JZB>100) OSSF" << endl;
141	>	dout << "Predicted (SB SF, JZB>100) OSSF" << endl;
142		pick_up_events((const char*)(sidebandcut&&cutOSSF&&cutnJets&&basiccut&&"((jzb[1]+0.06*pt-2.84727)>100)"));
143	<	cout << "Predicted (SB SF, JZB<-100) OSSF" << endl;
143	>	dout << "Predicted (SB SF, JZB<-100) OSSF" << endl;
144		pick_up_events((const char*)(sidebandcut&&cutOSSF&&cutnJets&&basiccut&&"((jzb[1]+0.06*pt-2.84727)<-100)"));
145		}
146
147	<	if(do_test) test();
147	>	if(do_run_check) run_check();
148	>	if(do_show_dibosons) show_dibosons(datajzb.str(),mcjzb.str());
149
150	+	if(do_plot_list) plot_list();
151
152	+	if(do_test) test();
153
154		return 0;
155		}

Diff Legend

–	Removed lines
+	Added lines
<	Changed lines
>	Changed lines