[ViewVC] Diff of: cvsroot/UserCode/cbrown/Development/Plotting/Modules/Plotting

Comparing UserCode/cbrown/Development/Plotting/Modules/Plotting_Functions.C (file contents):
Revision 1.97 by buchmann, Thu May 30 07:59:33 2013 UTC vs.
Revision 1.98 by buchmann, Fri May 31 09:45:15 2013 UTC

#	Line 109 \| Line 109 \| void find_one_peak_combination(TCut spec
109		delete tempcan;
110		}
111
112	–	void Old_FindPeakPileUpCorrection() {
113	–	bool StoreResultsInSeparateFile=false;
114	–	PeakLibrary::StoreHistos=true;
115	–
116	–
117	–	dout << "Finding peak position as a function of PU" << endl;
118	–	float MCPeak, MCPeakError, DataPeak, DataPeakError;
119	–	stringstream result, datajzb, mcjzb;
120	–	bool doPUreweighting=true;
121	–	bool SwitchOffNJetsCut=false;
122	–
123	–
124	–	TCanvas *can = new TCanvas("can","can");
125	–	float mcSigma,mcSigmaError,dataSigma,dataSigmaError;
126	–
127	–	float NumVtxBin[6] = {0,5,10,15,20,30};
128	–
129	–	stringstream NowCut;
130	–
131	–	TGraphErrors *gMCPeak = new TGraphErrors();
132	–	gMCPeak->SetTitle("gMCPeak");
133	–	gMCPeak->SetName("gMCPeak");
134	–	TGraphErrors *gDataPeak = new TGraphErrors();
135	–	gDataPeak->SetTitle("gDataPeak");
136	–	gDataPeak->SetName("gDataPeak");
137	–	TGraphErrors *gMCWidth = new TGraphErrors();
138	–	gMCWidth->SetTitle("gMCWidth");
139	–	gMCWidth->SetName("gMCWidth");
140	–	TGraphErrors *gDataWidth = new TGraphErrors();
141	–	gDataWidth->SetTitle("gDataWidth");
142	–	gDataWidth->SetName("gDataWidth");
143	–
144	–
145	–	for(int i=0;i<5;i++) {
146	–	NowCut.str("");
147	–	NowCut << "numVtx>=" << NumVtxBin[i] << "&&numVtx<" << NumVtxBin[i+1];
148	–	find_one_peak_combination(TCut(NowCut.str().c_str()),SwitchOffNJetsCut,MCPeak,MCPeakError, DataPeak,DataPeakError,mcSigma,mcSigmaError, dataSigma,dataSigmaError,result,doPUreweighting,"");
149	–	cout << " " << MCPeak << " +/- " << MCPeakError << endl;
150	–	cout << " " << DataPeak << " +/- " << DataPeakError << endl;
151	–
152	–	gMCPeak->SetPoint(i,0.5*(NumVtxBin[i]+NumVtxBin[i+1]),MCPeak);
153	–	gMCPeak->SetPointError(i,0.5*(NumVtxBin[i+1]-NumVtxBin[i]),MCPeakError);
154	–
155	–	gDataPeak->SetPoint(i,0.5*(NumVtxBin[i]+NumVtxBin[i+1]),DataPeak);
156	–	gDataPeak->SetPointError(i,0.5*(NumVtxBin[i+1]-NumVtxBin[i]),DataPeakError);
157	–
158	–	gMCWidth->SetPoint(i,0.5*(NumVtxBin[i]+NumVtxBin[i+1]),mcSigma);
159	–	gMCWidth->SetPointError(i,0.5*(NumVtxBin[i+1]-NumVtxBin[i]),mcSigmaError);
160	–
161	–	gDataWidth->SetPoint(i,0.5*(NumVtxBin[i]+NumVtxBin[i+1]),dataSigma);
162	–	gDataWidth->SetPointError(i,0.5*(NumVtxBin[i+1]-NumVtxBin[i]),dataSigmaError);
163	–	}
164	–
165	–	can->cd();
166	–	gMCPeak->GetXaxis()->SetTitle("N(Vertices)");
167	–	gMCPeak->GetYaxis()->SetTitle("Peak position");
168	–	gMCPeak->GetXaxis()->CenterTitle();
169	–	gMCPeak->GetYaxis()->CenterTitle();
170	–	gDataPeak->GetXaxis()->SetTitle("N(Vertices)");
171	–	gDataPeak->GetYaxis()->SetTitle("Peak position");
172	–	gDataPeak->GetXaxis()->CenterTitle();
173	–	gDataPeak->GetYaxis()->CenterTitle();
174	–
175	–	gDataWidth->GetXaxis()->SetTitle("N(Vertices)");
176	–	gDataWidth->GetYaxis()->SetTitle("#sigma_{JZB}");
177	–	gDataWidth->GetXaxis()->CenterTitle();
178	–	gDataWidth->GetYaxis()->CenterTitle();
179	–	gMCWidth->GetXaxis()->SetTitle("N(Vertices)");
180	–	gMCWidth->GetYaxis()->SetTitle("#sigma_{JZB}");
181	–	gMCWidth->GetXaxis()->CenterTitle();
182	–	gMCWidth->GetYaxis()->CenterTitle();
183	–
184	–	gMCPeak->SetFillColor(TColor::GetColor("#2E9AFE"));
185	–	gDataPeak->SetFillColor(TColor::GetColor("#2E9AFE"));
186	–	gDataWidth->SetFillColor(TColor::GetColor("#2E9AFE"));
187	–	gMCWidth->SetFillColor(TColor::GetColor("#2E9AFE"));
188	–
189	–	gMCPeak->Draw("A2");
190	–	DrawMCPrelim();
191	–	CompleteSave(can,"PUStudy/MCPeak");
192	–	can->cd();
193	–	gMCWidth->Draw("A2");
194	–	DrawMCPrelim();
195	–	CompleteSave(can,"PUStudy/MCWidth");
196	–
197	–	can->cd();
198	–	gDataPeak->Draw("A2");
199	–
200	–	DrawPrelim();
201	–	CompleteSave(can,"PUStudy/DataPeak");
202	–	can->cd();
203	–	gDataWidth->Draw("A2");
204	–	DrawPrelim();
205	–	CompleteSave(can,"PUStudy/DataWidth");
206	–
207	–	if(StoreResultsInSeparateFile) {
208	–	TFile *PeakFinding = new TFile("PeakFindingWithPeakCorrection.root","RECREATE");
209	–	gMCPeak->Write();
210	–	gDataPeak->Write();
211	–	gMCWidth->Write();
212	–	gDataWidth->Write();
213	–	PeakFinding->Close();
214	–	}
215	–
216	–	gDataPeak->Fit("pol1","L");
217	–	TF1 DataFit = (TF1)gDataPeak->GetFunction("pol1");
218	–	gMCPeak->Fit("pol1","L");
219	–	TF1 MCFit = (TF1)gMCPeak->GetFunction("pol1");
220	–
221	–
222	–	can->cd();
223	–	gMCPeak->Draw("A2");
224	–	DrawMCPrelim();
225	–	CompleteSave(can,"PUStudy/MCPeak_WithFit");
226	–	can->cd();
227	–
228	–	can->cd();
229	–	gDataPeak->Draw("A2");
230	–	DrawPrelim();
231	–	CompleteSave(can,"PUStudy/DataPeak_WithFit");
232	–
233	–
234	–
235	–	dout << "The suggested correction for Data is : " << DataFit->GetParameter(1) << endl;
236	–	dout << "The suggested correction for MC is : " << MCFit->GetParameter(1) << endl;
237	–
238	–	TLegend *leg = make_legend();
239	–	can->cd();
240	–	int j=0;
241	–	for(int i=PeakLibrary::DataHistogram.size()-5;i<PeakLibrary::DataHistogram.size();i++) {
242	–	PeakLibrary::DataHistogram[i]->Rebin(5);
243	–	PeakLibrary::DataHistogram[i]->SetLineColor(j+1);
244	–	if(i==PeakLibrary::DataHistogram.size()-5) PeakLibrary::DataHistogram[i]->DrawNormalized("histo");
245	–	else PeakLibrary::DataHistogram[i]->DrawNormalized("histo,same");
246	–	stringstream name;
247	–	name << NumVtxBin[j] << " < N_{vtx} < " << NumVtxBin[j+1];
248	–	leg->AddEntry(PeakLibrary::DataHistogram[i],name.str().c_str(),"l");
249	–	j++;
250	–	}
251	–
252	–	leg->Draw();
253	–	DrawPrelim();
254	–	CompleteSave(can,"PUStudy/DataOverview");
255	–
256	–	TLegend *leg2 = make_legend();
257	–	can->cd();
258	–	j=0;
259	–	for(int i=PeakLibrary::MCHistogram.size()-5;i<PeakLibrary::MCHistogram.size();i++) {
260	–	PeakLibrary::MCHistogram[i]->Rebin(5);
261	–	PeakLibrary::MCHistogram[i]->SetLineColor(j+1);
262	–	if(i==PeakLibrary::MCHistogram.size()-5) PeakLibrary::MCHistogram[i]->DrawNormalized("histo");
263	–	else PeakLibrary::MCHistogram[i]->DrawNormalized("histo,same");
264	–	stringstream name;
265	–	name << NumVtxBin[j] << " < N_{vtx} < " << NumVtxBin[j+1];
266	–	leg->AddEntry(PeakLibrary::MCHistogram[i],name.str().c_str(),"l");
267	–	j++;
268	–	}
269	–
270	–	leg->Draw();
271	–	DrawPrelim();
272	–	CompleteSave(can,"PUStudy/MCOverview");
273	–
274	–	delete can;
275	–	PeakLibrary::StoreHistos=false;
276	–	}
112
113		void FindPeakPileUpCorrection() {
114		bool StoreResultsInSeparateFile=true;
#	Line 290 \| Line 125 \| void FindPeakPileUpCorrection() {
125		TCanvas *can = new TCanvas("can","can");
126		float mcSigma,mcSigmaError,dataSigma,dataSigmaError;
127
128	<	const int nbins=11;
129	<	float NumVtxBin[nbins] = {0,7,9,11,13,15,17,19,21,24,50};
128	>	const int nbins=10;
129	>	float NumVtxBin[nbins] = {0,7,9,11,13,15,17,19,21,50};
130
131		stringstream NowCut;
132
#	Line 432 \| Line 267 \| void FindPeakPileUpCorrection() {
267		gDataWidth->SetFillColor(TColor::GetColor("#2E9AFE"));
268		gMCWidth->SetFillColor(TColor::GetColor("#2E9AFE"));
269
270	+	can->cd();
271		gMCPeak->Draw("A2");
272		DrawMCPrelim();
273		CompleteSave(can,"PUStudy/MCPeak");
274	+
275		can->cd();
276	<	gMCWidth->Draw("A2");
276	>	gStyle->SetOptFit(0);
277	>	gMCPeak->Fit("pol1");
278	>	gMCPeak->GetYaxis()->SetRangeUser(0,15);
279	>	TF1 MCFit = (TF1)gMCPeak->GetFunction("pol1");
280	>
281	>	gMCPeak->Draw("A2");
282	>	MCFit->SetLineColor(kBlue);
283	>	MCFit->SetLineWidth(2);
284	>	MCFit->Draw("same");
285	>	stringstream mresult;
286	>	mresult << "Slope: " << std::setprecision(2) << MCFit->GetParameter(1) << " +/- " << MCFit->GetParError(1);
287	>	TText* mrestitle = write_text(0.7,0.2,mresult.str());
288	>	mrestitle->Draw();
289	>	DrawMCPrelim();
290	>	CompleteSave(can,"PUStudy/MCPeak_Fit");
291	>
292	>	gMCPeak->Draw("A2");
293		DrawMCPrelim();
294		CompleteSave(can,"PUStudy/MCWidth");
295
296	+
297		can->cd();
298		gDataPeak->Draw("A2");
445	–
299		DrawPrelim();
300		CompleteSave(can,"PUStudy/DataPeak");
301	+
302	+	can->cd();
303	+	gStyle->SetOptFit(0);
304	+	gDataPeak->GetYaxis()->SetRangeUser(0,15);
305	+	gDataPeak->Fit("pol1");
306	+
307	+	TF1 DataFit = (TF1)gDataPeak->GetFunction("pol1");
308	+	DataFit->SetLineColor(kBlue);
309	+	DataFit->SetLineWidth(2);
310	+	DataFit->Draw("same");
311	+	stringstream dresult;
312	+	dresult << "Slope: " << std::setprecision(2) << DataFit->GetParameter(1) << " +/- " << DataFit->GetParError(1);
313	+	TText* drestitle = write_text(0.7,0.2,dresult.str());
314	+	drestitle->Draw();
315	+	DrawPrelim();
316	+	CompleteSave(can,"PUStudy/DataPeak_Fit");
317	+
318		can->cd();
319		gDataWidth->Draw("A2");
320		DrawPrelim();
#	Line 459 \| Line 329 \| void FindPeakPileUpCorrection() {
329		PeakFinding->Close();
330		}
331
462	–	gDataPeak->Fit("pol1");
463	–	TF1 DataFit = (TF1)gDataPeak->GetFunction("pol1");
464	–	gMCPeak->Fit("pol1");
465	–	TF1 MCFit = (TF1)gMCPeak->GetFunction("pol1");
332
333		dout << "The suggested correction for Data is : " << DataFit->GetParameter(1) << endl;
334		dout << "The suggested correction for MC is : " << MCFit->GetParameter(1) << endl;
#	Line 755 \| Line 621 \| void make_special_obs_pred_mll_plot(stri
621		dleg->SetY1(0.7);
622
623
624	+	diff->GetYaxis()->SetTitle("true - predicted");
625	+	if(MCorData==data) diff->GetYaxis()->SetTitle("observed - predicted");
626		diff->Draw("e1");
627		BpredSysBand->Draw("e2");
628		diff->Draw("e1,same");
#	Line 1038 \| Line 906 \| float make_data_comparison_plot(string v
906		}
907
908		void make_OFSF_plots(string variable, string addcut, int nbins, float min, float max, bool logscale,
909	<	string xlabel, string filename, float legendPosition=0.55) {
909	>	string xlabel, string filename, float tymax=-1, float legendPosition=0.55) {
910
911		string mllcuts[] = { "mll>20","mll>15&&mll<70", "mll>75&&mll<105", "mll>120" };
912		string mllcutname[] = { "m_{ll} > 20 GeV", "20 < m_{ll} < 70 GeV", "70 < m_{ll} < 110 GeV", "m_{ll} > 120 GeV" };
913		string plotname[] = {"_all","_low","_peak","_high"};
914		float ymax;
915	+	if(tymax>0) ymax=tymax;
916
917		int start = 0;
918		if ( !PlottingSetup::openBox ) start = 3;
919
920		for ( int i=start; i<4; ++i ) {
921		if ( addcut != "" ) mllcuts[i] += "&&"+addcut;
922	<	if ( i==start ) {
922	>	if ( i==start && !(tymax>0) ) {
923		ymax = make_one_OFSF_plot(variable, mllcuts[i], mllcutname[i], nbins, min, max, -1, logscale, xlabel,
924		filename+plotname[i], legendPosition );
925		} else {
#	Line 2438 \| Line 2307 \| void ProduceOFSFPlots(string mcjzb, stri
2307		Prefix="Aachen/";
2308		}
2309
2310	<	make_OFSF_plots("met[4]", BaseMetCut, 30, 100., 400., true, "E_{T}^{miss}", Prefix+"MET");
2311	<	make_OFSF_plots("pfJetGoodNumBtag30", BaseMetCut, 5, -0.5, 4.5, true, "N_{b-jets}", Prefix+"nbjets");
2312	<	make_OFSF_plots("st", BaseMetCut, 20, 0, 2000., false, "#tilde{S}_{T}", Prefix+"FrankT");
2313	<	make_OFSF_plots("pt1+pt2+met[4]+Sum$(pfJetGoodPt)", BaseMetCut, 20, 0, 2000., false, "#tilde{S}_{T}", Prefix+"POOR_FrankT");
2314	<	/* make_OFSF_plots("mt2j", BaseMetCut, 35, 0, 700., false, "M_{T2J}", Prefix+"MT2J");
2315	<
2316	<	make_OFSF_plots("ml1b", BaseMetCut, 20, 0, 800., false, "m_{l_{1}b}", Prefix+"ML1B");
2317	<	make_OFSF_plots("ml2b", BaseMetCut, 20, 0, 400., false, "m_{l_{2}b}", Prefix+"ML2B");
2318	<	make_OFSF_plots("min(ml1b,ml2b)", BaseMetCut, 20, 0, 400., false, "m_{lb}", Prefix+"MLB");
2319	<
2320	<	make_OFSF_plots(datajzb, BaseMetCut, 55, -100, 450., false, "JZB", Prefix+"JZB_In_Signal_Region");
2321	<	make_OFSF_plots("mt2", BaseMetCut, 15, 0, 150., false, "M_{T2}", Prefix+"MT2");
2322	<	make_OFSF_plots("d2", BaseMetCut, 20, 0, 200., true, "D2", Prefix+"D2");
2323	<	make_OFSF_plots("pfJetGoodNum40", BaseMetCut, 10, -0.5, 9.5, true, "N_{jets}", Prefix+"nJets");
2324	<	*/
2325	<	make_OFSF_plots("mll", BaseMetCut, 60, 20., 320., false, "m_{ll}", Prefix+"mll");
2326	<	make_OFSF_plots("leptonNum", BaseMetCut, 3, 2, 5., false, "N(leptons)", Prefix+"NLeptons");
2327	<
2328	<	make_OFSF_plots("pfJetGoodNum40", BaseMetCut, 7, 3, 10, true, "#(jets)", Prefix+"njets");
2329	<	// make_OFSF_plots("pfJetGoodNum40", BaseMetCut+"&&pfJetGoodNumBtag30==0", 7, 3, 10, true, "#(jets)", Prefix+"njets_btagVeto");
2330	<	// make_OFSF_plots("pfJetGoodNum40", BaseMetCut"&&pfJetGoodNumBtag30>0", 7, 3, 10, true, "#(jets)", Prefix+"njets_AtLeastOneBJet30");
2331	<
2332	<	make_OFSF_plots("pfJetGoodPtBtag[0]", BaseMetCut+"&&pfJetGoodNumBtag30>0", 20, 0, 400, true, "p_{T}(leading b-jet)", Prefix+"ptb1");
2333	<	make_OFSF_plots("pfJetGoodPt", BaseMetCut, 30, 0., 400., true, "p_{T}^{J}", Prefix+"ptJet");
2334	<
2335	<	// make_OFSF_plots("iso1", BaseMetCut, 20, 0, 0.3, true, "lepton 1 isolation", Prefix+"iso1");
2336	<	// make_OFSF_plots("iso2", BaseMetCut, 20, 0, 0.3, true, "lepton 2 isolation", Prefix+"iso2");
2337	<	make_OFSF_plots("pt", BaseMetCut, 30, 0., 300., true, "p_{ll}", Prefix+"pt");
2338	<	make_OFSF_plots("pt1", BaseMetCut, 30, 0., 300., true, "p_{T,1}", Prefix+"pt1");
2339	<	make_OFSF_plots("pt2", BaseMetCut, 22, 0., 220., true, "p_{T,2}", Prefix+"pt2");
2340	<	// make_OFSF_plots("eta1", BaseMetCut, 10, -2.5, 2.5, false, "#eta_{1}", Prefix+"eta1", 0.15);
2341	<	// make_OFSF_plots("eta2", BaseMetCut, 10, -2.5, 2.5, false, "#eta_{2}", Prefix+"eta2", 0.15);
2342	<	// make_OFSF_plots("phi1", BaseMetCut, 10, -TMath::Pi(), TMath::Pi(), false, "#phi_{1}", Prefix+"phi1", 0.2);
2343	<	// make_OFSF_plots("phi2", BaseMetCut, 10, -TMath::Pi(), TMath::Pi(), false, "#phi_{2}", Prefix+"phi2", 0.2);
2344	<	// make_OFSF_plots("pfJetGoodPt[0]/pfJetGoodPt[1]", BaseMetCut, 20, 1, 10, true, "pt_{j}^{1}/pt_{j}^{2}", Prefix+"jpt1pt2", 0.2);
2345	<	// make_OFSF_plots("TMath::Abs(pfJetDphiMet[0])", BaseMetCut, 16, 0, 3.2, false, "\|#Delta#phi(jet1,MET)\|", Prefix+"dphij1met", 0.2);
2310	>	// make_OFSF_plots("met[4]", BaseMetCut, 30, 100., 400., true, "E_{T}^{miss}", Prefix+"MET",1500);
2311	>	// make_OFSF_plots("pfJetGoodNumBtag30", BaseMetCut, 6, -0.5, 5.5, true, "N_{b-jets}", Prefix+"nbjets",5000);
2312	>	make_OFSF_plots("max(eta1,eta2)", BaseMetCut, 20, -2.4, 2.4, true, "max(#eta_{l1},#eta_{l2})", Prefix+"MaxEta",-1);
2313	>	// make_OFSF_plots("st", BaseMetCut, 20, 0, 2000., false, "#tilde{S}_{T}", Prefix+"FrankT",120);
2314	>	// make_OFSF_plots("pt1+pt2+met[4]+Sum$(pfJetGoodPt)", BaseMetCut, 20, 0, 2000., false, "#tilde{S}_{T}", Prefix+"POOR_FrankT",120);
2315	>
2316	>
2317	>	// make_OFSF_plots("mt2j", BaseMetCut, 35, 0, 700., false, "M_{T2J}", Prefix+"MT2J");
2318	>	//
2319	>	// make_OFSF_plots("ml1b", BaseMetCut, 20, 0, 800., false, "m_{l_{1}b}", Prefix+"ML1B");
2320	>	// make_OFSF_plots("ml2b", BaseMetCut, 20, 0, 400., false, "m_{l_{2}b}", Prefix+"ML2B");
2321	>	// make_OFSF_plots("min(ml1b,ml2b)", BaseMetCut, 20, 0, 400., false, "m_{lb}", Prefix+"MLB");
2322	>	//
2323	>	// make_OFSF_plots(datajzb, BaseMetCut, 55, -100, 450., false, "JZB", Prefix+"JZB_In_Signal_Region");
2324	>	// make_OFSF_plots("mt2", BaseMetCut, 15, 0, 150., false, "M_{T2}", Prefix+"MT2");
2325	>	// make_OFSF_plots("d2", BaseMetCut, 20, 0, 200., true, "D2", Prefix+"D2");
2326	>	make_OFSF_plots("pfJetGoodNum40", BaseMetCut, 10, -0.5, 9.5, true, "N_{jets}", Prefix+"nJets");
2327	>	//
2328	>	// make_OFSF_plots("mll", BaseMetCut, 60, 20., 320., false, "m_{ll}", Prefix+"mll");
2329	>	// make_OFSF_plots("leptonNum", BaseMetCut, 3, 2, 5., false, "N(leptons)", Prefix+"NLeptons");
2330	>	//
2331	>	// make_OFSF_plots("pfJetGoodNum40", BaseMetCut, 7, 3, 10, true, "#(jets)", Prefix+"njets");
2332	>	// make_OFSF_plots("pfJetGoodNum40", BaseMetCut+"&&pfJetGoodNumBtag30==0", 7, 3, 10, true, "#(jets)", Prefix+"njets_btagVeto");
2333	>	// make_OFSF_plots("pfJetGoodNum40", BaseMetCut"&&pfJetGoodNumBtag30>0", 7, 3, 10, true, "#(jets)", Prefix+"njets_AtLeastOneBJet30");
2334	>	//
2335	>	// make_OFSF_plots("pfJetGoodPtBtag[0]", BaseMetCut+"&&pfJetGoodNumBtag30>0", 20, 0, 400, true, "p_{T}(leading b-jet)", Prefix+"ptb1");
2336	>	make_OFSF_plots("pfJetGoodPt", BaseMetCut, 30, 0., 400., true, "p_{T}^{J}", Prefix+"ptJet");
2337	>	//
2338	>	// make_OFSF_plots("iso1", BaseMetCut, 20, 0, 0.3, true, "lepton 1 isolation", Prefix+"iso1");
2339	>	// make_OFSF_plots("iso2", BaseMetCut, 20, 0, 0.3, true, "lepton 2 isolation", Prefix+"iso2");
2340	>	make_OFSF_plots("pt", BaseMetCut, 30, 0., 300., true, "p_{ll}", Prefix+"pt");
2341	>	make_OFSF_plots("pt1", BaseMetCut, 30, 0., 300., true, "p_{T,1}", Prefix+"pt1");
2342	>	make_OFSF_plots("pt2", BaseMetCut, 22, 0., 220., true, "p_{T,2}", Prefix+"pt2");
2343	>	// make_OFSF_plots("eta1", BaseMetCut, 10, -2.5, 2.5, false, "#eta_{1}", Prefix+"eta1", -1, 0.15);
2344	>	// make_OFSF_plots("eta2", BaseMetCut, 10, -2.5, 2.5, false, "#eta_{2}", Prefix+"eta2", -1, 0.15);
2345	>	// make_OFSF_plots("phi1", BaseMetCut, 10, -TMath::Pi(), TMath::Pi(), false, "#phi_{1}", Prefix+"phi1", -1, 0.2);
2346	>	// make_OFSF_plots("phi2", BaseMetCut, 10, -TMath::Pi(), TMath::Pi(), false, "#phi_{2}", Prefix+"phi2", -1, 0.2);
2347	>	// make_OFSF_plots("pfJetGoodPt[0]/pfJetGoodPt[1]", BaseMetCut, 20, 1, 10, true, "pt_{j}^{1}/pt_{j}^{2}", Prefix+"jpt1pt2", -1, 0.2);
2348	>	make_OFSF_plots("TMath::Abs(pfJetDphiMet[0])", BaseMetCut, 16, 0, 3.2, false, "\|#Delta#phi(jet1,MET)\|", Prefix+"dphij1met", -1, 0.2);
2349
2350	<	// make_OFSF_plots("TMath::Abs(dphi)", BaseMetCut, 16, 0, 3.2, false, "\|#Delta#phi(l1,l2)\|", Prefix+"dphi", 0.2);
2351	<	// make_OFSF_plots("TMath::Abs(dphiMet1)", BaseMetCut, 16, 0, 3.2, false, "\|#Delta#phi(l1,MET)\|", Prefix+"dphiMet1", 0.2);
2352	<	// make_OFSF_plots("TMath::Abs(dphiMet2)", BaseMetCut, 16, 0, 3.2, false, "\|#Delta#phi(l2,MET)\|", Prefix+"dphiMet2", 0.2);
2350	>	// make_OFSF_plots("TMath::Abs(dphi)", BaseMetCut, 16, 0, 3.2, false, "\|#Delta#phi(l1,l2)\|", Prefix+"dphi", -1, 0.2);
2351	>	// make_OFSF_plots("TMath::Abs(dphiMet1)", BaseMetCut, 16, 0, 3.2, false, "\|#Delta#phi(l1,MET)\|", Prefix+"dphiMet1", -1, 0.2);
2352	>	// make_OFSF_plots("TMath::Abs(dphiMet2)", BaseMetCut, 16, 0, 3.2, false, "\|#Delta#phi(l2,MET)\|", Prefix+"dphiMet2", -1, 0.2);
2353		// make_OFSF_plots("TMath::Min(TMath::Abs(dphiMet1), TMath::Abs(dphiMet2))", BaseMetCut, 16, 0, 3.2, false, "Min(\|#Delta#phi(l,MET)\|)", Prefix+"dphilc");
2354		// make_OFSF_plots("TMath::Min(TMath::Abs(pfJetDphiMet[0]), TMath::Min(TMath::Abs(pfJetDphiMet[1]), TMath::Abs(pfJetDphiMet[2])))", BaseMetCut, 16, 0, 3.2, false, "Min(\|#Delta#phi(jet,MET)\|)", Prefix+"dphijc");
2355		// make_OFSF_plots("TMath::Min((TMath::Pi()-TMath::Abs(dphiMet1)), (TMath::Pi() - TMath::Abs(dphiMet2)))", BaseMetCut, 16, 0, 3.2, false, "Min(#pi - \|#Delta#phi(l,MET)\|)", Prefix+"dphilco");
#	Line 2515 \| Line 2387 \| void do_kinematic_plots(string mcjzb, st
2387		bool dolog=true;
2388		bool nolog=false;
2389
2390	<	bool doOFSF = true;
2390	>	bool doOFSF = false;
2391		bool doKin = false;
2392		bool doDataComp = false;
2393		bool MakeTwoThreeJetComparison = false;
#	Line 6373 \| Line 6245 \| void make_SF_over_OF_plot(string datajzb
6245		leg->AddEntry(gra,"MC","fp");
6246		leg->Draw();
6247
6248	<	TLine *g = new TLine(0,1.01,200,1.01);
6248	>	TLine *g = new TLine(0,1.0,200,1.0);
6249		g->SetLineColor(kBlue);
6250		g->Draw();
6251

Diff Legend

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

Comparing UserCode/cbrown/Development/Plotting/Modules/Plotting_Functions.C (file contents): Revision 1.97 by buchmann, Thu May 30 07:59:33 2013 UTC vs. Revision 1.98 by buchmann, Fri May 31 09:45:15 2013 UTC

Diff Legend

Comparing UserCode/cbrown/Development/Plotting/Modules/Plotting_Functions.C (file contents):
Revision 1.97 by buchmann, Thu May 30 07:59:33 2013 UTC vs.
Revision 1.98 by buchmann, Fri May 31 09:45:15 2013 UTC