Inheritance diagram for HcalHBHEMuonHighEtaAnalyzer:

Public Member Functions
	HcalHBHEMuonHighEtaAnalyzer (const edm::ParameterSet &)

Public Member Functions inherited from edm::one::EDAnalyzer< edm::one::WatchRuns, edm::one::SharedResources >
	EDAnalyzer ()=default

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

Public Member Functions inherited from edm::one::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDAnalyzerBase () override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase &&)=default

	EDConsumerBase (EDConsumerBase const &)=delete

ESProxyIndex const *	esGetTokenIndices (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::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase &	operator= (EDConsumerBase &&)=default

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Static Public Member Functions inherited from edm::one::EDAnalyzerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Private Member Functions
double	activeLength (const DetId &)

void	analyze (edm::Event const &, edm::EventSetup const &) override

bool	analyzeHadron (edm::Event const &, math::XYZPoint &)

bool	analyzeMuon (edm::Event const &, math::XYZPoint &)

bool	analyzeTracks (const reco::Track *, math::XYZPoint &, int, std::vector< spr::propagatedTrackID > &, bool)

void	beginJob () override

void	beginRun (edm::Run const &, edm::EventSetup const &) override

void	clearVectors ()

int	depth16HE (int, int)

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

void	fillTrackParameters (const reco::Track *, math::XYZPoint)

double	gainFactor (const edm::ESHandle< HcalDbService > &, const HcalDetId &)

bool	goodCell (const HcalDetId &, const reco::Track , const CaloGeometry , const MagneticField *)

bool	isGoodVertex (const reco::Vertex &)

int	matchId (const HcalDetId &, const HcalDetId &)

double	respCorr (const DetId &)

Private Attributes
std::vector< HcalDDDRecConstants::HcalActiveLength >	actHB

std::vector< HcalDDDRecConstants::HcalActiveLength >	actHE

const bool	analyzeMuon_

edm::Handle< EcalRecHitCollection >	barrelRecHitsHandle_

edm::ESHandle< MagneticField >	bFieldH_

const bool	collapseDepth_

edm::ESHandle< HcalDbService >	conditions_

std::map< DetId, double >	corrValue_

std::vector< double >	ecal3x3Energy_

edm::ESHandle< EcalChannelStatus >	ecalChStatus_

std::vector< unsigned int >	ecalDetId_

std::vector< double >	ecalEnergy_

std::vector< unsigned int >	ehcalDetId_

std::vector< double >	emaxNearP_

const double	emaxNearPThr_

edm::Handle< EcalRecHitCollection >	endcapRecHitsHandle_

std::vector< double >	energyMuon_

std::vector< double >	etaGlob_

const double	etaMin_

unsigned int	eventNumber_

const std::string	fileInCorr_

const bool	getCharge_

unsigned int	goodVertex_

edm::Handle< HBHERecHitCollection >	hbhe_

std::vector< double >	hcal1x1Energy_

std::vector< int >	hcal_ieta_

std::vector< int >	hcal_iphi_

std::vector< double >	hcalActiveLength_

std::vector< double >	hcalActiveLengthHot_

std::vector< double >	hcalDepthActiveLength_ [depthMax_]

std::vector< double >	hcalDepthActiveLengthHot_ [depthMax_]

std::vector< double >	hcalDepthChargeHot_ [depthMax_]

std::vector< double >	hcalDepthChargeHotBG_ [depthMax_]

std::vector< double >	hcalDepthEnergy_ [depthMax_]

std::vector< double >	hcalDepthEnergyCorr_ [depthMax_]

std::vector< double >	hcalDepthEnergyHot_ [depthMax_]

std::vector< double >	hcalDepthEnergyHotCorr_ [depthMax_]

std::vector< bool >	hcalDepthMatch_ [depthMax_]

std::vector< bool >	hcalDepthMatchHot_ [depthMax_]

std::vector< unsigned int >	hcalDetId_

std::vector< double >	hcalEnergy_

std::vector< bool >	hcalHot_

const HcalDDDRecConstants *	hdc_

std::vector< double >	hoEnergy_

const bool	ignoreHECorr_

const bool	isItPlan1_

const bool	isItPreRecHit_

std::vector< double >	isolationR03_

std::vector< double >	isolationR04_

int	kount_

const edm::InputTag	labelEBRecHit_

const edm::InputTag	labelEERecHit_

const std::string	labelGenTrack_

const edm::InputTag	labelHBHERecHit_

const std::string	labelMuon_

const std::string	labelVtx_

std::vector< bool >	matchedId_

int	maxDepth_

std::vector< bool >	mediumMuon_

bool	mergedDepth_

edm::ESHandle< CaloGeometry >	pG_

std::vector< double >	phiGlob_

std::vector< double >	pMuon_

std::vector< double >	ptGlob_

HcalRespCorrs *	respCorrs_

unsigned int	runNumber_

spr::trackSelectionParameters	selectionParameter_

edm::ESHandle< EcalSeverityLevelAlgo >	sevlv_

edm::ESHandle< CaloTopology >	theCaloTopology_

const HcalTopology *	theHBHETopology_

const std::string	theTrackQuality_

edm::EDGetTokenT< EcalRecHitCollection >	tok_EB_

edm::EDGetTokenT< EcalRecHitCollection >	tok_EE_

edm::EDGetTokenT< reco::TrackCollection >	tok_genTrack_

edm::EDGetTokenT< HBHERecHitCollection >	tok_HBHE_

edm::EDGetTokenT< reco::MuonCollection >	tok_Muon_

edm::EDGetTokenT< reco::VertexCollection >	tok_Vtx_

std::vector< double >	trackDz_

std::vector< int >	trackLayerCrossed_

std::vector< int >	trackMissedInnerHits_

std::vector< int >	trackMissedOuterHits_

std::vector< int >	trackOuterHit_

TTree *	tree_

const bool	unCorrect_

bool	useMyCorr_

const int	useRaw_

const int	verbosity_

const bool	writeRespCorr_

Static Private Attributes
static const int	depthMax_ = 7

Additional Inherited Members
Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Protected Member Functions inherited from edm::EDConsumerBase
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

Detailed Description

Definition at line 76 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Constructor & Destructor Documentation

◆ HcalHBHEMuonHighEtaAnalyzer()

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

explicit

Definition at line 167 of file HcalHBHEMuonHighEtaAnalyzer.cc.

     : 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")),
       labelGenTrack_(iConfig.getParameter<std::string>("labelTrack")),
       etaMin_(iConfig.getParameter<double>("etaMin")),
       emaxNearPThr_(iConfig.getParameter<double>("emaxNearPThreshold")),
       analyzeMuon_(iConfig.getParameter<bool>("analyzeMuon")),
       unCorrect_(iConfig.getParameter<bool>("unCorrect")),
       collapseDepth_(iConfig.getParameter<bool>("collapseDepth")),
       isItPlan1_(iConfig.getParameter<bool>("isItPlan1")),
       getCharge_(iConfig.getParameter<bool>("getCharge")),
       useRaw_(iConfig.getParameter<int>("useRaw")),
       verbosity_(iConfig.getParameter<int>("verbosity")),
       theTrackQuality_(iConfig.getUntrackedParameter<std::string>("trackQuality")),
       fileInCorr_(iConfig.getUntrackedParameter<std::string>("fileInCorr", "")),
       ignoreHECorr_(iConfig.getUntrackedParameter<bool>("ignoreHECorr", false)),
       isItPreRecHit_(iConfig.getUntrackedParameter<bool>("isItPreRecHit", false)),
       writeRespCorr_(iConfig.getUntrackedParameter<bool>("writeRespCorr", false)),
       hdc_(nullptr),
       theHBHETopology_(nullptr),
       respCorrs_(nullptr),
       tree_(nullptr) {
   usesResource(TFileService::kSharedResource);
   //now do what ever initialization is needed
   kount_ = 0;
   maxDepth_ = iConfig.getUntrackedParameter<int>("maxDepth", 7);
   if (maxDepth_ > depthMax_)
     maxDepth_ = depthMax_;
   else if (maxDepth_ < 1)
     maxDepth_ = 4;
  
   reco::TrackBase::TrackQuality trackQuality = reco::TrackBase::qualityByName(theTrackQuality_);
   selectionParameter_.minPt = iConfig.getUntrackedParameter<double>("minTrackPt");
   selectionParameter_.minQuality = trackQuality;
   selectionParameter_.maxDxyPV = iConfig.getUntrackedParameter<double>("maxDxyPV");
   selectionParameter_.maxDzPV = iConfig.getUntrackedParameter<double>("maxDzPV");
   selectionParameter_.maxChi2 = iConfig.getUntrackedParameter<double>("maxChi2");
   selectionParameter_.maxDpOverP = iConfig.getUntrackedParameter<double>("maxDpOverP");
   selectionParameter_.minOuterHit = selectionParameter_.minLayerCrossed = 0;
   selectionParameter_.maxInMiss = selectionParameter_.maxOutMiss = 2;
  
   mergedDepth_ = (!isItPreRecHit_) || (collapseDepth_);
   tok_EB_ = consumes<EcalRecHitCollection>(labelEBRecHit_);
   tok_EE_ = consumes<EcalRecHitCollection>(labelEERecHit_);
   tok_HBHE_ = consumes<HBHERecHitCollection>(labelHBHERecHit_);
   tok_Vtx_ = consumes<reco::VertexCollection>(labelVtx_);
   tok_Muon_ = consumes<reco::MuonCollection>(labelMuon_);
   tok_genTrack_ = consumes<reco::TrackCollection>(labelGenTrack_);
   edm::LogVerbatim("HBHEMuon") << "Labels used: Track " << labelGenTrack_ << " Vtx " << labelVtx_ << " EB "
                                << labelEBRecHit_ << " EE " << labelEERecHit_ << " HBHE " << labelHBHERecHit_ << " MU "
                                << labelMuon_;
  
   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_;
 }

Member Function Documentation

◆ activeLength()

double HcalHBHEMuonHighEtaAnalyzer::activeLength ( const DetId & hid )

private

Definition at line 933 of file HcalHBHEMuonHighEtaAnalyzer.cc.

                                                                  {
   HcalDetId id(hid);
   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 (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;
       }
     }
   } 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;
       }
     }
   }
   return lx;
 }

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

Referenced by analyzeTracks(), and beginRun().

◆ analyze()

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

overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 342 of file HcalHBHEMuonHighEtaAnalyzer.cc.

                                                                                              {
   ++kount_;
   clearVectors();
   runNumber_ = iEvent.id().run();
   eventNumber_ = iEvent.id().event();
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HBHEMuon") << "Run " << runNumber_ << " Event " << eventNumber_;
 #endif
  
   // get handles to calogeometry and calotopology
   iSetup.get<CaloGeometryRecord>().get(pG_);
  
   iSetup.get<IdealMagneticFieldRecord>().get(bFieldH_);
   iSetup.get<EcalChannelStatusRcd>().get(ecalChStatus_);
   iSetup.get<EcalSeverityLevelAlgoRcd>().get(sevlv_);
   iSetup.get<CaloTopologyRecord>().get(theCaloTopology_);
   iSetup.get<HcalDbRecord>().get(conditions_);
  
   // Relevant blocks from iEvent
   edm::Handle<reco::VertexCollection> vtx;
   iEvent.getByToken(tok_Vtx_, vtx);
  
   iEvent.getByToken(tok_EB_, barrelRecHitsHandle_);
   iEvent.getByToken(tok_EE_, endcapRecHitsHandle_);
   iEvent.getByToken(tok_HBHE_, hbhe_);
  
   // require a good vertex
   math::XYZPoint pvx;
   goodVertex_ = 0;
   if (!vtx.isValid()) {
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("HBHEMuon") << "No Good Vertex found == Reject\n";
 #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 (barrelRecHitsHandle_.isValid() && endcapRecHitsHandle_.isValid() && hbhe_.isValid()) {
     accept = analyzeMuon_ ? analyzeMuon(iEvent, pvx) : analyzeHadron(iEvent, pvx);
   }
   if (accept) {
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("HBHEMuon") << "Total of " << hcal_ieta_.size() << " propagated points";
     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];
     if ((verbosity_ / 100) % 10 > 0) {
       edm::LogVerbatim("HBHEMuon") << "Sizes:: ptGlob:" << ptGlob_.size() << " etaGlob:" << etaGlob_.size()
                                    << " phiGlob:" << phiGlob_.size() << " energyMuon:" << energyMuon_.size()
                                    << " pMuon:" << pMuon_.size() << " mediumMuon: " << mediumMuon_.size()
                                    << " isolation:" << isolationR04_.size() << ":" << isolationR03_.size()
                                    << " e|h|ho energy: " << ecalEnergy_.size() << ":" << hcalEnergy_.size() << ":"
                                    << hoEnergy_.size();
       edm::LogVerbatim("HBHEMuon") << "        matchedId:" << matchedId_.size() << " hcalHot:" << hcalHot_.size()
                                    << " 3x3|1x1 energy:" << ecal3x3Energy_.size() << ":" << hcal1x1Energy_.size()
                                    << " detId:" << ecalDetId_.size() << ":" << hcalDetId_.size() << ":"
                                    << ehcalDetId_.size() << " eta|phi:" << hcal_ieta_.size() << ":"
                                    << hcal_iphi_.size();
       edm::LogVerbatim("HBHEMuon") << "        activeLength:" << hcalActiveLength_.size() << ":"
                                    << hcalActiveLengthHot_.size() << " emaxNearP:" << emaxNearP_.size()
                                    << " trackDz: " << trackDz_.size() << " tracks:" << trackLayerCrossed_.size() << ":"
                                    << trackOuterHit_.size() << ":" << trackMissedInnerHits_.size() << ":"
                                    << trackMissedOuterHits_.size();
       for (unsigned int i = 0; i < depthMax_; ++i)
         edm::LogVerbatim("HBHEMuon")
             << "Depth " << i
             << " Energy|Length|EnergyHot|LengthHot|Charge|ChargeBG|EnergyCorr|EnergyHotCorr|Match|MatchHot:"
             << hcalDepthEnergy_[i].size() << ":" << hcalDepthActiveLength_[i].size() << ":"
             << hcalDepthEnergyHot_[i].size() << ":" << hcalDepthActiveLengthHot_[i].size() << ":"
             << hcalDepthChargeHot_[i].size() << ":" << hcalDepthChargeHotBG_[i].size() << ":"
             << hcalDepthEnergyCorr_[i].size() << ":" << hcalDepthEnergyHotCorr_[i].size() << ":"
             << hcalDepthMatch_[i].size() << ":" << hcalDepthMatchHot_[i].size();
     }
 #endif
     tree_->Fill();
   }
 }

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

◆ analyzeHadron()

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

private

Definition at line 557 of file HcalHBHEMuonHighEtaAnalyzer.cc.

                                                                                             {
   //Get track collection
   edm::Handle<reco::TrackCollection> trkCollection;
   iEvent.getByToken(tok_genTrack_, trkCollection);
   bool accept = false;
  
   if (!trkCollection.isValid()) {
     const CaloGeometry* geo = pG_.product();
     const MagneticField* bField = bFieldH_.product();
     std::vector<spr::propagatedTrackID> trkCaloDets;
     spr::propagateCALO(trkCollection, geo, bField, theTrackQuality_, trkCaloDets, false);
     int nTrack(0);
     std::vector<spr::propagatedTrackID>::const_iterator trkDetItr;
     for (trkDetItr = trkCaloDets.begin(), nTrack = 0; trkDetItr != trkCaloDets.end(); trkDetItr++, nTrack++) {
       const reco::Track* pTrack = &(*(trkDetItr->trkItr));
       if (std::abs(pTrack->eta()) > etaMin_) {
         accept = analyzeTracks(pTrack, leadPV, nTrack, trkCaloDets, true);
       }
     }
   }
   return accept;
 }

References funct::abs(), accept(), analyzeTracks(), Calorimetry_cff::bField, bFieldH_, reco::TrackBase::eta(), etaMin_, iEvent, edm::HandleBase::isValid(), pG_, edm::ESHandle< T >::product(), spr::propagateCALO(), theTrackQuality_, and tok_genTrack_.

Referenced by analyze().

◆ analyzeMuon()

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

private

Definition at line 492 of file HcalHBHEMuonHighEtaAnalyzer.cc.

                                                                                           {
   edm::Handle<reco::MuonCollection> _Muon;
   iEvent.getByToken(tok_Muon_, _Muon);
   bool accept = false;
  
   if (_Muon.isValid()) {
     int nTrack(0);
     std::vector<spr::propagatedTrackID> trkCaloDets;
     for (reco::MuonCollection::const_iterator RecMuon = _Muon->begin(); RecMuon != _Muon->end(); ++RecMuon) {
       if (RecMuon->innerTrack().isNonnull()) {
         const reco::Track* pTrack = (RecMuon->innerTrack()).get();
         if (std::abs(pTrack->eta()) > etaMin_) {
           if (analyzeTracks(pTrack, leadPV, nTrack, trkCaloDets, false)) {
             accept = true;
             ptGlob_.emplace_back((RecMuon)->pt());
             etaGlob_.emplace_back(RecMuon->eta());
             phiGlob_.emplace_back(RecMuon->phi());
             energyMuon_.push_back(RecMuon->energy());
             pMuon_.emplace_back(RecMuon->p());
             bool mediumMuon = (((RecMuon->isPFMuon()) && (RecMuon->isGlobalMuon() || RecMuon->isTrackerMuon())) &&
                                (RecMuon->innerTrack()->validFraction() > 0.49));
             if (mediumMuon) {
               double chiGlobal =
                   ((RecMuon->globalTrack().isNonnull()) ? RecMuon->globalTrack()->normalizedChi2() : 999);
               bool goodGlob =
                   (RecMuon->isGlobalMuon() && chiGlobal < 3 && RecMuon->combinedQuality().chi2LocalPosition < 12 &&
                    RecMuon->combinedQuality().trkKink < 20);
               mediumMuon = muon::segmentCompatibility(*RecMuon) > (goodGlob ? 0.303 : 0.451);
             }
             mediumMuon_.emplace_back(mediumMuon);
             bool isoR03 =
                 ((RecMuon->pfIsolationR03().sumChargedHadronPt +
                   std::max(0.,
                            RecMuon->pfIsolationR03().sumNeutralHadronEt + RecMuon->pfIsolationR03().sumPhotonEt -
                                (0.5 * RecMuon->pfIsolationR03().sumPUPt))) /
                  RecMuon->pt());
             bool isoR04 =
                 ((RecMuon->pfIsolationR04().sumChargedHadronPt +
                   std::max(0.,
                            RecMuon->pfIsolationR04().sumNeutralHadronEt + RecMuon->pfIsolationR04().sumPhotonEt -
                                (0.5 * RecMuon->pfIsolationR04().sumPUPt))) /
                  RecMuon->pt());
             isolationR03_.emplace_back(isoR03);
             isolationR04_.emplace_back(isoR04);
  
             ecalEnergy_.emplace_back(RecMuon->calEnergy().emS9);
             hcalEnergy_.emplace_back(RecMuon->calEnergy().hadS9);
             hoEnergy_.emplace_back(RecMuon->calEnergy().hoS9);
 #ifdef EDM_ML_DEBUG
             if ((verbosity_ / 100) % 10 > 0)
               edm::LogVerbatim("HBHEMuon")
                   << "Muon[" << ptGlob_.size() << "] pt:eta:phi:p " << ptGlob_.back() << ":" << etaGlob_.back() << ":"
                   << phiGlob_.back() << ":" << energyMuon_.back() << ":" << pMuon_.back() << ":"
                   << " Medium:i3:i4 " << mediumMuon_.back() << ":" << isolationR03_.back() << ":"
                   << isolationR04_.back() << ":"
                   << " Energy EC:HC:HO " << ecalEnergy_.back() << ":" << hcalEnergy_.back() << ":" << hoEnergy_.back();
 #endif
           }
         }
       }
     }
   }
   return accept;
 }

References funct::abs(), accept(), analyzeTracks(), ecalEnergy_, energyMuon_, reco::TrackBase::eta(), etaGlob_, etaMin_, get, hcalEnergy_, hoEnergy_, iEvent, isolationR03_, isolationR04_, edm::HandleBase::isValid(), SiStripPI::max, mediumMuon_, phiGlob_, pMuon_, DiDispStaMuonMonitor_cfi::pt, ptGlob_, muon::segmentCompatibility(), tok_Muon_, and verbosity_.

Referenced by analyze().

◆ analyzeTracks()

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

private

Definition at line 580 of file HcalHBHEMuonHighEtaAnalyzer.cc.

                                                                {
   const CaloGeometry* geo = pG_.product();
   const MagneticField* bField = bFieldH_.product();
   const EcalChannelStatus* theEcalChStatus = ecalChStatus_.product();
   const CaloTopology* caloTopology = theCaloTopology_.product();
   bool accept(false);
  
   if (spr::goodTrack(pTrack, leadPV, selectionParameter_, false)) {
     spr::propagatedTrackID trackID = spr::propagateCALO(pTrack, geo, bField, false);
  
     if (trackID.okECAL && trackID.okHCAL) {
       double emaxNearP = (ifHadron) ? spr::chargeIsolationEcal(nTrack, trkCaloDets, geo, caloTopology, 15, 15) : 0;
       if (emaxNearP < emaxNearPThr_) {
         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;
         }
  
         HcalDetId check;
         std::pair<bool, HcalDetId> info = spr::propagateHCALBack(pTrack, geo, bField, false);
         if (info.first)
           check = info.second;
  
         const DetId isoCell(trackID.detIdECAL);
         std::pair<double, bool> e3x3 = spr::eECALmatrix(isoCell,
                                                         barrelRecHitsHandle_,
                                                         endcapRecHitsHandle_,
                                                         *theEcalChStatus,
                                                         geo,
                                                         caloTopology,
                                                         sevlv_.product(),
                                                         1,
                                                         1,
                                                         -100.0,
                                                         -100.0,
                                                         -500.0,
                                                         500.0,
                                                         false);
         eEcal = e3x3.first;
 #ifdef EDM_ML_DEBUG
         if (verbosity_ % 10 > 0)
           edm::LogVerbatim("HBHEMuon") << "Propagate Track to ECAL: " << e3x3.second << ":" << trackID.okECAL << " E "
                                        << eEcal;
 #endif
  
         DetId closestCell(trackID.detIdHCAL);
         HcalDetId hcidt(closestCell.rawId());
         if ((hcidt.ieta() == check.ieta()) && (hcidt.iphi() == check.iphi()))
           tmpmatch = true;
 #ifdef EDM_ML_DEBUG
         if (verbosity_ % 10 > 0)
           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,
                             depthMax_,
                             -100.0,
                             -100.0,
                             -100.0,
                             -100.0,
                             -500.0,
                             500.0,
                             useRaw_,
                             depth16HE(ieta, iphi),
                             false);
         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));
           double ene = ehdepth[i].first;
           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
                 if (verbosity_ % 10 > 0) {
                   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) % 10 > 0)
                 edm::LogVerbatim("HBHEMuon")
                     << hcid0 << " E " << ene << ":" << enec << " L " << actL << " Match " << tmpC;
 #endif
             }
           }
         }
 #ifdef EDM_ML_DEBUG
         if ((verbosity_ / 10) % 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 << "] " << matchDepth[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,
                               depthMax_,
                               -100.0,
                               -100.0,
                               -100.0,
                               -100.0,
                               -500.0,
                               500.0,
                               useRaw_,
                               depth16HE(ieta, iphi),
                               false);
           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));
             double ene = ehdepth[i].first;
             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
                   if (verbosity_ % 10 > 0) {
                     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
                   if (verbosity_ % 10 > 0)
                     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) % 10 > 0)
                   edm::LogVerbatim("HBHEMuon") << hcid0 << " depth " << depth << " E " << ene << ":" << enec << " C "
                                                << chg << " L " << actL << " Match " << tmpC;
 #endif
               }
             }
           }
         }
 #ifdef EDM_ML_DEBUG
         edm::LogVerbatim("HBHEMuon") << "Propagate Track to HCAL: " << trackID.okHCAL << " Match " << tmpmatch
                                      << " Hot " << isHot << " Energy " << eHcal;
 #endif
  
         accept = true;
         ecalDetId_.emplace_back((trackID.detIdECAL)());
         hcalDetId_.emplace_back((trackID.detIdHCAL)());
         ehcalDetId_.emplace_back((trackID.detIdEHCAL)());
         emaxNearP_.emplace_back(emaxNearP);
         matchedId_.emplace_back(tmpmatch);
         ecal3x3Energy_.emplace_back(eEcal);
         hcal1x1Energy_.emplace_back(eHcal);
         hcal_ieta_.emplace_back(ieta);
         hcal_iphi_.emplace_back(iphi);
         for (int i = 0; i < maxDepth_; ++i) {
           hcalDepthEnergy_[i].emplace_back(eHcalDepth[i]);
           hcalDepthActiveLength_[i].emplace_back(activeL[i]);
           hcalDepthEnergyHot_[i].emplace_back(eHcalDepthHot[i]);
           hcalDepthActiveLengthHot_[i].emplace_back(activeHotL[i]);
           hcalDepthEnergyCorr_[i].emplace_back(eHcalDepthC[i]);
           hcalDepthEnergyHotCorr_[i].emplace_back(eHcalDepthHotC[i]);
           hcalDepthChargeHot_[i].emplace_back(cHcalDepthHot[i]);
           hcalDepthChargeHotBG_[i].emplace_back(cHcalDepthHotBG[i]);
           hcalDepthMatch_[i].emplace_back(matchDepth[i]);
           hcalDepthMatchHot_[i].emplace_back(matchDepthHot[i]);
         }
         hcalActiveLength_.emplace_back(activeLengthTot);
         hcalHot_.emplace_back(isHot);
         hcalActiveLengthHot_.emplace_back(activeLengthHotTot);
 #ifdef EDM_ML_DEBUG
         if ((verbosity_ / 100) % 10 > 0) {
           edm::LogVerbatim("HBHEMuon") << "Track " << std::hex << ecalDetId_.back() << ":" << hcalDetId_.back() << ":"
                                        << ehcalDetId_.back() << std::dec << ":" << emaxNearP_.back() << ":"
                                        << matchedId_.back() << ":" << ecal3x3Energy_.back() << ":"
                                        << hcal1x1Energy_.back() << ":" << hcal_ieta_.back() << ":" << hcal_iphi_.back()
                                        << ":" << hcalActiveLength_.back() << ":" << hcalHot_.back() << ":"
                                        << hcalActiveLengthHot_.back();
           for (int i = 0; i < maxDepth_; ++i) {
             edm::LogVerbatim("HBHEMuon") << "Depth[" << i << "] " << hcalDepthEnergy_[i].back() << ":"
                                          << hcalDepthActiveLength_[i].back() << ":" << hcalDepthEnergyHot_[i].back()
                                          << ":" << hcalDepthActiveLengthHot_[i].back() << ":"
                                          << hcalDepthEnergyCorr_[i].back() << ":" << hcalDepthEnergyHotCorr_[i].back()
                                          << ":" << hcalDepthChargeHot_[i].back() << ":"
                                          << hcalDepthChargeHotBG_[i].back() << ":" << hcalDepthMatch_[i].back() << ":"
                                          << hcalDepthMatchHot_[i].back();
           }
         }
 #endif
         fillTrackParameters(pTrack, leadPV);
       }
     }
   }
   return accept;
 }

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

Referenced by analyzeHadron(), and analyzeMuon().

◆ beginJob()

void HcalHBHEMuonHighEtaAnalyzer::beginJob ( void )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 280 of file HcalHBHEMuonHighEtaAnalyzer.cc.

                                            {
   edm::Service<TFileService> fs;
   tree_ = fs->make<TTree>("HBHEMuonHighEta", "HBHEMuonHighEta");
   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("MediumMuon", &mediumMuon_);
   tree_->Branch("IsolationR04", &isolationR04_);
   tree_->Branch("IsolationR03", &isolationR03_);
   tree_->Branch("ecal_3into3", &ecalEnergy_);
   tree_->Branch("hcal_3into3", &hcalEnergy_);
   tree_->Branch("ho_3into3", &hoEnergy_);
   tree_->Branch("emaxNearP", &emaxNearP_);
  
   tree_->Branch("Run_No", &runNumber_);
   tree_->Branch("Event_No", &eventNumber_);
   tree_->Branch("GoodVertex", &goodVertex_);
   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("trackDz", &trackDz_);
   tree_->Branch("trackLayerCrossed", &trackLayerCrossed_);
   tree_->Branch("trackOuterHit", &trackOuterHit_);
   tree_->Branch("trackMissedInnerHits", &trackMissedInnerHits_);
   tree_->Branch("trackMissedOuterHits", &trackMissedOuterHits_);
 }

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

◆ beginRun()

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

overrideprivate

Definition at line 432 of file HcalHBHEMuonHighEtaAnalyzer.cc.

                                                                                           {
   edm::ESHandle<HcalDDDRecConstants> pHRNDC;
   iSetup.get<HcalRecNumberingRecord>().get(pHRNDC);
   hdc_ = pHRNDC.product();
   actHB.clear();
   actHE.clear();
   actHB = hdc_->getThickActive(0);
   actHE = hdc_->getThickActive(1);
 #ifdef EDM_ML_DEBUG
   if (verbosity_ % 10 > 0) {
     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));
       ++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));
       ++k2;
     }
   }
 #endif
  
   edm::ESHandle<HcalTopology> htopo;
   iSetup.get<HcalRecNumberingRecord>().get(htopo);
   theHBHETopology_ = htopo.product();
  
   edm::ESHandle<HcalRespCorrs> resp;
   iSetup.get<HcalRespCorrsRcd>().get(resp);
   respCorrs_ = new HcalRespCorrs(*resp.product());
   respCorrs_->setTopo(theHBHETopology_);
  
   // Write correction factors for all HB/HE events
   if (writeRespCorr_) {
     edm::ESHandle<CaloGeometry> pG;
     iSetup.get<CaloGeometryRecord>().get(pG);
     const CaloGeometry* geo = pG.product();
     const HcalGeometry* gHcal = (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();
       }
     }
   }
 }

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

◆ clearVectors()

void HcalHBHEMuonHighEtaAnalyzer::clearVectors ( )

private

clearing vectots

Definition at line 876 of file HcalHBHEMuonHighEtaAnalyzer.cc.

                                                {
   eventNumber_ = -99999;
   runNumber_ = -99999;
   goodVertex_ = -99999;
  
   mediumMuon_.clear();
   ptGlob_.clear();
   etaGlob_.clear();
   phiGlob_.clear();
   energyMuon_.clear();
   pMuon_.clear();
   isolationR04_.clear();
   isolationR03_.clear();
   ecalEnergy_.clear();
   hcalEnergy_.clear();
   hoEnergy_.clear();
  
   matchedId_.clear();
   hcalHot_.clear();
   ecal3x3Energy_.clear();
   hcal1x1Energy_.clear();
   ecalDetId_.clear();
   hcalDetId_.clear();
   ehcalDetId_.clear();
   hcal_ieta_.clear();
   hcal_iphi_.clear();
   for (int i = 0; i < depthMax_; ++i) {
     hcalDepthEnergy_[i].clear();
     hcalDepthActiveLength_[i].clear();
     hcalDepthEnergyHot_[i].clear();
     hcalDepthActiveLengthHot_[i].clear();
     hcalDepthChargeHot_[i].clear();
     hcalDepthChargeHotBG_[i].clear();
     hcalDepthEnergyCorr_[i].clear();
     hcalDepthEnergyHotCorr_[i].clear();
     hcalDepthMatch_[i].clear();
     hcalDepthMatchHot_[i].clear();
   }
   hcalActiveLength_.clear();
   hcalActiveLengthHot_.clear();
  
   emaxNearP_.clear();
   trackDz_.clear();
   trackLayerCrossed_.clear();
   trackOuterHit_.clear();
   trackMissedInnerHits_.clear();
   trackMissedOuterHits_.clear();
 }

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

Referenced by analyze().

◆ depth16HE()

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

private

Definition at line 1000 of file HcalHBHEMuonHighEtaAnalyzer.cc.

                                                              {
   // 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
   if (verbosity_ % 10 > 0)
     edm::LogVerbatim("HBHEMuon") << "Plan1 " << isItPlan1_ << " PreRecHit " << isItPreRecHit_ << " phi " << iphi
                                  << " depth " << depth;
 #endif
   return depth;
 }

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

Referenced by analyzeTracks().

◆ endRun()

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

inlineoverrideprivate

Definition at line 86 of file HcalHBHEMuonHighEtaAnalyzer.cc.

86 {}

◆ fillDescriptions()

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

static

Definition at line 248 of file HcalHBHEMuonHighEtaAnalyzer.cc.

                                                                                              {
   edm::ParameterSetDescription desc;
   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");
   desc.add<std::string>("labelTrack", "generalTracks");
   desc.add<double>("etaMin", 2.0);
   desc.add<double>("emaxNearPThreshold", 10.0);
   desc.add<bool>("analyzeMuon", true);
   desc.add<bool>("unCorrect", false);
   desc.add<bool>("collapseDepth", false);
   desc.add<bool>("isItPlan1", false);
   desc.add<bool>("getCharge", false);
   desc.add<int>("useRaw", 0);
   desc.add<int>("verbosity", 0);
   desc.addUntracked<std::string>("fileInCorr", "");
   desc.addUntracked<std::string>("trackQuality", "highPurity");
   desc.addUntracked<double>("minTrackPt", 1.0);
   desc.addUntracked<double>("maxDxyPV", 0.02);
   desc.addUntracked<double>("maxDzPV", 100.0);
   desc.addUntracked<double>("maxChi2", 5.0);
   desc.addUntracked<double>("maxDpOverP", 0.1);
   desc.addUntracked<bool>("ignoreHECorr", false);
   desc.addUntracked<bool>("isItPreRecHit", false);
   desc.addUntracked<bool>("writeRespCorr", false);
   desc.addUntracked<int>("maxDepth", 7);
   descriptions.add("hcalHBHEMuonHighEta", desc);
 }

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), edm::ParameterSetDescription::addUntracked(), HLT_2018_cff::InputTag, and AlCaHLTBitMon_QueryRunRegistry::string.

◆ fillTrackParameters()

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

private

Definition at line 1026 of file HcalHBHEMuonHighEtaAnalyzer.cc.

                                                                                                   {
   trackDz_.emplace_back(pTrack->dz(leadPV));
   const reco::HitPattern& hitp = pTrack->hitPattern();
   trackLayerCrossed_.emplace_back(hitp.trackerLayersWithMeasurement());
   trackOuterHit_.emplace_back(hitp.stripTOBLayersWithMeasurement() + hitp.stripTECLayersWithMeasurement());
   trackMissedInnerHits_.emplace_back(hitp.trackerLayersWithoutMeasurement(reco::HitPattern::MISSING_INNER_HITS));
   trackMissedOuterHits_.emplace_back(hitp.trackerLayersWithoutMeasurement(reco::HitPattern::MISSING_OUTER_HITS));
 }

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

Referenced by analyzeTracks().

◆ gainFactor()

double HcalHBHEMuonHighEtaAnalyzer::gainFactor	(	const edm::ESHandle< HcalDbService > &	conditions,
		const HcalDetId &	id
	)

private

Definition at line 992 of file HcalHBHEMuonHighEtaAnalyzer.cc.

                                                                                                                 {
   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;
 }

References PedestalClient_cfi::gain, HcalDbService::getHcalCalibrations(), and HcalCalibrations::respcorrgain().

Referenced by analyzeTracks().

◆ goodCell()

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

private

Definition at line 1016 of file HcalHBHEMuonHighEtaAnalyzer.cc.

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

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

Referenced by analyzeTracks().

◆ isGoodVertex()

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

private

Definition at line 968 of file HcalHBHEMuonHighEtaAnalyzer.cc.

                                                                     {
   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;
 }

References badGlobalMuonTaggersAOD_cff::vtx.

Referenced by analyze().

◆ matchId()

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

private

Definition at line 926 of file HcalHBHEMuonHighEtaAnalyzer.cc.

                                                                                    {
   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;
 }

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

Referenced by analyzeTracks().

◆ respCorr()

double HcalHBHEMuonHighEtaAnalyzer::respCorr ( const DetId & id )

private

Definition at line 980 of file HcalHBHEMuonHighEtaAnalyzer.cc.

                                                             {
   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;
 }

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

Referenced by analyzeTracks().

Member Data Documentation

◆ actHB

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

private

Definition at line 162 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by activeLength(), and beginRun().

◆ actHE

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

private

Definition at line 162 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by activeLength(), and beginRun().

◆ analyzeMuon_

const bool HcalHBHEMuonHighEtaAnalyzer::analyzeMuon_

private

Definition at line 105 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyze().

◆ barrelRecHitsHandle_

edm::Handle<EcalRecHitCollection> HcalHBHEMuonHighEtaAnalyzer::barrelRecHitsHandle_

private

Definition at line 131 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyze(), and analyzeTracks().

◆ bFieldH_

edm::ESHandle<MagneticField> HcalHBHEMuonHighEtaAnalyzer::bFieldH_

private

Definition at line 125 of file HcalHBHEMuonHighEtaAnalyzer.cc.

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

◆ collapseDepth_

const bool HcalHBHEMuonHighEtaAnalyzer::collapseDepth_

private

Definition at line 105 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyzeTracks(), and HcalHBHEMuonHighEtaAnalyzer().

◆ conditions_

edm::ESHandle<HcalDbService> HcalHBHEMuonHighEtaAnalyzer::conditions_

private

Definition at line 129 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyze(), and analyzeTracks().

◆ corrValue_

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

private

Definition at line 163 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by HcalHBHEMuonHighEtaAnalyzer(), and respCorr().

◆ depthMax_

const int HcalHBHEMuonHighEtaAnalyzer::depthMax_ = 7

staticprivate

Definition at line 136 of file HcalHBHEMuonHighEtaAnalyzer.cc.

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

◆ ecal3x3Energy_

std::vector<double> HcalHBHEMuonHighEtaAnalyzer::ecal3x3Energy_

private

Definition at line 144 of file HcalHBHEMuonHighEtaAnalyzer.cc.

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

◆ ecalChStatus_

edm::ESHandle<EcalChannelStatus> HcalHBHEMuonHighEtaAnalyzer::ecalChStatus_

private

Definition at line 126 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyze(), and analyzeTracks().

◆ ecalDetId_

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

private

Definition at line 145 of file HcalHBHEMuonHighEtaAnalyzer.cc.

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

◆ ecalEnergy_

std::vector<double> HcalHBHEMuonHighEtaAnalyzer::ecalEnergy_

private

Definition at line 142 of file HcalHBHEMuonHighEtaAnalyzer.cc.

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

◆ ehcalDetId_

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

private

Definition at line 145 of file HcalHBHEMuonHighEtaAnalyzer.cc.

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

◆ emaxNearP_

std::vector<double> HcalHBHEMuonHighEtaAnalyzer::emaxNearP_

private

Definition at line 158 of file HcalHBHEMuonHighEtaAnalyzer.cc.

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

◆ emaxNearPThr_

const double HcalHBHEMuonHighEtaAnalyzer::emaxNearPThr_

private

Definition at line 104 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyzeTracks().

◆ endcapRecHitsHandle_

edm::Handle<EcalRecHitCollection> HcalHBHEMuonHighEtaAnalyzer::endcapRecHitsHandle_

private

Definition at line 132 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyze(), and analyzeTracks().

◆ energyMuon_

std::vector<double> HcalHBHEMuonHighEtaAnalyzer::energyMuon_

private

Definition at line 140 of file HcalHBHEMuonHighEtaAnalyzer.cc.

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

◆ etaGlob_

std::vector<double> HcalHBHEMuonHighEtaAnalyzer::etaGlob_

private

Definition at line 140 of file HcalHBHEMuonHighEtaAnalyzer.cc.

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

◆ etaMin_

const double HcalHBHEMuonHighEtaAnalyzer::etaMin_

private

Definition at line 104 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyzeHadron(), and analyzeMuon().

◆ eventNumber_

unsigned int HcalHBHEMuonHighEtaAnalyzer::eventNumber_

private

Definition at line 138 of file HcalHBHEMuonHighEtaAnalyzer.cc.

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

◆ fileInCorr_

const std::string HcalHBHEMuonHighEtaAnalyzer::fileInCorr_

private

Definition at line 107 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by HcalHBHEMuonHighEtaAnalyzer().

◆ getCharge_

const bool HcalHBHEMuonHighEtaAnalyzer::getCharge_

private

Definition at line 105 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyzeTracks(), and HcalHBHEMuonHighEtaAnalyzer().

◆ goodVertex_

unsigned int HcalHBHEMuonHighEtaAnalyzer::goodVertex_

private

Definition at line 138 of file HcalHBHEMuonHighEtaAnalyzer.cc.

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

◆ hbhe_

edm::Handle<HBHERecHitCollection> HcalHBHEMuonHighEtaAnalyzer::hbhe_

private

Definition at line 133 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyze(), and analyzeTracks().

◆ hcal1x1Energy_

std::vector<double> HcalHBHEMuonHighEtaAnalyzer::hcal1x1Energy_

private

Definition at line 144 of file HcalHBHEMuonHighEtaAnalyzer.cc.

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

◆ hcal_ieta_

std::vector<int> HcalHBHEMuonHighEtaAnalyzer::hcal_ieta_

private

Definition at line 146 of file HcalHBHEMuonHighEtaAnalyzer.cc.

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

◆ hcal_iphi_

std::vector<int> HcalHBHEMuonHighEtaAnalyzer::hcal_iphi_

private

Definition at line 146 of file HcalHBHEMuonHighEtaAnalyzer.cc.

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

◆ hcalActiveLength_

std::vector<double> HcalHBHEMuonHighEtaAnalyzer::hcalActiveLength_

private

Definition at line 157 of file HcalHBHEMuonHighEtaAnalyzer.cc.

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

◆ hcalActiveLengthHot_

std::vector<double> HcalHBHEMuonHighEtaAnalyzer::hcalActiveLengthHot_

private

Definition at line 157 of file HcalHBHEMuonHighEtaAnalyzer.cc.

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

◆ hcalDepthActiveLength_

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

private

Definition at line 148 of file HcalHBHEMuonHighEtaAnalyzer.cc.

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

◆ hcalDepthActiveLengthHot_

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

private

Definition at line 150 of file HcalHBHEMuonHighEtaAnalyzer.cc.

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

◆ hcalDepthChargeHot_

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

private

Definition at line 151 of file HcalHBHEMuonHighEtaAnalyzer.cc.

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

◆ hcalDepthChargeHotBG_

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

private

Definition at line 152 of file HcalHBHEMuonHighEtaAnalyzer.cc.

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

◆ hcalDepthEnergy_

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

private

Definition at line 147 of file HcalHBHEMuonHighEtaAnalyzer.cc.

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

◆ hcalDepthEnergyCorr_

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

private

Definition at line 153 of file HcalHBHEMuonHighEtaAnalyzer.cc.

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

◆ hcalDepthEnergyHot_

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

private

Definition at line 149 of file HcalHBHEMuonHighEtaAnalyzer.cc.

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

◆ hcalDepthEnergyHotCorr_

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

private

Definition at line 154 of file HcalHBHEMuonHighEtaAnalyzer.cc.

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

◆ hcalDepthMatch_

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

private

Definition at line 155 of file HcalHBHEMuonHighEtaAnalyzer.cc.

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

◆ hcalDepthMatchHot_

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

private

Definition at line 156 of file HcalHBHEMuonHighEtaAnalyzer.cc.

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

◆ hcalDetId_

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

private

Definition at line 145 of file HcalHBHEMuonHighEtaAnalyzer.cc.

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

◆ hcalEnergy_

std::vector<double> HcalHBHEMuonHighEtaAnalyzer::hcalEnergy_

private

Definition at line 142 of file HcalHBHEMuonHighEtaAnalyzer.cc.

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

◆ hcalHot_

std::vector<bool> HcalHBHEMuonHighEtaAnalyzer::hcalHot_

private

Definition at line 143 of file HcalHBHEMuonHighEtaAnalyzer.cc.

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

◆ hdc_

const HcalDDDRecConstants* HcalHBHEMuonHighEtaAnalyzer::hdc_

private

Definition at line 113 of file HcalHBHEMuonHighEtaAnalyzer.cc.

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

◆ hoEnergy_

std::vector<double> HcalHBHEMuonHighEtaAnalyzer::hoEnergy_

private

Definition at line 142 of file HcalHBHEMuonHighEtaAnalyzer.cc.

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

◆ ignoreHECorr_

const bool HcalHBHEMuonHighEtaAnalyzer::ignoreHECorr_

private

Definition at line 108 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyzeTracks(), and HcalHBHEMuonHighEtaAnalyzer().

◆ isItPlan1_

const bool HcalHBHEMuonHighEtaAnalyzer::isItPlan1_

private

Definition at line 105 of file HcalHBHEMuonHighEtaAnalyzer.cc.

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

◆ isItPreRecHit_

const bool HcalHBHEMuonHighEtaAnalyzer::isItPreRecHit_

private

Definition at line 108 of file HcalHBHEMuonHighEtaAnalyzer.cc.

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

◆ isolationR03_

std::vector<double> HcalHBHEMuonHighEtaAnalyzer::isolationR03_

private

Definition at line 141 of file HcalHBHEMuonHighEtaAnalyzer.cc.

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

◆ isolationR04_

std::vector<double> HcalHBHEMuonHighEtaAnalyzer::isolationR04_

private

Definition at line 141 of file HcalHBHEMuonHighEtaAnalyzer.cc.

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

◆ kount_

int HcalHBHEMuonHighEtaAnalyzer::kount_

private

Definition at line 110 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyze(), and HcalHBHEMuonHighEtaAnalyzer().

◆ labelEBRecHit_

const edm::InputTag HcalHBHEMuonHighEtaAnalyzer::labelEBRecHit_

private

Definition at line 102 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by HcalHBHEMuonHighEtaAnalyzer().

◆ labelEERecHit_

const edm::InputTag HcalHBHEMuonHighEtaAnalyzer::labelEERecHit_

private

Definition at line 102 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by HcalHBHEMuonHighEtaAnalyzer().

◆ labelGenTrack_

const std::string HcalHBHEMuonHighEtaAnalyzer::labelGenTrack_

private

Definition at line 103 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by HcalHBHEMuonHighEtaAnalyzer().

◆ labelHBHERecHit_

const edm::InputTag HcalHBHEMuonHighEtaAnalyzer::labelHBHERecHit_

private

Definition at line 102 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by HcalHBHEMuonHighEtaAnalyzer().

◆ labelMuon_

const std::string HcalHBHEMuonHighEtaAnalyzer::labelMuon_

private

Definition at line 103 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by HcalHBHEMuonHighEtaAnalyzer().

◆ labelVtx_

const std::string HcalHBHEMuonHighEtaAnalyzer::labelVtx_

private

Definition at line 103 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by HcalHBHEMuonHighEtaAnalyzer().

◆ matchedId_

std::vector<bool> HcalHBHEMuonHighEtaAnalyzer::matchedId_

private

Definition at line 143 of file HcalHBHEMuonHighEtaAnalyzer.cc.

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

◆ maxDepth_

int HcalHBHEMuonHighEtaAnalyzer::maxDepth_

private

Definition at line 110 of file HcalHBHEMuonHighEtaAnalyzer.cc.

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

◆ mediumMuon_

std::vector<bool> HcalHBHEMuonHighEtaAnalyzer::mediumMuon_

private

Definition at line 139 of file HcalHBHEMuonHighEtaAnalyzer.cc.

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

◆ mergedDepth_

bool HcalHBHEMuonHighEtaAnalyzer::mergedDepth_

private

Definition at line 109 of file HcalHBHEMuonHighEtaAnalyzer.cc.

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

◆ pG_

edm::ESHandle<CaloGeometry> HcalHBHEMuonHighEtaAnalyzer::pG_

private

Definition at line 124 of file HcalHBHEMuonHighEtaAnalyzer.cc.

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

◆ phiGlob_

std::vector<double> HcalHBHEMuonHighEtaAnalyzer::phiGlob_

private

Definition at line 140 of file HcalHBHEMuonHighEtaAnalyzer.cc.

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

◆ pMuon_

std::vector<double> HcalHBHEMuonHighEtaAnalyzer::pMuon_

private

Definition at line 140 of file HcalHBHEMuonHighEtaAnalyzer.cc.

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

◆ ptGlob_

std::vector<double> HcalHBHEMuonHighEtaAnalyzer::ptGlob_

private

Definition at line 140 of file HcalHBHEMuonHighEtaAnalyzer.cc.

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

◆ respCorrs_

HcalRespCorrs* HcalHBHEMuonHighEtaAnalyzer::respCorrs_

private

Definition at line 115 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by beginRun(), and respCorr().

◆ runNumber_

unsigned int HcalHBHEMuonHighEtaAnalyzer::runNumber_

private

Definition at line 138 of file HcalHBHEMuonHighEtaAnalyzer.cc.

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

◆ selectionParameter_

spr::trackSelectionParameters HcalHBHEMuonHighEtaAnalyzer::selectionParameter_

private

Definition at line 111 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyzeTracks(), and HcalHBHEMuonHighEtaAnalyzer().

◆ sevlv_

edm::ESHandle<EcalSeverityLevelAlgo> HcalHBHEMuonHighEtaAnalyzer::sevlv_

private

Definition at line 127 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyze(), and analyzeTracks().

◆ theCaloTopology_

edm::ESHandle<CaloTopology> HcalHBHEMuonHighEtaAnalyzer::theCaloTopology_

private

Definition at line 128 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyze(), and analyzeTracks().

◆ theHBHETopology_

const HcalTopology* HcalHBHEMuonHighEtaAnalyzer::theHBHETopology_

private

Definition at line 114 of file HcalHBHEMuonHighEtaAnalyzer.cc.

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

◆ theTrackQuality_

const std::string HcalHBHEMuonHighEtaAnalyzer::theTrackQuality_

private

Definition at line 107 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyzeHadron(), and HcalHBHEMuonHighEtaAnalyzer().

◆ tok_EB_

edm::EDGetTokenT<EcalRecHitCollection> HcalHBHEMuonHighEtaAnalyzer::tok_EB_

private

Definition at line 118 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyze(), and HcalHBHEMuonHighEtaAnalyzer().

◆ tok_EE_

edm::EDGetTokenT<EcalRecHitCollection> HcalHBHEMuonHighEtaAnalyzer::tok_EE_

private

Definition at line 119 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyze(), and HcalHBHEMuonHighEtaAnalyzer().

◆ tok_genTrack_

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

private

Definition at line 122 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyzeHadron(), and HcalHBHEMuonHighEtaAnalyzer().

◆ tok_HBHE_

edm::EDGetTokenT<HBHERecHitCollection> HcalHBHEMuonHighEtaAnalyzer::tok_HBHE_

private

Definition at line 120 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyze(), and HcalHBHEMuonHighEtaAnalyzer().

◆ tok_Muon_

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

private

Definition at line 121 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyzeMuon(), and HcalHBHEMuonHighEtaAnalyzer().

◆ tok_Vtx_

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

private

Definition at line 117 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyze(), and HcalHBHEMuonHighEtaAnalyzer().

◆ trackDz_

std::vector<double> HcalHBHEMuonHighEtaAnalyzer::trackDz_

private

Definition at line 158 of file HcalHBHEMuonHighEtaAnalyzer.cc.

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

◆ trackLayerCrossed_

std::vector<int> HcalHBHEMuonHighEtaAnalyzer::trackLayerCrossed_

private

Definition at line 159 of file HcalHBHEMuonHighEtaAnalyzer.cc.

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

◆ trackMissedInnerHits_

std::vector<int> HcalHBHEMuonHighEtaAnalyzer::trackMissedInnerHits_

private

Definition at line 160 of file HcalHBHEMuonHighEtaAnalyzer.cc.

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

◆ trackMissedOuterHits_

std::vector<int> HcalHBHEMuonHighEtaAnalyzer::trackMissedOuterHits_

private

Definition at line 160 of file HcalHBHEMuonHighEtaAnalyzer.cc.

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

◆ trackOuterHit_

std::vector<int> HcalHBHEMuonHighEtaAnalyzer::trackOuterHit_

private

Definition at line 159 of file HcalHBHEMuonHighEtaAnalyzer.cc.

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

◆ tree_

TTree* HcalHBHEMuonHighEtaAnalyzer::tree_

private

Definition at line 137 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyze(), and beginJob().

◆ unCorrect_

const bool HcalHBHEMuonHighEtaAnalyzer::unCorrect_

private

Definition at line 105 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyzeTracks(), and HcalHBHEMuonHighEtaAnalyzer().

◆ useMyCorr_

bool HcalHBHEMuonHighEtaAnalyzer::useMyCorr_

private

Definition at line 109 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by HcalHBHEMuonHighEtaAnalyzer(), and respCorr().

◆ useRaw_

const int HcalHBHEMuonHighEtaAnalyzer::useRaw_

private

Definition at line 106 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by analyzeTracks(), and HcalHBHEMuonHighEtaAnalyzer().

◆ verbosity_

const int HcalHBHEMuonHighEtaAnalyzer::verbosity_

private

Definition at line 106 of file HcalHBHEMuonHighEtaAnalyzer.cc.

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

◆ writeRespCorr_

const bool HcalHBHEMuonHighEtaAnalyzer::writeRespCorr_

private

Definition at line 108 of file HcalHBHEMuonHighEtaAnalyzer.cc.

Referenced by beginRun().

Public Member Functions

Static Public Member Functions

Private Member Functions

Private Attributes

Static Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ HcalHBHEMuonHighEtaAnalyzer()

Member Function Documentation

◆ activeLength()

◆ analyze()

◆ analyzeHadron()

◆ analyzeMuon()

◆ analyzeTracks()

◆ beginJob()

◆ beginRun()

◆ clearVectors()

◆ depth16HE()

◆ endRun()

◆ fillDescriptions()

◆ fillTrackParameters()

◆ gainFactor()

◆ goodCell()

◆ isGoodVertex()

◆ matchId()

◆ respCorr()

Member Data Documentation

◆ actHB

◆ actHE

◆ analyzeMuon_

◆ barrelRecHitsHandle_

◆ bFieldH_

◆ collapseDepth_

◆ conditions_

◆ corrValue_

◆ depthMax_

◆ ecal3x3Energy_

◆ ecalChStatus_

◆ ecalDetId_

◆ ecalEnergy_

◆ ehcalDetId_

◆ emaxNearP_

◆ emaxNearPThr_

◆ endcapRecHitsHandle_

◆ energyMuon_

◆ etaGlob_

◆ etaMin_

◆ eventNumber_

◆ fileInCorr_

◆ getCharge_

◆ goodVertex_

◆ hbhe_

◆ hcal1x1Energy_

◆ hcal_ieta_

◆ hcal_iphi_

◆ hcalActiveLength_

◆ hcalActiveLengthHot_

◆ hcalDepthActiveLength_

◆ hcalDepthActiveLengthHot_

◆ hcalDepthChargeHot_

◆ hcalDepthChargeHotBG_

◆ hcalDepthEnergy_

◆ hcalDepthEnergyCorr_

◆ hcalDepthEnergyHot_

◆ hcalDepthEnergyHotCorr_

◆ hcalDepthMatch_

◆ hcalDepthMatchHot_

◆ hcalDetId_

◆ hcalEnergy_

◆ hcalHot_

◆ hdc_

◆ hoEnergy_

◆ ignoreHECorr_

◆ isItPlan1_

◆ isItPreRecHit_

◆ isolationR03_

◆ isolationR04_

◆ kount_

◆ labelEBRecHit_