132 h_zmm_mass = fs->
make<TH1D>(
"zmm_mass",
"zmumu mass", 100, 0., 200.);
133 h_zmm2HLT_mass = fs->make<TH1D>(
"zmm2HLT_mass",
"zmumu 2HLT mass", 100, 0., 200.);
134 h_zmm1HLTplus_mass = fs->make<TH1D>(
"zmm1HLTplus_mass",
"zmumu 1HLT plus mass", 100, 0., 200.);
135 h_zmmNotIsoplus_mass = fs->make<TH1D>(
"zmmNotIsoplus_mass",
"zmumu a least One Not Iso plus mass", 100, 0., 200.);
136 h_zmsplus_mass = fs->make<TH1D>(
"zmsplus_mass",
"zmusta plus mass", 100, 0., 200.);
137 h_zmtplus_mass = fs->make<TH1D>(
"zmtplus_mass",
"zmutrack plus mass", 100, 0., 200.);
138 h_zmm1HLTminus_mass = fs->make<TH1D>(
"zmm1HLTminus_mass",
"zmumu 1HLT minus mass", 100, 0., 200.);
139 h_zmmNotIsominus_mass = fs->make<TH1D>(
"zmmNotIsominus_mass",
"zmumu a least One Not Iso minus mass", 100, 0., 200.);
140 h_zmsminus_mass = fs->make<TH1D>(
"zmsminus_mass",
"zmusta minus mass", 100, 0., 200.);
141 h_zmtminus_mass = fs->make<TH1D>(
"zmtminus_mass",
"zmutrack minus mass", 100, 0., 200.);
143 cout <<
"primo" << endl;
147 fs->
mkdir(
"etaIntervals");
149 fs->
mkdir(
"ptIntervals");
154 double etaRangeTmp[7] = {-2., -1.2, -0.8, 0., 0.8, 1.2, 2.};
155 double ptRangeTmp[5] = {20., 40., 60., 80., 100.};
158 for (
unsigned int i = 0;
i <=
ptBins;
i++)
162 cout <<
"eta istograms creation " << endl;
165 cout <<
" bin eta plus " <<
i << endl;
170 ap =
"zmumu2HLTplus_etaRange" + std::to_string(
i);
171 bp =
"zmumu2HLT plus mass eta Range " + std::to_string(range0) +
" to " + std::to_string(range1);
172 cout << ap <<
" " << bp << endl;
174 ap =
"zmumu1HLTplus_etaRange" + std::to_string(
i);
175 bp =
"zmumu1HLT plus mass eta Range " + std::to_string(range0) +
" to " + std::to_string(range1);
176 cout << ap <<
" " << bp << endl;
178 ap =
"zmustaplus_etaRange" + std::to_string(
i);
179 bp =
"zmusta plus mass eta Range " + std::to_string(range0) +
" to " + std::to_string(range1);
180 cout << ap <<
" " << bp << endl;
181 hmustaplus_eta.push_back(etaDirectory.
make<TH1D>(ap.c_str(), bp.c_str(), 50, 0., 200.));
182 ap =
"zmutrackplus_etaRange" + std::to_string(
i);
183 bp =
"zmutrack plus mass eta Range " + std::to_string(range0) +
" to " + std::to_string(range1);
184 cout << ap <<
" " << bp << endl;
186 ap =
"zmumuNotIsoplus_etaRange" + std::to_string(
i);
187 bp =
"zmumuNotIso plus mass eta Range " + std::to_string(range0) +
" to " + std::to_string(range1);
188 cout << ap <<
" " << bp << endl;
191 cout <<
" bin eta minus " <<
i << endl;
193 am =
"zmumu2HLTminus_etaRange" + std::to_string(
i);
194 bm =
"zmumu2HLT minus mass eta Range " + std::to_string(range0) +
" to " + std::to_string(range1);
195 cout << am <<
" " << bm << endl;
197 am =
"zmumu1HLTminus_etaRange" + std::to_string(
i);
198 bm =
"zmumu1HLT minus mass eta Range " + std::to_string(range0) +
" to " + std::to_string(range1);
199 cout << am <<
" " << bm << endl;
201 am =
"zmustaminus_etaRange" + std::to_string(
i);
202 bm =
"zmusta minus mass eta Range " + std::to_string(range0) +
" to " + std::to_string(range1);
203 cout << am <<
" " << bm << endl;
205 am =
"zmutrackminus_etaRange" + std::to_string(
i);
206 bm =
"zmutrack minus mass eta Range " + std::to_string(range0) +
" to " + std::to_string(range1);
207 cout << am <<
" " << bm << endl;
209 am =
"zmumuNotIsominus_etaRange" + std::to_string(
i);
210 bm =
"zmumuNotIso minus mass eta Range " + std::to_string(range0) +
" to " + std::to_string(range1);
211 cout << am <<
" " << bm << endl;
216 cout <<
"pt istograms creation " << endl;
218 for (
unsigned int i = 0;
i <
ptBins;
i++) {
222 cout <<
" bin pt plus " <<
i << endl;
224 ap1 =
"zmumu2HLTplus_ptRange" + std::to_string(
i);
225 bp1 =
"zmumu2HLT plus mass pt Range " + std::to_string(range0) +
" to " + std::to_string(range1);
226 cout << ap1 <<
" " << bp1 << endl;
228 ap1 =
"zmumu1HLTplus_ptRange" + std::to_string(
i);
229 bp1 =
"zmumu1HLT plus mass pt Range " + std::to_string(range0) +
" to " + std::to_string(range1);
230 cout << ap1 <<
" " << bp1 << endl;
232 ap1 =
"zmustaplus_ptRange" + std::to_string(
i);
233 bp1 =
"zmusta plus mass pt Range " + std::to_string(range0) +
" to " + std::to_string(range1);
234 cout << ap1 <<
" " << bp1 << endl;
235 hmustaplus_pt.push_back(ptDirectory.
make<TH1D>(ap1.c_str(), bp1.c_str(), 50, 0., 200.));
236 ap1 =
"zmutrackplus_ptRange" + std::to_string(
i);
237 bp1 =
"zmutrack plus mass pt Range " + std::to_string(range0) +
" to " + std::to_string(range1);
238 cout << ap1 <<
" " << bp1 << endl;
239 hmutrackplus_pt.push_back(ptDirectory.
make<TH1D>(ap1.c_str(), bp1.c_str(), 100, 0., 200.));
240 ap1 =
"zmumuNotIsoplus_ptRange" + std::to_string(
i);
241 bp1 =
"zmumuNotIso plus mass pt Range " + std::to_string(range0) +
" to " + std::to_string(range1);
242 cout << ap1 <<
" " << bp1 << endl;
245 cout <<
" bin pt minus " <<
i << endl;
247 am1 =
"zmumu2HLTminus_ptRange" + std::to_string(
i);
248 bm1 =
"zmumu2HLT minus mass pt Range " + std::to_string(range0) +
" to " + std::to_string(range1);
249 cout << am1 <<
" " << bm1 << endl;
251 am1 =
"zmumu1HLTminus_ptRange" + std::to_string(
i);
252 bm1 =
"zmumu1HLT minus mass pt Range " + std::to_string(range0) +
" to " + std::to_string(range1);
253 cout << am1 <<
" " << bm1 << endl;
255 am1 =
"zmustaminus_ptRange" + std::to_string(
i);
256 bm1 =
"zmusta minus mass pt Range " + std::to_string(range0) +
" to " + std::to_string(range1);
257 cout << am1 <<
" " << bm1 << endl;
258 hmustaminus_pt.push_back(ptDirectory.
make<TH1D>(am1.c_str(), bm1.c_str(), 50, 0., 200.));
259 am1 =
"zmutrackminus_ptRange" + std::to_string(
i);
260 bm1 =
"zmutrack minus mass pt Range " + std::to_string(range0) +
" to " + std::to_string(range1);
261 cout << am1 <<
" " << bm1 << endl;
263 am1 =
"zmumuNotIsominus_ptRange" + std::to_string(
i);
264 bm1 =
"zmumuNotIso minus mass pt Range " + std::to_string(range0) +
" to " + std::to_string(range1);
265 cout << am1 <<
" " << bm1 << endl;
T getParameter(std::string const &) const
vector< TH1D * > hmumuNotIsoplus_pt
T getUntrackedParameter(std::string const &, T const &) const
vector< TH1D * > hmumu1HLTplus_eta
vector< TH1D * > hmutrackplus_pt
EDGetTokenT< GenParticleMatch > zMuMuMatchMapToken_
EDGetTokenT< CandidateView > muonsToken_
vector< TH1D * > hmumu2HLTplus_eta
vector< TH1D * > hmumu2HLTminus_eta
EDGetTokenT< GenParticleMatch > zMuStandAloneMatchMapToken_
EDGetTokenT< GenParticleCollection > genParticlesToken_
vector< TH1D * > hmumuNotIsominus_eta
TH1D * h_zmmNotIsoplus_mass
vector< TH1D * > hmumuNotIsominus_pt
T * make(const Args &...args) const
make new ROOT object
TH1D * h_zmm1HLTplus_mass
EDGetTokenT< CandidateView > tracksToken_
TH1D * h_zmm1HLTminus_mass
EDGetTokenT< GenParticleMatch > zMuTrackMatchMapToken_
EDGetTokenT< CandidateView > zMuTrackToken_
vector< TH1D * > hmutrackminus_eta
vector< TH1D * > hmumu1HLTplus_pt
vector< TH1D * > hmustaplus_eta
TH1D * h_zmmNotIsominus_mass
vector< TH1D * > hmustaplus_pt
vector< TH1D * > hmustaminus_eta
vector< TH1D * > hmumu2HLTplus_pt
T * make(const Args &...args) const
make new ROOT object
vector< TH1D * > hmutrackplus_eta
EDGetTokenT< CandidateView > zMuStandAloneToken_
vector< TH1D * > hmustaminus_pt
vector< TH1D * > hmumu1HLTminus_eta
TFileDirectory mkdir(const std::string &dir, const std::string &descr="")
create a new subdirectory
vector< TH1D * > hmumu1HLTminus_pt
int nGlobalMuonsMatched_passed
EDGetTokenT< VertexCollection > primaryVerticesToken_
vector< TH1D * > hmumu2HLTminus_pt
vector< TH1D * > hmumuNotIsoplus_eta
EDGetTokenT< CandidateView > zMuMuToken_
vector< TH1D * > hmutrackminus_pt