39 if (!
tfile.isAvailable())
41 <<
"please add it to config file";
42 iniE =
tfile->make<TH1D>(
"iniE",
"Incident Energy (GeV)", 4000, 0., em1);
43 iEta =
tfile->make<TH1D>(
"iEta",
"Eta at incidence ", 300, 0., 3.);
44 iPhi =
tfile->make<TH1D>(
"iPhi",
"Phi at incidence ", 300, -1., 1.);
45 edepS =
tfile->make<TH1D>(
"edepS",
"Energy deposit == Total", 4000, 0., em1);
46 edecS =
tfile->make<TH1D>(
"edecS",
"Energy deposit == ECal ", 300, -2., 28.);
47 edhcS =
tfile->make<TH1D>(
"edhcS",
"Energy deposit == HCal ", 4000, 0., em2);
48 edepN =
tfile->make<TH1D>(
"edepN",
"Etot/Ebeam ", 200, -2.5, 2.5);
49 edecN =
tfile->make<TH1D>(
"edecN",
"Eecal/Ebeam ", 200, -2.5, 2.5);
50 edhcN =
tfile->make<TH1D>(
"edhcN",
"Ehcal/Ebeam ", 200, -2.5, 2.5);
51 emhcN =
tfile->make<TH1D>(
"emhcN",
"Ehcal/Ebeam MIP in Ecal", 200, -2.5, 2.5);
52 edehS =
tfile->make<TH2D>(
"edehS",
"Hcal vs Ecal", 100, 0., em1, 100, 0., em2);
55 tree_ =
tfile->make<TTree>(
"TB06Sim",
"TB06Sim");
94 edm::LogInfo(
"HcalTBSim") <<
"HcalTB06Histo:::fillEdep: Simulated Total " << etots <<
" ECal " << eecals <<
" HCal " 102 if (eecals <=
mip_) {
105 edehS->Fill(eecals, ehcals);
T getParameter(std::string const &) const
void fillEdep(double etots, double eecals, double ehcals)
T getUntrackedParameter(std::string const &, T const &) const
void fillTree(std::vector< double > &ecalo, std::vector< double > &etrig)
Log< level::Info, false > LogInfo
void fillPrimary(double energy, double eta, double phi)
HcalTB06Histo(const edm::ParameterSet &ps)