ViewVC Help

View File | Revision Log | Show Annotations | Root Listing

root/cvsroot/UserCode/cbrown/Development/Plotting/Modules/GeneralToolBox.C

Comparing UserCode/cbrown/Development/Plotting/Modules/GeneralToolBox.C (file contents):
Revision 1.3 by buchmann, Fri Mar 2 08:37:58 2012 UTC vs.
Revision 1.41 by buchmann, Fri Apr 26 07:00:24 2013 UTC

#	Line 17 | Line 17
17		#include <netdb.h>
18		#include <stdio.h>
19		#include <string.h>
20	+	#include <assert.h>
21
22		#include <TFile.h>
23		#include <TTree.h>
#	Line 38 | Line 39
39		#include <TPaveText.h>
40		#include <TRandom.h>
41		#include <TGraphErrors.h>
42	+	#include <TROOT.h>
43	+	#include <TPolyLine.h>
44	+
45		#ifndef Verbosity
46		#define Verbosity 0
47		#endif
#	Line 58 | Line 62 | namespace PlottingSetup {
62		vector<float> global_ratio_binning;
63		int publicmode=0;
64		int PaperMode=0; // PaperMode=true will suppress "Preliminary" in DrawPrelim()
65	+	int Approved=0; // Approved=true will only plot approved plots
66	+	bool is2012=true;
67	+	bool is53reco=true;
68	+	bool openBox = true;
69	+	vector<TH1F*> FakeHistoHeap;
70	+	bool DrawMetSignalRegionMllLines=false; // whether to draw the lines in mll plots for the MET signal region
71	+	float ConsideredZWidth=20;
72		}
73
74		bool dopng=false;
#	Line 84 | Line 95 | void write_error(string funcname, string
95		void write_info(string funcname, string text);
96		string get_directory();
97		bool Contains(string wholestring, string findme);
98	+	TH1F* CollapseStack(THStack stack,TString hname);
99		//-------------------------------------------------------------------------------------
100
101		template<typename U>
#	Line 119 | Line 131 | string topdir(string child) {
131		return tempdirectory.substr(0,ichar);
132		}
133		}
134	+	return "TOPDIRFAILURE";
135		}
136
137		template < typename CHAR_TYPE,
#	Line 250 | Line 263 | void ensure_directory_exists(string this
263		}
264		}
265
266	+	void DeleteStack(THStack mcstack) {
267	+	TH1F* h;
268	+	TIter nextHist(mcstack.GetHists());
269	+	while ( h = (TH1F*)nextHist() ) {
270	+	h->Delete();
271	+	}
272	+	delete h;
273	+	}
274	+
275	+
276	+	void DeadEnd(string message) {
277	+	dout << endl;
278	+	dout << "\033[1;31m     _                _                  _  " << endl;
279	+	dout << "\033[1;31m  __| | ___  __ _  __| |   ___ _ __   __| | " << endl;
280	+	dout << "\033[1;31m / _` |/ _ \\/ _` |/ _` |  / _ \\ '_ \\ / _` | " << endl;
281	+	dout << "\033[1;31m| (_| |  __/ (_| | (_| | |  __/ | | | (_| | " << endl;
282	+	dout << "\033[1;31m \\__,_|\\___|\\__,_|\\__,_|  \\___|_| |_|\\__,_| " << endl;
283	+	dout << "\033[1;31m                                            " << endl;
284	+	dout << "" << endl;
285	+	dout << " A totally fatal error has occurred: " << endl;
286	+	dout << " " << message << endl;
287	+	dout << endl;
288	+	dout << "If you do not know how to deal with this error message please call Customer Support Mumbai (or Zurich...) \033[0m" << endl;
289	+	dout << "   ... aborting now ... " << endl;
290	+	assert(0);
291	+	}
292	+
293	+
294	+	void SanityChecks() {
295	+	//this function gets called whenever samples are initialized (see ActiveSamples.C)
296	+	//whatever is necessary for the code to run should be put here.
297	+
298	+	if(gROOT->GetVersionInt()<53200) {
299	+	DeadEnd("ROOT version outdated, you are using "+any2string(gROOT->GetVersionInt())+" but the earliest recommended version is 5.32.00");
300	+	}
301	+
302	+
303	+	}
304	+
305		void initialize_log() {
306		if(!PlottingSetup::publicmode) {
307		dout << "____________________________________________________________" << endl;
#	Line 278 | Line 330 | void initialize_log() {
330		dout << "        PUBLIC MODE     " << endl;
331		}
332		time_t rawtime;
281	–	struct tm * timeinfo;
333		time ( &rawtime );
334		dout << "       Analysis run on " << asctime (localtime ( &rawtime ));
335		dout << "____________________________________________________________" << endl;
#	Line 289 | Line 340 | void extract_cbaf_dir(string curpath) {
340		int position=curpath.find("/Plotting");
341		if(position<0) position=curpath.find("/DistributedModelCalculations");
342		if(position<0) position=curpath.find("/various_assignments");
343	+	if(position<0) position=curpath.find("/exchange");
344		PlottingSetup::cbafbasedir=curpath.substr(0,position);
345		}
346
#	Line 296 | Line 348 | void set_directory(string basedir="") {
348		if(basedir.substr(0,1)=="/") basedir=basedir.substr(1,basedir.length()-1);
349		if(basedir.substr(basedir.length()-1,1)!="/") basedir+="/";
350		char currentpath[1024];
351	<	char *path = getcwd(currentpath,1024);
351	>	getcwd(currentpath,1024);
352		PlottingSetup::basedirectory=(string)currentpath+"/Plots/"+basedir;
353		ensure_directory_exists(PlottingSetup::basedirectory);
354		initialize_log();
#	Line 321 | Line 373 | string extract_directory(string savethis
373		string extract_root_dir(string name) {
374		int position = -1;
375		if(name.substr(0,1)=="/") name=name.substr(1,name.length()-1);
376	<	for(int ipos=0;ipos<name.length();ipos++) {
376	>	for(int ipos=0;ipos<(int)name.length();ipos++) {
377		if(name.substr(ipos,1)=="/") position=ipos;
378		}
379		if(position==-1) return "";
#	Line 331 | Line 383 | string extract_root_dir(string name) {
383		string extract_root_filename(string name) {
384		int position = -1;
385		if(name.substr(0,1)=="/") name=name.substr(1,name.length()-1);
386	<	for(int ipos=0;ipos<name.length();ipos++) {
386	>	for(int ipos=0;ipos<(int)name.length();ipos++) {
387		if(name.substr(ipos,1)=="/") position=ipos;
388		}
389		return name.substr(position+1,name.length()-position-1);
#	Line 389 | Line 441 | void CompleteSave(TVirtualPad *can, stri
441
442
443		void setlumi(float l) {
444	<	generaltoolboxlumi=l;
444	>	generaltoolboxlumi=l;
445		}
446
447		int write_first_line(vector<vector<string> > &entries) {
#	Line 414 | Line 466 | void write_entry(string entry,int width,
466		int currwidth=entry.size();
467		while(currwidth<width) {
468		entry=" "+entry;
469	<	if(entry.size()<width) entry=entry+" ";
469	>	if((int)entry.size()<width) entry=entry+" ";
470		currwidth=entry.size();
471		}
472		bool do_special=false;
#	Line 428 | Line 480 | void write_entry(string entry,int width,
480		void make_nice_table(vector<vector <string> > &entries) {
481		int ncolumns=write_first_line(entries);
482		int cellwidth=(int)(((float)(60-(ncolumns+1)))/(ncolumns));
483	<	for(int iline=0;iline<entries.size();iline++) {
483	>	for(int iline=0;iline<(int)entries.size();iline++) {
484		vector<string> currline = entries[iline];
485		dout << " |";
486	<	for(int ientry=0;ientry<currline.size();ientry++) {
486	>	for(int ientry=0;ientry<(int)currline.size();ientry++) {
487		write_entry(currline[ientry],cellwidth);
488		}
489		dout << endl;
#	Line 443 | Line 495 | void make_nice_table(vector<vector <stri
495		void make_nice_jzb_table(vector<vector <string> > &entries) {
496		int ncolumns=write_first_line(entries);
497		int cellwidth=(int)(((float)(60-(ncolumns+1)))/(ncolumns));
498	<	for(int iline=0;iline<entries.size();iline++) {
498	>	for(int iline=0;iline<(int)entries.size();iline++) {
499		vector<string> currline = entries[iline];
500		dout << " |";
501	<	for(int ientry=0;ientry<currline.size();ientry++) {
501	>	for(int ientry=0;ientry<(int)currline.size();ientry++) {
502		write_entry(currline[ientry],cellwidth,iline,ientry);
503		}
504		dout << endl;
#	Line 547 | Line 599 | void DrawPrelim(float writelumi=generalt
599		string prel=" Preliminary";
600		if(PlottingSetup::PaperMode) prel="";
601		//prelimtext << "CMS Preliminary 2011 , #sqrt{s}= 7 TeV, L= O(1) fb^{-1}"; //temporary replacement
602	+	//  if(PlottingSetup::is53reco) prel += " 53X";
603	+	string energy="7 TeV";
604	+	if(PlottingSetup::is2012) energy="8 TeV";
605		if(writelumi == 0) {
606	<	if(isMC) prelimtext << "CMS Simulation, #sqrt{s} = 7 TeV";
607	<	else prelimtext << "CMS" << prel << ", #sqrt{s} = 7 TeV";
606	>	if(isMC) prelimtext << "CMS Simulation, #sqrt{s} = " << energy;
607	>	else prelimtext << "CMS" << prel << ", #sqrt{s} = " << energy;
608		} else {
609	<	if(isMC) prelimtext << "CMS Simulation, #sqrt{s} = 7 TeV, L_{int} = " << std::setprecision(2) <<writelumi<<" "<<barn<<"^{-1}";
610	<	else prelimtext << "CMS" << prel << ", #sqrt{s} = 7 TeV, L_{int} = " << std::setprecision(2) <<writelumi<<" "<<barn<<"^{-1}";
609	>	if(isMC) prelimtext << "CMS Simulation, #sqrt{s} = " << energy << ", #scale[0.6]{#int}L dt = " << std::setprecision(3) <<writelumi<<" "<<barn<<"^{-1}";
610	>	else prelimtext << "CMS" << prel << ", #sqrt{s} = " << energy << ", #scale[0.6]{#int}L dt = " << std::setprecision(3) <<writelumi<<" "<<barn<<"^{-1}";
611		}
612		TPaveText *eventSelectionPaveText = new TPaveText(0.27, 0.93,0.77, 1.0,"blNDC");
613		eventSelectionPaveText->SetFillStyle(4000);
614		eventSelectionPaveText->SetBorderSize(0);
615		eventSelectionPaveText->SetFillColor(kWhite);
616		eventSelectionPaveText->SetTextFont(42);
617	+	//eventSelectionPaveText->SetTextFont(62); // paper style
618		eventSelectionPaveText->SetTextSize(0.042);
619		eventSelectionPaveText->AddText(prelimtext.str().c_str());
620		eventSelectionPaveText->Draw();
#	Line 596 | Line 652 | TGraph* make_nice_ratio(int nbins,float
652		bincontent[i]=0;
653		errors[i]=0;
654		}
655	<	float currlimit=binning[0];
655	>	//  float currlimit=binning[0];
656		int currtoplim=1;
657		for(int ibin=1;ibin<=histo->GetNbinsX();ibin++)
658		{
#	Line 629 | Line 685 | float computeProductError(float a, float
685		return val;
686		}
687
688	+	TGraphAsymmErrors *SysAndStatRatio(TH1F *h1,TH1F *h2, TH1F *sys, int id, vector<float>binning, bool precise=false) {
689	+	int absJZBbinsNumber = binning.size()-1;
690	+	TGraphAsymmErrors* graph = new TGraphAsymmErrors(absJZBbinsNumber);
691	+	for(int i=0;i<absJZBbinsNumber;i++)
692	+	{
693	+	float xCenter=h1->GetBinCenter(i+1);
694	+	float xWidth=(h1->GetBinWidth(i+1))*0.5;
695	+	float nominatorError = h1->GetBinError(i+1);
696	+	float nominator=h1->GetBinContent(i+1);
697	+	float denominatorError=h2->GetBinError(i+1);
698	+	float denominator=h2->GetBinContent(i+1);
699	+	float errorN = 0;
700	+	float errorP = computeRatioError(nominator,nominatorError,denominator,denominatorError);
701	+	float syserr=sys->GetBinContent(i+1);
702	+	float errorsys=0;
703	+	if(id==1) // (is data)
704	+	{
705	+	if(!precise) errorP = computeRatioError(nominator,statErrorP(nominator),denominator,statErrorP(denominator));
706	+	else         errorP = computeRatioError(nominator,nominatorError,denominator,denominatorError);
707	+	errorsys=computeRatioError(nominator,syserr,denominator,0);
708	+	errorP = TMath::Sqrt(errorsys*errorsys+errorP*errorP);
709	+	errorN = errorP; // symmetrize using statErrorP
710	+	} else {
711	+	errorN = computeRatioError(nominator,nominatorError,denominator,denominatorError);
712	+	errorsys=computeRatioError(nominator,syserr,denominator,0);
713	+	errorN = TMath::Sqrt(errorsys*errorsys+errorP*errorP);
714	+	errorP = errorN; // symmetrize using statErrorP
715	+	}
716	+	if(denominator!=0) {
717	+	graph->SetPoint(i, xCenter, nominator/denominator);
718	+	graph->SetPointError(i,xWidth,xWidth,errorN,errorP);
719	+	}
720	+	else {
721	+	graph->SetPoint(i, xCenter, -999);
722	+	graph->SetPointError(i,xWidth,xWidth,errorN,errorP);
723	+	}
724	+	}
725	+	return graph;
726	+	}
727	+
728	+
729		TGraphAsymmErrors *histRatio(TH1F *h1,TH1F *h2, int id, vector<float>binning, bool precise=false)
730		{
731		int absJZBbinsNumber = binning.size()-1;
732		TGraphAsymmErrors* graph = new TGraphAsymmErrors(absJZBbinsNumber);
733
734	<	for(unsigned int i=0;i<absJZBbinsNumber;i++)
734	>	for(int i=0;i<absJZBbinsNumber;i++)
735		{
736		float xCenter=h1->GetBinCenter(i+1);
737		float xWidth=(h1->GetBinWidth(i+1))*0.5;
#	Line 746 | Line 843 | void splitupcut(string incut, vector<str
843		int paranthesis_open=0;
844		int substr_start=0;
845		string currchar="";
846	<	for (int ichar=0;ichar<incut.length();ichar++)
846	>	for (int ichar=0;ichar<(int)incut.length();ichar++)
847		{
848		currchar=incut.substr(ichar,1);
849		//        if(currchar=="(") paranthesis_open++;
#	Line 759 | Line 856 | void splitupcut(string incut, vector<str
856		partvector.push_back(incut.substr(substr_start,incut.length()-substr_start));
857		if(Verbosity>1) {
858		dout << "[ splitupcut() ] : The cut vector now contains the following elements: "<< endl;
859	<	for (int ipart=0;ipart<partvector.size();ipart++)
859	>	for (int ipart=0;ipart<(int)partvector.size();ipart++)
860		{
861		dout << "     - " << partvector[ipart] << endl;
862		}
#	Line 768 | Line 865 | void splitupcut(string incut, vector<str
865
866		int atleastvalue(string expression, int &morethanlessthan) // takes in an expression such as ">2" or ">=3" and returns e.g. 3 (in both examples)
867		{
868	<	int retval=0;
868	>	//  int retval=0;
869		if(expression.substr(0,1)==">"&&expression.substr(1,1)=="=") {
870		//    dout << "The expression " << expression << " is saying that we have at least " << atoi(expression.substr(2,1).c_str()) << " jets" << endl;
871		morethanlessthan=1;
#	Line 794 | Line 891 | int atleastvalue(string expression, int
891		morethanlessthan=-1;
892		return 1+atoi(expression.substr(2,1).c_str());
893		}
894	+
895	+	return -1234567;
896		}
897
898		int do_jet_cut(string incut, int *nJets) {
#	Line 829 | Line 928 | int do_jet_cut(string incut, int *nJets)
928		dout << "Is of type >" << endl;
929		return 0;
930		}
931	+	return -12345;
932		}
933
934		string interpret_cut(string incut, bool &isJetCut, int *permittednJets)
#	Line 862 | Line 962 | string interpret_cut(string incut, bool
962
963		if(incut=="id1==0") return "ee";
964		if(incut=="id1==1") return "#mu#mu";
965	<	if(incut=="abs(mll-91.2)<20") return "|m_{l^{+}l^{-}}-m_{Z}|<20";
965	>	if(incut=="abs(mll-91)<10") return "|m_{l^{+}l^{-}}-m_{Z}|<10";
966		if(incut=="pfJetGoodID[0]") return "";
967		if(incut=="pfJetGoodID[1]") return "";
968		if((int)incut.find("pfJetGoodNum")>-1) {
#	Line 894 | Line 994 | string interpret_cuts(vector<string> &cu
994		int finalJetCut=-1;
995		int permittednJets[20];
996		for(int ijet=0;ijet<20;ijet++) permittednJets[ijet]=1;
997	<	int morethanlessthan=0;//-1: less than, 0: exactly, 1: more than
998	<	for(int icut=0;icut<cutparts.size();icut++)
997	>	//    int morethanlessthan=0;//-1: less than, 0: exactly, 1: more than
998	>	for(int icut=0;icut<(int)cutparts.size();icut++)
999		{
1000		if(icut==0) nicecut<<interpret_cut(cutparts[icut],isJetCut,permittednJets);
1001		else {
#	Line 960 | Line 1060 | TH1F * makehistofromfunction(TF1 *f1,TH1
1060		histo->SetTitle(histoname.str().c_str());
1061		histo->SetName(histoname.str().c_str());
1062		int nbins=histo->GetNbinsX();
1063	<	float low=histo->GetBinLowEdge(1);
1064	<	float hi=histo->GetBinLowEdge(histo->GetNbinsX())+histo->GetBinWidth(histo->GetNbinsX());
1063	>	//  float low=histo->GetBinLowEdge(1);
1064	>	//  float hi=histo->GetBinLowEdge(histo->GetNbinsX())+histo->GetBinWidth(histo->GetNbinsX());
1065
1066		for(int i=0;i<=nbins;i++) {
1067		histo->SetBinContent(i,(f1->Integral(histo->GetBinLowEdge(i),histo->GetBinLowEdge(i)+histo->GetBinWidth(i)))/histo->GetBinWidth(i));
#	Line 1138 | Line 1238 | bool Contains(string wholestring, string
1238		else return false;
1239		}
1240
1241	+
1242		//////////////////////////////////////////////////////////////////////////////
1243		//
1244		// http://stackoverflow.com/questions/669438/how-to-get-memory-usage-at-run-time-in-c
#	Line 1150 | Line 1251 | bool Contains(string wholestring, string
1251		/* usage:
1252		double vm2, rss2;
1253		process_mem_usage(vm2, rss2);
1254	<	cout << "Memory usage: VM: " << vm << "; RSS: " << rss << endl;
1254	>	cout << "Memory usage: VM: " << vm2 << "; RSS: " << rss2 << endl;
1255		*/
1256
1257		void process_mem_usage(double& vm_usage, double& resident_set)
#	Line 1209 | Line 1310 | TGraphAsymmErrors* produce_ratio_graph(T
1310		}
1311
1312
1313	+	TLatex* WriteSelection(int njets) {
1314	+	string sel="Loose selection";
1315	+	if(njets==3) sel="Low E_{T}^{miss} selection";
1316	+	assert(njets==2||njets==3);
1317	+	TLatex *sele = new TLatex(0.97,0.135,sel.c_str());
1318	+	sele->SetNDC(true);
1319	+	sele->SetTextColor(TColor::GetColor("#848484"));
1320	+	sele->SetTextFont(42);
1321	+	sele->SetTextAlign(32);
1322	+	sele->SetTextSize(0.03);
1323	+	sele->SetTextAngle(270);
1324	+	return sele;
1325	+	}
1326	+
1327		Double_t MarcosChi2TestX(const TH1* h1, const TH1* h2, Double_t &chi2, Int_t &ndf, Int_t &igood, Option_t *option)
1328		{
1329
#	Line 1231 | Line 1346 | Double_t MarcosChi2TestX(const TH1* h1,
1346		return prob;
1347		}
1348
1349	<	void save_with_ratio(TH1F *nominator, TH1F *denominator, TVirtualPad *canvas, string savemeas, bool do_bpred_ratio=false, bool extendrange=false, string yaxistitle="ratio") {
1349	>	void Save_With_Ratio(TH1F *nominator, TH1F *denominator, TVirtualPad *orig_canvas, string savemeas, bool do_bpred_ratio=false, bool extendrange=false, string yaxistitle="ratio",TH1F *syshisto=NULL) {
1350		//this function saves the pad being passed as well as a new one including the ratio.
1351	<	CompleteSave(canvas,savemeas);
1352	<
1351	>
1352	>	orig_canvas->cd();
1353	>	orig_canvas->Update();
1354	>	CompleteSave(orig_canvas,savemeas);
1355	>	TVirtualPad *canvas = (TVirtualPad*) orig_canvas->Clone("TempCanvas");//otherwise the Ratio_main_canvas will own our pad and destroy it upon deletion
1356	>
1357		float bottommargin=gStyle->GetPadBottomMargin();
1358		float canvas_height=gStyle->GetCanvasDefH();
1359		float canvas_width=gStyle->GetCanvasDefW();
#	Line 1245 | Line 1364 | void save_with_ratio(TH1F *nominator, TH
1364
1365		float xstretchfactor=((1-ratiospace)*(1-gStyle->GetPadTopMargin()))/((1)*ratiospace);
1366
1367	<	TCanvas *main_canvas = new TCanvas("main_canvas","main_canvas",(Int_t)canvas_width,(Int_t)(canvas_height*(1+ratiospace)));
1367	>	TCanvas *Ratio_main_canvas = new TCanvas("Ratio_main_canvas","Ratio_main_canvas",(Int_t)canvas_width,(Int_t)(canvas_height*(1+ratiospace)));
1368		TPad *mainpad = new TPad("mainpad","mainpad",0,1-(1.0/(1+ratiospace)),1,1);//top (main) pad
1369		TPad *coverpad = new TPad("coverpad","coverpad",gStyle->GetPadLeftMargin()-0.008,1-(1.0/(1+ratiospace))-0.04,1,1-(1.0/(1+ratiospace))+0.103);//pad covering up the x scale
1370		TPad *bottompad = new TPad("bottompad", "Ratio Pad",0,0,1,(1-(1-bottommargin)/(1+ratiospace))-0.015); //bottom pad
1371
1372	<	main_canvas->Range(0,0,1,1);
1373	<	main_canvas->SetBorderSize(0);
1374	<	main_canvas->SetFrameFillColor(0);
1256	<
1372	>	Ratio_main_canvas->Range(0,0,1,1);
1373	>	Ratio_main_canvas->SetBorderSize(0);
1374	>	Ratio_main_canvas->SetFrameFillColor(0);
1375		mainpad->Draw();
1376		mainpad->cd();
1377		mainpad->Range(0,0,1,1);
#	Line 1263 | Line 1381 | void save_with_ratio(TH1F *nominator, TH
1381		canvas->Range(0,0,1,1);
1382		canvas->Draw("same");
1383		mainpad->Modified();
1384	<	main_canvas->cd();
1384	>	Ratio_main_canvas->cd();
1385		coverpad->Draw();
1386		coverpad->cd();
1387		coverpad->Range(0,0,1,1);
#	Line 1271 | Line 1389 | void save_with_ratio(TH1F *nominator, TH
1389		coverpad->SetBorderSize(0);
1390		coverpad->SetFrameFillColor(0);
1391		coverpad->Modified();
1392	<	main_canvas->cd();
1392	>	Ratio_main_canvas->cd();
1393		bottompad->SetTopMargin(ratiotopmargin);
1394		bottompad->SetBottomMargin(ratiobottommargin);
1395		bottompad->Draw();
#	Line 1280 | Line 1398 | void save_with_ratio(TH1F *nominator, TH
1398		bottompad->SetFillColor(kWhite);
1399		TH1F *ratio = (TH1F*)nominator->Clone(GetNumericHistoName().c_str());
1400		ratio->Divide(denominator);
1283	–
1284	–
1401		TGraphAsymmErrors *eratio;
1402	+	TH1F *SystDown;
1403	+	TH1F *SystUp;
1404	+
1405		if(!do_bpred_ratio) eratio = produce_ratio_graph(ratio);
1406		else {
1407		bool using_data=false;
1408		if((int)savemeas.find("Data")>0) using_data=true;
1409		eratio = histRatio(nominator,denominator,using_data,PlottingSetup::global_ratio_binning,false);
1410	+	if(syshisto!=0) {
1411	+	SystDown = (TH1F*)syshisto->Clone("SystDown");
1412	+	SystUp = (TH1F*)syshisto->Clone("SystUp");
1413	+	for(int i=1;i<=syshisto->GetNbinsX();i++) {
1414	+	SystDown->SetBinContent(i,1-syshisto->GetBinContent(i));
1415	+	SystUp->SetBinContent(i,1+syshisto->GetBinContent(i));
1416	+	}
1417	+	SystDown->SetLineColor(TColor::GetColor("#006DE1"));
1418	+	SystUp->SetLineColor(TColor::GetColor("#006DE1"));
1419	+	}
1420		for(int i=1;i<=ratio->GetNbinsX();i++) {
1421		ratio->SetBinContent(i,0);
1422		ratio->SetBinError(i,0);
1423		}
1424	+
1425	+	}
1426	+
1427	+	if(syshisto && !do_bpred_ratio) {
1428	+	SystDown = (TH1F*)syshisto->Clone("SystDown");
1429	+	SystUp = (TH1F*)syshisto->Clone("SystUp");
1430	+	for(int i=1;i<=syshisto->GetNbinsX();i++) {
1431	+	SystDown->SetBinContent(i,1-syshisto->GetBinContent(i));
1432	+	SystUp->SetBinContent(i,1+syshisto->GetBinContent(i));
1433	+	}
1434	+	SystDown->SetLineColor(TColor::GetColor("#006DE1"));
1435	+	SystUp->SetLineColor(TColor::GetColor("#006DE1"));
1436		}
1437		eratio->SetFillColor(TColor::GetColor("#00ADE1"));
1297	–
1438		ratio->SetTitle("");
1439	<	ratio->GetYaxis()->SetRangeUser(0.5,1.5);
1439	>	ratio->GetYaxis()->SetRangeUser(0.0,2.0);
1440		if(do_bpred_ratio) ratio->GetYaxis()->SetRangeUser(0.0,2.0);
1441		if(extendrange) ratio->GetYaxis()->SetRangeUser(0.0,4.0);
1442		ratio->GetXaxis()->SetTitle(nominator->GetXaxis()->GetTitle());
#	Line 1312 | Line 1452 | void save_with_ratio(TH1F *nominator, TH
1452		ratio->GetYaxis()->SetNdivisions(502,false);
1453		ratio->SetFillColor(TColor::GetColor("#58D3F7"));
1454		ratio->SetMarkerSize(0);
1455	<	ratio->Draw("e2");
1456	<	eratio->Draw("2");
1455	>	ratio->Draw("e20");
1456	>
1457	>	TGraphAsymmErrors *ratio_center = (TGraphAsymmErrors*)eratio->Clone("ratio_center");
1458	>	for(int i=0;i<ratio_center->GetN();i++) {
1459	>	ratio_center->SetPointError(i,ratio_center->GetErrorXlow(i),ratio_center->GetErrorXhigh(i),0.005,0.005);
1460	>	}
1461	>	ratio_center->SetFillColor(TColor::GetColor("#006381"));
1462	>
1463	>	if(syshisto!=0) {
1464	>	//    sysratio->Draw("2");
1465	>	//    eratio->Draw("2,same");
1466	>	eratio->Draw("20");
1467	>	SystDown->Draw("histo,same");
1468	>	SystUp->Draw("histo,same");
1469	>	} else {
1470	>	eratio->Draw("20");
1471	>	}
1472	>	ratio_center->Draw("20");
1473		ratio->Draw("same,axis");
1474		TLine *oneline = new TLine(ratio->GetXaxis()->GetBinLowEdge(1),1,ratio->GetXaxis()->GetBinLowEdge(ratio->GetNbinsX())+ratio->GetXaxis()->GetBinWidth(ratio->GetNbinsX()),1);
1475		oneline->SetLineStyle(2);
1476		oneline->SetLineColor(kBlue);
1477		oneline->Draw("same");
1478	<
1479	<	main_canvas->cd();
1480	<	main_canvas->Modified();
1481	<	main_canvas->cd();
1482	<	main_canvas->SetSelected(main_canvas);
1478	>	if(PlottingSetup::DrawMetSignalRegionMllLines) {
1479	>	cout << "Drawing extra lines in ratio this time around ... " << endl;
1480	>	float RatioYMax=2.0;
1481	>	if(extendrange) RatioYMax=4.0;
1482	>
1483	>	TLine *SRline   = new TLine(70,0,70,RatioYMax);
1484	>	TLine *ZLowLine = new TLine(91.2-PlottingSetup::ConsideredZWidth,0,91.2-PlottingSetup::ConsideredZWidth,RatioYMax);
1485	>	TLine *ZHiLine  = new TLine(91.2+PlottingSetup::ConsideredZWidth,0,91.2+PlottingSetup::ConsideredZWidth,RatioYMax);
1486	>
1487	>	SRline->SetLineStyle(2);
1488	>	ZLowLine->SetLineStyle(2);
1489	>	ZHiLine->SetLineStyle(2);
1490	>	SRline->SetLineColor(kGray+2);
1491	>	ZLowLine->SetLineColor(kGray+2);
1492	>	ZHiLine->SetLineColor(kGray+2);
1493	>	SRline->Draw();
1494	>	ZLowLine->Draw();
1495	>	ZHiLine->Draw();
1496	>	}
1497	>
1498	>	Ratio_main_canvas->cd();
1499	>	Ratio_main_canvas->Modified();
1500	>	Ratio_main_canvas->cd();
1501	>	Ratio_main_canvas->SetSelected(Ratio_main_canvas);
1502
1503	<	CompleteSave(main_canvas,savemeas+"_withratio");
1503	>	CompleteSave(Ratio_main_canvas,savemeas+"_withratio");
1504		bottompad->cd();
1505
1506		Double_t chi2;
1507		Int_t ndf,igood;
1508	<	Double_t res=0.0;
1508	>	//  Double_t res=0.0;
1509		Double_t chi2prob = MarcosChi2TestX(nominator,denominator,chi2,ndf,igood,"UW");
1510
1511		stringstream Chi2text;
#	Line 1345 | Line 1520 | void save_with_ratio(TH1F *nominator, TH
1520		chi2probbox->SetTextSize(0.08);
1521		chi2probbox->SetTextAlign(11);
1522		chi2probbox->Draw();
1523	<	CompleteSave(main_canvas,savemeas+"_withratio_and_Chi2");
1523	>
1524	>	stringstream CompRatio;
1525	>	CompRatio << "Ratio: " << nominator->Integral() << " / " << denominator->Integral() << " \n " << nominator->Integral()/denominator->Integral();
1526	>	TText *CompleteRatio = write_text(0.98,0.04,CompRatio.str());
1527	>	CompleteRatio->SetTextSize(0.08);
1528	>	CompleteRatio->SetTextAlign(31);
1529	>	CompleteRatio->Draw();
1530	>
1531	>	CompleteSave(Ratio_main_canvas,savemeas+"_withratio_and_Chi2");
1532
1533		//  float KS   = nominator->KolmogorovTest(denominator);
1534		//  stringstream KStext;
1535		//  Chi2text << "KS = " << KS << endl;
1536		//cout << "Found : " << KStext.str() << endl;
1537
1538	+	delete eratio;
1539	+	delete ratio_center;
1540	+	delete ratio;
1541	+	delete Ratio_main_canvas;
1542	+	}
1543
1544	<	delete main_canvas;
1544	>	void Save_With_Ratio_And_Line(TH1F *nominator, TH1F *denominator, TVirtualPad *orig_canvas, string savemeas, bool do_bpred_ratio=false, bool extendrange=false, string yaxistitle="ratio",TH1F *syshisto=NULL) {
1545	>	PlottingSetup::ConsideredZWidth=10.0;
1546	>	PlottingSetup::DrawMetSignalRegionMllLines=true;
1547	>	Save_With_Ratio(nominator, denominator, orig_canvas, savemeas, do_bpred_ratio, extendrange, yaxistitle,syshisto);
1548	>	PlottingSetup::DrawMetSignalRegionMllLines=false;
1549	>	}
1550	>
1551	>	void Save_With_Ratio_And_Line(TH1F *nominator, THStack denominator, TVirtualPad *canvas, string savemeas, bool do_bpred_ratio=false) {
1552	>	TH1F *denominator_histo = (TH1F*) CollapseStack(denominator,"TemporaryStack");
1553	>	Save_With_Ratio_And_Line(nominator, denominator_histo, canvas, savemeas, do_bpred_ratio);
1554	>	delete denominator_histo;
1555		}
1556
1557	+	void save_with_ratio_and_sys_band(float ConsideredZWidth, TH1F *nominator, TH1F *denominator, TVirtualPad *orig_canvas, string savemeas, bool do_bpred_ratio=false, bool extendrange=false, string yaxistitle="ratio",TH1F *syshisto=NULL) {
1558	+	//this function saves the pad being passed as well as a new one including the SysRatio.
1559	+	orig_canvas->cd();
1560	+	orig_canvas->Update();
1561	+	CompleteSave(orig_canvas,savemeas);
1562	+
1563	+	TVirtualPad *canvas = (TVirtualPad*) orig_canvas->Clone("TempCanvas");//otherwise the Ratio_main_canvas will own our pad and destroy it upon deletion
1564
1565	<	TH1F* CollapseStack(THStack stack) {
1565	>	float bottommargin=gStyle->GetPadBottomMargin();
1566	>	float canvas_height=gStyle->GetCanvasDefH();
1567	>	float canvas_width=gStyle->GetCanvasDefW();
1568	>	float ratiospace=0.25;// space the ratio should take up (relative to original pad)
1569	>
1570	>	float ratiobottommargin=0.3;
1571	>	float ratiotopmargin=0.1;
1572	>
1573	>	float xstretchfactor=((1-ratiospace)*(1-gStyle->GetPadTopMargin()))/((1)*ratiospace);
1574	>
1575	>	TCanvas *Ratio_main_canvas = new TCanvas("Ratio_main_canvas","Ratio_main_canvas",(Int_t)canvas_width,(Int_t)(canvas_height*(1+ratiospace)));
1576	>	TPad *mainpad = new TPad("mainpad","mainpad",0,1-(1.0/(1+ratiospace)),1,1);//top (main) pad
1577	>	TPad *coverpad = new TPad("coverpad","coverpad",gStyle->GetPadLeftMargin()-0.008,1-(1.0/(1+ratiospace))-0.04,1,1-(1.0/(1+ratiospace))+0.103);//pad covering up the x scale
1578	>	TPad *bottompad = new TPad("bottompad", "Ratio Pad",0,0,1,(1-(1-bottommargin)/(1+ratiospace))-0.015); //bottom pad
1579	>
1580	>	Ratio_main_canvas->Range(0,0,1,1);
1581	>	Ratio_main_canvas->SetBorderSize(0);
1582	>	Ratio_main_canvas->SetFrameFillColor(0);
1583	>
1584	>	mainpad->Draw();
1585	>	mainpad->cd();
1586	>	mainpad->Range(0,0,1,1);
1587	>	mainpad->SetFillColor(kWhite);
1588	>	mainpad->SetBorderSize(0);
1589	>	mainpad->SetFrameFillColor(0);
1590	>	canvas->Range(0,0,1,1);
1591	>	canvas->Draw("same");
1592	>	mainpad->Modified();
1593	>	Ratio_main_canvas->cd();
1594	>	coverpad->Draw();
1595	>	coverpad->cd();
1596	>	coverpad->Range(0,0,1,1);
1597	>	coverpad->SetFillColor(kWhite);
1598	>	coverpad->SetBorderSize(0);
1599	>	coverpad->SetFrameFillColor(0);
1600	>	coverpad->Modified();
1601	>	Ratio_main_canvas->cd();
1602	>	bottompad->SetTopMargin(ratiotopmargin);
1603	>	bottompad->SetBottomMargin(ratiobottommargin);
1604	>	bottompad->Draw();
1605	>	bottompad->cd();
1606	>	bottompad->Range(0,0,1,1);
1607	>	bottompad->SetFillColor(kWhite);
1608	>	TH1F *ratio = (TH1F*)nominator->Clone(GetNumericHistoName().c_str());
1609	>	ratio->Divide(denominator);
1610	>
1611	>
1612	>	TGraphAsymmErrors *eratio;
1613	>	TGraphAsymmErrors *SysEnvelope = new TGraphAsymmErrors();
1614	>
1615	>	if(syshisto && !do_bpred_ratio) {
1616	>	for(int i=1;i<=syshisto->GetNbinsX();i++) {
1617	>	float dx=syshisto->GetBinWidth(i)/2.0;
1618	>	float dy=syshisto->GetBinContent(i);
1619	>	if(dy<0.01) dy=0.07;
1620	>	SysEnvelope->SetPoint(i-1,syshisto->GetBinCenter(i),1.0);
1621	>	SysEnvelope->SetPointError(i-1,dx,dx,dy,dy);
1622	>	}
1623	>	SysEnvelope->SetFillColor(TColor::GetColor("#006DE1"));
1624	>	}
1625	>
1626	>	ratio->SetTitle("");
1627	>	ratio->GetYaxis()->SetRangeUser(0.0,2.0);
1628	>	if(do_bpred_ratio) ratio->GetYaxis()->SetRangeUser(0.0,2.0);
1629	>	if(extendrange) ratio->GetYaxis()->SetRangeUser(0.0,4.0);
1630	>	ratio->GetXaxis()->SetTitle(nominator->GetXaxis()->GetTitle());
1631	>	ratio->GetXaxis()->CenterTitle();
1632	>	ratio->GetYaxis()->SetTitle(yaxistitle.c_str());
1633	>	ratio->GetYaxis()->SetTitleOffset(0.4);
1634	>	ratio->GetYaxis()->CenterTitle();
1635	>	ratio->SetStats(0);
1636	>	ratio->GetXaxis()->SetLabelSize(xstretchfactor*ratio->GetXaxis()->GetLabelSize());
1637	>	ratio->GetYaxis()->SetLabelSize(xstretchfactor*ratio->GetYaxis()->GetLabelSize());
1638	>	ratio->GetXaxis()->SetTitleSize(xstretchfactor*gStyle->GetTitleSize());
1639	>	ratio->GetYaxis()->SetTitleSize(xstretchfactor*gStyle->GetTitleSize());
1640	>	ratio->GetYaxis()->SetNdivisions(502,false);
1641	>	ratio->SetFillColor(TColor::GetColor("#58D3F7"));
1642	>	ratio->Draw("e0");
1643	>
1644	>	if(syshisto!=0) {
1645	>	SysEnvelope->SetFillColor(TColor::GetColor("#FE9A2E"));
1646	>	SysEnvelope->Draw("2,0,same");
1647	>	ratio->Draw("e0,same");
1648	>	} else {
1649	>	eratio->Draw("0");
1650	>	}
1651	>	ratio->Draw("same,axis");
1652	>	TLine *oneline = new TLine(ratio->GetXaxis()->GetBinLowEdge(1),1,ratio->GetXaxis()->GetBinLowEdge(ratio->GetNbinsX())+ratio->GetXaxis()->GetBinWidth(ratio->GetNbinsX()),1);
1653	>	oneline->SetLineStyle(2);
1654	>	oneline->SetLineColor(kBlue);
1655	>	oneline->Draw("same");
1656	>
1657	>
1658	>
1659	>	if(ConsideredZWidth>0.1) {
1660	>	cout << "Drawing extra lines in ratio this time around ... " << endl;
1661	>	float RatioYMax=2.0;
1662	>	if(extendrange) RatioYMax=4.0;
1663	>
1664	>	TLine *SRline   = new TLine(70,0,70,RatioYMax);
1665	>	TLine *ZLowLine = new TLine(91.2-ConsideredZWidth,0,91.2-ConsideredZWidth,RatioYMax);
1666	>	TLine *ZHiLine  = new TLine(91.2+ConsideredZWidth,0,91.2+ConsideredZWidth,RatioYMax);
1667	>
1668	>	SRline->SetLineStyle(2);
1669	>	ZLowLine->SetLineStyle(2);
1670	>	ZHiLine->SetLineStyle(2);
1671	>	SRline->SetLineColor(kGray+2);
1672	>	ZLowLine->SetLineColor(kGray+2);
1673	>	ZHiLine->SetLineColor(kGray+2);
1674	>	SRline->Draw();
1675	>	ZLowLine->Draw();
1676	>	ZHiLine->Draw();
1677	>	}
1678	>
1679	>	Ratio_main_canvas->cd();
1680	>	Ratio_main_canvas->Modified();
1681	>	Ratio_main_canvas->cd();
1682	>	Ratio_main_canvas->SetSelected(Ratio_main_canvas);
1683	>
1684	>	CompleteSave(Ratio_main_canvas,savemeas+"_withSysRatio");
1685	>	bottompad->cd();
1686	>
1687	>	Double_t chi2;
1688	>	Int_t ndf,igood;
1689	>	Double_t chi2prob = MarcosChi2TestX(nominator,denominator,chi2,ndf,igood,"UW");
1690	>
1691	>	stringstream Chi2text;
1692	>	Chi2text << "#chi^{2} / ndf: " << chi2 << " / " << ndf;
1693	>	stringstream Chi2probtext;
1694	>	Chi2probtext << "#chi^{2} prob: " << chi2prob;
1695	>	TText* chi2box = write_text(0.02,0.11,Chi2text.str());
1696	>	chi2box->SetTextSize(0.08);
1697	>	chi2box->SetTextAlign(11);
1698	>	chi2box->Draw();
1699	>	TText* chi2probbox = write_text(0.02,0.04,Chi2probtext.str());
1700	>	chi2probbox->SetTextSize(0.08);
1701	>	chi2probbox->SetTextAlign(11);
1702	>	chi2probbox->Draw();
1703	>
1704	>	stringstream CompRatio;
1705	>	CompRatio << "Ratio: " << nominator->Integral() << " / " << denominator->Integral() << " \n " << nominator->Integral()/denominator->Integral();
1706	>	TText *CompleteRatio = write_text(0.98,0.04,CompRatio.str());
1707	>	CompleteRatio->SetTextSize(0.08);
1708	>	CompleteRatio->SetTextAlign(31);
1709	>	CompleteRatio->Draw();
1710	>
1711	>	CompleteSave(Ratio_main_canvas,savemeas+"_withSysRatio_and_Chi2");
1712	>	delete Ratio_main_canvas;
1713	>	delete ratio;
1714	>	}
1715	>
1716	>	TH1F* CollapseStack(THStack stack,TString hname="CollapsedStack") {
1717		TH1F *bhist = ((TH1F*)((stack.GetHists())->At(0)));
1718	<	TH1F *basehisto = (TH1F*)bhist->Clone("base");
1718	>	TH1F *basehisto = (TH1F*)bhist->Clone(hname);
1719		TIter next(stack.GetHists());
1720		TH1F *h;
1721		int counter=0;
#	Line 1371 | Line 1727 | TH1F* CollapseStack(THStack stack) {
1727		return basehisto;
1728		}
1729
1730	<	void save_with_ratio(TH1F *nominator, THStack denominator, TVirtualPad *canvas, string savemeas, bool do_bpred_ratio=false) {
1731	<	save_with_ratio(nominator, CollapseStack(denominator), canvas, savemeas, do_bpred_ratio);
1730	>	void Save_With_Ratio(TH1F *nominator, THStack denominator, TVirtualPad *canvas, string savemeas, bool do_bpred_ratio=false) {
1731	>	TH1F *denominator_histo = (TH1F*) CollapseStack(denominator);
1732	>	Save_With_Ratio(nominator, denominator_histo, canvas, savemeas, do_bpred_ratio);
1733	>	delete denominator_histo;
1734		}
1735
1736		void flag_this_change(string function, int line, int checked=0) {
#	Line 1388 | Line 1746 | void flag_this_change(string function, i
1746		if(checked==2) write_info(function,peakmodificationwarning.str());
1747		}
1748
1749	<	void write_analysis_type(bool isonpeak) {
1749	>	void write_analysis_type(bool isonpeak, bool dobtag) {
1750		//http://www.network-science.de/ascii/, ogre
1751		if(!PlottingSetup::publicmode) {
1752		if(isonpeak) {
1753	+	//font: big
1754		dout << "\033[1;34m" << endl;
1755	<	dout <<  "    //\\\\                                  _                          _           _     " << endl;
1756	<	dout <<  "   //  \\\\     ___  _ __  _ __   ___  __ _| | __     __ _ _ __   __ _| |_   _ ___(_)___" << endl;
1757	<	dout <<  "  //    \\\\   / _ \\| '_ \\| '_ \\ / _ \\/ _` | |/ /    / _` | '_ \\ / _` | | | | / __| / __|" << endl;
1758	<	dout <<  " //      \\\\ | (_) | | | | |_) |  __/ (_| |   <    | (_| | | | | (_| | | |_| \\__ \\ \\__ \\" << endl;
1759	<	dout <<  "//        \\\\ \\___/|_| |_| .__/ \\___|\\__,_|_|\\_\\    \\__,_|_| |_|\\__,_|_|\\__, |___/_|___/" << endl;
1760	<	dout <<  "                        |_|                                            |___/" << endl;
1755	>	dout <<  "                             _ __________                                                         " << endl;
1756	>	dout <<  "                            | |___  /  _ \\                                                       " << endl;
1757	>	dout <<  "       /\\           ___     | |  / /| |_) |                                                       " << endl;
1758	>	dout <<  "      /  \\         /  _|_   | | / / |  _ <                                                         " << endl;
1759	>	dout <<  "     /    \\        | (_| |__| |/ /__| |_) |                                                       " << endl;
1760	>	dout <<  "                   \\___|\\____//_____|____/                                                       " << endl;
1761		} else {
1762	+	//font: big
1763		dout << "\033[1;33m" << endl;
1764	<	dout <<  "                    __  __                  _                          _           _     " << endl;
1765	<	dout <<  "              ___  / _|/ _|_ __   ___  __ _| | __     __ _ _ __   __ _| |_   _ ___(_)___ " << endl;
1766	<	dout <<  "    /\\       / _ \\| |_| |_| '_ \\ / _ \\/ _` | |/ /    / _` | '_ \\ / _` | | | | / __| / __|" << endl;
1767	<	dout <<  " __/  \\__   | (_) |  _|  _| |_) |  __/ (_| |   <    | (_| | | | | (_| | | |_| \\__ \\ \\__ \\" << endl;
1768	<	dout <<  "             \\___/|_| |_| | .__/ \\___|\\__,_|_|\\_\\    \\__,_|_| |_|\\__,_|_|\\__, |___/_|___/" << endl;
1769	<	dout <<  "                          |_|                                            |___/           " << endl;
1764	>	dout <<  "                    _      _ __________                                                             " << endl;
1765	>	dout <<  "                   (_)    | |___  /  _ \\                                                            " << endl;
1766	>	dout <<  "       /\\           _     | |  / /| |_) |                                                           " << endl;
1767	>	dout <<  "    /\\/  \\         | |_   | | / / |  _ <                                                            " << endl;
1768	>	dout <<  "   /      \\        | | |__| |/ /__| |_) |                                                           " << endl;
1769	>	dout <<  "                   |_|\\____//_____|____/                                                            " << endl;
1770		}
1771	<	dout << "\033[0m" << endl;
1771	>
1772	>	if(dobtag) {
1773	>	//font: big
1774	>	dout <<  "\033[1;32m                  \\     /                                                                " << endl;
1775	>	dout <<  "             \\    o ^ o    /                                                             " << endl;
1776	>	dout <<  "               \\ (     ) /                                                               " << endl;
1777	>	dout <<  "    ____________(%%%%%%%)____________        _          _ __________                     " << endl;
1778	>	dout <<  "   (     /   /  )%%%%%%%(  \\   \\     )      | |        | |___  /  _ \\                    " << endl;
1779	>	dout <<  "   (___/___/__/           \\__\\___\\___)      | |__      | |  / /| |_) |                   " << endl;
1780	>	dout <<  "      (     /  /(%%%%%%%)\\  \\     )         | '_ \\ _   | | / / |  _ <                    " << endl;
1781	>	dout <<  "       (__/___/ (%%%%%%%) \\___\\__)          | |_) | |__| |/ /__| |_) |                   " << endl;
1782	>	dout <<  "               /(       )\\                  |_.__/ \\____//_____|____/                    " << endl;
1783	>	dout <<  "             /   (%%%%%)   \\                                                             " << endl;
1784	>	dout <<  "                  (%%%)                                                                  " << endl;
1785	>	dout <<  "                    !                                                                    " << endl;
1786	>	}
1787	>	} else {
1788	>	//we're in public! don't advertise what we're up to :-)
1789	>	dout << "~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~" << endl;
1790	>	dout << "                              Starting the analysis                     " << endl;
1791	>	dout << "                                     i:" << !isonpeak << " , b " << dobtag << endl;
1792	>	dout << "~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~" << endl;
1793		}
1794	+	dout << "\033[0m" << endl;
1795	+
1796		}
1797
1798
#	Line 1418 | Line 1801 | vector<string> StringSplit(string str, s
1801		vector<string> results;
1802
1803		int cutAt;
1804	<	while( (cutAt = str.find_first_of(delim)) != str.npos )
1804	>	while( (cutAt = str.find_first_of(delim)) != (int)str.npos )
1805		{
1806		if(cutAt > 0)
1807		{
#	Line 1435 | Line 1818 | vector<string> StringSplit(string str, s
1818
1819		void manually_set_jzb_cuts(vector<float> &jzb_cut,string jzbcut_string) {
1820		vector<string> jzbcutvector = StringSplit(jzbcut_string,",");
1821	<	for(int i=0;i<jzbcutvector.size();i++) {
1821	>	for(int i=0;i<(int)jzbcutvector.size();i++) {
1822		jzb_cut.push_back(atoi(jzbcutvector[i].c_str()));
1823		dout << "Added a JZB cut manually at " << atoi(jzbcutvector[i].c_str()) << endl;
1824		}
#	Line 1449 | Line 1832 | void process_directory(TString directory
1832		dp = opendir (directory);
1833		if (dp != NULL)
1834		{
1835	<	while (ep = readdir (dp))
1835	>	while ((ep = readdir (dp)))
1836		{
1837		string filename=(string)ep->d_name;
1838	<	if(filename.find(".root")!=-1)
1838	>	if((int)filename.find(".root")!=-1)
1839		{
1840		files.push_back(string(directory)+filename);
1841		}
#	Line 1463 | Line 1846 | void process_directory(TString directory
1846		perror ("Couldn't open the directory");
1847		}
1848
1849	+	void ClearHisto(TH1* histo) {
1850	+	for(int ix=0;ix<=histo->GetNbinsX()+1;ix++) {//also reset over- and underflow
1851	+	histo->SetBinContent(ix,0);
1852	+	histo->SetBinContent(ix,0);
1853	+	}
1854	+	}
1855	+
1856	+	void AddSquared(TH1F *histo, TH1F *addthis, float weight) {
1857	+	for(int ix=1;ix<=addthis->GetNbinsX()+1;ix++) {//also include over- and underflow
1858	+	float bincenter=addthis->GetBinCenter(ix);
1859	+	int nBinHisto= histo->FindBin(bincenter);
1860	+	float old_content=histo->GetBinContent(nBinHisto);
1861	+	histo->SetBinContent(nBinHisto,old_content+((addthis->GetBinContent(ix))*(addthis->GetBinContent(ix))*weight));
1862	+	}
1863	+	}
1864	+
1865	+	void SQRT(TH1* histo) {
1866	+	for(int ix=0;ix<=histo->GetNbinsX()+1;ix++) {//also include over- and underflow
1867	+	histo->SetBinContent(ix,TMath::Sqrt(histo->GetBinContent(ix)));
1868	+	}
1869	+	}
1870	+
1871		/*string GetCompleteHostname() {
1872		struct addrinfo hints, *info, *p;
1873		int gai_result;
#	Line 1494 | Line 1899 | const char* concatenate (const char* a,
1899		return bla.str().c_str();
1900		}
1901
1902	+
1903	+	string ReplaceCharacter(string originalstring,string replacethis,string replacewiththis)
1904	+	{
1905	+	int pos = originalstring.find(replacethis);
1906	+	while(pos != -1) {
1907	+	originalstring.replace(pos, replacewiththis.length(), replacewiththis);
1908	+	pos = originalstring.find(replacethis);
1909	+	}
1910	+	return originalstring;
1911	+	}
1912	+
1913	+	string removefunnystring(string name) {
1914	+	name=ReplaceCharacter(name,"[","_");
1915	+	name=ReplaceCharacter(name,"]","_");
1916	+	name=ReplaceCharacter(name,"{","_");
1917	+	name=ReplaceCharacter(name,"}","_");
1918	+	name=ReplaceCharacter(name,"(","_");
1919	+	name=ReplaceCharacter(name,")","_");
1920	+	name=ReplaceCharacter(name," ","_");
1921	+	name=ReplaceCharacter(name,".","_");
1922	+	name=ReplaceCharacter(name,",","_");
1923	+	name=ReplaceCharacter(name,";","_");
1924	+	name=ReplaceCharacter(name,":","_");
1925	+	name=ReplaceCharacter(name,"'","_");
1926	+	name=ReplaceCharacter(name,"$","_");
1927	+	name=ReplaceCharacter(name,"@","_");
1928	+	name=ReplaceCharacter(name,"#","_");
1929	+	name=ReplaceCharacter(name,"+","_");
1930	+	name=ReplaceCharacter(name,"-","_");
1931	+	name=ReplaceCharacter(name,"/","_");
1932	+	return name;
1933	+	}
1934	+
1935	+	TPolyLine* GetFitUncertaintyShape(TF1 *fit, float low, float high,float x0,float x1) {
1936	+	int nPoints=1000;
1937	+	double x[nPoints+1];
1938	+	double y[nPoints];
1939	+
1940	+	float par1=fit->GetParameter(0);
1941	+	float dpar1=fit->GetParError(0);
1942	+	float par2=fit->GetParameter(1);
1943	+	float dpar2=fit->GetParError(1);
1944	+
1945	+	float step=(x1-x0)/(nPoints/2.0 -1);
1946	+	for(int i=0;i<nPoints/2.0;i++) {
1947	+	x[i]=x0+i*step;
1948	+	y[i]=(par1+x[i]*par2)    -    sqrt(dpar1*dpar1+dpar2*dpar2*x[i]*x[i]);
1949	+	x[nPoints-1-i]=x[i];
1950	+	y[nPoints-1-i]=(par1+x[i]*par2)    +    sqrt(dpar1*dpar1+dpar2*dpar2*x[i]*x[i]);
1951	+	if(y[i]<low) y[i]=low;
1952	+	if(y[i]>high) y[i]=high;
1953	+	if(y[nPoints-1-i]<low) y[nPoints-1-i]=low;
1954	+	if(y[nPoints-1-i]>high) y[nPoints-1-i]=high;
1955	+	}
1956	+
1957	+	x[nPoints]=x[0];
1958	+	y[nPoints]=y[0];
1959	+
1960	+	TPolyLine *l = new TPolyLine(nPoints+1,x,y);
1961	+	l->SetFillColor(TColor::GetColor("#A2A2FA"));
1962	+	l->SetLineColor(TColor::GetColor("#A2A2FA"));
1963	+	l->SetLineWidth(1);
1964	+	return l;
1965	+	}
1966	+
1967	+	//code for ReplaceAll copied from
1968	+	//http://stackoverflow.com/questions/5343190/c-how-do-i-replace-all-instances-of-of-a-string-with-another-string
1969	+	string ReplaceAll(string str, string from, string to) {
1970	+	size_t start_pos = 0;
1971	+	while((start_pos = str.find(from, start_pos)) != std::string::npos) {
1972	+	str.replace(start_pos, from.length(), to);
1973	+	start_pos += to.length(); // ...
1974	+	}
1975	+	return str;
1976	+	}
1977	+
1978	+	TPolyLine* GetFitUncertaintyShape(TH1 *histo, string fitname, float low, float high) {
1979	+	TF1 *fit = (TF1*)histo->GetFunction(fitname.c_str());
1980	+	if(!fit) {
1981	+	cout << "CANNOT PRODUCE FIT UNCERTAINTY AS THERE IS NO FIT CALLED " << fitname << " IN " << histo->GetName() << endl;
1982	+	return NULL;
1983	+	}
1984	+	float x0=histo->GetBinLowEdge(1);
1985	+	float x1=histo->GetBinLowEdge(histo->GetNbinsX())+histo->GetBinWidth(histo->GetNbinsX());
1986	+	return GetFitUncertaintyShape(fit, low, high,x0,x1);
1987	+	}
1988	+
1989	+	TPolyLine* GetFitUncertaintyShape(TGraphErrors *gr, string fitname, float low, float high, float minx=-999, float maxx=-999) {
1990	+	TF1 *fit = (TF1*)gr->GetFunction(fitname.c_str());
1991	+	if(!fit) {
1992	+	cout << "CANNOT PRODUCE FIT UNCERTAINTY AS THERE IS NO FIT CALLED " << fitname << " IN " << gr->GetName() << endl;
1993	+	return NULL;
1994	+	}
1995	+	double x,y;
1996	+	gr->GetPoint(0,x,y);
1997	+	float min=x;
1998	+	gr->GetPoint(gr->GetN()-1,x,y);
1999	+	float max=x;
2000	+	if(minx>-999 && maxx>-999 && minx<maxx) {
2001	+	max=maxx;
2002	+	min=minx;
2003	+	}
2004	+	return GetFitUncertaintyShape(fit, low, high,min,max);
2005	+	}
2006	+
2007		stringstream all_bugs;
2008
2009		void bug_tracker(string function, int line, string description) {
#	Line 1531 | Line 2041 | void bug_tracker(string function, int li
2041		//TODO: Write a bug summary at the end.
2042
2043		float pSTAR(float mglu, float mlsp, float mchi2) {
2044	<	float mz=91.2;
2044	>	float mz=91.0;
2045		float res=((mchi2*mchi2)-(mlsp+mz)*(mlsp+mz));
2046		res*=((mchi2*mchi2)-(mlsp-mz)*(mlsp-mz));
2047		res=TMath::Sqrt(res)/(2*mchi2);
#	Line 1539 | Line 2049 | float pSTAR(float mglu, float mlsp, floa
2049		}
2050
2051		float generalizedpSTAR(float massmother, float massdaughter1, float massdaughter2) {
1542	–	float mz=91.2;
2052		float res=((massmother*massmother)-(massdaughter1+massdaughter2)*(massdaughter1+massdaughter2));
2053		res*=((massmother*massmother)-(massdaughter1-massdaughter2)*(massdaughter1-massdaughter2));
2054		res=TMath::Sqrt(res)/(2*massmother);
2055		return res;
2056		}
2057	+
2058	+	void FindMinMax(TGraphErrors *gr, float &min, float &max) {
2059	+	double x,y;
2060	+	gr->GetPoint(0,x,y);
2061	+	min=500*y;
2062	+	max=0.0002*y;
2063	+	for(int i=0;i<gr->GetN();i++) {
2064	+	gr->GetPoint(i,x,y);
2065	+	cout << "Point at " << x << "/" << y << " with error " << gr->GetErrorY(i) << endl;
2066	+	float potmin=y-1.5*gr->GetErrorY(i); // taking 1.5 times the error to allow for enough space to present the plot
2067	+	float potmax=y+1.5*gr->GetErrorY(i); // taking 1.5 times the error to allow for enough space to present the plot
2068	+	if(potmin<min) min=potmin;
2069	+	if(potmax>max) max=potmax;
2070	+	}
2071	+	}
2072	+
2073	+
2074	+	void CleanLegends() {
2075	+	for(int ihist=(int)PlottingSetup::FakeHistoHeap.size()-1;ihist>=0;ihist--) {
2076	+	PlottingSetup::FakeHistoHeap[ihist]->Delete();
2077	+	PlottingSetup::FakeHistoHeap.pop_back();
2078	+	}
2079	+	}
2080	+
2081	+	string DigitsAfterComma(float number, int digits) {
2082	+	float temp=number*pow(10.0,digits);
2083	+	temp=int(temp)/pow(10.0,digits);
2084	+	return any2string(temp);
2085	+	}

Diff Legend

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