Skip to content

Commit 49696eb

Browse files
author
wteo
committed
update to latest lumi and sideband
1 parent 1aa0e6a commit 49696eb

6 files changed

+4075
-61
lines changed

Diff for: ‎ElHad_300_75.cpp

+979
Large diffs are not rendered by default.

Diff for: ‎ElHad_300_80.cpp

+979
Large diffs are not rendered by default.

Diff for: ‎EventCalculator.C

+153-61
Original file line numberDiff line numberDiff line change
@@ -748,52 +748,49 @@ float EventCalculator::getHLTMHTeff(float offMET, float offHT, uint nElectrons,
748748
//eff = gr->Eval(offMET);
749749

750750
//for 2011 full result
751-
if(offHT>400 && offMET>=150 && offMET<200){
751+
//updated with pixelLumiCalc numbers
752+
if(offHT>400 && offMET>=150 && offMET<250){
752753
if( is0L && mindphin > 4 ){
753-
eff = 0.832;
754+
eff = 0.850;
754755
}
755756
else if( isSingleE && mindphin > 4){
756-
eff = 0.931;
757+
eff = 0.955;
757758
}
758759
else if( isSingleMu && mindphin > 4){
759-
eff = 0.983;
760+
eff = 0.990;
760761
}
761762
else if ( is0L && mindphin <= 4){
762-
eff = 0.910;
763+
eff = 0.912;
763764
}
764-
765-
766-
}
767-
if(offHT>400 && offMET>=200 && offMET<250){
768-
// eff = 0.94;
769-
//eff = 0.859362; //after HT cut
770-
771-
//errors are wrong!
772-
773-
////assign a 3% error in the SB region
774-
//if ( theHLTEffType_ ==kHLTEffup){
775-
// eff = eff*1.03; assert(0);
776-
//}
777-
//else if ( theHLTEffType_ ==kHLTEffdown){
778-
// eff = eff*0.97; assert(0);
779-
//}
780-
781-
//after HT>400, ==0b
782-
if( is0L && mindphin > 4 ){
783-
eff = 0.841;
784-
}
785-
else if( isSingleE && mindphin > 4){
786-
eff = 0.991;
787-
}
788-
else if( isSingleMu && mindphin > 4){
789-
eff = 1.0;
790-
}
791-
else if ( is0L && mindphin <= 4){
792-
eff = 0.936;
793-
}
794-
795-
796765
}
766+
//if(offHT>400 && offMET>=200 && offMET<250){
767+
// // eff = 0.94;
768+
// //eff = 0.859362; //after HT cut
769+
//
770+
// //errors are wrong!
771+
//
772+
// ////assign a 3% error in the SB region
773+
// //if ( theHLTEffType_ ==kHLTEffup){
774+
// // eff = eff*1.03; assert(0);
775+
// //}
776+
// //else if ( theHLTEffType_ ==kHLTEffdown){
777+
// // eff = eff*0.97; assert(0);
778+
// //}
779+
//
780+
// //after HT>400, ==0b
781+
// if( is0L && mindphin > 4 ){
782+
// eff = 0.841;
783+
// }
784+
// else if( isSingleE && mindphin > 4){
785+
// eff = 0.991;
786+
// }
787+
// else if( isSingleMu && mindphin > 4){
788+
// eff = 1.0;
789+
// }
790+
// else if ( is0L && mindphin <= 4){
791+
// eff = 0.936;
792+
// }
793+
//}
797794
else if(offHT>400 && offMET>250){
798795
// eff = 0.998;
799796
//eff = 0.975299; //after HT cut
@@ -807,18 +804,19 @@ float EventCalculator::getHLTMHTeff(float offMET, float offHT, uint nElectrons,
807804
//}
808805

809806
//after HT>400, ==0b
807+
//updated with pixelLumiCalc numbers
810808
if( is0L && mindphin > 4 ){
811-
eff = 0.982;
809+
eff = 0.981;
812810
}
813811
else if( isSingleE && mindphin > 4){
814812
eff = 1.0;
815813
}
816814
else if( isSingleMu && mindphin > 4){
817-
eff = 0.998;
815+
eff = 0.999;
818816
}
819817

820818
else if ( is0L && mindphin <= 4){
821-
eff = 0.982; // same as nominal for now
819+
eff = 0.981; // same as nominal for now
822820
}
823821

824822

@@ -900,7 +898,10 @@ float EventCalculator::getHLTHTeff(float offHT) {
900898
}
901899

902900
//2011 dataset - weighted average of HT efficiencies
903-
eff = 0.076*1.0 + 0.464*gr1->Eval(offHT) + 0.460*gr2->Eval(offHT);
901+
double frac_HT2XX = 0.071;
902+
double frac_HT300 = 0.462;
903+
double frac_HT350 = 0.467;
904+
eff = frac_HT2XX*1.0 + frac_HT300*gr1->Eval(offHT) + frac_HT350*gr2->Eval(offHT);
904905

905906
return eff;
906907
}
@@ -3401,18 +3402,18 @@ double EventCalculator::getHLTMHTeffBNN(float offMET, float offHT, uint nElectro
34013402
bool is0L = (nElectrons==0 && nMuons==0);
34023403

34033404
std::vector<double> vlumi;
3404-
vlumi.push_back(0.0067);//ht260_mht60
3405-
vlumi.push_back(0.0409);//ht250_mht60_v2
3406-
vlumi.push_back(0.1686);//ht250_mht60_v3
3407-
vlumi.push_back(0.1391);//ht250_mht70
3408-
vlumi.push_back(0.0982+0.4225);//ht300_mht75_v7
3409-
vlumi.push_back(0.0043);//ht300_mht75_v8
3410-
vlumi.push_back(0.2657);//ht300_mht80_v1
3411-
vlumi.push_back(0.7804);//ht300_mht80_v2
3412-
vlumi.push_back(0.6044);//ht300_mht90
3413-
vlumi.push_back(1.3405);//ht350_mht95
3414-
vlumi.push_back(0.8125);//ht350_mht110
3415-
double totallumi = 4.6838; //this is exactly the sum of the above lumi numbers
3405+
vlumi.push_back(0.0063);//ht260_mht60
3406+
vlumi.push_back(0.0407);//ht250_mht60_v2
3407+
vlumi.push_back(0.1687);//ht250_mht60_v3
3408+
vlumi.push_back(0.1364);//ht250_mht70
3409+
vlumi.push_back(0.1005+0.4366);//ht300_mht75_v7
3410+
vlumi.push_back(0.0044);//ht300_mht75_v8
3411+
vlumi.push_back(0.2773);//ht300_mht80_v1
3412+
vlumi.push_back(0.8313);//ht300_mht80_v2
3413+
vlumi.push_back(0.6522);//ht300_mht90
3414+
vlumi.push_back(1.4421);//ht350_mht95
3415+
vlumi.push_back(0.8855);//ht350_mht110
3416+
double totallumi = 4.982; //this is exactly the sum of the above lumi numbers
34163417

34173418
// BNNs:
34183419
std::vector<double> pfmet;
@@ -3841,16 +3842,62 @@ double EventCalculator::getHLTMHTeffBNN(float offMET, float offHT, uint nElectro
38413842
else if( offHT>400 && isSingleE && mindphin > 4){
38423843

38433844
std::vector<double> vlumi_sl;
3844-
vlumi_sl.push_back(0.6044);//ht300_mht90
3845-
vlumi_sl.push_back(1.3405);//ht350_mht95
3846-
vlumi_sl.push_back(0.8125);//ht350_mht110
3847-
double totallumi_sl = 2.7574; //this is exactly the sum of the above lumi numbers
3845+
vlumi_sl.push_back(0.0063);//ht260_mht60
3846+
vlumi_sl.push_back(0.0407+0.1687);//ht250_mht60
3847+
vlumi_sl.push_back(0.1364);//ht250_mht70
3848+
vlumi_sl.push_back(0.1005+0.0044+0.4366);//ht300_mht75
3849+
vlumi_sl.push_back(0.2773+0.8313);//ht300_mht80
3850+
vlumi_sl.push_back(0.6522);//ht300_mht90
3851+
vlumi_sl.push_back(1.4421);//ht350_mht95
3852+
vlumi_sl.push_back(0.8855);//ht350_mht110
3853+
double totallumi_sl = 4.982; //this is exactly the sum of the above lumi numbers
38483854

38493855
std::vector<double> trigeffs;
38503856
std::vector<double> trigeffserr;
38513857
double stddeveff = 0, efferr = 0;
3858+
3859+
trigeffs.push_back(ElHad_300_75( pfmet ));
3860+
stddeveff = 0;
3861+
for(uint i = 0; i < 100; ++i){
3862+
effn[i] = ElHad_300_75(pfmet, i, i);
3863+
stddeveff += (effn[i] - trigeffs.back())*(effn[i] - trigeffs.back());
3864+
}
3865+
stddeveff = sqrt(stddeveff / 99); trigeffserr.push_back(stddeveff);
3866+
3867+
trigeffs.push_back(ElHad_300_75( pfmet ));
3868+
stddeveff = 0;
3869+
for(uint i = 0; i < 100; ++i){
3870+
effn[i] = ElHad_300_75(pfmet, i, i);
3871+
stddeveff += (effn[i] - trigeffs.back())*(effn[i] - trigeffs.back());
3872+
}
3873+
stddeveff = sqrt(stddeveff / 99); trigeffserr.push_back(stddeveff);
3874+
3875+
trigeffs.push_back(ElHad_300_75( pfmet ));
3876+
stddeveff = 0;
3877+
for(uint i = 0; i < 100; ++i){
3878+
effn[i] = ElHad_300_75(pfmet, i, i);
3879+
stddeveff += (effn[i] - trigeffs.back())*(effn[i] - trigeffs.back());
3880+
}
3881+
stddeveff = sqrt(stddeveff / 99); trigeffserr.push_back(stddeveff);
3882+
3883+
trigeffs.push_back(ElHad_300_75( pfmet ));
3884+
stddeveff = 0;
3885+
for(uint i = 0; i < 100; ++i){
3886+
effn[i] = ElHad_300_75(pfmet, i, i);
3887+
stddeveff += (effn[i] - trigeffs.back())*(effn[i] - trigeffs.back());
3888+
}
3889+
stddeveff = sqrt(stddeveff / 99); trigeffserr.push_back(stddeveff);
3890+
3891+
trigeffs.push_back(ElHad_300_80( pfmet ));
3892+
stddeveff = 0;
3893+
for(uint i = 0; i < 100; ++i){
3894+
effn[i] = ElHad_300_80(pfmet, i, i);
3895+
stddeveff += (effn[i] - trigeffs.back())*(effn[i] - trigeffs.back());
3896+
}
3897+
stddeveff = sqrt(stddeveff / 99); trigeffserr.push_back(stddeveff);
38523898

38533899
trigeffs.push_back(HT300MHT90elNew( pfmet ));
3900+
stddeveff = 0;
38543901
for(uint i = 0; i < 100; ++i){
38553902
effn[i] = HT300MHT90elNew(pfmet, i, i);
38563903
stddeveff += (effn[i] - trigeffs.back())*(effn[i] - trigeffs.back());
@@ -3902,15 +3949,60 @@ double EventCalculator::getHLTMHTeffBNN(float offMET, float offHT, uint nElectro
39023949
else if( offHT>400 && isSingleMu && mindphin > 4){
39033950

39043951
std::vector<double> vlumi_sl;
3905-
vlumi_sl.push_back(0.6044);//ht300_mht90
3906-
vlumi_sl.push_back(1.3405);//ht350_mht95
3907-
vlumi_sl.push_back(0.8125);//ht350_mht110
3908-
double totallumi_sl = 2.7574; //this is exactly the sum of the above lumi numbers
3952+
vlumi_sl.push_back(0.0063);//ht260_mht60
3953+
vlumi_sl.push_back(0.0407+0.1687);//ht250_mht60
3954+
vlumi_sl.push_back(0.1364);//ht250_mht70
3955+
vlumi_sl.push_back(0.1005+0.0044+0.4366);//ht300_mht75
3956+
vlumi_sl.push_back(0.2773+0.8313);//ht300_mht80
3957+
vlumi_sl.push_back(0.6522);//ht300_mht90
3958+
vlumi_sl.push_back(1.4421);//ht350_mht95
3959+
vlumi_sl.push_back(0.8855);//ht350_mht110
3960+
double totallumi_sl = 4.982; //this is exactly the sum of the above lumi numbers
39093961

39103962
std::vector<double> trigeffs;
39113963
std::vector<double> trigeffserr;
39123964
double stddeveff = 0, efferr = 0;
39133965

3966+
trigeffs.push_back(MuHad_300_75( pfmet ));
3967+
stddeveff = 0;
3968+
for(uint i = 0; i < 100; ++i){
3969+
effn[i] = MuHad_300_75(pfmet, i, i);
3970+
stddeveff += (effn[i] - trigeffs.back())*(effn[i] - trigeffs.back());
3971+
}
3972+
stddeveff = sqrt(stddeveff / 99); trigeffserr.push_back(stddeveff);
3973+
3974+
trigeffs.push_back(MuHad_300_75( pfmet ));
3975+
stddeveff = 0;
3976+
for(uint i = 0; i < 100; ++i){
3977+
effn[i] = MuHad_300_75(pfmet, i, i);
3978+
stddeveff += (effn[i] - trigeffs.back())*(effn[i] - trigeffs.back());
3979+
}
3980+
stddeveff = sqrt(stddeveff / 99); trigeffserr.push_back(stddeveff);
3981+
3982+
trigeffs.push_back(MuHad_300_75( pfmet ));
3983+
stddeveff = 0;
3984+
for(uint i = 0; i < 100; ++i){
3985+
effn[i] = MuHad_300_75(pfmet, i, i);
3986+
stddeveff += (effn[i] - trigeffs.back())*(effn[i] - trigeffs.back());
3987+
}
3988+
stddeveff = sqrt(stddeveff / 99); trigeffserr.push_back(stddeveff);
3989+
3990+
trigeffs.push_back(MuHad_300_75( pfmet ));
3991+
stddeveff = 0;
3992+
for(uint i = 0; i < 100; ++i){
3993+
effn[i] = MuHad_300_75(pfmet, i, i);
3994+
stddeveff += (effn[i] - trigeffs.back())*(effn[i] - trigeffs.back());
3995+
}
3996+
stddeveff = sqrt(stddeveff / 99); trigeffserr.push_back(stddeveff);
3997+
3998+
trigeffs.push_back(MuHad_300_80( pfmet ));
3999+
stddeveff = 0;
4000+
for(uint i = 0; i < 100; ++i){
4001+
effn[i] = MuHad_300_80(pfmet, i, i);
4002+
stddeveff += (effn[i] - trigeffs.back())*(effn[i] - trigeffs.back());
4003+
}
4004+
stddeveff = sqrt(stddeveff / 99); trigeffserr.push_back(stddeveff);
4005+
39144006
trigeffs.push_back(HT300MHT90muNew( pfmet ));
39154007
for(uint i = 0; i < 100; ++i){
39164008
effn[i] = HT300MHT90muNew(pfmet, i, i);

Diff for: ‎EventCalculator.h

+6
Original file line numberDiff line numberDiff line change
@@ -55,6 +55,12 @@ double HT300_MHT80_v2_HT300_v9_170065_173211_sel2(std::vector<double>& inputvars
5555
double HT300_MHT90_v2_HT300_v9_173212_176544_sel2(std::vector<double>& inputvars, int first=0, int last=100-1);
5656
double HT350_MHT90_v1_HT350_v8_176545_178410_sel2(std::vector<double>& inputvars, int first=0, int last=100-1);
5757
double HT350_MHT110_v3_HT350_v11_178411_180252_sel2(std::vector<double>& inputvars, int first=0, int last=100-1);
58+
59+
60+
double ElHad_300_75(std::vector<double>& inputvars, int first=0, int last=100-1);
61+
double MuHad_300_75(std::vector<double>& inputvars, int first=0, int last=100-1);
62+
double ElHad_300_80(std::vector<double>& inputvars, int first=0, int last=100-1);
63+
double MuHad_300_80(std::vector<double>& inputvars, int first=0, int last=100-1);
5864
double HT300MHT90elNew(std::vector<double>& inputvars, int first=0, int last=100-1);
5965
double HT300MHT90muNew(std::vector<double>& inputvars, int first=0, int last=100-1);
6066
double HT350MHT90elNew(std::vector<double>& inputvars, int first=0, int last=100-1);

0 commit comments

Comments
 (0)