27 std::vector<std::string> subDirVec;
31 bool found_reco_dir =
false;
32 bool found_miniaod_dir =
false;
37 found_reco_dir =
true;
40 found_miniaod_dir =
true;
43 if (found_miniaod_dir && found_reco_dir) {
76 MonitorElement *mPtCorrOverReco_Eta_20_40_Reco = iget_.
get(rundir_reco +
"/" +
"PtCorrOverReco_Eta_20_40");
77 MonitorElement *mPtCorrOverReco_Eta_200_600_Reco = iget_.
get(rundir_reco +
"/" +
"PtCorrOverReco_Eta_200_600");
78 MonitorElement *mPtCorrOverReco_Eta_1500_3500_Reco = iget_.
get(rundir_reco +
"/" +
"PtCorrOverReco_Eta_1500_3500");
79 MonitorElement *mPtCorrOverGen_GenEta_40_200_Reco = iget_.
get(rundir_reco +
"/" +
"PtCorrOverGen_GenEta_40_200");
81 iget_.
get(rundir_reco +
"/" +
"PtCorrOverGen_GenEta_600_1500");
87 MonitorElement *mchargedHadronMultiplicity_Reco = iget_.
get(rundir_reco +
"/" +
"chargedHadronMultiplicity");
88 MonitorElement *mneutralHadronMultiplicity_Reco = iget_.
get(rundir_reco +
"/" +
"neutralHadronMultiplicity");
89 MonitorElement *mphotonMultiplicity_Reco = iget_.
get(rundir_reco +
"/" +
"photonMultiplicity");
90 MonitorElement *mphotonEnergyFraction_Reco = iget_.
get(rundir_reco +
"/" +
"photonEnergyFraction");
91 MonitorElement *mneutralHadronEnergyFraction_Reco = iget_.
get(rundir_reco +
"/" +
"neutralHadronEnergyFraction");
92 MonitorElement *mchargedHadronEnergyFraction_Reco = iget_.
get(rundir_reco +
"/" +
"chargedHadronEnergyFraction");
94 std::vector<MonitorElement *> ME_Reco;
95 ME_Reco.push_back(mGenPt_Reco);
96 ME_Reco.push_back(mGenPhi_Reco);
97 ME_Reco.push_back(mGenEta_Reco);
98 ME_Reco.push_back(mPt_Reco);
99 ME_Reco.push_back(mPhi_Reco);
100 ME_Reco.push_back(mEta_Reco);
101 ME_Reco.push_back(mCorrJetPt_Reco);
102 ME_Reco.push_back(mCorrJetPhi_Reco);
103 ME_Reco.push_back(mCorrJetEta_Reco);
104 ME_Reco.push_back(mPtCorrOverReco_Eta_20_40_Reco);
105 ME_Reco.push_back(mPtCorrOverReco_Eta_200_600_Reco);
106 ME_Reco.push_back(mPtCorrOverReco_Eta_1500_3500_Reco);
107 ME_Reco.push_back(mPtCorrOverGen_GenEta_40_200_Reco);
108 ME_Reco.push_back(mPtCorrOverGen_GenEta_600_1500_Reco);
109 ME_Reco.push_back(mDeltaEta_Reco);
110 ME_Reco.push_back(mDeltaPhi_Reco);
111 ME_Reco.push_back(mDeltaPt_Reco);
112 ME_Reco.push_back(mMjj_Reco);
113 ME_Reco.push_back(mNJets40_Reco);
114 ME_Reco.push_back(mchargedHadronMultiplicity_Reco);
115 ME_Reco.push_back(mneutralHadronMultiplicity_Reco);
116 ME_Reco.push_back(mphotonMultiplicity_Reco);
117 ME_Reco.push_back(mphotonEnergyFraction_Reco);
118 ME_Reco.push_back(mneutralHadronEnergyFraction_Reco);
119 ME_Reco.push_back(mchargedHadronEnergyFraction_Reco);
127 MonitorElement *mCorrJetPt_MiniAOD = iget_.
get(rundir_miniaod +
"/" +
"CorrJetPt");
128 MonitorElement *mCorrJetPhi_MiniAOD = iget_.
get(rundir_miniaod +
"/" +
"CorrJetPhi");
129 MonitorElement *mCorrJetEta_MiniAOD = iget_.
get(rundir_miniaod +
"/" +
"CorrJetEta");
130 MonitorElement *mPtCorrOverReco_Eta_20_40_MiniAOD = iget_.
get(rundir_miniaod +
"/" +
"PtCorrOverReco_Eta_20_40");
132 iget_.
get(rundir_miniaod +
"/" +
"PtCorrOverReco_Eta_200_600");
134 iget_.
get(rundir_miniaod +
"/" +
"PtCorrOverReco_Eta_1500_3500");
136 iget_.
get(rundir_miniaod +
"/" +
"PtCorrOverGen_GenEta_40_200");
138 iget_.
get(rundir_miniaod +
"/" +
"PtCorrOverGen_GenEta_600_1500");
139 MonitorElement *mDeltaEta_MiniAOD = iget_.
get(rundir_miniaod +
"/" +
"DeltaEta");
140 MonitorElement *mDeltaPhi_MiniAOD = iget_.
get(rundir_miniaod +
"/" +
"DeltaPhi");
144 MonitorElement *mchargedHadronMultiplicity_MiniAOD = iget_.
get(rundir_miniaod +
"/" +
"chargedHadronMultiplicity");
145 MonitorElement *mneutralHadronMultiplicity_MiniAOD = iget_.
get(rundir_miniaod +
"/" +
"neutralHadronMultiplicity");
146 MonitorElement *mphotonMultiplicity_MiniAOD = iget_.
get(rundir_miniaod +
"/" +
"photonMultiplicity");
147 MonitorElement *mphotonEnergyFraction_MiniAOD = iget_.
get(rundir_miniaod +
"/" +
"photonEnergyFraction");
149 iget_.
get(rundir_miniaod +
"/" +
"neutralHadronEnergyFraction");
151 iget_.
get(rundir_miniaod +
"/" +
"chargedHadronEnergyFraction");
153 std::vector<MonitorElement *> ME_MiniAOD;
154 ME_MiniAOD.push_back(mGenPt_MiniAOD);
155 ME_MiniAOD.push_back(mGenPhi_MiniAOD);
156 ME_MiniAOD.push_back(mGenEta_MiniAOD);
157 ME_MiniAOD.push_back(mPt_MiniAOD);
158 ME_MiniAOD.push_back(mPhi_MiniAOD);
159 ME_MiniAOD.push_back(mEta_MiniAOD);
160 ME_MiniAOD.push_back(mCorrJetPt_MiniAOD);
161 ME_MiniAOD.push_back(mCorrJetPhi_MiniAOD);
162 ME_MiniAOD.push_back(mCorrJetEta_MiniAOD);
163 ME_MiniAOD.push_back(mPtCorrOverReco_Eta_20_40_MiniAOD);
164 ME_MiniAOD.push_back(mPtCorrOverReco_Eta_200_600_MiniAOD);
165 ME_MiniAOD.push_back(mPtCorrOverReco_Eta_1500_3500_MiniAOD);
166 ME_MiniAOD.push_back(mPtCorrOverGen_GenEta_40_200_MiniAOD);
167 ME_MiniAOD.push_back(mPtCorrOverGen_GenEta_600_1500_MiniAOD);
168 ME_MiniAOD.push_back(mDeltaEta_MiniAOD);
169 ME_MiniAOD.push_back(mDeltaPhi_MiniAOD);
170 ME_MiniAOD.push_back(mDeltaPt_MiniAOD);
171 ME_MiniAOD.push_back(mMjj_MiniAOD);
172 ME_MiniAOD.push_back(mNJets40_MiniAOD);
173 ME_MiniAOD.push_back(mchargedHadronMultiplicity_MiniAOD);
174 ME_MiniAOD.push_back(mneutralHadronMultiplicity_MiniAOD);
175 ME_MiniAOD.push_back(mphotonMultiplicity_MiniAOD);
176 ME_MiniAOD.push_back(mphotonEnergyFraction_MiniAOD);
177 ME_MiniAOD.push_back(mneutralHadronEnergyFraction_MiniAOD);
178 ME_MiniAOD.push_back(mchargedHadronEnergyFraction_MiniAOD);
195 float etarange[91] = {-6.0, -5.8, -5.6, -5.4, -5.2, -5.0, -4.8, -4.6, -4.4, -4.2, -4.0, -3.8, -3.6,
196 -3.4, -3.2, -3.0, -2.9, -2.8, -2.7, -2.6, -2.5, -2.4, -2.3, -2.2, -2.1, -2.0,
197 -1.9, -1.8, -1.7, -1.6, -1.5, -1.4, -1.3, -1.2, -1.1, -1.0, -0.9, -0.8, -0.7,
198 -0.6, -0.5, -0.4, -0.3, -0.2, -0.1, 0.0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6,
199 0.7, 0.8, 0.9, 1.0, 1.1, 1.2, 1.3, 1.4, 1.5, 1.6, 1.7, 1.8, 1.9,
200 2.0, 2.1, 2.2, 2.3, 2.4, 2.5, 2.6, 2.7, 2.8, 2.9, 3.0, 3.2, 3.4,
201 3.6, 3.8, 4.0, 4.2, 4.4, 4.6, 4.8, 5.0, 5.2, 5.4, 5.6, 5.8, 6.0};
204 ibook_.
book1D(
"PtCorrOverReco_Eta_20_40_MiniAOD_over_RECO",
"20<genpt<40", 90, etarange);
206 ibook_.
book1D(
"PtCorrOverReco_Eta_200_600_MiniAOD_over_RECO",
"200<genpt<600", 90, etarange);
208 ibook_.
book1D(
"PtCorrOverReco_Eta_1500_3500_MiniAOD_over_RECO",
"1500<genpt<3500", 90, etarange);
210 ibook_.
book1D(
"PtCorrOverGen_GenEta_40_200_MiniAOD_over_RECO",
"40<genpt<200", 90, etarange);
212 ibook_.
book1D(
"PtCorrOverGen_GenEta_600_1500_MiniAOD_over_RECO",
"600<genpt<1500", 90, etarange);
219 "chargedHadronMultiplicity_MiniAOD_over_RECO", (TH1F *)mchargedHadronMultiplicity_Reco->
getRootObject());
221 "neutralHadronMultiplicity_MiniAOD_over_RECO", (TH1F *)mneutralHadronMultiplicity_Reco->
getRootObject());
223 ibook_.
book1D(
"photonMultiplicity_MiniAOD_over_RECO", (TH1F *)mphotonMultiplicity_Reco->
getRootObject());
225 "chargedHadronEnergyFraction_MiniAOD_over_RECO", (TH1F *)mchargedHadronEnergyFraction_Reco->
getRootObject());
227 "neutralHadronEnergyFraction_MiniAOD_over_RECO", (TH1F *)mneutralHadronEnergyFraction_Reco->
getRootObject());
229 ibook_.
book1D(
"photonEnergyFraction_MiniAOD_over_RECO", (TH1F *)mphotonEnergyFraction_Reco->
getRootObject());
231 std::vector<MonitorElement *> ME_MiniAOD_over_Reco;
257 for (
unsigned int j = 0;
j < ME_MiniAOD_over_Reco.size();
j++) {
263 if (monMiniAOD_over_RECO && monMiniAOD_over_RECO->
getRootObject()) {
264 for (
int i = 0;
i <= (monMiniAOD_over_RECO->
getNbinsX() + 1);
i++) {