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.6 by buchmann, Wed Jul 20 08:51:33 2011 UTC vs.
Revision 1.26 by buchmann, Tue Sep 13 09:26:10 2011 UTC

#	Line 8 | Line 8 | to compile: just run make ... that shoul
8		#include <sys/stat.h>
9		#include "Modules/GeneralToolBox.C"
10		#include "Modules/SampleClass.C"
11	+	#include "Modules/Setup.C"
12		#include "Modules/PeakFinder.C"
12	–	#include "Modules/HistoDB.C"
13		#include "Modules/Poisson_Calculator.C"
14		#include "Modules/setTDRStyle.C"
15		#include "Modules/ActiveSamples.C"
16		#include "Modules/UpperLimitsWithShape.C"
17	<	#include "Modules/Setup.C"
17	>	//#include "Modules/external/roostats_cl95.C"
18	>	//#include "Modules/external/cl95cms.c"
19		#include "Modules/Plotting_Functions.C" //also included for peak finding etc.
20	+	#include "Modules/ResultModule.C"
21	+	#include "Modules/ExperimentalModule.C"
22		#include "Modules/StudyModule.C" //also included for peak finding etc.
23	+	#include "Modules/LimitCalculation.C"
24	+	#include "Modules/Systematics.C"
25	+	#include "Modules/CrossSectionReader.C"
26	+	#include "Modules/SUSYScan.C"
27
28		#include <TCut.h>
29		#include <TROOT.h>
#	Line 44 | Line 51 | using namespace PlottingSetup;
51		int main()
52		{
53		gROOT->SetStyle("Plain");
54	<	setTDRStyle();
54	>	bool do_fat_line=false; // if you want to have HistLineWidth=1 and FuncWidth=1 as it was before instead of 2
55	>	setTDRStyle(do_fat_line);
56		gStyle->SetTextFont(42);
57		bool showList=false;
58		set_treename("events");//you can set the treename here to be used; options are "events" (for reco) for "PFevents" (for particle flow)
59	<	define_samples(showList,allsamples,signalsamples);
59	>	define_samples(showList,allsamples,signalsamples,scansample);
60		setlumi(luminosity);
61		do_png(true);
62		do_pdf(false);
63		do_eps(false);
64		do_C(false);
65	<	set_directory("summer_mc");//here you can give it a name, for instance "1234pb" or whatever you like :-)
65	>	set_directory(directoryname);//Indicate the directory name where you'd like to save the output files in Setup.C
66		setessentialcut(essential);//this sets the essential cut; this one is used in the draw command so it is AUTOMATICALLY applied everywhere. IMPORTANT: Do NOT store weights here!
67		stringstream resultsummary;
68
69		//what to do :
70		bool do_all=false;/// DONE
71	<	bool do_peak_finding=false; /// DONE
71	>	bool do_peak_finding=true; /// DONE
72		bool calculate_pred_and_observed=false; /// DONE
73		bool study_sidebands=false; /// FOR REAL
74		bool do_test=false;/// DONE    (just any test you wish)
75		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
76		bool do_find_sideband_definition=false;
77	<	bool do_pick_up_events=true; ///DONE
78	<
77	>	bool do_pick_up_events=false; ///DONE
78	>	bool do_ttbar_with_shapes=false;
79	>	bool do_upper_limits=false;
80	>	bool do_run_check=false;
81	>	bool do_kinematic_dist_of_pred_and_obs=false;//former do_plot_list
82	>
83	>	bool do_show_dibosons=false;
84	>
85	>	bool overlay_signal=false;
86	>
87	>	bool do_efficiency_scan_in_susy_space=false; /// NOW OUTSOURCED TO T3
88	>	bool requireZ=true;
89	>
90	>	bool do_beautiful_ratio_plots=false;
91	>
92		//**** part 1 : peak finding
93		float MCPeak=0,MCPeakError=0,DataPeak=0,DataPeakError=0,MCSigma=10,DataSigma=10;
94		method=Kostasmethod;//Kostasmethod;//dogaus3sigma;// options: dogaus,doKM,dogaus2sigma,dogaus3sigma
95	<	if(do_peak_finding||calculate_pred_and_observed||do_all) find_peaks(MCPeak,MCPeakError, DataPeak, DataPeakError,MCSigma,DataSigma,resultsummary);
95	>	if(do_peak_finding||do_show_dibosons||calculate_pred_and_observed||do_all) find_peaks(MCPeak,MCPeakError, DataPeak, DataPeakError,MCSigma,DataSigma,resultsummary);
96
97		stringstream datajzb;
98		if(DataPeak>0) datajzb<<"("<<jzbvariabledata<<"-"<<TMath::Abs(DataPeak)<<")";
#	Line 82 | Line 103 | int main()
103
104		dout << "With peak correction, we get : " << endl;
105		dout << "    Data : " << datajzb.str() << endl;
106	<	dout << "    MC : " << datajzb.str() << endl;
106	>	dout << "    MC : " << mcjzb.str() << endl;
107
108		if(do_find_sideband_definition) find_sideband_definition();
109
#	Line 98 | Line 119 | int main()
119
120		if(study_sidebands) look_at_sidebands(mcjzb.str(),datajzb.str());
121
122	+	vector<float>jzb_limit_bins;
123	+	jzb_limit_bins.push_back(50);
124	+	//jzb_limit_bins.push_back(75);jzb_limit_bins.push_back(100);
125	+	//jzb_limit_bins.push_back(150);jzb_limit_bins.push_back(200);jzb_limit_bins.push_back(500);
126	+
127	+
128	+	if(do_ttbar_with_shapes) prepare_ttbar_limits(mcjzb.str(),datajzb.str(),DataPeakError,MCPeakError,jzb_limit_bins);
129	+
130	+	if(do_upper_limits) calculate_upper_limits(mcjzb.str(),datajzb.str());
131	+
132		if(do_pick_up_events) {
133
134		dout << "Observed: " << endl;
#	Line 118 | Line 149 | int main()
149		pick_up_events((const char*)(sidebandcut&&cutOSSF&&cutnJets&&basiccut&&"((jzb[1]+0.06*pt-2.84727)<-100)"));
150		}
151
152	+	if(do_run_check) run_check();
153	+	if(do_show_dibosons) show_dibosons(datajzb.str(),mcjzb.str());
154	+
155	+	if(do_kinematic_dist_of_pred_and_obs) kinematic_dist_of_pred_and_obs();
156	+
157		if(do_test) test();
158
159	+	vector<float> ratio_binning;
160	+	ratio_binning.push_back(0);
161	+	ratio_binning.push_back(5);
162	+	ratio_binning.push_back(10);
163	+	ratio_binning.push_back(20);
164	+	ratio_binning.push_back(50);
165	+	ratio_binning.push_back(100);
166	+	ratio_binning.push_back(200);
167	+	ratio_binning.push_back(350);
168	+	//ratio_binning.push_back(500);
169	+	if(do_beautiful_ratio_plots) Poisson_ratio_plots(mcjzb.str(),datajzb.str(),ratio_binning,MCPeakError,DataPeakError);
170
171	+	write_warning(__FUNCTION__,"Need to remove the next line:");
172	+	//  do_new_prediction_plots(mcjzb.str(),datajzb.str(),DataSigma,MCSigma,overlay_signal);
173
174	+	//  ttbar_sidebands_comparison(mcjzb.str());
175	+
176	+	lm0_illustration();
177		return 0;
178		}
179

Diff Legend

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