194 const int ntptarray = 23;
195 const int njptarray = 14;
196 const int netaarray = 19;
197 Double_t jetetabins[netaarray] = {-2.5,-2.0,-1.75,-1.5,-1.25,-1.0,-0.75,-0.5,-0.25,0.0,0.25,0.5,0.75,1.0,1.25,1.5,1.75,2.0,2.5};
198 Double_t jetptxbins[njptarray] = {0.,10.,20.,30., 40., 50., 60., 70., 80, 90., 100., 120., 140., 230.};
199 Double_t jetpttbins[ntptarray] = {0.,1.,2.,4.,6.,8.,10.,15.,20.,25.,30., 35.,40., 45., 50., 60., 70., 80, 90., 100., 120., 140., 230.};
205 histoid =
"jet_pt_"+
flavor; histotitle =
flavor+
" jet p_{T} [GeV/c]";
206 histocontainer_[histoid]=fs->
make<TH1D>(histoid.c_str(),histotitle.c_str(),njptarray-1,jetptxbins);
207 histoid =
"jet_pt_uncorr_"+
flavor; histotitle =
flavor+
" jet uncorrected p_{T} [GeV/c]";
208 histocontainer_[histoid]=fs->
make<TH1D>(histoid.c_str(),histotitle.c_str(),njptarray-1,jetptxbins);
209 histoid =
"jet_eta_"+
flavor;
211 histoid =
"jet_eta_"+
flavor+
"030";
213 histoid =
"jet_eta_"+
flavor+
"3050";
215 histoid =
"jet_eta_"+
flavor+
"50";
217 histoid =
"jet_pt_"+
flavor+
"_tagged";
218 histocontainer_[
"jet_pt_"+
flavor+
"_tagged"]=fs->
make<TH1D>(histoid.c_str(),
"jet_tagged p_{T} [GeV/c]",njptarray-1,jetptxbins);
219 histoid =
"jet_pt_uncorr_"+
flavor+
"_tagged";
220 histocontainer_[
"jet_pt_uncorr_"+
flavor+
"_tagged"]=fs->
make<TH1D>(histoid.c_str(),
"jet_tagged p_{T} [GeV/c]",njptarray-1,jetptxbins);
221 histoid =
"jet_eta_"+
flavor+
"_tagged";
225 histoid =
"jet_pt_"+
flavor+
"_taggedmincut";
226 histocontainer_[
"jet_pt_"+
flavor+
"_taggedmincut"]=fs->
make<TH1D>(histoid.c_str(),
"jet_tagged p_{T} [GeV/c]",njptarray-1,jetptxbins);
227 histoid =
"jet_pt_uncorr_"+
flavor+
"_taggedmincut";
228 histocontainer_[
"jet_pt_uncorr_"+
flavor+
"_taggedmincut"]=fs->
make<TH1D>(histoid.c_str(),
"jet_tagged p_{T} [GeV/c]",njptarray-1,jetptxbins);
229 histoid =
"jet_eta_"+
flavor+
"_taggedmincut";
230 histocontainer_[
"jet_eta_"+
flavor+
"_taggedmincut"]=fs->
make<TH1D>(histoid.c_str(),
"jet_tagged #eta",netaarray-1,jetetabins);
233 histoid =
"tracks_in_jet_"+
flavor; histotitle =
"traks per jet "+
flavor;
235 histoid =
"tracks_in_jet_"+
flavor+
"030"; histotitle =
"traks per jet "+
flavor+
"pt_{T} < 30[GeV/c]";
237 histoid =
"tracks_in_jet_"+
flavor+
"3050"; histotitle =
"traks per jet "+
flavor+
"30 < pt_{T} < 50[GeV/c]";
239 histoid =
"tracks_in_jet_"+
flavor+
"50"; histotitle =
"traks per jet "+
flavor+
"pt_{T} > 50[GeV/c]";
243 histoid=
"pt_tracks_in_jet_"+
flavor; histotitle =
"track p_{T} [GeV/c] "+
flavor+
" jets";
244 histocontainer_[histoid]=fs->
make<TH1D>(histoid.c_str(),histotitle.c_str(),ntptarray-1,jetpttbins);
245 histoid=
"pt_tracks_in_jet_"+
flavor+
"030"; histotitle =
"track p_{T} [GeV/c] "+
flavor+
" jets";
246 histocontainer_[histoid]=fs->
make<TH1D>(histoid.c_str(),histotitle.c_str(),ntptarray-1,jetpttbins);
247 histoid=
"pt_tracks_in_jet_"+
flavor+
"3050"; histotitle =
"track p_{T} [GeV/c] "+
flavor+
" jets";
248 histocontainer_[histoid]=fs->
make<TH1D>(histoid.c_str(),histotitle.c_str(),ntptarray-1,jetpttbins);
249 histoid=
"pt_tracks_in_jet_"+
flavor+
"50"; histotitle =
"track p_{T} [GeV/c] "+
flavor+
" jets";
250 histocontainer_[histoid]=fs->
make<TH1D>(histoid.c_str(),histotitle.c_str(),ntptarray-1,jetpttbins);
254 histoid=
"pt_tracks_in_jet_"+
flavor+
"_center"; histotitle =
"track p_{T} [GeV/c] "+
flavor+
" jets";
255 histocontainer_[histoid]=fs->
make<TH1D>(histoid.c_str(),histotitle.c_str(),ntptarray-1,jetpttbins);
256 histoid=
"pt_tracks_in_jet_"+
flavor+
"_center030"; histotitle =
"track p_{T} [GeV/c] "+
flavor+
" jets";
257 histocontainer_[histoid]=fs->
make<TH1D>(histoid.c_str(),histotitle.c_str(),ntptarray-1,jetpttbins);
258 histoid=
"pt_tracks_in_jet_"+
flavor+
"_center3050"; histotitle =
"track p_{T} [GeV/c] "+
flavor+
" jets";
259 histocontainer_[histoid]=fs->
make<TH1D>(histoid.c_str(),histotitle.c_str(),ntptarray-1,jetpttbins);
260 histoid=
"pt_tracks_in_jet_"+
flavor+
"_center50"; histotitle =
"track p_{T} [GeV/c] "+
flavor+
" jets";
261 histocontainer_[histoid]=fs->
make<TH1D>(histoid.c_str(),histotitle.c_str(),ntptarray-1,jetpttbins);
262 histoid=
"pt_tracks_in_jet_"+
flavor+
"_sides"; histotitle =
"track p_{T} [GeV/c] "+
flavor+
" jets";
263 histocontainer_[histoid]=fs->
make<TH1D>(histoid.c_str(),histotitle.c_str(),ntptarray-1,jetpttbins);
264 histoid=
"pt_tracks_in_jet_"+
flavor+
"_sides030"; histotitle =
"track p_{T} [GeV/c] "+
flavor+
" jets";
265 histocontainer_[histoid]=fs->
make<TH1D>(histoid.c_str(),histotitle.c_str(),ntptarray-1,jetpttbins);
266 histoid=
"pt_tracks_in_jet_"+
flavor+
"_sides3050"; histotitle =
"track p_{T} [GeV/c] "+
flavor+
" jets";
267 histocontainer_[histoid]=fs->
make<TH1D>(histoid.c_str(),histotitle.c_str(),ntptarray-1,jetpttbins);
268 histoid=
"pt_tracks_in_jet_"+
flavor+
"_sides50"; histotitle =
"track p_{T} [GeV/c] "+
flavor+
" jets";
269 histocontainer_[histoid]=fs->
make<TH1D>(histoid.c_str(),histotitle.c_str(),ntptarray-1,jetpttbins);
272 std::map<int,std::string> tagbl;
273 tagbl[0]=
"0";tagbl[1]=
"1";tagbl[2]=
"2";
274 for (
size_t i=0;
i < tagbl.size();
i++){
275 histoid =
"jet_pt_"+
flavor+
"_taggable"+tagbl[
i]; histotitle =
flavor+
" jet_taggable p_{T} [GeV/c]";
276 histocontainer_[
"jet_pt_"+
flavor+
"_taggable"+tagbl[
i]]=fs->
make<TH1D>(histoid.c_str(),histotitle.c_str(),njptarray-1,jetptxbins);
277 histoid =
"jet_pt_uncorr_"+
flavor+
"_taggable"+tagbl[
i]; histotitle =
flavor+
" jet_taggable p_{T} uncorr [GeV/c]";
278 histocontainer_[
"jet_pt_uncorr_"+
flavor+
"_taggable"+tagbl[
i]]=fs->
make<TH1D>(histoid.c_str(),histotitle.c_str(),njptarray-1,jetptxbins);
279 histoid =
"jet_eta_"+
flavor+
"_taggable"+tagbl[
i]; histotitle =
flavor+
" jet_taggable #eta";
280 histocontainer_[
"jet_eta_"+
flavor+
"_taggable"+tagbl[
i]]=fs->
make<TH1D>(histoid.c_str(),histotitle.c_str(),netaarray-1,jetetabins);
281 histoid =
"jet_eta_"+
flavor+
"_taggable"+tagbl[
i]+
"_030";
282 histocontainer_[
"jet_eta_"+
flavor+
"_taggable"+tagbl[
i]+
"_030"]=fs->
make<TH1D>(histoid.c_str(),
"jet #eta",netaarray-1,jetetabins);
283 histoid =
"jet_eta_"+
flavor+
"_taggable"+tagbl[
i]+
"_3050";
284 histocontainer_[
"jet_eta_"+
flavor+
"_taggable"+tagbl[
i]+
"_3050"]=fs->
make<TH1D>(histoid.c_str(),
"jet #eta",netaarray-1,jetetabins);
285 histoid =
"jet_eta_"+
flavor+
"_taggable"+tagbl[
i]+
"_50";
286 histocontainer_[
"jet_eta_"+
flavor+
"_taggable"+tagbl[
i]+
"_50"]=fs->
make<TH1D>(histoid.c_str(),
"jet #eta",netaarray-1,jetetabins);
289 histoid =
"jet_scatter_pt_"+
flavor;
290 h2_[
"jet_scatter_pt_"+
flavor]=fs->
make<TH2D>(histoid.c_str(),
"jet p_{T} uncorr(y) vs corr(x) [GeV/c]",
291 njptarray-1,jetptxbins,njptarray-1,jetptxbins);
std::map< std::string, TH1D * > histocontainer_
T * make(const Args &...args) const
make new ROOT object
std::map< std::string, TH2D * > h2_