ScansICHEP2012/code/getSStopLimit.C

{
float lumiErr   = 0.045;  // 4.5%
float trgErr    = 0.030;  // 3%
float isoErr    = 0.100;  // 5% per lepton
float btagErr   = 0.0714; // from cvs (Claudio)
float jesErr     = 0.05;  // from cvs (Claudio)
 float totErr = sqrt(lumiErr*lumiErr + trgErr*trgErr + isoErr*isoErr + btagErr*btagErr + jesErr*jesErr);
 float nev = observedLimitOnNumberOfEvents(3,100.*totErr,true);

 //float lumi=2870.; // in invpb
 float lumi = 3950.;
float acc=1.78445e-3; // from cvs (Claudio)
float xsecLimit = nev / (acc*lumi);
 float nevExpected = expectedLimitOnNumberOfEvents(3,100.*totErr,true);
 float xsecLimitExpected = nevExpected / (acc*lumi);


 cout << nev<< " " << xsecLimit << endl;
 cout << nevExpected << " " << xsecLimitExpected << endl;

}
Revision:	1.1
Committed:	Thu Jun 21 11:15:30 2012 UTC (12 years, 10 months ago) by claudioc
Content type:	text/plain
Branch:	MAIN
CVS Tags:	HEAD
Log Message:	3.95 invfb
#	User	Rev	Content
1	claudioc	1.1	{
2			float lumiErr = 0.045; // 4.5%
3			float trgErr = 0.030; // 3%
4			float isoErr = 0.100; // 5% per lepton
5			float btagErr = 0.0714; // from cvs (Claudio)
6			float jesErr = 0.05; // from cvs (Claudio)
7			float totErr = sqrt(lumiErrlumiErr + trgErrtrgErr + isoErrisoErr + btagErrbtagErr + jesErr*jesErr);
8			float nev = observedLimitOnNumberOfEvents(3,100.*totErr,true);
9
10			//float lumi=2870.; // in invpb
11			float lumi = 3950.;
12			float acc=1.78445e-3; // from cvs (Claudio)
13			float xsecLimit = nev / (acc*lumi);
14			float nevExpected = expectedLimitOnNumberOfEvents(3,100.*totErr,true);
15			float xsecLimitExpected = nevExpected / (acc*lumi);
16
17
18			cout << nev<< " " << xsecLimit << endl;
19			cout << nevExpected << " " << xsecLimitExpected << endl;
20
21			}