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.2 by buchmann, Thu Jul 14 06:25:39 2011 UTC vs.
Revision 1.16 by fronga, Mon Aug 8 14:21:55 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/Var_Studies.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 49 | 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);
56	>	define_samples(showList,allsamples,signalsamples,scansample);
57		setlumi(luminosity);
58		do_png(true);
59	<	do_pdf(false);
59	>	do_pdf(true);
60		do_eps(false);
61		do_C(false);
62	<	set_directory("comparison__corrected__fixed");//here you can give it a name, for instance "1234pb" or whatever you like :-)
62	>	set_directory("summer_mc");//here you can give it a name, for instance "1234pb" or whatever you like :-)
63		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!
64		stringstream resultsummary;
65
66		//what to do :
67		bool do_all=false;/// DONE
68	<	bool do_peak_finding=false; /// DONE
68	>	bool do_peak_finding=true; /// DONE
69		bool calculate_pred_and_observed=false; /// DONE
70	<	bool study_sidebands=true; /// FOR REAL
70	>	bool study_sidebands=false; /// FOR REAL
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	<
73	>	bool do_find_sideband_definition=false;
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 79 | 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;
100	<
101	<	//jzb response correction
102	<	if(do_jzb_correction) find_correction_factors(jzbvariabledata,jzbvariablemc);
103	<	cout << "Please update Modules/Setup.C to reflect the following values (round them first ... )" <<endl;
104	<	cout << "Corrected JZB variable definition " << endl << "  Data: " << jzbvariabledata << "  and MC: " << jzbvariablemc << endl;
105	<	cout << endl;
106	<	cout << "If you're feeling lazy, copy & paste this : " << endl;
107	<	cout << "  string jzbvariabledata=\""<<jzbvariabledata<<"\";"<<endl;
108	<	cout << "  string jzbvariablemc=\""  <<jzbvariablemc<<"\";"<<endl;
97	>	dout << "With peak correction, we get : " << endl;
98	>	dout << "    Data : " << datajzb.str() << endl;
99	>	dout << "    MC : " << mcjzb.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	>	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	<	if(do_test) test();
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	+	dout << "Observed: " << endl;
130	+	pick_up_events((const char*)(cutmass&&cutOSSF&&cutnJets&&basiccut&&"((jzb[1]+0.06*pt-2.84727)>100)"));
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	+	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	+	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	+	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	+	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	+	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	+	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_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