CMS 3D CMS Logo
List of all members | Public Member Functions | Static Public Member Functions | Private Member Functions | Private Attributes | Static Private Attributes
HcalHBHEMuonHighEtaAnalyzer Class Reference
Inheritance diagram for HcalHBHEMuonHighEtaAnalyzer:
edm::one::EDAnalyzer< edm::one::WatchRuns, edm::one::SharedResources > edm::one::EDAnalyzerBase edm::EDConsumerBase

Public Member Functions

 HcalHBHEMuonHighEtaAnalyzer (const edm::ParameterSet &)
 
- Public Member Functions inherited from edm::one::EDAnalyzer< edm::one::WatchRuns, edm::one::SharedResources >
 EDAnalyzer ()=default
 
 EDAnalyzer (const EDAnalyzer &)=delete
 
SerialTaskQueue * globalLuminosityBlocksQueue () final
 
SerialTaskQueue * globalRunsQueue () final
 
const EDAnalyzer & operator= (const EDAnalyzer &)=delete
 
bool wantsGlobalLuminosityBlocks () const final
 
bool wantsGlobalRuns () const final
 
bool wantsInputProcessBlocks () const final
 
bool wantsProcessBlocks () const final
 
- Public Member Functions inherited from edm::one::EDAnalyzerBase
void callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
 
 EDAnalyzerBase ()
 
ModuleDescription const & moduleDescription () const
 
bool wantsStreamLuminosityBlocks () const
 
bool wantsStreamRuns () const
 
 ~EDAnalyzerBase () override
 
- Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo > consumesInfo () const
 
void convertCurrentProcessAlias (std::string const &processName)
 Convert "@currentProcess" in InputTag process names to the actual current process name. More...
 
 EDConsumerBase ()
 
 EDConsumerBase (EDConsumerBase const &)=delete
 
 EDConsumerBase (EDConsumerBase &&)=default
 
ESResolverIndex const * esGetTokenIndices (edm::Transition iTrans) const
 
std::vector< ESResolverIndex > const & esGetTokenIndicesVector (edm::Transition iTrans) const
 
std::vector< ESRecordIndex > const & esGetTokenRecordIndicesVector (edm::Transition iTrans) const
 
ProductResolverIndexAndSkipBit indexFrom (EDGetToken, BranchType, TypeID const &) const
 
void itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
 
void itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
 
std::vector< ProductResolverIndexAndSkipBit > const & itemsToGetFrom (BranchType iType) const
 
void labelsForToken (EDGetToken iToken, Labels &oLabels) const
 
void modulesWhoseProductsAreConsumed (std::array< std::vector< ModuleDescription const *> *, NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const *> const &labelsToDesc, std::string const &processName) const
 
EDConsumerBase const & operator= (EDConsumerBase const &)=delete
 
EDConsumerBase & operator= (EDConsumerBase &&)=default
 
bool registeredToConsume (ProductResolverIndex, bool, BranchType) const
 
void selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)
 
ProductResolverIndexAndSkipBit uncheckedIndexFrom (EDGetToken) const
 
void updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
 
void updateLookup (eventsetup::ESRecordsToProductResolverIndices const &)
 
virtual ~EDConsumerBase () noexcept(false)
 

Static Public Member Functions

static void fillDescriptions (edm::ConfigurationDescriptions &descriptions)
 
- Static Public Member Functions inherited from edm::one::EDAnalyzerBase
static const std::string & baseType ()
 
static void fillDescriptions (ConfigurationDescriptions &descriptions)
 
static void prevalidate (ConfigurationDescriptions &descriptions)
 

Private Member Functions

double activeLength (const DetId &)
 
void analyze (edm::Event const &, edm::EventSetup const &) override
 
bool analyzeHadron (edm::Event const &, math::XYZPoint &)
 
bool analyzeMuon (edm::Event const &, math::XYZPoint &)
 
bool analyzeTracks (const reco::Track *, math::XYZPoint &, int, std::vector< spr::propagatedTrackID > &, bool)
 
void beginJob () override
 
void beginRun (edm::Run const &, edm::EventSetup const &) override
 
void clearVectors ()
 
int depth16HE (int, int)
 
void endRun (edm::Run const &, edm::EventSetup const &) override
 
void fillTrackParameters (const reco::Track *, math::XYZPoint)
 
double gainFactor (const HcalDbService *, const HcalDetId &)
 
bool goodCell (const HcalDetId &, const reco::Track *, const CaloGeometry *, const MagneticField *)
 
bool isGoodVertex (const reco::Vertex &)
 
int matchId (const HcalDetId &, const HcalDetId &)
 
double respCorr (const DetId &)
 

Private Attributes

std::vector< HcalDDDRecConstants::HcalActiveLength > actHB
 
std::vector< HcalDDDRecConstants::HcalActiveLength > actHE
 
const bool analyzeMuon_
 
edm::Handle< EcalRecHitCollection > barrelRecHitsHandle_
 
const MagneticField * bField_
 
const CaloTopology * caloTopology_
 
const bool collapseDepth_
 
const HcalDbService * conditions_
 
std::map< DetId, double > corrValue_
 
std::vector< double > ecal3x3Energy_
 
std::vector< unsigned int > ecalDetId_
 
std::vector< double > ecalEnergy_
 
std::vector< unsigned int > ehcalDetId_
 
std::vector< double > emaxNearP_
 
const double emaxNearPThr_
 
edm::Handle< EcalRecHitCollection > endcapRecHitsHandle_
 
std::vector< double > energyMuon_
 
std::vector< double > etaGlob_
 
const double etaMin_
 
unsigned int eventNumber_
 
const std::string fileInCorr_
 
const CaloGeometry * geo_
 
const bool getCharge_
 
unsigned int goodVertex_
 
edm::Handle< HBHERecHitCollection > hbhe_
 
std::vector< double > hcal1x1Energy_
 
std::vector< int > hcal_ieta_
 
std::vector< int > hcal_iphi_
 
std::vector< double > hcalActiveLength_
 
std::vector< double > hcalActiveLengthHot_
 
std::vector< double > hcalDepthActiveLength_ [depthMax_]
 
std::vector< double > hcalDepthActiveLengthHot_ [depthMax_]
 
std::vector< double > hcalDepthChargeHot_ [depthMax_]
 
std::vector< double > hcalDepthChargeHotBG_ [depthMax_]
 
std::vector< double > hcalDepthEnergy_ [depthMax_]
 
std::vector< double > hcalDepthEnergyCorr_ [depthMax_]
 
std::vector< double > hcalDepthEnergyHot_ [depthMax_]
 
std::vector< double > hcalDepthEnergyHotCorr_ [depthMax_]
 
std::vector< bool > hcalDepthMatch_ [depthMax_]
 
std::vector< bool > hcalDepthMatchHot_ [depthMax_]
 
std::vector< unsigned int > hcalDetId_
 
std::vector< double > hcalEnergy_
 
std::vector< bool > hcalHot_
 
const HcalDDDRecConstants * hdc_
 
std::vector< double > hoEnergy_
 
const bool ignoreHECorr_
 
const bool isItPlan1_
 
const bool isItPreRecHit_
 
std::vector< double > isolationR03_
 
std::vector< double > isolationR04_
 
int kount_
 
const edm::InputTag labelEBRecHit_
 
const edm::InputTag labelEERecHit_
 
const std::string labelGenTrack_
 
const edm::InputTag labelHBHERecHit_
 
const std::string labelMuon_
 
const std::string labelVtx_
 
std::vector< bool > matchedId_
 
const int maxDepth_
 
std::vector< bool > mediumMuon_
 
bool mergedDepth_
 
std::vector< double > phiGlob_
 
std::vector< double > pMuon_
 
std::vector< double > ptGlob_
 
HcalRespCorrs * respCorrs_
 
unsigned int runNumber_
 
spr::trackSelectionParameters selectionParameter_
 
const EcalSeverityLevelAlgo * sevlv_
 
const EcalChannelStatus * theEcalChStatus_
 
const HcalTopology * theHBHETopology_
 
const std::string theTrackQuality_
 
const edm::ESGetToken< EcalChannelStatus, EcalChannelStatusRcd > tok_chan_
 
const edm::ESGetToken< HcalDbService, HcalDbRecord > tok_dbservice_
 
const edm::ESGetToken< HcalDDDRecConstants, HcalRecNumberingRecord > tok_ddrec_
 
const edm::EDGetTokenT< EcalRecHitCollection > tok_EB_
 
const edm::EDGetTokenT< EcalRecHitCollection > tok_EE_
 
const edm::EDGetTokenT< reco::TrackCollection > tok_genTrack_
 
const edm::ESGetToken< CaloGeometry, CaloGeometryRecord > tok_geom_
 
const edm::EDGetTokenT< HBHERecHitCollection > tok_HBHE_
 
const edm::ESGetToken< HcalTopology, HcalRecNumberingRecord > tok_htopo_
 
const edm::ESGetToken< MagneticField, IdealMagneticFieldRecord > tok_magField_
 
const edm::EDGetTokenT< reco::MuonCollection > tok_Muon_
 
const edm::ESGetToken< HcalRespCorrs, HcalRespCorrsRcd > tok_respcorr_
 
const edm::ESGetToken< EcalSeverityLevelAlgo, EcalSeverityLevelAlgoRcd > tok_sevlv_
 
const edm::ESGetToken< CaloTopology, CaloTopologyRecord > tok_topo_
 
const edm::EDGetTokenT< reco::VertexCollection > tok_Vtx_
 
std::vector< double > trackDz_
 
std::vector< int > trackLayerCrossed_
 
std::vector< int > trackMissedInnerHits_
 
std::vector< int > trackMissedOuterHits_
 
std::vector< int > trackOuterHit_
 
TTree * tree_
 
const bool unCorrect_
 
bool useMyCorr_
 
const int useRaw_
 
const int verbosity_
 
const bool writeRespCorr_
 

Static Private Attributes

static const int depthMax_ = 7
 

Additional Inherited Members

- Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase ModuleType
 
- Public Types inherited from edm::EDConsumerBase
typedef ProductLabels Labels
 
- Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > consumes (edm::InputTag const &tag)
 
template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B > consumes (edm::InputTag tag) noexcept
 
EDGetToken consumes (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken consumes (TypeToGet const &id, edm::InputTag const &tag)
 
ConsumesCollector consumesCollector ()
 Use a ConsumesCollector to gather consumes information from helper functions. More...
 
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto esConsumes ()
 
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto esConsumes (ESInputTag const &tag)
 
template<Transition Tr = Transition::Event>
constexpr auto esConsumes ()
 
template<Transition Tr = Transition::Event>
auto esConsumes (ESInputTag tag)
 
template<Transition Tr = Transition::Event>
ESGetTokenGeneric esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
 Used with EventSetupRecord::doGet. More...
 
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > mayConsume (edm::InputTag const &tag)
 
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 
void resetItemsToGetFrom (BranchType iType)
 

Detailed Description

Definition at line 76 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Constructor & Destructor Documentation

◆ HcalHBHEMuonHighEtaAnalyzer()

HcalHBHEMuonHighEtaAnalyzer::HcalHBHEMuonHighEtaAnalyzer ( const edm::ParameterSet &  iConfig)
explicit

Definition at line 178 of file HcalHBHEMuonHighEtaAnalyzer.cc.

References collapseDepth_, corrValue_, fileInCorr_, getCharge_, edm::ParameterSet::getUntrackedParameter(), l1ctLayer2EG_cff::id, ignoreHECorr_, timingPdfMaker::infile, isItPlan1_, isItPreRecHit_, kount_, TFileService::kSharedResource, labelEBRecHit_, labelEERecHit_, labelGenTrack_, labelHBHERecHit_, labelMuon_, labelVtx_, spr::trackSelectionParameters::maxChi2, spr::trackSelectionParameters::maxDpOverP, spr::trackSelectionParameters::maxDxyPV, spr::trackSelectionParameters::maxDzPV, spr::trackSelectionParameters::maxInMiss, spr::trackSelectionParameters::maxOutMiss, mergedDepth_, spr::trackSelectionParameters::minLayerCrossed, spr::trackSelectionParameters::minOuterHit, spr::trackSelectionParameters::minPt, spr::trackSelectionParameters::minQuality, reco::TrackBase::qualityByName(), selectionParameter_, theTrackQuality_, StandaloneTrackMonitor_cfi::trackQuality, unCorrect_, useMyCorr_, and useRaw_.

179  : labelEBRecHit_(iConfig.getParameter<edm::InputTag>("labelEBRecHit")),
180  labelEERecHit_(iConfig.getParameter<edm::InputTag>("labelEERecHit")),
181  labelHBHERecHit_(iConfig.getParameter<edm::InputTag>("labelHBHERecHit")),
182  labelVtx_(iConfig.getParameter<std::string>("labelVertex")),
183  labelMuon_(iConfig.getParameter<std::string>("labelMuon")),
184  labelGenTrack_(iConfig.getParameter<std::string>("labelTrack")),
185  etaMin_(iConfig.getParameter<double>("etaMin")),
186  emaxNearPThr_(iConfig.getParameter<double>("emaxNearPThreshold")),
187  analyzeMuon_(iConfig.getParameter<bool>("analyzeMuon")),
188  unCorrect_(iConfig.getParameter<bool>("unCorrect")),
189  collapseDepth_(iConfig.getParameter<bool>("collapseDepth")),
190  isItPlan1_(iConfig.getParameter<bool>("isItPlan1")),
191  getCharge_(iConfig.getParameter<bool>("getCharge")),
192  useRaw_(iConfig.getParameter<int>("useRaw")),
193  verbosity_(iConfig.getParameter<int>("verbosity")),
194  theTrackQuality_(iConfig.getUntrackedParameter<std::string>("trackQuality")),
195  fileInCorr_(iConfig.getUntrackedParameter<std::string>("fileInCorr", "")),
196  ignoreHECorr_(iConfig.getUntrackedParameter<bool>("ignoreHECorr", false)),
197  isItPreRecHit_(iConfig.getUntrackedParameter<bool>("isItPreRecHit", false)),
198  writeRespCorr_(iConfig.getUntrackedParameter<bool>("writeRespCorr", false)),
199  maxDepth_(iConfig.getUntrackedParameter<int>("maxDepth", 7)),
200  tok_Vtx_(consumes<reco::VertexCollection>(labelVtx_)),
201  tok_EB_(consumes<EcalRecHitCollection>(labelEBRecHit_)),
202  tok_EE_(consumes<EcalRecHitCollection>(labelEERecHit_)),
203  tok_HBHE_(consumes<HBHERecHitCollection>(labelHBHERecHit_)),
204  tok_Muon_(consumes<reco::MuonCollection>(labelMuon_)),
205  tok_genTrack_(consumes<reco::TrackCollection>(labelGenTrack_)),
206  tok_ddrec_(esConsumes<HcalDDDRecConstants, HcalRecNumberingRecord, edm::Transition::BeginRun>()),
207  tok_htopo_(esConsumes<HcalTopology, HcalRecNumberingRecord, edm::Transition::BeginRun>()),
208  tok_respcorr_(esConsumes<HcalRespCorrs, HcalRespCorrsRcd, edm::Transition::BeginRun>()),
209  tok_geom_(esConsumes<CaloGeometry, CaloGeometryRecord, edm::Transition::BeginRun>()),
210  tok_magField_(esConsumes<MagneticField, IdealMagneticFieldRecord>()),
211  tok_chan_(esConsumes<EcalChannelStatus, EcalChannelStatusRcd>()),
212  tok_sevlv_(esConsumes<EcalSeverityLevelAlgo, EcalSeverityLevelAlgoRcd>()),
213  tok_topo_(esConsumes<CaloTopology, CaloTopologyRecord>()),
214  tok_dbservice_(esConsumes<HcalDbService, HcalDbRecord>()),
215  hdc_(nullptr),
216  theHBHETopology_(nullptr),
217  respCorrs_(nullptr),
218  tree_(nullptr) {
219  usesResource(TFileService::kSharedResource);
220  //now do what ever initialization is needed
221  kount_ = 0;
222 
223  reco::TrackBase::TrackQuality trackQuality = reco::TrackBase::qualityByName(theTrackQuality_);
224  selectionParameter_.minPt = iConfig.getUntrackedParameter<double>("minTrackPt");
225  selectionParameter_.minQuality = trackQuality;
226  selectionParameter_.maxDxyPV = iConfig.getUntrackedParameter<double>("maxDxyPV");
227  selectionParameter_.maxDzPV = iConfig.getUntrackedParameter<double>("maxDzPV");
228  selectionParameter_.maxChi2 = iConfig.getUntrackedParameter<double>("maxChi2");
229  selectionParameter_.maxDpOverP = iConfig.getUntrackedParameter<double>("maxDpOverP");
230  selectionParameter_.minOuterHit = selectionParameter_.minLayerCrossed = 0;
231  selectionParameter_.maxInMiss = selectionParameter_.maxOutMiss = 2;
232 
233  mergedDepth_ = (!isItPreRecHit_) || (collapseDepth_);
234  edm::LogVerbatim("HBHEMuon") << "Labels used: Track " << labelGenTrack_ << " Vtx " << labelVtx_ << " EB "
235  << labelEBRecHit_ << " EE " << labelEERecHit_ << " HBHE " << labelHBHERecHit_ << " MU "
236  << labelMuon_;
237 
238  if (!fileInCorr_.empty()) {
239  std::ifstream infile(fileInCorr_.c_str());
240  if (infile.is_open()) {
241  while (true) {
242  unsigned int id;
243  double cfac;
244  infile >> id >> cfac;
245  if (!infile.good())
246  break;
247  corrValue_[DetId(id)] = cfac;
248  }
249  infile.close();
250  }
251  }
252  useMyCorr_ = (!corrValue_.empty());
253  edm::LogVerbatim("HBHEMuon") << "Flags used: UseRaw " << useRaw_ << " GetCharge " << getCharge_ << " UnCorrect "
254  << unCorrect_ << " IgnoreHECorr " << ignoreHECorr_ << " CollapseDepth " << collapseDepth_
255  << ":" << mergedDepth_ << " IsItPlan1 " << isItPlan1_ << " IsItPreRecHit "
256  << isItPreRecHit_ << " UseMyCorr " << useMyCorr_;
257 }
TFileService::kSharedResource
static const std::string kSharedResource
Definition: TFileService.h:76
edm::LogVerbatim
Log< level::Info, true > LogVerbatim
Definition: MessageLogger.h:134
HcalHBHEMuonHighEtaAnalyzer::tok_EB_
const edm::EDGetTokenT< EcalRecHitCollection > tok_EB_
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:112
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
Definition: ParameterSet.h:303
HcalHBHEMuonHighEtaAnalyzer::selectionParameter_
spr::trackSelectionParameters selectionParameter_
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:130
HcalHBHEMuonHighEtaAnalyzer::tok_HBHE_
const edm::EDGetTokenT< HBHERecHitCollection > tok_HBHE_
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:114
reco::TrackBase::TrackQuality
TrackQuality
track quality
Definition: TrackBase.h:150
HcalHBHEMuonHighEtaAnalyzer::tok_topo_
const edm::ESGetToken< CaloTopology, CaloTopologyRecord > tok_topo_
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:125
HcalHBHEMuonHighEtaAnalyzer::tok_htopo_
const edm::ESGetToken< HcalTopology, HcalRecNumberingRecord > tok_htopo_
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:119
HcalHBHEMuonHighEtaAnalyzer::tok_magField_
const edm::ESGetToken< MagneticField, IdealMagneticFieldRecord > tok_magField_
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:122
StandaloneTrackMonitor_cfi.trackQuality
trackQuality
Definition: StandaloneTrackMonitor_cfi.py:16
l1ctLayer2EG_cff.id
id
Definition: l1ctLayer2EG_cff.py:85
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
HcalHBHEMuonHighEtaAnalyzer::tok_Muon_
const edm::EDGetTokenT< reco::MuonCollection > tok_Muon_
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:115
HcalHBHEMuonHighEtaAnalyzer::tok_Vtx_
const edm::EDGetTokenT< reco::VertexCollection > tok_Vtx_
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:111
HcalHBHEMuonHighEtaAnalyzer::tok_chan_
const edm::ESGetToken< EcalChannelStatus, EcalChannelStatusRcd > tok_chan_
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:123
HcalHBHEMuonHighEtaAnalyzer::tok_respcorr_
const edm::ESGetToken< HcalRespCorrs, HcalRespCorrsRcd > tok_respcorr_
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:120
HcalHBHEMuonHighEtaAnalyzer::tok_geom_
const edm::ESGetToken< CaloGeometry, CaloGeometryRecord > tok_geom_
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:121
HcalHBHEMuonHighEtaAnalyzer::tok_ddrec_
const edm::ESGetToken< HcalDDDRecConstants, HcalRecNumberingRecord > tok_ddrec_
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:118
DetId
Definition: DetId.h:17
reco::TrackBase::qualityByName
static TrackQuality qualityByName(const std::string &name)
Definition: TrackBase.cc:126
HcalHBHEMuonHighEtaAnalyzer::tok_genTrack_
const edm::EDGetTokenT< reco::TrackCollection > tok_genTrack_
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:116
HcalHBHEMuonHighEtaAnalyzer::tok_dbservice_
const edm::ESGetToken< HcalDbService, HcalDbRecord > tok_dbservice_
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:126
spr::trackSelectionParameters::minQuality
reco::TrackBase::TrackQuality minQuality
Definition: TrackSelection.h:21
HcalHBHEMuonHighEtaAnalyzer::tok_EE_
const edm::EDGetTokenT< EcalRecHitCollection > tok_EE_
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:113
HcalHBHEMuonHighEtaAnalyzer::tok_sevlv_
const edm::ESGetToken< EcalSeverityLevelAlgo, EcalSeverityLevelAlgoRcd > tok_sevlv_
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:124

Member Function Documentation

◆ activeLength()

double HcalHBHEMuonHighEtaAnalyzer::activeLength ( const DetId &  hid)
private

Definition at line 929 of file HcalHBHEMuonHighEtaAnalyzer.cc.

References actHB, actHE, hcalRecHitTable_cff::depth, spr::find(), HcalBarrel, hdc_, mps_fire::i, l1ctLayer2EG_cff::id, hcalRecHitTable_cff::ieta, hcalRecHitTable_cff::iphi, mergedDepth_, HcalDDDRecConstants::unmergeDepthDetId(), and ecaldqm::zside().

Referenced by analyzeTracks(), and beginRun().

929  {
930  HcalDetId id(hid);
931  int ieta = id.ietaAbs();
932  int zside = id.zside();
933  int iphi = id.iphi();
934  std::vector<int> dpths;
935  if (mergedDepth_) {
936  std::vector<HcalDetId> ids;
937  hdc_->unmergeDepthDetId(id, ids);
938  for (auto idh : ids)
939  dpths.emplace_back(idh.depth());
940  } else {
941  dpths.emplace_back(id.depth());
942  }
943  double lx(0);
944  if (id.subdet() == HcalBarrel) {
945  for (unsigned int i = 0; i < actHB.size(); ++i) {
946  if ((ieta == actHB[i].ieta) && (zside == actHB[i].zside) &&
947  (std::find(dpths.begin(), dpths.end(), actHB[i].depth) != dpths.end()) &&
948  (std::find(actHB[i].iphis.begin(), actHB[i].iphis.end(), iphi) != actHB[i].iphis.end())) {
949  lx += actHB[i].thick;
950  }
951  }
952  } else {
953  for (unsigned int i = 0; i < actHE.size(); ++i) {
954  if ((ieta == actHE[i].ieta) && (zside == actHE[i].zside) &&
955  (std::find(dpths.begin(), dpths.end(), actHE[i].depth) != dpths.end()) &&
956  (std::find(actHE[i].iphis.begin(), actHE[i].iphis.end(), iphi) != actHE[i].iphis.end())) {
957  lx += actHE[i].thick;
958  }
959  }
960  }
961  return lx;
962 }
mps_fire.i
i
Definition: mps_fire.py:429
HcalHBHEMuonHighEtaAnalyzer::actHE
std::vector< HcalDDDRecConstants::HcalActiveLength > actHE
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:173
ecaldqm::zside
int zside(DetId const &)
Definition: EcalDQMCommonUtils.cc:189
l1ctLayer2EG_cff.id
id
Definition: l1ctLayer2EG_cff.py:85
spr::find
void find(edm::Handle< EcalRecHitCollection > &hits, DetId thisDet, std::vector< EcalRecHitCollection::const_iterator > &hit, bool debug=false)
Definition: FindCaloHit.cc:19
hcalRecHitTable_cff.ieta
ieta
Definition: hcalRecHitTable_cff.py:15
hcalRecHitTable_cff.iphi
iphi
Definition: hcalRecHitTable_cff.py:16
hcalRecHitTable_cff.depth
depth
Definition: hcalRecHitTable_cff.py:17
HcalHBHEMuonHighEtaAnalyzer::actHB
std::vector< HcalDDDRecConstants::HcalActiveLength > actHB
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:173
HcalDDDRecConstants::unmergeDepthDetId
void unmergeDepthDetId(const HcalDetId &id, std::vector< HcalDetId > &ids) const
Definition: HcalDDDRecConstants.cc:696

◆ analyze()

void HcalHBHEMuonHighEtaAnalyzer::analyze ( edm::Event const &  iEvent,
edm::EventSetup const &  iSetup 
)
overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 358 of file HcalHBHEMuonHighEtaAnalyzer.cc.

References accept(), analyzeHadron(), analyzeMuon(), analyzeMuon_, barrelRecHitsHandle_, bField_, caloTopology_, clearVectors(), conditions_, depthMax_, ecal3x3Energy_, ecalDetId_, ecalEnergy_, ehcalDetId_, emaxNearP_, endcapRecHitsHandle_, energyMuon_, etaGlob_, eventNumber_, edm::EventSetup::getData(), goodVertex_, hbhe_, hcal1x1Energy_, hcal_ieta_, hcal_iphi_, hcalActiveLength_, hcalActiveLengthHot_, hcalDepthActiveLength_, hcalDepthActiveLengthHot_, hcalDepthChargeHot_, hcalDepthChargeHotBG_, hcalDepthEnergy_, hcalDepthEnergyCorr_, hcalDepthEnergyHot_, hcalDepthEnergyHotCorr_, hcalDepthMatch_, hcalDepthMatchHot_, hcalDetId_, hcalEnergy_, hcalHot_, hoEnergy_, mps_fire::i, iEvent, isGoodVertex(), isolationR03_, isolationR04_, kount_, matchedId_, mediumMuon_, phiGlob_, pMuon_, ptGlob_, runNumber_, sevlv_, theEcalChStatus_, tok_chan_, tok_dbservice_, tok_EB_, tok_EE_, tok_HBHE_, tok_magField_, tok_sevlv_, tok_topo_, tok_Vtx_, trackDz_, trackLayerCrossed_, trackMissedInnerHits_, trackMissedOuterHits_, trackOuterHit_, tree_, verbosity_, and L1BJetProducer_cff::vtx.

358  {
359  ++kount_;
360  clearVectors();
361  runNumber_ = iEvent.id().run();
362  eventNumber_ = iEvent.id().event();
363 #ifdef EDM_ML_DEBUG
364  edm::LogVerbatim("HBHEMuon") << "Run " << runNumber_ << " Event " << eventNumber_;
365 #endif
366 
367  // get handles to calogeometry and calotopology
368  bField_ = &iSetup.getData(tok_magField_);
369  theEcalChStatus_ = &iSetup.getData(tok_chan_);
370  sevlv_ = &iSetup.getData(tok_sevlv_);
371  caloTopology_ = &iSetup.getData(tok_topo_);
372  conditions_ = &iSetup.getData(tok_dbservice_);
373 
374  // Relevant blocks from iEvent
375  const edm::Handle<reco::VertexCollection>& vtx = iEvent.getHandle(tok_Vtx_);
376 
377  iEvent.getByToken(tok_EB_, barrelRecHitsHandle_);
378  iEvent.getByToken(tok_EE_, endcapRecHitsHandle_);
379  iEvent.getByToken(tok_HBHE_, hbhe_);
380 
381  // require a good vertex
382  math::XYZPoint pvx;
383  goodVertex_ = 0;
384  if (!vtx.isValid()) {
385 #ifdef EDM_ML_DEBUG
386  edm::LogVerbatim("HBHEMuon") << "No Good Vertex found == Reject\n";
387 #endif
388  return;
389  }
390 
391  reco::VertexCollection::const_iterator firstGoodVertex = vtx->end();
392  for (reco::VertexCollection::const_iterator it = vtx->begin(); it != vtx->end(); it++) {
393  if (isGoodVertex(*it)) {
394  if (firstGoodVertex == vtx->end())
395  firstGoodVertex = it;
396  ++goodVertex_;
397  }
398  }
399  if (firstGoodVertex != vtx->end())
400  pvx = firstGoodVertex->position();
401 
402  bool accept(false);
403  if (barrelRecHitsHandle_.isValid() && endcapRecHitsHandle_.isValid() && hbhe_.isValid()) {
404  accept = analyzeMuon_ ? analyzeMuon(iEvent, pvx) : analyzeHadron(iEvent, pvx);
405  }
406  if (accept) {
407 #ifdef EDM_ML_DEBUG
408  edm::LogVerbatim("HBHEMuon") << "Total of " << hcal_ieta_.size() << " propagated points";
409  for (unsigned int i = 0; i < hcal_ieta_.size(); ++i)
410  edm::LogVerbatim("HBHEMuon") << "[" << i << "] ieta/iphi for entry to "
411  << "HCAL has value of " << hcal_ieta_[i] << ":" << hcal_iphi_[i];
412  if ((verbosity_ / 100) % 10 > 0) {
413  edm::LogVerbatim("HBHEMuon") << "Sizes:: ptGlob:" << ptGlob_.size() << " etaGlob:" << etaGlob_.size()
414  << " phiGlob:" << phiGlob_.size() << " energyMuon:" << energyMuon_.size()
415  << " pMuon:" << pMuon_.size() << " mediumMuon: " << mediumMuon_.size()
416  << " isolation:" << isolationR04_.size() << ":" << isolationR03_.size()
417  << " e|h|ho energy: " << ecalEnergy_.size() << ":" << hcalEnergy_.size() << ":"
418  << hoEnergy_.size();
419  edm::LogVerbatim("HBHEMuon") << " matchedId:" << matchedId_.size() << " hcalHot:" << hcalHot_.size()
420  << " 3x3|1x1 energy:" << ecal3x3Energy_.size() << ":" << hcal1x1Energy_.size()
421  << " detId:" << ecalDetId_.size() << ":" << hcalDetId_.size() << ":"
422  << ehcalDetId_.size() << " eta|phi:" << hcal_ieta_.size() << ":"
423  << hcal_iphi_.size();
424  edm::LogVerbatim("HBHEMuon") << " activeLength:" << hcalActiveLength_.size() << ":"
425  << hcalActiveLengthHot_.size() << " emaxNearP:" << emaxNearP_.size()
426  << " trackDz: " << trackDz_.size() << " tracks:" << trackLayerCrossed_.size() << ":"
427  << trackOuterHit_.size() << ":" << trackMissedInnerHits_.size() << ":"
428  << trackMissedOuterHits_.size();
429  for (unsigned int i = 0; i < depthMax_; ++i)
430  edm::LogVerbatim("HBHEMuon")
431  << "Depth " << i
432  << " Energy|Length|EnergyHot|LengthHot|Charge|ChargeBG|EnergyCorr|EnergyHotCorr|Match|MatchHot:"
433  << hcalDepthEnergy_[i].size() << ":" << hcalDepthActiveLength_[i].size() << ":"
434  << hcalDepthEnergyHot_[i].size() << ":" << hcalDepthActiveLengthHot_[i].size() << ":"
435  << hcalDepthChargeHot_[i].size() << ":" << hcalDepthChargeHotBG_[i].size() << ":"
436  << hcalDepthEnergyCorr_[i].size() << ":" << hcalDepthEnergyHotCorr_[i].size() << ":"
437  << hcalDepthMatch_[i].size() << ":" << hcalDepthMatchHot_[i].size();
438  }
439 #endif
440  tree_->Fill();
441  }
442 }
edm::LogVerbatim
Log< level::Info, true > LogVerbatim
Definition: MessageLogger.h:134
HcalHBHEMuonHighEtaAnalyzer::tok_EB_
const edm::EDGetTokenT< EcalRecHitCollection > tok_EB_
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:112
HcalHBHEMuonHighEtaAnalyzer::theEcalChStatus_
const EcalChannelStatus * theEcalChStatus_
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:137
mps_fire.i
i
Definition: mps_fire.py:429
HcalHBHEMuonHighEtaAnalyzer::hbhe_
edm::Handle< HBHERecHitCollection > hbhe_
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:144
HcalHBHEMuonHighEtaAnalyzer::tok_HBHE_
const edm::EDGetTokenT< HBHERecHitCollection > tok_HBHE_
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:114
HcalHBHEMuonHighEtaAnalyzer::hcalDepthActiveLengthHot_
std::vector< double > hcalDepthActiveLengthHot_[depthMax_]
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:161
HcalHBHEMuonHighEtaAnalyzer::hcalDepthEnergyCorr_
std::vector< double > hcalDepthEnergyCorr_[depthMax_]
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:164
HcalHBHEMuonHighEtaAnalyzer::hcalDepthEnergy_
std::vector< double > hcalDepthEnergy_[depthMax_]
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:158
HcalHBHEMuonHighEtaAnalyzer::tok_topo_
const edm::ESGetToken< CaloTopology, CaloTopologyRecord > tok_topo_
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:125
L1BJetProducer_cff.vtx
vtx
Definition: L1BJetProducer_cff.py:10
HcalHBHEMuonHighEtaAnalyzer::sevlv_
const EcalSeverityLevelAlgo * sevlv_
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:138
HcalHBHEMuonHighEtaAnalyzer::hcalDepthMatchHot_
std::vector< bool > hcalDepthMatchHot_[depthMax_]
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:167
HcalHBHEMuonHighEtaAnalyzer::tok_magField_
const edm::ESGetToken< MagneticField, IdealMagneticFieldRecord > tok_magField_
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:122
HcalHBHEMuonHighEtaAnalyzer::analyzeMuon
bool analyzeMuon(edm::Event const &, math::XYZPoint &)
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:497
HcalHBHEMuonHighEtaAnalyzer::hcalDetId_
std::vector< unsigned int > hcalDetId_
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:156
accept
bool accept(const edm::Event &event, const edm::TriggerResults &triggerTable, const std::string &triggerPath)
Definition: TopDQMHelpers.h:31
HcalHBHEMuonHighEtaAnalyzer::hcalDepthChargeHotBG_
std::vector< double > hcalDepthChargeHotBG_[depthMax_]
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:163
HcalHBHEMuonHighEtaAnalyzer::isGoodVertex
bool isGoodVertex(const reco::Vertex &)
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:964
iEvent
int iEvent
Definition: GenABIO.cc:224
HcalHBHEMuonHighEtaAnalyzer::tok_Vtx_
const edm::EDGetTokenT< reco::VertexCollection > tok_Vtx_
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:111
HcalHBHEMuonHighEtaAnalyzer::tok_chan_
const edm::ESGetToken< EcalChannelStatus, EcalChannelStatusRcd > tok_chan_
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:123
HcalHBHEMuonHighEtaAnalyzer::hcalDepthEnergyHot_
std::vector< double > hcalDepthEnergyHot_[depthMax_]
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:160
HcalHBHEMuonHighEtaAnalyzer::hcalDepthActiveLength_
std::vector< double > hcalDepthActiveLength_[depthMax_]
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:159
HcalHBHEMuonHighEtaAnalyzer::hcalDepthEnergyHotCorr_
std::vector< double > hcalDepthEnergyHotCorr_[depthMax_]
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:165
math::XYZPoint
XYZPointD XYZPoint
point in space with cartesian internal representation
Definition: Point3D.h:12
HcalHBHEMuonHighEtaAnalyzer::tok_dbservice_
const edm::ESGetToken< HcalDbService, HcalDbRecord > tok_dbservice_
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:126
HcalHBHEMuonHighEtaAnalyzer::ehcalDetId_
std::vector< unsigned int > ehcalDetId_
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:156
HcalHBHEMuonHighEtaAnalyzer::barrelRecHitsHandle_
edm::Handle< EcalRecHitCollection > barrelRecHitsHandle_
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:142
HcalHBHEMuonHighEtaAnalyzer::hcalDepthChargeHot_
std::vector< double > hcalDepthChargeHot_[depthMax_]
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:162
HcalHBHEMuonHighEtaAnalyzer::tok_EE_
const edm::EDGetTokenT< EcalRecHitCollection > tok_EE_
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:113
HcalHBHEMuonHighEtaAnalyzer::hcalDepthMatch_
std::vector< bool > hcalDepthMatch_[depthMax_]
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:166
HcalHBHEMuonHighEtaAnalyzer::ecalDetId_
std::vector< unsigned int > ecalDetId_
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:156
HcalHBHEMuonHighEtaAnalyzer::endcapRecHitsHandle_
edm::Handle< EcalRecHitCollection > endcapRecHitsHandle_
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:143
HcalHBHEMuonHighEtaAnalyzer::tok_sevlv_
const edm::ESGetToken< EcalSeverityLevelAlgo, EcalSeverityLevelAlgoRcd > tok_sevlv_
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:124
HcalHBHEMuonHighEtaAnalyzer::analyzeHadron
bool analyzeHadron(edm::Event const &, math::XYZPoint &)
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:560

◆ analyzeHadron()

bool HcalHBHEMuonHighEtaAnalyzer::analyzeHadron ( edm::Event const &  iEvent,
math::XYZPoint &  leadPV 
)
private

Definition at line 560 of file HcalHBHEMuonHighEtaAnalyzer.cc.

References funct::abs(), accept(), analyzeTracks(), bField_, reco::TrackBase::eta(), etaMin_, geo_, iEvent, edm::HandleBase::isValid(), l1ctLayer1_cff::nTrack, spr::propagateCALO(), theTrackQuality_, and tok_genTrack_.

Referenced by analyze().

560  {
561  //Get track collection
562  edm::Handle<reco::TrackCollection> trkCollection = iEvent.getHandle(tok_genTrack_);
563  bool accept = false;
564 
565  if (!trkCollection.isValid()) {
566  std::vector<spr::propagatedTrackID> trkCaloDets;
567  spr::propagateCALO(trkCollection, geo_, bField_, theTrackQuality_, trkCaloDets, false);
568  int nTrack(0);
569  std::vector<spr::propagatedTrackID>::const_iterator trkDetItr;
570  for (trkDetItr = trkCaloDets.begin(), nTrack = 0; trkDetItr != trkCaloDets.end(); trkDetItr++, nTrack++) {
571  const reco::Track* pTrack = &(*(trkDetItr->trkItr));
572  if (std::abs(pTrack->eta()) > etaMin_) {
573  accept = analyzeTracks(pTrack, leadPV, nTrack, trkCaloDets, true);
574  }
575  }
576  }
577  return accept;
578 }
spr::propagateCALO
std::vector< spr::propagatedTrackID > propagateCALO(edm::Handle< reco::TrackCollection > &trkCollection, const CaloGeometry *geo, const MagneticField *bField, const std::string &theTrackQuality, bool debug=false)
Definition: CaloPropagateTrack.cc:130
l1ctLayer1_cff.nTrack
nTrack
Definition: l1ctLayer1_cff.py:60
accept
bool accept(const edm::Event &event, const edm::TriggerResults &triggerTable, const std::string &triggerPath)
Definition: TopDQMHelpers.h:31
iEvent
int iEvent
Definition: GenABIO.cc:224
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
reco::TrackBase::eta
double eta() const
pseudorapidity of momentum vector
Definition: TrackBase.h:652
HcalHBHEMuonHighEtaAnalyzer::analyzeTracks
bool analyzeTracks(const reco::Track *, math::XYZPoint &, int, std::vector< spr::propagatedTrackID > &, bool)
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:580
HcalHBHEMuonHighEtaAnalyzer::tok_genTrack_
const edm::EDGetTokenT< reco::TrackCollection > tok_genTrack_
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:116
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:70

◆ analyzeMuon()

bool HcalHBHEMuonHighEtaAnalyzer::analyzeMuon ( edm::Event const &  iEvent,
math::XYZPoint &  leadPV 
)
private

Definition at line 497 of file HcalHBHEMuonHighEtaAnalyzer.cc.

References funct::abs(), accept(), analyzeTracks(), ecalEnergy_, energyMuon_, reco::TrackBase::eta(), etaGlob_, etaMin_, hcalEnergy_, hoEnergy_, iEvent, isolationR03_, isolationR04_, edm::HandleBase::isValid(), SiStripPI::max, mediumMuon_, l1ctLayer1_cff::nTrack, phiGlob_, pMuon_, edm::Handle< T >::product(), ptGlob_, muon::segmentCompatibility(), tok_Muon_, and verbosity_.

Referenced by analyze().

497  {
498  const edm::Handle<reco::MuonCollection>& _Muon = iEvent.getHandle(tok_Muon_);
499  bool accept = false;
500 
501  if (_Muon.isValid()) {
502  int nTrack(0);
503  std::vector<spr::propagatedTrackID> trkCaloDets;
504  for (const auto& RecMuon : (*(_Muon.product()))) {
505  if (RecMuon.innerTrack().isNonnull()) {
506  const reco::Track* pTrack = (RecMuon.innerTrack()).get();
507  if (std::abs(pTrack->eta()) > etaMin_) {
508  if (analyzeTracks(pTrack, leadPV, nTrack, trkCaloDets, false)) {
509  accept = true;
510  ptGlob_.emplace_back(RecMuon.pt());
511  etaGlob_.emplace_back(RecMuon.eta());
512  phiGlob_.emplace_back(RecMuon.phi());
513  energyMuon_.push_back(RecMuon.energy());
514  pMuon_.emplace_back(RecMuon.p());
515  bool mediumMuon = (((RecMuon.isPFMuon()) && (RecMuon.isGlobalMuon() || RecMuon.isTrackerMuon())) &&
516  (RecMuon.innerTrack()->validFraction() > 0.49));
517  if (mediumMuon) {
518  double chiGlobal = ((RecMuon.globalTrack().isNonnull()) ? RecMuon.globalTrack()->normalizedChi2() : 999);
519  bool goodGlob =
520  (RecMuon.isGlobalMuon() && chiGlobal < 3 && RecMuon.combinedQuality().chi2LocalPosition < 12 &&
521  RecMuon.combinedQuality().trkKink < 20);
522  mediumMuon = muon::segmentCompatibility(RecMuon) > (goodGlob ? 0.303 : 0.451);
523  }
524  mediumMuon_.emplace_back(mediumMuon);
525  bool isoR03 =
526  ((RecMuon.pfIsolationR03().sumChargedHadronPt +
527  std::max(0.,
528  RecMuon.pfIsolationR03().sumNeutralHadronEt + RecMuon.pfIsolationR03().sumPhotonEt -
529  (0.5 * RecMuon.pfIsolationR03().sumPUPt))) /
530  RecMuon.pt());
531  bool isoR04 =
532  ((RecMuon.pfIsolationR04().sumChargedHadronPt +
533  std::max(0.,
534  RecMuon.pfIsolationR04().sumNeutralHadronEt + RecMuon.pfIsolationR04().sumPhotonEt -
535  (0.5 * RecMuon.pfIsolationR04().sumPUPt))) /
536  RecMuon.pt());
537  isolationR03_.emplace_back(isoR03);
538  isolationR04_.emplace_back(isoR04);
539 
540  ecalEnergy_.emplace_back(RecMuon.calEnergy().emS9);
541  hcalEnergy_.emplace_back(RecMuon.calEnergy().hadS9);
542  hoEnergy_.emplace_back(RecMuon.calEnergy().hoS9);
543 #ifdef EDM_ML_DEBUG
544  if ((verbosity_ / 100) % 10 > 0)
545  edm::LogVerbatim("HBHEMuon")
546  << "Muon[" << ptGlob_.size() << "] pt:eta:phi:p " << ptGlob_.back() << ":" << etaGlob_.back() << ":"
547  << phiGlob_.back() << ":" << energyMuon_.back() << ":" << pMuon_.back() << ":"
548  << " Medium:i3:i4 " << mediumMuon_.back() << ":" << isolationR03_.back() << ":"
549  << isolationR04_.back() << ":"
550  << " Energy EC:HC:HO " << ecalEnergy_.back() << ":" << hcalEnergy_.back() << ":" << hoEnergy_.back();
551 #endif
552  }
553  }
554  }
555  }
556  }
557  return accept;
558 }
l1ctLayer1_cff.nTrack
nTrack
Definition: l1ctLayer1_cff.py:60
edm::Handle::product
T const * product() const
Definition: Handle.h:70
accept
bool accept(const edm::Event &event, const edm::TriggerResults &triggerTable, const std::string &triggerPath)
Definition: TopDQMHelpers.h:31
iEvent
int iEvent
Definition: GenABIO.cc:224
HcalHBHEMuonHighEtaAnalyzer::tok_Muon_
const edm::EDGetTokenT< reco::MuonCollection > tok_Muon_
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:115
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
muon::segmentCompatibility
float segmentCompatibility(const reco::Muon &muon, reco::Muon::ArbitrationType arbitrationType=reco::Muon::SegmentAndTrackArbitration)
Definition: MuonSelectors.cc:61
reco::TrackBase::eta
double eta() const
pseudorapidity of momentum vector
Definition: TrackBase.h:652
HcalHBHEMuonHighEtaAnalyzer::analyzeTracks
bool analyzeTracks(const reco::Track *, math::XYZPoint &, int, std::vector< spr::propagatedTrackID > &, bool)
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:580
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:70

◆ analyzeTracks()

bool HcalHBHEMuonHighEtaAnalyzer::analyzeTracks ( const reco::Track *  pTrack,
math::XYZPoint &  leadPV,
int  nTrack,
std::vector< spr::propagatedTrackID > &  trkCaloDets,
bool  ifHadron 
)
private

Definition at line 580 of file HcalHBHEMuonHighEtaAnalyzer.cc.

References funct::abs(), accept(), activeLength(), barrelRecHitsHandle_, bField_, caloTopology_, spr::chargeIsolationEcal(), RPCNoise_example::check, chg, collapseDepth_, conditions_, alignCSCRings::corr, TauDecayModes::dec, hcalRecHitTable_cff::depth, depth16HE(), depthMax_, spr::propagatedTrackID::detIdECAL, spr::propagatedTrackID::detIdEHCAL, spr::propagatedTrackID::detIdHCAL, ecal3x3Energy_, ecalDetId_, spr::eECALmatrix(), ehcalDetId_, spr::eHCALmatrix(), emaxNearP_, emaxNearPThr_, endcapRecHitsHandle_, spr::energyHCALCell(), fillTrackParameters(), PedestalClient_cfi::gain, gainFactor(), geo_, getCharge_, goodCell(), spr::goodTrack(), hbhe_, hcal1x1Energy_, hcal_ieta_, hcal_iphi_, hcalActiveLength_, hcalActiveLengthHot_, HcalBarrel, hcalDepthActiveLength_, hcalDepthActiveLengthHot_, hcalDepthChargeHot_, hcalDepthChargeHotBG_, hcalDepthEnergy_, hcalDepthEnergyCorr_, hcalDepthEnergyHot_, hcalDepthEnergyHotCorr_, hcalDepthMatch_, hcalDepthMatchHot_, hcalDetId_, HcalEndcap, hcalHot_, hdc_, mps_fire::i, hcalRecHitTable_cff::ieta, HcalDetId::ieta(), ignoreHECorr_, info(), hcalRecHitTable_cff::iphi, HcalDetId::iphi(), isItPlan1_, isItPreRecHit_, dqmdumpme::k, matchedId_, matchId(), maxDepth_, HcalDDDRecConstants::mergedDepthDetId(), l1ctLayer1_cff::nTrack, spr::propagatedTrackID::okECAL, spr::propagatedTrackID::okHCAL, spr::propagateCALO(), spr::propagateHCALBack(), DetId::rawId(), respCorr(), edm::second(), selectionParameter_, sevlv_, HcalDetId::subdet(), theEcalChStatus_, theHBHETopology_, unCorrect_, useRaw_, HcalTopology::validHcal(), and verbosity_.

Referenced by analyzeHadron(), and analyzeMuon().

584  {
585  bool accept(false);
586 
587  if (spr::goodTrack(pTrack, leadPV, selectionParameter_, false)) {
588  spr::propagatedTrackID trackID = spr::propagateCALO(pTrack, geo_, bField_, false);
589 
590  if (trackID.okECAL && trackID.okHCAL) {
591  double emaxNearP = (ifHadron) ? spr::chargeIsolationEcal(nTrack, trkCaloDets, geo_, caloTopology_, 15, 15) : 0;
592  if (emaxNearP < emaxNearPThr_) {
593  double eEcal(0), eHcal(0), activeLengthTot(0), activeLengthHotTot(0);
594  double eHcalDepth[depthMax_], eHcalDepthHot[depthMax_];
595  double eHcalDepthC[depthMax_], eHcalDepthHotC[depthMax_];
596  double cHcalDepthHot[depthMax_], cHcalDepthHotBG[depthMax_];
597  double activeL[depthMax_], activeHotL[depthMax_];
598  bool matchDepth[depthMax_], matchDepthHot[depthMax_];
599  HcalDetId eHcalDetId[depthMax_];
600  unsigned int isHot(0);
601  bool tmpmatch(false);
602  int ieta(-1000), iphi(-1000);
603  for (int i = 0; i < depthMax_; ++i) {
604  eHcalDepth[i] = eHcalDepthHot[i] = 0;
605  eHcalDepthC[i] = eHcalDepthHotC[i] = 0;
606  cHcalDepthHot[i] = cHcalDepthHotBG[i] = 0;
607  activeL[i] = activeHotL[i] = 0;
608  matchDepth[i] = matchDepthHot[i] = true;
609  }
610 
611  HcalDetId check;
612  std::pair<bool, HcalDetId> info = spr::propagateHCALBack(pTrack, geo_, bField_, false);
613  if (info.first)
614  check = info.second;
615 
616  const DetId isoCell(trackID.detIdECAL);
617  std::pair<double, bool> e3x3 = spr::eECALmatrix(isoCell,
618  barrelRecHitsHandle_,
619  endcapRecHitsHandle_,
620  *theEcalChStatus_,
621  geo_,
622  caloTopology_,
623  sevlv_,
624  1,
625  1,
626  -100.0,
627  -100.0,
628  -500.0,
629  500.0,
630  false);
631  eEcal = e3x3.first;
632 #ifdef EDM_ML_DEBUG
633  if (verbosity_ % 10 > 0)
634  edm::LogVerbatim("HBHEMuon") << "Propagate Track to ECAL: " << e3x3.second << ":" << trackID.okECAL << " E "
635  << eEcal;
636 #endif
637 
638  DetId closestCell(trackID.detIdHCAL);
639  HcalDetId hcidt(closestCell.rawId());
640  if ((hcidt.ieta() == check.ieta()) && (hcidt.iphi() == check.iphi()))
641  tmpmatch = true;
642 #ifdef EDM_ML_DEBUG
643  if (verbosity_ % 10 > 0)
644  edm::LogVerbatim("HBHEMuon") << "Front " << hcidt << " Back " << info.first << ":" << check << " Match "
645  << tmpmatch;
646 #endif
647 
648  HcalSubdetector subdet = hcidt.subdet();
649  ieta = hcidt.ieta();
650  iphi = hcidt.iphi();
651  bool hborhe = (std::abs(ieta) == 16);
652 
653  eHcal = spr::eHCALmatrix(theHBHETopology_,
654  closestCell,
655  hbhe_,
656  0,
657  0,
658  false,
659  true,
660  -100.0,
661  -100.0,
662  -100.0,
663  -100.0,
664  -500.,
665  500.,
666  useRaw_);
667  std::vector<std::pair<double, int>> ehdepth;
668  spr::energyHCALCell((HcalDetId)closestCell,
669  hbhe_,
670  ehdepth,
671  depthMax_,
672  -100.0,
673  -100.0,
674  -100.0,
675  -100.0,
676  -500.0,
677  500.0,
678  useRaw_,
679  depth16HE(ieta, iphi),
680  false);
681  for (int i = 0; i < depthMax_; ++i)
682  eHcalDetId[i] = HcalDetId();
683  for (unsigned int i = 0; i < ehdepth.size(); ++i) {
684  HcalSubdetector subdet0 =
685  (hborhe) ? ((ehdepth[i].second >= depth16HE(ieta, iphi)) ? HcalEndcap : HcalBarrel) : subdet;
686  HcalDetId hcid0(subdet0, ieta, iphi, ehdepth[i].second);
687  double actL = activeLength(DetId(hcid0));
688  double ene = ehdepth[i].first;
689  bool tmpC(false);
690  if (ene > 0.0) {
691  if (!(theHBHETopology_->validHcal(hcid0))) {
692  edm::LogWarning("HBHEMuon") << "(1) Invalid ID " << hcid0 << " with E = " << ene;
693  edm::LogWarning("HBHEMuon") << HcalDetId(closestCell) << " with " << ehdepth.size() << " depths:";
694  for (const auto& ehd : ehdepth)
695  edm::LogWarning("HBHEMuon") << " " << ehd.second << ":" << ehd.first;
696  } else {
697  tmpC = goodCell(hcid0, pTrack, geo_, bField_);
698  double enec(ene);
699  if (unCorrect_) {
700  double corr = (ignoreHECorr_ && (subdet0 == HcalEndcap)) ? 1.0 : respCorr(DetId(hcid0));
701  if (corr != 0)
702  ene /= corr;
703 #ifdef EDM_ML_DEBUG
704  if (verbosity_ % 10 > 0) {
705  HcalDetId id = (isItPlan1_ && isItPreRecHit_) ? hdc_->mergedDepthDetId(hcid0) : hcid0;
706  edm::LogVerbatim("HBHEMuon") << hcid0 << ":" << id << " Corr " << corr;
707  }
708 #endif
709  }
710  int depth = ehdepth[i].second - 1;
711  if (collapseDepth_) {
712  HcalDetId id = hdc_->mergedDepthDetId(hcid0);
713  depth = id.depth() - 1;
714  }
715  eHcalDepth[depth] += ene;
716  eHcalDepthC[depth] += enec;
717  activeL[depth] += actL;
718  activeLengthTot += actL;
719  matchDepth[depth] = (matchDepth[depth] && tmpC);
720 #ifdef EDM_ML_DEBUG
721  if ((verbosity_ / 10) % 10 > 0)
722  edm::LogVerbatim("HBHEMuon")
723  << hcid0 << " E " << ene << ":" << enec << " L " << actL << " Match " << tmpC;
724 #endif
725  }
726  }
727  }
728 #ifdef EDM_ML_DEBUG
729  if ((verbosity_ / 10) % 10 > 0) {
730  edm::LogVerbatim("HBHEMuon") << hcidt << " Match " << tmpmatch << " Depths " << ehdepth.size();
731  for (unsigned int k = 0; k < ehdepth.size(); ++k)
732  edm::LogVerbatim("HBHEMuon") << " [" << k << ":" << ehdepth[k].second << "] " << matchDepth[k];
733  }
734 #endif
735  HcalDetId hotCell;
736  spr::eHCALmatrix(geo_, theHBHETopology_, closestCell, hbhe_, 1, 1, hotCell, false, useRaw_, false);
737  isHot = matchId(closestCell, hotCell);
738  if (hotCell != HcalDetId()) {
739  subdet = HcalDetId(hotCell).subdet();
740  ieta = HcalDetId(hotCell).ieta();
741  iphi = HcalDetId(hotCell).iphi();
742  hborhe = (std::abs(ieta) == 16);
743  std::vector<std::pair<double, int>> ehdepth;
744  spr::energyHCALCell(hotCell,
745  hbhe_,
746  ehdepth,
747  depthMax_,
748  -100.0,
749  -100.0,
750  -100.0,
751  -100.0,
752  -500.0,
753  500.0,
754  useRaw_,
755  depth16HE(ieta, iphi),
756  false);
757  for (int i = 0; i < depthMax_; ++i)
758  eHcalDetId[i] = HcalDetId();
759  for (unsigned int i = 0; i < ehdepth.size(); ++i) {
760  HcalSubdetector subdet0 =
761  (hborhe) ? ((ehdepth[i].second >= depth16HE(ieta, iphi)) ? HcalEndcap : HcalBarrel) : subdet;
762  HcalDetId hcid0(subdet0, ieta, iphi, ehdepth[i].second);
763  double actL = activeLength(DetId(hcid0));
764  double ene = ehdepth[i].first;
765  bool tmpC(false);
766  if (ene > 0.0) {
767  if (!(theHBHETopology_->validHcal(hcid0))) {
768  edm::LogWarning("HBHEMuon") << "(2) Invalid ID " << hcid0 << " with E = " << ene;
769  edm::LogWarning("HBHEMuon") << HcalDetId(hotCell) << " with " << ehdepth.size() << " depths:";
770  for (const auto& ehd : ehdepth)
771  edm::LogWarning("HBHEMuon") << " " << ehd.second << ":" << ehd.first;
772  } else {
773  tmpC = goodCell(hcid0, pTrack, geo_, bField_);
774  double chg(ene), enec(ene);
775  if (unCorrect_) {
776  double corr = (ignoreHECorr_ && (subdet0 == HcalEndcap)) ? 1.0 : respCorr(DetId(hcid0));
777  if (corr != 0)
778  ene /= corr;
779 #ifdef EDM_ML_DEBUG
780  if (verbosity_ % 10 > 0) {
781  HcalDetId id = (isItPlan1_ && isItPreRecHit_) ? hdc_->mergedDepthDetId(hcid0) : hcid0;
782  edm::LogVerbatim("HBHEMuon")
783  << hcid0 << ":" << id << " Corr " << corr << " E " << ene << ":" << enec;
784  }
785 #endif
786  }
787  if (getCharge_) {
788  double gain = gainFactor(conditions_, hcid0);
789  if (gain != 0)
790  chg /= gain;
791 #ifdef EDM_ML_DEBUG
792  if (verbosity_ % 10 > 0)
793  edm::LogVerbatim("HBHEMuon") << hcid0 << " Gain " << gain << " C " << chg;
794 #endif
795  }
796  int depth = ehdepth[i].second - 1;
797  if (collapseDepth_) {
798  HcalDetId id = hdc_->mergedDepthDetId(hcid0);
799  depth = id.depth() - 1;
800  }
801  eHcalDepthHot[depth] += ene;
802  eHcalDepthHotC[depth] += enec;
803  cHcalDepthHot[depth] += chg;
804  activeHotL[depth] += actL;
805  activeLengthHotTot += actL;
806  matchDepthHot[depth] = (matchDepthHot[depth] && tmpC);
807 #ifdef EDM_ML_DEBUG
808  if ((verbosity_ / 10) % 10 > 0)
809  edm::LogVerbatim("HBHEMuon") << hcid0 << " depth " << depth << " E " << ene << ":" << enec << " C "
810  << chg << " L " << actL << " Match " << tmpC;
811 #endif
812  }
813  }
814  }
815  }
816 #ifdef EDM_ML_DEBUG
817  edm::LogVerbatim("HBHEMuon") << "Propagate Track to HCAL: " << trackID.okHCAL << " Match " << tmpmatch
818  << " Hot " << isHot << " Energy " << eHcal;
819 #endif
820 
821  accept = true;
822  ecalDetId_.emplace_back((trackID.detIdECAL)());
823  hcalDetId_.emplace_back((trackID.detIdHCAL)());
824  ehcalDetId_.emplace_back((trackID.detIdEHCAL)());
825  emaxNearP_.emplace_back(emaxNearP);
826  matchedId_.emplace_back(tmpmatch);
827  ecal3x3Energy_.emplace_back(eEcal);
828  hcal1x1Energy_.emplace_back(eHcal);
829  hcal_ieta_.emplace_back(ieta);
830  hcal_iphi_.emplace_back(iphi);
831  for (int i = 0; i < maxDepth_; ++i) {
832  hcalDepthEnergy_[i].emplace_back(eHcalDepth[i]);
833  hcalDepthActiveLength_[i].emplace_back(activeL[i]);
834  hcalDepthEnergyHot_[i].emplace_back(eHcalDepthHot[i]);
835  hcalDepthActiveLengthHot_[i].emplace_back(activeHotL[i]);
836  hcalDepthEnergyCorr_[i].emplace_back(eHcalDepthC[i]);
837  hcalDepthEnergyHotCorr_[i].emplace_back(eHcalDepthHotC[i]);
838  hcalDepthChargeHot_[i].emplace_back(cHcalDepthHot[i]);
839  hcalDepthChargeHotBG_[i].emplace_back(cHcalDepthHotBG[i]);
840  hcalDepthMatch_[i].emplace_back(matchDepth[i]);
841  hcalDepthMatchHot_[i].emplace_back(matchDepthHot[i]);
842  }
843  hcalActiveLength_.emplace_back(activeLengthTot);
844  hcalHot_.emplace_back(isHot);
845  hcalActiveLengthHot_.emplace_back(activeLengthHotTot);
846 #ifdef EDM_ML_DEBUG
847  if ((verbosity_ / 100) % 10 > 0) {
848  edm::LogVerbatim("HBHEMuon") << "Track " << std::hex << ecalDetId_.back() << ":" << hcalDetId_.back() << ":"
849  << ehcalDetId_.back() << std::dec << ":" << emaxNearP_.back() << ":"
850  << matchedId_.back() << ":" << ecal3x3Energy_.back() << ":"
851  << hcal1x1Energy_.back() << ":" << hcal_ieta_.back() << ":" << hcal_iphi_.back()
852  << ":" << hcalActiveLength_.back() << ":" << hcalHot_.back() << ":"
853  << hcalActiveLengthHot_.back();
854  for (int i = 0; i < maxDepth_; ++i) {
855  edm::LogVerbatim("HBHEMuon") << "Depth[" << i << "] " << hcalDepthEnergy_[i].back() << ":"
856  << hcalDepthActiveLength_[i].back() << ":" << hcalDepthEnergyHot_[i].back()
857  << ":" << hcalDepthActiveLengthHot_[i].back() << ":"
858  << hcalDepthEnergyCorr_[i].back() << ":" << hcalDepthEnergyHotCorr_[i].back()
859  << ":" << hcalDepthChargeHot_[i].back() << ":"
860  << hcalDepthChargeHotBG_[i].back() << ":" << hcalDepthMatch_[i].back() << ":"
861  << hcalDepthMatchHot_[i].back();
862  }
863  }
864 #endif
865  fillTrackParameters(pTrack, leadPV);
866  }
867  }
868  }
869  return accept;
870 }
HcalHBHEMuonHighEtaAnalyzer::matchId
int matchId(const HcalDetId &, const HcalDetId &)
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:922
edm::LogVerbatim
Log< level::Info, true > LogVerbatim
Definition: MessageLogger.h:134
spr::propagateCALO
std::vector< spr::propagatedTrackID > propagateCALO(edm::Handle< reco::TrackCollection > &trkCollection, const CaloGeometry *geo, const MagneticField *bField, const std::string &theTrackQuality, bool debug=false)
Definition: CaloPropagateTrack.cc:130
HcalHBHEMuonHighEtaAnalyzer::theEcalChStatus_
const EcalChannelStatus * theEcalChStatus_
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:137
info
static const TGPicture * info(bool iBackgroundIsBlack)
Definition: FWCollectionSummaryWidget.cc:153
mps_fire.i
i
Definition: mps_fire.py:429
RPCNoise_example.check
check
Definition: RPCNoise_example.py:71
l1ctLayer1_cff.nTrack
nTrack
Definition: l1ctLayer1_cff.py:60
HcalHBHEMuonHighEtaAnalyzer::selectionParameter_
spr::trackSelectionParameters selectionParameter_
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:130
HcalHBHEMuonHighEtaAnalyzer::hbhe_
edm::Handle< HBHERecHitCollection > hbhe_
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:144
HcalDDDRecConstants::mergedDepthDetId
HcalDetId mergedDepthDetId(const HcalDetId &id) const
Definition: HcalDDDRecConstants.cc:672
HcalHBHEMuonHighEtaAnalyzer::hcalDepthActiveLengthHot_
std::vector< double > hcalDepthActiveLengthHot_[depthMax_]
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:161
HcalHBHEMuonHighEtaAnalyzer::hcalDepthEnergyCorr_
std::vector< double > hcalDepthEnergyCorr_[depthMax_]
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:164
chg
const float chg[109]
Definition: CoreSimTrack.cc:5
HcalHBHEMuonHighEtaAnalyzer::hcalDepthEnergy_
std::vector< double > hcalDepthEnergy_[depthMax_]
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:158
HcalHBHEMuonHighEtaAnalyzer::sevlv_
const EcalSeverityLevelAlgo * sevlv_
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:138
HcalHBHEMuonHighEtaAnalyzer::hcalDepthMatchHot_
std::vector< bool > hcalDepthMatchHot_[depthMax_]
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:167
HcalHBHEMuonHighEtaAnalyzer::hcalDetId_
std::vector< unsigned int > hcalDetId_
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:156
HcalHBHEMuonHighEtaAnalyzer::goodCell
bool goodCell(const HcalDetId &, const reco::Track *, const CaloGeometry *, const MagneticField *)
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:1012
accept
bool accept(const edm::Event &event, const edm::TriggerResults &triggerTable, const std::string &triggerPath)
Definition: TopDQMHelpers.h:31
HcalHBHEMuonHighEtaAnalyzer::hcalDepthChargeHotBG_
std::vector< double > hcalDepthChargeHotBG_[depthMax_]
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:163
edm::second
U second(std::pair< T, U > const &p)
Definition: ParameterSet.cc:222
HcalHBHEMuonHighEtaAnalyzer::gainFactor
double gainFactor(const HcalDbService *, const HcalDetId &)
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:988
HcalDetId::subdet
constexpr HcalSubdetector subdet() const
get the subdetector
Definition: HcalDetId.h:138
hcalRecHitTable_cff.ieta
ieta
Definition: hcalRecHitTable_cff.py:15
spr::chargeIsolationEcal
double chargeIsolationEcal(unsigned int trkIndex, std::vector< spr::propagatedTrackID > &vdetIds, const CaloGeometry *geo, const CaloTopology *caloTopology, int ieta, int iphi, bool debug=false)
Definition: ChargeIsolation.cc:12
alignCSCRings.corr
dictionary corr
Definition: alignCSCRings.py:124
spr::propagateHCALBack
spr::propagatedTrackDirection propagateHCALBack(unsigned int thisTrk, edm::Handle< edm::SimTrackContainer > &SimTk, edm::Handle< edm::SimVertexContainer > &SimVtx, const CaloGeometry *geo, const MagneticField *bField, bool debug=false)
Definition: CaloPropagateTrack.cc:583
spr::goodTrack
bool goodTrack(const reco::Track *pTrack, math::XYZPoint leadPV, trackSelectionParameters parameters, bool debug=false)
Definition: TrackSelection.cc:9
HcalDetId::ieta
constexpr int ieta() const
get the cell ieta
Definition: HcalDetId.h:155
HcalHBHEMuonHighEtaAnalyzer::hcalDepthEnergyHot_
std::vector< double > hcalDepthEnergyHot_[depthMax_]
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:160
HcalSubdetector
HcalSubdetector
Definition: HcalAssistant.h:31
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
PedestalClient_cfi.gain
gain
Definition: PedestalClient_cfi.py:37
hcalRecHitTable_cff.iphi
iphi
Definition: hcalRecHitTable_cff.py:16
HcalHBHEMuonHighEtaAnalyzer::fillTrackParameters
void fillTrackParameters(const reco::Track *, math::XYZPoint)
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:1022
HcalHBHEMuonHighEtaAnalyzer::hcalDepthActiveLength_
std::vector< double > hcalDepthActiveLength_[depthMax_]
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:159
DetId
Definition: DetId.h:17
HcalHBHEMuonHighEtaAnalyzer::hcalDepthEnergyHotCorr_
std::vector< double > hcalDepthEnergyHotCorr_[depthMax_]
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:165
HcalHBHEMuonHighEtaAnalyzer::ehcalDetId_
std::vector< unsigned int > ehcalDetId_
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:156
HcalTopology::validHcal
bool validHcal(const HcalDetId &id) const
Definition: HcalTopology.cc:230
hcalRecHitTable_cff.depth
depth
Definition: hcalRecHitTable_cff.py:17
HcalHBHEMuonHighEtaAnalyzer::barrelRecHitsHandle_
edm::Handle< EcalRecHitCollection > barrelRecHitsHandle_
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:142
HcalHBHEMuonHighEtaAnalyzer::hcalDepthChargeHot_
std::vector< double > hcalDepthChargeHot_[depthMax_]
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:162
edm::LogWarning
Log< level::Warning, false > LogWarning
Definition: MessageLogger.h:128
HcalHBHEMuonHighEtaAnalyzer::hcalDepthMatch_
std::vector< bool > hcalDepthMatch_[depthMax_]
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:166
HcalHBHEMuonHighEtaAnalyzer::ecalDetId_
std::vector< unsigned int > ecalDetId_
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:156
spr::energyHCALCell
void energyHCALCell(HcalDetId detId, edm::Handle< T > &hits, std::vector< std::pair< double, int > > &energyCell, int maxDepth=1, double hbThr=-100, double heThr=-100, double hfThr=-100, double hoThr=-100, double tMin=-500, double tMax=500, int useRaw=0, int depthHE=3, bool debug=false)
HcalHBHEMuonHighEtaAnalyzer::endcapRecHitsHandle_
edm::Handle< EcalRecHitCollection > endcapRecHitsHandle_
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:143
HcalDetId::iphi
constexpr int iphi() const
get the cell iphi
Definition: HcalDetId.h:157
spr::eECALmatrix
double eECALmatrix(const DetId &detId, edm::Handle< T > &hitsEB, edm::Handle< T > &hitsEE, const CaloGeometry *geo, const CaloTopology *caloTopology, int ieta, int iphi, double ebThr=-100, double eeThr=-100, double tMin=-500, double tMax=500, bool debug=false)
spr::eHCALmatrix
double eHCALmatrix(const HcalTopology *topology, const DetId &det, edm::Handle< T > &hits, int ieta, int iphi, bool includeHO=false, bool algoNew=true, double hbThr=-100, double heThr=-100, double hfThr=-100, double hoThr=-100, double tMin=-500, double tMax=500, int useRaw=0, bool debug=false)

◆ beginJob()

void HcalHBHEMuonHighEtaAnalyzer::beginJob ( )
overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 296 of file HcalHBHEMuonHighEtaAnalyzer.cc.

References ecal3x3Energy_, ecalDetId_, ecalEnergy_, ehcalDetId_, emaxNearP_, energyMuon_, etaGlob_, eventNumber_, compareTotals::fs, goodVertex_, hcal1x1Energy_, hcal_ieta_, hcal_iphi_, hcalActiveLength_, hcalActiveLengthHot_, hcalDepthActiveLength_, hcalDepthActiveLengthHot_, hcalDepthChargeHot_, hcalDepthChargeHotBG_, hcalDepthEnergy_, hcalDepthEnergyCorr_, hcalDepthEnergyHot_, hcalDepthEnergyHotCorr_, hcalDepthMatch_, hcalDepthMatchHot_, hcalDetId_, hcalEnergy_, hcalHot_, hoEnergy_, isolationR03_, isolationR04_, dqmdumpme::k, matchedId_, maxDepth_, mediumMuon_, Skims_PA_cff::name, phiGlob_, pMuon_, ptGlob_, runNumber_, trackDz_, trackLayerCrossed_, trackMissedInnerHits_, trackMissedOuterHits_, trackOuterHit_, and tree_.

296  {
297  edm::Service<TFileService> fs;
298  tree_ = fs->make<TTree>("HBHEMuonHighEta", "HBHEMuonHighEta");
299  tree_->Branch("pt_of_muon", &ptGlob_);
300  tree_->Branch("eta_of_muon", &etaGlob_);
301  tree_->Branch("phi_of_muon", &phiGlob_);
302  tree_->Branch("energy_of_muon", &energyMuon_);
303  tree_->Branch("p_of_muon", &pMuon_);
304  tree_->Branch("MediumMuon", &mediumMuon_);
305  tree_->Branch("IsolationR04", &isolationR04_);
306  tree_->Branch("IsolationR03", &isolationR03_);
307  tree_->Branch("ecal_3into3", &ecalEnergy_);
308  tree_->Branch("hcal_3into3", &hcalEnergy_);
309  tree_->Branch("ho_3into3", &hoEnergy_);
310  tree_->Branch("emaxNearP", &emaxNearP_);
311 
312  tree_->Branch("Run_No", &runNumber_);
313  tree_->Branch("Event_No", &eventNumber_);
314  tree_->Branch("GoodVertex", &goodVertex_);
315  tree_->Branch("matchedId", &matchedId_);
316  tree_->Branch("hcal_cellHot", &hcalHot_);
317  tree_->Branch("ecal_3x3", &ecal3x3Energy_);
318  tree_->Branch("hcal_1x1", &hcal1x1Energy_);
319  tree_->Branch("ecal_detID", &ecalDetId_);
320  tree_->Branch("hcal_detID", &hcalDetId_);
321  tree_->Branch("ehcal_detID", &ehcalDetId_);
322  tree_->Branch("hcal_ieta", &hcal_ieta_);
323  tree_->Branch("hcal_iphi", &hcal_iphi_);
324 
325  char name[100];
326  for (int k = 0; k < maxDepth_; ++k) {
327  sprintf(name, "hcal_edepth%d", (k + 1));
328  tree_->Branch(name, &hcalDepthEnergy_[k]);
329  sprintf(name, "hcal_activeL%d", (k + 1));
330  tree_->Branch(name, &hcalDepthActiveLength_[k]);
331  sprintf(name, "hcal_edepthHot%d", (k + 1));
332  tree_->Branch(name, &hcalDepthEnergyHot_[k]);
333  sprintf(name, "hcal_activeHotL%d", (k + 1));
334  tree_->Branch(name, &hcalDepthActiveLengthHot_[k]);
335  sprintf(name, "hcal_cdepthHot%d", (k + 1));
336  tree_->Branch(name, &hcalDepthChargeHot_[k]);
337  sprintf(name, "hcal_cdepthHotBG%d", (k + 1));
338  tree_->Branch(name, &hcalDepthChargeHotBG_[k]);
339  sprintf(name, "hcal_edepthCorrect%d", (k + 1));
340  tree_->Branch(name, &hcalDepthEnergyCorr_[k]);
341  sprintf(name, "hcal_edepthHotCorrect%d", (k + 1));
342  tree_->Branch(name, &hcalDepthEnergyHotCorr_[k]);
343  sprintf(name, "hcal_depthMatch%d", (k + 1));
344  tree_->Branch(name, &hcalDepthMatch_[k]);
345  sprintf(name, "hcal_depthMatchHot%d", (k + 1));
346  tree_->Branch(name, &hcalDepthMatchHot_[k]);
347  }
348  tree_->Branch("activeLength", &hcalActiveLength_);
349  tree_->Branch("activeLengthHot", &hcalActiveLengthHot_);
350  tree_->Branch("trackDz", &trackDz_);
351  tree_->Branch("trackLayerCrossed", &trackLayerCrossed_);
352  tree_->Branch("trackOuterHit", &trackOuterHit_);
353  tree_->Branch("trackMissedInnerHits", &trackMissedInnerHits_);
354  tree_->Branch("trackMissedOuterHits", &trackMissedOuterHits_);
355 }
HcalHBHEMuonHighEtaAnalyzer::hcalDepthActiveLengthHot_
std::vector< double > hcalDepthActiveLengthHot_[depthMax_]
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:161
HcalHBHEMuonHighEtaAnalyzer::hcalDepthEnergyCorr_
std::vector< double > hcalDepthEnergyCorr_[depthMax_]
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:164
HcalHBHEMuonHighEtaAnalyzer::hcalDepthEnergy_
std::vector< double > hcalDepthEnergy_[depthMax_]
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:158
HcalHBHEMuonHighEtaAnalyzer::hcalDepthMatchHot_
std::vector< bool > hcalDepthMatchHot_[depthMax_]
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:167
HcalHBHEMuonHighEtaAnalyzer::hcalDetId_
std::vector< unsigned int > hcalDetId_
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:156
HcalHBHEMuonHighEtaAnalyzer::hcalDepthChargeHotBG_
std::vector< double > hcalDepthChargeHotBG_[depthMax_]
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:163
HcalHBHEMuonHighEtaAnalyzer::hcalDepthEnergyHot_
std::vector< double > hcalDepthEnergyHot_[depthMax_]
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:160
HcalHBHEMuonHighEtaAnalyzer::hcalDepthActiveLength_
std::vector< double > hcalDepthActiveLength_[depthMax_]
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:159
HcalHBHEMuonHighEtaAnalyzer::hcalDepthEnergyHotCorr_
std::vector< double > hcalDepthEnergyHotCorr_[depthMax_]
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:165
HcalHBHEMuonHighEtaAnalyzer::ehcalDetId_
std::vector< unsigned int > ehcalDetId_
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:156
HcalHBHEMuonHighEtaAnalyzer::hcalDepthChargeHot_
std::vector< double > hcalDepthChargeHot_[depthMax_]
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:162
HcalHBHEMuonHighEtaAnalyzer::hcalDepthMatch_
std::vector< bool > hcalDepthMatch_[depthMax_]
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:166
HcalHBHEMuonHighEtaAnalyzer::ecalDetId_
std::vector< unsigned int > ecalDetId_
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:156
Skims_PA_cff.name
name
Definition: Skims_PA_cff.py:17

◆ beginRun()

void HcalHBHEMuonHighEtaAnalyzer::beginRun ( edm::Run const &  iRun,
edm::EventSetup const &  iSetup 
)
overrideprivate

Definition at line 445 of file HcalHBHEMuonHighEtaAnalyzer.cc.

References actHB, actHE, activeLength(), geo_, edm::EventSetup::getData(), CaloGeometry::getSubdetectorGeometry(), HcalDDDRecConstants::getThickActive(), HcalGeometry::getValidDetIds(), reco::JetExtendedAssociation::getValue(), HcalCondObjectContainer< Item >::getValues(), DetId::Hcal, HcalBarrel, HcalEndcap, hdc_, mergedDepth_, HcalDDDRecConstants::mergedDepthDetId(), DetId::rawId(), respCorrs_, edm::RunBase::run(), HcalCondObjectContainerBase::setTopo(), theHBHETopology_, tok_ddrec_, tok_geom_, tok_htopo_, tok_respcorr_, verbosity_, and writeRespCorr_.

445  {
446  hdc_ = &iSetup.getData(tok_ddrec_);
447  actHB.clear();
448  actHE.clear();
449  actHB = hdc_->getThickActive(0);
450  actHE = hdc_->getThickActive(1);
451 #ifdef EDM_ML_DEBUG
452  if (verbosity_ % 10 > 0) {
453  unsigned int k1(0), k2(0);
454  edm::LogVerbatim("HBHEMuon") << actHB.size() << " Active Length for HB";
455  for (const auto& act : actHB) {
456  edm::LogVerbatim("HBHEMuon") << "[" << k1 << "] ieta " << act.ieta << " depth " << act.depth << " zside "
457  << act.zside << " type " << act.stype << " phi " << act.iphis.size() << ":"
458  << act.iphis[0] << " L " << act.thick;
459  HcalDetId hcid1(HcalBarrel, (act.ieta) * (act.zside), act.iphis[0], act.depth);
460  HcalDetId hcid2 = mergedDepth_ ? hdc_->mergedDepthDetId(hcid1) : hcid1;
461  edm::LogVerbatim("HBHEMuon") << hcid1 << " | " << hcid2 << " L " << activeLength(DetId(hcid2));
462  ++k1;
463  }
464  edm::LogVerbatim("HBHEMuon") << actHE.size() << " Active Length for HE";
465  for (const auto& act : actHE) {
466  edm::LogVerbatim("HBHEMuon") << "[" << k2 << "] ieta " << act.ieta << " depth " << act.depth << " zside "
467  << act.zside << " type " << act.stype << " phi " << act.iphis.size() << ":"
468  << act.iphis[0] << " L " << act.thick;
469  HcalDetId hcid1(HcalEndcap, (act.ieta) * (act.zside), act.iphis[0], act.depth);
470  HcalDetId hcid2 = mergedDepth_ ? hdc_->mergedDepthDetId(hcid1) : hcid1;
471  edm::LogVerbatim("HBHEMuon") << hcid1 << " | " << hcid2 << " L " << activeLength(DetId(hcid2));
472  ++k2;
473  }
474  }
475 #endif
476 
477  theHBHETopology_ = &iSetup.getData(tok_htopo_);
478  const HcalRespCorrs* resp = &iSetup.getData(tok_respcorr_);
479  respCorrs_ = new HcalRespCorrs(*resp);
480  respCorrs_->setTopo(theHBHETopology_);
481  geo_ = &iSetup.getData(tok_geom_);
482 
483  // Write correction factors for all HB/HE events
484  if (writeRespCorr_) {
485  const HcalGeometry* gHcal = (const HcalGeometry*)(geo_->getSubdetectorGeometry(DetId::Hcal, HcalBarrel));
486  const std::vector<DetId>& ids = gHcal->getValidDetIds(DetId::Hcal, 0);
487  edm::LogVerbatim("HBHEMuon") << "\nTable of Correction Factors for Run " << iRun.run() << "\n";
488  for (auto const& id : ids) {
489  if ((id.det() == DetId::Hcal) && ((id.subdetId() == HcalBarrel) || (id.subdetId() == HcalEndcap))) {
490  edm::LogVerbatim("HBHEMuon") << HcalDetId(id) << " " << id.rawId() << " "
491  << (respCorrs_->getValues(id))->getValue();
492  }
493  }
494  }
495 }
edm::LogVerbatim
Log< level::Info, true > LogVerbatim
Definition: MessageLogger.h:134
HcalGeometry::getValidDetIds
const std::vector< DetId > & getValidDetIds(DetId::Detector det=DetId::Detector(0), int subdet=0) const override
Get a list of valid detector ids (for the given subdetector)
Definition: HcalGeometry.cc:76
HcalDDDRecConstants::mergedDepthDetId
HcalDetId mergedDepthDetId(const HcalDetId &id) const
Definition: HcalDDDRecConstants.cc:672
HcalHBHEMuonHighEtaAnalyzer::actHE
std::vector< HcalDDDRecConstants::HcalActiveLength > actHE
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:173
HcalHBHEMuonHighEtaAnalyzer::tok_htopo_
const edm::ESGetToken< HcalTopology, HcalRecNumberingRecord > tok_htopo_
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:119
HcalCondObjectContainer::getValues
const Item * getValues(DetId fId, bool throwOnFail=true) const
Definition: HcalCondObjectContainer.h:159
HcalHBHEMuonHighEtaAnalyzer::tok_respcorr_
const edm::ESGetToken< HcalRespCorrs, HcalRespCorrsRcd > tok_respcorr_
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:120
HcalHBHEMuonHighEtaAnalyzer::tok_geom_
const edm::ESGetToken< CaloGeometry, CaloGeometryRecord > tok_geom_
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:121
HcalHBHEMuonHighEtaAnalyzer::tok_ddrec_
const edm::ESGetToken< HcalDDDRecConstants, HcalRecNumberingRecord > tok_ddrec_
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:118
DetId
Definition: DetId.h:17
DetId::rawId
constexpr uint32_t rawId() const
get the raw id
Definition: DetId.h:57
HcalHBHEMuonHighEtaAnalyzer::actHB
std::vector< HcalDDDRecConstants::HcalActiveLength > actHB
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:173
HcalDDDRecConstants::getThickActive
std::vector< HcalActiveLength > getThickActive(const int &type) const
Definition: HcalDDDRecConstants.cc:490
reco::JetExtendedAssociation::getValue
const JetExtendedData & getValue(const Container &, const reco::JetBaseRef &)
get value for the association. Throw exception if no association found
Definition: JetExtendedAssociation.cc:50
CaloGeometry::getSubdetectorGeometry
const CaloSubdetectorGeometry * getSubdetectorGeometry(const DetId &id) const
access the subdetector geometry for the given subdetector directly
Definition: CaloGeometry.cc:34
HcalCondObjectContainerBase::setTopo
void setTopo(const HcalTopology *topo)
Definition: HcalCondObjectContainerBase.cc:17

◆ clearVectors()

void HcalHBHEMuonHighEtaAnalyzer::clearVectors ( )
private

clearing vectots

Definition at line 872 of file HcalHBHEMuonHighEtaAnalyzer.cc.

References depthMax_, ecal3x3Energy_, ecalDetId_, ecalEnergy_, ehcalDetId_, emaxNearP_, energyMuon_, etaGlob_, eventNumber_, goodVertex_, hcal1x1Energy_, hcal_ieta_, hcal_iphi_, hcalActiveLength_, hcalActiveLengthHot_, hcalDepthActiveLength_, hcalDepthActiveLengthHot_, hcalDepthChargeHot_, hcalDepthChargeHotBG_, hcalDepthEnergy_, hcalDepthEnergyCorr_, hcalDepthEnergyHot_, hcalDepthEnergyHotCorr_, hcalDepthMatch_, hcalDepthMatchHot_, hcalDetId_, hcalEnergy_, hcalHot_, hoEnergy_, mps_fire::i, isolationR03_, isolationR04_, matchedId_, mediumMuon_, phiGlob_, pMuon_, ptGlob_, runNumber_, trackDz_, trackLayerCrossed_, trackMissedInnerHits_, trackMissedOuterHits_, and trackOuterHit_.

Referenced by analyze().

872  {
874  eventNumber_ = -99999;
875  runNumber_ = -99999;
876  goodVertex_ = -99999;
877 
878  mediumMuon_.clear();
879  ptGlob_.clear();
880  etaGlob_.clear();
881  phiGlob_.clear();
882  energyMuon_.clear();
883  pMuon_.clear();
884  isolationR04_.clear();
885  isolationR03_.clear();
886  ecalEnergy_.clear();
887  hcalEnergy_.clear();
888  hoEnergy_.clear();
889 
890  matchedId_.clear();
891  hcalHot_.clear();
892  ecal3x3Energy_.clear();
893  hcal1x1Energy_.clear();
894  ecalDetId_.clear();
895  hcalDetId_.clear();
896  ehcalDetId_.clear();
897  hcal_ieta_.clear();
898  hcal_iphi_.clear();
899  for (int i = 0; i < depthMax_; ++i) {
900  hcalDepthEnergy_[i].clear();
901  hcalDepthActiveLength_[i].clear();
902  hcalDepthEnergyHot_[i].clear();
903  hcalDepthActiveLengthHot_[i].clear();
904  hcalDepthChargeHot_[i].clear();
905  hcalDepthChargeHotBG_[i].clear();
906  hcalDepthEnergyCorr_[i].clear();
907  hcalDepthEnergyHotCorr_[i].clear();
908  hcalDepthMatch_[i].clear();
909  hcalDepthMatchHot_[i].clear();
910  }
911  hcalActiveLength_.clear();
912  hcalActiveLengthHot_.clear();
913 
914  emaxNearP_.clear();
915  trackDz_.clear();
916  trackLayerCrossed_.clear();
917  trackOuterHit_.clear();
918  trackMissedInnerHits_.clear();
919  trackMissedOuterHits_.clear();
920 }
mps_fire.i
i
Definition: mps_fire.py:429
HcalHBHEMuonHighEtaAnalyzer::hcalDepthActiveLengthHot_
std::vector< double > hcalDepthActiveLengthHot_[depthMax_]
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:161
HcalHBHEMuonHighEtaAnalyzer::hcalDepthEnergyCorr_
std::vector< double > hcalDepthEnergyCorr_[depthMax_]
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:164
HcalHBHEMuonHighEtaAnalyzer::hcalDepthEnergy_
std::vector< double > hcalDepthEnergy_[depthMax_]
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:158
HcalHBHEMuonHighEtaAnalyzer::hcalDepthMatchHot_
std::vector< bool > hcalDepthMatchHot_[depthMax_]
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:167
HcalHBHEMuonHighEtaAnalyzer::hcalDetId_
std::vector< unsigned int > hcalDetId_
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:156
HcalHBHEMuonHighEtaAnalyzer::hcalDepthChargeHotBG_
std::vector< double > hcalDepthChargeHotBG_[depthMax_]
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:163
HcalHBHEMuonHighEtaAnalyzer::hcalDepthEnergyHot_
std::vector< double > hcalDepthEnergyHot_[depthMax_]
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:160
HcalHBHEMuonHighEtaAnalyzer::hcalDepthActiveLength_
std::vector< double > hcalDepthActiveLength_[depthMax_]
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:159
HcalHBHEMuonHighEtaAnalyzer::hcalDepthEnergyHotCorr_
std::vector< double > hcalDepthEnergyHotCorr_[depthMax_]
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:165
HcalHBHEMuonHighEtaAnalyzer::ehcalDetId_
std::vector< unsigned int > ehcalDetId_
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:156
HcalHBHEMuonHighEtaAnalyzer::hcalDepthChargeHot_
std::vector< double > hcalDepthChargeHot_[depthMax_]
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:162
HcalHBHEMuonHighEtaAnalyzer::hcalDepthMatch_
std::vector< bool > hcalDepthMatch_[depthMax_]
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:166
HcalHBHEMuonHighEtaAnalyzer::ecalDetId_
std::vector< unsigned int > ecalDetId_
Definition: HcalHBHEMuonHighEtaAnalyzer.cc:156

◆ depth16HE()

int HcalHBHEMuonHighEtaAnalyzer::depth16HE ( int  ieta,
int  iphi 
)
private

Definition at line 996 of file HcalHBHEMuonHighEtaAnalyzer.cc.

References HcalTopology::dddConstants(), hcalRecHitTable_cff::depth, HcalDDDRecConstants::getMinDepth(), hcalRecHitTable_cff::ieta, hcalRecHitTable_cff::iphi, isItPlan1_, isItPreRecHit_, theHBHETopology_, verbosity_, and ecaldqm::zside().

Referenced by analyzeTracks().

996  {
997  // Transition between HB/HE is special
998  // For Run 1 or for Plan1 standard reconstruction it is 3
999  // For runs beyond 2018 or in Plan1 for HEP17 it is 4
1000  int zside = (ieta > 0) ? 1 : -1;
1001  int depth = theHBHETopology_->dddConstants()->getMinDepth(1, 16, iphi, zside);
1002  if (isItPlan1_ && (!isItPreRecHit_))
1003  depth = 3;
1004 #ifdef EDM_ML_DEBUG
1005  if (verbosity_ % 10 > 0)
1006  edm::LogVerbatim("HBHEMuon") << "Plan1 " << isItPlan1_ << " PreRecHit " << isItPreRecHit_ << " phi " << iphi
1007  << " depth " << depth;
1008 #endif
1009  return depth;
1010 }
HcalDDDRecConstants::getMinDepth
int getMinDepth(const int &itype, const int &ieta, const int &iphi, const int &zside) const
Definition: HcalDDDRecConstants.cc:378
edm::LogVerbatim
Log< level::Info, true > LogVerbatim
Definition: MessageLogger.h:134
ecaldqm::zside
int zside(DetId const &)
Definition: EcalDQMCommonUtils.cc:189
hcalRecHitTable_cff.ieta
ieta
Definition: hcalRecHitTable_cff.py:15
hcalRecHitTable_cff.iphi
iphi
Definition: hcalRecHitTable_cff.py:16
hcalRecHitTable_cff.depth
depth
Definition: hcalRecHitTable_cff.py:17
HcalTopology::dddConstants
const HcalDDDRecConstants * dddConstants() const
Definition: HcalTopology.h:164

◆ endRun()

void HcalHBHEMuonHighEtaAnalyzer::endRun ( edm::Run const &  ,
edm::EventSetup const &   
)
inlineoverrideprivate

Definition at line 86 of file HcalHBHEMuonHighEtaAnalyzer.cc.

86 {}

◆ fillDescriptions()

void HcalHBHEMuonHighEtaAnalyzer::fillDescriptions ( edm::ConfigurationDescriptions &  descriptions)
static

Definition at line 264 of file HcalHBHEMuonHighEtaAnalyzer.cc.

References edm::ConfigurationDescriptions::add(), submitPVResolutionJobs::desc, ProducerED_cfi::InputTag, and AlCaHLTBitMon_QueryRunRegistry::string.

264  {
265  edm::ParameterSetDescription desc;
266  desc.add<edm::InputTag>("labelEBRecHit", edm::InputTag("ecalRecHit", "EcalRecHitsEB"));
267  desc.add<edm::InputTag>("labelEERecHit", edm::InputTag("ecalRecHit", "EcalRecHitsEE"));
268  desc.add<edm::InputTag>("labelHBHERecHit", edm::InputTag("hbhereco"));
269  desc.add<std::string>("labelVertex", "offlinePrimaryVertices");
270  desc.add<std::string>("labelMuon", "muons");
271  desc.add<std::string>("labelTrack", "generalTracks");
272  desc.add<double>("etaMin", 2.0);
273  desc.add<double>("emaxNearPThreshold", 10.0);
274  desc.add<bool>("analyzeMuon", true);
275  desc.add<bool>("unCorrect", true);
276  desc.add<bool>("collapseDepth", false);
277  desc.add<bool>("isItPlan1", false);
278  desc.add<bool>("getCharge", true);
279  desc.add<int>("useRaw", 0);
280  desc.add<int>("verbosity", 0);
281  desc.addUntracked<std::string>("fileInCorr", "");
282  desc.addUntracked<std::string>("trackQuality", "highPurity");
283  desc.addUntracked<double>("minTrackPt", 1.0);
284  desc.addUntracked<double>("maxDxyPV", 0.02);
285  desc.addUntracked<double>("maxDzPV", 100.0);
286  desc.addUntracked<double>("maxChi2", 5.0);
287  desc.addUntracked<double>("maxDpOverP", 0.1);
288  desc.addUntracked<bool>("ignoreHECorr", false);
289  desc.addUntracked<bool>("isItPreRecHit", false);
290  desc.addUntracked<bool>("writeRespCorr", false);
291  desc.addUntracked<int>("maxDepth", 7);
292  descriptions.add("hcalHBHEMuonHighEta", desc);
293 }
ProducerED_cfi.InputTag
InputTag
Definition: ProducerED_cfi.py:5
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
edm::ConfigurationDescriptions::add
void add(std::string const &label, ParameterSetDescription const &psetDescription)
Definition: ConfigurationDescriptions.cc:57

◆ fillTrackParameters()

void HcalHBHEMuonHighEtaAnalyzer::fillTrackParameters ( const reco::Track *  pTrack,
math::XYZPoint  leadPV 
)
private

Definition at line 1022 of file HcalHBHEMuonHighEtaAnalyzer.cc.

References reco::TrackBase::dz(), reco::TrackBase::hitPattern(), reco::HitPattern::MISSING_INNER_HITS, reco::HitPattern::MISSING_OUTER_HITS, reco::HitPattern::stripTECLayersWithMeasurement(), reco::HitPattern::stripTOBLayersWithMeasurement(), trackDz_, reco::HitPattern::trackerLayersWithMeasurement(), reco::HitPattern::trackerLayersWithoutMeasurement(), trackLayerCrossed_, trackMissedInnerHits_, trackMissedOuterHits_, and trackOuterHit_.

Referenced by analyzeTracks().

1022  {
1023  trackDz_.emplace_back(pTrack->dz(leadPV));
1024  const reco::HitPattern& hitp = pTrack->hitPattern();
1025  trackLayerCrossed_.emplace_back(hitp.trackerLayersWithMeasurement());
1026  trackOuterHit_.emplace_back(hitp.stripTOBLayersWithMeasurement() + hitp.stripTECLayersWithMeasurement());
1027  trackMissedInnerHits_.emplace_back(hitp.trackerLayersWithoutMeasurement(reco::HitPattern::MISSING_INNER_HITS));
1028  trackMissedOuterHits_.emplace_back(hitp.trackerLayersWithoutMeasurement(reco::HitPattern::MISSING_OUTER_HITS));
1029 }
reco::HitPattern::trackerLayersWithMeasurement
int trackerLayersWithMeasurement() const
Definition: HitPattern.cc:534
reco::HitPattern::stripTECLayersWithMeasurement
int stripTECLayersWithMeasurement() const
Definition: HitPattern.cc:628
reco::TrackBase::dz
double dz() const
dz parameter (= dsz/cos(lambda)). This is the track z0 w.r.t (0,0,0) only if the refPoint is close to...
Definition: TrackBase.h:622
reco::HitPattern::stripTOBLayersWithMeasurement
int stripTOBLayersWithMeasurement() const
Definition: HitPattern.cc:618
reco::HitPattern::trackerLayersWithoutMeasurement
int trackerLayersWithoutMeasurement(HitCategory category) const
Definition: HitPattern.cc:553
reco::TrackBase::hitPattern
const HitPattern & hitPattern() const
Access the hit pattern, indicating in which Tracker layers the track has hits.
Definition: TrackBase.h:504

◆ gainFactor()

double HcalHBHEMuonHighEtaAnalyzer::gainFactor ( const HcalDbService *  conditions,
const HcalDetId &  id 
)
private

Definition at line 988 of file HcalHBHEMuonHighEtaAnalyzer.cc.

References submitPVValidationJobs::conditions, PedestalClient_cfi::gain, and HcalCalibrations::respcorrgain().

Referenced by analyzeTracks().

988  {
989  double gain(0.0);
990  const HcalCalibrations& calibs = conditions->getHcalCalibrations(id);
991  for (int capid = 0; capid < 4; ++capid)
992  gain += (0.25 * calibs.respcorrgain(capid));
993  return gain;
994 }
PedestalClient_cfi.gain
gain
Definition: PedestalClient_cfi.py:37
HcalCalibrations::respcorrgain
constexpr double respcorrgain(int fCapId) const
get response corrected gain for capid=0..3
Definition: HcalCalibrations.h:42

◆ goodCell()

bool HcalHBHEMuonHighEtaAnalyzer::goodCell ( const HcalDetId &  hcid,
const reco::Track *  pTrack,
const CaloGeometry *  geo,
const MagneticField *  bField 
)
private

Definition at line 1012 of file HcalHBHEMuonHighEtaAnalyzer.cc.

References Calorimetry_cff::bField, HcalDDDRecConstants::getRZ(), HcalEndcap, hdc_, match(), spr::propagateHCAL(), and HcalDetId::subdet().

Referenced by analyzeTracks().

1015  {
1016  std::pair<double, double> rz = hdc_->getRZ(hcid);
1017  bool typeRZ = (hcid.subdet() == HcalEndcap) ? false : true;
1018  bool match = spr::propagateHCAL(pTrack, geo, bField, typeRZ, rz, false);
1019  return match;
1020 }
HcalDDDRecConstants::getRZ
double getRZ(const int &subdet, const int &ieta, const int &depth) const
Definition: HcalDDDRecConstants.cc:441
spr::propagateHCAL
std::pair< math::XYZPoint, bool > propagateHCAL(const reco::Track *, const MagneticField *, bool debug=false)
Definition: CaloPropagateTrack.cc:720
Calorimetry_cff.bField
bField
Definition: Calorimetry_cff.py:284
HcalDetId::subdet
constexpr HcalSubdetector subdet() const
get the subdetector
Definition: HcalDetId.h:138
match
std::pair< typename Association::data_type::first_type, double > match(Reference key, Association association, bool bestMatchByMaxValue)
Generic matching function.
Definition: Utils.h:10

◆ isGoodVertex()

bool HcalHBHEMuonHighEtaAnalyzer::isGoodVertex ( const reco::Vertex &  vtx)
private

Definition at line 964 of file HcalHBHEMuonHighEtaAnalyzer.cc.

References L1BJetProducer_cff::vtx.

Referenced by analyze().

964  {
965  if (vtx.isFake())
966  return false;
967  if (vtx.ndof() < 4)
968  return false;
969  if (vtx.position().Rho() > 2.)
970  return false;
971  if (fabs(vtx.position().Z()) > 24.)
972  return false;
973  return true;
974 }
L1BJetProducer_cff.vtx
vtx
Definition: L1BJetProducer_cff.py:10

◆ matchId()

int HcalHBHEMuonHighEtaAnalyzer::matchId ( const HcalDetId &  id1,
const HcalDetId &  id2 
)
private

Definition at line 922 of file HcalHBHEMuonHighEtaAnalyzer.cc.

References globals_cff::id1, globals_cff::id2, and match().

Referenced by analyzeTracks().

922  {
923  HcalDetId kd1(id1.subdet(), id1.ieta(), id1.iphi(), 1);
924  HcalDetId kd2(id1.subdet(), id2.ieta(), id2.iphi(), 1);
925  int match = ((kd1 == kd2) ? 1 : 0);
926  return match;
927 }
globals_cff.id1
id1
Definition: globals_cff.py:49
globals_cff.id2
id2
Definition: globals_cff.py:50
match
std::pair< typename Association::data_type::first_type, double > match(Reference key, Association association, bool bestMatchByMaxValue)
Generic matching function.
Definition: Utils.h:10

◆ respCorr()

double HcalHBHEMuonHighEtaAnalyzer::respCorr ( const DetId &  id)
private

Definition at line 976 of file HcalHBHEMuonHighEtaAnalyzer.cc.

References corrValue_, reco::JetExtendedAssociation::getValue(), HcalCondObjectContainer< Item >::getValues(), respCorrs_, and useMyCorr_.

Referenced by analyzeTracks().

976  {
977  double cfac(1.0);
978  if (useMyCorr_) {
979  auto itr = corrValue_.find(id);
980  if (itr != corrValue_.end())
981  cfac = itr->second;
982  } else if (respCorrs_ != nullptr) {
983  cfac = (respCorrs_->getValues(id))->getValue();
984  }
985  return cfac;
986 }
HcalCondObjectContainer::getValues
const Item * getValues(DetId fId, bool throwOnFail=true) const
Definition: HcalCondObjectContainer.h:159
reco::JetExtendedAssociation::getValue
const JetExtendedData & getValue(const Container &, const reco::JetBaseRef &)
get value for the association. Throw exception if no association found
Definition: JetExtendedAssociation.cc:50

Member Data Documentation

◆ actHB

std::vector<HcalDDDRecConstants::HcalActiveLength> HcalHBHEMuonHighEtaAnalyzer::actHB
private

Definition at line 173 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by activeLength(), and beginRun().

◆ actHE

std::vector<HcalDDDRecConstants::HcalActiveLength> HcalHBHEMuonHighEtaAnalyzer::actHE
private

Definition at line 173 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by activeLength(), and beginRun().

◆ analyzeMuon_

const bool HcalHBHEMuonHighEtaAnalyzer::analyzeMuon_
private

Definition at line 105 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyze().

◆ barrelRecHitsHandle_

edm::Handle<EcalRecHitCollection> HcalHBHEMuonHighEtaAnalyzer::barrelRecHitsHandle_
private

Definition at line 142 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyze(), and analyzeTracks().

◆ bField_

const MagneticField* HcalHBHEMuonHighEtaAnalyzer::bField_
private

Definition at line 136 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyze(), analyzeHadron(), and analyzeTracks().

◆ caloTopology_

const CaloTopology* HcalHBHEMuonHighEtaAnalyzer::caloTopology_
private

Definition at line 139 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyze(), and analyzeTracks().

◆ collapseDepth_

const bool HcalHBHEMuonHighEtaAnalyzer::collapseDepth_
private

Definition at line 105 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyzeTracks(), and HcalHBHEMuonHighEtaAnalyzer().

◆ conditions_

const HcalDbService* HcalHBHEMuonHighEtaAnalyzer::conditions_
private

Definition at line 140 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyze(), and analyzeTracks().

◆ corrValue_

std::map<DetId, double> HcalHBHEMuonHighEtaAnalyzer::corrValue_
private

Definition at line 174 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by HcalHBHEMuonHighEtaAnalyzer(), and respCorr().

◆ depthMax_

const int HcalHBHEMuonHighEtaAnalyzer::depthMax_ = 7
staticprivate

Definition at line 147 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyze(), analyzeTracks(), and clearVectors().

◆ ecal3x3Energy_

std::vector<double> HcalHBHEMuonHighEtaAnalyzer::ecal3x3Energy_
private

Definition at line 155 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyze(), analyzeTracks(), beginJob(), and clearVectors().

◆ ecalDetId_

std::vector<unsigned int> HcalHBHEMuonHighEtaAnalyzer::ecalDetId_
private

Definition at line 156 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyze(), analyzeTracks(), beginJob(), and clearVectors().

◆ ecalEnergy_

std::vector<double> HcalHBHEMuonHighEtaAnalyzer::ecalEnergy_
private

Definition at line 153 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyze(), analyzeMuon(), beginJob(), and clearVectors().

◆ ehcalDetId_

std::vector<unsigned int> HcalHBHEMuonHighEtaAnalyzer::ehcalDetId_
private

Definition at line 156 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyze(), analyzeTracks(), beginJob(), and clearVectors().

◆ emaxNearP_

std::vector<double> HcalHBHEMuonHighEtaAnalyzer::emaxNearP_
private

Definition at line 169 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyze(), analyzeTracks(), beginJob(), and clearVectors().

◆ emaxNearPThr_

const double HcalHBHEMuonHighEtaAnalyzer::emaxNearPThr_
private

Definition at line 104 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyzeTracks().

◆ endcapRecHitsHandle_

edm::Handle<EcalRecHitCollection> HcalHBHEMuonHighEtaAnalyzer::endcapRecHitsHandle_
private

Definition at line 143 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyze(), and analyzeTracks().

◆ energyMuon_

std::vector<double> HcalHBHEMuonHighEtaAnalyzer::energyMuon_
private

Definition at line 151 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyze(), analyzeMuon(), beginJob(), and clearVectors().

◆ etaGlob_

std::vector<double> HcalHBHEMuonHighEtaAnalyzer::etaGlob_
private

Definition at line 151 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyze(), analyzeMuon(), beginJob(), and clearVectors().

◆ etaMin_

const double HcalHBHEMuonHighEtaAnalyzer::etaMin_
private

Definition at line 104 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyzeHadron(), and analyzeMuon().

◆ eventNumber_

unsigned int HcalHBHEMuonHighEtaAnalyzer::eventNumber_
private

Definition at line 149 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyze(), beginJob(), and clearVectors().

◆ fileInCorr_

const std::string HcalHBHEMuonHighEtaAnalyzer::fileInCorr_
private

Definition at line 107 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by HcalHBHEMuonHighEtaAnalyzer().

◆ geo_

const CaloGeometry* HcalHBHEMuonHighEtaAnalyzer::geo_
private

Definition at line 134 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyzeHadron(), analyzeTracks(), and beginRun().

◆ getCharge_

const bool HcalHBHEMuonHighEtaAnalyzer::getCharge_
private

Definition at line 105 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyzeTracks(), and HcalHBHEMuonHighEtaAnalyzer().

◆ goodVertex_

unsigned int HcalHBHEMuonHighEtaAnalyzer::goodVertex_
private

Definition at line 149 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyze(), beginJob(), and clearVectors().

◆ hbhe_

edm::Handle<HBHERecHitCollection> HcalHBHEMuonHighEtaAnalyzer::hbhe_
private

Definition at line 144 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyze(), and analyzeTracks().

◆ hcal1x1Energy_

std::vector<double> HcalHBHEMuonHighEtaAnalyzer::hcal1x1Energy_
private

Definition at line 155 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyze(), analyzeTracks(), beginJob(), and clearVectors().

◆ hcal_ieta_

std::vector<int> HcalHBHEMuonHighEtaAnalyzer::hcal_ieta_
private

Definition at line 157 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyze(), analyzeTracks(), beginJob(), and clearVectors().

◆ hcal_iphi_

std::vector<int> HcalHBHEMuonHighEtaAnalyzer::hcal_iphi_
private

Definition at line 157 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyze(), analyzeTracks(), beginJob(), and clearVectors().

◆ hcalActiveLength_

std::vector<double> HcalHBHEMuonHighEtaAnalyzer::hcalActiveLength_
private

Definition at line 168 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyze(), analyzeTracks(), beginJob(), and clearVectors().

◆ hcalActiveLengthHot_

std::vector<double> HcalHBHEMuonHighEtaAnalyzer::hcalActiveLengthHot_
private

Definition at line 168 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyze(), analyzeTracks(), beginJob(), and clearVectors().

◆ hcalDepthActiveLength_

std::vector<double> HcalHBHEMuonHighEtaAnalyzer::hcalDepthActiveLength_[depthMax_]
private

Definition at line 159 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyze(), analyzeTracks(), beginJob(), and clearVectors().

◆ hcalDepthActiveLengthHot_

std::vector<double> HcalHBHEMuonHighEtaAnalyzer::hcalDepthActiveLengthHot_[depthMax_]
private

Definition at line 161 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyze(), analyzeTracks(), beginJob(), and clearVectors().

◆ hcalDepthChargeHot_

std::vector<double> HcalHBHEMuonHighEtaAnalyzer::hcalDepthChargeHot_[depthMax_]
private

Definition at line 162 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyze(), analyzeTracks(), beginJob(), and clearVectors().

◆ hcalDepthChargeHotBG_

std::vector<double> HcalHBHEMuonHighEtaAnalyzer::hcalDepthChargeHotBG_[depthMax_]
private

Definition at line 163 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyze(), analyzeTracks(), beginJob(), and clearVectors().

◆ hcalDepthEnergy_

std::vector<double> HcalHBHEMuonHighEtaAnalyzer::hcalDepthEnergy_[depthMax_]
private

Definition at line 158 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyze(), analyzeTracks(), beginJob(), and clearVectors().

◆ hcalDepthEnergyCorr_

std::vector<double> HcalHBHEMuonHighEtaAnalyzer::hcalDepthEnergyCorr_[depthMax_]
private

Definition at line 164 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyze(), analyzeTracks(), beginJob(), and clearVectors().

◆ hcalDepthEnergyHot_

std::vector<double> HcalHBHEMuonHighEtaAnalyzer::hcalDepthEnergyHot_[depthMax_]
private

Definition at line 160 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyze(), analyzeTracks(), beginJob(), and clearVectors().

◆ hcalDepthEnergyHotCorr_

std::vector<double> HcalHBHEMuonHighEtaAnalyzer::hcalDepthEnergyHotCorr_[depthMax_]
private

Definition at line 165 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyze(), analyzeTracks(), beginJob(), and clearVectors().

◆ hcalDepthMatch_

std::vector<bool> HcalHBHEMuonHighEtaAnalyzer::hcalDepthMatch_[depthMax_]
private

Definition at line 166 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyze(), analyzeTracks(), beginJob(), and clearVectors().

◆ hcalDepthMatchHot_

std::vector<bool> HcalHBHEMuonHighEtaAnalyzer::hcalDepthMatchHot_[depthMax_]
private

Definition at line 167 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyze(), analyzeTracks(), beginJob(), and clearVectors().

◆ hcalDetId_

std::vector<unsigned int> HcalHBHEMuonHighEtaAnalyzer::hcalDetId_
private

Definition at line 156 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyze(), analyzeTracks(), beginJob(), and clearVectors().

◆ hcalEnergy_

std::vector<double> HcalHBHEMuonHighEtaAnalyzer::hcalEnergy_
private

Definition at line 153 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyze(), analyzeMuon(), beginJob(), and clearVectors().

◆ hcalHot_

std::vector<bool> HcalHBHEMuonHighEtaAnalyzer::hcalHot_
private

Definition at line 154 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyze(), analyzeTracks(), beginJob(), and clearVectors().

◆ hdc_

const HcalDDDRecConstants* HcalHBHEMuonHighEtaAnalyzer::hdc_
private

Definition at line 132 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by activeLength(), analyzeTracks(), beginRun(), and goodCell().

◆ hoEnergy_

std::vector<double> HcalHBHEMuonHighEtaAnalyzer::hoEnergy_
private

Definition at line 153 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyze(), analyzeMuon(), beginJob(), and clearVectors().

◆ ignoreHECorr_

const bool HcalHBHEMuonHighEtaAnalyzer::ignoreHECorr_
private

Definition at line 108 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyzeTracks(), and HcalHBHEMuonHighEtaAnalyzer().

◆ isItPlan1_

const bool HcalHBHEMuonHighEtaAnalyzer::isItPlan1_
private

Definition at line 105 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyzeTracks(), depth16HE(), and HcalHBHEMuonHighEtaAnalyzer().

◆ isItPreRecHit_

const bool HcalHBHEMuonHighEtaAnalyzer::isItPreRecHit_
private

Definition at line 108 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyzeTracks(), depth16HE(), and HcalHBHEMuonHighEtaAnalyzer().

◆ isolationR03_

std::vector<double> HcalHBHEMuonHighEtaAnalyzer::isolationR03_
private

Definition at line 152 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyze(), analyzeMuon(), beginJob(), and clearVectors().

◆ isolationR04_

std::vector<double> HcalHBHEMuonHighEtaAnalyzer::isolationR04_
private

Definition at line 152 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyze(), analyzeMuon(), beginJob(), and clearVectors().

◆ kount_

int HcalHBHEMuonHighEtaAnalyzer::kount_
private

Definition at line 129 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyze(), and HcalHBHEMuonHighEtaAnalyzer().

◆ labelEBRecHit_

const edm::InputTag HcalHBHEMuonHighEtaAnalyzer::labelEBRecHit_
private

Definition at line 102 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by HcalHBHEMuonHighEtaAnalyzer().

◆ labelEERecHit_

const edm::InputTag HcalHBHEMuonHighEtaAnalyzer::labelEERecHit_
private

Definition at line 102 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by HcalHBHEMuonHighEtaAnalyzer().

◆ labelGenTrack_

const std::string HcalHBHEMuonHighEtaAnalyzer::labelGenTrack_
private

Definition at line 103 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by HcalHBHEMuonHighEtaAnalyzer().

◆ labelHBHERecHit_

const edm::InputTag HcalHBHEMuonHighEtaAnalyzer::labelHBHERecHit_
private

Definition at line 102 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by HcalHBHEMuonHighEtaAnalyzer().

◆ labelMuon_

const std::string HcalHBHEMuonHighEtaAnalyzer::labelMuon_
private

Definition at line 103 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by HcalHBHEMuonHighEtaAnalyzer().

◆ labelVtx_

const std::string HcalHBHEMuonHighEtaAnalyzer::labelVtx_
private

Definition at line 103 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by HcalHBHEMuonHighEtaAnalyzer().

◆ matchedId_

std::vector<bool> HcalHBHEMuonHighEtaAnalyzer::matchedId_
private

Definition at line 154 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyze(), analyzeTracks(), beginJob(), and clearVectors().

◆ maxDepth_

const int HcalHBHEMuonHighEtaAnalyzer::maxDepth_
private

Definition at line 109 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyzeTracks(), and beginJob().

◆ mediumMuon_

std::vector<bool> HcalHBHEMuonHighEtaAnalyzer::mediumMuon_
private

Definition at line 150 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyze(), analyzeMuon(), beginJob(), and clearVectors().

◆ mergedDepth_

bool HcalHBHEMuonHighEtaAnalyzer::mergedDepth_
private

Definition at line 128 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by activeLength(), beginRun(), and HcalHBHEMuonHighEtaAnalyzer().

◆ phiGlob_

std::vector<double> HcalHBHEMuonHighEtaAnalyzer::phiGlob_
private

Definition at line 151 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyze(), analyzeMuon(), beginJob(), and clearVectors().

◆ pMuon_

std::vector<double> HcalHBHEMuonHighEtaAnalyzer::pMuon_
private

Definition at line 151 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyze(), analyzeMuon(), beginJob(), and clearVectors().

◆ ptGlob_

std::vector<double> HcalHBHEMuonHighEtaAnalyzer::ptGlob_
private

Definition at line 151 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyze(), analyzeMuon(), beginJob(), and clearVectors().

◆ respCorrs_

HcalRespCorrs* HcalHBHEMuonHighEtaAnalyzer::respCorrs_
private

Definition at line 135 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by beginRun(), and respCorr().

◆ runNumber_

unsigned int HcalHBHEMuonHighEtaAnalyzer::runNumber_
private

Definition at line 149 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyze(), beginJob(), and clearVectors().

◆ selectionParameter_

spr::trackSelectionParameters HcalHBHEMuonHighEtaAnalyzer::selectionParameter_
private

Definition at line 130 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyzeTracks(), and HcalHBHEMuonHighEtaAnalyzer().

◆ sevlv_

const EcalSeverityLevelAlgo* HcalHBHEMuonHighEtaAnalyzer::sevlv_
private

Definition at line 138 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyze(), and analyzeTracks().

◆ theEcalChStatus_

const EcalChannelStatus* HcalHBHEMuonHighEtaAnalyzer::theEcalChStatus_
private

Definition at line 137 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyze(), and analyzeTracks().

◆ theHBHETopology_

const HcalTopology* HcalHBHEMuonHighEtaAnalyzer::theHBHETopology_
private

Definition at line 133 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyzeTracks(), beginRun(), and depth16HE().

◆ theTrackQuality_

const std::string HcalHBHEMuonHighEtaAnalyzer::theTrackQuality_
private

Definition at line 107 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyzeHadron(), and HcalHBHEMuonHighEtaAnalyzer().

◆ tok_chan_

const edm::ESGetToken<EcalChannelStatus, EcalChannelStatusRcd> HcalHBHEMuonHighEtaAnalyzer::tok_chan_
private

Definition at line 123 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyze().

◆ tok_dbservice_

const edm::ESGetToken<HcalDbService, HcalDbRecord> HcalHBHEMuonHighEtaAnalyzer::tok_dbservice_
private

Definition at line 126 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyze().

◆ tok_ddrec_

const edm::ESGetToken<HcalDDDRecConstants, HcalRecNumberingRecord> HcalHBHEMuonHighEtaAnalyzer::tok_ddrec_
private

Definition at line 118 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by beginRun().

◆ tok_EB_

const edm::EDGetTokenT<EcalRecHitCollection> HcalHBHEMuonHighEtaAnalyzer::tok_EB_
private

Definition at line 112 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyze().

◆ tok_EE_

const edm::EDGetTokenT<EcalRecHitCollection> HcalHBHEMuonHighEtaAnalyzer::tok_EE_
private

Definition at line 113 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyze().

◆ tok_genTrack_

const edm::EDGetTokenT<reco::TrackCollection> HcalHBHEMuonHighEtaAnalyzer::tok_genTrack_
private

Definition at line 116 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyzeHadron().

◆ tok_geom_

const edm::ESGetToken<CaloGeometry, CaloGeometryRecord> HcalHBHEMuonHighEtaAnalyzer::tok_geom_
private

Definition at line 121 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by beginRun().

◆ tok_HBHE_

const edm::EDGetTokenT<HBHERecHitCollection> HcalHBHEMuonHighEtaAnalyzer::tok_HBHE_
private

Definition at line 114 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyze().

◆ tok_htopo_

const edm::ESGetToken<HcalTopology, HcalRecNumberingRecord> HcalHBHEMuonHighEtaAnalyzer::tok_htopo_
private

Definition at line 119 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by beginRun().

◆ tok_magField_

const edm::ESGetToken<MagneticField, IdealMagneticFieldRecord> HcalHBHEMuonHighEtaAnalyzer::tok_magField_
private

Definition at line 122 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyze().

◆ tok_Muon_

const edm::EDGetTokenT<reco::MuonCollection> HcalHBHEMuonHighEtaAnalyzer::tok_Muon_
private

Definition at line 115 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyzeMuon().

◆ tok_respcorr_

const edm::ESGetToken<HcalRespCorrs, HcalRespCorrsRcd> HcalHBHEMuonHighEtaAnalyzer::tok_respcorr_
private

Definition at line 120 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by beginRun().

◆ tok_sevlv_

const edm::ESGetToken<EcalSeverityLevelAlgo, EcalSeverityLevelAlgoRcd> HcalHBHEMuonHighEtaAnalyzer::tok_sevlv_
private

Definition at line 124 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyze().

◆ tok_topo_

const edm::ESGetToken<CaloTopology, CaloTopologyRecord> HcalHBHEMuonHighEtaAnalyzer::tok_topo_
private

Definition at line 125 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyze().

◆ tok_Vtx_

const edm::EDGetTokenT<reco::VertexCollection> HcalHBHEMuonHighEtaAnalyzer::tok_Vtx_
private

Definition at line 111 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyze().

◆ trackDz_

std::vector<double> HcalHBHEMuonHighEtaAnalyzer::trackDz_
private

Definition at line 169 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyze(), beginJob(), clearVectors(), and fillTrackParameters().

◆ trackLayerCrossed_

std::vector<int> HcalHBHEMuonHighEtaAnalyzer::trackLayerCrossed_
private

Definition at line 170 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyze(), beginJob(), clearVectors(), and fillTrackParameters().

◆ trackMissedInnerHits_

std::vector<int> HcalHBHEMuonHighEtaAnalyzer::trackMissedInnerHits_
private

Definition at line 171 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyze(), beginJob(), clearVectors(), and fillTrackParameters().

◆ trackMissedOuterHits_

std::vector<int> HcalHBHEMuonHighEtaAnalyzer::trackMissedOuterHits_
private

Definition at line 171 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyze(), beginJob(), clearVectors(), and fillTrackParameters().

◆ trackOuterHit_

std::vector<int> HcalHBHEMuonHighEtaAnalyzer::trackOuterHit_
private

Definition at line 170 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyze(), beginJob(), clearVectors(), and fillTrackParameters().

◆ tree_

TTree* HcalHBHEMuonHighEtaAnalyzer::tree_
private

Definition at line 148 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ unCorrect_

const bool HcalHBHEMuonHighEtaAnalyzer::unCorrect_
private

Definition at line 105 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyzeTracks(), and HcalHBHEMuonHighEtaAnalyzer().

◆ useMyCorr_

bool HcalHBHEMuonHighEtaAnalyzer::useMyCorr_
private

Definition at line 128 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by HcalHBHEMuonHighEtaAnalyzer(), and respCorr().

◆ useRaw_

const int HcalHBHEMuonHighEtaAnalyzer::useRaw_
private

Definition at line 106 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyzeTracks(), and HcalHBHEMuonHighEtaAnalyzer().

◆ verbosity_

const int HcalHBHEMuonHighEtaAnalyzer::verbosity_
private

Definition at line 106 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyze(), analyzeMuon(), analyzeTracks(), beginRun(), and depth16HE().

◆ writeRespCorr_

const bool HcalHBHEMuonHighEtaAnalyzer::writeRespCorr_
private

Definition at line 108 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by beginRun().

Generated for CMSSW Reference Manual by doxygen 1.8.14