HcalHBHEMuonAnalyzer Class Reference

Inheritance diagram for HcalHBHEMuonAnalyzer:

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
ESResolverIndex const *	esGetTokenIndices (edm::Transition iTrans) const
std::vector< ESResolverIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const
std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (edm::Transition iTrans) const
ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesWhoseProductsAreConsumed (std::array< std::vector< ModuleDescription const *> *, NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const *> const &labelsToDesc, std::string const &processName) const
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
EDConsumerBase &	operator= (EDConsumerBase &&)=default
bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const
void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)
ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
void	updateLookup (eventsetup::ESRecordsToProductResolverIndices const &)
virtual	~EDConsumerBase () noexcept(false)

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

Private Member Functions
double	activeLength (const DetId &, bool)
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::ESGetToken< EcalPFRecHitThresholds, EcalPFRecHitThresholdsRcd >	tok_ecalPFRecHitThresholds_
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 bool	usePFThresh_
const int	useRaw_
const int	verbosity_
const bool	writeRespCorr_

Static Private Attributes
static const int	depthMax_ = 7

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

Detailed Description

Definition at line 74 of file HcalHBHEMuonAnalyzer.cc.

Constructor & Destructor Documentation

HcalHBHEMuonAnalyzer()

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

explicit

Definition at line 179 of file HcalHBHEMuonAnalyzer.cc.

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

    : hlTriggerResults_(iConfig.getParameter<edm::InputTag>("hlTriggerResults")),
      labelEBRecHit_(iConfig.getParameter<edm::InputTag>("labelEBRecHit")),
      labelEERecHit_(iConfig.getParameter<edm::InputTag>("labelEERecHit")),
      labelHBHERecHit_(iConfig.getParameter<edm::InputTag>("labelHBHERecHit")),
      labelVtx_(iConfig.getParameter<std::string>("labelVertex")),
      labelMuon_(iConfig.getParameter<std::string>("labelMuon")),
      fileInCorr_(iConfig.getUntrackedParameter<std::string>("fileInCorr", "")),
      triggers_(iConfig.getParameter<std::vector<std::string>>("triggers")),
      pMinMuon_(iConfig.getParameter<double>("pMinMuon")),
      verbosity_(iConfig.getUntrackedParameter<int>("verbosity", 0)),
      useRaw_(iConfig.getParameter<int>("useRaw")),
      unCorrect_(iConfig.getParameter<bool>("unCorrect")),
      collapseDepth_(iConfig.getParameter<bool>("collapseDepth")),
      isItPlan1_(iConfig.getParameter<bool>("isItPlan1")),
      ignoreHECorr_(iConfig.getUntrackedParameter<bool>("ignoreHECorr", false)),
      isItPreRecHit_(iConfig.getUntrackedParameter<bool>("isItPreRecHit", false)),
      getCharge_(iConfig.getParameter<bool>("getCharge")),
      writeRespCorr_(iConfig.getUntrackedParameter<bool>("writeRespCorr", false)),
      maxDepth_(iConfig.getUntrackedParameter<int>("maxDepth", 7)),
      modnam_(iConfig.getUntrackedParameter<std::string>("moduleName", "")),
      procnm_(iConfig.getUntrackedParameter<std::string>("processName", "")),
      usePFThresh_(iConfig.getParameter<bool>("usePFThreshold")),
      tok_trigRes_(consumes<edm::TriggerResults>(hlTriggerResults_)),
      tok_Vtx_((modnam_.empty()) ? consumes<reco::VertexCollection>(labelVtx_)
                                 : consumes<reco::VertexCollection>(edm::InputTag(modnam_, labelVtx_, procnm_))),
      tok_EB_(consumes<EcalRecHitCollection>(labelEBRecHit_)),
      tok_EE_(consumes<EcalRecHitCollection>(labelEERecHit_)),
      tok_HBHE_(consumes<HBHERecHitCollection>(labelHBHERecHit_)),
      tok_Muon_((modnam_.empty()) ? consumes<reco::MuonCollection>(labelMuon_)
                                  : consumes<reco::MuonCollection>(edm::InputTag(modnam_, labelMuon_, procnm_))),
      tok_ddrec_(esConsumes<HcalDDDRecConstants, HcalRecNumberingRecord, edm::Transition::BeginRun>()),
      tok_htopo_(esConsumes<HcalTopology, HcalRecNumberingRecord, edm::Transition::BeginRun>()),
      tok_respcorr_(esConsumes<HcalRespCorrs, HcalRespCorrsRcd, edm::Transition::BeginRun>()),
      tok_geom_(esConsumes<CaloGeometry, CaloGeometryRecord, edm::Transition::BeginRun>()),
      tok_magField_(esConsumes<MagneticField, IdealMagneticFieldRecord>()),
      tok_chan_(esConsumes<EcalChannelStatus, EcalChannelStatusRcd>()),
      tok_sevlv_(esConsumes<EcalSeverityLevelAlgo, EcalSeverityLevelAlgoRcd>()),
      tok_topo_(esConsumes<CaloTopology, CaloTopologyRecord>()),
      tok_dbservice_(esConsumes<HcalDbService, HcalDbRecord>()),
      tok_ecalPFRecHitThresholds_(esConsumes<EcalPFRecHitThresholds, EcalPFRecHitThresholdsRcd>()),
      hdc_(nullptr),
      theHBHETopology_(nullptr),
      respCorrs_(nullptr) {
  usesResource(TFileService::kSharedResource);
  //now do what ever initialization is needed
  kount_ = 0;
  mergedDepth_ = (!isItPreRecHit_) || (collapseDepth_);

  if (modnam_.empty()) {
    edm::LogVerbatim("HBHEMuon") << "Labels used: Trig " << hlTriggerResults_ << " Vtx " << labelVtx_ << " EB "
                                 << labelEBRecHit_ << " EE " << labelEERecHit_ << " HBHE " << labelHBHERecHit_ << " MU "
                                 << labelMuon_;
  } else {
    edm::LogVerbatim("HBHEMuon") << "Labels used Trig " << hlTriggerResults_ << "\n  Vtx  "
                                 << edm::InputTag(modnam_, labelVtx_, procnm_) << "\n  EB   " << labelEBRecHit_
                                 << "\n  EE   " << labelEERecHit_ << "\n  HBHE " << labelHBHERecHit_ << "\n  MU   "
                                 << edm::InputTag(modnam_, labelMuon_, procnm_);
  }

  if (!fileInCorr_.empty()) {
    std::ifstream infile(fileInCorr_.c_str());
    if (infile.is_open()) {
      while (true) {
        unsigned int id;
        double cfac;
        infile >> id >> cfac;
        if (!infile.good())
          break;
        corrValue_[DetId(id)] = cfac;
      }
      infile.close();
    }
  }
  useMyCorr_ = (!corrValue_.empty());
  edm::LogVerbatim("HBHEMuon") << "Flags used: UseRaw " << useRaw_ << " GetCharge " << getCharge_ << " UnCorrect "
                               << unCorrect_ << " IgnoreHECorr " << ignoreHECorr_ << " CollapseDepth " << collapseDepth_
                               << ":" << mergedDepth_ << " IsItPlan1 " << isItPlan1_ << " IsItPreRecHit "
                               << isItPreRecHit_ << " UseMyCorr " << useMyCorr_ << " pMinMuon " << pMinMuon_
                               << " usePFThresh " << usePFThresh_;
}

Member Function Documentation

activeLength()

double HcalHBHEMuonAnalyzer::activeLength ( const DetId &  id_, bool  debug  )

private

Definition at line 1209 of file HcalHBHEMuonAnalyzer.cc.

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

Referenced by analyze(), and beginRun().

                                                                      {
  HcalDetId id(id_);
  int ieta = id.ietaAbs();
  int zside = id.zside();
  int iphi = id.iphi();
  std::vector<int> dpths;
  if (mergedDepth_) {
    std::vector<HcalDetId> ids;
    hdc_->unmergeDepthDetId(id, ids);
    for (auto idh : ids)
      dpths.emplace_back(idh.depth());
  } else {
    dpths.emplace_back(id.depth());
  }
  double lx(0);
  if (debug) {
    std::ostringstream st1;
    st1 << "Subdet:zside:ieta:iphi:ndepth " << id.subdet() << " : " << zside << " : " << ieta << " : " << iphi << " : "
        << dpths.size() << " depths";
    for (unsigned int k = 0; k < dpths.size(); ++k)
      st1 << " : " << dpths[k];
    edm::LogVerbatim("HBHEMuon") << st1.str();
  }
  if (id.subdet() == HcalBarrel) {
    for (unsigned int i = 0; i < actHB.size(); ++i) {
      if ((ieta == actHB[i].ieta) && (zside == actHB[i].zside) &&
          (std::find(dpths.begin(), dpths.end(), actHB[i].depth) != dpths.end()) &&
          (std::find(actHB[i].iphis.begin(), actHB[i].iphis.end(), iphi) != actHB[i].iphis.end())) {
        lx += actHB[i].thick;
        if (debug) {
          edm::LogVerbatim("HBHEMuon") << "Eta: " << (ieta == actHB[i].ieta) << " zside " << (zside == actHB[i].zside)
                                       << " Depth " << actHB[i].depth << ":"
                                       << (std::find(dpths.begin(), dpths.end(), actHB[i].depth) != dpths.end())
                                       << " Phi: "
                                       << (std::find(actHB[i].iphis.begin(), actHB[i].iphis.end(), iphi) !=
                                           actHB[i].iphis.end())
                                       << " Length: " << actHB[i].thick << " : " << lx;
        }
      }
    }
  } else {
    for (unsigned int i = 0; i < actHE.size(); ++i) {
      if ((ieta == actHE[i].ieta) && (zside == actHE[i].zside) &&
          (std::find(dpths.begin(), dpths.end(), actHE[i].depth) != dpths.end()) &&
          (std::find(actHE[i].iphis.begin(), actHE[i].iphis.end(), iphi) != actHE[i].iphis.end())) {
        lx += actHE[i].thick;
        if (debug) {
          edm::LogVerbatim("HBHEMuon") << "Eta: " << (ieta == actHE[i].ieta) << " zside " << (zside == actHE[i].zside)
                                       << " Depth " << actHE[i].depth << ":"
                                       << (std::find(dpths.begin(), dpths.end(), actHE[i].depth) != dpths.end())
                                       << " Phi: "
                                       << (std::find(actHE[i].iphis.begin(), actHE[i].iphis.end(), iphi) !=
                                           actHE[i].iphis.end())
                                       << " Length: " << actHE[i].thick << " : " << lx;
        }
      }
    }
  }
  return lx;
}

analyze()

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

overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 266 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, hcalRecHitTable_cff::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, hcalRecHitTable_cff::ieta, HcalDetId::ieta(), iEvent, ignoreHECorr_, info(), susybsm::HSCParticleType::innerTrack, innerTrack_, innerTracketa_, innerTrackphi_, innerTrackpt_, hcalRecHitTable_cff::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_ecalPFRecHitThresholds_, tok_EE_, tok_HBHE_, tok_magField_, tok_Muon_, tok_sevlv_, tok_topo_, tok_trigRes_, tok_Vtx_, trackerLayer_, tree_, edm::TriggerNames::triggerIndex(), unCorrect_, usePFThresh_, useRaw_, HcalTopology::validHcal(), verbosity_, and L1BJetProducer_cff::vtx.

                                                                                      {
  ++kount_;
  clearVectors();
  std::vector<bool> muon_is_good, muon_global, muon_tracker;
  std::vector<bool> muon_is_tight, muon_is_medium;
  std::vector<double> ptGlob, etaGlob, phiGlob, energyMuon, pMuon;
  std::vector<float> muon_trkKink, muon_chi2LocalPosition, muon_segComp;
  std::vector<int> trackerLayer, numPixelLayers, tight_PixelHits;
  std::vector<bool> innerTrack, outerTrack, globalTrack;
  std::vector<double> chiTracker, dxyTracker, dzTracker;
  std::vector<double> innerTrackpt, innerTracketa, innerTrackphi;
  std::vector<double> tight_validFraction, outerTrackChi;
  std::vector<double> outerTrackPt, outerTrackEta, outerTrackPhi;
  std::vector<int> outerTrackHits, outerTrackRHits;
  std::vector<double> globalTrckPt, globalTrckEta, globalTrckPhi;
  std::vector<int> globalMuonHits, matchedStat;
  std::vector<double> chiGlobal, tight_LongPara, tight_TransImpara;
  std::vector<double> isolationR04, isolationR03;
  std::vector<double> ecalEnergy, hcalEnergy, hoEnergy;
  std::vector<bool> matchedId, hcalHot;
  std::vector<double> ecal3x3Energy, hcal1x1Energy;
  std::vector<unsigned int> ecalDetId, hcalDetId, ehcalDetId;
  std::vector<int> hcal_ieta, hcal_iphi;
  std::vector<double> hcalDepthEnergy[depthMax_];
  std::vector<double> hcalDepthActiveLength[depthMax_];
  std::vector<double> hcalDepthEnergyHot[depthMax_];
  std::vector<double> hcalDepthActiveLengthHot[depthMax_];
  std::vector<double> hcalDepthChargeHot[depthMax_];
  std::vector<double> hcalDepthChargeHotBG[depthMax_];
  std::vector<double> hcalDepthEnergyCorr[depthMax_];
  std::vector<double> hcalDepthEnergyHotCorr[depthMax_];
  std::vector<bool> hcalDepthMatch[depthMax_];
  std::vector<bool> hcalDepthMatchHot[depthMax_];
  std::vector<double> hcalActiveLength, hcalActiveLengthHot;
  runNumber_ = iEvent.id().run();
  eventNumber_ = iEvent.id().event();
  lumiNumber_ = iEvent.id().luminosityBlock();
  bxNumber_ = iEvent.bunchCrossing();
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HBHEMuon") << "Run " << runNumber_ << " Event " << eventNumber_ << " Lumi " << lumiNumber_ << " BX "
                               << bxNumber_ << std::endl;
#endif
  const edm::Handle<edm::TriggerResults>& _Triggers = iEvent.getHandle(tok_trigRes_);
#ifdef EDM_ML_DEBUG
  if ((verbosity_ / 10000) % 10 > 0)
    edm::LogVerbatim("HBHEMuon") << "Size of all triggers " << all_triggers_.size();
#endif
  int Ntriggers = static_cast<int>(all_triggers_.size());
#ifdef EDM_ML_DEBUG
  if ((verbosity_ / 10000) % 10 > 0)
    edm::LogVerbatim("HBHEMuon") << "Size of HLT MENU: " << _Triggers->size();
#endif
  if (_Triggers.isValid()) {
    const edm::TriggerNames& triggerNames_ = iEvent.triggerNames(*_Triggers);
    std::vector<int> index;
    for (int i = 0; i < Ntriggers; i++) {
      index.push_back(triggerNames_.triggerIndex(all_triggers_[i]));
      int triggerSize = static_cast<int>(_Triggers->size());
#ifdef EDM_ML_DEBUG
      if ((verbosity_ / 10000) % 10 > 0)
        edm::LogVerbatim("HBHEMuon") << "outside loop " << index[i] << "\ntriggerSize " << triggerSize;
#endif
      if (index[i] < triggerSize) {
        hltresults_.push_back(_Triggers->accept(index[i]));
#ifdef EDM_ML_DEBUG
        if ((verbosity_ / 10000) % 10 > 0)
          edm::LogVerbatim("HBHEMuon") << "Trigger_info " << triggerSize << " triggerSize " << index[i]
                                       << " trigger_index " << hltresults_.at(i) << " hltresult";
#endif
      } else {
        if ((verbosity_ / 10000) % 10 > 0)
          edm::LogVerbatim("HBHEMuon") << "Requested HLT path \""
                                       << "\" does not exist";
      }
    }
  }

  // get handles to calogeometry and calotopology
  const MagneticField* bField = &iSetup.getData(tok_magField_);
  const EcalChannelStatus* theEcalChStatus = &iSetup.getData(tok_chan_);
  const EcalSeverityLevelAlgo* sevlv = &iSetup.getData(tok_sevlv_);
  const CaloTopology* caloTopology = &iSetup.getData(tok_topo_);
  const HcalDbService* conditions = &iSetup.getData(tok_dbservice_);
  const EcalPFRecHitThresholds* eThresholds = &iSetup.getData(tok_ecalPFRecHitThresholds_);

  // Relevant blocks from iEvent
  const edm::Handle<reco::VertexCollection>& vtx = iEvent.getHandle(tok_Vtx_);

  edm::Handle<EcalRecHitCollection> barrelRecHitsHandle = iEvent.getHandle(tok_EB_);
  edm::Handle<EcalRecHitCollection> endcapRecHitsHandle = iEvent.getHandle(tok_EE_);

  edm::Handle<HBHERecHitCollection> hbhe = iEvent.getHandle(tok_HBHE_);

  const edm::Handle<reco::MuonCollection>& _Muon = iEvent.getHandle(tok_Muon_);

  // require a good vertex
  math::XYZPoint pvx;
  goodVertex_ = 0;
  if (!vtx.isValid()) {
#ifdef EDM_ML_DEBUG
    edm::LogVerbatim("HBHEMuon") << "No Good Vertex found == Reject";
#endif
    return;
  }
  reco::VertexCollection::const_iterator firstGoodVertex = vtx->end();
  for (reco::VertexCollection::const_iterator it = vtx->begin(); it != vtx->end(); it++) {
    if (isGoodVertex(*it)) {
      if (firstGoodVertex == vtx->end())
        firstGoodVertex = it;
      ++goodVertex_;
    }
  }
  if (firstGoodVertex != vtx->end())
    pvx = firstGoodVertex->position();

  bool accept(false);
  if (_Muon.isValid() && barrelRecHitsHandle.isValid() && endcapRecHitsHandle.isValid() && hbhe.isValid()) {
    for (const auto& RecMuon : (*(_Muon.product()))) {
      muon_is_good.push_back(RecMuon.isPFMuon());
      muon_global.push_back(RecMuon.isGlobalMuon());
      muon_tracker.push_back(RecMuon.isTrackerMuon());
      ptGlob.push_back(RecMuon.pt());
      etaGlob.push_back(RecMuon.eta());
      phiGlob.push_back(RecMuon.phi());
      energyMuon.push_back(RecMuon.energy());
      pMuon.push_back(RecMuon.p());
#ifdef EDM_ML_DEBUG
      edm::LogVerbatim("HBHEMuon") << "Energy:" << RecMuon.energy() << " P:" << RecMuon.p();
#endif
      muon_is_tight.push_back(muon::isTightMuon(RecMuon, *firstGoodVertex));
      muon_is_medium.push_back(muon::isMediumMuon(RecMuon));
      muon_trkKink.push_back(RecMuon.combinedQuality().trkKink);
      muon_chi2LocalPosition.push_back(RecMuon.combinedQuality().chi2LocalPosition);
      muon_segComp.push_back(muon::segmentCompatibility(RecMuon));
      // acessing tracker hits info
      if (RecMuon.track().isNonnull()) {
        trackerLayer.push_back(RecMuon.track()->hitPattern().trackerLayersWithMeasurement());
      } else {
        trackerLayer.push_back(-1);
      }
      if (RecMuon.innerTrack().isNonnull()) {
        innerTrack.push_back(true);
        numPixelLayers.push_back(RecMuon.innerTrack()->hitPattern().pixelLayersWithMeasurement());
        chiTracker.push_back(RecMuon.innerTrack()->normalizedChi2());
        dxyTracker.push_back(fabs(RecMuon.innerTrack()->dxy(pvx)));
        dzTracker.push_back(fabs(RecMuon.innerTrack()->dz(pvx)));
        innerTrackpt.push_back(RecMuon.innerTrack()->pt());
        innerTracketa.push_back(RecMuon.innerTrack()->eta());
        innerTrackphi.push_back(RecMuon.innerTrack()->phi());
        tight_PixelHits.push_back(RecMuon.innerTrack()->hitPattern().numberOfValidPixelHits());
        tight_validFraction.push_back(RecMuon.innerTrack()->validFraction());
      } else {
        innerTrack.push_back(false);
        numPixelLayers.push_back(0);
        chiTracker.push_back(0);
        dxyTracker.push_back(0);
        dzTracker.push_back(0);
        innerTrackpt.push_back(0);
        innerTracketa.push_back(0);
        innerTrackphi.push_back(0);
        tight_PixelHits.push_back(0);
        tight_validFraction.push_back(-99);
      }
      // outer track info
      if (RecMuon.outerTrack().isNonnull()) {
        outerTrack.push_back(true);
        outerTrackPt.push_back(RecMuon.outerTrack()->pt());
        outerTrackEta.push_back(RecMuon.outerTrack()->eta());
        outerTrackPhi.push_back(RecMuon.outerTrack()->phi());
        outerTrackChi.push_back(RecMuon.outerTrack()->normalizedChi2());
        outerTrackHits.push_back(RecMuon.outerTrack()->numberOfValidHits());
        outerTrackRHits.push_back(RecMuon.outerTrack()->recHitsSize());
      } else {
        outerTrack.push_back(false);
        outerTrackPt.push_back(0);
        outerTrackEta.push_back(0);
        outerTrackPhi.push_back(0);
        outerTrackChi.push_back(0);
        outerTrackHits.push_back(0);
        outerTrackRHits.push_back(0);
      }
      // Tight Muon cuts
      if (RecMuon.globalTrack().isNonnull()) {
        globalTrack.push_back(true);
        chiGlobal.push_back(RecMuon.globalTrack()->normalizedChi2());
        globalMuonHits.push_back(RecMuon.globalTrack()->hitPattern().numberOfValidMuonHits());
        matchedStat.push_back(RecMuon.numberOfMatchedStations());
        globalTrckPt.push_back(RecMuon.globalTrack()->pt());
        globalTrckEta.push_back(RecMuon.globalTrack()->eta());
        globalTrckPhi.push_back(RecMuon.globalTrack()->phi());
        tight_TransImpara.push_back(fabs(RecMuon.muonBestTrack()->dxy(pvx)));
        tight_LongPara.push_back(fabs(RecMuon.muonBestTrack()->dz(pvx)));
      } else {
        globalTrack.push_back(false);
        chiGlobal.push_back(0);
        globalMuonHits.push_back(0);
        matchedStat.push_back(0);
        globalTrckPt.push_back(0);
        globalTrckEta.push_back(0);
        globalTrckPhi.push_back(0);
        tight_TransImpara.push_back(0);
        tight_LongPara.push_back(0);
      }

      isolationR04.push_back(
          ((RecMuon.pfIsolationR04().sumChargedHadronPt +
            std::max(0.,
                     RecMuon.pfIsolationR04().sumNeutralHadronEt + RecMuon.pfIsolationR04().sumPhotonEt -
                         (0.5 * RecMuon.pfIsolationR04().sumPUPt))) /
           RecMuon.pt()));

      isolationR03.push_back(
          ((RecMuon.pfIsolationR03().sumChargedHadronPt +
            std::max(0.,
                     RecMuon.pfIsolationR03().sumNeutralHadronEt + RecMuon.pfIsolationR03().sumPhotonEt -
                         (0.5 * RecMuon.pfIsolationR03().sumPUPt))) /
           RecMuon.pt()));

      ecalEnergy.push_back(RecMuon.calEnergy().emS9);
      hcalEnergy.push_back(RecMuon.calEnergy().hadS9);
      hoEnergy.push_back(RecMuon.calEnergy().hoS9);

      double eEcal(0), eHcal(0), activeLengthTot(0), activeLengthHotTot(0);
      double eHcalDepth[depthMax_], eHcalDepthHot[depthMax_];
      double eHcalDepthC[depthMax_], eHcalDepthHotC[depthMax_];
      double cHcalDepthHot[depthMax_], cHcalDepthHotBG[depthMax_];
      double activeL[depthMax_], activeHotL[depthMax_];
      bool matchDepth[depthMax_], matchDepthHot[depthMax_];
      HcalDetId eHcalDetId[depthMax_];
      unsigned int isHot(0);
      bool tmpmatch(false);
      int ieta(-1000), iphi(-1000);
      for (int i = 0; i < depthMax_; ++i) {
        eHcalDepth[i] = eHcalDepthHot[i] = 0;
        eHcalDepthC[i] = eHcalDepthHotC[i] = 0;
        cHcalDepthHot[i] = cHcalDepthHotBG[i] = 0;
        activeL[i] = activeHotL[i] = 0;
        matchDepth[i] = matchDepthHot[i] = true;
      }
      if (RecMuon.innerTrack().isNonnull()) {
        const reco::Track* pTrack = (RecMuon.innerTrack()).get();
        spr::propagatedTrackID trackID = spr::propagateCALO(pTrack, geo_, bField, (((verbosity_ / 100) % 10 > 0)));
        if ((RecMuon.p() > pMinMuon_) && (trackID.okHCAL))
          accept = true;

        ecalDetId.push_back((trackID.detIdECAL)());
        hcalDetId.push_back((trackID.detIdHCAL)());
        ehcalDetId.push_back((trackID.detIdEHCAL)());

        HcalDetId check;
        std::pair<bool, HcalDetId> info = spr::propagateHCALBack(pTrack, geo_, bField, (((verbosity_ / 100) % 10 > 0)));
        if (info.first) {
          check = info.second;
        }

        bool okE = trackID.okECAL;
        if (okE) {
          const DetId isoCell(trackID.detIdECAL);
          std::pair<double, bool> e3x3 = (usePFThresh_ ? spr::eECALmatrix(isoCell,
                                                                          barrelRecHitsHandle,
                                                                          endcapRecHitsHandle,
                                                                          *theEcalChStatus,
                                                                          geo_,
                                                                          caloTopology,
                                                                          sevlv,
                                                                          eThresholds,
                                                                          1,
                                                                          1,
                                                                          false)
                                                       : spr::eECALmatrix(isoCell,
                                                                          barrelRecHitsHandle,
                                                                          endcapRecHitsHandle,
                                                                          *theEcalChStatus,
                                                                          geo_,
                                                                          caloTopology,
                                                                          sevlv,
                                                                          1,
                                                                          1,
                                                                          -100.0,
                                                                          -100.0,
                                                                          -500.0,
                                                                          500.0,
                                                                          false));
          eEcal = e3x3.first;
          okE = e3x3.second;
        }
#ifdef EDM_ML_DEBUG
        edm::LogVerbatim("HBHEMuon") << "Propagate Track to ECAL: " << okE << ":" << trackID.okECAL << " E " << eEcal;
#endif

        if (trackID.okHCAL) {
          DetId closestCell(trackID.detIdHCAL);
          HcalDetId hcidt(closestCell.rawId());
          if ((hcidt.ieta() == check.ieta()) && (hcidt.iphi() == check.iphi()))
            tmpmatch = true;
#ifdef EDM_ML_DEBUG
          edm::LogVerbatim("HBHEMuon") << "Front " << hcidt << " Back " << info.first << ":" << check << " Match "
                                       << tmpmatch;
#endif

          HcalSubdetector subdet = hcidt.subdet();
          ieta = hcidt.ieta();
          iphi = hcidt.iphi();
          bool hborhe = (std::abs(ieta) == 16);

          eHcal = spr::eHCALmatrix(theHBHETopology_,
                                   closestCell,
                                   hbhe,
                                   0,
                                   0,
                                   false,
                                   true,
                                   -100.0,
                                   -100.0,
                                   -100.0,
                                   -100.0,
                                   -500.,
                                   500.,
                                   useRaw_);
          std::vector<std::pair<double, int>> ehdepth;
          spr::energyHCALCell((HcalDetId)closestCell,
                              hbhe,
                              ehdepth,
                              maxDepth_,
                              -100.0,
                              -100.0,
                              -100.0,
                              -100.0,
                              -500.0,
                              500.0,
                              useRaw_,
                              depth16HE(ieta, iphi),
                              (((verbosity_ / 1000) % 10) > 0));
          for (int i = 0; i < depthMax_; ++i)
            eHcalDetId[i] = HcalDetId();
          for (unsigned int i = 0; i < ehdepth.size(); ++i) {
            HcalSubdetector subdet0 =
                (hborhe) ? ((ehdepth[i].second >= depth16HE(ieta, iphi)) ? HcalEndcap : HcalBarrel) : subdet;
            HcalDetId hcid0(subdet0, ieta, iphi, ehdepth[i].second);
            double actL = activeLength(DetId(hcid0), ((verbosity_ / 100000) % 10));
            double ene = ehdepth[i].first;
#ifdef EDM_ML_DEBUG
            if ((verbosity_ / 100000) % 10)
              edm::LogVerbatim("HBHEMuon")
                  << "DetId " << subdet0 << ":" << ieta << ":" << iphi << ":" << ehdepth[i].second << " E " << ene;
#endif
            bool tmpC(false);
            if (ene > 0.0) {
              if (!(theHBHETopology_->validHcal(hcid0))) {
                edm::LogWarning("HBHEMuon") << "(1) Invalid ID " << hcid0 << " with E = " << ene;
                edm::LogWarning("HBHEMuon") << HcalDetId(closestCell) << " with " << ehdepth.size() << " depths:";
                for (const auto& ehd : ehdepth)
                  edm::LogWarning("HBHEMuon") << " " << ehd.second << ":" << ehd.first;
              } else {
                tmpC = goodCell(hcid0, pTrack, geo_, bField);
                double enec(ene);
                if (unCorrect_) {
                  double corr = (ignoreHECorr_ && (subdet0 == HcalEndcap)) ? 1.0 : respCorr(DetId(hcid0));
                  if (corr != 0)
                    ene /= corr;
#ifdef EDM_ML_DEBUG
                  HcalDetId id = (isItPlan1_ && isItPreRecHit_) ? hdc_->mergedDepthDetId(hcid0) : hcid0;
                  edm::LogVerbatim("HBHEMuon") << hcid0 << ":" << id << " Corr " << corr;
#endif
                }
                int depth = ehdepth[i].second - 1;
                if (collapseDepth_) {
                  HcalDetId id = hdc_->mergedDepthDetId(hcid0);
                  depth = id.depth() - 1;
                }
                eHcalDepth[depth] += ene;
                eHcalDepthC[depth] += enec;
                activeL[depth] += actL;
                activeLengthTot += actL;
                matchDepth[depth] = (matchDepth[depth] && tmpC);
#ifdef EDM_ML_DEBUG
                if ((verbosity_ % 10) > 0)
                  edm::LogVerbatim("HBHEMuon")
                      << hcid0 << " E " << ene << ":" << enec << " L " << actL << " Match " << tmpC;
#endif
              }
            }
          }
#ifdef EDM_ML_DEBUG
          if ((verbosity_ % 10) > 0) {
            edm::LogVerbatim("HBHEMuon") << hcidt << " Match " << tmpmatch << " Depths " << ehdepth.size();
            for (unsigned int k = 0; k < ehdepth.size(); ++k)
              edm::LogVerbatim("HBHEMuon") << " [" << k << ":" << ehdepth[k].second << "] Match " << matchDepth[k]
                                           << " E " << eHcalDepth[k] << ":" << eHcalDepthC[k] << " L " << activeL[k];
          }
#endif
          HcalDetId hotCell;
          spr::eHCALmatrix(geo_, theHBHETopology_, closestCell, hbhe, 1, 1, hotCell, false, useRaw_, false);
          isHot = matchId(closestCell, hotCell);
          if (hotCell != HcalDetId()) {
            subdet = HcalDetId(hotCell).subdet();
            ieta = HcalDetId(hotCell).ieta();
            iphi = HcalDetId(hotCell).iphi();
            hborhe = (std::abs(ieta) == 16);
            std::vector<std::pair<double, int>> ehdepth;
            spr::energyHCALCell(hotCell,
                                hbhe,
                                ehdepth,
                                maxDepth_,
                                -100.0,
                                -100.0,
                                -100.0,
                                -100.0,
                                -500.0,
                                500.0,
                                useRaw_,
                                depth16HE(ieta, iphi),
                                false);  //(((verbosity_/1000)%10)>0    ));
            for (int i = 0; i < depthMax_; ++i)
              eHcalDetId[i] = HcalDetId();
            for (unsigned int i = 0; i < ehdepth.size(); ++i) {
              HcalSubdetector subdet0 =
                  (hborhe) ? ((ehdepth[i].second >= depth16HE(ieta, iphi)) ? HcalEndcap : HcalBarrel) : subdet;
              HcalDetId hcid0(subdet0, ieta, iphi, ehdepth[i].second);
              double actL = activeLength(DetId(hcid0), ((verbosity_ / 100000) % 10));
              double ene = ehdepth[i].first;
#ifdef EDM_ML_DEBUG
              if ((verbosity_ / 100000) % 10)
                edm::LogVerbatim("HBHEMuon")
                    << "DetId " << subdet0 << ":" << ieta << ":" << iphi << ":" << ehdepth[i].second << " E " << ene;
#endif
              bool tmpC(false);
              if (ene > 0.0) {
                if (!(theHBHETopology_->validHcal(hcid0))) {
                  edm::LogWarning("HBHEMuon") << "(2) Invalid ID " << hcid0 << " with E = " << ene;
                  edm::LogWarning("HBHEMuon") << HcalDetId(hotCell) << " with " << ehdepth.size() << " depths:";
                  for (const auto& ehd : ehdepth)
                    edm::LogWarning("HBHEMuon") << " " << ehd.second << ":" << ehd.first;
                } else {
                  tmpC = goodCell(hcid0, pTrack, geo_, bField);
                  double chg(ene), enec(ene);
                  if (unCorrect_) {
                    double corr = (ignoreHECorr_ && (subdet0 == HcalEndcap)) ? 1.0 : respCorr(DetId(hcid0));
                    if (corr != 0)
                      ene /= corr;
#ifdef EDM_ML_DEBUG
                    HcalDetId id = (isItPlan1_ && isItPreRecHit_) ? hdc_->mergedDepthDetId(hcid0) : hcid0;
                    edm::LogVerbatim("HBHEMuon")
                        << hcid0 << ":" << id << " Corr " << corr << " E " << ene << ":" << enec;
#endif
                  }
                  if (getCharge_) {
                    double gain = gainFactor(conditions, hcid0);
                    if (gain != 0)
                      chg /= gain;
#ifdef EDM_ML_DEBUG
                    edm::LogVerbatim("HBHEMuon") << hcid0 << " Gain " << gain << " C " << chg;
#endif
                  }
                  int depth = ehdepth[i].second - 1;
                  if (collapseDepth_) {
                    HcalDetId id = hdc_->mergedDepthDetId(hcid0);
                    depth = id.depth() - 1;
                  }
                  eHcalDepthHot[depth] += ene;
                  eHcalDepthHotC[depth] += enec;
                  cHcalDepthHot[depth] += chg;
                  activeHotL[depth] += actL;
                  activeLengthHotTot += actL;
                  matchDepthHot[depth] = (matchDepthHot[depth] && tmpC);
#ifdef EDM_ML_DEBUG
                  if ((verbosity_ % 10) > 0)
                    edm::LogVerbatim("HBHEMuon") << hcid0 << " depth " << depth << " E " << ene << ":" << enec << " C "
                                                 << chg << " L " << actL << " Match " << tmpC;
#endif
                }
              }
            }

            HcalDetId oppCell(subdet, -ieta, iphi, HcalDetId(hotCell).depth());
            std::vector<std::pair<double, int>> ehdeptho;
            spr::energyHCALCell(oppCell,
                                hbhe,
                                ehdeptho,
                                maxDepth_,
                                -100.0,
                                -100.0,
                                -100.0,
                                -100.0,
                                -500.0,
                                500.0,
                                useRaw_,
                                depth16HE(-ieta, iphi),
                                false);  //(((verbosity_/1000)%10)>0));
            for (unsigned int i = 0; i < ehdeptho.size(); ++i) {
              HcalSubdetector subdet0 =
                  (hborhe) ? ((ehdeptho[i].second >= depth16HE(-ieta, iphi)) ? HcalEndcap : HcalBarrel) : subdet;
              HcalDetId hcid0(subdet0, -ieta, iphi, ehdeptho[i].second);
              double ene = ehdeptho[i].first;
              if (ene > 0.0) {
                if (!(theHBHETopology_->validHcal(hcid0))) {
                  edm::LogWarning("HBHEMuon") << "(3) Invalid ID " << hcid0 << " with E = " << ene;
                  edm::LogWarning("HBHEMuon") << oppCell << " with " << ehdeptho.size() << " depths:";
                  for (const auto& ehd : ehdeptho)
                    edm::LogWarning("HBHEMuon") << " " << ehd.second << ":" << ehd.first;
                } else {
                  double chg(ene);
                  if (unCorrect_) {
                    double corr = (ignoreHECorr_ && (subdet0 == HcalEndcap)) ? 1.0 : respCorr(DetId(hcid0));
                    if (corr != 0)
                      ene /= corr;
#ifdef EDM_ML_DEBUG
                    HcalDetId id = (isItPlan1_ && isItPreRecHit_) ? hdc_->mergedDepthDetId(hcid0) : hcid0;
                    edm::LogVerbatim("HBHEMuon")
                        << hcid0 << ":" << id << " Corr " << corr << " E " << ene << ":" << ehdeptho[i].first;
#endif
                  }
                  if (getCharge_) {
                    double gain = gainFactor(conditions, hcid0);
                    if (gain != 0)
                      chg /= gain;
#ifdef EDM_ML_DEBUG
                    edm::LogVerbatim("HBHEMuon") << hcid0 << " Gain " << gain << " C " << chg;
#endif
                  }
                  int depth = ehdeptho[i].second - 1;
                  if (collapseDepth_) {
                    HcalDetId id = hdc_->mergedDepthDetId(hcid0);
                    depth = id.depth() - 1;
                  }
                  cHcalDepthHotBG[depth] += chg;
#ifdef EDM_ML_DEBUG
                  if ((verbosity_ % 10) > 0)
                    edm::LogVerbatim("HBHEMuon") << hcid0 << " Depth " << depth << " E " << ene << " C " << chg;
#endif
                }
              }
            }
          }
        }
#ifdef EDM_ML_DEBUG
        edm::LogVerbatim("HBHEMuon") << "Propagate Track to HCAL: " << trackID.okHCAL << " Match " << tmpmatch
                                     << " Hot " << isHot << " Energy " << eHcal << std::endl;
#endif

      } else {
        ecalDetId.push_back(0);
        hcalDetId.push_back(0);
        ehcalDetId.push_back(0);
      }

      matchedId.push_back(tmpmatch);
      ecal3x3Energy.push_back(eEcal);
      hcal1x1Energy.push_back(eHcal);
      hcal_ieta.push_back(ieta);
      hcal_iphi.push_back(iphi);
      for (int i = 0; i < depthMax_; ++i) {
        hcalDepthEnergy[i].push_back(eHcalDepth[i]);
        hcalDepthActiveLength[i].push_back(activeL[i]);
        hcalDepthEnergyHot[i].push_back(eHcalDepthHot[i]);
        hcalDepthActiveLengthHot[i].push_back(activeHotL[i]);
        hcalDepthEnergyCorr[i].push_back(eHcalDepthC[i]);
        hcalDepthEnergyHotCorr[i].push_back(eHcalDepthHotC[i]);
        hcalDepthChargeHot[i].push_back(cHcalDepthHot[i]);
        hcalDepthChargeHotBG[i].push_back(cHcalDepthHotBG[i]);
        hcalDepthMatch[i].push_back(matchDepth[i]);
        hcalDepthMatchHot[i].push_back(matchDepthHot[i]);
      }
      hcalActiveLength.push_back(activeLengthTot);
      hcalHot.push_back(isHot);
      hcalActiveLengthHot.push_back(activeLengthHotTot);
    }
  }
  if (accept) {
#ifdef EDM_ML_DEBUG
    for (unsigned int i = 0; i < hcal_ieta.size(); ++i)
      edm::LogVerbatim("HBHEMuon") << "[" << i << "] ieta/iphi for entry to "
                                   << "HCAL has value of " << hcal_ieta[i] << ":" << hcal_iphi[i];
#endif
    for (unsigned int k = 0; k < muon_is_good.size(); ++k) {
      muon_is_good_ = muon_is_good[k];
      muon_global_ = muon_global[k];
      muon_tracker_ = muon_tracker[k];
      muon_is_tight_ = muon_is_tight[k];
      muon_is_medium_ = muon_is_medium[k];
      ptGlob_ = ptGlob[k];
      etaGlob_ = etaGlob[k];
      phiGlob_ = phiGlob[k];
      energyMuon_ = energyMuon[k];
      pMuon_ = pMuon[k];
      muon_trkKink_ = muon_trkKink[k];
      muon_chi2LocalPosition_ = muon_chi2LocalPosition[k];
      muon_segComp_ = muon_segComp[k];
      trackerLayer_ = trackerLayer[k];
      numPixelLayers_ = numPixelLayers[k];
      tight_PixelHits_ = tight_PixelHits[k];
      innerTrack_ = innerTrack[k];
      outerTrack_ = outerTrack[k];
      globalTrack_ = globalTrack[k];
      chiTracker_ = chiTracker[k];
      dxyTracker_ = dxyTracker[k];
      dzTracker_ = dzTracker[k];
      innerTrackpt_ = innerTrackpt[k];
      innerTracketa_ = innerTracketa[k];
      innerTrackphi_ = innerTrackphi[k];
      tight_validFraction_ = tight_validFraction[k];
      outerTrackChi_ = outerTrackChi[k];
      outerTrackPt_ = outerTrackPt[k];
      outerTrackEta_ = outerTrackEta[k];
      outerTrackPhi_ = outerTrackPhi[k];
      outerTrackHits_ = outerTrackHits[k];
      outerTrackRHits_ = outerTrackRHits[k];
      globalTrckPt_ = globalTrckPt[k];
      globalTrckEta_ = globalTrckEta[k];
      globalTrckPhi_ = globalTrckPhi[k];
      globalMuonHits_ = globalMuonHits[k];
      matchedStat_ = matchedStat[k];
      chiGlobal_ = chiGlobal[k];
      tight_LongPara_ = tight_LongPara[k];
      tight_TransImpara_ = tight_TransImpara[k];
      isolationR04_ = isolationR04[k];
      isolationR03_ = isolationR03[k];
      ecalEnergy_ = ecalEnergy[k];
      hcalEnergy_ = hcalEnergy[k];
      hoEnergy_ = hoEnergy[k];
      matchedId_ = matchedId[k];
      hcalHot_ = hcalHot[k];
      ecal3x3Energy_ = ecal3x3Energy[k];
      hcal1x1Energy_ = hcal1x1Energy[k];
      ecalDetId_ = ecalDetId[k];
      hcalDetId_ = hcalDetId[k];
      ehcalDetId_ = ehcalDetId[k];
      hcal_ieta_ = hcal_ieta[k];
      hcal_iphi_ = hcal_iphi[k];
      for (int i = 0; i < depthMax_; ++i) {
        hcalDepthEnergy_[i] = hcalDepthEnergy[i][k];
        hcalDepthActiveLength_[i] = hcalDepthActiveLength[i][k];
        hcalDepthEnergyHot_[i] = hcalDepthEnergyHot[i][k];
        hcalDepthActiveLengthHot_[i] = hcalDepthActiveLengthHot[i][k];
        hcalDepthChargeHot_[i] = hcalDepthChargeHot[i][k];
        hcalDepthChargeHotBG_[i] = hcalDepthChargeHotBG[i][k];
        hcalDepthEnergyCorr_[i] = hcalDepthEnergyCorr[i][k];
        hcalDepthEnergyHotCorr_[i] = hcalDepthEnergyHotCorr[i][k];
        hcalDepthMatch_[i] = hcalDepthMatch[i][k];
        hcalDepthMatchHot_[i] = hcalDepthMatchHot[i][k];
      }
      hcalActiveLength_ = hcalActiveLength[k];
      hcalActiveLengthHot_ = hcalActiveLengthHot[k];
      tree_->Fill();
    }
  }
}

beginJob()

void HcalHBHEMuonAnalyzer::beginJob ( )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 915 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_.

                                    {
  edm::Service<TFileService> fs;
  tree_ = fs->make<TTree>("TREE", "TREE");
  tree_->Branch("Event_No", &eventNumber_);
  tree_->Branch("Run_No", &runNumber_);
  tree_->Branch("LumiNumber", &lumiNumber_);
  tree_->Branch("BXNumber", &bxNumber_);
  tree_->Branch("GoodVertex", &goodVertex_);
  tree_->Branch("PF_Muon", &muon_is_good_);
  tree_->Branch("Global_Muon", &muon_global_);
  tree_->Branch("Tracker_muon", &muon_tracker_);
  tree_->Branch("MuonIsTight", &muon_is_tight_);
  tree_->Branch("MuonIsMedium", &muon_is_medium_);
  tree_->Branch("pt_of_muon", &ptGlob_);
  tree_->Branch("eta_of_muon", &etaGlob_);
  tree_->Branch("phi_of_muon", &phiGlob_);
  tree_->Branch("energy_of_muon", &energyMuon_);
  tree_->Branch("p_of_muon", &pMuon_);
  tree_->Branch("muon_trkKink", &muon_trkKink_);
  tree_->Branch("muon_chi2LocalPosition", &muon_chi2LocalPosition_);
  tree_->Branch("muon_segComp", &muon_segComp_);

  tree_->Branch("TrackerLayer", &trackerLayer_);
  tree_->Branch("NumPixelLayers", &numPixelLayers_);
  tree_->Branch("InnerTrackPixelHits", &tight_PixelHits_);
  tree_->Branch("innerTrack", &innerTrack_);
  tree_->Branch("chiTracker", &chiTracker_);
  tree_->Branch("DxyTracker", &dxyTracker_);
  tree_->Branch("DzTracker", &dzTracker_);
  tree_->Branch("innerTrackpt", &innerTrackpt_);
  tree_->Branch("innerTracketa", &innerTracketa_);
  tree_->Branch("innerTrackphi", &innerTrackphi_);
  tree_->Branch("tight_validFraction", &tight_validFraction_);

  tree_->Branch("OuterTrack", &outerTrack_);
  tree_->Branch("OuterTrackChi", &outerTrackChi_);
  tree_->Branch("OuterTrackPt", &outerTrackPt_);
  tree_->Branch("OuterTrackEta", &outerTrackEta_);
  tree_->Branch("OuterTrackPhi", &outerTrackPhi_);
  tree_->Branch("OuterTrackHits", &outerTrackHits_);
  tree_->Branch("OuterTrackRHits", &outerTrackRHits_);

  tree_->Branch("GlobalTrack", &globalTrack_);
  tree_->Branch("GlobalTrckPt", &globalTrckPt_);
  tree_->Branch("GlobalTrckEta", &globalTrckEta_);
  tree_->Branch("GlobalTrckPhi", &globalTrckPhi_);
  tree_->Branch("Global_Muon_Hits", &globalMuonHits_);
  tree_->Branch("MatchedStations", &matchedStat_);
  tree_->Branch("GlobTrack_Chi", &chiGlobal_);
  tree_->Branch("Tight_LongitudinalImpactparameter", &tight_LongPara_);
  tree_->Branch("Tight_TransImpactparameter", &tight_TransImpara_);

  tree_->Branch("IsolationR04", &isolationR04_);
  tree_->Branch("IsolationR03", &isolationR03_);
  tree_->Branch("ecal_3into3", &ecalEnergy_);
  tree_->Branch("hcal_3into3", &hcalEnergy_);
  tree_->Branch("tracker_3into3", &hoEnergy_);

  tree_->Branch("matchedId", &matchedId_);
  tree_->Branch("hcal_cellHot", &hcalHot_);

  tree_->Branch("ecal_3x3", &ecal3x3Energy_);
  tree_->Branch("hcal_1x1", &hcal1x1Energy_);
  tree_->Branch("ecal_detID", &ecalDetId_);
  tree_->Branch("hcal_detID", &hcalDetId_);
  tree_->Branch("ehcal_detID", &ehcalDetId_);
  tree_->Branch("hcal_ieta", &hcal_ieta_);
  tree_->Branch("hcal_iphi", &hcal_iphi_);

  char name[100];
  for (int k = 0; k < maxDepth_; ++k) {
    sprintf(name, "hcal_edepth%d", (k + 1));
    tree_->Branch(name, &hcalDepthEnergy_[k]);
    sprintf(name, "hcal_activeL%d", (k + 1));
    tree_->Branch(name, &hcalDepthActiveLength_[k]);
    sprintf(name, "hcal_edepthHot%d", (k + 1));
    tree_->Branch(name, &hcalDepthEnergyHot_[k]);
    sprintf(name, "hcal_activeHotL%d", (k + 1));
    tree_->Branch(name, &hcalDepthActiveLengthHot_[k]);
    sprintf(name, "hcal_cdepthHot%d", (k + 1));
    tree_->Branch(name, &hcalDepthChargeHot_[k]);
    sprintf(name, "hcal_cdepthHotBG%d", (k + 1));
    tree_->Branch(name, &hcalDepthChargeHotBG_[k]);
    sprintf(name, "hcal_edepthCorrect%d", (k + 1));
    tree_->Branch(name, &hcalDepthEnergyCorr_[k]);
    sprintf(name, "hcal_edepthHotCorrect%d", (k + 1));
    tree_->Branch(name, &hcalDepthEnergyHotCorr_[k]);
    sprintf(name, "hcal_depthMatch%d", (k + 1));
    tree_->Branch(name, &hcalDepthMatch_[k]);
    sprintf(name, "hcal_depthMatchHot%d", (k + 1));
    tree_->Branch(name, &hcalDepthMatchHot_[k]);
  }

  tree_->Branch("activeLength", &hcalActiveLength_);
  tree_->Branch("activeLengthHot", &hcalActiveLengthHot_);

  tree_->Branch("hltresults", &hltresults_);
  tree_->Branch("all_triggers", &all_triggers_);
}

beginRun()

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

overrideprivate

Definition at line 1016 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_, verbosity_, and writeRespCorr_.

                                                                                   {
  hdc_ = &iSetup.getData(tok_ddrec_);
  actHB.clear();
  actHE.clear();
  actHB = hdc_->getThickActive(0);
  actHE = hdc_->getThickActive(1);
#ifdef EDM_ML_DEBUG
  unsigned int k1(0), k2(0);
  edm::LogVerbatim("HBHEMuon") << actHB.size() << " Active Length for HB";
  for (const auto& act : actHB) {
    edm::LogVerbatim("HBHEMuon") << "[" << k1 << "] ieta " << act.ieta << " depth " << act.depth << " zside "
                                 << act.zside << " type " << act.stype << " phi " << act.iphis.size() << ":"
                                 << act.iphis[0] << " L " << act.thick;
    HcalDetId hcid1(HcalBarrel, (act.ieta) * (act.zside), act.iphis[0], act.depth);
    HcalDetId hcid2 = mergedDepth_ ? hdc_->mergedDepthDetId(hcid1) : hcid1;
    edm::LogVerbatim("HBHEMuon") << hcid1 << " | " << hcid2 << " L "
                                 << activeLength(DetId(hcid2), ((verbosity_ / 100000) % 10));
    ++k1;
  }
  edm::LogVerbatim("HBHEMuon") << actHE.size() << " Active Length for HE";
  for (const auto& act : actHE) {
    edm::LogVerbatim("HBHEMuon") << "[" << k2 << "] ieta " << act.ieta << " depth " << act.depth << " zside "
                                 << act.zside << " type " << act.stype << " phi " << act.iphis.size() << ":"
                                 << act.iphis[0] << " L " << act.thick;
    HcalDetId hcid1(HcalEndcap, (act.ieta) * (act.zside), act.iphis[0], act.depth);
    HcalDetId hcid2 = mergedDepth_ ? hdc_->mergedDepthDetId(hcid1) : hcid1;
    edm::LogVerbatim("HBHEMuon") << hcid1 << " | " << hcid2 << " L "
                                 << activeLength(DetId(hcid2), ((verbosity_ / 100000) % 10));
    ++k2;
  }
#endif

  bool changed = true;
  all_triggers_.clear();
  if (hltConfig_.init(iRun, iSetup, "HLT", changed)) {
    // if init returns TRUE, initialisation has succeeded!
#ifdef EDM_ML_DEBUG
    edm::LogVerbatim("HBHEMuon") << "HLT config with process name HLT successfully extracted";
#endif
    unsigned int ntriggers = hltConfig_.size();
    for (unsigned int t = 0; t < ntriggers; ++t) {
      std::string hltname(hltConfig_.triggerName(t));
      for (unsigned int ik = 0; ik < triggers_.size(); ++ik) {
        if (hltname.find(triggers_[ik]) != std::string::npos) {
          all_triggers_.push_back(hltname);
          break;
        }
      }
    }  //loop over ntriggers
    edm::LogVerbatim("HBHEMuon") << "All triggers size in begin run " << all_triggers_.size();
  } else {
    edm::LogError("HBHEMuon") << "Error! HLT config extraction with process name HLT failed";
  }

  theHBHETopology_ = &iSetup.getData(tok_htopo_);
  const HcalRespCorrs* resp = &iSetup.getData(tok_respcorr_);
  respCorrs_ = new HcalRespCorrs(*resp);
  respCorrs_->setTopo(theHBHETopology_);
  geo_ = &iSetup.getData(tok_geom_);

  // Write correction factors for all HB/HE events
  if (writeRespCorr_) {
    const HcalGeometry* gHcal = static_cast<const HcalGeometry*>(geo_->getSubdetectorGeometry(DetId::Hcal, HcalBarrel));
    const std::vector<DetId>& ids = gHcal->getValidDetIds(DetId::Hcal, 0);
    edm::LogVerbatim("HBHEMuon") << "\nTable of Correction Factors for Run " << iRun.run() << "\n";
    for (auto const& id : ids) {
      if ((id.det() == DetId::Hcal) && ((id.subdetId() == HcalBarrel) || (id.subdetId() == HcalEndcap))) {
        edm::LogVerbatim("HBHEMuon") << HcalDetId(id) << " " << id.rawId() << " "
                                     << (respCorrs_->getValues(id))->getValue();
      }
    }
  }
}

clearVectors()

void HcalHBHEMuonAnalyzer::clearVectors ( )

private

clearing vectots

Definition at line 1119 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().

                                        {
  eventNumber_ = -99999;
  runNumber_ = -99999;
  lumiNumber_ = -99999;
  bxNumber_ = -99999;
  goodVertex_ = -99999;

  muon_is_good_ = false;
  muon_global_ = false;
  muon_tracker_ = false;
  ptGlob_ = 0;
  etaGlob_ = 0;
  phiGlob_ = 0;
  energyMuon_ = 0;
  pMuon_ = 0;
  muon_trkKink_ = 0;
  muon_chi2LocalPosition_ = 0;
  muon_segComp_ = 0;
  muon_is_tight_ = false;
  muon_is_medium_ = false;

  trackerLayer_ = 0;
  numPixelLayers_ = 0;
  tight_PixelHits_ = 0;
  innerTrack_ = false;
  chiTracker_ = 0;
  dxyTracker_ = 0;
  dzTracker_ = 0;
  innerTrackpt_ = 0;
  innerTracketa_ = 0;
  innerTrackphi_ = 0;
  tight_validFraction_ = 0;

  outerTrack_ = false;
  outerTrackPt_ = 0;
  outerTrackEta_ = 0;
  outerTrackPhi_ = 0;
  outerTrackHits_ = 0;
  outerTrackRHits_ = 0;
  outerTrackChi_ = 0;

  globalTrack_ = false;
  globalTrckPt_ = 0;
  globalTrckEta_ = 0;
  globalTrckPhi_ = 0;
  globalMuonHits_ = 0;
  matchedStat_ = 0;
  chiGlobal_ = 0;
  tight_LongPara_ = 0;
  tight_TransImpara_ = 0;

  isolationR04_ = 0;
  isolationR03_ = 0;
  ecalEnergy_ = 0;
  hcalEnergy_ = 0;
  hoEnergy_ = 0;
  matchedId_ = false;
  hcalHot_ = false;
  ecal3x3Energy_ = 0;
  hcal1x1Energy_ = 0;
  ecalDetId_ = 0;
  hcalDetId_ = 0;
  ehcalDetId_ = 0;
  hcal_ieta_ = 0;
  hcal_iphi_ = 0;
  for (int i = 0; i < maxDepth_; ++i) {
    hcalDepthEnergy_[i] = 0;
    hcalDepthActiveLength_[i] = 0;
    hcalDepthEnergyHot_[i] = 0;
    hcalDepthActiveLengthHot_[i] = 0;
    hcalDepthChargeHot_[i] = 0;
    hcalDepthChargeHotBG_[i] = 0;
    hcalDepthEnergyCorr_[i] = 0;
    hcalDepthEnergyHotCorr_[i] = 0;
    hcalDepthMatch_[i] = false;
    hcalDepthMatchHot_[i] = false;
  }
  hcalActiveLength_ = 0;
  hcalActiveLengthHot_ = 0;
  hltresults_.clear();
}

depth16HE()

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

private

Definition at line 1302 of file HcalHBHEMuonAnalyzer.cc.

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

Referenced by analyze().

                                                      {
  // Transition between HB/HE is special
  // For Run 1 or for Plan1 standard reconstruction it is 3
  // For runs beyond 2018 or in Plan1 for HEP17 it is 4
  int zside = (ieta > 0) ? 1 : -1;
  int depth = theHBHETopology_->dddConstants()->getMinDepth(1, 16, iphi, zside);
  if (isItPlan1_ && (!isItPreRecHit_))
    depth = 3;
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HBHEMuon") << "Plan1 " << isItPlan1_ << " PreRecHit " << isItPreRecHit_ << " phi " << iphi
                               << " depth " << depth;
#endif
  return depth;
}

endRun()

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

inlineoverrideprivate

Definition at line 84 of file HcalHBHEMuonAnalyzer.cc.

{}

fillDescriptions()

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

static

Definition at line 1091 of file HcalHBHEMuonAnalyzer.cc.

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

                                                                                      {
  edm::ParameterSetDescription desc;
  desc.add<edm::InputTag>("hlTriggerResults", edm::InputTag("TriggerResults", "", "HLT"));
  desc.add<edm::InputTag>("labelEBRecHit", edm::InputTag("ecalRecHit", "EcalRecHitsEB"));
  desc.add<edm::InputTag>("labelEERecHit", edm::InputTag("ecalRecHit", "EcalRecHitsEE"));
  desc.add<edm::InputTag>("labelHBHERecHit", edm::InputTag("hbhereco"));
  desc.add<std::string>("labelVertex", "offlinePrimaryVertices");
  desc.add<std::string>("labelMuon", "muons");
  std::vector<std::string> trig = {};
  desc.add<std::vector<std::string>>("triggers", trig);
  desc.add<double>("pMinMuon", 5.0);
  desc.addUntracked<int>("verbosity", 0);
  desc.add<int>("useRaw", 0);
  desc.add<bool>("unCorrect", true);
  desc.add<bool>("getCharge", true);
  desc.add<bool>("collapseDepth", false);
  desc.add<bool>("isItPlan1", false);
  desc.addUntracked<bool>("ignoreHECorr", false);
  desc.addUntracked<bool>("isItPreRecHit", false);
  desc.addUntracked<std::string>("moduleName", "");
  desc.addUntracked<std::string>("processName", "");
  desc.addUntracked<int>("maxDepth", 7);
  desc.addUntracked<std::string>("fileInCorr", "");
  desc.addUntracked<bool>("writeRespCorr", false);
  desc.add<bool>("usePFThreshold", true);
  descriptions.add("hcalHBHEMuon", desc);
}

gainFactor()

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

private

Definition at line 1294 of file HcalHBHEMuonAnalyzer.cc.

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

Referenced by analyze().

                                                                                            {
  double gain(0.0);
  const HcalCalibrations& calibs = conditions->getHcalCalibrations(id);
  for (int capid = 0; capid < 4; ++capid)
    gain += (0.25 * calibs.respcorrgain(capid));
  return gain;
}

goodCell()

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

private

Definition at line 1317 of file HcalHBHEMuonAnalyzer.cc.

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

Referenced by analyze().

                                                                 {
  std::pair<double, double> rz = hdc_->getRZ(hcid);
  bool typeRZ = (hcid.subdet() == HcalEndcap) ? false : true;
  bool match = spr::propagateHCAL(pTrack, geo, bField, typeRZ, rz, (((verbosity_ / 10000) % 10) > 0));
  return match;
}

isGoodVertex()

bool HcalHBHEMuonAnalyzer::isGoodVertex ( const reco::Vertex &  vtx )

private

Definition at line 1270 of file HcalHBHEMuonAnalyzer.cc.

References L1BJetProducer_cff::vtx.

Referenced by analyze().

                                                             {
  if (vtx.isFake())
    return false;
  if (vtx.ndof() < 4)
    return false;
  if (vtx.position().Rho() > 2.)
    return false;
  if (fabs(vtx.position().Z()) > 24.)
    return false;
  return true;
}

matchId()

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

private

Definition at line 1202 of file HcalHBHEMuonAnalyzer.cc.

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

Referenced by analyze().

                                                                            {
  HcalDetId kd1(id1.subdet(), id1.ieta(), id1.iphi(), 1);
  HcalDetId kd2(id1.subdet(), id2.ieta(), id2.iphi(), 1);
  int match = ((kd1 == kd2) ? 1 : 0);
  return match;
}

respCorr()

double HcalHBHEMuonAnalyzer::respCorr ( const DetId &  id )

private

Definition at line 1282 of file HcalHBHEMuonAnalyzer.cc.

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

Referenced by analyze().

                                                     {
  double cfac(1.0);
  if (useMyCorr_) {
    auto itr = corrValue_.find(id);
    if (itr != corrValue_.end())
      cfac = itr->second;
  } else if (respCorrs_ != nullptr) {
    cfac = (respCorrs_->getValues(id))->getValue();
  }
  return cfac;
}

Member Data Documentation

actHB

std::vector<HcalDDDRecConstants::HcalActiveLength> HcalHBHEMuonAnalyzer::actHB

private

Definition at line 174 of file HcalHBHEMuonAnalyzer.cc.

Referenced by activeLength(), and beginRun().

actHE

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

private

Definition at line 174 of file HcalHBHEMuonAnalyzer.cc.

Referenced by activeLength(), and beginRun().

all_triggers_

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

private

Definition at line 171 of file HcalHBHEMuonAnalyzer.cc.

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

bxNumber_

unsigned int HcalHBHEMuonAnalyzer::bxNumber_

private

Definition at line 138 of file HcalHBHEMuonAnalyzer.cc.

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

chiGlobal_

double HcalHBHEMuonAnalyzer::chiGlobal_

private

Definition at line 153 of file HcalHBHEMuonAnalyzer.cc.

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

chiTracker_

double HcalHBHEMuonAnalyzer::chiTracker_

private

Definition at line 146 of file HcalHBHEMuonAnalyzer.cc.

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

collapseDepth_

const bool HcalHBHEMuonAnalyzer::collapseDepth_

private

Definition at line 102 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze(), and HcalHBHEMuonAnalyzer().

corrValue_

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

private

Definition at line 175 of file HcalHBHEMuonAnalyzer.cc.

Referenced by HcalHBHEMuonAnalyzer(), and respCorr().

depthMax_

const int HcalHBHEMuonAnalyzer::depthMax_ = 7

staticprivate

Definition at line 136 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze().

dxyTracker_

double HcalHBHEMuonAnalyzer::dxyTracker_

private

Definition at line 146 of file HcalHBHEMuonAnalyzer.cc.

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

dzTracker_

double HcalHBHEMuonAnalyzer::dzTracker_

private

Definition at line 146 of file HcalHBHEMuonAnalyzer.cc.

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

ecal3x3Energy_

double HcalHBHEMuonAnalyzer::ecal3x3Energy_

private

Definition at line 157 of file HcalHBHEMuonAnalyzer.cc.

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

ecalDetId_

unsigned int HcalHBHEMuonAnalyzer::ecalDetId_

private

Definition at line 158 of file HcalHBHEMuonAnalyzer.cc.

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

ecalEnergy_

double HcalHBHEMuonAnalyzer::ecalEnergy_

private

Definition at line 155 of file HcalHBHEMuonAnalyzer.cc.

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

ehcalDetId_

unsigned int HcalHBHEMuonAnalyzer::ehcalDetId_

private

Definition at line 158 of file HcalHBHEMuonAnalyzer.cc.

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

energyMuon_

double HcalHBHEMuonAnalyzer::energyMuon_

private

Definition at line 142 of file HcalHBHEMuonAnalyzer.cc.

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

etaGlob_

double HcalHBHEMuonAnalyzer::etaGlob_

private

Definition at line 142 of file HcalHBHEMuonAnalyzer.cc.

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

eventNumber_

unsigned int HcalHBHEMuonAnalyzer::eventNumber_

private

Definition at line 138 of file HcalHBHEMuonAnalyzer.cc.

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

fileInCorr_

const std::string HcalHBHEMuonAnalyzer::fileInCorr_

private

Definition at line 98 of file HcalHBHEMuonAnalyzer.cc.

Referenced by HcalHBHEMuonAnalyzer().

geo_

const CaloGeometry* HcalHBHEMuonAnalyzer::geo_

private

Definition at line 129 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze(), and beginRun().

getCharge_

const bool HcalHBHEMuonAnalyzer::getCharge_

private

Definition at line 104 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze(), and HcalHBHEMuonAnalyzer().

globalMuonHits_

int HcalHBHEMuonAnalyzer::globalMuonHits_

private

Definition at line 152 of file HcalHBHEMuonAnalyzer.cc.

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

globalTrack_

bool HcalHBHEMuonAnalyzer::globalTrack_

private

Definition at line 145 of file HcalHBHEMuonAnalyzer.cc.

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

globalTrckEta_

double HcalHBHEMuonAnalyzer::globalTrckEta_

private

Definition at line 151 of file HcalHBHEMuonAnalyzer.cc.

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

globalTrckPhi_

double HcalHBHEMuonAnalyzer::globalTrckPhi_

private

Definition at line 151 of file HcalHBHEMuonAnalyzer.cc.

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

globalTrckPt_

double HcalHBHEMuonAnalyzer::globalTrckPt_

private

Definition at line 151 of file HcalHBHEMuonAnalyzer.cc.

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

goodVertex_

unsigned int HcalHBHEMuonAnalyzer::goodVertex_

private

Definition at line 139 of file HcalHBHEMuonAnalyzer.cc.

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

hcal1x1Energy_

double HcalHBHEMuonAnalyzer::hcal1x1Energy_

private

Definition at line 157 of file HcalHBHEMuonAnalyzer.cc.

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

hcal_ieta_

int HcalHBHEMuonAnalyzer::hcal_ieta_

private

Definition at line 159 of file HcalHBHEMuonAnalyzer.cc.

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

hcal_iphi_

int HcalHBHEMuonAnalyzer::hcal_iphi_

private

Definition at line 159 of file HcalHBHEMuonAnalyzer.cc.

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

hcalActiveLength_

double HcalHBHEMuonAnalyzer::hcalActiveLength_

private

Definition at line 170 of file HcalHBHEMuonAnalyzer.cc.

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

hcalActiveLengthHot_

double HcalHBHEMuonAnalyzer::hcalActiveLengthHot_

private

Definition at line 170 of file HcalHBHEMuonAnalyzer.cc.

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

hcalDepthActiveLength_

double HcalHBHEMuonAnalyzer::hcalDepthActiveLength_[depthMax_]

private

Definition at line 161 of file HcalHBHEMuonAnalyzer.cc.

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

hcalDepthActiveLengthHot_

double HcalHBHEMuonAnalyzer::hcalDepthActiveLengthHot_[depthMax_]

private

Definition at line 163 of file HcalHBHEMuonAnalyzer.cc.

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

hcalDepthChargeHot_

double HcalHBHEMuonAnalyzer::hcalDepthChargeHot_[depthMax_]

private

Definition at line 164 of file HcalHBHEMuonAnalyzer.cc.

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

hcalDepthChargeHotBG_

double HcalHBHEMuonAnalyzer::hcalDepthChargeHotBG_[depthMax_]

private

Definition at line 165 of file HcalHBHEMuonAnalyzer.cc.

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

hcalDepthEnergy_

double HcalHBHEMuonAnalyzer::hcalDepthEnergy_[depthMax_]

private

Definition at line 160 of file HcalHBHEMuonAnalyzer.cc.

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

hcalDepthEnergyCorr_

double HcalHBHEMuonAnalyzer::hcalDepthEnergyCorr_[depthMax_]

private

Definition at line 166 of file HcalHBHEMuonAnalyzer.cc.

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

hcalDepthEnergyHot_

double HcalHBHEMuonAnalyzer::hcalDepthEnergyHot_[depthMax_]

private

Definition at line 162 of file HcalHBHEMuonAnalyzer.cc.

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

hcalDepthEnergyHotCorr_

double HcalHBHEMuonAnalyzer::hcalDepthEnergyHotCorr_[depthMax_]

private

Definition at line 167 of file HcalHBHEMuonAnalyzer.cc.

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

hcalDepthMatch_

bool HcalHBHEMuonAnalyzer::hcalDepthMatch_[depthMax_]

private

Definition at line 168 of file HcalHBHEMuonAnalyzer.cc.

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

hcalDepthMatchHot_

bool HcalHBHEMuonAnalyzer::hcalDepthMatchHot_[depthMax_]

private

Definition at line 169 of file HcalHBHEMuonAnalyzer.cc.

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

hcalDetId_

unsigned int HcalHBHEMuonAnalyzer::hcalDetId_

private

Definition at line 158 of file HcalHBHEMuonAnalyzer.cc.

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

hcalEnergy_

double HcalHBHEMuonAnalyzer::hcalEnergy_

private

Definition at line 155 of file HcalHBHEMuonAnalyzer.cc.

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

hcalHot_

bool HcalHBHEMuonAnalyzer::hcalHot_

private

Definition at line 156 of file HcalHBHEMuonAnalyzer.cc.

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

hdc_

const HcalDDDRecConstants* HcalHBHEMuonAnalyzer::hdc_

private

Definition at line 127 of file HcalHBHEMuonAnalyzer.cc.

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

hltConfig_

HLTConfigProvider HcalHBHEMuonAnalyzer::hltConfig_

private

Definition at line 95 of file HcalHBHEMuonAnalyzer.cc.

Referenced by beginRun().

hltresults_

std::vector<int> HcalHBHEMuonAnalyzer::hltresults_

private

Definition at line 172 of file HcalHBHEMuonAnalyzer.cc.

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

hlTriggerResults_

const edm::InputTag HcalHBHEMuonAnalyzer::hlTriggerResults_

private

Definition at line 96 of file HcalHBHEMuonAnalyzer.cc.

Referenced by HcalHBHEMuonAnalyzer().

hoEnergy_

double HcalHBHEMuonAnalyzer::hoEnergy_

private

Definition at line 155 of file HcalHBHEMuonAnalyzer.cc.

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

ignoreHECorr_

const bool HcalHBHEMuonAnalyzer::ignoreHECorr_

private

Definition at line 103 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze(), and HcalHBHEMuonAnalyzer().

innerTrack_

bool HcalHBHEMuonAnalyzer::innerTrack_

private

Definition at line 145 of file HcalHBHEMuonAnalyzer.cc.

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

innerTracketa_

double HcalHBHEMuonAnalyzer::innerTracketa_

private

Definition at line 147 of file HcalHBHEMuonAnalyzer.cc.

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

innerTrackphi_

double HcalHBHEMuonAnalyzer::innerTrackphi_

private

Definition at line 147 of file HcalHBHEMuonAnalyzer.cc.

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

innerTrackpt_

double HcalHBHEMuonAnalyzer::innerTrackpt_

private

Definition at line 147 of file HcalHBHEMuonAnalyzer.cc.

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

isItPlan1_

const bool HcalHBHEMuonAnalyzer::isItPlan1_

private

Definition at line 102 of file HcalHBHEMuonAnalyzer.cc.

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

isItPreRecHit_

const bool HcalHBHEMuonAnalyzer::isItPreRecHit_

private

Definition at line 103 of file HcalHBHEMuonAnalyzer.cc.

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

isolationR03_

double HcalHBHEMuonAnalyzer::isolationR03_

private

Definition at line 154 of file HcalHBHEMuonAnalyzer.cc.

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

isolationR04_

double HcalHBHEMuonAnalyzer::isolationR04_

private

Definition at line 154 of file HcalHBHEMuonAnalyzer.cc.

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

kount_

int HcalHBHEMuonAnalyzer::kount_

private

Definition at line 133 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze(), and HcalHBHEMuonAnalyzer().

labelEBRecHit_

const edm::InputTag HcalHBHEMuonAnalyzer::labelEBRecHit_

private

Definition at line 97 of file HcalHBHEMuonAnalyzer.cc.

Referenced by HcalHBHEMuonAnalyzer().

labelEERecHit_

const edm::InputTag HcalHBHEMuonAnalyzer::labelEERecHit_

private

Definition at line 97 of file HcalHBHEMuonAnalyzer.cc.

Referenced by HcalHBHEMuonAnalyzer().

labelHBHERecHit_

const edm::InputTag HcalHBHEMuonAnalyzer::labelHBHERecHit_

private

Definition at line 97 of file HcalHBHEMuonAnalyzer.cc.

Referenced by HcalHBHEMuonAnalyzer().

labelMuon_

const std::string HcalHBHEMuonAnalyzer::labelMuon_

private

Definition at line 98 of file HcalHBHEMuonAnalyzer.cc.

Referenced by HcalHBHEMuonAnalyzer().

labelVtx_

const std::string HcalHBHEMuonAnalyzer::labelVtx_

private

Definition at line 98 of file HcalHBHEMuonAnalyzer.cc.

Referenced by HcalHBHEMuonAnalyzer().

lumiNumber_

unsigned int HcalHBHEMuonAnalyzer::lumiNumber_

private

Definition at line 138 of file HcalHBHEMuonAnalyzer.cc.

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

matchedId_

bool HcalHBHEMuonAnalyzer::matchedId_

private

Definition at line 156 of file HcalHBHEMuonAnalyzer.cc.

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

matchedStat_

int HcalHBHEMuonAnalyzer::matchedStat_

private

Definition at line 152 of file HcalHBHEMuonAnalyzer.cc.

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

maxDepth_

const int HcalHBHEMuonAnalyzer::maxDepth_

private

Definition at line 105 of file HcalHBHEMuonAnalyzer.cc.

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

mergedDepth_

bool HcalHBHEMuonAnalyzer::mergedDepth_

private

Definition at line 132 of file HcalHBHEMuonAnalyzer.cc.

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

modnam_

const std::string HcalHBHEMuonAnalyzer::modnam_

private

Definition at line 106 of file HcalHBHEMuonAnalyzer.cc.

Referenced by HcalHBHEMuonAnalyzer().

muon_chi2LocalPosition_

float HcalHBHEMuonAnalyzer::muon_chi2LocalPosition_

private

Definition at line 143 of file HcalHBHEMuonAnalyzer.cc.

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

muon_global_

bool HcalHBHEMuonAnalyzer::muon_global_

private

Definition at line 140 of file HcalHBHEMuonAnalyzer.cc.

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

muon_is_good_

bool HcalHBHEMuonAnalyzer::muon_is_good_

private

Definition at line 140 of file HcalHBHEMuonAnalyzer.cc.

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

muon_is_medium_

bool HcalHBHEMuonAnalyzer::muon_is_medium_

private

Definition at line 141 of file HcalHBHEMuonAnalyzer.cc.

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

muon_is_tight_

bool HcalHBHEMuonAnalyzer::muon_is_tight_

private

Definition at line 141 of file HcalHBHEMuonAnalyzer.cc.

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

muon_segComp_

float HcalHBHEMuonAnalyzer::muon_segComp_

private

Definition at line 143 of file HcalHBHEMuonAnalyzer.cc.

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

muon_tracker_

bool HcalHBHEMuonAnalyzer::muon_tracker_

private

Definition at line 140 of file HcalHBHEMuonAnalyzer.cc.

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

muon_trkKink_

float HcalHBHEMuonAnalyzer::muon_trkKink_

private

Definition at line 143 of file HcalHBHEMuonAnalyzer.cc.

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

numPixelLayers_

int HcalHBHEMuonAnalyzer::numPixelLayers_

private

Definition at line 144 of file HcalHBHEMuonAnalyzer.cc.

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

outerTrack_

bool HcalHBHEMuonAnalyzer::outerTrack_

private

Definition at line 145 of file HcalHBHEMuonAnalyzer.cc.

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

outerTrackChi_

double HcalHBHEMuonAnalyzer::outerTrackChi_

private

Definition at line 148 of file HcalHBHEMuonAnalyzer.cc.

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

outerTrackEta_

double HcalHBHEMuonAnalyzer::outerTrackEta_

private

Definition at line 149 of file HcalHBHEMuonAnalyzer.cc.

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

outerTrackHits_

int HcalHBHEMuonAnalyzer::outerTrackHits_

private

Definition at line 150 of file HcalHBHEMuonAnalyzer.cc.

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

outerTrackPhi_

double HcalHBHEMuonAnalyzer::outerTrackPhi_

private

Definition at line 149 of file HcalHBHEMuonAnalyzer.cc.

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

outerTrackPt_

double HcalHBHEMuonAnalyzer::outerTrackPt_

private

Definition at line 149 of file HcalHBHEMuonAnalyzer.cc.

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

outerTrackRHits_

int HcalHBHEMuonAnalyzer::outerTrackRHits_

private

Definition at line 150 of file HcalHBHEMuonAnalyzer.cc.

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

phiGlob_

double HcalHBHEMuonAnalyzer::phiGlob_

private

Definition at line 142 of file HcalHBHEMuonAnalyzer.cc.

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

pMinMuon_

const double HcalHBHEMuonAnalyzer::pMinMuon_

private

Definition at line 100 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze(), and HcalHBHEMuonAnalyzer().

pMuon_

double HcalHBHEMuonAnalyzer::pMuon_

private

Definition at line 142 of file HcalHBHEMuonAnalyzer.cc.

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

procnm_

const std::string HcalHBHEMuonAnalyzer::procnm_

private

Definition at line 106 of file HcalHBHEMuonAnalyzer.cc.

Referenced by HcalHBHEMuonAnalyzer().

ptGlob_

double HcalHBHEMuonAnalyzer::ptGlob_

private

Definition at line 142 of file HcalHBHEMuonAnalyzer.cc.

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

respCorrs_

HcalRespCorrs* HcalHBHEMuonAnalyzer::respCorrs_

private

Definition at line 130 of file HcalHBHEMuonAnalyzer.cc.

Referenced by beginRun(), and respCorr().

runNumber_

unsigned int HcalHBHEMuonAnalyzer::runNumber_

private

Definition at line 138 of file HcalHBHEMuonAnalyzer.cc.

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

theHBHETopology_

const HcalTopology* HcalHBHEMuonAnalyzer::theHBHETopology_

private

Definition at line 128 of file HcalHBHEMuonAnalyzer.cc.

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

tight_LongPara_

double HcalHBHEMuonAnalyzer::tight_LongPara_

private

Definition at line 153 of file HcalHBHEMuonAnalyzer.cc.

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

tight_PixelHits_

int HcalHBHEMuonAnalyzer::tight_PixelHits_

private

Definition at line 144 of file HcalHBHEMuonAnalyzer.cc.

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

tight_TransImpara_

double HcalHBHEMuonAnalyzer::tight_TransImpara_

private

Definition at line 153 of file HcalHBHEMuonAnalyzer.cc.

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

tight_validFraction_

double HcalHBHEMuonAnalyzer::tight_validFraction_

private

Definition at line 148 of file HcalHBHEMuonAnalyzer.cc.

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

tok_chan_

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

private

Definition at line 121 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze().

tok_dbservice_

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

private

Definition at line 124 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze().

tok_ddrec_

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

private

Definition at line 116 of file HcalHBHEMuonAnalyzer.cc.

Referenced by beginRun().

tok_EB_

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

private

Definition at line 111 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze().

tok_ecalPFRecHitThresholds_

const edm::ESGetToken<EcalPFRecHitThresholds, EcalPFRecHitThresholdsRcd> HcalHBHEMuonAnalyzer::tok_ecalPFRecHitThresholds_

private

Definition at line 125 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze().

tok_EE_

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

private

Definition at line 112 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze().

tok_geom_

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

private

Definition at line 119 of file HcalHBHEMuonAnalyzer.cc.

Referenced by beginRun().

tok_HBHE_

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

private

Definition at line 113 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze().

tok_htopo_

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

private

Definition at line 117 of file HcalHBHEMuonAnalyzer.cc.

Referenced by beginRun().

tok_magField_

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

private

Definition at line 120 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze().

tok_Muon_

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

private

Definition at line 114 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze().

tok_respcorr_

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

private

Definition at line 118 of file HcalHBHEMuonAnalyzer.cc.

Referenced by beginRun().

tok_sevlv_

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

private

Definition at line 122 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze().

tok_topo_

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

private

Definition at line 123 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze().

tok_trigRes_

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

private

Definition at line 109 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze().

tok_Vtx_

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

private

Definition at line 110 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze().

trackerLayer_

int HcalHBHEMuonAnalyzer::trackerLayer_

private

Definition at line 144 of file HcalHBHEMuonAnalyzer.cc.

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

tree_

TTree* HcalHBHEMuonAnalyzer::tree_

private

Definition at line 137 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze(), and beginJob().

triggers_

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

private

Definition at line 99 of file HcalHBHEMuonAnalyzer.cc.

Referenced by beginRun().

unCorrect_

const bool HcalHBHEMuonAnalyzer::unCorrect_

private

Definition at line 102 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze(), and HcalHBHEMuonAnalyzer().

useMyCorr_

bool HcalHBHEMuonAnalyzer::useMyCorr_

private

Definition at line 132 of file HcalHBHEMuonAnalyzer.cc.

Referenced by HcalHBHEMuonAnalyzer(), and respCorr().

usePFThresh_

const bool HcalHBHEMuonAnalyzer::usePFThresh_

private

Definition at line 107 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze(), and HcalHBHEMuonAnalyzer().

useRaw_

const int HcalHBHEMuonAnalyzer::useRaw_

private

Definition at line 101 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze(), and HcalHBHEMuonAnalyzer().

verbosity_

const int HcalHBHEMuonAnalyzer::verbosity_

private

Definition at line 101 of file HcalHBHEMuonAnalyzer.cc.

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

writeRespCorr_

const bool HcalHBHEMuonAnalyzer::writeRespCorr_

private

Definition at line 104 of file HcalHBHEMuonAnalyzer.cc.

Referenced by beginRun().