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

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 const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

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 const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

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

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 edm::ESHandle< HcalDbService > &, 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_

std::vector< double >	chiGlobal_

std::vector< double >	chiTracker_

const bool	collapseDepth_

std::map< DetId, double >	corrValue_

std::vector< double >	dxyTracker_

std::vector< double >	dzTracker_

std::vector< double >	ecal3x3Energy_

std::vector< unsigned int >	ecalDetId_

std::vector< double >	ecalEnergy_

std::vector< unsigned int >	ehcalDetId_

std::vector< double >	energyMuon_

std::vector< double >	etaGlob_

unsigned int	eventNumber_

const std::string	fileInCorr_

edm::Service< TFileService >	fs

const bool	getCharge_

std::vector< int >	globalMuonHits_

std::vector< bool >	globalTrack_

std::vector< double >	globalTrckEta_

std::vector< double >	globalTrckPhi_

std::vector< double >	globalTrckPt_

unsigned int	goodVertex_

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_

HLTConfigProvider	hltConfig_

std::vector< int >	hltresults_

const edm::InputTag	hlTriggerResults_

std::vector< double >	hoEnergy_

const bool	ignoreHECorr_

std::vector< bool >	innerTrack_

std::vector< double >	innerTracketa_

std::vector< double >	innerTrackphi_

std::vector< double >	innerTrackpt_

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 edm::InputTag	labelHBHERecHit_

const std::string	labelMuon_

const std::string	labelVtx_

unsigned int	lumiNumber_

std::vector< bool >	matchedId_

std::vector< int >	matchedStat_

int	maxDepth_

bool	mergedDepth_

std::vector< float >	muon_chi2LocalPosition

std::vector< bool >	muon_global_

std::vector< bool >	muon_is_good_

std::vector< bool >	muon_is_medium_

std::vector< bool >	muon_is_tight_

std::vector< float >	muon_segComp

std::vector< bool >	muon_tracker_

std::vector< float >	muon_trkKink

std::vector< int >	numPixelLayers_

std::vector< bool >	outerTrack_

std::vector< double >	outerTrackChi_

std::vector< double >	outerTrackEta_

std::vector< int >	outerTrackHits_

std::vector< double >	outerTrackPhi_

std::vector< double >	outerTrackPt_

std::vector< int >	outerTrackRHits_

std::vector< double >	phiGlob_

std::vector< double >	pMuon_

std::vector< double >	ptGlob_

HcalRespCorrs *	respCorrs_

unsigned int	runNumber_

const HcalTopology *	theHBHETopology_

std::vector< double >	tight_LongPara_

std::vector< int >	tight_PixelHits_

std::vector< double >	tight_TransImpara_

std::vector< double >	tight_validFraction_

edm::EDGetTokenT< EcalRecHitCollection >	tok_EB_

edm::EDGetTokenT< EcalRecHitCollection >	tok_EE_

edm::EDGetTokenT< HBHERecHitCollection >	tok_HBHE_

edm::EDGetTokenT< reco::MuonCollection >	tok_Muon_

edm::EDGetTokenT< edm::TriggerResults >	tok_trigRes_

edm::EDGetTokenT< reco::VertexCollection >	tok_Vtx_

std::vector< int >	trackerLayer_

TTree *	tree_

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

const bool	unCorrect_

bool	useMyCorr_

const int	useRaw_

const int	verbosity_

const bool	writeRespCorr_

Static Private Attributes
static const int	depthMax_ = 7

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

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

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

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

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

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

void	consumesMany (const TypeToGet &id)

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

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

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

template<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)

Detailed Description

Definition at line 72 of file HcalHBHEMuonAnalyzer.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 161 of file HcalHBHEMuonAnalyzer.cc.

References collapseDepth_, corrValue_, depthMax_, fileInCorr_, getCharge_, edm::ParameterSet::getUntrackedParameter(), hlTriggerResults_, triggerObjects_cff::id, ignoreHECorr_, isItPlan1_, isItPreRecHit_, kount_, TFileService::kSharedResource, labelEBRecHit_, labelEERecHit_, labelHBHERecHit_, labelMuon_, labelVtx_, maxDepth_, mergedDepth_, AlCaHLTBitMon_QueryRunRegistry::string, tok_EB_, tok_EE_, tok_HBHE_, tok_Muon_, tok_trigRes_, tok_Vtx_, unCorrect_, useMyCorr_, 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")),
       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)),
       hdc_(nullptr),
       theHBHETopology_(nullptr),
       respCorrs_(nullptr) {
   usesResource(TFileService::kSharedResource);
   //now do what ever initialization is needed
   kount_ = 0;
   maxDepth_ = iConfig.getUntrackedParameter<int>("maxDepth", 4);
   if (maxDepth_ > depthMax_)
     maxDepth_ = depthMax_;
   else if (maxDepth_ < 1)
     maxDepth_ = 4;
   std::string modnam = iConfig.getUntrackedParameter<std::string>("moduleName", "");
   std::string procnm = iConfig.getUntrackedParameter<std::string>("processName", "");
 
   mergedDepth_ = (!isItPreRecHit_) || (collapseDepth_);
   tok_trigRes_ = consumes<edm::TriggerResults>(hlTriggerResults_);
   tok_EB_ = consumes<EcalRecHitCollection>(labelEBRecHit_);
   tok_EE_ = consumes<EcalRecHitCollection>(labelEERecHit_);
   tok_HBHE_ = consumes<HBHERecHitCollection>(labelHBHERecHit_);
   if (modnam.empty()) {
     tok_Vtx_ = consumes<reco::VertexCollection>(labelVtx_);
     tok_Muon_ = consumes<reco::MuonCollection>(labelMuon_);
     edm::LogVerbatim("HBHEMuon") << "Labels used: Trig " << hlTriggerResults_ << " Vtx " << labelVtx_ << " EB "
                                  << labelEBRecHit_ << " EE " << labelEERecHit_ << " HBHE " << labelHBHERecHit_ << " MU "
                                  << labelMuon_;
   } else {
     tok_Vtx_ = consumes<reco::VertexCollection>(edm::InputTag(modnam, labelVtx_, procnm));
     tok_Muon_ = consumes<reco::MuonCollection>(edm::InputTag(modnam, labelMuon_, procnm));
     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_;
 }

Member Function Documentation

double HcalHBHEMuonAnalyzer::activeLength ( const DetId & id_ )

private

Definition at line 1086 of file HcalHBHEMuonAnalyzer.cc.

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

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

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

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

overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 239 of file HcalHBHEMuonAnalyzer.cc.

References funct::abs(), accept(), edm::HLTGlobalStatus::accept(), activeLength(), all_triggers_, ecalTB2006H4_GenSimDigiReco_cfg::bField, edm::EventBase::bunchCrossing(), bxNumber_, trackerTree::check(), chg, chiGlobal_, chiTracker_, clearVectors(), collapseDepth_, corr, egammaForCoreTracking_cff::depth, HcalDetId::depth(), depth16HE(), depthMax_, spr::propagatedTrackID::detIdECAL, spr::propagatedTrackID::detIdEHCAL, spr::propagatedTrackID::detIdHCAL, dxyTracker_, dzTracker_, ecal3x3Energy_, ecalDetId_, ecalEnergy_, spr::eECALmatrix(), ehcalDetId_, spr::eHCALmatrix(), spr::energyHCALCell(), energyMuon_, etaGlob_, edm::EventID::event(), eventNumber_, muonCSCDigis_cfi::gain, gainFactor(), edm::EventSetup::get(), edm::Event::getByToken(), getCharge_, globalMuonHits_, 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, edm::EventBase::id(), HcalDetId::ieta(), ignoreHECorr_, info(), innerTrack_, innerTracketa_, innerTrackphi_, innerTrackpt_, createfilelist::int, HcalDetId::iphi(), isGoodVertex(), isItPlan1_, isItPreRecHit_, muon::isMediumMuon(), isolationR03_, isolationR04_, muon::isTightMuon(), edm::HandleBase::isValid(), gen::k, kount_, edm::EventID::luminosityBlock(), 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_, pMuon_, edm::ESHandle< T >::product(), spr::propagateCALO(), spr::propagateHCALBack(), EnergyCorrector::pt, ptGlob_, DetId::rawId(), respCorr(), edm::EventID::run(), runNumber_, edm::second(), muon::segmentCompatibility(), edm::HLTGlobalStatus::size(), HcalDetId::subdet(), theHBHETopology_, tight_LongPara_, tight_PixelHits_, tight_TransImpara_, tight_validFraction_, tok_EB_, tok_EE_, tok_HBHE_, tok_Muon_, tok_trigRes_, tok_Vtx_, trackerLayer_, tree_, edm::TriggerNames::triggerIndex(), edm::Event::triggerNames(), unCorrect_, useRaw_, HcalTopology::validHcal(), verbosity_, and extraflags_cff::vtx.

                                                                                       {
   ++kount_;
   clearVectors();
   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
   edm::Handle<edm::TriggerResults> _Triggers;
   iEvent.getByToken(tok_trigRes_, _Triggers);
 #ifdef EDM_ML_DEBUG
   if ((verbosity_ / 10000) % 10 > 0)
     edm::LogVerbatim("HBHEMuon") << "Size of all triggers " << all_triggers_.size() << std::endl;
 #endif
   int Ntriggers = all_triggers_.size();
 #ifdef EDM_ML_DEBUG
   if ((verbosity_ / 10000) % 10 > 0)
     edm::LogVerbatim("HBHEMuon") << "Size of HLT MENU: " << _Triggers->size() << std::endl;
 #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 = int(_Triggers->size());
 #ifdef EDM_ML_DEBUG
       if ((verbosity_ / 10000) % 10 > 0)
         edm::LogVerbatim("HBHEMuon") << "outside loop " << index[i] << "\ntriggerSize " << triggerSize << std::endl;
 #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" << std::endl;
 #endif
       } else {
         if ((verbosity_ / 10000) % 10 > 0)
           edm::LogVerbatim("HBHEMuon") << "Requested HLT path \""
                                        << "\" does not exist\n";
       }
     }
   }
 
   // get handles to calogeometry and calotopology
   edm::ESHandle<CaloGeometry> pG;
   iSetup.get<CaloGeometryRecord>().get(pG);
   const CaloGeometry* geo = pG.product();
 
   edm::ESHandle<MagneticField> bFieldH;
   iSetup.get<IdealMagneticFieldRecord>().get(bFieldH);
   const MagneticField* bField = bFieldH.product();
 
   edm::ESHandle<EcalChannelStatus> ecalChStatus;
   iSetup.get<EcalChannelStatusRcd>().get(ecalChStatus);
   const EcalChannelStatus* theEcalChStatus = ecalChStatus.product();
 
   edm::ESHandle<EcalSeverityLevelAlgo> sevlv;
   iSetup.get<EcalSeverityLevelAlgoRcd>().get(sevlv);
 
   edm::ESHandle<CaloTopology> theCaloTopology;
   iSetup.get<CaloTopologyRecord>().get(theCaloTopology);
   const CaloTopology* caloTopology = theCaloTopology.product();
 
   edm::ESHandle<HcalDbService> conditions;
   iSetup.get<HcalDbRecord>().get(conditions);
 
   // Relevant blocks from iEvent
   edm::Handle<reco::VertexCollection> vtx;
   iEvent.getByToken(tok_Vtx_, vtx);
 
   edm::Handle<EcalRecHitCollection> barrelRecHitsHandle;
   iEvent.getByToken(tok_EB_, barrelRecHitsHandle);
   edm::Handle<EcalRecHitCollection> endcapRecHitsHandle;
   iEvent.getByToken(tok_EE_, endcapRecHitsHandle);
 
   edm::Handle<HBHERecHitCollection> hbhe;
   iEvent.getByToken(tok_HBHE_, hbhe);
 
   edm::Handle<reco::MuonCollection> _Muon;
   iEvent.getByToken(tok_Muon_, _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\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 (_Muon.isValid() && barrelRecHitsHandle.isValid() && endcapRecHitsHandle.isValid() && hbhe.isValid()) {
     for (reco::MuonCollection::const_iterator RecMuon = _Muon->begin(); RecMuon != _Muon->end(); ++RecMuon) {
       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() << std::endl;
 #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() > 10.0) && (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 = spr::eECALmatrix(isoCell,
                                                           barrelRecHitsHandle,
                                                           endcapRecHitsHandle,
                                                           *theEcalChStatus,
                                                           geo,
                                                           caloTopology,
                                                           sevlv.product(),
                                                           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));
             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
                   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 << "] " << 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,
                                 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));
               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
                     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
     tree_->Fill();
   }
 }

void HcalHBHEMuonAnalyzer::beginJob ( void )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 786 of file HcalHBHEMuonAnalyzer.cc.

References all_triggers_, bxNumber_, chiGlobal_, chiTracker_, dxyTracker_, dzTracker_, ecal3x3Energy_, ecalDetId_, ecalEnergy_, ehcalDetId_, energyMuon_, etaGlob_, eventNumber_, 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_, gen::k, lumiNumber_, TFileService::make(), matchedId_, matchedStat_, maxDepth_, muon_chi2LocalPosition, muon_global_, muon_is_good_, muon_is_medium_, muon_is_tight_, muon_segComp, muon_tracker_, muon_trkKink, dataset::name, numPixelLayers_, outerTrack_, outerTrackChi_, outerTrackEta_, outerTrackHits_, outerTrackPhi_, outerTrackPt_, outerTrackRHits_, phiGlob_, pMuon_, ptGlob_, runNumber_, tight_LongPara_, tight_PixelHits_, tight_TransImpara_, tight_validFraction_, trackerLayer_, and tree_.

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

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

overrideprivate

Definition at line 886 of file HcalHBHEMuonAnalyzer.cc.

References actHB, actHE, activeLength(), all_triggers_, edm::EventSetup::get(), CaloGeometry::getSubdetectorGeometry(), HcalDDDRecConstants::getThickActive(), HcalGeometry::getValidDetIds(), reco::JetExtendedAssociation::getValue(), HcalCondObjectContainer< Item >::getValues(), DetId::Hcal, HcalBarrel, HcalEndcap, hdc_, hltConfig_, lumiCalc2::hltname, HLTConfigProvider::init(), mergedDepth_, HcalDDDRecConstants::mergedDepthDetId(), edm::ESHandle< T >::product(), DetId::rawId(), respCorrs_, edm::RunBase::run(), HcalCondObjectContainerBase::setTopo(), HLTConfigProvider::size(), AlCaHLTBitMon_QueryRunRegistry::string, lumiQTWidget::t, theHBHETopology_, HLTConfigProvider::triggerName(), triggers_, and writeRespCorr_.

                                                                                    {
   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
   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
 
   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" << std::endl;
 #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 < 6; ++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() << std::endl;
   } else {
     edm::LogError("HBHEMuon") << "Error! HLT config extraction with process "
                               << "name HLT failed";
   }
 
   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();
       }
     }
   }
 }

void HcalHBHEMuonAnalyzer::clearVectors ( )

private

clearing vectots

Definition at line 996 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(), and endRun().

                                         {
   eventNumber_ = -99999;
   runNumber_ = -99999;
   lumiNumber_ = -99999;
   bxNumber_ = -99999;
   goodVertex_ = -99999;
 
   muon_is_good_.clear();
   muon_global_.clear();
   muon_tracker_.clear();
   ptGlob_.clear();
   etaGlob_.clear();
   phiGlob_.clear();
   energyMuon_.clear();
   pMuon_.clear();
   muon_trkKink.clear();
   muon_chi2LocalPosition.clear();
   muon_segComp.clear();
   muon_is_tight_.clear();
   muon_is_medium_.clear();
 
   trackerLayer_.clear();
   numPixelLayers_.clear();
   tight_PixelHits_.clear();
   innerTrack_.clear();
   chiTracker_.clear();
   dxyTracker_.clear();
   dzTracker_.clear();
   innerTrackpt_.clear();
   innerTracketa_.clear();
   innerTrackphi_.clear();
   tight_validFraction_.clear();
 
   outerTrack_.clear();
   outerTrackPt_.clear();
   outerTrackEta_.clear();
   outerTrackPhi_.clear();
   outerTrackHits_.clear();
   outerTrackRHits_.clear();
   outerTrackChi_.clear();
 
   globalTrack_.clear();
   globalTrckPt_.clear();
   globalTrckEta_.clear();
   globalTrckPhi_.clear();
   globalMuonHits_.clear();
   matchedStat_.clear();
   chiGlobal_.clear();
   tight_LongPara_.clear();
   tight_TransImpara_.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 < maxDepth_; ++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();
   hltresults_.clear();
 }

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

private

Definition at line 1153 of file HcalHBHEMuonAnalyzer.cc.

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

Referenced by analyze(), and endRun().

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

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

inlineoverrideprivate

Definition at line 82 of file HcalHBHEMuonAnalyzer.cc.

References activeLength(), ecalTB2006H4_GenSimDigiReco_cfg::bField, clearVectors(), depth16HE(), gainFactor(), goodCell(), isGoodVertex(), matchId(), respCorr(), and extraflags_cff::vtx.

82 {}

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

static

Definition at line 970 of file HcalHBHEMuonAnalyzer.cc.

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

                                                                                       {
   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 = {"HLT_IsoMu17", "HLT_IsoMu20", "HLT_IsoMu24", "HLT_IsoMu27", "HLT_Mu45", "HLT_Mu50"};
   desc.add<std::vector<std::string>>("triggers", trig);
   desc.addUntracked<int>("verbosity", 0);
   desc.add<int>("useRaw", 0);
   desc.add<bool>("unCorrect", false);
   desc.add<bool>("getCharge", false);
   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", 4);
   desc.addUntracked<std::string>("fileInCorr", "");
   desc.addUntracked<bool>("writeRespCorr", false);
   descriptions.add("hcalHBHEMuon", desc);
 }

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

private

Definition at line 1145 of file HcalHBHEMuonAnalyzer.cc.

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

Referenced by analyze(), and endRun().

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

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

private

Definition at line 1168 of file HcalHBHEMuonAnalyzer.cc.

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

Referenced by analyze(), and endRun().

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

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

private

Definition at line 1121 of file HcalHBHEMuonAnalyzer.cc.

References reco::Vertex::isFake(), reco::Vertex::ndof(), and reco::Vertex::position().

Referenced by analyze(), and endRun().

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

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

private

Definition at line 1079 of file HcalHBHEMuonAnalyzer.cc.

References HcalDetId::ieta(), HcalDetId::iphi(), match(), and HcalDetId::subdet().

Referenced by analyze(), and endRun().

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

double HcalHBHEMuonAnalyzer::respCorr ( const DetId & id )

private

Definition at line 1133 of file HcalHBHEMuonAnalyzer.cc.

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

Referenced by analyze(), and endRun().

                                                      {
   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

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

private

Definition at line 156 of file HcalHBHEMuonAnalyzer.cc.

Referenced by activeLength(), and beginRun().

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

private

Definition at line 156 of file HcalHBHEMuonAnalyzer.cc.

Referenced by activeLength(), and beginRun().

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

private

Definition at line 153 of file HcalHBHEMuonAnalyzer.cc.

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

unsigned int HcalHBHEMuonAnalyzer::bxNumber_

private

Definition at line 120 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<double> HcalHBHEMuonAnalyzer::chiGlobal_

private

Definition at line 135 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<double> HcalHBHEMuonAnalyzer::chiTracker_

private

Definition at line 128 of file HcalHBHEMuonAnalyzer.cc.

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

const bool HcalHBHEMuonAnalyzer::collapseDepth_

private

Definition at line 100 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze(), and HcalHBHEMuonAnalyzer().

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

private

Definition at line 157 of file HcalHBHEMuonAnalyzer.cc.

Referenced by HcalHBHEMuonAnalyzer(), and respCorr().

const int HcalHBHEMuonAnalyzer::depthMax_ = 7

staticprivate

Definition at line 118 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze(), and HcalHBHEMuonAnalyzer().

std::vector<double> HcalHBHEMuonAnalyzer::dxyTracker_

private

Definition at line 128 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<double> HcalHBHEMuonAnalyzer::dzTracker_

private

Definition at line 128 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<double> HcalHBHEMuonAnalyzer::ecal3x3Energy_

private

Definition at line 139 of file HcalHBHEMuonAnalyzer.cc.

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

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

private

Definition at line 140 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<double> HcalHBHEMuonAnalyzer::ecalEnergy_

private

Definition at line 137 of file HcalHBHEMuonAnalyzer.cc.

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

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

private

Definition at line 140 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<double> HcalHBHEMuonAnalyzer::energyMuon_

private

Definition at line 124 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<double> HcalHBHEMuonAnalyzer::etaGlob_

private

Definition at line 124 of file HcalHBHEMuonAnalyzer.cc.

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

unsigned int HcalHBHEMuonAnalyzer::eventNumber_

private

Definition at line 120 of file HcalHBHEMuonAnalyzer.cc.

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

const std::string HcalHBHEMuonAnalyzer::fileInCorr_

private

Definition at line 97 of file HcalHBHEMuonAnalyzer.cc.

Referenced by HcalHBHEMuonAnalyzer().

edm::Service<TFileService> HcalHBHEMuonAnalyzer::fs

private

Definition at line 94 of file HcalHBHEMuonAnalyzer.cc.

Referenced by beginJob().

const bool HcalHBHEMuonAnalyzer::getCharge_

private

Definition at line 102 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze(), and HcalHBHEMuonAnalyzer().

std::vector<int> HcalHBHEMuonAnalyzer::globalMuonHits_

private

Definition at line 134 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<bool> HcalHBHEMuonAnalyzer::globalTrack_

private

Definition at line 127 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<double> HcalHBHEMuonAnalyzer::globalTrckEta_

private

Definition at line 133 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<double> HcalHBHEMuonAnalyzer::globalTrckPhi_

private

Definition at line 133 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<double> HcalHBHEMuonAnalyzer::globalTrckPt_

private

Definition at line 133 of file HcalHBHEMuonAnalyzer.cc.

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

unsigned int HcalHBHEMuonAnalyzer::goodVertex_

private

Definition at line 121 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<double> HcalHBHEMuonAnalyzer::hcal1x1Energy_

private

Definition at line 139 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<int> HcalHBHEMuonAnalyzer::hcal_ieta_

private

Definition at line 141 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<int> HcalHBHEMuonAnalyzer::hcal_iphi_

private

Definition at line 141 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<double> HcalHBHEMuonAnalyzer::hcalActiveLength_

private

Definition at line 152 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<double> HcalHBHEMuonAnalyzer::hcalActiveLengthHot_

private

Definition at line 152 of file HcalHBHEMuonAnalyzer.cc.

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

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

private

Definition at line 143 of file HcalHBHEMuonAnalyzer.cc.

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

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

private

Definition at line 145 of file HcalHBHEMuonAnalyzer.cc.

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

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

private

Definition at line 146 of file HcalHBHEMuonAnalyzer.cc.

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

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

private

Definition at line 147 of file HcalHBHEMuonAnalyzer.cc.

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

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

private

Definition at line 142 of file HcalHBHEMuonAnalyzer.cc.

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

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

private

Definition at line 148 of file HcalHBHEMuonAnalyzer.cc.

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

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

private

Definition at line 144 of file HcalHBHEMuonAnalyzer.cc.

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

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

private

Definition at line 149 of file HcalHBHEMuonAnalyzer.cc.

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

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

private

Definition at line 150 of file HcalHBHEMuonAnalyzer.cc.

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

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

private

Definition at line 151 of file HcalHBHEMuonAnalyzer.cc.

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

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

private

Definition at line 140 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<double> HcalHBHEMuonAnalyzer::hcalEnergy_

private

Definition at line 137 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<bool> HcalHBHEMuonAnalyzer::hcalHot_

private

Definition at line 138 of file HcalHBHEMuonAnalyzer.cc.

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

const HcalDDDRecConstants* HcalHBHEMuonAnalyzer::hdc_

private

Definition at line 106 of file HcalHBHEMuonAnalyzer.cc.

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

HLTConfigProvider HcalHBHEMuonAnalyzer::hltConfig_

private

Definition at line 93 of file HcalHBHEMuonAnalyzer.cc.

Referenced by beginRun().

std::vector<int> HcalHBHEMuonAnalyzer::hltresults_

private

Definition at line 154 of file HcalHBHEMuonAnalyzer.cc.

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

const edm::InputTag HcalHBHEMuonAnalyzer::hlTriggerResults_

private

Definition at line 95 of file HcalHBHEMuonAnalyzer.cc.

Referenced by HcalHBHEMuonAnalyzer().

std::vector<double> HcalHBHEMuonAnalyzer::hoEnergy_

private

Definition at line 137 of file HcalHBHEMuonAnalyzer.cc.

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

const bool HcalHBHEMuonAnalyzer::ignoreHECorr_

private

Definition at line 101 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze(), and HcalHBHEMuonAnalyzer().

std::vector<bool> HcalHBHEMuonAnalyzer::innerTrack_

private

Definition at line 127 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<double> HcalHBHEMuonAnalyzer::innerTracketa_

private

Definition at line 129 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<double> HcalHBHEMuonAnalyzer::innerTrackphi_

private

Definition at line 129 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<double> HcalHBHEMuonAnalyzer::innerTrackpt_

private

Definition at line 129 of file HcalHBHEMuonAnalyzer.cc.

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

const bool HcalHBHEMuonAnalyzer::isItPlan1_

private

Definition at line 100 of file HcalHBHEMuonAnalyzer.cc.

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

const bool HcalHBHEMuonAnalyzer::isItPreRecHit_

private

Definition at line 101 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<double> HcalHBHEMuonAnalyzer::isolationR03_

private

Definition at line 136 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<double> HcalHBHEMuonAnalyzer::isolationR04_

private

Definition at line 136 of file HcalHBHEMuonAnalyzer.cc.

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

int HcalHBHEMuonAnalyzer::kount_

private

Definition at line 104 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze(), and HcalHBHEMuonAnalyzer().

const edm::InputTag HcalHBHEMuonAnalyzer::labelEBRecHit_

private

Definition at line 96 of file HcalHBHEMuonAnalyzer.cc.

Referenced by HcalHBHEMuonAnalyzer().

const edm::InputTag HcalHBHEMuonAnalyzer::labelEERecHit_

private

Definition at line 96 of file HcalHBHEMuonAnalyzer.cc.

Referenced by HcalHBHEMuonAnalyzer().

const edm::InputTag HcalHBHEMuonAnalyzer::labelHBHERecHit_

private

Definition at line 96 of file HcalHBHEMuonAnalyzer.cc.

Referenced by HcalHBHEMuonAnalyzer().

const std::string HcalHBHEMuonAnalyzer::labelMuon_

private

Definition at line 97 of file HcalHBHEMuonAnalyzer.cc.

Referenced by HcalHBHEMuonAnalyzer().

const std::string HcalHBHEMuonAnalyzer::labelVtx_

private

Definition at line 97 of file HcalHBHEMuonAnalyzer.cc.

Referenced by HcalHBHEMuonAnalyzer().

unsigned int HcalHBHEMuonAnalyzer::lumiNumber_

private

Definition at line 120 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<bool> HcalHBHEMuonAnalyzer::matchedId_

private

Definition at line 138 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<int> HcalHBHEMuonAnalyzer::matchedStat_

private

Definition at line 134 of file HcalHBHEMuonAnalyzer.cc.

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

int HcalHBHEMuonAnalyzer::maxDepth_

private

Definition at line 104 of file HcalHBHEMuonAnalyzer.cc.

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

bool HcalHBHEMuonAnalyzer::mergedDepth_

private

Definition at line 103 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<float> HcalHBHEMuonAnalyzer::muon_chi2LocalPosition

private

Definition at line 125 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<bool> HcalHBHEMuonAnalyzer::muon_global_

private

Definition at line 122 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<bool> HcalHBHEMuonAnalyzer::muon_is_good_

private

Definition at line 122 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<bool> HcalHBHEMuonAnalyzer::muon_is_medium_

private

Definition at line 123 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<bool> HcalHBHEMuonAnalyzer::muon_is_tight_

private

Definition at line 123 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<float> HcalHBHEMuonAnalyzer::muon_segComp

private

Definition at line 125 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<bool> HcalHBHEMuonAnalyzer::muon_tracker_

private

Definition at line 122 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<float> HcalHBHEMuonAnalyzer::muon_trkKink

private

Definition at line 125 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<int> HcalHBHEMuonAnalyzer::numPixelLayers_

private

Definition at line 126 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<bool> HcalHBHEMuonAnalyzer::outerTrack_

private

Definition at line 127 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<double> HcalHBHEMuonAnalyzer::outerTrackChi_

private

Definition at line 130 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<double> HcalHBHEMuonAnalyzer::outerTrackEta_

private

Definition at line 131 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<int> HcalHBHEMuonAnalyzer::outerTrackHits_

private

Definition at line 132 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<double> HcalHBHEMuonAnalyzer::outerTrackPhi_

private

Definition at line 131 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<double> HcalHBHEMuonAnalyzer::outerTrackPt_

private

Definition at line 131 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<int> HcalHBHEMuonAnalyzer::outerTrackRHits_

private

Definition at line 132 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<double> HcalHBHEMuonAnalyzer::phiGlob_

private

Definition at line 124 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<double> HcalHBHEMuonAnalyzer::pMuon_

private

Definition at line 124 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<double> HcalHBHEMuonAnalyzer::ptGlob_

private

Definition at line 124 of file HcalHBHEMuonAnalyzer.cc.

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

HcalRespCorrs* HcalHBHEMuonAnalyzer::respCorrs_

private

Definition at line 108 of file HcalHBHEMuonAnalyzer.cc.

Referenced by beginRun(), and respCorr().

unsigned int HcalHBHEMuonAnalyzer::runNumber_

private

Definition at line 120 of file HcalHBHEMuonAnalyzer.cc.

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

const HcalTopology* HcalHBHEMuonAnalyzer::theHBHETopology_

private

Definition at line 107 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<double> HcalHBHEMuonAnalyzer::tight_LongPara_

private

Definition at line 135 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<int> HcalHBHEMuonAnalyzer::tight_PixelHits_

private

Definition at line 126 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<double> HcalHBHEMuonAnalyzer::tight_TransImpara_

private

Definition at line 135 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<double> HcalHBHEMuonAnalyzer::tight_validFraction_

private

Definition at line 130 of file HcalHBHEMuonAnalyzer.cc.

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

edm::EDGetTokenT<EcalRecHitCollection> HcalHBHEMuonAnalyzer::tok_EB_

private

Definition at line 112 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze(), and HcalHBHEMuonAnalyzer().

edm::EDGetTokenT<EcalRecHitCollection> HcalHBHEMuonAnalyzer::tok_EE_

private

Definition at line 113 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze(), and HcalHBHEMuonAnalyzer().

edm::EDGetTokenT<HBHERecHitCollection> HcalHBHEMuonAnalyzer::tok_HBHE_

private

Definition at line 114 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze(), and HcalHBHEMuonAnalyzer().

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

private

Definition at line 115 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze(), and HcalHBHEMuonAnalyzer().

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

private

Definition at line 110 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze(), and HcalHBHEMuonAnalyzer().

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

private

Definition at line 111 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze(), and HcalHBHEMuonAnalyzer().

std::vector<int> HcalHBHEMuonAnalyzer::trackerLayer_

private

Definition at line 126 of file HcalHBHEMuonAnalyzer.cc.

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

TTree* HcalHBHEMuonAnalyzer::tree_

private

Definition at line 119 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze(), and beginJob().

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

private

Definition at line 98 of file HcalHBHEMuonAnalyzer.cc.

Referenced by beginRun().

const bool HcalHBHEMuonAnalyzer::unCorrect_

private

Definition at line 100 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze(), and HcalHBHEMuonAnalyzer().

bool HcalHBHEMuonAnalyzer::useMyCorr_

private

Definition at line 103 of file HcalHBHEMuonAnalyzer.cc.

Referenced by HcalHBHEMuonAnalyzer(), and respCorr().

const int HcalHBHEMuonAnalyzer::useRaw_

private

Definition at line 99 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze(), and HcalHBHEMuonAnalyzer().

const int HcalHBHEMuonAnalyzer::verbosity_

private

Definition at line 99 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze(), and goodCell().

const bool HcalHBHEMuonAnalyzer::writeRespCorr_

private

Definition at line 102 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

Member Function Documentation

Member Data Documentation