Plotting/Modules/ResultLibraryClass.C

#include <iostream>
#include "YieldClass.C"
#define ResultLibraryClassLoaded

using namespace std;

string ampersand=" & \t";

class tripple
{
public:
  string name;
  float first;
  float second;
};

class LeptonResult 
{
public:
  float eemm;
  float ee;
  float mm;
  float em;
};

std::ostream &operator<<(std::ostream &ostr, LeptonResult l)
{
  return ostr << l.eemm << " (" << l.ee << " , " << l.mm << ")";
}

class LibraryForOneCut {
public:
  Yield Zbkg; // Filled in ResultModule
  Yield Flavorsym; // Filled in ResultModule
  Yield total; // Filled in ResultModule
  LeptonResult observed; // Filled in ResultModule
  vector<tripple> SignalYield;
  float UpperLimit;//Filled in ResultModule
  float jzbcut; // Filled in ResultModule
  void Print();
};

void LibraryForOneCut::Print() {
  cout << "Library for jzb cut at " << this->jzbcut << " : " << endl;
  cout << "Z prediction : " << Zbkg.value << " +/- " << Zbkg.syserror << " (sys) +/- " << Zbkg.staterror << " (stat)" << endl;
  cout << "Flavorsym prediction : " << Flavorsym.value << " +/- " << Flavorsym.syserror << " (sys) +/- " << Flavorsym.staterror << " (stat)" << endl;
  cout << "Total prediction: " << total.value << " +/- " << total.syserror << " (sys) +/- " << total.staterror << " (stat)" << endl;
  cout << "Observed : " << observed.eemm << " (" << observed.ee << " , " << observed.mm << ")" << endl;
  cout << "Upper Limits: " << UpperLimit << endl;
  
}

class ResultLibrary {
public:
  vector<LibraryForOneCut> predictions;
  void Add(LibraryForOneCut);
  int Find(float);
  void Print();
};

void ResultLibrary::Add(LibraryForOneCut alpha) {
  this->predictions.push_back(alpha);
}

int ResultLibrary::Find(float jzbcut) {
  for(int i=0;i<this->predictions.size();i++) if(this->predictions[i].jzbcut==jzbcut) return i;
  // if cut wasn't found, create this library and return the new index!
  LibraryForOneCut additional;
  additional.jzbcut=jzbcut;
  this->Add(additional);
  cout << "A library for cut " << jzbcut << " was added " << endl;
  return (this->predictions).size()-1;
}

void ResultLibrary::Print() {
  cout << " \\hline" << endl;
  cout << " " << ampersand << " ";
  for(int i=0;i<this->predictions.size();i++) cout << ampersand << " JZB > " << predictions[i].jzbcut << " GeV " << " ";
  cout << " \\\\ \\hline" << endl;
  cout << " " << ampersand << " Z bkgd ";
  for(int i=0;i<this->predictions.size();i++) cout << ampersand << " " << predictions[i].Zbkg;
  cout << " \\\\" << endl;
  cout << " " << ampersand << " flavor-symmetric ";
  for(int i=0;i<this->predictions.size();i++) cout << ampersand << " " << predictions[i].Flavorsym;
  cout << " \\\\" << endl;
  cout << " " << ampersand << " Total ";
  for(int i=0;i<this->predictions.size();i++) cout << ampersand << " " << predictions[i].total;
  cout << " \\\\ \\hline" << endl;
  cout << " " << ampersand << " Observed ";
  for(int i=0;i<this->predictions.size();i++) cout << ampersand << " " << predictions[i].observed;
  cout << " \\\\ \\hline" << endl;
  cout << " " << ampersand << " Upper Limit";
  for(int i=0;i<this->predictions.size();i++) cout << ampersand << " " << predictions[i].UpperLimit;
  cout << " \\\\ \\hline";
  for(int isignal=0;isignal<this->predictions[0].SignalYield.size();isignal++) {
    cout << endl;
    cout << " " << ampersand << " " << predictions[0].SignalYield[isignal].name;
    for(int i=0;i<this->predictions.size();i++) {
      if(predictions[i].SignalYield.size()<=isignal) cout << "X";
      else cout << ampersand << " " << predictions[i].SignalYield[isignal].first << " +/- " << predictions[i].SignalYield[isignal].second << "  \\\\";
    }
  }
  cout << "  \\hline " << endl;
  
}

/*
int main() {
  LibraryForOneCut tester;
  LeptonResult observed;
  observed.ee=2;
  observed.mm=3;
  observed.eemm=5;
  Yield Zbkg;Zbkg.value=3;Zbkg.staterror=2;Zbkg.syserror=1;
  Yield Flavorsym;Flavorsym.value=6;Flavorsym.staterror=5;Flavorsym.syserror=4;
  Yield total;total.value=9;total.staterror=8;total.syserror=7;
  float UpperLimit=15;
  tester.observed=observed;
  tester.jzbcut=11;
  tester.Zbkg=Zbkg;
  tester.Flavorsym=Flavorsym;
  tester.total=total;
  tester.UpperLimit=UpperLimit;
  ResultLibrary hello;
  vector<tripple> signals;
  tripple signal;
  signal.name="LM4";
  signal.first=2;
  signal.second=3;
  (tester.SignalYield).push_back(signal);
  signal.name="LM8";
  signal.first=4;
  signal.second=5;
  tester.SignalYield.push_back(signal);
  
  
  hello.Add(tester);
  tester.SignalYield.push_back(signal);
  
  tester.observed.ee=12;
  tester.observed.mm=13;
  tester.observed.eemm=25;
  tester.jzbcut=50;
  hello.Add(tester);
  tester.jzbcut=150;
  hello.Add(tester);
  tester.Print();
  hello.Print(); 
}
*/
Revision:	1.2
Committed:	Mon Dec 5 14:26:03 2011 UTC (13 years, 5 months ago) by buchmann
Content type:	text/plain
Branch:	MAIN
CVS Tags:	cbaf_4_98ifb_paper, beforeFR20120418, cbaf_4p7ifb, HEAD
Changes since 1.1:	+1 -1 lines
Log Message:	Updated range for JZB plots
#	User	Rev	Content
1	buchmann	1.1	#include <iostream>
2			#include "YieldClass.C"
3			#define ResultLibraryClassLoaded
4
5			using namespace std;
6
7	buchmann	1.2	string ampersand=" & \t";
8	buchmann	1.1
9			class tripple
10			{
11			public:
12			string name;
13			float first;
14			float second;
15			};
16
17			class LeptonResult
18			{
19			public:
20			float eemm;
21			float ee;
22			float mm;
23			float em;
24			};
25
26			std::ostream &operator<<(std::ostream &ostr, LeptonResult l)
27			{
28			return ostr << l.eemm << " (" << l.ee << " , " << l.mm << ")";
29			}
30
31			class LibraryForOneCut {
32			public:
33			Yield Zbkg; // Filled in ResultModule
34			Yield Flavorsym; // Filled in ResultModule
35			Yield total; // Filled in ResultModule
36			LeptonResult observed; // Filled in ResultModule
37			vector<tripple> SignalYield;
38			float UpperLimit;//Filled in ResultModule
39			float jzbcut; // Filled in ResultModule
40			void Print();
41			};
42
43			void LibraryForOneCut::Print() {
44			cout << "Library for jzb cut at " << this->jzbcut << " : " << endl;
45			cout << "Z prediction : " << Zbkg.value << " +/- " << Zbkg.syserror << " (sys) +/- " << Zbkg.staterror << " (stat)" << endl;
46			cout << "Flavorsym prediction : " << Flavorsym.value << " +/- " << Flavorsym.syserror << " (sys) +/- " << Flavorsym.staterror << " (stat)" << endl;
47			cout << "Total prediction: " << total.value << " +/- " << total.syserror << " (sys) +/- " << total.staterror << " (stat)" << endl;
48			cout << "Observed : " << observed.eemm << " (" << observed.ee << " , " << observed.mm << ")" << endl;
49			cout << "Upper Limits: " << UpperLimit << endl;
50
51			}
52
53			class ResultLibrary {
54			public:
55			vector<LibraryForOneCut> predictions;
56			void Add(LibraryForOneCut);
57			int Find(float);
58			void Print();
59			};
60
61			void ResultLibrary::Add(LibraryForOneCut alpha) {
62			this->predictions.push_back(alpha);
63			}
64
65			int ResultLibrary::Find(float jzbcut) {
66			for(int i=0;i<this->predictions.size();i++) if(this->predictions[i].jzbcut==jzbcut) return i;
67			// if cut wasn't found, create this library and return the new index!
68			LibraryForOneCut additional;
69			additional.jzbcut=jzbcut;
70			this->Add(additional);
71			cout << "A library for cut " << jzbcut << " was added " << endl;
72			return (this->predictions).size()-1;
73			}
74
75			void ResultLibrary::Print() {
76			cout << " \\hline" << endl;
77			cout << " " << ampersand << " ";
78			for(int i=0;i<this->predictions.size();i++) cout << ampersand << " JZB > " << predictions[i].jzbcut << " GeV " << " ";
79			cout << " \\\\ \\hline" << endl;
80			cout << " " << ampersand << " Z bkgd ";
81			for(int i=0;i<this->predictions.size();i++) cout << ampersand << " " << predictions[i].Zbkg;
82			cout << " \\\\" << endl;
83			cout << " " << ampersand << " flavor-symmetric ";
84			for(int i=0;i<this->predictions.size();i++) cout << ampersand << " " << predictions[i].Flavorsym;
85			cout << " \\\\" << endl;
86			cout << " " << ampersand << " Total ";
87			for(int i=0;i<this->predictions.size();i++) cout << ampersand << " " << predictions[i].total;
88			cout << " \\\\ \\hline" << endl;
89			cout << " " << ampersand << " Observed ";
90			for(int i=0;i<this->predictions.size();i++) cout << ampersand << " " << predictions[i].observed;
91			cout << " \\\\ \\hline" << endl;
92			cout << " " << ampersand << " Upper Limit";
93			for(int i=0;i<this->predictions.size();i++) cout << ampersand << " " << predictions[i].UpperLimit;
94			cout << " \\\\ \\hline";
95			for(int isignal=0;isignal<this->predictions[0].SignalYield.size();isignal++) {
96			cout << endl;
97			cout << " " << ampersand << " " << predictions[0].SignalYield[isignal].name;
98			for(int i=0;i<this->predictions.size();i++) {
99			if(predictions[i].SignalYield.size()<=isignal) cout << "X";
100			else cout << ampersand << " " << predictions[i].SignalYield[isignal].first << " +/- " << predictions[i].SignalYield[isignal].second << " \\\\";
101			}
102			}
103			cout << " \\hline " << endl;
104
105			}
106
107			/*
108			int main() {
109			LibraryForOneCut tester;
110			LeptonResult observed;
111			observed.ee=2;
112			observed.mm=3;
113			observed.eemm=5;
114			Yield Zbkg;Zbkg.value=3;Zbkg.staterror=2;Zbkg.syserror=1;
115			Yield Flavorsym;Flavorsym.value=6;Flavorsym.staterror=5;Flavorsym.syserror=4;
116			Yield total;total.value=9;total.staterror=8;total.syserror=7;
117			float UpperLimit=15;
118			tester.observed=observed;
119			tester.jzbcut=11;
120			tester.Zbkg=Zbkg;
121			tester.Flavorsym=Flavorsym;
122			tester.total=total;
123			tester.UpperLimit=UpperLimit;
124			ResultLibrary hello;
125			vector<tripple> signals;
126			tripple signal;
127			signal.name="LM4";
128			signal.first=2;
129			signal.second=3;
130			(tester.SignalYield).push_back(signal);
131			signal.name="LM8";
132			signal.first=4;
133			signal.second=5;
134			tester.SignalYield.push_back(signal);
135
136
137			hello.Add(tester);
138			tester.SignalYield.push_back(signal);
139
140			tester.observed.ee=12;
141			tester.observed.mm=13;
142			tester.observed.eemm=25;
143			tester.jzbcut=50;
144			hello.Add(tester);
145			tester.jzbcut=150;
146			hello.Add(tester);
147			tester.Print();
148			hello.Print();
149			}
150			*/