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

Public Member Functions

 HcalHBHEMuonAnalyzer (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
 
ESProxyIndex const * esGetTokenIndices (edm::Transition iTrans) const
 
std::vector< ESProxyIndex > 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
 
bool registeredToConsumeMany (TypeID const &, 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::ESRecordsToProxyIndices 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
 
void beginJob () override
 
void beginRun (edm::Run const &, edm::EventSetup const &) override
 
void clearVectors ()
 
int depth16HE (int ieta, int iphi)
 
void endRun (edm::Run const &, edm::EventSetup const &) override
 
double gainFactor (const HcalDbService *dbserv, const HcalDetId &id)
 
bool goodCell (const HcalDetId &hcid, const reco::Track *pTrack, const CaloGeometry *geo, const MagneticField *bField)
 
bool isGoodVertex (const reco::Vertex &vtx)
 
int matchId (const HcalDetId &, const HcalDetId &)
 
double respCorr (const DetId &id)
 

Private Attributes

std::vector< HcalDDDRecConstants::HcalActiveLength > actHB
 
std::vector< HcalDDDRecConstants::HcalActiveLength > actHE
 
std::vector< std::string > all_triggers_
 
unsigned int bxNumber_
 
double chiGlobal_
 
double chiTracker_
 
const bool collapseDepth_
 
std::map< DetId, double > corrValue_
 
double dxyTracker_
 
double dzTracker_
 
double ecal3x3Energy_
 
unsigned int ecalDetId_
 
double ecalEnergy_
 
unsigned int ehcalDetId_
 
double energyMuon_
 
double etaGlob_
 
unsigned int eventNumber_
 
const std::string fileInCorr_
 
const CaloGeometry * geo_
 
const bool getCharge_
 
int globalMuonHits_
 
bool globalTrack_
 
double globalTrckEta_
 
double globalTrckPhi_
 
double globalTrckPt_
 
unsigned int goodVertex_
 
double hcal1x1Energy_
 
int hcal_ieta_
 
int hcal_iphi_
 
double hcalActiveLength_
 
double hcalActiveLengthHot_
 
double hcalDepthActiveLength_ [depthMax_]
 
double hcalDepthActiveLengthHot_ [depthMax_]
 
double hcalDepthChargeHot_ [depthMax_]
 
double hcalDepthChargeHotBG_ [depthMax_]
 
double hcalDepthEnergy_ [depthMax_]
 
double hcalDepthEnergyCorr_ [depthMax_]
 
double hcalDepthEnergyHot_ [depthMax_]
 
double hcalDepthEnergyHotCorr_ [depthMax_]
 
bool hcalDepthMatch_ [depthMax_]
 
bool hcalDepthMatchHot_ [depthMax_]
 
unsigned int hcalDetId_
 
double hcalEnergy_
 
bool hcalHot_
 
const HcalDDDRecConstants * hdc_
 
HLTConfigProvider hltConfig_
 
std::vector< int > hltresults_
 
const edm::InputTag hlTriggerResults_
 
double hoEnergy_
 
const bool ignoreHECorr_
 
bool innerTrack_
 
double innerTracketa_
 
double innerTrackphi_
 
double innerTrackpt_
 
const bool isItPlan1_
 
const bool isItPreRecHit_
 
double isolationR03_
 
double isolationR04_
 
int kount_
 
const edm::InputTag labelEBRecHit_
 
const edm::InputTag labelEERecHit_
 
const edm::InputTag labelHBHERecHit_
 
const std::string labelMuon_
 
const std::string labelVtx_
 
unsigned int lumiNumber_
 
bool matchedId_
 
int matchedStat_
 
const int maxDepth_
 
bool mergedDepth_
 
const std::string modnam_
 
float muon_chi2LocalPosition_
 
bool muon_global_
 
bool muon_is_good_
 
bool muon_is_medium_
 
bool muon_is_tight_
 
float muon_segComp_
 
bool muon_tracker_
 
float muon_trkKink_
 
int numPixelLayers_
 
bool outerTrack_
 
double outerTrackChi_
 
double outerTrackEta_
 
int outerTrackHits_
 
double outerTrackPhi_
 
double outerTrackPt_
 
int outerTrackRHits_
 
double phiGlob_
 
const double pMinMuon_
 
double pMuon_
 
const std::string procnm_
 
double ptGlob_
 
HcalRespCorrs * respCorrs_
 
unsigned int runNumber_
 
const HcalTopology * theHBHETopology_
 
double tight_LongPara_
 
int tight_PixelHits_
 
double tight_TransImpara_
 
double tight_validFraction_
 
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::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< edm::TriggerResults > tok_trigRes_
 
const edm::EDGetTokenT< reco::VertexCollection > tok_Vtx_
 
int trackerLayer_
 
TTree * tree_
 
const std::vector< std::string > triggers_
 
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 ProductType , BranchType B = InEvent>
void consumesMany ()
 
void consumesMany (const TypeToGet &id)
 
template<BranchType B>
void consumesMany (const TypeToGet &id)
 
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 72 of file HcalHBHEMuonAnalyzer.cc.

Constructor & Destructor Documentation

◆ HcalHBHEMuonAnalyzer()

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

Definition at line 175 of file HcalHBHEMuonAnalyzer.cc.

References collapseDepth_, corrValue_, fileInCorr_, getCharge_, hlTriggerResults_, l1ctLayer2EG_cff::id, ignoreHECorr_, timingPdfMaker::infile, HLT_2022v15_cff::InputTag, isItPlan1_, isItPreRecHit_, kount_, TFileService::kSharedResource, labelEBRecHit_, labelEERecHit_, labelHBHERecHit_, labelMuon_, labelVtx_, mergedDepth_, modnam_, pMinMuon_, procnm_, unCorrect_, useMyCorr_, and useRaw_.

176  : hlTriggerResults_(iConfig.getParameter<edm::InputTag>("hlTriggerResults")),
177  labelEBRecHit_(iConfig.getParameter<edm::InputTag>("labelEBRecHit")),
178  labelEERecHit_(iConfig.getParameter<edm::InputTag>("labelEERecHit")),
179  labelHBHERecHit_(iConfig.getParameter<edm::InputTag>("labelHBHERecHit")),
180  labelVtx_(iConfig.getParameter<std::string>("labelVertex")),
181  labelMuon_(iConfig.getParameter<std::string>("labelMuon")),
182  fileInCorr_(iConfig.getUntrackedParameter<std::string>("fileInCorr", "")),
183  triggers_(iConfig.getParameter<std::vector<std::string>>("triggers")),
184  pMinMuon_(iConfig.getParameter<double>("pMinMuon")),
185  verbosity_(iConfig.getUntrackedParameter<int>("verbosity", 0)),
186  useRaw_(iConfig.getParameter<int>("useRaw")),
187  unCorrect_(iConfig.getParameter<bool>("unCorrect")),
188  collapseDepth_(iConfig.getParameter<bool>("collapseDepth")),
189  isItPlan1_(iConfig.getParameter<bool>("isItPlan1")),
190  ignoreHECorr_(iConfig.getUntrackedParameter<bool>("ignoreHECorr", false)),
191  isItPreRecHit_(iConfig.getUntrackedParameter<bool>("isItPreRecHit", false)),
192  getCharge_(iConfig.getParameter<bool>("getCharge")),
193  writeRespCorr_(iConfig.getUntrackedParameter<bool>("writeRespCorr", false)),
194  maxDepth_(iConfig.getUntrackedParameter<int>("maxDepth", 7)),
195  modnam_(iConfig.getUntrackedParameter<std::string>("moduleName", "")),
196  procnm_(iConfig.getUntrackedParameter<std::string>("processName", "")),
197  tok_trigRes_(consumes<edm::TriggerResults>(hlTriggerResults_)),
198  tok_Vtx_((modnam_.empty()) ? consumes<reco::VertexCollection>(labelVtx_)
199  : consumes<reco::VertexCollection>(edm::InputTag(modnam_, labelVtx_, procnm_))),
200  tok_EB_(consumes<EcalRecHitCollection>(labelEBRecHit_)),
201  tok_EE_(consumes<EcalRecHitCollection>(labelEERecHit_)),
202  tok_HBHE_(consumes<HBHERecHitCollection>(labelHBHERecHit_)),
203  tok_Muon_((modnam_.empty()) ? consumes<reco::MuonCollection>(labelMuon_)
204  : consumes<reco::MuonCollection>(edm::InputTag(modnam_, labelMuon_, procnm_))),
205  tok_ddrec_(esConsumes<HcalDDDRecConstants, HcalRecNumberingRecord, edm::Transition::BeginRun>()),
206  tok_htopo_(esConsumes<HcalTopology, HcalRecNumberingRecord, edm::Transition::BeginRun>()),
207  tok_respcorr_(esConsumes<HcalRespCorrs, HcalRespCorrsRcd, edm::Transition::BeginRun>()),
208  tok_geom_(esConsumes<CaloGeometry, CaloGeometryRecord, edm::Transition::BeginRun>()),
209  tok_magField_(esConsumes<MagneticField, IdealMagneticFieldRecord>()),
210  tok_chan_(esConsumes<EcalChannelStatus, EcalChannelStatusRcd>()),
211  tok_sevlv_(esConsumes<EcalSeverityLevelAlgo, EcalSeverityLevelAlgoRcd>()),
212  tok_topo_(esConsumes<CaloTopology, CaloTopologyRecord>()),
213  tok_dbservice_(esConsumes<HcalDbService, HcalDbRecord>()),
214  hdc_(nullptr),
215  theHBHETopology_(nullptr),
216  respCorrs_(nullptr) {
217  usesResource(TFileService::kSharedResource);
218  //now do what ever initialization is needed
219  kount_ = 0;
220  mergedDepth_ = (!isItPreRecHit_) || (collapseDepth_);
221 
222  if (modnam_.empty()) {
223  edm::LogVerbatim("HBHEMuon") << "Labels used: Trig " << hlTriggerResults_ << " Vtx " << labelVtx_ << " EB "
224  << labelEBRecHit_ << " EE " << labelEERecHit_ << " HBHE " << labelHBHERecHit_ << " MU "
225  << labelMuon_;
226  } else {
227  edm::LogVerbatim("HBHEMuon") << "Labels used Trig " << hlTriggerResults_ << "\n Vtx "
228  << edm::InputTag(modnam_, labelVtx_, procnm_) << "\n EB " << labelEBRecHit_
229  << "\n EE " << labelEERecHit_ << "\n HBHE " << labelHBHERecHit_ << "\n MU "
230  << edm::InputTag(modnam_, labelMuon_, procnm_);
231  }
232 
233  if (!fileInCorr_.empty()) {
234  std::ifstream infile(fileInCorr_.c_str());
235  if (infile.is_open()) {
236  while (true) {
237  unsigned int id;
238  double cfac;
239  infile >> id >> cfac;
240  if (!infile.good())
241  break;
242  corrValue_[DetId(id)] = cfac;
243  }
244  infile.close();
245  }
246  }
247  useMyCorr_ = (!corrValue_.empty());
248  edm::LogVerbatim("HBHEMuon") << "Flags used: UseRaw " << useRaw_ << " GetCharge " << getCharge_ << " UnCorrect "
249  << unCorrect_ << " IgnoreHECorr " << ignoreHECorr_ << " CollapseDepth " << collapseDepth_
250  << ":" << mergedDepth_ << " IsItPlan1 " << isItPlan1_ << " IsItPreRecHit "
251  << isItPreRecHit_ << " UseMyCorr " << useMyCorr_ << " pMinMuon " << pMinMuon_;
252 }
TFileService::kSharedResource
static const std::string kSharedResource
Definition: TFileService.h:76
edm::LogVerbatim
Log< level::Info, true > LogVerbatim
Definition: MessageLogger.h:128
HcalHBHEMuonAnalyzer::labelMuon_
const std::string labelMuon_
Definition: HcalHBHEMuonAnalyzer.cc:96
HcalHBHEMuonAnalyzer::theHBHETopology_
const HcalTopology * theHBHETopology_
Definition: HcalHBHEMuonAnalyzer.cc:124
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
Definition: ParameterSet.h:303
HcalHBHEMuonAnalyzer::tok_respcorr_
const edm::ESGetToken< HcalRespCorrs, HcalRespCorrsRcd > tok_respcorr_
Definition: HcalHBHEMuonAnalyzer.cc:115
HcalHBHEMuonAnalyzer::tok_htopo_
const edm::ESGetToken< HcalTopology, HcalRecNumberingRecord > tok_htopo_
Definition: HcalHBHEMuonAnalyzer.cc:114
HcalHBHEMuonAnalyzer::tok_Muon_
const edm::EDGetTokenT< reco::MuonCollection > tok_Muon_
Definition: HcalHBHEMuonAnalyzer.cc:111
HcalHBHEMuonAnalyzer::tok_sevlv_
const edm::ESGetToken< EcalSeverityLevelAlgo, EcalSeverityLevelAlgoRcd > tok_sevlv_
Definition: HcalHBHEMuonAnalyzer.cc:119
HLT_2022v15_cff.InputTag
InputTag
Definition: HLT_2022v15_cff.py:60426
HcalHBHEMuonAnalyzer::hlTriggerResults_
const edm::InputTag hlTriggerResults_
Definition: HcalHBHEMuonAnalyzer.cc:94
HcalHBHEMuonAnalyzer::tok_topo_
const edm::ESGetToken< CaloTopology, CaloTopologyRecord > tok_topo_
Definition: HcalHBHEMuonAnalyzer.cc:120
l1ctLayer2EG_cff.id
id
Definition: l1ctLayer2EG_cff.py:65
l1t::VertexCollection
std::vector< Vertex > VertexCollection
Definition: Vertex.h:31
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
reco::MuonCollection
std::vector< Muon > MuonCollection
collection of Muon objects
Definition: MuonFwd.h:9
edm::EDConsumerBase::consumes
EDGetTokenT< ProductType > consumes(edm::InputTag const &tag)
Definition: EDConsumerBase.h:161
HcalHBHEMuonAnalyzer::labelEERecHit_
const edm::InputTag labelEERecHit_
Definition: HcalHBHEMuonAnalyzer.cc:95
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
HcalHBHEMuonAnalyzer::labelEBRecHit_
const edm::InputTag labelEBRecHit_
Definition: HcalHBHEMuonAnalyzer.cc:95
HcalHBHEMuonAnalyzer::tok_dbservice_
const edm::ESGetToken< HcalDbService, HcalDbRecord > tok_dbservice_
Definition: HcalHBHEMuonAnalyzer.cc:121
HcalHBHEMuonAnalyzer::tok_HBHE_
const edm::EDGetTokenT< HBHERecHitCollection > tok_HBHE_
Definition: HcalHBHEMuonAnalyzer.cc:110
edm::Transition
Transition
Definition: Transition.h:12
HcalHBHEMuonAnalyzer::labelVtx_
const std::string labelVtx_
Definition: HcalHBHEMuonAnalyzer.cc:96
HcalHBHEMuonAnalyzer::labelHBHERecHit_
const edm::InputTag labelHBHERecHit_
Definition: HcalHBHEMuonAnalyzer.cc:95
DetId
Definition: DetId.h:17
HcalHBHEMuonAnalyzer::tok_trigRes_
const edm::EDGetTokenT< edm::TriggerResults > tok_trigRes_
Definition: HcalHBHEMuonAnalyzer.cc:106
HcalHBHEMuonAnalyzer::corrValue_
std::map< DetId, double > corrValue_
Definition: HcalHBHEMuonAnalyzer.cc:171
HcalHBHEMuonAnalyzer::tok_geom_
const edm::ESGetToken< CaloGeometry, CaloGeometryRecord > tok_geom_
Definition: HcalHBHEMuonAnalyzer.cc:116
HcalHBHEMuonAnalyzer::tok_EE_
const edm::EDGetTokenT< EcalRecHitCollection > tok_EE_
Definition: HcalHBHEMuonAnalyzer.cc:109
HcalHBHEMuonAnalyzer::tok_EB_
const edm::EDGetTokenT< EcalRecHitCollection > tok_EB_
Definition: HcalHBHEMuonAnalyzer.cc:108
HcalHBHEMuonAnalyzer::hdc_
const HcalDDDRecConstants * hdc_
Definition: HcalHBHEMuonAnalyzer.cc:123
reco
fixed size matrix
Definition: AlignmentAlgorithmBase.h:46
edm
HLT enums.
Definition: AlignableModifier.h:19
HcalHBHEMuonAnalyzer::tok_chan_
const edm::ESGetToken< EcalChannelStatus, EcalChannelStatusRcd > tok_chan_
Definition: HcalHBHEMuonAnalyzer.cc:118
HcalHBHEMuonAnalyzer::tok_Vtx_
const edm::EDGetTokenT< reco::VertexCollection > tok_Vtx_
Definition: HcalHBHEMuonAnalyzer.cc:107
HcalHBHEMuonAnalyzer::tok_magField_
const edm::ESGetToken< MagneticField, IdealMagneticFieldRecord > tok_magField_
Definition: HcalHBHEMuonAnalyzer.cc:117
HcalHBHEMuonAnalyzer::fileInCorr_
const std::string fileInCorr_
Definition: HcalHBHEMuonAnalyzer.cc:96
HcalHBHEMuonAnalyzer::tok_ddrec_
const edm::ESGetToken< HcalDDDRecConstants, HcalRecNumberingRecord > tok_ddrec_
Definition: HcalHBHEMuonAnalyzer.cc:113
HcalHBHEMuonAnalyzer::triggers_
const std::vector< std::string > triggers_
Definition: HcalHBHEMuonAnalyzer.cc:97

Member Function Documentation

◆ activeLength()

double HcalHBHEMuonAnalyzer::activeLength ( const DetId &  id_)
private

Definition at line 1176 of file HcalHBHEMuonAnalyzer.cc.

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

Referenced by analyze(), and beginRun().

1176  {
1177  HcalDetId id(id_);
1178  int ieta = id.ietaAbs();
1179  int zside = id.zside();
1180  int iphi = id.iphi();
1181  std::vector<int> dpths;
1182  if (mergedDepth_) {
1183  std::vector<HcalDetId> ids;
1184  hdc_->unmergeDepthDetId(id, ids);
1185  for (auto idh : ids)
1186  dpths.emplace_back(idh.depth());
1187  } else {
1188  dpths.emplace_back(id.depth());
1189  }
1190  double lx(0);
1191  if (id.subdet() == HcalBarrel) {
1192  for (unsigned int i = 0; i < actHB.size(); ++i) {
1193  if ((ieta == actHB[i].ieta) && (zside == actHB[i].zside) &&
1194  (std::find(dpths.begin(), dpths.end(), actHB[i].depth) != dpths.end()) &&
1195  (std::find(actHB[i].iphis.begin(), actHB[i].iphis.end(), iphi) != actHB[i].iphis.end())) {
1196  lx += actHB[i].thick;
1197  }
1198  }
1199  } else {
1200  for (unsigned int i = 0; i < actHE.size(); ++i) {
1201  if ((ieta == actHE[i].ieta) && (zside == actHE[i].zside) &&
1202  (std::find(dpths.begin(), dpths.end(), actHE[i].depth) != dpths.end()) &&
1203  (std::find(actHE[i].iphis.begin(), actHE[i].iphis.end(), iphi) != actHE[i].iphis.end())) {
1204  lx += actHE[i].thick;
1205  }
1206  }
1207  }
1208  return lx;
1209 }
mps_fire.i
i
Definition: mps_fire.py:429
ecaldqm::zside
int zside(DetId const &)
Definition: EcalDQMCommonUtils.cc:189
HcalHBHEMuonAnalyzer::actHE
std::vector< HcalDDDRecConstants::HcalActiveLength > actHE
Definition: HcalHBHEMuonAnalyzer.cc:170
l1ctLayer2EG_cff.id
id
Definition: l1ctLayer2EG_cff.py:65
spr::find
void find(edm::Handle< EcalRecHitCollection > &hits, DetId thisDet, std::vector< EcalRecHitCollection::const_iterator > &hit, bool debug=false)
Definition: FindCaloHit.cc:19
HcalHBHEMuonAnalyzer::actHB
std::vector< HcalDDDRecConstants::HcalActiveLength > actHB
Definition: HcalHBHEMuonAnalyzer.cc:170
HcalHBHEMuonAnalyzer::hdc_
const HcalDDDRecConstants * hdc_
Definition: HcalHBHEMuonAnalyzer.cc:123
HcalDDDRecConstants::unmergeDepthDetId
void unmergeDepthDetId(const HcalDetId &id, std::vector< HcalDetId > &ids) const
Definition: HcalDDDRecConstants.cc:677

◆ analyze()

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

Implements edm::one::EDAnalyzerBase.

Definition at line 259 of file HcalHBHEMuonAnalyzer.cc.

References funct::abs(), accept(), edm::HLTGlobalStatus::accept(), activeLength(), all_triggers_, Calorimetry_cff::bField, bxNumber_, RPCNoise_example::check, chg, chiGlobal_, chiTracker_, clearVectors(), collapseDepth_, submitPVValidationJobs::conditions, alignCSCRings::corr, LEDCalibrationChannels::depth, depth16HE(), depthMax_, spr::propagatedTrackID::detIdECAL, spr::propagatedTrackID::detIdEHCAL, spr::propagatedTrackID::detIdHCAL, dxyTracker_, dzTracker_, ecal3x3Energy_, ecalDetId_, ecalEnergy_, spr::eECALmatrix(), ehcalDetId_, spr::eHCALmatrix(), spr::energyHCALCell(), energyMuon_, etaGlob_, eventNumber_, PedestalClient_cfi::gain, gainFactor(), geo_, getCharge_, edm::EventSetup::getData(), globalMuonHits_, cutBasedMuonId_MuonPOG_V0_cff::globalTrack, globalTrack_, globalTrckEta_, globalTrckPhi_, globalTrckPt_, goodCell(), goodVertex_, photonIsolationHIProducer_cfi::hbhe, hcal1x1Energy_, hcal_ieta_, hcal_iphi_, hcalActiveLength_, hcalActiveLengthHot_, HcalBarrel, hcalDepthActiveLength_, hcalDepthActiveLengthHot_, hcalDepthChargeHot_, hcalDepthChargeHotBG_, hcalDepthEnergy_, hcalDepthEnergyCorr_, hcalDepthEnergyHot_, hcalDepthEnergyHotCorr_, hcalDepthMatch_, hcalDepthMatchHot_, hcalDetId_, HcalEndcap, hcalEnergy_, hcalHot_, hdc_, hltresults_, hoEnergy_, mps_fire::i, LEDCalibrationChannels::ieta, HcalDetId::ieta(), iEvent, ignoreHECorr_, info(), susybsm::HSCParticleType::innerTrack, innerTrack_, innerTracketa_, innerTrackphi_, innerTrackpt_, LEDCalibrationChannels::iphi, HcalDetId::iphi(), isGoodVertex(), isItPlan1_, isItPreRecHit_, muon::isMediumMuon(), muons2muons_cfi::isolationR03, isolationR03_, muons2muons_cfi::isolationR04, isolationR04_, muon::isTightMuon(), edm::HandleBase::isValid(), dqmdumpme::k, kount_, lumiNumber_, matchedId_, matchedStat_, matchId(), SiStripPI::max, maxDepth_, HcalDDDRecConstants::mergedDepthDetId(), muon_chi2LocalPosition_, muon_global_, muon_is_good_, muon_is_medium_, muon_is_tight_, muon_segComp_, muon_tracker_, muon_trkKink_, numPixelLayers_, spr::propagatedTrackID::okECAL, spr::propagatedTrackID::okHCAL, outerTrack_, outerTrackChi_, outerTrackEta_, outerTrackHits_, outerTrackPhi_, outerTrackPt_, outerTrackRHits_, phiGlob_, pMinMuon_, pMuon_, edm::Handle< T >::product(), spr::propagateCALO(), spr::propagateHCALBack(), ptGlob_, DetId::rawId(), respCorr(), runNumber_, edm::second(), muon::segmentCompatibility(), edm::HLTGlobalStatus::size(), HcalDetId::subdet(), theHBHETopology_, tight_LongPara_, tight_PixelHits_, tight_TransImpara_, tight_validFraction_, tok_chan_, tok_dbservice_, tok_EB_, tok_EE_, tok_HBHE_, tok_magField_, tok_Muon_, tok_sevlv_, tok_topo_, tok_trigRes_, tok_Vtx_, trackerLayer_, tree_, edm::TriggerNames::triggerIndex(), unCorrect_, useRaw_, HcalTopology::validHcal(), verbosity_, and extraflags_cff::vtx.

259  {
260  ++kount_;
261  clearVectors();
262  std::vector<bool> muon_is_good, muon_global, muon_tracker;
263  std::vector<bool> muon_is_tight, muon_is_medium;
264  std::vector<double> ptGlob, etaGlob, phiGlob, energyMuon, pMuon;
265  std::vector<float> muon_trkKink, muon_chi2LocalPosition, muon_segComp;
266  std::vector<int> trackerLayer, numPixelLayers, tight_PixelHits;
267  std::vector<bool> innerTrack, outerTrack, globalTrack;
268  std::vector<double> chiTracker, dxyTracker, dzTracker;
269  std::vector<double> innerTrackpt, innerTracketa, innerTrackphi;
270  std::vector<double> tight_validFraction, outerTrackChi;
271  std::vector<double> outerTrackPt, outerTrackEta, outerTrackPhi;
272  std::vector<int> outerTrackHits, outerTrackRHits;
273  std::vector<double> globalTrckPt, globalTrckEta, globalTrckPhi;
274  std::vector<int> globalMuonHits, matchedStat;
275  std::vector<double> chiGlobal, tight_LongPara, tight_TransImpara;
276  std::vector<double> isolationR04, isolationR03;
277  std::vector<double> ecalEnergy, hcalEnergy, hoEnergy;
278  std::vector<bool> matchedId, hcalHot;
279  std::vector<double> ecal3x3Energy, hcal1x1Energy;
280  std::vector<unsigned int> ecalDetId, hcalDetId, ehcalDetId;
281  std::vector<int> hcal_ieta, hcal_iphi;
282  std::vector<double> hcalDepthEnergy[depthMax_];
283  std::vector<double> hcalDepthActiveLength[depthMax_];
284  std::vector<double> hcalDepthEnergyHot[depthMax_];
285  std::vector<double> hcalDepthActiveLengthHot[depthMax_];
286  std::vector<double> hcalDepthChargeHot[depthMax_];
287  std::vector<double> hcalDepthChargeHotBG[depthMax_];
288  std::vector<double> hcalDepthEnergyCorr[depthMax_];
289  std::vector<double> hcalDepthEnergyHotCorr[depthMax_];
290  std::vector<bool> hcalDepthMatch[depthMax_];
291  std::vector<bool> hcalDepthMatchHot[depthMax_];
292  std::vector<double> hcalActiveLength, hcalActiveLengthHot;
293  runNumber_ = iEvent.id().run();
294  eventNumber_ = iEvent.id().event();
295  lumiNumber_ = iEvent.id().luminosityBlock();
296  bxNumber_ = iEvent.bunchCrossing();
297 #ifdef EDM_ML_DEBUG
298  edm::LogVerbatim("HBHEMuon") << "Run " << runNumber_ << " Event " << eventNumber_ << " Lumi " << lumiNumber_ << " BX "
299  << bxNumber_ << std::endl;
300 #endif
301  const edm::Handle<edm::TriggerResults>& _Triggers = iEvent.getHandle(tok_trigRes_);
302 #ifdef EDM_ML_DEBUG
303  if ((verbosity_ / 10000) % 10 > 0)
304  edm::LogVerbatim("HBHEMuon") << "Size of all triggers " << all_triggers_.size();
305 #endif
306  int Ntriggers = static_cast<int>(all_triggers_.size());
307 #ifdef EDM_ML_DEBUG
308  if ((verbosity_ / 10000) % 10 > 0)
309  edm::LogVerbatim("HBHEMuon") << "Size of HLT MENU: " << _Triggers->size();
310 #endif
311  if (_Triggers.isValid()) {
312  const edm::TriggerNames& triggerNames_ = iEvent.triggerNames(*_Triggers);
313  std::vector<int> index;
314  for (int i = 0; i < Ntriggers; i++) {
315  index.push_back(triggerNames_.triggerIndex(all_triggers_[i]));
316  int triggerSize = static_cast<int>(_Triggers->size());
317 #ifdef EDM_ML_DEBUG
318  if ((verbosity_ / 10000) % 10 > 0)
319  edm::LogVerbatim("HBHEMuon") << "outside loop " << index[i] << "\ntriggerSize " << triggerSize;
320 #endif
321  if (index[i] < triggerSize) {
322  hltresults_.push_back(_Triggers->accept(index[i]));
323 #ifdef EDM_ML_DEBUG
324  if ((verbosity_ / 10000) % 10 > 0)
325  edm::LogVerbatim("HBHEMuon") << "Trigger_info " << triggerSize << " triggerSize " << index[i]
326  << " trigger_index " << hltresults_.at(i) << " hltresult";
327 #endif
328  } else {
329  if ((verbosity_ / 10000) % 10 > 0)
330  edm::LogVerbatim("HBHEMuon") << "Requested HLT path \""
331  << "\" does not exist";
332  }
333  }
334  }
335 
336  // get handles to calogeometry and calotopology
337  const MagneticField* bField = &iSetup.getData(tok_magField_);
338  const EcalChannelStatus* theEcalChStatus = &iSetup.getData(tok_chan_);
339  const EcalSeverityLevelAlgo* sevlv = &iSetup.getData(tok_sevlv_);
340  const CaloTopology* caloTopology = &iSetup.getData(tok_topo_);
341  const HcalDbService* conditions = &iSetup.getData(tok_dbservice_);
342 
343  // Relevant blocks from iEvent
344  const edm::Handle<reco::VertexCollection>& vtx = iEvent.getHandle(tok_Vtx_);
345 
346  edm::Handle<EcalRecHitCollection> barrelRecHitsHandle = iEvent.getHandle(tok_EB_);
347  edm::Handle<EcalRecHitCollection> endcapRecHitsHandle = iEvent.getHandle(tok_EE_);
348 
349  edm::Handle<HBHERecHitCollection> hbhe = iEvent.getHandle(tok_HBHE_);
350 
351  const edm::Handle<reco::MuonCollection>& _Muon = iEvent.getHandle(tok_Muon_);
352 
353  // require a good vertex
354  math::XYZPoint pvx;
355  goodVertex_ = 0;
356  if (!vtx.isValid()) {
357 #ifdef EDM_ML_DEBUG
358  edm::LogVerbatim("HBHEMuon") << "No Good Vertex found == Reject";
359 #endif
360  return;
361  }
362  reco::VertexCollection::const_iterator firstGoodVertex = vtx->end();
363  for (reco::VertexCollection::const_iterator it = vtx->begin(); it != vtx->end(); it++) {
364  if (isGoodVertex(*it)) {
365  if (firstGoodVertex == vtx->end())
366  firstGoodVertex = it;
367  ++goodVertex_;
368  }
369  }
370  if (firstGoodVertex != vtx->end())
371  pvx = firstGoodVertex->position();
372 
373  bool accept(false);
374  if (_Muon.isValid() && barrelRecHitsHandle.isValid() && endcapRecHitsHandle.isValid() && hbhe.isValid()) {
375  for (const auto& RecMuon : (*(_Muon.product()))) {
376  muon_is_good.push_back(RecMuon.isPFMuon());
377  muon_global.push_back(RecMuon.isGlobalMuon());
378  muon_tracker.push_back(RecMuon.isTrackerMuon());
379  ptGlob.push_back(RecMuon.pt());
380  etaGlob.push_back(RecMuon.eta());
381  phiGlob.push_back(RecMuon.phi());
382  energyMuon.push_back(RecMuon.energy());
383  pMuon.push_back(RecMuon.p());
384 #ifdef EDM_ML_DEBUG
385  edm::LogVerbatim("HBHEMuon") << "Energy:" << RecMuon.energy() << " P:" << RecMuon.p();
386 #endif
387  muon_is_tight.push_back(muon::isTightMuon(RecMuon, *firstGoodVertex));
388  muon_is_medium.push_back(muon::isMediumMuon(RecMuon));
389  muon_trkKink.push_back(RecMuon.combinedQuality().trkKink);
390  muon_chi2LocalPosition.push_back(RecMuon.combinedQuality().chi2LocalPosition);
391  muon_segComp.push_back(muon::segmentCompatibility(RecMuon));
392  // acessing tracker hits info
393  if (RecMuon.track().isNonnull()) {
394  trackerLayer.push_back(RecMuon.track()->hitPattern().trackerLayersWithMeasurement());
395  } else {
396  trackerLayer.push_back(-1);
397  }
398  if (RecMuon.innerTrack().isNonnull()) {
399  innerTrack.push_back(true);
400  numPixelLayers.push_back(RecMuon.innerTrack()->hitPattern().pixelLayersWithMeasurement());
401  chiTracker.push_back(RecMuon.innerTrack()->normalizedChi2());
402  dxyTracker.push_back(fabs(RecMuon.innerTrack()->dxy(pvx)));
403  dzTracker.push_back(fabs(RecMuon.innerTrack()->dz(pvx)));
404  innerTrackpt.push_back(RecMuon.innerTrack()->pt());
405  innerTracketa.push_back(RecMuon.innerTrack()->eta());
406  innerTrackphi.push_back(RecMuon.innerTrack()->phi());
407  tight_PixelHits.push_back(RecMuon.innerTrack()->hitPattern().numberOfValidPixelHits());
408  tight_validFraction.push_back(RecMuon.innerTrack()->validFraction());
409  } else {
410  innerTrack.push_back(false);
411  numPixelLayers.push_back(0);
412  chiTracker.push_back(0);
413  dxyTracker.push_back(0);
414  dzTracker.push_back(0);
415  innerTrackpt.push_back(0);
416  innerTracketa.push_back(0);
417  innerTrackphi.push_back(0);
418  tight_PixelHits.push_back(0);
419  tight_validFraction.push_back(-99);
420  }
421  // outer track info
422  if (RecMuon.outerTrack().isNonnull()) {
423  outerTrack.push_back(true);
424  outerTrackPt.push_back(RecMuon.outerTrack()->pt());
425  outerTrackEta.push_back(RecMuon.outerTrack()->eta());
426  outerTrackPhi.push_back(RecMuon.outerTrack()->phi());
427  outerTrackChi.push_back(RecMuon.outerTrack()->normalizedChi2());
428  outerTrackHits.push_back(RecMuon.outerTrack()->numberOfValidHits());
429  outerTrackRHits.push_back(RecMuon.outerTrack()->recHitsSize());
430  } else {
431  outerTrack.push_back(false);
432  outerTrackPt.push_back(0);
433  outerTrackEta.push_back(0);
434  outerTrackPhi.push_back(0);
435  outerTrackChi.push_back(0);
436  outerTrackHits.push_back(0);
437  outerTrackRHits.push_back(0);
438  }
439  // Tight Muon cuts
440  if (RecMuon.globalTrack().isNonnull()) {
441  globalTrack.push_back(true);
442  chiGlobal.push_back(RecMuon.globalTrack()->normalizedChi2());
443  globalMuonHits.push_back(RecMuon.globalTrack()->hitPattern().numberOfValidMuonHits());
444  matchedStat.push_back(RecMuon.numberOfMatchedStations());
445  globalTrckPt.push_back(RecMuon.globalTrack()->pt());
446  globalTrckEta.push_back(RecMuon.globalTrack()->eta());
447  globalTrckPhi.push_back(RecMuon.globalTrack()->phi());
448  tight_TransImpara.push_back(fabs(RecMuon.muonBestTrack()->dxy(pvx)));
449  tight_LongPara.push_back(fabs(RecMuon.muonBestTrack()->dz(pvx)));
450  } else {
451  globalTrack.push_back(false);
452  chiGlobal.push_back(0);
453  globalMuonHits.push_back(0);
454  matchedStat.push_back(0);
455  globalTrckPt.push_back(0);
456  globalTrckEta.push_back(0);
457  globalTrckPhi.push_back(0);
458  tight_TransImpara.push_back(0);
459  tight_LongPara.push_back(0);
460  }
461 
462  isolationR04.push_back(
463  ((RecMuon.pfIsolationR04().sumChargedHadronPt +
464  std::max(0.,
465  RecMuon.pfIsolationR04().sumNeutralHadronEt + RecMuon.pfIsolationR04().sumPhotonEt -
466  (0.5 * RecMuon.pfIsolationR04().sumPUPt))) /
467  RecMuon.pt()));
468 
469  isolationR03.push_back(
470  ((RecMuon.pfIsolationR03().sumChargedHadronPt +
471  std::max(0.,
472  RecMuon.pfIsolationR03().sumNeutralHadronEt + RecMuon.pfIsolationR03().sumPhotonEt -
473  (0.5 * RecMuon.pfIsolationR03().sumPUPt))) /
474  RecMuon.pt()));
475 
476  ecalEnergy.push_back(RecMuon.calEnergy().emS9);
477  hcalEnergy.push_back(RecMuon.calEnergy().hadS9);
478  hoEnergy.push_back(RecMuon.calEnergy().hoS9);
479 
480  double eEcal(0), eHcal(0), activeLengthTot(0), activeLengthHotTot(0);
481  double eHcalDepth[depthMax_], eHcalDepthHot[depthMax_];
482  double eHcalDepthC[depthMax_], eHcalDepthHotC[depthMax_];
483  double cHcalDepthHot[depthMax_], cHcalDepthHotBG[depthMax_];
484  double activeL[depthMax_], activeHotL[depthMax_];
485  bool matchDepth[depthMax_], matchDepthHot[depthMax_];
486  HcalDetId eHcalDetId[depthMax_];
487  unsigned int isHot(0);
488  bool tmpmatch(false);
489  int ieta(-1000), iphi(-1000);
490  for (int i = 0; i < depthMax_; ++i) {
491  eHcalDepth[i] = eHcalDepthHot[i] = 0;
492  eHcalDepthC[i] = eHcalDepthHotC[i] = 0;
493  cHcalDepthHot[i] = cHcalDepthHotBG[i] = 0;
494  activeL[i] = activeHotL[i] = 0;
495  matchDepth[i] = matchDepthHot[i] = true;
496  }
497  if (RecMuon.innerTrack().isNonnull()) {
498  const reco::Track* pTrack = (RecMuon.innerTrack()).get();
499  spr::propagatedTrackID trackID = spr::propagateCALO(pTrack, geo_, bField, (((verbosity_ / 100) % 10 > 0)));
500  if ((RecMuon.p() > pMinMuon_) && (trackID.okHCAL))
501  accept = true;
502 
503  ecalDetId.push_back((trackID.detIdECAL)());
504  hcalDetId.push_back((trackID.detIdHCAL)());
505  ehcalDetId.push_back((trackID.detIdEHCAL)());
506 
507  HcalDetId check;
508  std::pair<bool, HcalDetId> info = spr::propagateHCALBack(pTrack, geo_, bField, (((verbosity_ / 100) % 10 > 0)));
509  if (info.first) {
510  check = info.second;
511  }
512 
513  bool okE = trackID.okECAL;
514  if (okE) {
515  const DetId isoCell(trackID.detIdECAL);
516  std::pair<double, bool> e3x3 = spr::eECALmatrix(isoCell,
517  barrelRecHitsHandle,
518  endcapRecHitsHandle,
519  *theEcalChStatus,
520  geo_,
521  caloTopology,
522  sevlv,
523  1,
524  1,
525  -100.0,
526  -100.0,
527  -500.0,
528  500.0,
529  false);
530  eEcal = e3x3.first;
531  okE = e3x3.second;
532  }
533 #ifdef EDM_ML_DEBUG
534  edm::LogVerbatim("HBHEMuon") << "Propagate Track to ECAL: " << okE << ":" << trackID.okECAL << " E " << eEcal;
535 #endif
536 
537  if (trackID.okHCAL) {
538  DetId closestCell(trackID.detIdHCAL);
539  HcalDetId hcidt(closestCell.rawId());
540  if ((hcidt.ieta() == check.ieta()) && (hcidt.iphi() == check.iphi()))
541  tmpmatch = true;
542 #ifdef EDM_ML_DEBUG
543  edm::LogVerbatim("HBHEMuon") << "Front " << hcidt << " Back " << info.first << ":" << check << " Match "
544  << tmpmatch;
545 #endif
546 
547  HcalSubdetector subdet = hcidt.subdet();
548  ieta = hcidt.ieta();
549  iphi = hcidt.iphi();
550  bool hborhe = (std::abs(ieta) == 16);
551 
552  eHcal = spr::eHCALmatrix(theHBHETopology_,
553  closestCell,
554  hbhe,
555  0,
556  0,
557  false,
558  true,
559  -100.0,
560  -100.0,
561  -100.0,
562  -100.0,
563  -500.,
564  500.,
565  useRaw_);
566  std::vector<std::pair<double, int>> ehdepth;
567  spr::energyHCALCell((HcalDetId)closestCell,
568  hbhe,
569  ehdepth,
570  maxDepth_,
571  -100.0,
572  -100.0,
573  -100.0,
574  -100.0,
575  -500.0,
576  500.0,
577  useRaw_,
578  depth16HE(ieta, iphi),
579  (((verbosity_ / 1000) % 10) > 0));
580  for (int i = 0; i < depthMax_; ++i)
581  eHcalDetId[i] = HcalDetId();
582  for (unsigned int i = 0; i < ehdepth.size(); ++i) {
583  HcalSubdetector subdet0 =
584  (hborhe) ? ((ehdepth[i].second >= depth16HE(ieta, iphi)) ? HcalEndcap : HcalBarrel) : subdet;
585  HcalDetId hcid0(subdet0, ieta, iphi, ehdepth[i].second);
586  double actL = activeLength(DetId(hcid0));
587  double ene = ehdepth[i].first;
588  bool tmpC(false);
589  if (ene > 0.0) {
590  if (!(theHBHETopology_->validHcal(hcid0))) {
591  edm::LogWarning("HBHEMuon") << "(1) Invalid ID " << hcid0 << " with E = " << ene;
592  edm::LogWarning("HBHEMuon") << HcalDetId(closestCell) << " with " << ehdepth.size() << " depths:";
593  for (const auto& ehd : ehdepth)
594  edm::LogWarning("HBHEMuon") << " " << ehd.second << ":" << ehd.first;
595  } else {
596  tmpC = goodCell(hcid0, pTrack, geo_, bField);
597  double enec(ene);
598  if (unCorrect_) {
599  double corr = (ignoreHECorr_ && (subdet0 == HcalEndcap)) ? 1.0 : respCorr(DetId(hcid0));
600  if (corr != 0)
601  ene /= corr;
602 #ifdef EDM_ML_DEBUG
603  HcalDetId id = (isItPlan1_ && isItPreRecHit_) ? hdc_->mergedDepthDetId(hcid0) : hcid0;
604  edm::LogVerbatim("HBHEMuon") << hcid0 << ":" << id << " Corr " << corr;
605 #endif
606  }
607  int depth = ehdepth[i].second - 1;
608  if (collapseDepth_) {
609  HcalDetId id = hdc_->mergedDepthDetId(hcid0);
610  depth = id.depth() - 1;
611  }
612  eHcalDepth[depth] += ene;
613  eHcalDepthC[depth] += enec;
614  activeL[depth] += actL;
615  activeLengthTot += actL;
616  matchDepth[depth] = (matchDepth[depth] && tmpC);
617 #ifdef EDM_ML_DEBUG
618  if ((verbosity_ % 10) > 0)
619  edm::LogVerbatim("HBHEMuon")
620  << hcid0 << " E " << ene << ":" << enec << " L " << actL << " Match " << tmpC;
621 #endif
622  }
623  }
624  }
625 #ifdef EDM_ML_DEBUG
626  if ((verbosity_ % 10) > 0) {
627  edm::LogVerbatim("HBHEMuon") << hcidt << " Match " << tmpmatch << " Depths " << ehdepth.size();
628  for (unsigned int k = 0; k < ehdepth.size(); ++k)
629  edm::LogVerbatim("HBHEMuon") << " [" << k << ":" << ehdepth[k].second << "] " << matchDepth[k];
630  }
631 #endif
632  HcalDetId hotCell;
633  spr::eHCALmatrix(geo_, theHBHETopology_, closestCell, hbhe, 1, 1, hotCell, false, useRaw_, false);
634  isHot = matchId(closestCell, hotCell);
635  if (hotCell != HcalDetId()) {
636  subdet = HcalDetId(hotCell).subdet();
637  ieta = HcalDetId(hotCell).ieta();
638  iphi = HcalDetId(hotCell).iphi();
639  hborhe = (std::abs(ieta) == 16);
640  std::vector<std::pair<double, int>> ehdepth;
641  spr::energyHCALCell(hotCell,
642  hbhe,
643  ehdepth,
644  maxDepth_,
645  -100.0,
646  -100.0,
647  -100.0,
648  -100.0,
649  -500.0,
650  500.0,
651  useRaw_,
652  depth16HE(ieta, iphi),
653  false); //(((verbosity_/1000)%10)>0 ));
654  for (int i = 0; i < depthMax_; ++i)
655  eHcalDetId[i] = HcalDetId();
656  for (unsigned int i = 0; i < ehdepth.size(); ++i) {
657  HcalSubdetector subdet0 =
658  (hborhe) ? ((ehdepth[i].second >= depth16HE(ieta, iphi)) ? HcalEndcap : HcalBarrel) : subdet;
659  HcalDetId hcid0(subdet0, ieta, iphi, ehdepth[i].second);
660  double actL = activeLength(DetId(hcid0));
661  double ene = ehdepth[i].first;
662  bool tmpC(false);
663  if (ene > 0.0) {
664  if (!(theHBHETopology_->validHcal(hcid0))) {
665  edm::LogWarning("HBHEMuon") << "(2) Invalid ID " << hcid0 << " with E = " << ene;
666  edm::LogWarning("HBHEMuon") << HcalDetId(hotCell) << " with " << ehdepth.size() << " depths:";
667  for (const auto& ehd : ehdepth)
668  edm::LogWarning("HBHEMuon") << " " << ehd.second << ":" << ehd.first;
669  } else {
670  tmpC = goodCell(hcid0, pTrack, geo_, bField);
671  double chg(ene), enec(ene);
672  if (unCorrect_) {
673  double corr = (ignoreHECorr_ && (subdet0 == HcalEndcap)) ? 1.0 : respCorr(DetId(hcid0));
674  if (corr != 0)
675  ene /= corr;
676 #ifdef EDM_ML_DEBUG
677  HcalDetId id = (isItPlan1_ && isItPreRecHit_) ? hdc_->mergedDepthDetId(hcid0) : hcid0;
678  edm::LogVerbatim("HBHEMuon")
679  << hcid0 << ":" << id << " Corr " << corr << " E " << ene << ":" << enec;
680 #endif
681  }
682  if (getCharge_) {
683  double gain = gainFactor(conditions, hcid0);
684  if (gain != 0)
685  chg /= gain;
686 #ifdef EDM_ML_DEBUG
687  edm::LogVerbatim("HBHEMuon") << hcid0 << " Gain " << gain << " C " << chg;
688 #endif
689  }
690  int depth = ehdepth[i].second - 1;
691  if (collapseDepth_) {
692  HcalDetId id = hdc_->mergedDepthDetId(hcid0);
693  depth = id.depth() - 1;
694  }
695  eHcalDepthHot[depth] += ene;
696  eHcalDepthHotC[depth] += enec;
697  cHcalDepthHot[depth] += chg;
698  activeHotL[depth] += actL;
699  activeLengthHotTot += actL;
700  matchDepthHot[depth] = (matchDepthHot[depth] && tmpC);
701 #ifdef EDM_ML_DEBUG
702  if ((verbosity_ % 10) > 0)
703  edm::LogVerbatim("HBHEMuon") << hcid0 << " depth " << depth << " E " << ene << ":" << enec << " C "
704  << chg << " L " << actL << " Match " << tmpC;
705 #endif
706  }
707  }
708  }
709 
710  HcalDetId oppCell(subdet, -ieta, iphi, HcalDetId(hotCell).depth());
711  std::vector<std::pair<double, int>> ehdeptho;
712  spr::energyHCALCell(oppCell,
713  hbhe,
714  ehdeptho,
715  maxDepth_,
716  -100.0,
717  -100.0,
718  -100.0,
719  -100.0,
720  -500.0,
721  500.0,
722  useRaw_,
723  depth16HE(-ieta, iphi),
724  false); //(((verbosity_/1000)%10)>0));
725  for (unsigned int i = 0; i < ehdeptho.size(); ++i) {
726  HcalSubdetector subdet0 =
727  (hborhe) ? ((ehdeptho[i].second >= depth16HE(-ieta, iphi)) ? HcalEndcap : HcalBarrel) : subdet;
728  HcalDetId hcid0(subdet0, -ieta, iphi, ehdeptho[i].second);
729  double ene = ehdeptho[i].first;
730  if (ene > 0.0) {
731  if (!(theHBHETopology_->validHcal(hcid0))) {
732  edm::LogWarning("HBHEMuon") << "(3) Invalid ID " << hcid0 << " with E = " << ene;
733  edm::LogWarning("HBHEMuon") << oppCell << " with " << ehdeptho.size() << " depths:";
734  for (const auto& ehd : ehdeptho)
735  edm::LogWarning("HBHEMuon") << " " << ehd.second << ":" << ehd.first;
736  } else {
737  double chg(ene);
738  if (unCorrect_) {
739  double corr = (ignoreHECorr_ && (subdet0 == HcalEndcap)) ? 1.0 : respCorr(DetId(hcid0));
740  if (corr != 0)
741  ene /= corr;
742 #ifdef EDM_ML_DEBUG
743  HcalDetId id = (isItPlan1_ && isItPreRecHit_) ? hdc_->mergedDepthDetId(hcid0) : hcid0;
744  edm::LogVerbatim("HBHEMuon")
745  << hcid0 << ":" << id << " Corr " << corr << " E " << ene << ":" << ehdeptho[i].first;
746 #endif
747  }
748  if (getCharge_) {
749  double gain = gainFactor(conditions, hcid0);
750  if (gain != 0)
751  chg /= gain;
752 #ifdef EDM_ML_DEBUG
753  edm::LogVerbatim("HBHEMuon") << hcid0 << " Gain " << gain << " C " << chg;
754 #endif
755  }
756  int depth = ehdeptho[i].second - 1;
757  if (collapseDepth_) {
758  HcalDetId id = hdc_->mergedDepthDetId(hcid0);
759  depth = id.depth() - 1;
760  }
761  cHcalDepthHotBG[depth] += chg;
762 #ifdef EDM_ML_DEBUG
763  if ((verbosity_ % 10) > 0)
764  edm::LogVerbatim("HBHEMuon") << hcid0 << " Depth " << depth << " E " << ene << " C " << chg;
765 #endif
766  }
767  }
768  }
769  }
770  }
771 #ifdef EDM_ML_DEBUG
772  edm::LogVerbatim("HBHEMuon") << "Propagate Track to HCAL: " << trackID.okHCAL << " Match " << tmpmatch
773  << " Hot " << isHot << " Energy " << eHcal << std::endl;
774 #endif
775 
776  } else {
777  ecalDetId.push_back(0);
778  hcalDetId.push_back(0);
779  ehcalDetId.push_back(0);
780  }
781 
782  matchedId.push_back(tmpmatch);
783  ecal3x3Energy.push_back(eEcal);
784  hcal1x1Energy.push_back(eHcal);
785  hcal_ieta.push_back(ieta);
786  hcal_iphi.push_back(iphi);
787  for (int i = 0; i < depthMax_; ++i) {
788  hcalDepthEnergy[i].push_back(eHcalDepth[i]);
789  hcalDepthActiveLength[i].push_back(activeL[i]);
790  hcalDepthEnergyHot[i].push_back(eHcalDepthHot[i]);
791  hcalDepthActiveLengthHot[i].push_back(activeHotL[i]);
792  hcalDepthEnergyCorr[i].push_back(eHcalDepthC[i]);
793  hcalDepthEnergyHotCorr[i].push_back(eHcalDepthHotC[i]);
794  hcalDepthChargeHot[i].push_back(cHcalDepthHot[i]);
795  hcalDepthChargeHotBG[i].push_back(cHcalDepthHotBG[i]);
796  hcalDepthMatch[i].push_back(matchDepth[i]);
797  hcalDepthMatchHot[i].push_back(matchDepthHot[i]);
798  }
799  hcalActiveLength.push_back(activeLengthTot);
800  hcalHot.push_back(isHot);
801  hcalActiveLengthHot.push_back(activeLengthHotTot);
802  }
803  }
804  if (accept) {
805 #ifdef EDM_ML_DEBUG
806  for (unsigned int i = 0; i < hcal_ieta.size(); ++i)
807  edm::LogVerbatim("HBHEMuon") << "[" << i << "] ieta/iphi for entry to "
808  << "HCAL has value of " << hcal_ieta[i] << ":" << hcal_iphi[i];
809 #endif
810  for (unsigned int k = 0; k < muon_is_good.size(); ++k) {
811  muon_is_good_ = muon_is_good[k];
812  muon_global_ = muon_global[k];
813  muon_tracker_ = muon_tracker[k];
814  muon_is_tight_ = muon_is_tight[k];
815  muon_is_medium_ = muon_is_medium[k];
816  ptGlob_ = ptGlob[k];
817  etaGlob_ = etaGlob[k];
818  phiGlob_ = phiGlob[k];
819  energyMuon_ = energyMuon[k];
820  pMuon_ = pMuon[k];
821  muon_trkKink_ = muon_trkKink[k];
822  muon_chi2LocalPosition_ = muon_chi2LocalPosition[k];
823  muon_segComp_ = muon_segComp[k];
824  trackerLayer_ = trackerLayer[k];
825  numPixelLayers_ = numPixelLayers[k];
826  tight_PixelHits_ = tight_PixelHits[k];
827  innerTrack_ = innerTrack[k];
828  outerTrack_ = outerTrack[k];
829  globalTrack_ = globalTrack[k];
830  chiTracker_ = chiTracker[k];
831  dxyTracker_ = dxyTracker[k];
832  dzTracker_ = dzTracker[k];
833  innerTrackpt_ = innerTrackpt[k];
834  innerTracketa_ = innerTracketa[k];
835  innerTrackphi_ = innerTrackphi[k];
836  tight_validFraction_ = tight_validFraction[k];
837  outerTrackChi_ = outerTrackChi[k];
838  outerTrackPt_ = outerTrackPt[k];
839  outerTrackEta_ = outerTrackEta[k];
840  outerTrackPhi_ = outerTrackPhi[k];
841  outerTrackHits_ = outerTrackHits[k];
842  outerTrackRHits_ = outerTrackRHits[k];
843  globalTrckPt_ = globalTrckPt[k];
844  globalTrckEta_ = globalTrckEta[k];
845  globalTrckPhi_ = globalTrckPhi[k];
846  globalMuonHits_ = globalMuonHits[k];
847  matchedStat_ = matchedStat[k];
848  chiGlobal_ = chiGlobal[k];
849  tight_LongPara_ = tight_LongPara[k];
850  tight_TransImpara_ = tight_TransImpara[k];
851  isolationR04_ = isolationR04[k];
852  isolationR03_ = isolationR03[k];
853  ecalEnergy_ = ecalEnergy[k];
854  hcalEnergy_ = hcalEnergy[k];
855  hoEnergy_ = hoEnergy[k];
856  matchedId_ = matchedId[k];
857  hcalHot_ = hcalHot[k];
858  ecal3x3Energy_ = ecal3x3Energy[k];
859  hcal1x1Energy_ = hcal1x1Energy[k];
860  ecalDetId_ = ecalDetId[k];
861  hcalDetId_ = hcalDetId[k];
862  ehcalDetId_ = ehcalDetId[k];
863  hcal_ieta_ = hcal_ieta[k];
864  hcal_iphi_ = hcal_iphi[k];
865  for (int i = 0; i < depthMax_; ++i) {
866  hcalDepthEnergy_[i] = hcalDepthEnergy[i][k];
867  hcalDepthActiveLength_[i] = hcalDepthActiveLength[i][k];
868  hcalDepthEnergyHot_[i] = hcalDepthEnergyHot[i][k];
869  hcalDepthActiveLengthHot_[i] = hcalDepthActiveLengthHot[i][k];
870  hcalDepthChargeHot_[i] = hcalDepthChargeHot[i][k];
871  hcalDepthChargeHotBG_[i] = hcalDepthChargeHotBG[i][k];
872  hcalDepthEnergyCorr_[i] = hcalDepthEnergyCorr[i][k];
873  hcalDepthEnergyHotCorr_[i] = hcalDepthEnergyHotCorr[i][k];
874  hcalDepthMatch_[i] = hcalDepthMatch[i][k];
875  hcalDepthMatchHot_[i] = hcalDepthMatchHot[i][k];
876  }
877  hcalActiveLength_ = hcalActiveLength[k];
878  hcalActiveLengthHot_ = hcalActiveLengthHot[k];
879  tree_->Fill();
880  }
881  }
882 }
HcalHBHEMuonAnalyzer::gainFactor
double gainFactor(const HcalDbService *dbserv, const HcalDetId &id)
Definition: HcalHBHEMuonAnalyzer.cc:1235
HcalHBHEMuonAnalyzer::hltresults_
std::vector< int > hltresults_
Definition: HcalHBHEMuonAnalyzer.cc:168
edm::HLTGlobalStatus::accept
bool accept() const
Has at least one path accepted the event?
Definition: HLTGlobalStatus.h:51
edm::LogVerbatim
Log< level::Info, true > LogVerbatim
Definition: MessageLogger.h:128
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
HcalHBHEMuonAnalyzer::theHBHETopology_
const HcalTopology * theHBHETopology_
Definition: HcalHBHEMuonAnalyzer.cc:124
HcalHBHEMuonAnalyzer::hcalDepthEnergyHot_
double hcalDepthEnergyHot_[depthMax_]
Definition: HcalHBHEMuonAnalyzer.cc:158
info
static const TGPicture * info(bool iBackgroundIsBlack)
Definition: FWCollectionSummaryWidget.cc:153
HcalHBHEMuonAnalyzer::all_triggers_
std::vector< std::string > all_triggers_
Definition: HcalHBHEMuonAnalyzer.cc:167
mps_fire.i
i
Definition: mps_fire.py:429
RPCNoise_example.check
check
Definition: RPCNoise_example.py:71
extraflags_cff.vtx
vtx
Definition: extraflags_cff.py:19
HcalDDDRecConstants::mergedDepthDetId
HcalDetId mergedDepthDetId(const HcalDetId &id) const
Definition: HcalDDDRecConstants.cc:653
muon::isMediumMuon
bool isMediumMuon(const reco::Muon &, bool run2016_hip_mitigation=false)
Definition: MuonSelectors.cc:899
chg
const float chg[109]
Definition: CoreSimTrack.cc:5
HcalHBHEMuonAnalyzer::tok_Muon_
const edm::EDGetTokenT< reco::MuonCollection > tok_Muon_
Definition: HcalHBHEMuonAnalyzer.cc:111
edm::Handle::product
T const * product() const
Definition: Handle.h:70
HcalHBHEMuonAnalyzer::tok_sevlv_
const edm::ESGetToken< EcalSeverityLevelAlgo, EcalSeverityLevelAlgoRcd > tok_sevlv_
Definition: HcalHBHEMuonAnalyzer.cc:119
HcalHBHEMuonAnalyzer::respCorr
double respCorr(const DetId &id)
Definition: HcalHBHEMuonAnalyzer.cc:1223
HcalHBHEMuonAnalyzer::hcalDepthChargeHot_
double hcalDepthChargeHot_[depthMax_]
Definition: HcalHBHEMuonAnalyzer.cc:160
HcalHBHEMuonAnalyzer::tok_topo_
const edm::ESGetToken< CaloTopology, CaloTopologyRecord > tok_topo_
Definition: HcalHBHEMuonAnalyzer.cc:120
accept
bool accept(const edm::Event &event, const edm::TriggerResults &triggerTable, const std::string &triggerPath)
Definition: TopDQMHelpers.h:31
cutBasedMuonId_MuonPOG_V0_cff.globalTrack
globalTrack
Definition: cutBasedMuonId_MuonPOG_V0_cff.py:135
HcalHBHEMuonAnalyzer::hcalDepthActiveLengthHot_
double hcalDepthActiveLengthHot_[depthMax_]
Definition: HcalHBHEMuonAnalyzer.cc:159
HcalHBHEMuonAnalyzer::hcalDepthEnergyCorr_
double hcalDepthEnergyCorr_[depthMax_]
Definition: HcalHBHEMuonAnalyzer.cc:162
Calorimetry_cff.bField
bField
Definition: Calorimetry_cff.py:284
edm::TriggerNames::triggerIndex
unsigned int triggerIndex(std::string_view name) const
Definition: TriggerNames.cc:52
HcalHBHEMuonAnalyzer::hcalDepthActiveLength_
double hcalDepthActiveLength_[depthMax_]
Definition: HcalHBHEMuonAnalyzer.cc:157
edm::second
U second(std::pair< T, U > const &p)
Definition: ParameterSet.cc:222
photonIsolationHIProducer_cfi.hbhe
hbhe
Definition: photonIsolationHIProducer_cfi.py:8
HcalHBHEMuonAnalyzer::tok_dbservice_
const edm::ESGetToken< HcalDbService, HcalDbRecord > tok_dbservice_
Definition: HcalHBHEMuonAnalyzer.cc:121
HcalDetId::subdet
constexpr HcalSubdetector subdet() const
get the subdetector
Definition: HcalDetId.h:138
HcalHBHEMuonAnalyzer::hcalDepthEnergyHotCorr_
double hcalDepthEnergyHotCorr_[depthMax_]
Definition: HcalHBHEMuonAnalyzer.cc:163
HcalHBHEMuonAnalyzer::hcalDepthMatchHot_
bool hcalDepthMatchHot_[depthMax_]
Definition: HcalHBHEMuonAnalyzer.cc:165
edm::HLTGlobalStatus::size
unsigned int size() const
Get number of paths stored.
Definition: HLTGlobalStatus.h:37
iEvent
int iEvent
Definition: GenABIO.cc:224
muons2muons_cfi.isolationR03
isolationR03
Definition: muons2muons_cfi.py:25
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
HcalDetId::ieta
constexpr int ieta() const
get the cell ieta
Definition: HcalDetId.h:155
HcalSubdetector
HcalSubdetector
Definition: HcalAssistant.h:31
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
HcalHBHEMuonAnalyzer::tok_HBHE_
const edm::EDGetTokenT< HBHERecHitCollection > tok_HBHE_
Definition: HcalHBHEMuonAnalyzer.cc:110
HcalHBHEMuonAnalyzer::geo_
const CaloGeometry * geo_
Definition: HcalHBHEMuonAnalyzer.cc:125
muon::segmentCompatibility
float segmentCompatibility(const reco::Muon &muon, reco::Muon::ArbitrationType arbitrationType=reco::Muon::SegmentAndTrackArbitration)
Definition: MuonSelectors.cc:61
PedestalClient_cfi.gain
gain
Definition: PedestalClient_cfi.py:37
HcalHBHEMuonAnalyzer::activeLength
double activeLength(const DetId &)
Definition: HcalHBHEMuonAnalyzer.cc:1176
HcalHBHEMuonAnalyzer::goodCell
bool goodCell(const HcalDetId &hcid, const reco::Track *pTrack, const CaloGeometry *geo, const MagneticField *bField)
Definition: HcalHBHEMuonAnalyzer.cc:1258
HcalHBHEMuonAnalyzer::hcalDepthChargeHotBG_
double hcalDepthChargeHotBG_[depthMax_]
Definition: HcalHBHEMuonAnalyzer.cc:161
DetId
Definition: DetId.h:17
HcalHBHEMuonAnalyzer::tok_trigRes_
const edm::EDGetTokenT< edm::TriggerResults > tok_trigRes_
Definition: HcalHBHEMuonAnalyzer.cc:106
HcalHBHEMuonAnalyzer::hcalDepthEnergy_
double hcalDepthEnergy_[depthMax_]
Definition: HcalHBHEMuonAnalyzer.cc:156
HcalHBHEMuonAnalyzer::depth16HE
int depth16HE(int ieta, int iphi)
Definition: HcalHBHEMuonAnalyzer.cc:1243
math::XYZPoint
XYZPointD XYZPoint
point in space with cartesian internal representation
Definition: Point3D.h:12
HcalHBHEMuonAnalyzer::hcalDepthMatch_
bool hcalDepthMatch_[depthMax_]
Definition: HcalHBHEMuonAnalyzer.cc:164
HcalHBHEMuonAnalyzer::tok_EE_
const edm::EDGetTokenT< EcalRecHitCollection > tok_EE_
Definition: HcalHBHEMuonAnalyzer.cc:109
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:70
HcalTopology::validHcal
bool validHcal(const HcalDetId &id) const
Definition: HcalTopology.cc:230
HcalHBHEMuonAnalyzer::matchId
int matchId(const HcalDetId &, const HcalDetId &)
Definition: HcalHBHEMuonAnalyzer.cc:1169
HcalHBHEMuonAnalyzer::tok_EB_
const edm::EDGetTokenT< EcalRecHitCollection > tok_EB_
Definition: HcalHBHEMuonAnalyzer.cc:108
HcalHBHEMuonAnalyzer::hdc_
const HcalDDDRecConstants * hdc_
Definition: HcalHBHEMuonAnalyzer.cc:123
HcalHBHEMuonAnalyzer::tok_chan_
const edm::ESGetToken< EcalChannelStatus, EcalChannelStatusRcd > tok_chan_
Definition: HcalHBHEMuonAnalyzer.cc:118
HcalHBHEMuonAnalyzer::tok_Vtx_
const edm::EDGetTokenT< reco::VertexCollection > tok_Vtx_
Definition: HcalHBHEMuonAnalyzer.cc:107
HcalHBHEMuonAnalyzer::tok_magField_
const edm::ESGetToken< MagneticField, IdealMagneticFieldRecord > tok_magField_
Definition: HcalHBHEMuonAnalyzer.cc:117
edm::LogWarning
Log< level::Warning, false > LogWarning
Definition: MessageLogger.h:122
muon::isTightMuon
bool isTightMuon(const reco::Muon &, const reco::Vertex &)
Definition: MuonSelectors.cc:880
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)
muons2muons_cfi.isolationR04
isolationR04
Definition: muons2muons_cfi.py:32
HcalDetId::iphi
constexpr int iphi() const
get the cell iphi
Definition: HcalDetId.h:157
HcalHBHEMuonAnalyzer::isGoodVertex
bool isGoodVertex(const reco::Vertex &vtx)
Definition: HcalHBHEMuonAnalyzer.cc:1211
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 HcalHBHEMuonAnalyzer::beginJob ( )
overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 885 of file HcalHBHEMuonAnalyzer.cc.

References all_triggers_, bxNumber_, chiGlobal_, chiTracker_, dxyTracker_, dzTracker_, ecal3x3Energy_, ecalDetId_, ecalEnergy_, ehcalDetId_, energyMuon_, etaGlob_, eventNumber_, compareTotals::fs, globalMuonHits_, globalTrack_, globalTrckEta_, globalTrckPhi_, globalTrckPt_, goodVertex_, hcal1x1Energy_, hcal_ieta_, hcal_iphi_, hcalActiveLength_, hcalActiveLengthHot_, hcalDepthActiveLength_, hcalDepthActiveLengthHot_, hcalDepthChargeHot_, hcalDepthChargeHotBG_, hcalDepthEnergy_, hcalDepthEnergyCorr_, hcalDepthEnergyHot_, hcalDepthEnergyHotCorr_, hcalDepthMatch_, hcalDepthMatchHot_, hcalDetId_, hcalEnergy_, hcalHot_, hltresults_, hoEnergy_, innerTrack_, innerTracketa_, innerTrackphi_, innerTrackpt_, isolationR03_, isolationR04_, dqmdumpme::k, lumiNumber_, matchedId_, matchedStat_, maxDepth_, muon_chi2LocalPosition_, muon_global_, muon_is_good_, muon_is_medium_, muon_is_tight_, muon_segComp_, muon_tracker_, muon_trkKink_, Skims_PA_cff::name, numPixelLayers_, outerTrack_, outerTrackChi_, outerTrackEta_, outerTrackHits_, outerTrackPhi_, outerTrackPt_, outerTrackRHits_, phiGlob_, pMuon_, ptGlob_, runNumber_, tight_LongPara_, tight_PixelHits_, tight_TransImpara_, tight_validFraction_, trackerLayer_, and tree_.

885  {
886  edm::Service<TFileService> fs;
887  tree_ = fs->make<TTree>("TREE", "TREE");
888  tree_->Branch("Event_No", &eventNumber_);
889  tree_->Branch("Run_No", &runNumber_);
890  tree_->Branch("LumiNumber", &lumiNumber_);
891  tree_->Branch("BXNumber", &bxNumber_);
892  tree_->Branch("GoodVertex", &goodVertex_);
893  tree_->Branch("PF_Muon", &muon_is_good_);
894  tree_->Branch("Global_Muon", &muon_global_);
895  tree_->Branch("Tracker_muon", &muon_tracker_);
896  tree_->Branch("MuonIsTight", &muon_is_tight_);
897  tree_->Branch("MuonIsMedium", &muon_is_medium_);
898  tree_->Branch("pt_of_muon", &ptGlob_);
899  tree_->Branch("eta_of_muon", &etaGlob_);
900  tree_->Branch("phi_of_muon", &phiGlob_);
901  tree_->Branch("energy_of_muon", &energyMuon_);
902  tree_->Branch("p_of_muon", &pMuon_);
903  tree_->Branch("muon_trkKink", &muon_trkKink_);
904  tree_->Branch("muon_chi2LocalPosition", &muon_chi2LocalPosition_);
905  tree_->Branch("muon_segComp", &muon_segComp_);
906 
907  tree_->Branch("TrackerLayer", &trackerLayer_);
908  tree_->Branch("NumPixelLayers", &numPixelLayers_);
909  tree_->Branch("InnerTrackPixelHits", &tight_PixelHits_);
910  tree_->Branch("innerTrack", &innerTrack_);
911  tree_->Branch("chiTracker", &chiTracker_);
912  tree_->Branch("DxyTracker", &dxyTracker_);
913  tree_->Branch("DzTracker", &dzTracker_);
914  tree_->Branch("innerTrackpt", &innerTrackpt_);
915  tree_->Branch("innerTracketa", &innerTracketa_);
916  tree_->Branch("innerTrackphi", &innerTrackphi_);
917  tree_->Branch("tight_validFraction", &tight_validFraction_);
918 
919  tree_->Branch("OuterTrack", &outerTrack_);
920  tree_->Branch("OuterTrackChi", &outerTrackChi_);
921  tree_->Branch("OuterTrackPt", &outerTrackPt_);
922  tree_->Branch("OuterTrackEta", &outerTrackEta_);
923  tree_->Branch("OuterTrackPhi", &outerTrackPhi_);
924  tree_->Branch("OuterTrackHits", &outerTrackHits_);
925  tree_->Branch("OuterTrackRHits", &outerTrackRHits_);
926 
927  tree_->Branch("GlobalTrack", &globalTrack_);
928  tree_->Branch("GlobalTrckPt", &globalTrckPt_);
929  tree_->Branch("GlobalTrckEta", &globalTrckEta_);
930  tree_->Branch("GlobalTrckPhi", &globalTrckPhi_);
931  tree_->Branch("Global_Muon_Hits", &globalMuonHits_);
932  tree_->Branch("MatchedStations", &matchedStat_);
933  tree_->Branch("GlobTrack_Chi", &chiGlobal_);
934  tree_->Branch("Tight_LongitudinalImpactparameter", &tight_LongPara_);
935  tree_->Branch("Tight_TransImpactparameter", &tight_TransImpara_);
936 
937  tree_->Branch("IsolationR04", &isolationR04_);
938  tree_->Branch("IsolationR03", &isolationR03_);
939  tree_->Branch("ecal_3into3", &ecalEnergy_);
940  tree_->Branch("hcal_3into3", &hcalEnergy_);
941  tree_->Branch("tracker_3into3", &hoEnergy_);
942 
943  tree_->Branch("matchedId", &matchedId_);
944  tree_->Branch("hcal_cellHot", &hcalHot_);
945 
946  tree_->Branch("ecal_3x3", &ecal3x3Energy_);
947  tree_->Branch("hcal_1x1", &hcal1x1Energy_);
948  tree_->Branch("ecal_detID", &ecalDetId_);
949  tree_->Branch("hcal_detID", &hcalDetId_);
950  tree_->Branch("ehcal_detID", &ehcalDetId_);
951  tree_->Branch("hcal_ieta", &hcal_ieta_);
952  tree_->Branch("hcal_iphi", &hcal_iphi_);
953 
954  char name[100];
955  for (int k = 0; k < maxDepth_; ++k) {
956  sprintf(name, "hcal_edepth%d", (k + 1));
957  tree_->Branch(name, &hcalDepthEnergy_[k]);
958  sprintf(name, "hcal_activeL%d", (k + 1));
959  tree_->Branch(name, &hcalDepthActiveLength_[k]);
960  sprintf(name, "hcal_edepthHot%d", (k + 1));
961  tree_->Branch(name, &hcalDepthEnergyHot_[k]);
962  sprintf(name, "hcal_activeHotL%d", (k + 1));
963  tree_->Branch(name, &hcalDepthActiveLengthHot_[k]);
964  sprintf(name, "hcal_cdepthHot%d", (k + 1));
965  tree_->Branch(name, &hcalDepthChargeHot_[k]);
966  sprintf(name, "hcal_cdepthHotBG%d", (k + 1));
967  tree_->Branch(name, &hcalDepthChargeHotBG_[k]);
968  sprintf(name, "hcal_edepthCorrect%d", (k + 1));
969  tree_->Branch(name, &hcalDepthEnergyCorr_[k]);
970  sprintf(name, "hcal_edepthHotCorrect%d", (k + 1));
971  tree_->Branch(name, &hcalDepthEnergyHotCorr_[k]);
972  sprintf(name, "hcal_depthMatch%d", (k + 1));
973  tree_->Branch(name, &hcalDepthMatch_[k]);
974  sprintf(name, "hcal_depthMatchHot%d", (k + 1));
975  tree_->Branch(name, &hcalDepthMatchHot_[k]);
976  }
977 
978  tree_->Branch("activeLength", &hcalActiveLength_);
979  tree_->Branch("activeLengthHot", &hcalActiveLengthHot_);
980 
981  tree_->Branch("hltresults", &hltresults_);
982  tree_->Branch("all_triggers", &all_triggers_);
983 }
HcalHBHEMuonAnalyzer::hltresults_
std::vector< int > hltresults_
Definition: HcalHBHEMuonAnalyzer.cc:168
HcalHBHEMuonAnalyzer::hcalDepthEnergyHot_
double hcalDepthEnergyHot_[depthMax_]
Definition: HcalHBHEMuonAnalyzer.cc:158
HcalHBHEMuonAnalyzer::all_triggers_
std::vector< std::string > all_triggers_
Definition: HcalHBHEMuonAnalyzer.cc:167
HcalHBHEMuonAnalyzer::hcalDepthChargeHot_
double hcalDepthChargeHot_[depthMax_]
Definition: HcalHBHEMuonAnalyzer.cc:160
HcalHBHEMuonAnalyzer::hcalDepthActiveLengthHot_
double hcalDepthActiveLengthHot_[depthMax_]
Definition: HcalHBHEMuonAnalyzer.cc:159
HcalHBHEMuonAnalyzer::hcalDepthEnergyCorr_
double hcalDepthEnergyCorr_[depthMax_]
Definition: HcalHBHEMuonAnalyzer.cc:162
HcalHBHEMuonAnalyzer::hcalDepthActiveLength_
double hcalDepthActiveLength_[depthMax_]
Definition: HcalHBHEMuonAnalyzer.cc:157
HcalHBHEMuonAnalyzer::hcalDepthEnergyHotCorr_
double hcalDepthEnergyHotCorr_[depthMax_]
Definition: HcalHBHEMuonAnalyzer.cc:163
HcalHBHEMuonAnalyzer::hcalDepthMatchHot_
bool hcalDepthMatchHot_[depthMax_]
Definition: HcalHBHEMuonAnalyzer.cc:165
HcalHBHEMuonAnalyzer::hcalDepthChargeHotBG_
double hcalDepthChargeHotBG_[depthMax_]
Definition: HcalHBHEMuonAnalyzer.cc:161
HcalHBHEMuonAnalyzer::hcalDepthEnergy_
double hcalDepthEnergy_[depthMax_]
Definition: HcalHBHEMuonAnalyzer.cc:156
HcalHBHEMuonAnalyzer::hcalDepthMatch_
bool hcalDepthMatch_[depthMax_]
Definition: HcalHBHEMuonAnalyzer.cc:164
Skims_PA_cff.name
name
Definition: Skims_PA_cff.py:17

◆ beginRun()

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

Definition at line 986 of file HcalHBHEMuonAnalyzer.cc.

References actHB, actHE, activeLength(), all_triggers_, geo_, edm::EventSetup::getData(), CaloGeometry::getSubdetectorGeometry(), HcalDDDRecConstants::getThickActive(), HcalGeometry::getValidDetIds(), reco::JetExtendedAssociation::getValue(), HcalCondObjectContainer< Item >::getValues(), DetId::Hcal, HcalBarrel, HcalEndcap, hdc_, hltConfig_, HLTConfigProvider::init(), mergedDepth_, HcalDDDRecConstants::mergedDepthDetId(), DetId::rawId(), respCorrs_, edm::RunBase::run(), HcalCondObjectContainerBase::setTopo(), HLTConfigProvider::size(), AlCaHLTBitMon_QueryRunRegistry::string, submitPVValidationJobs::t, theHBHETopology_, tok_ddrec_, tok_geom_, tok_htopo_, tok_respcorr_, HLTConfigProvider::triggerName(), triggers_, and writeRespCorr_.

986  {
987  hdc_ = &iSetup.getData(tok_ddrec_);
988  actHB.clear();
989  actHE.clear();
990  actHB = hdc_->getThickActive(0);
991  actHE = hdc_->getThickActive(1);
992 #ifdef EDM_ML_DEBUG
993  unsigned int k1(0), k2(0);
994  edm::LogVerbatim("HBHEMuon") << actHB.size() << " Active Length for HB";
995  for (const auto& act : actHB) {
996  edm::LogVerbatim("HBHEMuon") << "[" << k1 << "] ieta " << act.ieta << " depth " << act.depth << " zside "
997  << act.zside << " type " << act.stype << " phi " << act.iphis.size() << ":"
998  << act.iphis[0] << " L " << act.thick;
999  HcalDetId hcid1(HcalBarrel, (act.ieta) * (act.zside), act.iphis[0], act.depth);
1000  HcalDetId hcid2 = mergedDepth_ ? hdc_->mergedDepthDetId(hcid1) : hcid1;
1001  edm::LogVerbatim("HBHEMuon") << hcid1 << " | " << hcid2 << " L " << activeLength(DetId(hcid2));
1002  ++k1;
1003  }
1004  edm::LogVerbatim("HBHEMuon") << actHE.size() << " Active Length for HE";
1005  for (const auto& act : actHE) {
1006  edm::LogVerbatim("HBHEMuon") << "[" << k2 << "] ieta " << act.ieta << " depth " << act.depth << " zside "
1007  << act.zside << " type " << act.stype << " phi " << act.iphis.size() << ":"
1008  << act.iphis[0] << " L " << act.thick;
1009  HcalDetId hcid1(HcalEndcap, (act.ieta) * (act.zside), act.iphis[0], act.depth);
1010  HcalDetId hcid2 = mergedDepth_ ? hdc_->mergedDepthDetId(hcid1) : hcid1;
1011  edm::LogVerbatim("HBHEMuon") << hcid1 << " | " << hcid2 << " L " << activeLength(DetId(hcid2));
1012  ++k2;
1013  }
1014 #endif
1015 
1016  bool changed = true;
1017  all_triggers_.clear();
1018  if (hltConfig_.init(iRun, iSetup, "HLT", changed)) {
1019  // if init returns TRUE, initialisation has succeeded!
1020 #ifdef EDM_ML_DEBUG
1021  edm::LogVerbatim("HBHEMuon") << "HLT config with process name HLT successfully extracted";
1022 #endif
1023  unsigned int ntriggers = hltConfig_.size();
1024  for (unsigned int t = 0; t < ntriggers; ++t) {
1025  std::string hltname(hltConfig_.triggerName(t));
1026  for (unsigned int ik = 0; ik < triggers_.size(); ++ik) {
1027  if (hltname.find(triggers_[ik]) != std::string::npos) {
1028  all_triggers_.push_back(hltname);
1029  break;
1030  }
1031  }
1032  } //loop over ntriggers
1033  edm::LogVerbatim("HBHEMuon") << "All triggers size in begin run " << all_triggers_.size();
1034  } else {
1035  edm::LogError("HBHEMuon") << "Error! HLT config extraction with process name HLT failed";
1036  }
1037 
1038  theHBHETopology_ = &iSetup.getData(tok_htopo_);
1039  const HcalRespCorrs* resp = &iSetup.getData(tok_respcorr_);
1040  respCorrs_ = new HcalRespCorrs(*resp);
1041  respCorrs_->setTopo(theHBHETopology_);
1042  geo_ = &iSetup.getData(tok_geom_);
1043 
1044  // Write correction factors for all HB/HE events
1045  if (writeRespCorr_) {
1046  const HcalGeometry* gHcal = static_cast<const HcalGeometry*>(geo_->getSubdetectorGeometry(DetId::Hcal, HcalBarrel));
1047  const std::vector<DetId>& ids = gHcal->getValidDetIds(DetId::Hcal, 0);
1048  edm::LogVerbatim("HBHEMuon") << "\nTable of Correction Factors for Run " << iRun.run() << "\n";
1049  for (auto const& id : ids) {
1050  if ((id.det() == DetId::Hcal) && ((id.subdetId() == HcalBarrel) || (id.subdetId() == HcalEndcap))) {
1051  edm::LogVerbatim("HBHEMuon") << HcalDetId(id) << " " << id.rawId() << " "
1052  << (respCorrs_->getValues(id))->getValue();
1053  }
1054  }
1055  }
1056 }
edm::LogVerbatim
Log< level::Info, true > LogVerbatim
Definition: MessageLogger.h:128
HLTConfigProvider::triggerName
const std::string & triggerName(unsigned int triggerIndex) const
Definition: HLTConfigProvider.h:71
HcalHBHEMuonAnalyzer::theHBHETopology_
const HcalTopology * theHBHETopology_
Definition: HcalHBHEMuonAnalyzer.cc:124
HcalHBHEMuonAnalyzer::all_triggers_
std::vector< std::string > all_triggers_
Definition: HcalHBHEMuonAnalyzer.cc:167
HcalHBHEMuonAnalyzer::tok_respcorr_
const edm::ESGetToken< HcalRespCorrs, HcalRespCorrsRcd > tok_respcorr_
Definition: HcalHBHEMuonAnalyzer.cc:115
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:653
HcalHBHEMuonAnalyzer::tok_htopo_
const edm::ESGetToken< HcalTopology, HcalRecNumberingRecord > tok_htopo_
Definition: HcalHBHEMuonAnalyzer.cc:114
HcalHBHEMuonAnalyzer::actHE
std::vector< HcalDDDRecConstants::HcalActiveLength > actHE
Definition: HcalHBHEMuonAnalyzer.cc:170
edm::LogError
Log< level::Error, false > LogError
Definition: MessageLogger.h:123
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
HcalCondObjectContainer::getValues
const Item * getValues(DetId fId, bool throwOnFail=true) const
Definition: HcalCondObjectContainer.h:159
HLTConfigProvider::size
unsigned int size() const
number of trigger paths in trigger table
Definition: HLTConfigProvider.h:64
HcalHBHEMuonAnalyzer::geo_
const CaloGeometry * geo_
Definition: HcalHBHEMuonAnalyzer.cc:125
HcalHBHEMuonAnalyzer::activeLength
double activeLength(const DetId &)
Definition: HcalHBHEMuonAnalyzer.cc:1176
DetId
Definition: DetId.h:17
HcalHBHEMuonAnalyzer::actHB
std::vector< HcalDDDRecConstants::HcalActiveLength > actHB
Definition: HcalHBHEMuonAnalyzer.cc:170
HcalHBHEMuonAnalyzer::tok_geom_
const edm::ESGetToken< CaloGeometry, CaloGeometryRecord > tok_geom_
Definition: HcalHBHEMuonAnalyzer.cc:116
DetId::rawId
constexpr uint32_t rawId() const
get the raw id
Definition: DetId.h:57
HLTConfigProvider::init
bool init(const edm::Run &iRun, const edm::EventSetup &iSetup, const std::string &processName, bool &changed)
d&#39;tor
Definition: HLTConfigProvider.cc:36
HcalHBHEMuonAnalyzer::hdc_
const HcalDDDRecConstants * hdc_
Definition: HcalHBHEMuonAnalyzer.cc:123
HcalHBHEMuonAnalyzer::hltConfig_
HLTConfigProvider hltConfig_
Definition: HcalHBHEMuonAnalyzer.cc:93
HcalDDDRecConstants::getThickActive
std::vector< HcalActiveLength > getThickActive(const int &type) const
Definition: HcalDDDRecConstants.cc:471
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
HcalHBHEMuonAnalyzer::tok_ddrec_
const edm::ESGetToken< HcalDDDRecConstants, HcalRecNumberingRecord > tok_ddrec_
Definition: HcalHBHEMuonAnalyzer.cc:113
HcalCondObjectContainerBase::setTopo
void setTopo(const HcalTopology *topo)
Definition: HcalCondObjectContainerBase.cc:17
HcalHBHEMuonAnalyzer::triggers_
const std::vector< std::string > triggers_
Definition: HcalHBHEMuonAnalyzer.cc:97

◆ clearVectors()

void HcalHBHEMuonAnalyzer::clearVectors ( )
private

clearing vectots

Definition at line 1086 of file HcalHBHEMuonAnalyzer.cc.

References bxNumber_, chiGlobal_, chiTracker_, dxyTracker_, dzTracker_, ecal3x3Energy_, ecalDetId_, ecalEnergy_, ehcalDetId_, energyMuon_, etaGlob_, eventNumber_, globalMuonHits_, globalTrack_, globalTrckEta_, globalTrckPhi_, globalTrckPt_, goodVertex_, hcal1x1Energy_, hcal_ieta_, hcal_iphi_, hcalActiveLength_, hcalActiveLengthHot_, hcalDepthActiveLength_, hcalDepthActiveLengthHot_, hcalDepthChargeHot_, hcalDepthChargeHotBG_, hcalDepthEnergy_, hcalDepthEnergyCorr_, hcalDepthEnergyHot_, hcalDepthEnergyHotCorr_, hcalDepthMatch_, hcalDepthMatchHot_, hcalDetId_, hcalEnergy_, hcalHot_, hltresults_, hoEnergy_, mps_fire::i, innerTrack_, innerTracketa_, innerTrackphi_, innerTrackpt_, isolationR03_, isolationR04_, lumiNumber_, matchedId_, matchedStat_, maxDepth_, muon_chi2LocalPosition_, muon_global_, muon_is_good_, muon_is_medium_, muon_is_tight_, muon_segComp_, muon_tracker_, muon_trkKink_, numPixelLayers_, outerTrack_, outerTrackChi_, outerTrackEta_, outerTrackHits_, outerTrackPhi_, outerTrackPt_, outerTrackRHits_, phiGlob_, pMuon_, ptGlob_, runNumber_, tight_LongPara_, tight_PixelHits_, tight_TransImpara_, tight_validFraction_, and trackerLayer_.

Referenced by analyze().

1086  {
1088  eventNumber_ = -99999;
1089  runNumber_ = -99999;
1090  lumiNumber_ = -99999;
1091  bxNumber_ = -99999;
1092  goodVertex_ = -99999;
1093 
1094  muon_is_good_ = false;
1095  muon_global_ = false;
1096  muon_tracker_ = false;
1097  ptGlob_ = 0;
1098  etaGlob_ = 0;
1099  phiGlob_ = 0;
1100  energyMuon_ = 0;
1101  pMuon_ = 0;
1102  muon_trkKink_ = 0;
1103  muon_chi2LocalPosition_ = 0;
1104  muon_segComp_ = 0;
1105  muon_is_tight_ = false;
1106  muon_is_medium_ = false;
1107 
1108  trackerLayer_ = 0;
1109  numPixelLayers_ = 0;
1110  tight_PixelHits_ = 0;
1111  innerTrack_ = false;
1112  chiTracker_ = 0;
1113  dxyTracker_ = 0;
1114  dzTracker_ = 0;
1115  innerTrackpt_ = 0;
1116  innerTracketa_ = 0;
1117  innerTrackphi_ = 0;
1118  tight_validFraction_ = 0;
1119 
1120  outerTrack_ = false;
1121  outerTrackPt_ = 0;
1122  outerTrackEta_ = 0;
1123  outerTrackPhi_ = 0;
1124  outerTrackHits_ = 0;
1125  outerTrackRHits_ = 0;
1126  outerTrackChi_ = 0;
1127 
1128  globalTrack_ = false;
1129  globalTrckPt_ = 0;
1130  globalTrckEta_ = 0;
1131  globalTrckPhi_ = 0;
1132  globalMuonHits_ = 0;
1133  matchedStat_ = 0;
1134  chiGlobal_ = 0;
1135  tight_LongPara_ = 0;
1136  tight_TransImpara_ = 0;
1137 
1138  isolationR04_ = 0;
1139  isolationR03_ = 0;
1140  ecalEnergy_ = 0;
1141  hcalEnergy_ = 0;
1142  hoEnergy_ = 0;
1143  matchedId_ = false;
1144  hcalHot_ = false;
1145  ecal3x3Energy_ = 0;
1146  hcal1x1Energy_ = 0;
1147  ecalDetId_ = 0;
1148  hcalDetId_ = 0;
1149  ehcalDetId_ = 0;
1150  hcal_ieta_ = 0;
1151  hcal_iphi_ = 0;
1152  for (int i = 0; i < maxDepth_; ++i) {
1153  hcalDepthEnergy_[i] = 0;
1154  hcalDepthActiveLength_[i] = 0;
1155  hcalDepthEnergyHot_[i] = 0;
1156  hcalDepthActiveLengthHot_[i] = 0;
1157  hcalDepthChargeHot_[i] = 0;
1158  hcalDepthChargeHotBG_[i] = 0;
1159  hcalDepthEnergyCorr_[i] = 0;
1160  hcalDepthEnergyHotCorr_[i] = 0;
1161  hcalDepthMatch_[i] = false;
1162  hcalDepthMatchHot_[i] = false;
1163  }
1164  hcalActiveLength_ = 0;
1165  hcalActiveLengthHot_ = 0;
1166  hltresults_.clear();
1167 }
HcalHBHEMuonAnalyzer::hltresults_
std::vector< int > hltresults_
Definition: HcalHBHEMuonAnalyzer.cc:168
HcalHBHEMuonAnalyzer::hcalDepthEnergyHot_
double hcalDepthEnergyHot_[depthMax_]
Definition: HcalHBHEMuonAnalyzer.cc:158
mps_fire.i
i
Definition: mps_fire.py:429
HcalHBHEMuonAnalyzer::hcalDepthChargeHot_
double hcalDepthChargeHot_[depthMax_]
Definition: HcalHBHEMuonAnalyzer.cc:160
HcalHBHEMuonAnalyzer::hcalDepthActiveLengthHot_
double hcalDepthActiveLengthHot_[depthMax_]
Definition: HcalHBHEMuonAnalyzer.cc:159
HcalHBHEMuonAnalyzer::hcalDepthEnergyCorr_
double hcalDepthEnergyCorr_[depthMax_]
Definition: HcalHBHEMuonAnalyzer.cc:162
HcalHBHEMuonAnalyzer::hcalDepthActiveLength_
double hcalDepthActiveLength_[depthMax_]
Definition: HcalHBHEMuonAnalyzer.cc:157
HcalHBHEMuonAnalyzer::hcalDepthEnergyHotCorr_
double hcalDepthEnergyHotCorr_[depthMax_]
Definition: HcalHBHEMuonAnalyzer.cc:163
HcalHBHEMuonAnalyzer::hcalDepthMatchHot_
bool hcalDepthMatchHot_[depthMax_]
Definition: HcalHBHEMuonAnalyzer.cc:165
HcalHBHEMuonAnalyzer::hcalDepthChargeHotBG_
double hcalDepthChargeHotBG_[depthMax_]
Definition: HcalHBHEMuonAnalyzer.cc:161
HcalHBHEMuonAnalyzer::hcalDepthEnergy_
double hcalDepthEnergy_[depthMax_]
Definition: HcalHBHEMuonAnalyzer.cc:156
HcalHBHEMuonAnalyzer::hcalDepthMatch_
bool hcalDepthMatch_[depthMax_]
Definition: HcalHBHEMuonAnalyzer.cc:164

◆ depth16HE()

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

Definition at line 1243 of file HcalHBHEMuonAnalyzer.cc.

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

Referenced by analyze().

1243  {
1244  // Transition between HB/HE is special
1245  // For Run 1 or for Plan1 standard reconstruction it is 3
1246  // For runs beyond 2018 or in Plan1 for HEP17 it is 4
1247  int zside = (ieta > 0) ? 1 : -1;
1248  int depth = theHBHETopology_->dddConstants()->getMinDepth(1, 16, iphi, zside);
1249  if (isItPlan1_ && (!isItPreRecHit_))
1250  depth = 3;
1251 #ifdef EDM_ML_DEBUG
1252  edm::LogVerbatim("HBHEMuon") << "Plan1 " << isItPlan1_ << " PreRecHit " << isItPreRecHit_ << " phi " << iphi
1253  << " depth " << depth;
1254 #endif
1255  return depth;
1256 }
HcalDDDRecConstants::getMinDepth
int getMinDepth(const int &itype, const int &ieta, const int &iphi, const int &zside) const
Definition: HcalDDDRecConstants.cc:359
edm::LogVerbatim
Log< level::Info, true > LogVerbatim
Definition: MessageLogger.h:128
HcalHBHEMuonAnalyzer::theHBHETopology_
const HcalTopology * theHBHETopology_
Definition: HcalHBHEMuonAnalyzer.cc:124
ecaldqm::zside
int zside(DetId const &)
Definition: EcalDQMCommonUtils.cc:189
HcalTopology::dddConstants
const HcalDDDRecConstants * dddConstants() const
Definition: HcalTopology.h:164

◆ endRun()

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

Definition at line 82 of file HcalHBHEMuonAnalyzer.cc.

82 {}

◆ fillDescriptions()

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

Definition at line 1059 of file HcalHBHEMuonAnalyzer.cc.

References edm::ConfigurationDescriptions::add(), submitPVResolutionJobs::desc, HLT_2022v15_cff::InputTag, AlCaHLTBitMon_QueryRunRegistry::string, and ClusterTask_cfi::trig.

1059  {
1060  edm::ParameterSetDescription desc;
1061  desc.add<edm::InputTag>("hlTriggerResults", edm::InputTag("TriggerResults", "", "HLT"));
1062  desc.add<edm::InputTag>("labelEBRecHit", edm::InputTag("ecalRecHit", "EcalRecHitsEB"));
1063  desc.add<edm::InputTag>("labelEERecHit", edm::InputTag("ecalRecHit", "EcalRecHitsEE"));
1064  desc.add<edm::InputTag>("labelHBHERecHit", edm::InputTag("hbhereco"));
1065  desc.add<std::string>("labelVertex", "offlinePrimaryVertices");
1066  desc.add<std::string>("labelMuon", "muons");
1067  std::vector<std::string> trig = {};
1068  desc.add<std::vector<std::string>>("triggers", trig);
1069  desc.add<double>("pMinMuon", 5.0);
1070  desc.addUntracked<int>("verbosity", 0);
1071  desc.add<int>("useRaw", 0);
1072  desc.add<bool>("unCorrect", true);
1073  desc.add<bool>("getCharge", true);
1074  desc.add<bool>("collapseDepth", false);
1075  desc.add<bool>("isItPlan1", false);
1076  desc.addUntracked<bool>("ignoreHECorr", false);
1077  desc.addUntracked<bool>("isItPreRecHit", false);
1078  desc.addUntracked<std::string>("moduleName", "");
1079  desc.addUntracked<std::string>("processName", "");
1080  desc.addUntracked<int>("maxDepth", 7);
1081  desc.addUntracked<std::string>("fileInCorr", "");
1082  desc.addUntracked<bool>("writeRespCorr", false);
1083  descriptions.add("hcalHBHEMuon", desc);
1084 }
HLT_2022v15_cff.InputTag
InputTag
Definition: HLT_2022v15_cff.py:60426
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
ClusterTask_cfi.trig
trig
Definition: ClusterTask_cfi.py:398
edm::ConfigurationDescriptions::add
void add(std::string const &label, ParameterSetDescription const &psetDescription)
Definition: ConfigurationDescriptions.cc:57

◆ gainFactor()

double HcalHBHEMuonAnalyzer::gainFactor ( const HcalDbService *  dbserv,
const HcalDetId &  id 
)
private

Definition at line 1235 of file HcalHBHEMuonAnalyzer.cc.

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

Referenced by analyze().

1235  {
1236  double gain(0.0);
1237  const HcalCalibrations& calibs = conditions->getHcalCalibrations(id);
1238  for (int capid = 0; capid < 4; ++capid)
1239  gain += (0.25 * calibs.respcorrgain(capid));
1240  return gain;
1241 }
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 HcalHBHEMuonAnalyzer::goodCell ( const HcalDetId &  hcid,
const reco::Track *  pTrack,
const CaloGeometry *  geo,
const MagneticField *  bField 
)
private

Definition at line 1258 of file HcalHBHEMuonAnalyzer.cc.

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

Referenced by analyze().

1261  {
1262  std::pair<double, double> rz = hdc_->getRZ(hcid);
1263  bool typeRZ = (hcid.subdet() == HcalEndcap) ? false : true;
1264  bool match = spr::propagateHCAL(pTrack, geo, bField, typeRZ, rz, (((verbosity_ / 10000) % 10) > 0));
1265  return match;
1266 }
HcalDDDRecConstants::getRZ
double getRZ(const int &subdet, const int &ieta, const int &depth) const
Definition: HcalDDDRecConstants.cc:422
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
HcalHBHEMuonAnalyzer::hdc_
const HcalDDDRecConstants * hdc_
Definition: HcalHBHEMuonAnalyzer.cc:123
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 HcalHBHEMuonAnalyzer::isGoodVertex ( const reco::Vertex &  vtx)
private

Definition at line 1211 of file HcalHBHEMuonAnalyzer.cc.

References extraflags_cff::vtx.

Referenced by analyze().

1211  {
1212  if (vtx.isFake())
1213  return false;
1214  if (vtx.ndof() < 4)
1215  return false;
1216  if (vtx.position().Rho() > 2.)
1217  return false;
1218  if (fabs(vtx.position().Z()) > 24.)
1219  return false;
1220  return true;
1221 }
extraflags_cff.vtx
vtx
Definition: extraflags_cff.py:19

◆ matchId()

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

Definition at line 1169 of file HcalHBHEMuonAnalyzer.cc.

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

Referenced by analyze().

1169  {
1170  HcalDetId kd1(id1.subdet(), id1.ieta(), id1.iphi(), 1);
1171  HcalDetId kd2(id1.subdet(), id2.ieta(), id2.iphi(), 1);
1172  int match = ((kd1 == kd2) ? 1 : 0);
1173  return match;
1174 }
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 HcalHBHEMuonAnalyzer::respCorr ( const DetId &  id)
private

Definition at line 1223 of file HcalHBHEMuonAnalyzer.cc.

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

Referenced by analyze().

1223  {
1224  double cfac(1.0);
1225  if (useMyCorr_) {
1226  auto itr = corrValue_.find(id);
1227  if (itr != corrValue_.end())
1228  cfac = itr->second;
1229  } else if (respCorrs_ != nullptr) {
1230  cfac = (respCorrs_->getValues(id))->getValue();
1231  }
1232  return cfac;
1233 }
HcalCondObjectContainer::getValues
const Item * getValues(DetId fId, bool throwOnFail=true) const
Definition: HcalCondObjectContainer.h:159
HcalHBHEMuonAnalyzer::corrValue_
std::map< DetId, double > corrValue_
Definition: HcalHBHEMuonAnalyzer.cc:171
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> HcalHBHEMuonAnalyzer::actHB
private

Definition at line 170 of file HcalHBHEMuonAnalyzer.cc.

Referenced by activeLength(), and beginRun().

◆ actHE

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

Definition at line 170 of file HcalHBHEMuonAnalyzer.cc.

Referenced by activeLength(), and beginRun().

◆ all_triggers_

std::vector<std::string> HcalHBHEMuonAnalyzer::all_triggers_
private

Definition at line 167 of file HcalHBHEMuonAnalyzer.cc.

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

◆ bxNumber_

unsigned int HcalHBHEMuonAnalyzer::bxNumber_
private

Definition at line 134 of file HcalHBHEMuonAnalyzer.cc.

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

◆ chiGlobal_

double HcalHBHEMuonAnalyzer::chiGlobal_
private

Definition at line 149 of file HcalHBHEMuonAnalyzer.cc.

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

◆ chiTracker_

double HcalHBHEMuonAnalyzer::chiTracker_
private

Definition at line 142 of file HcalHBHEMuonAnalyzer.cc.

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

◆ collapseDepth_

const bool HcalHBHEMuonAnalyzer::collapseDepth_
private

Definition at line 100 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze(), and HcalHBHEMuonAnalyzer().

◆ corrValue_

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

Definition at line 171 of file HcalHBHEMuonAnalyzer.cc.

Referenced by HcalHBHEMuonAnalyzer(), and respCorr().

◆ depthMax_

const int HcalHBHEMuonAnalyzer::depthMax_ = 7
staticprivate

Definition at line 132 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze().

◆ dxyTracker_

double HcalHBHEMuonAnalyzer::dxyTracker_
private

Definition at line 142 of file HcalHBHEMuonAnalyzer.cc.

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

◆ dzTracker_

double HcalHBHEMuonAnalyzer::dzTracker_
private

Definition at line 142 of file HcalHBHEMuonAnalyzer.cc.

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

◆ ecal3x3Energy_

double HcalHBHEMuonAnalyzer::ecal3x3Energy_
private

Definition at line 153 of file HcalHBHEMuonAnalyzer.cc.

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

◆ ecalDetId_

unsigned int HcalHBHEMuonAnalyzer::ecalDetId_
private

Definition at line 154 of file HcalHBHEMuonAnalyzer.cc.

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

◆ ecalEnergy_

double HcalHBHEMuonAnalyzer::ecalEnergy_
private

Definition at line 151 of file HcalHBHEMuonAnalyzer.cc.

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

◆ ehcalDetId_

unsigned int HcalHBHEMuonAnalyzer::ehcalDetId_
private

Definition at line 154 of file HcalHBHEMuonAnalyzer.cc.

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

◆ energyMuon_

double HcalHBHEMuonAnalyzer::energyMuon_
private

Definition at line 138 of file HcalHBHEMuonAnalyzer.cc.

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

◆ etaGlob_

double HcalHBHEMuonAnalyzer::etaGlob_
private

Definition at line 138 of file HcalHBHEMuonAnalyzer.cc.

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

◆ eventNumber_

unsigned int HcalHBHEMuonAnalyzer::eventNumber_
private

Definition at line 134 of file HcalHBHEMuonAnalyzer.cc.

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

◆ fileInCorr_

const std::string HcalHBHEMuonAnalyzer::fileInCorr_
private

Definition at line 96 of file HcalHBHEMuonAnalyzer.cc.

Referenced by HcalHBHEMuonAnalyzer().

◆ geo_

const CaloGeometry* HcalHBHEMuonAnalyzer::geo_
private

Definition at line 125 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze(), and beginRun().

◆ getCharge_

const bool HcalHBHEMuonAnalyzer::getCharge_
private

Definition at line 102 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze(), and HcalHBHEMuonAnalyzer().

◆ globalMuonHits_

int HcalHBHEMuonAnalyzer::globalMuonHits_
private

Definition at line 148 of file HcalHBHEMuonAnalyzer.cc.

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

◆ globalTrack_

bool HcalHBHEMuonAnalyzer::globalTrack_
private

Definition at line 141 of file HcalHBHEMuonAnalyzer.cc.

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

◆ globalTrckEta_

double HcalHBHEMuonAnalyzer::globalTrckEta_
private

Definition at line 147 of file HcalHBHEMuonAnalyzer.cc.

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

◆ globalTrckPhi_

double HcalHBHEMuonAnalyzer::globalTrckPhi_
private

Definition at line 147 of file HcalHBHEMuonAnalyzer.cc.

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

◆ globalTrckPt_

double HcalHBHEMuonAnalyzer::globalTrckPt_
private

Definition at line 147 of file HcalHBHEMuonAnalyzer.cc.

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

◆ goodVertex_

unsigned int HcalHBHEMuonAnalyzer::goodVertex_
private

Definition at line 135 of file HcalHBHEMuonAnalyzer.cc.

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

◆ hcal1x1Energy_

double HcalHBHEMuonAnalyzer::hcal1x1Energy_
private

Definition at line 153 of file HcalHBHEMuonAnalyzer.cc.

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

◆ hcal_ieta_

int HcalHBHEMuonAnalyzer::hcal_ieta_
private

Definition at line 155 of file HcalHBHEMuonAnalyzer.cc.

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

◆ hcal_iphi_

int HcalHBHEMuonAnalyzer::hcal_iphi_
private

Definition at line 155 of file HcalHBHEMuonAnalyzer.cc.

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

◆ hcalActiveLength_

double HcalHBHEMuonAnalyzer::hcalActiveLength_
private

Definition at line 166 of file HcalHBHEMuonAnalyzer.cc.

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

◆ hcalActiveLengthHot_

double HcalHBHEMuonAnalyzer::hcalActiveLengthHot_
private

Definition at line 166 of file HcalHBHEMuonAnalyzer.cc.

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

◆ hcalDepthActiveLength_

double HcalHBHEMuonAnalyzer::hcalDepthActiveLength_[depthMax_]
private

Definition at line 157 of file HcalHBHEMuonAnalyzer.cc.

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

◆ hcalDepthActiveLengthHot_

double HcalHBHEMuonAnalyzer::hcalDepthActiveLengthHot_[depthMax_]
private

Definition at line 159 of file HcalHBHEMuonAnalyzer.cc.

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

◆ hcalDepthChargeHot_

double HcalHBHEMuonAnalyzer::hcalDepthChargeHot_[depthMax_]
private

Definition at line 160 of file HcalHBHEMuonAnalyzer.cc.

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

◆ hcalDepthChargeHotBG_

double HcalHBHEMuonAnalyzer::hcalDepthChargeHotBG_[depthMax_]
private

Definition at line 161 of file HcalHBHEMuonAnalyzer.cc.

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

◆ hcalDepthEnergy_

double HcalHBHEMuonAnalyzer::hcalDepthEnergy_[depthMax_]
private

Definition at line 156 of file HcalHBHEMuonAnalyzer.cc.

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

◆ hcalDepthEnergyCorr_

double HcalHBHEMuonAnalyzer::hcalDepthEnergyCorr_[depthMax_]
private

Definition at line 162 of file HcalHBHEMuonAnalyzer.cc.

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

◆ hcalDepthEnergyHot_

double HcalHBHEMuonAnalyzer::hcalDepthEnergyHot_[depthMax_]
private

Definition at line 158 of file HcalHBHEMuonAnalyzer.cc.

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

◆ hcalDepthEnergyHotCorr_

double HcalHBHEMuonAnalyzer::hcalDepthEnergyHotCorr_[depthMax_]
private

Definition at line 163 of file HcalHBHEMuonAnalyzer.cc.

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

◆ hcalDepthMatch_

bool HcalHBHEMuonAnalyzer::hcalDepthMatch_[depthMax_]
private

Definition at line 164 of file HcalHBHEMuonAnalyzer.cc.

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

◆ hcalDepthMatchHot_

bool HcalHBHEMuonAnalyzer::hcalDepthMatchHot_[depthMax_]
private

Definition at line 165 of file HcalHBHEMuonAnalyzer.cc.

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

◆ hcalDetId_

unsigned int HcalHBHEMuonAnalyzer::hcalDetId_
private

Definition at line 154 of file HcalHBHEMuonAnalyzer.cc.

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

◆ hcalEnergy_

double HcalHBHEMuonAnalyzer::hcalEnergy_
private

Definition at line 151 of file HcalHBHEMuonAnalyzer.cc.

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

◆ hcalHot_

bool HcalHBHEMuonAnalyzer::hcalHot_
private

Definition at line 152 of file HcalHBHEMuonAnalyzer.cc.

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

◆ hdc_

const HcalDDDRecConstants* HcalHBHEMuonAnalyzer::hdc_
private

Definition at line 123 of file HcalHBHEMuonAnalyzer.cc.

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

◆ hltConfig_

HLTConfigProvider HcalHBHEMuonAnalyzer::hltConfig_
private

Definition at line 93 of file HcalHBHEMuonAnalyzer.cc.

Referenced by beginRun().

◆ hltresults_

std::vector<int> HcalHBHEMuonAnalyzer::hltresults_
private

Definition at line 168 of file HcalHBHEMuonAnalyzer.cc.

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

◆ hlTriggerResults_

const edm::InputTag HcalHBHEMuonAnalyzer::hlTriggerResults_
private

Definition at line 94 of file HcalHBHEMuonAnalyzer.cc.

Referenced by HcalHBHEMuonAnalyzer().

◆ hoEnergy_

double HcalHBHEMuonAnalyzer::hoEnergy_
private

Definition at line 151 of file HcalHBHEMuonAnalyzer.cc.

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

◆ ignoreHECorr_

const bool HcalHBHEMuonAnalyzer::ignoreHECorr_
private

Definition at line 101 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze(), and HcalHBHEMuonAnalyzer().

◆ innerTrack_

bool HcalHBHEMuonAnalyzer::innerTrack_
private

Definition at line 141 of file HcalHBHEMuonAnalyzer.cc.

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

◆ innerTracketa_

double HcalHBHEMuonAnalyzer::innerTracketa_
private

Definition at line 143 of file HcalHBHEMuonAnalyzer.cc.

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

◆ innerTrackphi_

double HcalHBHEMuonAnalyzer::innerTrackphi_
private

Definition at line 143 of file HcalHBHEMuonAnalyzer.cc.

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

◆ innerTrackpt_

double HcalHBHEMuonAnalyzer::innerTrackpt_
private

Definition at line 143 of file HcalHBHEMuonAnalyzer.cc.

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

◆ isItPlan1_

const bool HcalHBHEMuonAnalyzer::isItPlan1_
private

Definition at line 100 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze(), depth16HE(), and HcalHBHEMuonAnalyzer().

◆ isItPreRecHit_

const bool HcalHBHEMuonAnalyzer::isItPreRecHit_
private

Definition at line 101 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze(), depth16HE(), and HcalHBHEMuonAnalyzer().

◆ isolationR03_

double HcalHBHEMuonAnalyzer::isolationR03_
private

Definition at line 150 of file HcalHBHEMuonAnalyzer.cc.

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

◆ isolationR04_

double HcalHBHEMuonAnalyzer::isolationR04_
private

Definition at line 150 of file HcalHBHEMuonAnalyzer.cc.

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

◆ kount_

int HcalHBHEMuonAnalyzer::kount_
private

Definition at line 129 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze(), and HcalHBHEMuonAnalyzer().

◆ labelEBRecHit_

const edm::InputTag HcalHBHEMuonAnalyzer::labelEBRecHit_
private

Definition at line 95 of file HcalHBHEMuonAnalyzer.cc.

Referenced by HcalHBHEMuonAnalyzer().

◆ labelEERecHit_

const edm::InputTag HcalHBHEMuonAnalyzer::labelEERecHit_
private

Definition at line 95 of file HcalHBHEMuonAnalyzer.cc.

Referenced by HcalHBHEMuonAnalyzer().

◆ labelHBHERecHit_

const edm::InputTag HcalHBHEMuonAnalyzer::labelHBHERecHit_
private

Definition at line 95 of file HcalHBHEMuonAnalyzer.cc.

Referenced by HcalHBHEMuonAnalyzer().

◆ labelMuon_

const std::string HcalHBHEMuonAnalyzer::labelMuon_
private

Definition at line 96 of file HcalHBHEMuonAnalyzer.cc.

Referenced by HcalHBHEMuonAnalyzer().

◆ labelVtx_

const std::string HcalHBHEMuonAnalyzer::labelVtx_
private

Definition at line 96 of file HcalHBHEMuonAnalyzer.cc.

Referenced by HcalHBHEMuonAnalyzer().

◆ lumiNumber_

unsigned int HcalHBHEMuonAnalyzer::lumiNumber_
private

Definition at line 134 of file HcalHBHEMuonAnalyzer.cc.

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

◆ matchedId_

bool HcalHBHEMuonAnalyzer::matchedId_
private

Definition at line 152 of file HcalHBHEMuonAnalyzer.cc.

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

◆ matchedStat_

int HcalHBHEMuonAnalyzer::matchedStat_
private

Definition at line 148 of file HcalHBHEMuonAnalyzer.cc.

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

◆ maxDepth_

const int HcalHBHEMuonAnalyzer::maxDepth_
private

Definition at line 103 of file HcalHBHEMuonAnalyzer.cc.

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

◆ mergedDepth_

bool HcalHBHEMuonAnalyzer::mergedDepth_
private

Definition at line 128 of file HcalHBHEMuonAnalyzer.cc.

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

◆ modnam_

const std::string HcalHBHEMuonAnalyzer::modnam_
private

Definition at line 104 of file HcalHBHEMuonAnalyzer.cc.

Referenced by HcalHBHEMuonAnalyzer().

◆ muon_chi2LocalPosition_

float HcalHBHEMuonAnalyzer::muon_chi2LocalPosition_
private

Definition at line 139 of file HcalHBHEMuonAnalyzer.cc.

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

◆ muon_global_

bool HcalHBHEMuonAnalyzer::muon_global_
private

Definition at line 136 of file HcalHBHEMuonAnalyzer.cc.

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

◆ muon_is_good_

bool HcalHBHEMuonAnalyzer::muon_is_good_
private

Definition at line 136 of file HcalHBHEMuonAnalyzer.cc.

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

◆ muon_is_medium_

bool HcalHBHEMuonAnalyzer::muon_is_medium_
private

Definition at line 137 of file HcalHBHEMuonAnalyzer.cc.

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

◆ muon_is_tight_

bool HcalHBHEMuonAnalyzer::muon_is_tight_
private

Definition at line 137 of file HcalHBHEMuonAnalyzer.cc.

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

◆ muon_segComp_

float HcalHBHEMuonAnalyzer::muon_segComp_
private

Definition at line 139 of file HcalHBHEMuonAnalyzer.cc.

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

◆ muon_tracker_

bool HcalHBHEMuonAnalyzer::muon_tracker_
private

Definition at line 136 of file HcalHBHEMuonAnalyzer.cc.

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

◆ muon_trkKink_

float HcalHBHEMuonAnalyzer::muon_trkKink_
private

Definition at line 139 of file HcalHBHEMuonAnalyzer.cc.

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

◆ numPixelLayers_

int HcalHBHEMuonAnalyzer::numPixelLayers_
private

Definition at line 140 of file HcalHBHEMuonAnalyzer.cc.

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

◆ outerTrack_

bool HcalHBHEMuonAnalyzer::outerTrack_
private

Definition at line 141 of file HcalHBHEMuonAnalyzer.cc.

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

◆ outerTrackChi_

double HcalHBHEMuonAnalyzer::outerTrackChi_
private

Definition at line 144 of file HcalHBHEMuonAnalyzer.cc.

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

◆ outerTrackEta_

double HcalHBHEMuonAnalyzer::outerTrackEta_
private

Definition at line 145 of file HcalHBHEMuonAnalyzer.cc.

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

◆ outerTrackHits_

int HcalHBHEMuonAnalyzer::outerTrackHits_
private

Definition at line 146 of file HcalHBHEMuonAnalyzer.cc.

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

◆ outerTrackPhi_

double HcalHBHEMuonAnalyzer::outerTrackPhi_
private

Definition at line 145 of file HcalHBHEMuonAnalyzer.cc.

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

◆ outerTrackPt_

double HcalHBHEMuonAnalyzer::outerTrackPt_
private

Definition at line 145 of file HcalHBHEMuonAnalyzer.cc.

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

◆ outerTrackRHits_

int HcalHBHEMuonAnalyzer::outerTrackRHits_
private

Definition at line 146 of file HcalHBHEMuonAnalyzer.cc.

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

◆ phiGlob_

double HcalHBHEMuonAnalyzer::phiGlob_
private

Definition at line 138 of file HcalHBHEMuonAnalyzer.cc.

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

◆ pMinMuon_

const double HcalHBHEMuonAnalyzer::pMinMuon_
private

Definition at line 98 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze(), and HcalHBHEMuonAnalyzer().

◆ pMuon_

double HcalHBHEMuonAnalyzer::pMuon_
private

Definition at line 138 of file HcalHBHEMuonAnalyzer.cc.

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

◆ procnm_

const std::string HcalHBHEMuonAnalyzer::procnm_
private

Definition at line 104 of file HcalHBHEMuonAnalyzer.cc.

Referenced by HcalHBHEMuonAnalyzer().

◆ ptGlob_

double HcalHBHEMuonAnalyzer::ptGlob_
private

Definition at line 138 of file HcalHBHEMuonAnalyzer.cc.

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

◆ respCorrs_

HcalRespCorrs* HcalHBHEMuonAnalyzer::respCorrs_
private

Definition at line 126 of file HcalHBHEMuonAnalyzer.cc.

Referenced by beginRun(), and respCorr().

◆ runNumber_

unsigned int HcalHBHEMuonAnalyzer::runNumber_
private

Definition at line 134 of file HcalHBHEMuonAnalyzer.cc.

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

◆ theHBHETopology_

const HcalTopology* HcalHBHEMuonAnalyzer::theHBHETopology_
private

Definition at line 124 of file HcalHBHEMuonAnalyzer.cc.

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

◆ tight_LongPara_

double HcalHBHEMuonAnalyzer::tight_LongPara_
private

Definition at line 149 of file HcalHBHEMuonAnalyzer.cc.

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

◆ tight_PixelHits_

int HcalHBHEMuonAnalyzer::tight_PixelHits_
private

Definition at line 140 of file HcalHBHEMuonAnalyzer.cc.

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

◆ tight_TransImpara_

double HcalHBHEMuonAnalyzer::tight_TransImpara_
private

Definition at line 149 of file HcalHBHEMuonAnalyzer.cc.

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

◆ tight_validFraction_

double HcalHBHEMuonAnalyzer::tight_validFraction_
private

Definition at line 144 of file HcalHBHEMuonAnalyzer.cc.

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

◆ tok_chan_

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

Definition at line 118 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze().

◆ tok_dbservice_

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

Definition at line 121 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze().

◆ tok_ddrec_

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

Definition at line 113 of file HcalHBHEMuonAnalyzer.cc.

Referenced by beginRun().

◆ tok_EB_

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

Definition at line 108 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze().

◆ tok_EE_

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

Definition at line 109 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze().

◆ tok_geom_

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

Definition at line 116 of file HcalHBHEMuonAnalyzer.cc.

Referenced by beginRun().

◆ tok_HBHE_

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

Definition at line 110 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze().

◆ tok_htopo_

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

Definition at line 114 of file HcalHBHEMuonAnalyzer.cc.

Referenced by beginRun().

◆ tok_magField_

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

Definition at line 117 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze().

◆ tok_Muon_

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

Definition at line 111 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze().

◆ tok_respcorr_

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

Definition at line 115 of file HcalHBHEMuonAnalyzer.cc.

Referenced by beginRun().

◆ tok_sevlv_

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

Definition at line 119 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze().

◆ tok_topo_

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

Definition at line 120 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze().

◆ tok_trigRes_

const edm::EDGetTokenT<edm::TriggerResults> HcalHBHEMuonAnalyzer::tok_trigRes_
private

Definition at line 106 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze().

◆ tok_Vtx_

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

Definition at line 107 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze().

◆ trackerLayer_

int HcalHBHEMuonAnalyzer::trackerLayer_
private

Definition at line 140 of file HcalHBHEMuonAnalyzer.cc.

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

◆ tree_

TTree* HcalHBHEMuonAnalyzer::tree_
private

Definition at line 133 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ triggers_

const std::vector<std::string> HcalHBHEMuonAnalyzer::triggers_
private

Definition at line 97 of file HcalHBHEMuonAnalyzer.cc.

Referenced by beginRun().

◆ unCorrect_

const bool HcalHBHEMuonAnalyzer::unCorrect_
private

Definition at line 100 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze(), and HcalHBHEMuonAnalyzer().

◆ useMyCorr_

bool HcalHBHEMuonAnalyzer::useMyCorr_
private

Definition at line 128 of file HcalHBHEMuonAnalyzer.cc.

Referenced by HcalHBHEMuonAnalyzer(), and respCorr().

◆ useRaw_

const int HcalHBHEMuonAnalyzer::useRaw_
private

Definition at line 99 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze(), and HcalHBHEMuonAnalyzer().

◆ verbosity_

const int HcalHBHEMuonAnalyzer::verbosity_
private

Definition at line 99 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze(), and goodCell().

◆ writeRespCorr_

const bool HcalHBHEMuonAnalyzer::writeRespCorr_
private

Definition at line 102 of file HcalHBHEMuonAnalyzer.cc.

Referenced by beginRun().

Generated for CMSSW Reference Manual by doxygen 1.8.14