HcalHBHEMuonAnalyzer Class Reference

Inheritance diagram for HcalHBHEMuonAnalyzer:

Public Member Functions
	HcalHBHEMuonAnalyzer (const edm::ParameterSet &)
Public Member Functions inherited from edm::one::EDAnalyzer< edm::one::WatchRuns, edm::one::SharedResources >
	EDAnalyzer ()=default
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
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)
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 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 164 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 1063 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 246 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 badGlobalMuonTaggersAOD_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 749 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 850 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 972 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 1127 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 83 of file HcalHBHEMuonAnalyzer.cc.

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

{}

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

static

Definition at line 944 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 1119 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 1142 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 1100 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 1055 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 1108 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 158 of file HcalHBHEMuonAnalyzer.cc.

Referenced by activeLength(), and beginRun().

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

private

Definition at line 158 of file HcalHBHEMuonAnalyzer.cc.

Referenced by activeLength(), and beginRun().

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

private

Definition at line 155 of file HcalHBHEMuonAnalyzer.cc.

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

unsigned int HcalHBHEMuonAnalyzer::bxNumber_

private

Definition at line 122 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<double> HcalHBHEMuonAnalyzer::chiGlobal_

private

Definition at line 137 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<double> HcalHBHEMuonAnalyzer::chiTracker_

private

Definition at line 130 of file HcalHBHEMuonAnalyzer.cc.

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

const bool HcalHBHEMuonAnalyzer::collapseDepth_

private

Definition at line 102 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze(), and HcalHBHEMuonAnalyzer().

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

private

Definition at line 159 of file HcalHBHEMuonAnalyzer.cc.

Referenced by HcalHBHEMuonAnalyzer(), and respCorr().

const int HcalHBHEMuonAnalyzer::depthMax_ = 7

staticprivate

Definition at line 120 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze(), and HcalHBHEMuonAnalyzer().

std::vector<double> HcalHBHEMuonAnalyzer::dxyTracker_

private

Definition at line 130 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<double> HcalHBHEMuonAnalyzer::dzTracker_

private

Definition at line 130 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<double> HcalHBHEMuonAnalyzer::ecal3x3Energy_

private

Definition at line 141 of file HcalHBHEMuonAnalyzer.cc.

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

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

private

Definition at line 142 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<double> HcalHBHEMuonAnalyzer::ecalEnergy_

private

Definition at line 139 of file HcalHBHEMuonAnalyzer.cc.

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

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

private

Definition at line 142 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<double> HcalHBHEMuonAnalyzer::energyMuon_

private

Definition at line 126 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<double> HcalHBHEMuonAnalyzer::etaGlob_

private

Definition at line 126 of file HcalHBHEMuonAnalyzer.cc.

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

unsigned int HcalHBHEMuonAnalyzer::eventNumber_

private

Definition at line 122 of file HcalHBHEMuonAnalyzer.cc.

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

const std::string HcalHBHEMuonAnalyzer::fileInCorr_

private

Definition at line 99 of file HcalHBHEMuonAnalyzer.cc.

Referenced by HcalHBHEMuonAnalyzer().

edm::Service<TFileService> HcalHBHEMuonAnalyzer::fs

private

Definition at line 96 of file HcalHBHEMuonAnalyzer.cc.

Referenced by beginJob().

const bool HcalHBHEMuonAnalyzer::getCharge_

private

Definition at line 104 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze(), and HcalHBHEMuonAnalyzer().

std::vector<int> HcalHBHEMuonAnalyzer::globalMuonHits_

private

Definition at line 136 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<bool> HcalHBHEMuonAnalyzer::globalTrack_

private

Definition at line 129 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<double> HcalHBHEMuonAnalyzer::globalTrckEta_

private

Definition at line 135 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<double> HcalHBHEMuonAnalyzer::globalTrckPhi_

private

Definition at line 135 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<double> HcalHBHEMuonAnalyzer::globalTrckPt_

private

Definition at line 135 of file HcalHBHEMuonAnalyzer.cc.

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

unsigned int HcalHBHEMuonAnalyzer::goodVertex_

private

Definition at line 123 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<double> HcalHBHEMuonAnalyzer::hcal1x1Energy_

private

Definition at line 141 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<int> HcalHBHEMuonAnalyzer::hcal_ieta_

private

Definition at line 143 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<int> HcalHBHEMuonAnalyzer::hcal_iphi_

private

Definition at line 143 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<double> HcalHBHEMuonAnalyzer::hcalActiveLength_

private

Definition at line 154 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<double> HcalHBHEMuonAnalyzer::hcalActiveLengthHot_

private

Definition at line 154 of file HcalHBHEMuonAnalyzer.cc.

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

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

private

Definition at line 145 of file HcalHBHEMuonAnalyzer.cc.

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

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

private

Definition at line 147 of file HcalHBHEMuonAnalyzer.cc.

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

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

private

Definition at line 148 of file HcalHBHEMuonAnalyzer.cc.

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

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

private

Definition at line 149 of file HcalHBHEMuonAnalyzer.cc.

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

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

private

Definition at line 144 of file HcalHBHEMuonAnalyzer.cc.

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

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

private

Definition at line 150 of file HcalHBHEMuonAnalyzer.cc.

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

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

private

Definition at line 146 of file HcalHBHEMuonAnalyzer.cc.

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

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

private

Definition at line 151 of file HcalHBHEMuonAnalyzer.cc.

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

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

private

Definition at line 152 of file HcalHBHEMuonAnalyzer.cc.

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

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

private

Definition at line 153 of file HcalHBHEMuonAnalyzer.cc.

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

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

private

Definition at line 142 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<double> HcalHBHEMuonAnalyzer::hcalEnergy_

private

Definition at line 139 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<bool> HcalHBHEMuonAnalyzer::hcalHot_

private

Definition at line 140 of file HcalHBHEMuonAnalyzer.cc.

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

const HcalDDDRecConstants* HcalHBHEMuonAnalyzer::hdc_

private

Definition at line 108 of file HcalHBHEMuonAnalyzer.cc.

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

HLTConfigProvider HcalHBHEMuonAnalyzer::hltConfig_

private

Definition at line 95 of file HcalHBHEMuonAnalyzer.cc.

Referenced by beginRun().

std::vector<int> HcalHBHEMuonAnalyzer::hltresults_

private

Definition at line 156 of file HcalHBHEMuonAnalyzer.cc.

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

const edm::InputTag HcalHBHEMuonAnalyzer::hlTriggerResults_

private

Definition at line 97 of file HcalHBHEMuonAnalyzer.cc.

Referenced by HcalHBHEMuonAnalyzer().

std::vector<double> HcalHBHEMuonAnalyzer::hoEnergy_

private

Definition at line 139 of file HcalHBHEMuonAnalyzer.cc.

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

const bool HcalHBHEMuonAnalyzer::ignoreHECorr_

private

Definition at line 103 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze(), and HcalHBHEMuonAnalyzer().

std::vector<bool> HcalHBHEMuonAnalyzer::innerTrack_

private

Definition at line 129 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<double> HcalHBHEMuonAnalyzer::innerTracketa_

private

Definition at line 131 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<double> HcalHBHEMuonAnalyzer::innerTrackphi_

private

Definition at line 131 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<double> HcalHBHEMuonAnalyzer::innerTrackpt_

private

Definition at line 131 of file HcalHBHEMuonAnalyzer.cc.

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

const bool HcalHBHEMuonAnalyzer::isItPlan1_

private

Definition at line 102 of file HcalHBHEMuonAnalyzer.cc.

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

const bool HcalHBHEMuonAnalyzer::isItPreRecHit_

private

Definition at line 103 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<double> HcalHBHEMuonAnalyzer::isolationR03_

private

Definition at line 138 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<double> HcalHBHEMuonAnalyzer::isolationR04_

private

Definition at line 138 of file HcalHBHEMuonAnalyzer.cc.

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

int HcalHBHEMuonAnalyzer::kount_

private

Definition at line 106 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze(), and HcalHBHEMuonAnalyzer().

const edm::InputTag HcalHBHEMuonAnalyzer::labelEBRecHit_

private

Definition at line 98 of file HcalHBHEMuonAnalyzer.cc.

Referenced by HcalHBHEMuonAnalyzer().

const edm::InputTag HcalHBHEMuonAnalyzer::labelEERecHit_

private

Definition at line 98 of file HcalHBHEMuonAnalyzer.cc.

Referenced by HcalHBHEMuonAnalyzer().

const edm::InputTag HcalHBHEMuonAnalyzer::labelHBHERecHit_

private

Definition at line 98 of file HcalHBHEMuonAnalyzer.cc.

Referenced by HcalHBHEMuonAnalyzer().

const std::string HcalHBHEMuonAnalyzer::labelMuon_

private

Definition at line 99 of file HcalHBHEMuonAnalyzer.cc.

Referenced by HcalHBHEMuonAnalyzer().

const std::string HcalHBHEMuonAnalyzer::labelVtx_

private

Definition at line 99 of file HcalHBHEMuonAnalyzer.cc.

Referenced by HcalHBHEMuonAnalyzer().

unsigned int HcalHBHEMuonAnalyzer::lumiNumber_

private

Definition at line 122 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<bool> HcalHBHEMuonAnalyzer::matchedId_

private

Definition at line 140 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<int> HcalHBHEMuonAnalyzer::matchedStat_

private

Definition at line 136 of file HcalHBHEMuonAnalyzer.cc.

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

int HcalHBHEMuonAnalyzer::maxDepth_

private

Definition at line 106 of file HcalHBHEMuonAnalyzer.cc.

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

bool HcalHBHEMuonAnalyzer::mergedDepth_

private

Definition at line 105 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<float> HcalHBHEMuonAnalyzer::muon_chi2LocalPosition

private

Definition at line 127 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<bool> HcalHBHEMuonAnalyzer::muon_global_

private

Definition at line 124 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<bool> HcalHBHEMuonAnalyzer::muon_is_good_

private

Definition at line 124 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<bool> HcalHBHEMuonAnalyzer::muon_is_medium_

private

Definition at line 125 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<bool> HcalHBHEMuonAnalyzer::muon_is_tight_

private

Definition at line 125 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<float> HcalHBHEMuonAnalyzer::muon_segComp

private

Definition at line 127 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<bool> HcalHBHEMuonAnalyzer::muon_tracker_

private

Definition at line 124 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<float> HcalHBHEMuonAnalyzer::muon_trkKink

private

Definition at line 127 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<int> HcalHBHEMuonAnalyzer::numPixelLayers_

private

Definition at line 128 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<bool> HcalHBHEMuonAnalyzer::outerTrack_

private

Definition at line 129 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<double> HcalHBHEMuonAnalyzer::outerTrackChi_

private

Definition at line 132 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<double> HcalHBHEMuonAnalyzer::outerTrackEta_

private

Definition at line 133 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<int> HcalHBHEMuonAnalyzer::outerTrackHits_

private

Definition at line 134 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<double> HcalHBHEMuonAnalyzer::outerTrackPhi_

private

Definition at line 133 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<double> HcalHBHEMuonAnalyzer::outerTrackPt_

private

Definition at line 133 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<int> HcalHBHEMuonAnalyzer::outerTrackRHits_

private

Definition at line 134 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<double> HcalHBHEMuonAnalyzer::phiGlob_

private

Definition at line 126 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<double> HcalHBHEMuonAnalyzer::pMuon_

private

Definition at line 126 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<double> HcalHBHEMuonAnalyzer::ptGlob_

private

Definition at line 126 of file HcalHBHEMuonAnalyzer.cc.

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

HcalRespCorrs* HcalHBHEMuonAnalyzer::respCorrs_

private

Definition at line 110 of file HcalHBHEMuonAnalyzer.cc.

Referenced by beginRun(), and respCorr().

unsigned int HcalHBHEMuonAnalyzer::runNumber_

private

Definition at line 122 of file HcalHBHEMuonAnalyzer.cc.

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

const HcalTopology* HcalHBHEMuonAnalyzer::theHBHETopology_

private

Definition at line 109 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<double> HcalHBHEMuonAnalyzer::tight_LongPara_

private

Definition at line 137 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<int> HcalHBHEMuonAnalyzer::tight_PixelHits_

private

Definition at line 128 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<double> HcalHBHEMuonAnalyzer::tight_TransImpara_

private

Definition at line 137 of file HcalHBHEMuonAnalyzer.cc.

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

std::vector<double> HcalHBHEMuonAnalyzer::tight_validFraction_

private

Definition at line 132 of file HcalHBHEMuonAnalyzer.cc.

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

edm::EDGetTokenT<EcalRecHitCollection> HcalHBHEMuonAnalyzer::tok_EB_

private

Definition at line 114 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze(), and HcalHBHEMuonAnalyzer().

edm::EDGetTokenT<EcalRecHitCollection> HcalHBHEMuonAnalyzer::tok_EE_

private

Definition at line 115 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze(), and HcalHBHEMuonAnalyzer().

edm::EDGetTokenT<HBHERecHitCollection> HcalHBHEMuonAnalyzer::tok_HBHE_

private

Definition at line 116 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze(), and HcalHBHEMuonAnalyzer().

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

private

Definition at line 117 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze(), and HcalHBHEMuonAnalyzer().

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

private

Definition at line 112 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze(), and HcalHBHEMuonAnalyzer().

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

private

Definition at line 113 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze(), and HcalHBHEMuonAnalyzer().

std::vector<int> HcalHBHEMuonAnalyzer::trackerLayer_

private

Definition at line 128 of file HcalHBHEMuonAnalyzer.cc.

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

TTree* HcalHBHEMuonAnalyzer::tree_

private

Definition at line 121 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze(), and beginJob().

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

private

Definition at line 100 of file HcalHBHEMuonAnalyzer.cc.

Referenced by beginRun().

const bool HcalHBHEMuonAnalyzer::unCorrect_

private

Definition at line 102 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze(), and HcalHBHEMuonAnalyzer().

bool HcalHBHEMuonAnalyzer::useMyCorr_

private

Definition at line 105 of file HcalHBHEMuonAnalyzer.cc.

Referenced by HcalHBHEMuonAnalyzer(), and respCorr().

const int HcalHBHEMuonAnalyzer::useRaw_

private

Definition at line 101 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze(), and HcalHBHEMuonAnalyzer().

const int HcalHBHEMuonAnalyzer::verbosity_

private

Definition at line 101 of file HcalHBHEMuonAnalyzer.cc.

Referenced by analyze(), and goodCell().

const bool HcalHBHEMuonAnalyzer::writeRespCorr_

private

Definition at line 104 of file HcalHBHEMuonAnalyzer.cc.

Referenced by beginRun().