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

Diff Legend

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