Definition at line 70 of file CaloSimHitStudy.cc.
References edep_, edepC_, edepEM_, edepHad_, edepT_, edepTk_, eMIP_, eneInc_, etaInc_, etot_, etotg_, g4Label, edm::ParameterSet::getUntrackedParameter(), hit_, hitHigh, hitLab, hitLow, hitMu, hitTk_, mps_fire::i, HLT_2018_cff::InputTag, edm::Service< T >::isAvailable(), TFileService::kSharedResource, TFileService::make(), muonLab, Skims_PA_cff::name, phiInc_, ptInc_, storeRL_, AlCaHLTBitMon_QueryRunRegistry::string, testNumber_, compare::tfile, time_, timeAll_, overlapproblemtsosanalyzer_cfi::title, tkHighLab, tkLowLab, tmax_, tofTk_, tok_evt_, toks_calo_, toks_tkHigh_, toks_tkLow_, toks_track_, and L1TOccupancyClient_cfi::ymax.
90 tkHighLab[0] =
"TrackerHitsPixelBarrelHighTof";
91 tkHighLab[1] =
"TrackerHitsPixelEndcapHighTof";
96 tkLowLab[0] =
"TrackerHitsPixelBarrelLowTof";
97 tkLowLab[1] =
"TrackerHitsPixelEndcapLowTof";
98 tkLowLab[2] =
"TrackerHitsTECLowTof";
99 tkLowLab[3] =
"TrackerHitsTIBLowTof";
100 tkLowLab[4] =
"TrackerHitsTIDLowTof";
101 tkLowLab[5] =
"TrackerHitsTOBLowTof";
104 for (
unsigned i = 0;
i != 4;
i++)
107 for (
unsigned i = 0; i != 3; i++)
110 for (
unsigned i = 0; i != 6; i++) {
116 <<
hitLab[2] <<
", " <<
hitLab[3] <<
" MaxEnergy: " << maxEnergy_
117 <<
" Tmax: " <<
tmax_ <<
" MIP Cut: " <<
eMIP_;
123 <<
"please add it to config file";
125 sprintf(title,
"Incident PT (GeV)");
127 ptInc_->GetXaxis()->SetTitle(title);
128 ptInc_->GetYaxis()->SetTitle(
"Events");
129 sprintf(title,
"Incident Energy (GeV)");
131 eneInc_->GetXaxis()->SetTitle(title);
132 eneInc_->GetYaxis()->SetTitle(
"Events");
133 sprintf(title,
"Incident #eta");
135 etaInc_->GetXaxis()->SetTitle(title);
136 etaInc_->GetYaxis()->SetTitle(
"Events");
137 sprintf(title,
"Incident #phi");
139 phiInc_->GetXaxis()->SetTitle(title);
140 phiInc_->GetYaxis()->SetTitle(
"Events");
141 std::string dets[9] = {
"EB",
"EB(APD)",
"EB(ATJ)",
"EE",
"ES",
"HB",
"HE",
"HO",
"HF"};
142 for (
int i = 0; i < 9; i++) {
143 sprintf(name,
"Hit%d", i);
144 sprintf(title,
"Number of hits in %s", dets[i].c_str());
146 hit_[
i]->GetXaxis()->SetTitle(title);
147 hit_[
i]->GetYaxis()->SetTitle(
"Events");
148 sprintf(name,
"Time%d", i);
149 sprintf(title,
"Time of the hit (ns) in %s", dets[i].c_str());
151 time_[
i]->GetXaxis()->SetTitle(title);
152 time_[
i]->GetYaxis()->SetTitle(
"Hits");
153 sprintf(name,
"TimeAll%d", i);
154 sprintf(title,
"Hit time (ns) in %s with no check on Track ID", dets[i].c_str());
156 timeAll_[
i]->GetXaxis()->SetTitle(title);
157 timeAll_[
i]->GetYaxis()->SetTitle(
"Hits");
158 double ymax = maxEnergy_;
159 if (i == 1 || i == 2 || i == 4)
161 else if (i > 4 && i < 8)
163 sprintf(name,
"Edep%d", i);
164 sprintf(title,
"Energy deposit (GeV) in %s", dets[i].c_str());
166 edep_[
i]->GetXaxis()->SetTitle(title);
167 edep_[
i]->GetYaxis()->SetTitle(
"Hits");
168 sprintf(name,
"EdepEM%d", i);
169 sprintf(title,
"Energy deposit (GeV) by EM particles in %s", dets[i].c_str());
171 edepEM_[
i]->GetXaxis()->SetTitle(title);
172 edepEM_[
i]->GetYaxis()->SetTitle(
"Hits");
173 sprintf(name,
"EdepHad%d", i);
174 sprintf(title,
"Energy deposit (GeV) by hadrons in %s", dets[i].c_str());
176 edepHad_[
i]->GetXaxis()->SetTitle(title);
177 edepHad_[
i]->GetYaxis()->SetTitle(
"Hits");
178 sprintf(name,
"Etot%d", i);
179 sprintf(title,
"Total energy deposit (GeV) in %s", dets[i].c_str());
181 etot_[
i]->GetXaxis()->SetTitle(title);
182 etot_[
i]->GetYaxis()->SetTitle(
"Events");
183 sprintf(name,
"EtotG%d", i);
184 sprintf(title,
"Total energy deposit (GeV) in %s (t < 100 ns)", dets[i].c_str());
186 etotg_[
i]->GetXaxis()->SetTitle(title);
187 etotg_[
i]->GetYaxis()->SetTitle(
"Events");
198 for (
int i = 0; i < 9; i++) {
200 if (i == 0 || i == 3 || i == 6)
202 sprintf(name,
"EdepCal%d", i);
203 sprintf(title,
"Energy deposit in %s", detx[i].c_str());
205 edepC_[
i]->GetXaxis()->SetTitle(title);
206 edepC_[
i]->GetYaxis()->SetTitle(
"Events");
207 sprintf(name,
"EdepCalT%d", i);
208 sprintf(title,
"Energy deposit (t < %f ns) in %s",
tmax_, detx[i].c_str());
210 edepT_[
i]->GetXaxis()->SetTitle(title);
211 edepT_[
i]->GetYaxis()->SetTitle(
"Events");
213 hitLow = tfile->
make<TH1F>(
"HitLow",
"Number of hits in Track (Low)", 1000, 0, 10000.);
214 hitLow->GetXaxis()->SetTitle(
"Number of hits in Track (Low)");
215 hitLow->GetYaxis()->SetTitle(
"Events");
216 hitHigh = tfile->
make<TH1F>(
"HitHigh",
"Number of hits in Track (High)", 1000, 0, 10000.);
217 hitHigh->GetXaxis()->SetTitle(
"Number of hits in Track (High)");
218 hitHigh->GetYaxis()->SetTitle(
"Events");
219 hitMu = tfile->
make<TH1F>(
"HitMu",
"Number of hits in Track (Muon)", 1000, 0, 5000.);
220 hitMu->GetXaxis()->SetTitle(
"Number of hits in Muon");
221 hitMu->GetYaxis()->SetTitle(
"Events");
223 "Pixel Endcap (High)",
228 "Pixel Barrel (Low)",
229 "Pixel Endcap (Low)",
237 for (
int i = 0; i < 15; i++) {
238 sprintf(name,
"HitTk%d", i);
239 sprintf(title,
"Number of hits in %s", dett[i].c_str());
241 hitTk_[
i]->GetXaxis()->SetTitle(title);
242 hitTk_[
i]->GetYaxis()->SetTitle(
"Events");
243 sprintf(name,
"TimeTk%d", i);
244 sprintf(title,
"Time of the hit (ns) in %s", dett[i].c_str());
246 tofTk_[
i]->GetXaxis()->SetTitle(title);
247 tofTk_[
i]->GetYaxis()->SetTitle(
"Hits");
248 sprintf(name,
"EdepTk%d", i);
249 sprintf(title,
"Energy deposit (GeV) in %s", dett[i].c_str());
251 edepTk_[
i]->GetXaxis()->SetTitle(title);
252 edepTk_[
i]->GetYaxis()->SetTitle(
"Hits");
static const std::string kSharedResource
edm::EDGetTokenT< edm::PSimHitContainer > toks_tkHigh_[6]
T getUntrackedParameter(std::string const &, T const &) const
T * make(const Args &...args) const
make new ROOT object
edm::EDGetTokenT< edm::HepMCProduct > tok_evt_
edm::EDGetTokenT< edm::PSimHitContainer > toks_tkLow_[6]
edm::EDGetTokenT< edm::PSimHitContainer > toks_track_[3]
edm::EDGetTokenT< edm::PCaloHitContainer > toks_calo_[4]