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 noexcept final

bool	wantsGlobalRuns () const noexcept final

bool	wantsInputProcessBlocks () const noexcept final

bool	wantsProcessBlocks () const noexcept 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 noexcept

bool	wantsStreamRuns () const noexcept

	~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 70 of file HcalHBHEMuonAnalyzer.cc.

Constructor & Destructor Documentation

◆ HcalHBHEMuonAnalyzer()

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

explicit

Definition at line 175 of file HcalHBHEMuonAnalyzer.cc.

References collapseDepth_, corrValue_, fileInCorr_, getCharge_, hlTriggerResults_, l1ctLayer2EG_cff::id, ignoreHECorr_, timingPdfMaker::infile, 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 1205 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 262 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_, electrons_cff::ecalEnergy, 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(), ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, 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 911 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 1012 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 1115 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 1298 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 80 of file HcalHBHEMuonAnalyzer.cc.

80 {}

◆ fillDescriptions()

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

static

Definition at line 1087 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 1290 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 1313 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 1266 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 1198 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 1278 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 170 of file HcalHBHEMuonAnalyzer.cc.

Referenced by activeLength(), and beginRun().

◆ actHE

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

private

Definition at line 170 of file HcalHBHEMuonAnalyzer.cc.

Referenced by activeLength(), and beginRun().

◆ all_triggers_

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

private

Definition at line 167 of file HcalHBHEMuonAnalyzer.cc.

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

◆ bxNumber_

unsigned int HcalHBHEMuonAnalyzer::bxNumber_

private

Definition at line 134 of file HcalHBHEMuonAnalyzer.cc.

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

◆ chiGlobal_

double HcalHBHEMuonAnalyzer::chiGlobal_

private

Definition at line 149 of file HcalHBHEMuonAnalyzer.cc.

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

◆ chiTracker_

double HcalHBHEMuonAnalyzer::chiTracker_

private

Definition at line 142 of file HcalHBHEMuonAnalyzer.cc.

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

◆ collapseDepth_

const bool HcalHBHEMuonAnalyzer::collapseDepth_

private

Definition at line 98 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze(), and HcalHBHEMuonAnalyzer().

◆ corrValue_

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

private

Definition at line 171 of file HcalHBHEMuonAnalyzer.cc.

Referenced by HcalHBHEMuonAnalyzer(), and respCorr().

◆ depthMax_

const int HcalHBHEMuonAnalyzer::depthMax_ = 7

staticprivate

Definition at line 132 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze().

◆ dxyTracker_

double HcalHBHEMuonAnalyzer::dxyTracker_

private

Definition at line 142 of file HcalHBHEMuonAnalyzer.cc.

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

◆ dzTracker_

double HcalHBHEMuonAnalyzer::dzTracker_

private

Definition at line 142 of file HcalHBHEMuonAnalyzer.cc.

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

◆ ecal3x3Energy_

double HcalHBHEMuonAnalyzer::ecal3x3Energy_

private

Definition at line 153 of file HcalHBHEMuonAnalyzer.cc.

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

◆ ecalDetId_

unsigned int HcalHBHEMuonAnalyzer::ecalDetId_

private

Definition at line 154 of file HcalHBHEMuonAnalyzer.cc.

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

◆ ecalEnergy_

double HcalHBHEMuonAnalyzer::ecalEnergy_

private

Definition at line 151 of file HcalHBHEMuonAnalyzer.cc.

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

◆ ehcalDetId_

unsigned int HcalHBHEMuonAnalyzer::ehcalDetId_

private

Definition at line 154 of file HcalHBHEMuonAnalyzer.cc.

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

◆ energyMuon_

double HcalHBHEMuonAnalyzer::energyMuon_

private

Definition at line 138 of file HcalHBHEMuonAnalyzer.cc.

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

◆ etaGlob_

double HcalHBHEMuonAnalyzer::etaGlob_

private

Definition at line 138 of file HcalHBHEMuonAnalyzer.cc.

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

◆ eventNumber_

unsigned int HcalHBHEMuonAnalyzer::eventNumber_

private

Definition at line 134 of file HcalHBHEMuonAnalyzer.cc.

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

◆ fileInCorr_

const std::string HcalHBHEMuonAnalyzer::fileInCorr_

private

Definition at line 94 of file HcalHBHEMuonAnalyzer.cc.

Referenced by HcalHBHEMuonAnalyzer().

◆ geo_

const CaloGeometry* HcalHBHEMuonAnalyzer::geo_

private

Definition at line 125 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze(), and beginRun().

◆ getCharge_

const bool HcalHBHEMuonAnalyzer::getCharge_

private

Definition at line 100 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze(), and HcalHBHEMuonAnalyzer().

◆ globalMuonHits_

int HcalHBHEMuonAnalyzer::globalMuonHits_

private

Definition at line 148 of file HcalHBHEMuonAnalyzer.cc.

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

◆ globalTrack_

bool HcalHBHEMuonAnalyzer::globalTrack_

private

Definition at line 141 of file HcalHBHEMuonAnalyzer.cc.

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

◆ globalTrckEta_

double HcalHBHEMuonAnalyzer::globalTrckEta_

private

Definition at line 147 of file HcalHBHEMuonAnalyzer.cc.

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

◆ globalTrckPhi_

double HcalHBHEMuonAnalyzer::globalTrckPhi_

private

Definition at line 147 of file HcalHBHEMuonAnalyzer.cc.

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

◆ globalTrckPt_

double HcalHBHEMuonAnalyzer::globalTrckPt_

private

Definition at line 147 of file HcalHBHEMuonAnalyzer.cc.

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

◆ goodVertex_

unsigned int HcalHBHEMuonAnalyzer::goodVertex_

private

Definition at line 135 of file HcalHBHEMuonAnalyzer.cc.

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

◆ hcal1x1Energy_

double HcalHBHEMuonAnalyzer::hcal1x1Energy_

private

Definition at line 153 of file HcalHBHEMuonAnalyzer.cc.

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

◆ hcal_ieta_

int HcalHBHEMuonAnalyzer::hcal_ieta_

private

Definition at line 155 of file HcalHBHEMuonAnalyzer.cc.

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

◆ hcal_iphi_

int HcalHBHEMuonAnalyzer::hcal_iphi_

private

Definition at line 155 of file HcalHBHEMuonAnalyzer.cc.

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

◆ hcalActiveLength_

double HcalHBHEMuonAnalyzer::hcalActiveLength_

private

Definition at line 166 of file HcalHBHEMuonAnalyzer.cc.

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

◆ hcalActiveLengthHot_

double HcalHBHEMuonAnalyzer::hcalActiveLengthHot_

private

Definition at line 166 of file HcalHBHEMuonAnalyzer.cc.

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

◆ hcalDepthActiveLength_

double HcalHBHEMuonAnalyzer::hcalDepthActiveLength_[depthMax_]

private

Definition at line 157 of file HcalHBHEMuonAnalyzer.cc.

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

◆ hcalDepthActiveLengthHot_

double HcalHBHEMuonAnalyzer::hcalDepthActiveLengthHot_[depthMax_]

private

Definition at line 159 of file HcalHBHEMuonAnalyzer.cc.

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

◆ hcalDepthChargeHot_

double HcalHBHEMuonAnalyzer::hcalDepthChargeHot_[depthMax_]

private

Definition at line 160 of file HcalHBHEMuonAnalyzer.cc.

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

◆ hcalDepthChargeHotBG_

double HcalHBHEMuonAnalyzer::hcalDepthChargeHotBG_[depthMax_]

private

Definition at line 161 of file HcalHBHEMuonAnalyzer.cc.

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

◆ hcalDepthEnergy_

double HcalHBHEMuonAnalyzer::hcalDepthEnergy_[depthMax_]

private

Definition at line 156 of file HcalHBHEMuonAnalyzer.cc.

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

◆ hcalDepthEnergyCorr_

double HcalHBHEMuonAnalyzer::hcalDepthEnergyCorr_[depthMax_]

private

Definition at line 162 of file HcalHBHEMuonAnalyzer.cc.

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

◆ hcalDepthEnergyHot_

double HcalHBHEMuonAnalyzer::hcalDepthEnergyHot_[depthMax_]

private

Definition at line 158 of file HcalHBHEMuonAnalyzer.cc.

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

◆ hcalDepthEnergyHotCorr_

double HcalHBHEMuonAnalyzer::hcalDepthEnergyHotCorr_[depthMax_]

private

Definition at line 163 of file HcalHBHEMuonAnalyzer.cc.

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

◆ hcalDepthMatch_

bool HcalHBHEMuonAnalyzer::hcalDepthMatch_[depthMax_]

private

Definition at line 164 of file HcalHBHEMuonAnalyzer.cc.

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

◆ hcalDepthMatchHot_

bool HcalHBHEMuonAnalyzer::hcalDepthMatchHot_[depthMax_]

private

Definition at line 165 of file HcalHBHEMuonAnalyzer.cc.

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

◆ hcalDetId_

unsigned int HcalHBHEMuonAnalyzer::hcalDetId_

private

Definition at line 154 of file HcalHBHEMuonAnalyzer.cc.

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

◆ hcalEnergy_

double HcalHBHEMuonAnalyzer::hcalEnergy_

private

Definition at line 151 of file HcalHBHEMuonAnalyzer.cc.

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

◆ hcalHot_

bool HcalHBHEMuonAnalyzer::hcalHot_

private

Definition at line 152 of file HcalHBHEMuonAnalyzer.cc.

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

◆ hdc_

const HcalDDDRecConstants* HcalHBHEMuonAnalyzer::hdc_

private

Definition at line 123 of file HcalHBHEMuonAnalyzer.cc.

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

◆ hltConfig_

HLTConfigProvider HcalHBHEMuonAnalyzer::hltConfig_

private

Definition at line 91 of file HcalHBHEMuonAnalyzer.cc.

Referenced by beginRun().

◆ hltresults_

std::vector<int> HcalHBHEMuonAnalyzer::hltresults_

private

Definition at line 168 of file HcalHBHEMuonAnalyzer.cc.

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

◆ hlTriggerResults_

const edm::InputTag HcalHBHEMuonAnalyzer::hlTriggerResults_

private

Definition at line 92 of file HcalHBHEMuonAnalyzer.cc.

Referenced by HcalHBHEMuonAnalyzer().

◆ hoEnergy_

double HcalHBHEMuonAnalyzer::hoEnergy_

private

Definition at line 151 of file HcalHBHEMuonAnalyzer.cc.

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

◆ ignoreHECorr_

const bool HcalHBHEMuonAnalyzer::ignoreHECorr_

private

Definition at line 99 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze(), and HcalHBHEMuonAnalyzer().

◆ innerTrack_

bool HcalHBHEMuonAnalyzer::innerTrack_

private

Definition at line 141 of file HcalHBHEMuonAnalyzer.cc.

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

◆ innerTracketa_

double HcalHBHEMuonAnalyzer::innerTracketa_

private

Definition at line 143 of file HcalHBHEMuonAnalyzer.cc.

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

◆ innerTrackphi_

double HcalHBHEMuonAnalyzer::innerTrackphi_

private

Definition at line 143 of file HcalHBHEMuonAnalyzer.cc.

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

◆ innerTrackpt_

double HcalHBHEMuonAnalyzer::innerTrackpt_

private

Definition at line 143 of file HcalHBHEMuonAnalyzer.cc.

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

◆ isItPlan1_

const bool HcalHBHEMuonAnalyzer::isItPlan1_

private

Definition at line 98 of file HcalHBHEMuonAnalyzer.cc.

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

◆ isItPreRecHit_

const bool HcalHBHEMuonAnalyzer::isItPreRecHit_

private

Definition at line 99 of file HcalHBHEMuonAnalyzer.cc.

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

◆ isolationR03_

double HcalHBHEMuonAnalyzer::isolationR03_

private

Definition at line 150 of file HcalHBHEMuonAnalyzer.cc.

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

◆ isolationR04_

double HcalHBHEMuonAnalyzer::isolationR04_

private

Definition at line 150 of file HcalHBHEMuonAnalyzer.cc.

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

◆ kount_

int HcalHBHEMuonAnalyzer::kount_

private

Definition at line 129 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze(), and HcalHBHEMuonAnalyzer().

◆ labelEBRecHit_

const edm::InputTag HcalHBHEMuonAnalyzer::labelEBRecHit_

private

Definition at line 93 of file HcalHBHEMuonAnalyzer.cc.

Referenced by HcalHBHEMuonAnalyzer().

◆ labelEERecHit_

const edm::InputTag HcalHBHEMuonAnalyzer::labelEERecHit_

private

Definition at line 93 of file HcalHBHEMuonAnalyzer.cc.

Referenced by HcalHBHEMuonAnalyzer().

◆ labelHBHERecHit_

const edm::InputTag HcalHBHEMuonAnalyzer::labelHBHERecHit_

private

Definition at line 93 of file HcalHBHEMuonAnalyzer.cc.

Referenced by HcalHBHEMuonAnalyzer().

◆ labelMuon_

const std::string HcalHBHEMuonAnalyzer::labelMuon_

private

Definition at line 94 of file HcalHBHEMuonAnalyzer.cc.

Referenced by HcalHBHEMuonAnalyzer().

◆ labelVtx_

const std::string HcalHBHEMuonAnalyzer::labelVtx_

private

Definition at line 94 of file HcalHBHEMuonAnalyzer.cc.

Referenced by HcalHBHEMuonAnalyzer().

◆ lumiNumber_

unsigned int HcalHBHEMuonAnalyzer::lumiNumber_

private

Definition at line 134 of file HcalHBHEMuonAnalyzer.cc.

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

◆ matchedId_

bool HcalHBHEMuonAnalyzer::matchedId_

private

Definition at line 152 of file HcalHBHEMuonAnalyzer.cc.

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

◆ matchedStat_

int HcalHBHEMuonAnalyzer::matchedStat_

private

Definition at line 148 of file HcalHBHEMuonAnalyzer.cc.

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

◆ maxDepth_

const int HcalHBHEMuonAnalyzer::maxDepth_

private

Definition at line 101 of file HcalHBHEMuonAnalyzer.cc.

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

◆ mergedDepth_

bool HcalHBHEMuonAnalyzer::mergedDepth_

private

Definition at line 128 of file HcalHBHEMuonAnalyzer.cc.

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

◆ modnam_

const std::string HcalHBHEMuonAnalyzer::modnam_

private

Definition at line 102 of file HcalHBHEMuonAnalyzer.cc.

Referenced by HcalHBHEMuonAnalyzer().

◆ muon_chi2LocalPosition_

float HcalHBHEMuonAnalyzer::muon_chi2LocalPosition_

private

Definition at line 139 of file HcalHBHEMuonAnalyzer.cc.

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

◆ muon_global_

bool HcalHBHEMuonAnalyzer::muon_global_

private

Definition at line 136 of file HcalHBHEMuonAnalyzer.cc.

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

◆ muon_is_good_

bool HcalHBHEMuonAnalyzer::muon_is_good_

private

Definition at line 136 of file HcalHBHEMuonAnalyzer.cc.

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

◆ muon_is_medium_

bool HcalHBHEMuonAnalyzer::muon_is_medium_

private

Definition at line 137 of file HcalHBHEMuonAnalyzer.cc.

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

◆ muon_is_tight_

bool HcalHBHEMuonAnalyzer::muon_is_tight_

private

Definition at line 137 of file HcalHBHEMuonAnalyzer.cc.

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

◆ muon_segComp_

float HcalHBHEMuonAnalyzer::muon_segComp_

private

Definition at line 139 of file HcalHBHEMuonAnalyzer.cc.

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

◆ muon_tracker_

bool HcalHBHEMuonAnalyzer::muon_tracker_

private

Definition at line 136 of file HcalHBHEMuonAnalyzer.cc.

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

◆ muon_trkKink_

float HcalHBHEMuonAnalyzer::muon_trkKink_

private

Definition at line 139 of file HcalHBHEMuonAnalyzer.cc.

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

◆ numPixelLayers_

int HcalHBHEMuonAnalyzer::numPixelLayers_

private

Definition at line 140 of file HcalHBHEMuonAnalyzer.cc.

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

◆ outerTrack_

bool HcalHBHEMuonAnalyzer::outerTrack_

private

Definition at line 141 of file HcalHBHEMuonAnalyzer.cc.

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

◆ outerTrackChi_

double HcalHBHEMuonAnalyzer::outerTrackChi_

private

Definition at line 144 of file HcalHBHEMuonAnalyzer.cc.

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

◆ outerTrackEta_

double HcalHBHEMuonAnalyzer::outerTrackEta_

private

Definition at line 145 of file HcalHBHEMuonAnalyzer.cc.

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

◆ outerTrackHits_

int HcalHBHEMuonAnalyzer::outerTrackHits_

private

Definition at line 146 of file HcalHBHEMuonAnalyzer.cc.

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

◆ outerTrackPhi_

double HcalHBHEMuonAnalyzer::outerTrackPhi_

private

Definition at line 145 of file HcalHBHEMuonAnalyzer.cc.

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

◆ outerTrackPt_

double HcalHBHEMuonAnalyzer::outerTrackPt_

private

Definition at line 145 of file HcalHBHEMuonAnalyzer.cc.

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

◆ outerTrackRHits_

int HcalHBHEMuonAnalyzer::outerTrackRHits_

private

Definition at line 146 of file HcalHBHEMuonAnalyzer.cc.

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

◆ phiGlob_

double HcalHBHEMuonAnalyzer::phiGlob_

private

Definition at line 138 of file HcalHBHEMuonAnalyzer.cc.

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

◆ pMinMuon_

const double HcalHBHEMuonAnalyzer::pMinMuon_

private

Definition at line 96 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze(), and HcalHBHEMuonAnalyzer().

◆ pMuon_

double HcalHBHEMuonAnalyzer::pMuon_

private

Definition at line 138 of file HcalHBHEMuonAnalyzer.cc.

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

◆ procnm_

const std::string HcalHBHEMuonAnalyzer::procnm_

private

Definition at line 102 of file HcalHBHEMuonAnalyzer.cc.

Referenced by HcalHBHEMuonAnalyzer().

◆ ptGlob_

double HcalHBHEMuonAnalyzer::ptGlob_

private

Definition at line 138 of file HcalHBHEMuonAnalyzer.cc.

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

◆ respCorrs_

HcalRespCorrs* HcalHBHEMuonAnalyzer::respCorrs_

private

Definition at line 126 of file HcalHBHEMuonAnalyzer.cc.

Referenced by beginRun(), and respCorr().

◆ runNumber_

unsigned int HcalHBHEMuonAnalyzer::runNumber_

private

Definition at line 134 of file HcalHBHEMuonAnalyzer.cc.

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

◆ theHBHETopology_

const HcalTopology* HcalHBHEMuonAnalyzer::theHBHETopology_

private

Definition at line 124 of file HcalHBHEMuonAnalyzer.cc.

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

◆ tight_LongPara_

double HcalHBHEMuonAnalyzer::tight_LongPara_

private

Definition at line 149 of file HcalHBHEMuonAnalyzer.cc.

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

◆ tight_PixelHits_

int HcalHBHEMuonAnalyzer::tight_PixelHits_

private

Definition at line 140 of file HcalHBHEMuonAnalyzer.cc.

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

◆ tight_TransImpara_

double HcalHBHEMuonAnalyzer::tight_TransImpara_

private

Definition at line 149 of file HcalHBHEMuonAnalyzer.cc.

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

◆ tight_validFraction_

double HcalHBHEMuonAnalyzer::tight_validFraction_

private

Definition at line 144 of file HcalHBHEMuonAnalyzer.cc.

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

◆ tok_chan_

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

private

Definition at line 117 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze().

◆ tok_dbservice_

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

private

Definition at line 120 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze().

◆ tok_ddrec_

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

private

Definition at line 112 of file HcalHBHEMuonAnalyzer.cc.

Referenced by beginRun().

◆ tok_EB_

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

private

Definition at line 107 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze().

◆ tok_ecalPFRecHitThresholds_

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

private

Definition at line 121 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze().

◆ tok_EE_

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

private

Definition at line 108 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze().

◆ tok_geom_

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

private

Definition at line 115 of file HcalHBHEMuonAnalyzer.cc.

Referenced by beginRun().

◆ tok_HBHE_

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

private

Definition at line 109 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze().

◆ tok_htopo_

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

private

Definition at line 113 of file HcalHBHEMuonAnalyzer.cc.

Referenced by beginRun().

◆ tok_magField_

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

private

Definition at line 116 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze().

◆ tok_Muon_

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

private

Definition at line 110 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze().

◆ tok_respcorr_

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

private

Definition at line 114 of file HcalHBHEMuonAnalyzer.cc.

Referenced by beginRun().

◆ tok_sevlv_

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

private

Definition at line 118 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze().

◆ tok_topo_

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

private

Definition at line 119 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze().

◆ tok_trigRes_

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

private

Definition at line 105 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze().

◆ tok_Vtx_

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

private

Definition at line 106 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze().

◆ trackerLayer_

int HcalHBHEMuonAnalyzer::trackerLayer_

private

Definition at line 140 of file HcalHBHEMuonAnalyzer.cc.

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

◆ tree_

TTree* HcalHBHEMuonAnalyzer::tree_

private

Definition at line 133 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ triggers_

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

private

Definition at line 95 of file HcalHBHEMuonAnalyzer.cc.

Referenced by beginRun().

◆ unCorrect_

const bool HcalHBHEMuonAnalyzer::unCorrect_

private

Definition at line 98 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze(), and HcalHBHEMuonAnalyzer().

◆ useMyCorr_

bool HcalHBHEMuonAnalyzer::useMyCorr_

private

Definition at line 128 of file HcalHBHEMuonAnalyzer.cc.

Referenced by HcalHBHEMuonAnalyzer(), and respCorr().

◆ usePFThresh_

const bool HcalHBHEMuonAnalyzer::usePFThresh_

private

Definition at line 103 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze(), and HcalHBHEMuonAnalyzer().

◆ useRaw_

const int HcalHBHEMuonAnalyzer::useRaw_

private

Definition at line 97 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze(), and HcalHBHEMuonAnalyzer().

◆ verbosity_

const int HcalHBHEMuonAnalyzer::verbosity_

private

Definition at line 97 of file HcalHBHEMuonAnalyzer.cc.

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

◆ writeRespCorr_

const bool HcalHBHEMuonAnalyzer::writeRespCorr_

private

Definition at line 100 of file HcalHBHEMuonAnalyzer.cc.

Referenced by beginRun().

Public Member Functions

Static Public Member Functions

Private Member Functions

Private Attributes

Static Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ HcalHBHEMuonAnalyzer()

Member Function Documentation

◆ activeLength()

◆ analyze()

◆ beginJob()

◆ beginRun()

◆ clearVectors()

◆ depth16HE()

◆ endRun()

◆ fillDescriptions()

◆ gainFactor()

◆ goodCell()

◆ isGoodVertex()

◆ matchId()

◆ respCorr()

Member Data Documentation

◆ actHB

◆ actHE

◆ all_triggers_

◆ bxNumber_

◆ chiGlobal_

◆ chiTracker_

◆ collapseDepth_

◆ corrValue_

◆ depthMax_

◆ dxyTracker_

◆ dzTracker_

◆ ecal3x3Energy_

◆ ecalDetId_

◆ ecalEnergy_

◆ ehcalDetId_

◆ energyMuon_

◆ etaGlob_

◆ eventNumber_

◆ fileInCorr_

◆ geo_

◆ getCharge_

◆ 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_

◆ hdc_

◆ hltConfig_

◆ hltresults_

◆ hlTriggerResults_

◆ hoEnergy_

◆ ignoreHECorr_

◆ innerTrack_

◆ innerTracketa_

◆ innerTrackphi_

◆ innerTrackpt_

◆ isItPlan1_