UserArea/StatTools/countingLimits.C

#include "cl95cms.C"

//const double lumi = 34.18, lumiErr = 0.052;
const double lumi = 34.74, lumiErr = 0.11;

void run(double qTcut, double qTcutErr, double accCut, double accCutErr, double effCut, double effCutErr, double bgExpect, double bgExpectErr, int bg){
  // Full efficiency is the product -> error propagation (assuing absolute errors) is:
  double uncertainty = qTcut*accCut*effCut * sqrt( qTcutErr/qTcut*qTcutErr/qTcut + accCutErr/accCut*accCutErr/accCut + effCutErr/effCut*effCutErr/effCut );
  double uncertNoAcc = qTcut*       effCut * sqrt( qTcutErr/qTcut*qTcutErr/qTcut +                                     effCutErr/effCut*effCutErr/effCut );

  cout<<"Efficiency= "<<qTcut*accCut*effCut<<" +- "<<uncertainty<<endl;

  double expectedLimitSigma = CLA (lumi, lumi*lumiErr, qTcut*accCut*effCut, uncertainty, bgExpect, bgExpectErr);
  double actualLimitSigma   = CL95(lumi, lumi*lumiErr, qTcut*accCut*effCut, uncertainty, bgExpect, bgExpectErr, bg);

  double expectedLimitSigmaAcc = CLA (lumi, lumi*lumiErr, qTcut*effCut, uncertNoAcc, bgExpect, bgExpectErr);
  double actualLimitSigmaAcc   = CL95(lumi, lumi*lumiErr, qTcut*effCut, uncertNoAcc, bgExpect, bgExpectErr, bg);

  cout<<"Found limit on sigma="<<actualLimitSigma<<" vs. expected limit on sigma="<<expectedLimitSigma<<" (pb)"<<endl;
  cout<<"Found limit on sigma x acc="<<actualLimitSigmaAcc<<" vs. expected limit on sigma x acc="<<expectedLimitSigmaAcc<<" (pb)"<<endl;
}

void countingLimits(void){
/*
  double qTCI05     = 0.3466;      // qT threshold efficiency (generated with boostedZ/cl95cms/cuts.C)
  double qTCI05err  = 0.006;       //  and the absolute statistical uncertainty
  double accCI05    = 0.589/0.839; // acceptance above qT threshold (from GP's table)
  double accCI05err = sqrt(0.007*0.007+0.02*0.02*accCI05*accCI05);   //  and the absolute statistical (0.009/sqrt(2)) [+] systematic uncertainty (from Joe's plots)
  double effCI05    = 0.839;       // selection efficiency above the qT threshold within the acceptance
  double effCI05err = sqrt(0.007*0.007+0.015*0.015*effCI05*effCI05); //  and the absolute statistical (0.009/sqrt(2)) [+] systematic uncertainty (my calcs)
  //double bgExpect300    = 5.3; // From GP's data-driven bg estimation
  //double bgExpectErr300 = 2.7;
  double bgExpect300    = 3.0;
  double bgExpectErr300 = 2.6;
  double bgMeasured300  = 5;

  cout<<"Contact Interaction M= 500 GeV/c"<<endl;
  run(qTCI05,qTCI05err,accCI05,accCI05err,effCI05,effCI05err,bgExpect300,bgExpectErr300,bgMeasured300);
  cout<<endl<<endl;
*/
/*
  double qTCI10     = 0.478943;    // qT threshold efficiency (generated with boostedZ/cl95cms/cuts.C)
  double qTCI10err  = 0.006;       //  and the relative statistical uncertainty
  double accCI10    = 0.582/0.840; // acceptance above qT threshold (from GP's table)
  double accCI10err = sqrt(0.006*0.006+0.02*accCI10*0.02*accCI10);   //  and the relative statistical (0.008/sqrt(2)) [+] systematic uncertainty (from Joe's plots)
  double effCI10    = 0.840;       // selection efficiency above the qT threshold within the acceptance
  double effCI10err = sqrt(0.006*0.006+0.015*effCI10*0.015*effCI10); //  and the relative statistical (0.008/sqrt(2)) [+] systematic uncertainty (my calcs)
  //double bgExpect350    = 2.4; // From GP's data-driven bg estimation
  //double bgExpectErr350 = 1.4;
  double bgExpect350    = 1.3;
  double bgExpectErr350 = 1.25;
  double bgMeasured350  = 1;

  cout<<"Contact Interaction M = 1000 GeV/c"<<endl;
  run(qTCI10,qTCI10err,accCI10,accCI10err,effCI10,effCI10err,bgExpect350,bgExpectErr350,bgMeasured350);
  cout<<endl<<endl;

  double qTCI15     = 0.544992;    // qT threshold efficiency (generated with boostedZ/cl95cms/cuts.C)
  double qTCI15err  = 0.006;       //  and the relative statistical uncertainty
  double accCI15    = 0.512/0.839; // acceptance above qT threshold (from GP's table)
  double accCI15err = sqrt(0.005*0.005+0.02*accCI15*0.02*accCI15);   //  and the relative statistical (0.007/sqrt(2)) [+] systematic uncertainty (from Joe's plots)
  double effCI15    = 0.839;       // selection efficiency above the qT threshold within the acceptance
  double effCI15err = sqrt(0.005*0.005+0.015*effCI15*0.015*effCI15); //  and the relative statistical (0.007/sqrt(2)) [+] systematic uncertainty (my calcs)
  //double bgExpect400    = 1.2; // From GP's data-driven bg estimation
  //double bgExpectErr400 = 0.7;
  double bgExpect400    = 0.57;
  double bgExpectErr400 = 0.58;
  double bgMeasured400  = 1;

  cout<<"Contact Interaction M = 1500 GeV/c"<<endl;
  run(qTCI15,qTCI15err,accCI15,accCI15err,effCI15,effCI15err,bgExpect400,bgExpectErr400,bgMeasured400);
  cout<<endl<<endl;

  double qTCI20     = 0.581699;    // qT threshold efficiency (generated with boostedZ/cl95cms/cuts.C)
  double qTCI20err  = 0.006;       //  and the relative statistical uncertainty
  double accCI20    = 0.443/0.832; // acceptance above qT threshold (from GP's table)
  double accCI20err = sqrt(0.005*0.005+0.02*accCI20*0.02*accCI20);   //  and the relative statistical (0.007/sqrt(2)) [+] systematic uncertainty (from Joe's plots)
  double effCI20    = 0.832;       // selection efficiency above the qT threshold within the acceptance
  double effCI20err = sqrt(0.005*0.005+0.015*effCI20*0.015*effCI20); //  and the relative statistical (0.007/sqrt(2)) [+] systematic uncertainty (my calcs)

  cout<<"Contact Interaction M = 2000 GeV/c"<<endl;
  run(qTCI20,qTCI20err,accCI20,accCI20err,effCI20,effCI20err,bgExpect400,bgExpectErr400,bgMeasured400);
  cout<<endl<<endl;
*/
  double qTGI05     = 0.434767;    // qT threshold efficiency (generated with boostedZ/cl95cms/cuts.C)
  double qTGI05err  = 0.007;       //  and the absolute statistical uncertainty
  double accGI05    = 0.704/0.831; // acceptance above qT threshold (from GP's table)
  double accGI05err = sqrt(0.007*0.007+0.02*0.02*accGI05*accGI05);   //  and the absolute statistical (0.009/sqrt(2)) [+] systematic uncertainty (from Joe's plots)
  double effGI05    = 0.831;       // selection efficiency above the qT threshold within the acceptance
  double effGI05err = sqrt(0.007*0.007+0.015*0.015*effGI05*effGI05); //  and the absolute statistical (0.009/sqrt(2)) [+] systematic uncertainty (my calcs)
  //double bgExpect200    = 29; // From GP's data-driven bg estimation
  //double bgExpectErr200 = 9;
  double bgExpect200    = 20.;
  double bgExpectErr200 = 10.1;
  double bgMeasured200  = 16;

  cout<<"Gauge Interaction M= 500 GeV/c"<<endl;
  run(qTGI05,qTGI05err,accGI05,accGI05err,effGI05,effGI05err,bgExpect200,bgExpectErr200,bgMeasured200);
  cout<<endl<<endl;
return;
/*
  double qTGI07     = 0.71325;    // qT threshold efficiency (generated with boostedZ/cl95cms/cuts.C)
  double qTGI07err  = 0.006;       //  and the absolute statistical uncertainty
  double accGI07    = 0.698/0.831; // acceptance above qT threshold (from GP's table)
  double accGI07err = sqrt(0.005*0.005+0.02*0.02*accGI07*accGI07);   //  and the absolute statistical (0.007/sqrt(2)) [+] systematic uncertainty (from Joe's plots)
  double effGI07    = 0.831;       // selection efficiency above the qT threshold within the acceptance
  double effGI07err = sqrt(0.005*0.005+0.015*0.015*effGI07*effGI07); //  and the absolute statistical (0.007/sqrt(2)) [+] systematic uncertainty (my calcs)
  double bgExpect200    = 29; // From GP's data-driven bg estimation
  double bgExpectErr200 = 9;
  double bgMeasured200  = 16;

  cout<<"Gauge Interaction M= 750 GeV/c"<<endl;
  run(qTGI07,qTGI07err,accGI07,accGI07err,effGI07,effGI07err,bgExpect200,bgExpectErr200,bgMeasured200);
  cout<<endl<<endl;
*/
  double qTGI10     = 0.552582;    // qT threshold efficiency (generated with boostedZ/cl95cms/cuts.C)
  double qTGI10err  = 0.007;       //  and the absolute statistical uncertainty
  double accGI10    = 0.701/0.865; // acceptance above qT threshold (from GP's table)
  double accGI10err = sqrt(0.005*0.005+0.02*0.02*accGI10*accGI10);   //  and the absolute statistical (0.007/sqrt(2)) [+] systematic uncertainty (from Joe's plots)
  double effGI10    = 0.865;       // selection efficiency above the qT threshold within the acceptance
  double effGI10err = sqrt(0.005*0.005+0.015*0.015*effGI10*effGI10); //  and the absolute statistical (0.007/sqrt(2)) [+] systematic uncertainty (my calcs)
  //double bgExpect350    = 2.4; // From GP's data-driven bg estimation
  //double bgExpectErr350 = 1.4;
  double bgExpect350    = 1.3;
  double bgExpectErr350 = 1.25;
  double bgMeasured350  = 1;

  cout<<"Gauge Interaction M= 1000 GeV/c"<<endl;
  run(qTGI10,qTGI10err,accGI10,accGI10err,effGI10,effGI10err,bgExpect350,bgExpectErr350,bgMeasured350);
  cout<<endl<<endl;
/*
  double qTGI12     = 0.650356;    // qT threshold efficiency (generated with boostedZ/cl95cms/cuts.C)
  double qTGI12err  = 0.006;       //  and the absolute statistical uncertainty
  double accGI12    = 0.673/0.864; // acceptance above qT threshold (from GP's table)
  double accGI12err = sqrt(0.005*0.005+0.02*0.02*accGI12*accGI12);   //  and the absolute statistical (0.007/sqrt(2)) [+] systematic uncertainty (from Joe's plots)
  double effGI12    = 0.864;       // selection efficiency above the qT threshold within the acceptance
  double effGI12err = sqrt(0.005*0.005+0.015*0.015*effGI12*effGI12); //  and the absolute statistical (0.007/sqrt(2)) [+] systematic uncertainty (my calcs)
  double bgExpect350    = 2.4; // From GP's data-driven bg estimation
  double bgExpectErr350 = 1.4;
  double bgMeasured350  = 1;

  cout<<"Gauge Interaction M= 1250 GeV/c"<<endl;
  run(qTGI12,qTGI12err,accGI12,accGI12err,effGI12,effGI12err,bgExpect350,bgExpectErr350,bgMeasured350);
  cout<<endl<<endl;
*/
  double qTGI15     = 0.607149;    // qT threshold efficiency (generated with boostedZ/cl95cms/cuts.C)
  double qTGI15err  = 0.006;       //  and the absolute statistical uncertainty
  double accGI15    = 0.721/0.879; // acceptance above qT threshold (from GP's table)
  double accGI15err = sqrt(0.016*0.016+0.02*0.02*accGI15*accGI15);   //  and the absolute statistical (0.007/sqrt(2)) [+] systematic uncertainty (from Joe's plots)
  double effGI15    = 0.879;       // selection efficiency above the qT threshold within the acceptance
  double effGI15err = sqrt(0.016*0.016+0.015*0.015*effGI15*effGI15); //  and the absolute statistical (0.007/sqrt(2)) [+] systematic uncertainty (my calcs)
  //double bgExpect400    = 1.2; // From GP's data-driven bg estimation
  //double bgExpectErr400 = 0.7;
  double bgExpect400    = 0.57;
  double bgExpectErr400 = 0.58;
  double bgMeasured400  = 1;

  cout<<"Gauge Interaction M= 1500 GeV/c"<<endl;
  run(qTGI15,qTGI15err,accGI15,accGI15err,effGI15,effGI15err,bgExpect400,bgExpectErr400,bgMeasured400);
  cout<<endl<<endl;
/*
  double qTZP05     = 0.320869;    // qT threshold efficiency (generated with boostedZ/cl95cms/cuts.C)
  double qTZP05err  = 0.0114643;   //  and the absolute statistical uncertainty
  double accZP05    = 0.479/0.734; // acceptance above qT threshold (from GP's table)
  double accZP05err = sqrt(0.013*0.013+0.02*0.02*accZP05*accZP05);   //  and the absolute statistical (0.019/sqrt(2)) [+] systematic uncertainty (from Joe's plots)
  double effZP05    = 0.734;       // selection efficiency above the qT threshold within the acceptance
  double effZP05err = sqrt(0.013*0.013+0.015*0.015*effZP05*effZP05); //  and the absolute statistical (0.019/sqrt(2)) [+] systematic uncertainty (my calcs)
//  double bgExpect200    = 29; // From GP's data-driven bg estimation
//  double bgExpectErr200 = 9;
//  double bgMeasured200  = 16;

  cout<<"Z' M= 500 GeV/c"<<endl;
  run(qTZP05,qTZP05err,accZP05,accZP05err,effZP05,effZP05err,bgExpect200,bgExpectErr200,bgMeasured200);
  cout<<endl<<endl;

  double qTZP10     = 0.475628;    // qT threshold efficiency (generated with boostedZ/cl95cms/cuts.C)
  double qTZP10err  = 0.00683797;  //  and the absolute statistical uncertainty
  double accZP10    = 0.499/0.822; // acceptance above qT threshold (from GP's table)
  double accZP10err = sqrt(0.006*0.006+0.02*0.02*accZP10*accZP10);   //  and the absolute statistical (0.008/sqrt(2)) [+] systematic uncertainty (from Joe's plots)
  double effZP10    = 0.822;       // selection efficiency above the qT threshold within the acceptance
  double effZP10err = sqrt(0.006*0.006+0.015*0.015*effZP10*effZP10); //  and the absolute statistical (0.008/sqrt(2)) [+] systematic uncertainty (my calcs)
//  double bgExpect350    = 2.4; // From GP's data-driven bg estimation
//  double bgExpectErr350 = 1.4;
//  double bgMeasured350  = 1;

  cout<<"Z' M= 1000 GeV/c"<<endl;
  run(qTZP10,qTZP10err,accZP10,accZP10err,effZP10,effZP10err,bgExpect350,bgExpectErr350,bgMeasured350);
  cout<<endl<<endl;

  double qTZP15     = 0.540415;    // qT threshold efficiency (generated with boostedZ/cl95cms/cuts.C)
  double qTZP15err  = 0.00802145;  //  and the absolute statistical uncertainty
  double accZP15    = 0.439/0.822; // acceptance above qT threshold (from GP's table)
  double accZP15err = sqrt(0.007*0.007+0.02*0.02*accZP15*accZP15);   //  and the absolute statistical (0.010/sqrt(2)) [+] systematic uncertainty (from Joe's plots)
  double effZP15    = 0.822;       // selection efficiency above the qT threshold within the acceptance
  double effZP15err = sqrt(0.007*0.007+0.015*0.015*effZP15*effZP15); //  and the absolute statistical (0.010/sqrt(2)) [+] systematic uncertainty (my calcs)
//  double bgExpect400    = 1.2; // From GP's data-driven bg estimation
//  double bgExpectErr400 = 0.7;
//  double bgMeasured400  = 1;

  cout<<"Z' M= 1500 GeV/c"<<endl;
  run(qTZP15,qTZP15err,accZP15,accZP15err,effZP15,effZP15err,bgExpect400,bgExpectErr400,bgMeasured400);
  cout<<endl<<endl;
*/
}
Revision:	1.1
Committed:	Tue Feb 15 14:03:43 2011 UTC (14 years, 2 months ago) by kkotov
Content type:	text/plain
Branch:	MAIN
CVS Tags:	V2012-H-02, V2012-01-00, V2011-01-01, V2011-01-00, AnnaDimuon, V01-00-01, V01-00-00, HEAD
Error occurred while calculating annotation data.
Log Message:	* empty log message *
#	Content
1	#include "cl95cms.C"
2
3	//const double lumi = 34.18, lumiErr = 0.052;
4	const double lumi = 34.74, lumiErr = 0.11;
5
6	void run(double qTcut, double qTcutErr, double accCut, double accCutErr, double effCut, double effCutErr, double bgExpect, double bgExpectErr, int bg){
7	// Full efficiency is the product -> error propagation (assuing absolute errors) is:
8	double uncertainty = qTcutaccCuteffCut * sqrt( qTcutErr/qTcutqTcutErr/qTcut + accCutErr/accCutaccCutErr/accCut + effCutErr/effCut*effCutErr/effCut );
9	double uncertNoAcc = qTcut* effCut * sqrt( qTcutErr/qTcutqTcutErr/qTcut + effCutErr/effCuteffCutErr/effCut );
10
11	cout<<"Efficiency= "<<qTcutaccCuteffCut<<" +- "<<uncertainty<<endl;
12
13	double expectedLimitSigma = CLA (lumi, lumilumiErr, qTcutaccCut*effCut, uncertainty, bgExpect, bgExpectErr);
14	double actualLimitSigma = CL95(lumi, lumilumiErr, qTcutaccCut*effCut, uncertainty, bgExpect, bgExpectErr, bg);
15
16	double expectedLimitSigmaAcc = CLA (lumi, lumilumiErr, qTcuteffCut, uncertNoAcc, bgExpect, bgExpectErr);
17	double actualLimitSigmaAcc = CL95(lumi, lumilumiErr, qTcuteffCut, uncertNoAcc, bgExpect, bgExpectErr, bg);
18
19	cout<<"Found limit on sigma="<<actualLimitSigma<<" vs. expected limit on sigma="<<expectedLimitSigma<<" (pb)"<<endl;
20	cout<<"Found limit on sigma x acc="<<actualLimitSigmaAcc<<" vs. expected limit on sigma x acc="<<expectedLimitSigmaAcc<<" (pb)"<<endl;
21	}
22
23	void countingLimits(void){
24	/*
25	double qTCI05 = 0.3466; // qT threshold efficiency (generated with boostedZ/cl95cms/cuts.C)
26	double qTCI05err = 0.006; // and the absolute statistical uncertainty
27	double accCI05 = 0.589/0.839; // acceptance above qT threshold (from GP's table)
28	double accCI05err = sqrt(0.0070.007+0.020.02accCI05accCI05); // and the absolute statistical (0.009/sqrt(2)) [+] systematic uncertainty (from Joe's plots)
29	double effCI05 = 0.839; // selection efficiency above the qT threshold within the acceptance
30	double effCI05err = sqrt(0.0070.007+0.0150.015effCI05effCI05); // and the absolute statistical (0.009/sqrt(2)) [+] systematic uncertainty (my calcs)
31	//double bgExpect300 = 5.3; // From GP's data-driven bg estimation
32	//double bgExpectErr300 = 2.7;
33	double bgExpect300 = 3.0;
34	double bgExpectErr300 = 2.6;
35	double bgMeasured300 = 5;
36
37	cout<<"Contact Interaction M= 500 GeV/c"<<endl;
38	run(qTCI05,qTCI05err,accCI05,accCI05err,effCI05,effCI05err,bgExpect300,bgExpectErr300,bgMeasured300);
39	cout<<endl<<endl;
40	*/
41	/*
42	double qTCI10 = 0.478943; // qT threshold efficiency (generated with boostedZ/cl95cms/cuts.C)
43	double qTCI10err = 0.006; // and the relative statistical uncertainty
44	double accCI10 = 0.582/0.840; // acceptance above qT threshold (from GP's table)
45	double accCI10err = sqrt(0.0060.006+0.02accCI100.02accCI10); // and the relative statistical (0.008/sqrt(2)) [+] systematic uncertainty (from Joe's plots)
46	double effCI10 = 0.840; // selection efficiency above the qT threshold within the acceptance
47	double effCI10err = sqrt(0.0060.006+0.015effCI100.015effCI10); // and the relative statistical (0.008/sqrt(2)) [+] systematic uncertainty (my calcs)
48	//double bgExpect350 = 2.4; // From GP's data-driven bg estimation
49	//double bgExpectErr350 = 1.4;
50	double bgExpect350 = 1.3;
51	double bgExpectErr350 = 1.25;
52	double bgMeasured350 = 1;
53
54	cout<<"Contact Interaction M = 1000 GeV/c"<<endl;
55	run(qTCI10,qTCI10err,accCI10,accCI10err,effCI10,effCI10err,bgExpect350,bgExpectErr350,bgMeasured350);
56	cout<<endl<<endl;
57
58	double qTCI15 = 0.544992; // qT threshold efficiency (generated with boostedZ/cl95cms/cuts.C)
59	double qTCI15err = 0.006; // and the relative statistical uncertainty
60	double accCI15 = 0.512/0.839; // acceptance above qT threshold (from GP's table)
61	double accCI15err = sqrt(0.0050.005+0.02accCI150.02accCI15); // and the relative statistical (0.007/sqrt(2)) [+] systematic uncertainty (from Joe's plots)
62	double effCI15 = 0.839; // selection efficiency above the qT threshold within the acceptance
63	double effCI15err = sqrt(0.0050.005+0.015effCI150.015effCI15); // and the relative statistical (0.007/sqrt(2)) [+] systematic uncertainty (my calcs)
64	//double bgExpect400 = 1.2; // From GP's data-driven bg estimation
65	//double bgExpectErr400 = 0.7;
66	double bgExpect400 = 0.57;
67	double bgExpectErr400 = 0.58;
68	double bgMeasured400 = 1;
69
70	cout<<"Contact Interaction M = 1500 GeV/c"<<endl;
71	run(qTCI15,qTCI15err,accCI15,accCI15err,effCI15,effCI15err,bgExpect400,bgExpectErr400,bgMeasured400);
72	cout<<endl<<endl;
73
74	double qTCI20 = 0.581699; // qT threshold efficiency (generated with boostedZ/cl95cms/cuts.C)
75	double qTCI20err = 0.006; // and the relative statistical uncertainty
76	double accCI20 = 0.443/0.832; // acceptance above qT threshold (from GP's table)
77	double accCI20err = sqrt(0.0050.005+0.02accCI200.02accCI20); // and the relative statistical (0.007/sqrt(2)) [+] systematic uncertainty (from Joe's plots)
78	double effCI20 = 0.832; // selection efficiency above the qT threshold within the acceptance
79	double effCI20err = sqrt(0.0050.005+0.015effCI200.015effCI20); // and the relative statistical (0.007/sqrt(2)) [+] systematic uncertainty (my calcs)
80
81	cout<<"Contact Interaction M = 2000 GeV/c"<<endl;
82	run(qTCI20,qTCI20err,accCI20,accCI20err,effCI20,effCI20err,bgExpect400,bgExpectErr400,bgMeasured400);
83	cout<<endl<<endl;
84	*/
85	double qTGI05 = 0.434767; // qT threshold efficiency (generated with boostedZ/cl95cms/cuts.C)
86	double qTGI05err = 0.007; // and the absolute statistical uncertainty
87	double accGI05 = 0.704/0.831; // acceptance above qT threshold (from GP's table)
88	double accGI05err = sqrt(0.0070.007+0.020.02accGI05accGI05); // and the absolute statistical (0.009/sqrt(2)) [+] systematic uncertainty (from Joe's plots)
89	double effGI05 = 0.831; // selection efficiency above the qT threshold within the acceptance
90	double effGI05err = sqrt(0.0070.007+0.0150.015effGI05effGI05); // and the absolute statistical (0.009/sqrt(2)) [+] systematic uncertainty (my calcs)
91	//double bgExpect200 = 29; // From GP's data-driven bg estimation
92	//double bgExpectErr200 = 9;
93	double bgExpect200 = 20.;
94	double bgExpectErr200 = 10.1;
95	double bgMeasured200 = 16;
96
97	cout<<"Gauge Interaction M= 500 GeV/c"<<endl;
98	run(qTGI05,qTGI05err,accGI05,accGI05err,effGI05,effGI05err,bgExpect200,bgExpectErr200,bgMeasured200);
99	cout<<endl<<endl;
100	return;
101	/*
102	double qTGI07 = 0.71325; // qT threshold efficiency (generated with boostedZ/cl95cms/cuts.C)
103	double qTGI07err = 0.006; // and the absolute statistical uncertainty
104	double accGI07 = 0.698/0.831; // acceptance above qT threshold (from GP's table)
105	double accGI07err = sqrt(0.0050.005+0.020.02accGI07accGI07); // and the absolute statistical (0.007/sqrt(2)) [+] systematic uncertainty (from Joe's plots)
106	double effGI07 = 0.831; // selection efficiency above the qT threshold within the acceptance
107	double effGI07err = sqrt(0.0050.005+0.0150.015effGI07effGI07); // and the absolute statistical (0.007/sqrt(2)) [+] systematic uncertainty (my calcs)
108	double bgExpect200 = 29; // From GP's data-driven bg estimation
109	double bgExpectErr200 = 9;
110	double bgMeasured200 = 16;
111
112	cout<<"Gauge Interaction M= 750 GeV/c"<<endl;
113	run(qTGI07,qTGI07err,accGI07,accGI07err,effGI07,effGI07err,bgExpect200,bgExpectErr200,bgMeasured200);
114	cout<<endl<<endl;
115	*/
116	double qTGI10 = 0.552582; // qT threshold efficiency (generated with boostedZ/cl95cms/cuts.C)
117	double qTGI10err = 0.007; // and the absolute statistical uncertainty
118	double accGI10 = 0.701/0.865; // acceptance above qT threshold (from GP's table)
119	double accGI10err = sqrt(0.0050.005+0.020.02accGI10accGI10); // and the absolute statistical (0.007/sqrt(2)) [+] systematic uncertainty (from Joe's plots)
120	double effGI10 = 0.865; // selection efficiency above the qT threshold within the acceptance
121	double effGI10err = sqrt(0.0050.005+0.0150.015effGI10effGI10); // and the absolute statistical (0.007/sqrt(2)) [+] systematic uncertainty (my calcs)
122	//double bgExpect350 = 2.4; // From GP's data-driven bg estimation
123	//double bgExpectErr350 = 1.4;
124	double bgExpect350 = 1.3;
125	double bgExpectErr350 = 1.25;
126	double bgMeasured350 = 1;
127
128	cout<<"Gauge Interaction M= 1000 GeV/c"<<endl;
129	run(qTGI10,qTGI10err,accGI10,accGI10err,effGI10,effGI10err,bgExpect350,bgExpectErr350,bgMeasured350);
130	cout<<endl<<endl;
131	/*
132	double qTGI12 = 0.650356; // qT threshold efficiency (generated with boostedZ/cl95cms/cuts.C)
133	double qTGI12err = 0.006; // and the absolute statistical uncertainty
134	double accGI12 = 0.673/0.864; // acceptance above qT threshold (from GP's table)
135	double accGI12err = sqrt(0.0050.005+0.020.02accGI12accGI12); // and the absolute statistical (0.007/sqrt(2)) [+] systematic uncertainty (from Joe's plots)
136	double effGI12 = 0.864; // selection efficiency above the qT threshold within the acceptance
137	double effGI12err = sqrt(0.0050.005+0.0150.015effGI12effGI12); // and the absolute statistical (0.007/sqrt(2)) [+] systematic uncertainty (my calcs)
138	double bgExpect350 = 2.4; // From GP's data-driven bg estimation
139	double bgExpectErr350 = 1.4;
140	double bgMeasured350 = 1;
141
142	cout<<"Gauge Interaction M= 1250 GeV/c"<<endl;
143	run(qTGI12,qTGI12err,accGI12,accGI12err,effGI12,effGI12err,bgExpect350,bgExpectErr350,bgMeasured350);
144	cout<<endl<<endl;
145	*/
146	double qTGI15 = 0.607149; // qT threshold efficiency (generated with boostedZ/cl95cms/cuts.C)
147	double qTGI15err = 0.006; // and the absolute statistical uncertainty
148	double accGI15 = 0.721/0.879; // acceptance above qT threshold (from GP's table)
149	double accGI15err = sqrt(0.0160.016+0.020.02accGI15accGI15); // and the absolute statistical (0.007/sqrt(2)) [+] systematic uncertainty (from Joe's plots)
150	double effGI15 = 0.879; // selection efficiency above the qT threshold within the acceptance
151	double effGI15err = sqrt(0.0160.016+0.0150.015effGI15effGI15); // and the absolute statistical (0.007/sqrt(2)) [+] systematic uncertainty (my calcs)
152	//double bgExpect400 = 1.2; // From GP's data-driven bg estimation
153	//double bgExpectErr400 = 0.7;
154	double bgExpect400 = 0.57;
155	double bgExpectErr400 = 0.58;
156	double bgMeasured400 = 1;
157
158	cout<<"Gauge Interaction M= 1500 GeV/c"<<endl;
159	run(qTGI15,qTGI15err,accGI15,accGI15err,effGI15,effGI15err,bgExpect400,bgExpectErr400,bgMeasured400);
160	cout<<endl<<endl;
161	/*
162	double qTZP05 = 0.320869; // qT threshold efficiency (generated with boostedZ/cl95cms/cuts.C)
163	double qTZP05err = 0.0114643; // and the absolute statistical uncertainty
164	double accZP05 = 0.479/0.734; // acceptance above qT threshold (from GP's table)
165	double accZP05err = sqrt(0.0130.013+0.020.02accZP05accZP05); // and the absolute statistical (0.019/sqrt(2)) [+] systematic uncertainty (from Joe's plots)
166	double effZP05 = 0.734; // selection efficiency above the qT threshold within the acceptance
167	double effZP05err = sqrt(0.0130.013+0.0150.015effZP05effZP05); // and the absolute statistical (0.019/sqrt(2)) [+] systematic uncertainty (my calcs)
168	// double bgExpect200 = 29; // From GP's data-driven bg estimation
169	// double bgExpectErr200 = 9;
170	// double bgMeasured200 = 16;
171
172	cout<<"Z' M= 500 GeV/c"<<endl;
173	run(qTZP05,qTZP05err,accZP05,accZP05err,effZP05,effZP05err,bgExpect200,bgExpectErr200,bgMeasured200);
174	cout<<endl<<endl;
175
176	double qTZP10 = 0.475628; // qT threshold efficiency (generated with boostedZ/cl95cms/cuts.C)
177	double qTZP10err = 0.00683797; // and the absolute statistical uncertainty
178	double accZP10 = 0.499/0.822; // acceptance above qT threshold (from GP's table)
179	double accZP10err = sqrt(0.0060.006+0.020.02accZP10accZP10); // and the absolute statistical (0.008/sqrt(2)) [+] systematic uncertainty (from Joe's plots)
180	double effZP10 = 0.822; // selection efficiency above the qT threshold within the acceptance
181	double effZP10err = sqrt(0.0060.006+0.0150.015effZP10effZP10); // and the absolute statistical (0.008/sqrt(2)) [+] systematic uncertainty (my calcs)
182	// double bgExpect350 = 2.4; // From GP's data-driven bg estimation
183	// double bgExpectErr350 = 1.4;
184	// double bgMeasured350 = 1;
185
186	cout<<"Z' M= 1000 GeV/c"<<endl;
187	run(qTZP10,qTZP10err,accZP10,accZP10err,effZP10,effZP10err,bgExpect350,bgExpectErr350,bgMeasured350);
188	cout<<endl<<endl;
189
190	double qTZP15 = 0.540415; // qT threshold efficiency (generated with boostedZ/cl95cms/cuts.C)
191	double qTZP15err = 0.00802145; // and the absolute statistical uncertainty
192	double accZP15 = 0.439/0.822; // acceptance above qT threshold (from GP's table)
193	double accZP15err = sqrt(0.0070.007+0.020.02accZP15accZP15); // and the absolute statistical (0.010/sqrt(2)) [+] systematic uncertainty (from Joe's plots)
194	double effZP15 = 0.822; // selection efficiency above the qT threshold within the acceptance
195	double effZP15err = sqrt(0.0070.007+0.0150.015effZP15effZP15); // and the absolute statistical (0.010/sqrt(2)) [+] systematic uncertainty (my calcs)
196	// double bgExpect400 = 1.2; // From GP's data-driven bg estimation
197	// double bgExpectErr400 = 0.7;
198	// double bgMeasured400 = 1;
199
200	cout<<"Z' M= 1500 GeV/c"<<endl;
201	run(qTZP15,qTZP15err,accZP15,accZP15err,effZP15,effZP15err,bgExpect400,bgExpectErr400,bgMeasured400);
202	cout<<endl<<endl;
203	*/
204	}