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.3 by buchmann, Thu Jul 14 15:50:58 2011 UTC vs.
Revision 1.5 by buchmann, Tue Jul 19 07:06:51 2011 UTC

#	Line 48 | Line 48 | int main()
48		gStyle->SetTextFont(42);
49		bool showList=false;
50		set_treename("events");//you can set the treename here to be used; options are "events" (for reco) for "PFevents" (for particle flow)
51	<	define_samples(showList,allsamples);
51	>	define_samples(showList,allsamples,signalsamples);
52		setlumi(luminosity);
53		do_png(true);
54		do_pdf(false);
#	Line 62 | Line 62 | int main()
62		bool do_all=false;/// DONE
63		bool do_peak_finding=false; /// DONE
64		bool calculate_pred_and_observed=false; /// DONE
65	<	bool study_sidebands=true; /// FOR REAL
65	>	bool study_sidebands=false; /// FOR REAL
66		bool do_test=false;/// DONE    (just any test you wish)
67		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
68		bool do_find_sideband_definition=false;
69	+	bool do_pick_up_events=true; ///DONE
70
71		//**** part 1 : peak finding
72		float MCPeak=0,MCPeakError=0,DataPeak=0,DataPeakError=0,MCSigma=10,DataSigma=10;
#	Line 97 | Line 98 | int main()
98
99		if(study_sidebands) look_at_sidebands(mcjzb.str(),datajzb.str());
100
101	+	if(do_pick_up_events) {
102	+
103	+	cout << "Observed: " << endl;
104	+	pick_up_events((const char*)(cutmass&&cutOSSF&&cutnJets&&basiccut&&"((jzb[1]+0.06*pt-2.84727)>100)"));
105	+	cout << "Predicted (JZB<-100) OSSF" << endl;
106	+	pick_up_events((const char*)(cutmass&&cutOSSF&&cutnJets&&basiccut&&"((jzb[1]+0.06*pt-2.84727)<-100)"));
107	+	cout << "Predicted (emu, JZB>100) OSOF" << endl;
108	+	pick_up_events((const char*)(cutmass&&cutOSOF&&cutnJets&&basiccut&&"((jzb[1]+0.06*pt-2.84727)>100)"));
109	+	cout << "Predicted (emu, JZB<-100) OSOF" << endl;
110	+	pick_up_events((const char*)(cutmass&&cutOSOF&&cutnJets&&basiccut&&"((jzb[1]+0.06*pt-2.84727)<-100)"));
111	+	cout << "Predicted (SB emu, JZB>100) OSOF" << endl;
112	+	pick_up_events((const char*)(sidebandcut&&cutOSOF&&cutnJets&&basiccut&&"((jzb[1]+0.06*pt-2.84727)>100)"));
113	+	cout << "Predicted (SB emu, JZB<-100) OSOF" << endl;
114	+	pick_up_events((const char*)(cutOSOF&&cutnJets&&basiccut&&sidebandcut&&"((jzb[1]+0.06*pt-2.84727)<-100)"));
115	+	cout << "Predicted (SB SF, JZB>100) OSSF" << endl;
116	+	pick_up_events((const char*)(sidebandcut&&cutOSSF&&cutnJets&&basiccut&&"((jzb[1]+0.06*pt-2.84727)>100)"));
117	+	cout << "Predicted (SB SF, JZB<-100) OSSF" << endl;
118	+	pick_up_events((const char*)(sidebandcut&&cutOSSF&&cutnJets&&basiccut&&"((jzb[1]+0.06*pt-2.84727)<-100)"));
119	+	}
120	+
121		if(do_test) test();
122
123

Diff Legend

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