| 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 |
– |
|
| 12 |
|
void HistosMET::FillEventHistograms(const edm::Handle<std::vector<pat::MET> > & aCol) |
| 13 |
|
{ |
| 14 |
|
assert(aCol->size() == 1); |
