9 |
|
|
10 |
|
namespace HbbAnalysis {//namespace |
11 |
|
|
12 |
< |
void HistosMET::Initialise(TFileDirectory & aDir, std::string aName, bool aDoGenMatched){ |
13 |
< |
|
14 |
< |
p_MET = aDir.make<TH1F>("p_MET",";MET (GeV.c^{-1});N_{entries}",200,0,200); |
15 |
< |
p_MEx = aDir.make<TH1F>("p_MEx",";MEX (GeV.c^{-1});N_{entries}",200,-100,100); |
16 |
< |
p_MEy = aDir.make<TH1F>("p_MEy",";MEY (GeV.c^{-1});N_{entries}",200,-100,100); |
17 |
< |
p_SumET = aDir.make<TH1F>("p_SumET",";SumET (GeV.c^{-1});N_{entries}",500,0,500); |
18 |
< |
p_phi = aDir.make<TH1F>("p_phi",";#phi;N_{entries}",64,-3.2,3.2); |
19 |
< |
p_mEtSig = aDir.make<TH1F>("p_mEtSig",";MET/#sqrt{SumET};N_{entries}",100,0,10); |
20 |
< |
|
21 |
< |
p_genMET_MET = aDir.make<TH1F>("p_genMET_MET",";MET (GeV.c^{-1});N_{entries}",200,0,200); |
22 |
< |
p_genMET_MEx = aDir.make<TH1F>("p_genMET_MEx",";MEX (GeV.c^{-1});N_{entries}",200,-100,100); |
23 |
< |
p_genMET_MEy = aDir.make<TH1F>("p_genMET_MEy",";MEY (GeV.c^{-1});N_{entries}",200,-100,100); |
24 |
< |
p_genMET_SumET = aDir.make<TH1F>("p_genMET_SumET",";SumET (GeV.c^{-1});N_{entries}",500,0,500); |
25 |
< |
p_genMET_phi = aDir.make<TH1F>("p_genMET_phi",";#phi;N_{entries}",64,-3.2,3.2); |
26 |
< |
p_genMET_mEtSig = aDir.make<TH1F>("p_genMET_mEtSig",";MET/#sqrt{SumET};N_{entries}",100,0,10); |
27 |
< |
|
28 |
< |
p_RecoOverGen_MET = aDir.make<TH1F>("p_RecoOverGen_MET",";MET^{reco}/MET^{gen};N_{entries}",500,0,100); |
29 |
< |
p_RecoOverGen_MEx = aDir.make<TH1F>("p_RecoOverGen_MEx",";MEx^{reco}/MEx^{gen};N_{entries}",500,-50,50); |
30 |
< |
p_RecoOverGen_MEy = aDir.make<TH1F>("p_RecoOverGen_MEy",";MEy^{reco}/MEy^{gen};N_{entries}",500,-50,50); |
31 |
< |
p_RecoOverGen_SumET = aDir.make<TH1F>("p_RecoOverGen_SumET",";SumET^{reco}/SumET^{gen};N_{entries}",500,0,10); |
32 |
< |
|
33 |
< |
} |
34 |
< |
|
35 |
< |
void HistosMET::FillEventHistograms(edm::Handle<std::vector<pat::MET> > aCol) |
12 |
> |
void HistosMET::FillEventHistograms(const edm::Handle<std::vector<pat::MET> > & aCol) |
13 |
|
{ |
14 |
|
assert(aCol->size() == 1); |
15 |
|
const pat::MET & lMet = *(aCol->begin()); |