Inizialize parameters for histo binning.
151 int etaNBinsPFCand =
v->getParameter<
int>(
"etaNBins");
152 double etaMinPFCand =
v->getParameter<
double>(
"etaMin");
153 double etaMaxPFCand =
v->getParameter<
double>(
"etaMax");
154 int phiNBinsPFCand =
v->getParameter<
int>(
"phiNBins");
155 double phiMinPFCand =
v->getParameter<
double>(
"phiMin");
156 double phiMaxPFCand =
v->getParameter<
double>(
"phiMax");
157 int nMinPFCand =
v->getParameter<
int>(
"nMin");
158 int nMaxPFCand =
v->getParameter<
int>(
"nMax");
159 int nbinsPFCand =
v->getParameter<
double>(
"nbins");
175 if(etaMinPFCand*etaMaxPFCand<0){
179 const int nbins_eta_endcap=2*(etaNBinsPFCand+1);
180 double eta_limits_endcap[nbins_eta_endcap];
181 for(
int i=0;
i<nbins_eta_endcap;
i++){
182 if(
i<(etaNBinsPFCand+1)){
183 eta_limits_endcap[
i]=etaMinPFCand+
i*(etaMaxPFCand-etaMinPFCand)/(
double)etaNBinsPFCand;
185 eta_limits_endcap[
i]= -etaMaxPFCand +(
i- (etaNBinsPFCand+1) )*(etaMaxPFCand-etaMinPFCand)/(
double)etaNBinsPFCand;
190 TH2F* hist_temp_pt =
new TH2F((
ptPFCand_nameRECO_[
ptPFCand_nameRECO_.size()-1]).c_str(),
"pt",nbins_eta_endcap-1, eta_limits_endcap, phiNBinsPFCand, phiMinPFCand, phiMaxPFCand);
214 const int nbins_eta=16;
216 double eta_limits[nbins_eta]={-2.650,-2.500,-2.322,-2.172,-2.043,-1.930,-1.830,-1.740,1.740,1.830,1.930,2.043,2.172,2.3122,2.500,2.650};
218 TH2F* hist_temp_HCAL =
new TH2F(
"IsoPfChHad_HCAL_profile_endcap",
"IsolatedPFChHa HCAL occupancy (outer endcap)",nbins_eta-1,eta_limits, 36, -
M_PI,
M_PI);
219 TH2F* hist_tempPt_HCAL=(TH2F*)hist_temp_HCAL->Clone(
"Isolated PFCHHadron HadPt (outer endcap)");
239 const int nbins_eta_hcal_total=28;
240 double eta_limits_hcal_total[nbins_eta_hcal_total]={-2.650,-2.322,-2.043,-1.830,-1.653,-1.479,-1.305,-1.131,-0.957,-0.783,-0.609,-0.435,-0.261,-0.087,
241 0.087, 0.261, 0.435, 0.609, 0.783, 0.957, 1.131, 1.305, 1.479, 1.653, 1.830, 2.043, 2.322, 2.650};
242 float eta_limits_hcal_total_f[nbins_eta_hcal_total];
243 float log_bin_spacing =
log(200.)/40.;
244 const int nbins_pt_total_hcal= 41;
245 double pt_limits_hcal[nbins_pt_total_hcal];
246 float pt_limits_hcal_f[nbins_pt_total_hcal];
247 for(
int i=0;
i<nbins_pt_total_hcal;
i++){
248 pt_limits_hcal[
i]=
exp(
i*log_bin_spacing);
249 pt_limits_hcal_f[
i]=
exp(
i*log_bin_spacing);
251 for(
int i=0;
i<nbins_eta_hcal_total;
i++){
252 eta_limits_hcal_total_f[
i]=(
float)eta_limits_hcal_total[
i];
254 m_HOverTrackP_trackPtVsEta= ibooker.
book2D(
"HOverTrackP_trackPtVsEta",
"HOverTrackP_trackPtVsEta",nbins_pt_total_hcal-1,pt_limits_hcal_f,nbins_eta_hcal_total-1,eta_limits_hcal_total_f);
308 int etaNBinsPFCand =
v->getParameter<
int>(
"etaNBins");
309 double etaMinPFCand =
v->getParameter<
double>(
"etaMin");
310 double etaMaxPFCand =
v->getParameter<
double>(
"etaMax");
311 int phiNBinsPFCand =
v->getParameter<
int>(
"phiNBins");
312 double phiMinPFCand =
v->getParameter<
double>(
"phiMin");
313 double phiMaxPFCand =
v->getParameter<
double>(
"phiMax");
314 int nMinPFCand =
v->getParameter<
int>(
"nMin");
315 int nMaxPFCand =
v->getParameter<
int>(
"nMax");
316 int nbinsPFCand =
v->getParameter<
double>(
"nbins");
336 if(etaMinPFCand*etaMaxPFCand<0){
342 const int nbins_eta_endcap=2*(etaNBinsPFCand+1);
343 double eta_limits_endcap[nbins_eta_endcap];
344 for(
int i=0;
i<nbins_eta_endcap;
i++){
345 if(
i<(etaNBinsPFCand+1)){
346 eta_limits_endcap[
i]=etaMinPFCand+
i*(etaMaxPFCand-etaMinPFCand)/(
double)etaNBinsPFCand;
348 eta_limits_endcap[
i]= -etaMaxPFCand +(
i- (etaNBinsPFCand+1) )*(etaMaxPFCand-etaMinPFCand)/(
double)etaNBinsPFCand;
MonitorElement * m_HOverTrackPVsTrackP_EndCap
std::vector< int > typePFCand_
std::vector< MonitorElement * > occupancyPFCand_puppiNolepWeight_
MonitorElement * m_HOverTrackPVsTrackPt_Barrel
std::vector< double > etaMinPFCandRECO_
std::vector< std::string > ptPFCand_name_puppiNolepWeight_
MonitorElement * m_HOverTrackP_Barrel_hPt_20_50
MonitorElement * bookProfile(Args &&...args)
MonitorElement * m_HOverTrackPVsEta_hPt_20_50
std::vector< std::string > multiplicityPFCand_nameRECO_
MonitorElement * m_HOverTrackP_trackPtVsEta
MonitorElement * m_HOverTrackPVsEta_hPt_10_20
MonitorElement * m_HOverTrackP_Barrel_hPt_1_10
std::vector< std::string > occupancyPFCand_name_
MonitorElement * m_HOverTrackPVsTrackP_Barrel
std::vector< int > typePFCandRECO_
MonitorElement * m_HOverTrackP_EndCap_hPt_50
std::vector< int > countsPFCandRECO_
std::vector< MonitorElement * > multiplicityPFCandRECO_
MonitorElement * m_HOverTrackPVsEta_hPt_1_10
MonitorElement * mProfileIsoPFChHad_HadPtEndcap
MonitorElement * mProfileIsoPFChHad_HcalOccupancyEndcap
std::vector< MonitorElement * > ptPFCandRECO_
MonitorElement * mProfileIsoPFChHad_EMPtCentral
void setCurrentFolder(std::string const &fullpath)
std::vector< MonitorElement * > occupancyPFCandRECO_
std::vector< MonitorElement * > occupancyPFCand_
MonitorElement * book1D(Args &&...args)
std::vector< double > etaMaxPFCandRECO_
std::vector< std::string > multiplicityPFCand_name_
std::vector< std::string > occupancyPFCand_nameRECO_
std::vector< double > etaMinPFCand_
std::map< std::string, MonitorElement * > map_of_MEs
MonitorElement * m_HOverTrackPVsTrackPt_EndCap
std::vector< MonitorElement * > ptPFCand_puppiNolepWeight_
MonitorElement * mProfileIsoPFChHad_EcalOccupancyEndcap
std::vector< std::string > occupancyPFCand_name_puppiNolepWeight_
MonitorElement * m_HOverTrackP_Barrel_hPt_10_20
MonitorElement * book2D(Args &&...args)
std::vector< int > countsPFCand_
MonitorElement * m_HOverTrackP_EndCap_hPt_1_10
MonitorElement * mProfileIsoPFChHad_TrackPt
std::vector< std::string > ptPFCand_nameRECO_
MonitorElement * mProfileIsoPFChHad_HcalOccupancyCentral
MonitorElement * mProfileIsoPFChHad_TrackOccupancy
edm::InputTag mInputCollection_
MonitorElement * m_HOverTrackP_Barrel_hPt_50
MonitorElement * mProfileIsoPFChHad_EcalOccupancyCentral
std::vector< MonitorElement * > multiplicityPFCand_
MonitorElement * mProfileIsoPFChHad_EMPtEndcap
MonitorElement * m_HOverTrackPVsEta_hPt_50
std::vector< double > etaMaxPFCand_
MonitorElement * m_HOverTrackP_EndCap_hPt_10_20
MonitorElement * m_HOverTrackP_EndCap_hPt_20_50
std::vector< edm::ParameterSet > diagnosticsParameters_
std::vector< MonitorElement * > ptPFCand_
std::vector< std::string > ptPFCand_name_
MonitorElement * mProfileIsoPFChHad_HadPtCentral