Skip to content

Commit

Permalink
Merge pull request #94 from kreczko/fixing_non_iso
Browse files Browse the repository at this point in the history
Fixes for issue #92
  • Loading branch information
senkin committed Feb 24, 2014
2 parents be9ee40 + e21efca commit 5b1348a
Show file tree
Hide file tree
Showing 15 changed files with 123 additions and 606 deletions.
1 change: 1 addition & 0 deletions .gitignore
Original file line number Diff line number Diff line change
Expand Up @@ -9,6 +9,7 @@ Makefile
*.py[cod]
*.log
*.so
*.root
Debug/
Profile/
Release/
Expand Down
23 changes: 14 additions & 9 deletions README.md
Original file line number Diff line number Diff line change
Expand Up @@ -5,23 +5,28 @@ Analysis software for TTbar+X differential cross section analysis (TOP-12-042)

## General Recipe

To setup the standalone version of the code:
To setup the stand-alone version of the code:

```
#get the code from the repository
# get the code from the repository
git clone https://github.com/BristolTopGroup/AnalysisSoftware.git AnalysisSoftware
cd AnalysisSoftware
#create the makefile using cmake (install it if you don't have it: https://github.com/Kitware/CMake)
# create the makefile using cmake (install it if you don't have it: https://github.com/Kitware/CMake)
cmake CMakeLists.txt
#compile and build
# for older versions of boost you will need to
#sed -i 's/lib64\/lib64/lib64/g' ./CMakeFiles/AnalysisSoftware.dir/build.make
#sed -i 's/lib64\/lib64/lib64/g' ./CMakeFiles/AnalysisSoftware.dir/link.txt
# as the library path is wrong (contains lib64/lib64)
# compile and build
make -j24
#adjust the toolsFolder environment variable
# adjust the toolsFolder environment variable
export toolsFolder=$PWD/
#test run the code:
# test run the code:
nohup ./AnalysisSoftware python/test_cfg.py &> test.log &
```

Expand All @@ -34,15 +39,15 @@ cmsrel CMSSW_6_0_0
cd CMSSW_6_0_0/src
cmsenv
#get the code from the repository
# get the code from the repository
git clone https://github.com/BristolTopGroup/AnalysisSoftware.git BristolAnalysis/Tools
#compile
# compile
scram b -j2
hash -r #or rehash in case that BAT cannot be found
#test run the code:
# test run the code:
nohup BAT BristolAnalysis/Tools/python/test_cfg.py &> test.log &
```

Expand Down
38 changes: 3 additions & 35 deletions interface/Analysers/TTbar_plus_X_analyser.h
Original file line number Diff line number Diff line change
Expand Up @@ -32,25 +32,18 @@ class TTbar_plus_X_analyser: public BAT::BasicAnalyser {
void muPlusJetsQcdAnalysis(const EventPtr);
void ePlusJetsSignalAnalysis(const EventPtr);
void muPlusJetsSignalAnalysis(const EventPtr);
//compare shift up/down for the systematic samples
//for the argument of 0-btag vs >=2 b-tag shape
// void ePlusJetsSystematicsAnalysis(const EventPtr);
// void muPlusJetsSystematicsAnalysis(const EventPtr);

private:
//signal selections
SelectionPointer topEplusJetsRefSelection_, topMuplusJetsRefSelection_;
//QCD selections with respect to reference selection
SelectionPointer qcdNonIsoElectronSelection_, qcdConversionSelection_;
SelectionPointer qcd_noniso_muon_plus_jets_selection_, qcd_noniso_muon_plus_jets_selection_ge4j_;
// SelectionPointer qcdPFRelIsoEPlusJetsSelection_, qcdPFRelIsoMuPlusJetsSelection_,
// qcd_noiso_muon_plus_jets_selection_;
SelectionPointer qcdPFRelIsoEPlusJetsSelection_, qcdPFRelIsoMuPlusJetsSelection_;

/**
* Analysers
*/
//no selection
// METAnalyserLocalPtr metAnalyserEPlusJetsNoSelection_, metAnalyserMuPlusJetsNoSelection_;
//signal regions
METAnalyserLocalPtr metAnalyserEPlusJetsRefSelection_, metAnalyserMuPlusJetsRefSelection_;
ElectronAnalyserLocalPtr electronAnalyserRefSelection_;
Expand All @@ -66,74 +59,49 @@ class TTbar_plus_X_analyser: public BAT::BasicAnalyser {
METAnalyserLocalPtr metAnalyserqcdConversionSelection_;
ElectronAnalyserLocalPtr qcdConversionsElectronAnalyser_;
//No iso
// ElectronAnalyserLocalPtr qcdEPlusjetsPFRelIsoElectronAnalyser_;
// MuonAnalyserLocalPtr qcdMuPlusjetsPFRelIsoMuonAnalyser_, qcdNoIsolationMuonAnalyser_;
ElectronAnalyserLocalPtr qcdEPlusjetsPFRelIsoElectronAnalyser_;
MuonAnalyserLocalPtr qcdMuPlusjetsPFRelIsoMuonAnalyser_;

std::vector<double> metBins_, ht_bins_, st_bins_, mt_bins_, wpt_bins_;
//MET analysers electron
// std::vector<Binned_Variable_analyser_ptr> no_selection_binned_MET_analyser_electron_eta_;
std::vector<Binned_Variable_analyser_ptr> ref_selection_binned_MET_analyser_electron_eta_;
std::vector<Binned_Variable_analyser_ptr> qcd_conversion_binned_MET_analyser_electron_eta_;
std::vector<Binned_Variable_analyser_ptr> qcd_noniso_binned_MET_analyser_electron_eta_;
// std::vector<Binned_Variable_analyser_ptr> qcd_PFRelIso_binned_MET_analyser_electron_isolation_; //for QCD rate estimation
// std::vector<Binned_Variable_analyser_ptr> qcd_PFRelIso_binned_MET_analyser_electron_rhocorrected_isolation_; //for QCD rate estimation
//MET analysers muon
// std::vector<Binned_Variable_analyser_ptr> no_selection_binned_MET_analyser_muon_eta_;
std::vector<Binned_Variable_analyser_ptr> ref_selection_binned_MET_analyser_muon_eta_;
std::vector<Binned_Variable_analyser_ptr> qcd_noniso_binned_MET_analyser_muon_eta_;
// std::vector<Binned_Variable_analyser_ptr> qcd_PFRelIso_binned_MET_analyser_muon_eta_;

//HT analysers electron
// Binned_Variable_analyser_ptr no_selection_binned_HT_analyser_electron_eta_;
Binned_Variable_analyser_ptr ref_selection_binned_HT_analyser_electron_eta_;
Binned_Variable_analyser_ptr qcd_conversion_binned_HT_analyser_electron_eta_;
Binned_Variable_analyser_ptr qcd_noniso_binned_HT_analyser_electron_eta_;
// Binned_Variable_analyser_ptr qcd_PFRelIso_binned_HT_analyser_electron_isolation_; //for QCD rate estimation
// Binned_Variable_analyser_ptr qcd_PFRelIso_binned_HT_analyser_electron_rhocorrected_isolation_; //for QCD rate estimation
//HT analysers muon
// Binned_Variable_analyser_ptr no_selection_binned_HT_analyser_muon_eta_;
Binned_Variable_analyser_ptr ref_selection_binned_HT_analyser_muon_eta_;
Binned_Variable_analyser_ptr qcd_noniso_binned_HT_analyser_muon_eta_;
// Binned_Variable_analyser_ptr qcd_PFRelIso_binned_HT_analyser_muon_eta_;

//ST analysers electron
// std::vector<Binned_Variable_analyser_ptr> no_selection_binned_ST_analyser_electron_eta_;
std::vector<Binned_Variable_analyser_ptr> ref_selection_binned_ST_analyser_electron_eta_;
std::vector<Binned_Variable_analyser_ptr> qcd_conversion_binned_ST_analyser_electron_eta_;
std::vector<Binned_Variable_analyser_ptr> qcd_noniso_binned_ST_analyser_electron_eta_;
// std::vector<Binned_Variable_analyser_ptr> qcd_PFRelIso_binned_ST_analyser_electron_isolation_; //for QCD rate estimation
// std::vector<Binned_Variable_analyser_ptr> qcd_PFRelIso_binned_ST_analyser_electron_rhocorrected_isolation_; //for QCD rate estimation
//ST analysers muon
// std::vector<Binned_Variable_analyser_ptr> no_selection_binned_ST_analyser_muon_eta_;
std::vector<Binned_Variable_analyser_ptr> ref_selection_binned_ST_analyser_muon_eta_;
std::vector<Binned_Variable_analyser_ptr> qcd_noniso_binned_ST_analyser_muon_eta_;
// std::vector<Binned_Variable_analyser_ptr> qcd_PFRelIso_binned_ST_analyser_muon_eta_;

//MT analysers electron
// std::vector<Binned_Variable_analyser_ptr> no_selection_binned_MT_analyser_electron_eta_;
std::vector<Binned_Variable_analyser_ptr> ref_selection_binned_MT_analyser_electron_eta_;
std::vector<Binned_Variable_analyser_ptr> qcd_conversion_binned_MT_analyser_electron_eta_;
std::vector<Binned_Variable_analyser_ptr> qcd_noniso_binned_MT_analyser_electron_eta_;
// std::vector<Binned_Variable_analyser_ptr> qcd_PFRelIso_binned_MT_analyser_electron_isolation_; //for QCD rate estimation
// std::vector<Binned_Variable_analyser_ptr> qcd_PFRelIso_binned_MT_analyser_electron_rhocorrected_isolation_; //for QCD rate estimation
//MT analysers muon
// std::vector<Binned_Variable_analyser_ptr> no_selection_binned_MT_analyser_muon_eta_;
std::vector<Binned_Variable_analyser_ptr> ref_selection_binned_MT_analyser_muon_eta_;
std::vector<Binned_Variable_analyser_ptr> qcd_noniso_binned_MT_analyser_muon_eta_;
// std::vector<Binned_Variable_analyser_ptr> qcd_PFRelIso_binned_MT_analyser_muon_eta_;

//WPT analysers electron
// std::vector<Binned_Variable_analyser_ptr> no_selection_binned_WPT_analyser_electron_eta_;
std::vector<Binned_Variable_analyser_ptr> ref_selection_binned_WPT_analyser_electron_eta_;
std::vector<Binned_Variable_analyser_ptr> qcd_conversion_binned_WPT_analyser_electron_eta_;
std::vector<Binned_Variable_analyser_ptr> qcd_noniso_binned_WPT_analyser_electron_eta_;
// std::vector<Binned_Variable_analyser_ptr> qcd_PFRelIso_binned_WPT_analyser_electron_isolation_; //for QCD rate estimation
// std::vector<Binned_Variable_analyser_ptr> qcd_PFRelIso_binned_WPT_analyser_electron_rhocorrected_isolation_; //for QCD rate estimation
//WPT analysers muon
// std::vector<Binned_Variable_analyser_ptr> no_selection_binned_WPT_analyser_muon_eta_;
std::vector<Binned_Variable_analyser_ptr> ref_selection_binned_WPT_analyser_muon_eta_;
std::vector<Binned_Variable_analyser_ptr> qcd_noniso_binned_WPT_analyser_muon_eta_;
// std::vector<Binned_Variable_analyser_ptr> qcd_PFRelIso_binned_WPT_analyser_muon_eta_;

JetAnalyserLocalPtr jetAnalyserEPlusJetsRefSelection_, jetAnalyserMuPlusJetsRefSelection_;

Expand Down
2 changes: 1 addition & 1 deletion interface/DataTypes.h
Original file line number Diff line number Diff line change
Expand Up @@ -173,7 +173,7 @@ const boost::array<std::string, DataType::NUMBER_OF_DATA_TYPES> names = { {
"QCD_Pt-1000_MuEnrichedPt5", //
"GJets_HT-40To100", //
"GJets_HT-100To200", //
"GJets_HT-200", //
"GJets_HT-200ToInf", //
//2012 PJ samples
"GJets_HT-200To400", //
"GJets_HT-400ToInf", //
Expand Down
2 changes: 1 addition & 1 deletion python/DataSetInfo.py
Original file line number Diff line number Diff line change
Expand Up @@ -33,7 +33,7 @@

datasetInfo['GJets_HT-40To100'] = {"cross-section": 23620., "NumberOfProcessedEvents":12730863}
datasetInfo['GJets_HT-100To200'] = {"cross-section": 3476., "NumberOfProcessedEvents":1536287}
datasetInfo['GJets_HT-200'] = {"cross-section": 485., "NumberOfProcessedEvents":9377168}
datasetInfo['GJets_HT-200ToInf'] = {"cross-section": 485., "NumberOfProcessedEvents":9377168}

datasetInfo['QCD_Pt-20to30_BCtoE'] = {"cross-section": 0.2355e9 * 0.00046, "NumberOfProcessedEvents":2002588}
datasetInfo['QCD_Pt-30to80_BCtoE'] = {"cross-section": 0.0593e9 * 0.00234, "NumberOfProcessedEvents":2030030}
Expand Down
2 changes: 1 addition & 1 deletion python/DataSetInfo_7TeV.py
Original file line number Diff line number Diff line change
Expand Up @@ -29,7 +29,7 @@

datasetInfo['GJets_HT-40To100'] = {"cross-section": 23620., "NumberOfProcessedEvents":12730169}
datasetInfo['GJets_HT-100To200'] = {"cross-section": 3476., "NumberOfProcessedEvents":1535665}
datasetInfo['GJets_HT-200'] = {"cross-section": 485., "NumberOfProcessedEvents":9366381}
datasetInfo['GJets_HT-200ToInf'] = {"cross-section": 485., "NumberOfProcessedEvents":9366381}

datasetInfo['QCD_Pt-20to30_BCtoE'] = {"cross-section": 0.2355e9 * 0.00046, "NumberOfProcessedEvents":2081553}
datasetInfo['QCD_Pt-30to80_BCtoE'] = {"cross-section": 0.0593e9 * 0.00234, "NumberOfProcessedEvents":2030029}
Expand Down
2 changes: 1 addition & 1 deletion python/DataSetInfo_V7.py
Original file line number Diff line number Diff line change
Expand Up @@ -33,7 +33,7 @@

datasetInfo['GJets_HT-40To100'] = {"cross-section": 23620., "NumberOfProcessedEvents":12658883}
datasetInfo['GJets_HT-100To200'] = {"cross-section": 3476., "NumberOfProcessedEvents":1520063}
datasetInfo['GJets_HT-200'] = {"cross-section": 485., "NumberOfProcessedEvents":9377168}
datasetInfo['GJets_HT-200ToInf'] = {"cross-section": 485., "NumberOfProcessedEvents":9377168}

datasetInfo['QCD_Pt-20to30_BCtoE'] = {"cross-section": 0.2355e9 * 0.00046, "NumberOfProcessedEvents":1796637}
datasetInfo['QCD_Pt-30to80_BCtoE'] = {"cross-section": 0.0593e9 * 0.00234, "NumberOfProcessedEvents":2030030}
Expand Down
6 changes: 3 additions & 3 deletions python/test_cfg.py
Original file line number Diff line number Diff line change
Expand Up @@ -5,7 +5,7 @@

#center of mass energy: 7TeV for 2010/2011 data/MC, 8TeV for 2012 data
#this value will be part of the output file name: DataType_CenterOfMassEnergyTeV_lumipb-1_....
centerOfMassEnergy = 7
centerOfMassEnergy = 8
nTuple_version = 10
#number of events to be processed
maxEvents = 10000
Expand Down Expand Up @@ -60,11 +60,11 @@
useHitFit = False
produceFitterASCIIoutput = False
inputFiles = [
# '/storage/TopQuarkGroup/data/2012/SingleElectron/nTuple_v10_Run2012D-22Jan2013-v1_AOD_GoldenJSON_LeptonPlus3Jets/*.root' # 53X 8TeV data
'/storage/TopQuarkGroup/data/2012/SingleElectron/nTuple_v10_Run2012D-22Jan2013-v1_AOD_GoldenJSON_LeptonPlus3Jets/*.root' # 53X 8TeV data
# '/storage/TopQuarkGroup/mc/8TeV/LeptonPlus3JetsSkim/TTJets_MassiveBinDECAY_TuneZ2star_8TeV-madgraph-tauola/nTuple_v10_Summer12_DR53X-PU_S10_START53_V7C-v1_LeptonPlus3Jets/*.root' #53X 8TeV MC
# '/storage/TopQuarkGroup/data/2011/ElectronHad/nTuple_v10_Run2011A-12Oct2013-v1_53X_GoldenJSON_LeptonPlus3Jets/*.root' # 53X 7TeV data
# ''53X 7TeV MC not available yet
'/storage/TopQuarkGroup/data/2011/ElectronHad/nTuple_v10_Run2011B-19Nov2011-v1_44X_GoldenJSON_LeptonPlus3Jets/*.root' #44X 7TeV data
# '/storage/TopQuarkGroup/data/2011/ElectronHad/nTuple_v10_Run2011B-19Nov2011-v1_44X_GoldenJSON_LeptonPlus3Jets/*.root' #44X 7TeV data
# '/storage/TopQuarkGroup/mc/7TeV/v10/LeptonPlus3JetsSkim/TTJets_TuneZ2_7TeV-madgraph-tauola/nTuple_v10_Fall11-PU_S6_START44_V9B-v1_LeptonPlus3Jets/*.root' #44X 7TeV MC
]

Expand Down
42 changes: 21 additions & 21 deletions src/Analysers/ElectronAnalyser.cpp
Original file line number Diff line number Diff line change
Expand Up @@ -22,9 +22,9 @@ void ElectronAnalyser::analyse(const EventPtr event) {
histMan_->H1D_BJetBinned("All_Electron_Eta")->Fill(electron->eta(), weight_);
histMan_->H1D_BJetBinned("All_Electron_AbsEta")->Fill(fabs(electron->eta()), weight_);
histMan_->H1D_BJetBinned("All_Electron_Phi")->Fill(electron->phi(), weight_);
histMan_->H1D_BJetBinned("All_Electron_pfIsolation_03")->Fill(electron->pfRelativeIsolation(0.3), weight_);
histMan_->H1D_BJetBinned("All_Electron_pfIsolation_04")->Fill(electron->pfRelativeIsolation(0.4), weight_);
histMan_->H1D_BJetBinned("All_Electron_pfIsolation_05")->Fill(electron->pfRelativeIsolation(0.5), weight_);
// histMan_->H1D_BJetBinned("All_Electron_pfIsolation_03")->Fill(electron->pfRelativeIsolation(0.3), weight_);
// histMan_->H1D_BJetBinned("All_Electron_pfIsolation_04")->Fill(electron->pfRelativeIsolation(0.4), weight_);
// histMan_->H1D_BJetBinned("All_Electron_pfIsolation_05")->Fill(electron->pfRelativeIsolation(0.5), weight_);
histMan_->H1D_BJetBinned("All_Electron_rhoCorrectedIso_03")->Fill(electron->pfRelativeIsolationRhoCorrected(), weight_);

histMan_->H1D_BJetBinned("All_Electron_sigma_ietaieta")->Fill(electron->sigmaIEtaIEta(), weight_);
Expand All @@ -44,15 +44,15 @@ void ElectronAnalyser::analyseElectron(const ElectronPointer electron, double we
return;
histMan_->H1D_BJetBinned("electron_eta")->Fill(electron->eta(), weight_);
histMan_->H1D_BJetBinned("electron_AbsEta")->Fill(fabs(electron->eta()), weight_);
histMan_->H1D_BJetBinned("electron_pfIsolation_03")->Fill(electron->pfRelativeIsolation(0.3), weight_);
histMan_->H1D_BJetBinned("electron_rhoCorrectedIso_03")->Fill(electron->pfRelativeIsolationRhoCorrected(), weight_);
// histMan_->H1D_BJetBinned("electron_pfIsolation_03")->Fill(electron->pfRelativeIsolation(0.3), weight_);

if (!ttbarPlusMETAnalysisSetup_) {
histMan_->H1D_BJetBinned("electron_pT")->Fill(electron->pt(), weight_);
histMan_->H1D_BJetBinned("electron_phi")->Fill(electron->phi(), weight_);

histMan_->H1D_BJetBinned("electron_pfIsolation_04")->Fill(electron->pfRelativeIsolation(0.4), weight_);
histMan_->H1D_BJetBinned("electron_pfIsolation_05")->Fill(electron->pfRelativeIsolation(0.5), weight_);
histMan_->H1D_BJetBinned("electron_rhoCorrectedIso_03")->Fill(electron->pfRelativeIsolationRhoCorrected(), weight_);
// histMan_->H1D_BJetBinned("electron_pfIsolation_04")->Fill(electron->pfRelativeIsolation(0.4), weight_);
// histMan_->H1D_BJetBinned("electron_pfIsolation_05")->Fill(electron->pfRelativeIsolation(0.5), weight_);

histMan_->H1D_BJetBinned("electron_sigma_ietaieta")->Fill(electron->sigmaIEtaIEta(), weight_);
histMan_->H1D_BJetBinned("electron_dPhi_in")->Fill(electron->dPhiIn(), weight_);
Expand Down Expand Up @@ -84,12 +84,12 @@ void ElectronAnalyser::createHistograms() {
histMan_->addH1D_BJetBinned("All_Electron_Eta", "Electron #eta; #eta(e); Events/(0.02)", 300, -3, 3);
histMan_->addH1D_BJetBinned("All_Electron_AbsEta", "Electron |#eta|; |#eta(e)|; Events/(0.01)", 300, 0, 3);
histMan_->addH1D_BJetBinned("All_Electron_Phi", "Electron #phi; #phi(e); Events/(0.02)", 400, -4, 4);
histMan_->addH1D_BJetBinned("All_Electron_pfIsolation_03",
"Electron relative pf isolation (DR=0.3); PF relative isolation; Events/(0.01)", 500, 0, 5);
histMan_->addH1D_BJetBinned("All_Electron_pfIsolation_04",
"Electron relative pf isolation (DR=0.4); PF relative isolation; Events/(0.01)", 500, 0, 5);
histMan_->addH1D_BJetBinned("All_Electron_pfIsolation_05",
"Electron relative pf isolation (DR=0.5); PF relative isolation; Events/(0.01)", 500, 0, 5);
// histMan_->addH1D_BJetBinned("All_Electron_pfIsolation_03",
// "Electron relative pf isolation (DR=0.3); PF relative isolation; Events/(0.01)", 500, 0, 5);
// histMan_->addH1D_BJetBinned("All_Electron_pfIsolation_04",
// "Electron relative pf isolation (DR=0.4); PF relative isolation; Events/(0.01)", 500, 0, 5);
// histMan_->addH1D_BJetBinned("All_Electron_pfIsolation_05",
// "Electron relative pf isolation (DR=0.5); PF relative isolation; Events/(0.01)", 500, 0, 5);
histMan_->addH1D_BJetBinned("All_Electron_rhoCorrectedIso_03",
"Electron rho-corrected isolation (DR=0.3); rho-corrected isolation; Events/(0.01)", 500, 0, 5);
histMan_->addH1D_BJetBinned("All_Electron_sigma_ietaieta",
Expand All @@ -110,12 +110,10 @@ void ElectronAnalyser::createHistograms() {

//b-jet binning useful for QCD estimates

histMan_->addH1D_BJetBinned("electron_pfIsolation_04",
"Electron relative pf isolation (DR=0.4); PF relative isolation; Events/(0.01)", 500, 0, 5);
histMan_->addH1D_BJetBinned("electron_pfIsolation_05",
"Electron relative pf isolation (DR=0.5); PF relative isolation; Events/(0.01)", 500, 0, 5);
histMan_->addH1D_BJetBinned("electron_rhoCorrectedIso_03",
"Electron rho-corrected isolation (DR=0.3); rho-corrected isolation; Events/(0.01)", 500, 0, 5);
// histMan_->addH1D_BJetBinned("electron_pfIsolation_04",
// "Electron relative pf isolation (DR=0.4); PF relative isolation; Events/(0.01)", 500, 0, 5);
// histMan_->addH1D_BJetBinned("electron_pfIsolation_05",
// "Electron relative pf isolation (DR=0.5); PF relative isolation; Events/(0.01)", 500, 0, 5);
histMan_->addH1D_BJetBinned("electron_sigma_ietaieta",
"Electron #sigma_{i#etai#eta}; #sigma_{i#etai#eta}; Events/(0.001)", 100, 0, 0.1);
histMan_->addH1D_BJetBinned("electron_dPhi_in", "Electron #Delta#Phi_{in}; #Delta#Phi_{in}; Events/(0.01)", 200, -1,
Expand All @@ -133,8 +131,10 @@ void ElectronAnalyser::createHistograms() {
if (singleElectronOnly_) {
histMan_->addH1D_BJetBinned("electron_eta", "Electron #eta; #eta(e); Events/(0.02)", 300, -3, 3);
histMan_->addH1D_BJetBinned("electron_AbsEta", "Electron |#eta|; |#eta(e)|; Events/(0.01)", 300, 0, 3);
histMan_->addH1D_BJetBinned("electron_pfIsolation_03",
"Electron relative pf isolation (DR=0.3); PF relative isolation; Events/(0.01)", 500, 0, 5);
histMan_->addH1D_BJetBinned("electron_rhoCorrectedIso_03",
"Electron rho-corrected isolation (DR=0.3); rho-corrected isolation; Events/(0.01)", 500, 0, 5);
// histMan_->addH1D_BJetBinned("electron_pfIsolation_03",
// "Electron relative pf isolation (DR=0.3); PF relative isolation; Events/(0.01)", 500, 0, 5);
}

}
Expand Down
Loading

0 comments on commit 5b1348a

Please sign in to comment.