Inheritance diagram for HcalIsoTrkAnalyzer:

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

	~HcalIsoTrkAnalyzer () override

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

	EDAnalyzer (const EDAnalyzer &)=delete

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

const EDAnalyzer &	operator= (const EDAnalyzer &)=delete

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

bool	wantsInputProcessBlocks () const final

bool	wantsProcessBlocks () const final

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

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDAnalyzerBase () override

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

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

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

std::vector< ESProxyIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const

std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (edm::Transition iTrans) const

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector < ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::array< std::vector< ModuleDescription const * > , NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const > const &labelsToDesc, std::string const &processName) const

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

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

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

void	updateLookup (eventsetup::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
void	analyze (edm::Event const &, edm::EventSetup const &) override

void	beginJob () override

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

double	dR (math::XYZTLorentzVector &, math::XYZTLorentzVector &)

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

double	eThreshold (const DetId &id, const CaloGeometry *geo) const

std::array< int, 3 >	fillTree (std::vector< math::XYZTLorentzVector > &vecL1, std::vector< math::XYZTLorentzVector > &vecL3, math::XYZPoint &leadPV, std::vector< spr::propagatedTrackDirection > &trkCaloDirections, std::vector< spr::propagatedTrackID > &trkCaloDets, const CaloGeometry geo, const CaloTopology topo, const HcalTopology theHBHETopology, const EcalChannelStatus theEcalChStatus, const EcalSeverityLevelAlgo theEcalSevlv, edm::Handle< EcalRecHitCollection > &barrelRecHitsHandle, edm::Handle< EcalRecHitCollection > &endcapRecHitsHandle, edm::Handle< HBHERecHitCollection > &hbhe, edm::Handle< CaloTowerCollection > &towerHandle, edm::Handle< reco::GenParticleCollection > &genParticles, const HcalRespCorrs respCorrs, const edm::Handle< reco::MuonCollection > &muonh)

DetId	newId (const DetId &)

bool	notaMuon (const reco::Track *pTrack0, const edm::Handle< reco::MuonCollection > &muonh)

double	rhoh (const edm::Handle< CaloTowerCollection > &)

void	storeEnergy (int indx, const HcalRespCorrs respCorrs, const std::vector< DetId > &ids, std::vector< double > &edet, double &eHcal, std::vector< unsigned int > detIds, std::vector< double > *hitEnergies)

double	trackP (const reco::Track *, const edm::Handle< reco::GenParticleCollection > &)

Private Attributes
double	a_charIsoR_

double	a_coneR1_

double	a_coneR2_

const double	a_coneR_

const double	a_mipR2_

const double	a_mipR3_

const double	a_mipR4_

const double	a_mipR5_

const double	a_mipR_

const bool	collapseDepth_

const int	dataType_

const std::vector< int >	debEvents_

bool	debug_

const double	eEcalMax_

const double	eIsolate1_

const double	eIsolate2_

std::vector< double >	etabins_

double	etadist_

double	etahalfdist_

edm::Service< TFileService >	fs

const double	hcalScale_

const HcalDDDRecConstants *	hdc_

const bool	hep17_

const double	hitEthrEB_

const double	hitEthrEE0_

const double	hitEthrEE1_

const double	hitEthrEE2_

const double	hitEthrEE3_

const double	hitEthrEEHi_

const double	hitEthrEELo_

HLTConfigProvider	hltConfig_

const bool	ignoreTrigger_

const std::string	l1Filter_

l1t::L1TGlobalUtil *	l1GtUtils_

const std::string	l1TrigName_

const std::string	l2Filter_

const std::string	l3Filter_

const std::string	labelEB_

const std::string	labelEE_

const std::string	labelGenTrack_

const std::string	labelHBHE_

const std::string	labelRecVtx_

const std::string	labelTower_

const double	maxRestrictionP_

const int	mode_

const std::vector< int >	newDepth_

unsigned int	nHigh_

unsigned int	nLow_

unsigned int	nRun_

std::vector< int >	oldDepth_

std::vector< int >	oldDet_

std::vector< int >	oldEta_

const std::vector< int >	oldID_

std::vector< double >	phibins_

double	phidist_

double	phihalfdist_

const int	prescaleHigh_

const int	prescaleLow_

const std::string	processName_

const double	pTrackHigh_

const double	pTrackLow_

const double	pTrackMin_

spr::trackSelectionParameters	selectionParameter_

const double	slopeRestrictionP_

int	t_allvertex

int	t_DataType

std::vector< unsigned int > *	t_DetIds

std::vector< unsigned int > *	t_DetIds1

std::vector< unsigned int > *	t_DetIds3

double	t_eAnnular

double	t_eHcal

double	t_eHcal10

double	t_eHcal30

double	t_emaxNearP

double	t_eMipDR

double	t_eMipDR2

double	t_eMipDR3

double	t_eMipDR4

double	t_eMipDR5

int	t_Event

unsigned int	t_EventNo

double	t_EventWeight

double	t_gentrackP

int	t_goodPV

double	t_hAnnular

std::vector< double > *	t_HitEnergies

std::vector< double > *	t_HitEnergies1

std::vector< double > *	t_HitEnergies3

std::vector< bool > *	t_hltbits

double	t_hmaxNearP

int	t_ieta

std::vector< int > *	t_ietaAll

std::vector< int > *	t_ietaGood

int	t_iphi

bool	t_L1Bit

double	t_l1eta

double	t_l1phi

double	t_l1pt

double	t_l3eta

double	t_l3phi

double	t_l3pt

double	t_mindR1

double	t_mindR2

int	t_nTrk

int	t_nVtx

double	t_p

double	t_phi

double	t_pt

bool	t_qltyFlag

bool	t_qltyMissFlag

bool	t_qltyPVFlag

double	t_rhoh

int	t_Run

unsigned int	t_RunNo

bool	t_selectTk

int	t_Tracks

int	t_TracksLoose

int	t_TracksProp

int	t_TracksSaved

int	t_TracksTight

std::vector< int > *	t_trackType

std::vector< bool > *	t_trgbits

bool	t_TrigPass

bool	t_TrigPassSel

const std::string	theTrackQuality_

const edm::InputTag	theTriggerResultsLabel_

edm::EDGetTokenT< BXVector < GlobalAlgBlk > >	tok_alg_

edm::ESGetToken< MagneticField, IdealMagneticFieldRecord >	tok_bFieldH_

edm::EDGetTokenT< reco::BeamSpot >	tok_bs_

edm::EDGetTokenT < CaloTowerCollection >	tok_cala_

edm::ESGetToken< CaloTopology, CaloTopologyRecord >	tok_caloTopology_

edm::ESGetToken < HcalDDDRecConstants, HcalRecNumberingRecord >	tok_ddrec_

edm::EDGetTokenT < EcalRecHitCollection >	tok_EB_

edm::ESGetToken < EcalChannelStatus, EcalChannelStatusRcd >	tok_ecalChStatus_

edm::EDGetTokenT < EcalRecHitCollection >	tok_EE_

edm::EDGetTokenT < GenEventInfoProduct >	tok_ew_

edm::EDGetTokenT < reco::TrackCollection >	tok_genTrack_

edm::ESGetToken< CaloGeometry, CaloGeometryRecord >	tok_geom_

edm::EDGetTokenT < HBHERecHitCollection >	tok_hbhe_

edm::ESGetToken< HcalTopology, HcalRecNumberingRecord >	tok_htopo_

edm::EDGetTokenT < reco::MuonCollection >	tok_Muon_

edm::EDGetTokenT < reco::GenParticleCollection >	tok_parts_

edm::EDGetTokenT < reco::VertexCollection >	tok_recVtx_

edm::ESGetToken< HcalRespCorrs, HcalRespCorrsRcd >	tok_resp_

edm::ESGetToken < EcalSeverityLevelAlgo, EcalSeverityLevelAlgoRcd >	tok_sevlv_

edm::EDGetTokenT < trigger::TriggerEvent >	tok_trigEvt_

edm::EDGetTokenT < edm::TriggerResults >	tok_trigRes_

TTree *	tree

TTree *	tree2

const edm::InputTag	triggerEvent_

const std::vector< std::string >	trigNames_

const bool	unCorrect_

const bool	useL1Trigger_

const int	useRaw_

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

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept

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

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

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

template<typename 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<Transition Tr = Transition::Event>
constexpr auto	esConsumes () noexcept

template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag) noexcept

template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...

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

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

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

void	resetItemsToGetFrom (BranchType iType)

Detailed Description

Definition at line 91 of file HcalIsoTrkAnalyzer.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 215 of file HcalIsoTrkAnalyzer.cc.

     : trigNames_(iConfig.getParameter<std::vector<std::string>>("triggers")),
       theTrackQuality_(iConfig.getParameter<std::string>("trackQuality")),
       processName_(iConfig.getParameter<std::string>("processName")),
       l1Filter_(iConfig.getParameter<std::string>("l1Filter")),
       l2Filter_(iConfig.getParameter<std::string>("l2Filter")),
       l3Filter_(iConfig.getParameter<std::string>("l3Filter")),
       a_coneR_(iConfig.getParameter<double>("coneRadius")),
       a_mipR_(iConfig.getParameter<double>("coneRadiusMIP")),
       a_mipR2_(iConfig.getParameter<double>("coneRadiusMIP2")),
       a_mipR3_(iConfig.getParameter<double>("coneRadiusMIP3")),
       a_mipR4_(iConfig.getParameter<double>("coneRadiusMIP4")),
       a_mipR5_(iConfig.getParameter<double>("coneRadiusMIP5")),
       pTrackMin_(iConfig.getParameter<double>("minimumTrackP")),
       eEcalMax_(iConfig.getParameter<double>("maximumEcalEnergy")),
       maxRestrictionP_(iConfig.getParameter<double>("maxTrackP")),
       slopeRestrictionP_(iConfig.getParameter<double>("slopeTrackP")),
       hcalScale_(iConfig.getUntrackedParameter<double>("hHcalScale", 1.0)),
       eIsolate1_(iConfig.getParameter<double>("isolationEnergyTight")),
       eIsolate2_(iConfig.getParameter<double>("isolationEnergyLoose")),
       pTrackLow_(iConfig.getParameter<double>("momentumLow")),
       pTrackHigh_(iConfig.getParameter<double>("momentumHigh")),
       prescaleLow_(iConfig.getParameter<int>("prescaleLow")),
       prescaleHigh_(iConfig.getParameter<int>("prescaleHigh")),
       useRaw_(iConfig.getUntrackedParameter<int>("useRaw", 0)),
       dataType_(iConfig.getUntrackedParameter<int>("dataType", 0)),
       mode_(iConfig.getUntrackedParameter<int>("outMode", 11)),
       ignoreTrigger_(iConfig.getUntrackedParameter<bool>("ignoreTriggers", false)),
       useL1Trigger_(iConfig.getUntrackedParameter<bool>("useL1Trigger", false)),
       unCorrect_(iConfig.getUntrackedParameter<bool>("unCorrect", false)),
       collapseDepth_(iConfig.getUntrackedParameter<bool>("collapseDepth", false)),
       hitEthrEB_(iConfig.getParameter<double>("EBHitEnergyThreshold")),
       hitEthrEE0_(iConfig.getParameter<double>("EEHitEnergyThreshold0")),
       hitEthrEE1_(iConfig.getParameter<double>("EEHitEnergyThreshold1")),
       hitEthrEE2_(iConfig.getParameter<double>("EEHitEnergyThreshold2")),
       hitEthrEE3_(iConfig.getParameter<double>("EEHitEnergyThreshold3")),
       hitEthrEELo_(iConfig.getParameter<double>("EEHitEnergyThresholdLow")),
       hitEthrEEHi_(iConfig.getParameter<double>("EEHitEnergyThresholdHigh")),
       triggerEvent_(iConfig.getParameter<edm::InputTag>("labelTriggerEvent")),
       theTriggerResultsLabel_(iConfig.getParameter<edm::InputTag>("labelTriggerResult")),
       labelGenTrack_(iConfig.getParameter<std::string>("labelTrack")),
       labelRecVtx_(iConfig.getParameter<std::string>("labelVertex")),
       labelEB_(iConfig.getParameter<std::string>("labelEBRecHit")),
       labelEE_(iConfig.getParameter<std::string>("labelEERecHit")),
       labelHBHE_(iConfig.getParameter<std::string>("labelHBHERecHit")),
       labelTower_(iConfig.getParameter<std::string>("labelCaloTower")),
       l1TrigName_(iConfig.getUntrackedParameter<std::string>("l1TrigName", "L1_SingleJet60")),
       oldID_(iConfig.getUntrackedParameter<std::vector<int>>("oldID")),
       newDepth_(iConfig.getUntrackedParameter<std::vector<int>>("newDepth")),
       hep17_(iConfig.getUntrackedParameter<bool>("hep17")),
       debEvents_(iConfig.getParameter<std::vector<int>>("debugEvents")),
       nRun_(0),
       nLow_(0),
       nHigh_(0),
       hdc_(nullptr) {
   usesResource(TFileService::kSharedResource);
 
   //now do whatever initialization is needed
   const double isolationRadius(28.9), innerR(10.0), outerR(30.0);
   reco::TrackBase::TrackQuality trackQuality_ = reco::TrackBase::qualityByName(theTrackQuality_);
   selectionParameter_.minPt = iConfig.getParameter<double>("minTrackPt");
   ;
   selectionParameter_.minQuality = trackQuality_;
   selectionParameter_.maxDxyPV = iConfig.getParameter<double>("maxDxyPV");
   selectionParameter_.maxDzPV = iConfig.getParameter<double>("maxDzPV");
   selectionParameter_.maxChi2 = iConfig.getParameter<double>("maxChi2");
   selectionParameter_.maxDpOverP = iConfig.getParameter<double>("maxDpOverP");
   selectionParameter_.minOuterHit = iConfig.getParameter<int>("minOuterHit");
   selectionParameter_.minLayerCrossed = iConfig.getParameter<int>("minLayerCrossed");
   selectionParameter_.maxInMiss = iConfig.getParameter<int>("maxInMiss");
   selectionParameter_.maxOutMiss = iConfig.getParameter<int>("maxOutMiss");
   a_charIsoR_ = a_coneR_ + isolationRadius;
   a_coneR1_ = a_coneR_ + innerR;
   a_coneR2_ = a_coneR_ + outerR;
   // Different isolation cuts are described in DN-2016/029
   // Tight cut uses 2 GeV; Loose cut uses 10 GeV
   // Eta dependent cut uses (maxRestrictionP_ * exp(|ieta|*log(2.5)/18))
   // with the factor for exponential slopeRestrictionP_ = log(2.5)/18
   // maxRestrictionP_ = 8 GeV as came from a study
   std::string labelBS = iConfig.getParameter<std::string>("labelBeamSpot");
   std::string modnam = iConfig.getUntrackedParameter<std::string>("moduleName", "");
   std::string prdnam = iConfig.getUntrackedParameter<std::string>("producerName", "");
   edm::InputTag algTag = iConfig.getParameter<edm::InputTag>("algInputTag");
   edm::InputTag extTag = iConfig.getParameter<edm::InputTag>("extInputTag");
   std::string labelMuon = iConfig.getParameter<std::string>("labelMuon");
 
   for (unsigned int k = 0; k < oldID_.size(); ++k) {
     oldDet_.emplace_back((oldID_[k] / 10000) % 10);
     oldEta_.emplace_back((oldID_[k] / 100) % 100);
     oldDepth_.emplace_back(oldID_[k] % 100);
   }
 
   l1GtUtils_ = new l1t::L1TGlobalUtil(iConfig, consumesCollector(), *this, algTag, extTag, l1t::UseEventSetupIn::Event);
   // define tokens for access
   tok_trigEvt_ = consumes<trigger::TriggerEvent>(triggerEvent_);
   tok_trigRes_ = consumes<edm::TriggerResults>(theTriggerResultsLabel_);
   tok_bs_ = consumes<reco::BeamSpot>(labelBS);
   tok_genTrack_ = consumes<reco::TrackCollection>(labelGenTrack_);
   tok_ew_ = consumes<GenEventInfoProduct>(edm::InputTag("generator"));
   tok_parts_ = consumes<reco::GenParticleCollection>(edm::InputTag("genParticles"));
   tok_cala_ = consumes<CaloTowerCollection>(labelTower_);
   tok_alg_ = consumes<BXVector<GlobalAlgBlk>>(algTag);
   tok_Muon_ = consumes<reco::MuonCollection>(labelMuon);
 
   if (modnam.empty()) {
     tok_recVtx_ = consumes<reco::VertexCollection>(labelRecVtx_);
     tok_EB_ = consumes<EcalRecHitCollection>(edm::InputTag("ecalRecHit", labelEB_));
     tok_EE_ = consumes<EcalRecHitCollection>(edm::InputTag("ecalRecHit", labelEE_));
     tok_hbhe_ = consumes<HBHERecHitCollection>(labelHBHE_);
     edm::LogVerbatim("HcalIsoTrack") << "Labels used " << triggerEvent_ << " " << theTriggerResultsLabel_ << " "
                                      << labelBS << " " << labelRecVtx_ << " " << labelGenTrack_ << " "
                                      << edm::InputTag("ecalRecHit", labelEB_) << " "
                                      << edm::InputTag("ecalRecHit", labelEE_) << " " << labelHBHE_ << " " << labelTower_
                                      << " " << labelMuon;
   } else {
     tok_recVtx_ = consumes<reco::VertexCollection>(edm::InputTag(modnam, labelRecVtx_, prdnam));
     tok_EB_ = consumes<EcalRecHitCollection>(edm::InputTag(modnam, labelEB_, prdnam));
     tok_EE_ = consumes<EcalRecHitCollection>(edm::InputTag(modnam, labelEE_, prdnam));
     tok_hbhe_ = consumes<HBHERecHitCollection>(edm::InputTag(modnam, labelHBHE_, prdnam));
     edm::LogVerbatim("HcalIsoTrack") << "Labels used " << triggerEvent_ << " " << theTriggerResultsLabel_ << " "
                                      << labelBS << " " << edm::InputTag(modnam, labelRecVtx_, prdnam) << " "
                                      << labelGenTrack_ << " " << edm::InputTag(modnam, labelEB_, prdnam) << " "
                                      << edm::InputTag(modnam, labelEE_, prdnam) << " "
                                      << edm::InputTag(modnam, labelHBHE_, prdnam) << " " << labelTower_ << " "
                                      << labelMuon;
   }
 
   tok_ddrec_ = esConsumes<HcalDDDRecConstants, HcalRecNumberingRecord, edm::Transition::BeginRun>();
   tok_bFieldH_ = esConsumes<MagneticField, IdealMagneticFieldRecord>();
   tok_ecalChStatus_ = esConsumes<EcalChannelStatus, EcalChannelStatusRcd>();
   tok_sevlv_ = esConsumes<EcalSeverityLevelAlgo, EcalSeverityLevelAlgoRcd>();
   tok_geom_ = esConsumes<CaloGeometry, CaloGeometryRecord>();
   tok_caloTopology_ = esConsumes<CaloTopology, CaloTopologyRecord>();
   tok_htopo_ = esConsumes<HcalTopology, HcalRecNumberingRecord>();
   tok_resp_ = esConsumes<HcalRespCorrs, HcalRespCorrsRcd>();
 
   edm::LogVerbatim("HcalIsoTrack")
       << "Parameters read from config file \n"
       << "\t minPt " << selectionParameter_.minPt << "\t theTrackQuality " << theTrackQuality_ << "\t minQuality "
       << selectionParameter_.minQuality << "\t maxDxyPV " << selectionParameter_.maxDxyPV << "\t maxDzPV "
       << selectionParameter_.maxDzPV << "\t maxChi2 " << selectionParameter_.maxChi2 << "\t maxDpOverP "
       << selectionParameter_.maxDpOverP << "\t minOuterHit " << selectionParameter_.minOuterHit << "\t minLayerCrossed "
       << selectionParameter_.minLayerCrossed << "\t maxInMiss " << selectionParameter_.maxInMiss << "\t maxOutMiss "
       << selectionParameter_.maxOutMiss << "\t a_coneR " << a_coneR_ << ":" << a_coneR1_ << ":" << a_coneR2_
       << "\t a_charIsoR " << a_charIsoR_ << "\t a_mipR " << a_mipR_ << "\t a_mipR2 " << a_mipR2_ << "\t a_mipR3 "
       << a_mipR3_ << "\t a_mipR4 " << a_mipR4_ << "\t a_mipR5 " << a_mipR5_ << "\n pTrackMin_ " << pTrackMin_
       << "\t eEcalMax_ " << eEcalMax_ << "\t maxRestrictionP_ " << maxRestrictionP_ << "\t slopeRestrictionP_ "
       << slopeRestrictionP_ << "\t eIsolateStrong_ " << eIsolate1_ << "\t eIsolateSoft_ " << eIsolate2_
       << "\t hcalScale_ " << hcalScale_ << "\n\t momentumLow_ " << pTrackLow_ << "\t prescaleLow_ " << prescaleLow_
       << "\t momentumHigh_ " << pTrackHigh_ << "\t prescaleHigh_ " << prescaleHigh_ << "\n\t useRaw_ " << useRaw_
       << "\t ignoreTrigger_ " << ignoreTrigger_ << "\n\t useL1Trigegr_ " << useL1Trigger_ << "\t dataType_      "
       << dataType_ << "\t mode_          " << mode_ << "\t unCorrect_     " << unCorrect_ << "\t collapseDepth_ "
       << collapseDepth_ << "\t L1TrigName_    " << l1TrigName_ << "\nThreshold for EB " << hitEthrEB_ << " EE "
       << hitEthrEE0_ << ":" << hitEthrEE1_ << ":" << hitEthrEE2_ << ":" << hitEthrEE3_ << ":" << hitEthrEELo_ << ":"
       << hitEthrEEHi_ << " and " << debEvents_.size() << " events to be debugged";
   edm::LogVerbatim("HcalIsoTrack") << "Process " << processName_ << " L1Filter:" << l1Filter_
                                    << " L2Filter:" << l2Filter_ << " L3Filter:" << l3Filter_;
   for (unsigned int k = 0; k < trigNames_.size(); ++k) {
     edm::LogVerbatim("HcalIsoTrack") << "Trigger[" << k << "] " << trigNames_[k];
   }
   edm::LogVerbatim("HcalIsoTrack") << oldID_.size() << " DetIDs to be corrected with HEP17 flag:" << hep17_;
   for (unsigned int k = 0; k < oldID_.size(); ++k)
     edm::LogVerbatim("HcalIsoTrack") << "[" << k << "] Det " << oldDet_[k] << " EtaAbs " << oldEta_[k] << " Depth "
                                      << oldDepth_[k] << ":" << newDepth_[k];
 
   for (int i = 0; i < 10; i++)
     phibins_.push_back(-M_PI + 0.1 * (2 * i + 1) * M_PI);
   for (int i = 0; i < 8; ++i)
     etabins_.push_back(-2.1 + 0.6 * i);
   etadist_ = etabins_[1] - etabins_[0];
   phidist_ = phibins_[1] - phibins_[0];
   etahalfdist_ = 0.5 * etadist_;
   phihalfdist_ = 0.5 * phidist_;
   edm::LogVerbatim("HcalIsoTrack") << "EtaDist " << etadist_ << " " << etahalfdist_ << " PhiDist " << phidist_ << " "
                                    << phihalfdist_;
   unsigned int k1(0), k2(0);
   for (auto phi : phibins_) {
     edm::LogVerbatim("HcalIsoTrack") << "phibin_[" << k1 << "] " << phi;
     ++k1;
   }
   for (auto eta : etabins_) {
     edm::LogVerbatim("HcalIsoTrack") << "etabin_[" << k2 << "] " << eta;
     ++k2;
   }
 }

HcalIsoTrkAnalyzer::~HcalIsoTrkAnalyzer ( )

inlineoverride

Definition at line 94 of file HcalIsoTrkAnalyzer.cc.

94 {}

Member Function Documentation

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

overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 401 of file HcalIsoTrkAnalyzer.cc.

                                                                                     {
   t_Run = iEvent.id().run();
   t_Event = iEvent.id().event();
   t_DataType = dataType_;
   debug_ = (debEvents_.empty())
                ? true
                : (std::find(debEvents_.begin(), debEvents_.end(), iEvent.id().event()) != debEvents_.end());
 #ifdef EDM_ML_DEBUG
   if (debug_)
     edm::LogVerbatim("HcalIsoTrack") << "Run " << t_Run << " Event " << t_Event << " type " << t_DataType
                                      << " Luminosity " << iEvent.luminosityBlock() << " Bunch "
                                      << iEvent.bunchCrossing();
 #endif
   //Get magnetic field and ECAL channel status
   const MagneticField* bField = &iSetup.getData(tok_bFieldH_);
   const EcalChannelStatus* theEcalChStatus = &iSetup.getData(tok_ecalChStatus_);
   const EcalSeverityLevelAlgo* theEcalSevlv = &iSetup.getData(tok_sevlv_);
 
   // get handles to calogeometry and calotopology
   const CaloGeometry* geo = &iSetup.getData(tok_geom_);
   const CaloTopology* caloTopology = &iSetup.getData(tok_caloTopology_);
   const HcalTopology* theHBHETopology = &iSetup.getData(tok_htopo_);
   const HcalRespCorrs* respCorrs = &iSetup.getData(tok_resp_);
 
   //=== genParticle information
   edm::Handle<reco::GenParticleCollection> genParticles;
   iEvent.getByToken(tok_parts_, genParticles);
 
   bool okC(true);
   //Get track collection
   edm::Handle<reco::TrackCollection> trkCollection;
   iEvent.getByToken(tok_genTrack_, trkCollection);
   if (!trkCollection.isValid()) {
     edm::LogWarning("HcalIsoTrack") << "Cannot access the collection " << labelGenTrack_;
     okC = false;
   }
 
   //Get muon collection
   edm::Handle<reco::MuonCollection> muonh;
   iEvent.getByToken(tok_Muon_, muonh);
 
   //event weight for FLAT sample
   t_EventWeight = 1.0;
   edm::Handle<GenEventInfoProduct> genEventInfo;
   iEvent.getByToken(tok_ew_, genEventInfo);
   if (genEventInfo.isValid())
     t_EventWeight = genEventInfo->weight();
 
   //Define the best vertex and the beamspot
   edm::Handle<reco::VertexCollection> recVtxs;
   iEvent.getByToken(tok_recVtx_, recVtxs);
   edm::Handle<reco::BeamSpot> beamSpotH;
   iEvent.getByToken(tok_bs_, beamSpotH);
   math::XYZPoint leadPV(0, 0, 0);
   t_goodPV = t_nVtx = 0;
   if (recVtxs.isValid() && !(recVtxs->empty())) {
     t_nVtx = recVtxs->size();
     for (unsigned int k = 0; k < recVtxs->size(); ++k) {
       if (!((*recVtxs)[k].isFake()) && ((*recVtxs)[k].ndof() > 4)) {
         if (t_goodPV == 0)
           leadPV = math::XYZPoint((*recVtxs)[k].x(), (*recVtxs)[k].y(), (*recVtxs)[k].z());
         t_goodPV++;
       }
     }
   }
   if (t_goodPV == 0 && beamSpotH.isValid()) {
     leadPV = beamSpotH->position();
   }
   t_allvertex = t_goodPV;
 #ifdef EDM_ML_DEBUG
   if (debug_) {
     edm::LogVerbatim("HcalIsoTrack") << "Primary Vertex " << leadPV << " out of " << t_goodPV << " vertex";
     if (beamSpotH.isValid())
       edm::LogVerbatim("HcalIsoTrack") << " Beam Spot " << beamSpotH->position();
   }
 #endif
   // RecHits
   edm::Handle<EcalRecHitCollection> barrelRecHitsHandle;
   iEvent.getByToken(tok_EB_, barrelRecHitsHandle);
   if (!barrelRecHitsHandle.isValid()) {
     edm::LogWarning("HcalIsoTrack") << "Cannot access the collection " << labelEB_;
     okC = false;
   }
   edm::Handle<EcalRecHitCollection> endcapRecHitsHandle;
   iEvent.getByToken(tok_EE_, endcapRecHitsHandle);
   if (!endcapRecHitsHandle.isValid()) {
     edm::LogWarning("HcalIsoTrack") << "Cannot access the collection " << labelEE_;
     okC = false;
   }
   edm::Handle<HBHERecHitCollection> hbhe;
   iEvent.getByToken(tok_hbhe_, hbhe);
   if (!hbhe.isValid()) {
     edm::LogWarning("HcalIsoTrack") << "Cannot access the collection " << labelHBHE_;
     okC = false;
   }
   edm::Handle<CaloTowerCollection> caloTower;
   iEvent.getByToken(tok_cala_, caloTower);
 
   //Propagate tracks to calorimeter surface)
   std::vector<spr::propagatedTrackDirection> trkCaloDirections;
   spr::propagateCALO(trkCollection, geo, bField, theTrackQuality_, trkCaloDirections, false);
   std::vector<spr::propagatedTrackID> trkCaloDets;
   spr::propagateCALO(trkCollection, geo, bField, theTrackQuality_, trkCaloDets, false);
   std::vector<math::XYZTLorentzVector> vecL1, vecL3;
   t_RunNo = iEvent.id().run();
   t_EventNo = iEvent.id().event();
   t_Tracks = trkCollection->size();
   t_TracksProp = trkCaloDirections.size();
   t_ietaAll->clear();
   t_ietaGood->clear();
   t_trackType->clear();
   t_trgbits->clear();
   t_hltbits->clear();
 #ifdef EDM_ML_DEBUG
   if (debug_)
     edm::LogVerbatim("HcalIsoTrack") << "# of propagated tracks " << t_TracksProp << " out of " << t_Tracks
                                      << " with Trigger " << ignoreTrigger_;
 #endif
 
   //Trigger
   t_trgbits->assign(trigNames_.size(), false);
   t_hltbits->assign(trigNames_.size(), false);
   t_TracksSaved = t_TracksLoose = t_TracksTight = 0;
   t_L1Bit = true;
   t_TrigPass = false;
 
   edm::Handle<edm::TriggerResults> triggerResults;
   if (!ignoreTrigger_) {
     //L1
     l1GtUtils_->retrieveL1(iEvent, iSetup, tok_alg_);
     const std::vector<std::pair<std::string, bool>>& finalDecisions = l1GtUtils_->decisionsFinal();
     for (const auto& decision : finalDecisions) {
       if (decision.first.find(l1TrigName_) != std::string::npos) {
         t_L1Bit = decision.second;
         break;
       }
     }
 #ifdef EDM_ML_DEBUG
     if (debug_)
       edm::LogVerbatim("HcalIsoTrack") << "Trigger Information for " << l1TrigName_ << " is " << t_L1Bit
                                        << " from a list of " << finalDecisions.size() << " decisions";
 #endif
 
     //HLT
     iEvent.getByToken(tok_trigRes_, triggerResults);
     if (triggerResults.isValid()) {
       const edm::TriggerNames& triggerNames = iEvent.triggerNames(*triggerResults);
       const std::vector<std::string>& names = triggerNames.triggerNames();
       if (!trigNames_.empty()) {
         for (unsigned int iHLT = 0; iHLT < triggerResults->size(); iHLT++) {
           int hlt = triggerResults->accept(iHLT);
           for (unsigned int i = 0; i < trigNames_.size(); ++i) {
             if (names[iHLT].find(trigNames_[i]) != std::string::npos) {
               t_trgbits->at(i) = (hlt > 0);
               t_hltbits->at(i) = (hlt > 0);
               if (hlt > 0)
                 t_TrigPass = true;
 #ifdef EDM_ML_DEBUG
               if (debug_)
                 edm::LogVerbatim("HcalIsoTrack")
                     << "This trigger " << names[iHLT] << " Flag " << hlt << ":" << t_trgbits->at(i);
 #endif
             }
           }
         }
       }
     }
 #ifdef EDM_ML_DEBUG
     if (debug_)
       edm::LogVerbatim("HcalIsoTrack") << "HLT Information shows " << t_TrigPass << ":" << trigNames_.empty() << ":"
                                        << okC;
 #endif
   }
 
   std::array<int, 3> ntksave{{0, 0, 0}};
   if (ignoreTrigger_ || useL1Trigger_) {
     t_l1pt = t_l1eta = t_l1phi = 0;
     t_l3pt = t_l3eta = t_l3phi = 0;
     if (ignoreTrigger_ || t_L1Bit)
       ntksave = fillTree(vecL1,
                          vecL3,
                          leadPV,
                          trkCaloDirections,
                          trkCaloDets,
                          geo,
                          caloTopology,
                          theHBHETopology,
                          theEcalChStatus,
                          theEcalSevlv,
                          barrelRecHitsHandle,
                          endcapRecHitsHandle,
                          hbhe,
                          caloTower,
                          genParticles,
                          respCorrs,
                          muonh);
     t_TracksSaved = ntksave[0];
     t_TracksLoose = ntksave[1];
     t_TracksTight = ntksave[2];
   } else {
     trigger::TriggerEvent triggerEvent;
     edm::Handle<trigger::TriggerEvent> triggerEventHandle;
     iEvent.getByToken(tok_trigEvt_, triggerEventHandle);
     if (!triggerEventHandle.isValid()) {
       edm::LogWarning("HcalIsoTrack") << "Error! Can't get the product " << triggerEvent_.label();
     } else if (okC) {
       triggerEvent = *(triggerEventHandle.product());
       const trigger::TriggerObjectCollection& TOC(triggerEvent.getObjects());
       bool done(false);
       if (triggerResults.isValid()) {
         std::vector<std::string> modules;
         const edm::TriggerNames& triggerNames = iEvent.triggerNames(*triggerResults);
         const std::vector<std::string>& names = triggerNames.triggerNames();
         for (unsigned int iHLT = 0; iHLT < triggerResults->size(); iHLT++) {
           bool ok = (t_TrigPass) || (trigNames_.empty());
           if (ok) {
             unsigned int triggerindx = hltConfig_.triggerIndex(names[iHLT]);
             const std::vector<std::string>& moduleLabels(hltConfig_.moduleLabels(triggerindx));
             std::vector<math::XYZTLorentzVector> vecL2;
             vecL1.clear();
             vecL3.clear();
             //loop over all trigger filters in event (i.e. filters passed)
             for (unsigned int ifilter = 0; ifilter < triggerEvent.sizeFilters(); ++ifilter) {
               std::vector<int> Keys;
               std::string label = triggerEvent.filterTag(ifilter).label();
               //loop over keys to objects passing this filter
               for (unsigned int imodule = 0; imodule < moduleLabels.size(); imodule++) {
                 if (label.find(moduleLabels[imodule]) != std::string::npos) {
 #ifdef EDM_ML_DEBUG
                   if (debug_)
                     edm::LogVerbatim("HcalIsoTrack") << "FilterName " << label;
 #endif
                   for (unsigned int ifiltrKey = 0; ifiltrKey < triggerEvent.filterKeys(ifilter).size(); ++ifiltrKey) {
                     Keys.push_back(triggerEvent.filterKeys(ifilter)[ifiltrKey]);
                     const trigger::TriggerObject& TO(TOC[Keys[ifiltrKey]]);
                     math::XYZTLorentzVector v4(TO.px(), TO.py(), TO.pz(), TO.energy());
                     if (label.find(l2Filter_) != std::string::npos) {
                       vecL2.push_back(v4);
                     } else if (label.find(l3Filter_) != std::string::npos) {
                       vecL3.push_back(v4);
                     } else if ((label.find(l1Filter_) != std::string::npos) || (l1Filter_.empty())) {
                       vecL1.push_back(v4);
                     }
 #ifdef EDM_ML_DEBUG
                     if (debug_)
                       edm::LogVerbatim("HcalIsoTrack")
                           << "key " << ifiltrKey << " : pt " << TO.pt() << " eta " << TO.eta() << " phi " << TO.phi()
                           << " mass " << TO.mass() << " Id " << TO.id();
 #endif
                   }
 #ifdef EDM_ML_DEBUG
                   if (debug_)
                     edm::LogVerbatim("HcalIsoTrack")
                         << "sizes " << vecL1.size() << ":" << vecL2.size() << ":" << vecL3.size();
 #endif
                 }
               }
             }
             math::XYZTLorentzVector mindRvec1;
             double mindR1(999);
             for (unsigned int i = 0; i < vecL2.size(); i++) {
               double dr = dR(vecL1[0], vecL2[i]);
 #ifdef EDM_ML_DEBUG
               if (debug_)
                 edm::LogVerbatim("HcalIsoTrack") << "lvl2[" << i << "] dR " << dr;
 #endif
               if (dr < mindR1) {
                 mindR1 = dr;
                 mindRvec1 = vecL2[i];
               }
             }
 #ifdef EDM_ML_DEBUG
             if (debug_)
               edm::LogVerbatim("HcalIsoTrack") << "L2 object closest to L1 " << mindRvec1 << " at Dr " << mindR1;
 #endif
 
             if (!vecL1.empty()) {
               t_l1pt = vecL1[0].pt();
               t_l1eta = vecL1[0].eta();
               t_l1phi = vecL1[0].phi();
             } else {
               t_l1pt = t_l1eta = t_l1phi = 0;
             }
             if (!vecL3.empty()) {
               t_l3pt = vecL3[0].pt();
               t_l3eta = vecL3[0].eta();
               t_l3phi = vecL3[0].phi();
             } else {
               t_l3pt = t_l3eta = t_l3phi = 0;
             }
             // Now fill in the tree for each selected track
             if (!done) {
               ntksave = fillTree(vecL1,
                                  vecL3,
                                  leadPV,
                                  trkCaloDirections,
                                  trkCaloDets,
                                  geo,
                                  caloTopology,
                                  theHBHETopology,
                                  theEcalChStatus,
                                  theEcalSevlv,
                                  barrelRecHitsHandle,
                                  endcapRecHitsHandle,
                                  hbhe,
                                  caloTower,
                                  genParticles,
                                  respCorrs,
                                  muonh);
               t_TracksSaved += ntksave[0];
               t_TracksLoose += ntksave[1];
               t_TracksTight += ntksave[2];
               done = true;
             }
           }
         }
       }
     }
   }
 #ifdef EDM_ML_DEBUG
   if (debug_)
     edm::LogVerbatim("HcalIsoTrack") << "Final results on selected tracks " << t_TracksSaved << ":" << t_TracksLoose
                                      << ":" << t_TracksTight;
 #endif
   t_TrigPassSel = (t_TracksSaved > 0);
   tree2->Fill();
 }

void HcalIsoTrkAnalyzer::beginJob ( void )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 730 of file HcalIsoTrkAnalyzer.cc.

References fs, TFileService::make(), t_allvertex, t_DataType, t_DetIds, t_DetIds1, t_DetIds3, t_eAnnular, t_eHcal, t_eHcal10, t_eHcal30, t_emaxNearP, t_eMipDR, t_eMipDR2, t_eMipDR3, t_eMipDR4, t_eMipDR5, t_Event, t_EventNo, t_EventWeight, t_gentrackP, t_goodPV, t_hAnnular, t_HitEnergies, t_HitEnergies1, t_HitEnergies3, t_hltbits, t_hmaxNearP, t_ieta, t_ietaAll, t_ietaGood, t_iphi, t_L1Bit, t_l1eta, t_l1phi, t_l1pt, t_l3eta, t_l3phi, t_l3pt, t_mindR1, t_mindR2, t_nTrk, t_nVtx, t_p, t_phi, t_pt, t_qltyFlag, t_qltyMissFlag, t_qltyPVFlag, t_rhoh, t_Run, t_RunNo, t_selectTk, t_Tracks, t_TracksLoose, t_TracksProp, t_TracksSaved, t_TracksTight, t_trackType, t_trgbits, t_TrigPass, t_TrigPassSel, tree, and tree2.

                                   {
   tree = fs->make<TTree>("CalibTree", "CalibTree");
 
   tree->Branch("t_Run", &t_Run, "t_Run/I");
   tree->Branch("t_Event", &t_Event, "t_Event/I");
   tree->Branch("t_DataType", &t_DataType, "t_DataType/I");
   tree->Branch("t_ieta", &t_ieta, "t_ieta/I");
   tree->Branch("t_iphi", &t_iphi, "t_iphi/I");
   tree->Branch("t_EventWeight", &t_EventWeight, "t_EventWeight/D");
   tree->Branch("t_nVtx", &t_nVtx, "t_nVtx/I");
   tree->Branch("t_nTrk", &t_nTrk, "t_nTrk/I");
   tree->Branch("t_goodPV", &t_goodPV, "t_goodPV/I");
   tree->Branch("t_l1pt", &t_l1pt, "t_l1pt/D");
   tree->Branch("t_l1eta", &t_l1eta, "t_l1eta/D");
   tree->Branch("t_l1phi", &t_l1phi, "t_l1phi/D");
   tree->Branch("t_l3pt", &t_l3pt, "t_l3pt/D");
   tree->Branch("t_l3eta", &t_l3eta, "t_l3eta/D");
   tree->Branch("t_l3phi", &t_l3phi, "t_l3phi/D");
   tree->Branch("t_p", &t_p, "t_p/D");
   tree->Branch("t_pt", &t_pt, "t_pt/D");
   tree->Branch("t_phi", &t_phi, "t_phi/D");
   tree->Branch("t_mindR1", &t_mindR1, "t_mindR1/D");
   tree->Branch("t_mindR2", &t_mindR2, "t_mindR2/D");
   tree->Branch("t_eMipDR", &t_eMipDR, "t_eMipDR/D");
   tree->Branch("t_eMipDR2", &t_eMipDR2, "t_eMipDR2/D");
   tree->Branch("t_eMipDR3", &t_eMipDR3, "t_eMipDR3/D");
   tree->Branch("t_eMipDR4", &t_eMipDR4, "t_eMipDR4/D");
   tree->Branch("t_eMipDR5", &t_eMipDR5, "t_eMipDR5/D");
   tree->Branch("t_eHcal", &t_eHcal, "t_eHcal/D");
   tree->Branch("t_eHcal10", &t_eHcal10, "t_eHcal10/D");
   tree->Branch("t_eHcal30", &t_eHcal30, "t_eHcal30/D");
   tree->Branch("t_hmaxNearP", &t_hmaxNearP, "t_hmaxNearP/D");
   tree->Branch("t_emaxNearP", &t_emaxNearP, "t_emaxNearP/D");
   tree->Branch("t_eAnnular", &t_eAnnular, "t_eAnnular/D");
   tree->Branch("t_hAnnular", &t_hAnnular, "t_hAnnular/D");
   tree->Branch("t_rhoh", &t_rhoh, "t_rhoh/D");
   tree->Branch("t_selectTk", &t_selectTk, "t_selectTk/O");
   tree->Branch("t_qltyFlag", &t_qltyFlag, "t_qltyFlag/O");
   tree->Branch("t_qltyMissFlag", &t_qltyMissFlag, "t_qltyMissFlag/O");
   tree->Branch("t_qltyPVFlag", &t_qltyPVFlag, "t_qltyPVFlag/O");
   tree->Branch("t_gentrackP", &t_gentrackP, "t_gentrackP/D");
 
   t_DetIds = new std::vector<unsigned int>();
   t_DetIds1 = new std::vector<unsigned int>();
   t_DetIds3 = new std::vector<unsigned int>();
   t_HitEnergies = new std::vector<double>();
   t_HitEnergies1 = new std::vector<double>();
   t_HitEnergies3 = new std::vector<double>();
   t_trgbits = new std::vector<bool>();
   tree->Branch("t_DetIds", "std::vector<unsigned int>", &t_DetIds);
   tree->Branch("t_HitEnergies", "std::vector<double>", &t_HitEnergies);
   tree->Branch("t_trgbits", "std::vector<bool>", &t_trgbits);
   tree->Branch("t_DetIds1", "std::vector<unsigned int>", &t_DetIds1);
   tree->Branch("t_DetIds3", "std::vector<unsigned int>", &t_DetIds3);
   tree->Branch("t_HitEnergies1", "std::vector<double>", &t_HitEnergies1);
   tree->Branch("t_HitEnergies3", "std::vector<double>", &t_HitEnergies3);
 
   tree2 = fs->make<TTree>("EventInfo", "Event Information");
 
   tree2->Branch("t_RunNo", &t_RunNo, "t_RunNo/i");
   tree2->Branch("t_EventNo", &t_EventNo, "t_EventNo/i");
   tree2->Branch("t_Tracks", &t_Tracks, "t_Tracks/I");
   tree2->Branch("t_TracksProp", &t_TracksProp, "t_TracksProp/I");
   tree2->Branch("t_TracksSaved", &t_TracksSaved, "t_TracksSaved/I");
   tree2->Branch("t_TracksLoose", &t_TracksLoose, "t_TracksLoose/I");
   tree2->Branch("t_TracksTight", &t_TracksTight, "t_TracksTight/I");
   tree2->Branch("t_TrigPass", &t_TrigPass, "t_TrigPass/O");
   tree2->Branch("t_TrigPassSel", &t_TrigPassSel, "t_TrigPassSel/O");
   tree2->Branch("t_L1Bit", &t_L1Bit, "t_L1Bit/O");
   tree2->Branch("t_allvertex", &t_allvertex, "t_allvertex/I");
   t_hltbits = new std::vector<bool>();
   t_ietaAll = new std::vector<int>();
   t_ietaGood = new std::vector<int>();
   t_trackType = new std::vector<int>();
   tree2->Branch("t_ietaAll", "std::vector<int>", &t_ietaAll);
   tree2->Branch("t_ietaGood", "std::vector<int>", &t_ietaGood);
   tree2->Branch("t_trackType", "std::vector<int>", &t_trackType);
   tree2->Branch("t_hltbits", "std::vector<bool>", &t_hltbits);
 }

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

overrideprivate

Definition at line 811 of file HcalIsoTrkAnalyzer.cc.

References edm::EventSetup::getData(), hdc_, hltConfig_, ignoreTrigger_, HLTConfigProvider::init(), dqmiodumpmetadata::n, nRun_, processName_, edm::RunBase::run(), HLTConfigProvider::size(), tok_ddrec_, HLTConfigProvider::triggerIndex(), and trigNames_.

                                                                                  {
   hdc_ = &iSetup.getData(tok_ddrec_);
 
   if (!ignoreTrigger_) {
     bool changed_(true);
     bool flag = hltConfig_.init(iRun, iSetup, processName_, changed_);
     edm::LogVerbatim("HcalIsoTrack") << "Run[" << nRun_ << "] " << iRun.run() << " process " << processName_
                                      << " init flag " << flag << " change flag " << changed_;
     // check if trigger names in (new) config
     if (changed_) {
 #ifdef EDM_ML_DEBUG
       edm::LogVerbatim("HcalIsoTrack") << "New trigger menu found !!!";
 #endif
       const unsigned int n(hltConfig_.size());
       for (unsigned itrig = 0; itrig < trigNames_.size(); itrig++) {
         unsigned int triggerindx = hltConfig_.triggerIndex(trigNames_[itrig]);
         if (triggerindx >= n) {
           edm::LogWarning("HcalIsoTrack") << trigNames_[itrig] << " " << triggerindx << " does not exist in "
                                           << "the current menu";
 #ifdef EDM_ML_DEBUG
         } else {
           edm::LogVerbatim("HcalIsoTrack") << trigNames_[itrig] << " " << triggerindx << " exists";
 #endif
         }
       }
     }
   }
 }

double HcalIsoTrkAnalyzer::dR	(	math::XYZTLorentzVector &	vec1,
		math::XYZTLorentzVector &	vec2
	)

private

Definition at line 1322 of file HcalIsoTrkAnalyzer.cc.

References reco::deltaR().

Referenced by analyze(), fillTree(), notaMuon(), and trackP().

                                                                                       {
   return reco::deltaR(vec1.eta(), vec1.phi(), vec2.eta(), vec2.phi());
 }

void HcalIsoTrkAnalyzer::endRun	(	edm::Run const &	iRun,
		edm::EventSetup const &
	)

overrideprivate

Definition at line 841 of file HcalIsoTrkAnalyzer.cc.

References nRun_, and edm::RunBase::run().

                                                                         {
   nRun_++;
   edm::LogVerbatim("HcalIsoTrack") << "endRun[" << nRun_ << "] " << iRun.run();
 }

double HcalIsoTrkAnalyzer::eThreshold	(	const DetId &	id,
		const CaloGeometry *	geo
	)		const

private

Definition at line 1374 of file HcalIsoTrkAnalyzer.cc.

References funct::abs(), EcalBarrel, PVValHelper::eta, PV3DBase< T, PVType, FrameType >::eta(), CaloGeometry::getPosition(), hitEthrEB_, hitEthrEE0_, hitEthrEE1_, hitEthrEE2_, hitEthrEE3_, hitEthrEEHi_, and hitEthrEELo_.

Referenced by fillTree().

                                                                                     {
   const GlobalPoint& pos = geo->getPosition(id);
   double eta = std::abs(pos.eta());
   double eThr(hitEthrEB_);
   if (id.subdetId() != EcalBarrel) {
     eThr = (((eta * hitEthrEE3_ + hitEthrEE2_) * eta + hitEthrEE1_) * eta + hitEthrEE0_);
     if (eThr < hitEthrEELo_)
       eThr = hitEthrEELo_;
     else if (eThr > hitEthrEEHi_)
       eThr = hitEthrEEHi_;
   }
   return eThr;
 }

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

static

Definition at line 846 of file HcalIsoTrkAnalyzer.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), edm::ParameterSetDescription::addUntracked(), submitPVResolutionJobs::desc, patZpeak::events, HLT_FULL_cff::InputTag, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                     {
   edm::ParameterSetDescription desc;
   std::vector<std::string> trig = {"HLT_PFJet40",
                                    "HLT_PFJet60",
                                    "HLT_PFJet80",
                                    "HLT_PFJet140",
                                    "HLT_PFJet200",
                                    "HLT_PFJet260",
                                    "HLT_PFJet320",
                                    "HLT_PFJet400",
                                    "HLT_PFJet450",
                                    "HLT_PFJet500"};
   desc.add<std::vector<std::string>>("triggers", trig);
   desc.add<std::string>("processName", "HLT");
   desc.add<std::string>("l1Filter", "");
   desc.add<std::string>("l2Filter", "L2Filter");
   desc.add<std::string>("l3Filter", "Filter");
   // following 10 parameters are parameters to select good tracks
   desc.add<std::string>("trackQuality", "highPurity");
   desc.add<double>("minTrackPt", 1.0);
   desc.add<double>("maxDxyPV", 0.02);
   desc.add<double>("maxDzPV", 0.02);
   desc.add<double>("maxChi2", 5.0);
   desc.add<double>("maxDpOverP", 0.1);
   desc.add<int>("minOuterHit", 4);
   desc.add<int>("minLayerCrossed", 8);
   desc.add<int>("maxInMiss", 0);
   desc.add<int>("maxOutMiss", 0);
   // Minimum momentum of selected isolated track and signal zone
   desc.add<double>("minimumTrackP", 10.0);
   desc.add<double>("coneRadius", 34.98);
   // signal zone in ECAL and MIP energy cutoff
   desc.add<double>("coneRadiusMIP", 14.0);
   desc.add<double>("coneRadiusMIP2", 18.0);
   desc.add<double>("coneRadiusMIP3", 20.0);
   desc.add<double>("coneRadiusMIP4", 22.0);
   desc.add<double>("coneRadiusMIP5", 24.0);
   desc.add<double>("maximumEcalEnergy", 2.0);
   // following 4 parameters are for isolation cuts and described in the code
   desc.add<double>("maxTrackP", 8.0);
   desc.add<double>("slopeTrackP", 0.05090504066);
   desc.add<double>("isolationEnergyTight", 2.0);
   desc.add<double>("isolationEnergyLoose", 10.0);
   // energy thershold for ECAL (from Egamma group)
   desc.add<double>("EBHitEnergyThreshold", 0.08);
   desc.add<double>("EEHitEnergyThreshold0", 0.30);
   desc.add<double>("EEHitEnergyThreshold1", 0.00);
   desc.add<double>("EEHitEnergyThreshold2", 0.00);
   desc.add<double>("EEHitEnergyThreshold3", 0.00);
   desc.add<double>("EEHitEnergyThresholdLow", 0.30);
   desc.add<double>("EEHitEnergyThresholdHigh", 0.30);
   // prescale factors
   desc.add<double>("momentumLow", 40.0);
   desc.add<double>("momentumHigh", 60.0);
   desc.add<int>("prescaleLow", 1);
   desc.add<int>("prescaleHigh", 1);
   // various labels for collections used in the code
   desc.add<edm::InputTag>("labelTriggerEvent", edm::InputTag("hltTriggerSummaryAOD", "", "HLT"));
   desc.add<edm::InputTag>("labelTriggerResult", edm::InputTag("TriggerResults", "", "HLT"));
   desc.add<std::string>("labelTrack", "generalTracks");
   desc.add<std::string>("labelVertex", "offlinePrimaryVertices");
   desc.add<std::string>("labelEBRecHit", "EcalRecHitsEB");
   desc.add<std::string>("labelEERecHit", "EcalRecHitsEE");
   desc.add<std::string>("labelHBHERecHit", "hbhereco");
   desc.add<std::string>("labelBeamSpot", "offlineBeamSpot");
   desc.add<std::string>("labelCaloTower", "towerMaker");
   desc.add<std::string>("labelMuon", "muons");
   desc.add<edm::InputTag>("algInputTag", edm::InputTag("gtStage2Digis"));
   desc.add<edm::InputTag>("extInputTag", edm::InputTag("gtStage2Digis"));
   desc.addUntracked<std::string>("moduleName", "");
   desc.addUntracked<std::string>("producerName", "");
   //  Various flags used for selecting tracks, choice of energy Method2/0
   //  Data type 0/1 for single jet trigger or others
   desc.addUntracked<int>("useRaw", 0);
   desc.addUntracked<bool>("ignoreTriggers", false);
   desc.addUntracked<bool>("useL1Trigger", false);
   desc.addUntracked<double>("hcalScale", 1.0);
   desc.addUntracked<int>("dataType", 0);
   desc.addUntracked<bool>("unCorrect", false);
   desc.addUntracked<bool>("collapseDepth", false);
   desc.addUntracked<std::string>("l1TrigName", "L1_SingleJet60");
   desc.addUntracked<int>("outMode", 11);
   std::vector<int> dummy;
   desc.addUntracked<std::vector<int>>("oldID", dummy);
   desc.addUntracked<std::vector<int>>("newDepth", dummy);
   desc.addUntracked<bool>("hep17", false);
   std::vector<int> events;
   desc.add<std::vector<int>>("debugEvents", events);
   descriptions.add("HcalIsoTrkAnalyzer", desc);
 }

std::array< int, 3 > HcalIsoTrkAnalyzer::fillTree	(	std::vector< math::XYZTLorentzVector > &	vecL1,
		std::vector< math::XYZTLorentzVector > &	vecL3,
		math::XYZPoint &	leadPV,
		std::vector< spr::propagatedTrackDirection > &	trkCaloDirections,
		std::vector< spr::propagatedTrackID > &	trkCaloDets,
		const CaloGeometry *	geo,
		const CaloTopology *	topo,
		const HcalTopology *	theHBHETopology,
		const EcalChannelStatus *	theEcalChStatus,
		const EcalSeverityLevelAlgo *	theEcalSevlv,
		edm::Handle< EcalRecHitCollection > &	barrelRecHitsHandle,
		edm::Handle< EcalRecHitCollection > &	endcapRecHitsHandle,
		edm::Handle< HBHERecHitCollection > &	hbhe,
		edm::Handle< CaloTowerCollection > &	towerHandle,
		edm::Handle< reco::GenParticleCollection > &	genParticles,
		const HcalRespCorrs *	respCorrs,
		const edm::Handle< reco::MuonCollection > &	muonh
	)

private

Definition at line 937 of file HcalIsoTrkAnalyzer.cc.

References a_charIsoR_, a_coneR1_, a_coneR2_, a_coneR_, a_mipR2_, a_mipR3_, a_mipR4_, a_mipR5_, a_mipR_, funct::abs(), accept(), spr::chargeIsolationCone(), spr::chargeIsolationEcal(), debug_, dR(), runTauDisplay::dr, spr::eCone_ecal(), spr::eCone_hcal(), spr::eECALmatrix(), eEcalMax_, spr::eHCALmatrix(), eIsolate1_, eIsolate2_, reco::TrackBase::eta(), eThreshold(), funct::exp(), spr::goodTrack(), HcalDetId::ieta(), HcalDetId::iphi(), edm::HandleBase::isValid(), isotrackApplyRegressor::k, spr::trackSelectionParameters::maxDxyPV, spr::trackSelectionParameters::maxDzPV, spr::trackSelectionParameters::maxInMiss, spr::trackSelectionParameters::maxOutMiss, maxRestrictionP_, newId(), nHigh_, nLow_, notaMuon(), BeamSpotPI::nTracks, oldID_, reco::TrackBase::p(), reco::TrackBase::phi(), prescaleHigh_, prescaleLow_, reco::TrackBase::pt(), pTrackHigh_, pTrackLow_, pTrackMin_, reco::TrackBase::px(), reco::TrackBase::py(), reco::TrackBase::pz(), rhoh(), selectionParameter_, slopeRestrictionP_, storeEnergy(), t_DetIds, t_DetIds1, t_DetIds3, t_eAnnular, t_eHcal, t_eHcal10, t_eHcal30, t_emaxNearP, t_eMipDR, t_eMipDR2, t_eMipDR3, t_eMipDR4, t_eMipDR5, t_EventNo, t_gentrackP, t_hAnnular, t_HitEnergies, t_HitEnergies1, t_HitEnergies3, t_hmaxNearP, t_ieta, t_ietaAll, t_ietaGood, t_iphi, t_mindR1, t_mindR2, t_nTrk, t_p, t_phi, t_pt, t_qltyFlag, t_qltyMissFlag, t_qltyPVFlag, t_rhoh, t_RunNo, t_selectTk, t_trackType, t_trgbits, trackP(), tree, and useRaw_.

Referenced by analyze(), and core.AutoFillTreeProducer.AutoFillTreeProducer::process().

                                                                                             {
   int nSave(0), nLoose(0), nTight(0);
   //Loop over tracks
   std::vector<spr::propagatedTrackDirection>::const_iterator trkDetItr;
   unsigned int nTracks(0), nselTracks(0);
   t_nTrk = trkCaloDirections.size();
   t_rhoh = (tower.isValid()) ? rhoh(tower) : 0;
   for (trkDetItr = trkCaloDirections.begin(), nTracks = 0; trkDetItr != trkCaloDirections.end();
        trkDetItr++, nTracks++) {
     const reco::Track* pTrack = &(*(trkDetItr->trkItr));
     math::XYZTLorentzVector v4(pTrack->px(), pTrack->py(), pTrack->pz(), pTrack->p());
     t_p = pTrack->p();
     t_pt = pTrack->pt();
     t_phi = pTrack->phi();
 
 #ifdef EDM_ML_DEBUG
     if (debug_)
       edm::LogVerbatim("HcalIsoTrack") << "This track : " << nTracks << " (pt|eta|phi|p) : " << t_pt << "|"
                                        << pTrack->eta() << "|" << t_phi << "|" << t_p;
 #endif
     t_mindR2 = 999;
     for (unsigned int k = 0; k < vecL3.size(); ++k) {
       double dr = dR(vecL3[k], v4);
       if (dr < t_mindR2) {
         t_mindR2 = dr;
       }
     }
     t_mindR1 = (!vecL1.empty()) ? dR(vecL1[0], v4) : 999;
 #ifdef EDM_ML_DEBUG
     if (debug_)
       edm::LogVerbatim("HcalIsoTrack") << "Closest L3 object at dr : " << t_mindR2 << " and from L1 " << t_mindR1;
 #endif
 
     t_ieta = t_iphi = 0;
     if (trkDetItr->okHCAL) {
       HcalDetId detId = (HcalDetId)(trkDetItr->detIdHCAL);
       t_ieta = detId.ieta();
       t_iphi = detId.iphi();
       if (t_p > 40.0 && t_p <= 60.0)
         t_ietaAll->emplace_back(t_ieta);
     }
     //Selection of good track
     t_selectTk = spr::goodTrack(pTrack, leadPV, selectionParameter_, false);
     spr::trackSelectionParameters oneCutParameters = selectionParameter_;
     oneCutParameters.maxDxyPV = 10;
     oneCutParameters.maxDzPV = 100;
     oneCutParameters.maxInMiss = 2;
     oneCutParameters.maxOutMiss = 2;
     bool qltyFlag = spr::goodTrack(pTrack, leadPV, oneCutParameters, false);
     oneCutParameters = selectionParameter_;
     oneCutParameters.maxDxyPV = 10;
     oneCutParameters.maxDzPV = 100;
     t_qltyMissFlag = spr::goodTrack(pTrack, leadPV, oneCutParameters, false);
     oneCutParameters = selectionParameter_;
     oneCutParameters.maxInMiss = 2;
     oneCutParameters.maxOutMiss = 2;
     t_qltyPVFlag = spr::goodTrack(pTrack, leadPV, oneCutParameters, false);
     double eIsolation = maxRestrictionP_ * exp(slopeRestrictionP_ * std::abs((double)t_ieta));
     if (eIsolation < eIsolate1_)
       eIsolation = eIsolate1_;
     if (eIsolation < eIsolate2_)
       eIsolation = eIsolate2_;
 #ifdef EDM_ML_DEBUG
     if (debug_)
       edm::LogVerbatim("HcalIsoTrack") << "qltyFlag|okECAL|okHCAL : " << qltyFlag << "|" << trkDetItr->okECAL << "|"
                                        << trkDetItr->okHCAL << " eIsolation " << eIsolation;
 #endif
     t_qltyFlag = (qltyFlag && trkDetItr->okECAL && trkDetItr->okHCAL);
     bool notMuon = (muonh.isValid()) ? notaMuon(pTrack, muonh) : true;
     if (t_qltyFlag && notMuon) {
       nselTracks++;
       int nNearTRKs(0);
       std::vector<DetId> eIds;
       std::vector<double> eHit;
       t_eMipDR = spr::eCone_ecal(geo,
                                  barrelRecHitsHandle,
                                  endcapRecHitsHandle,
                                  trkDetItr->pointHCAL,
                                  trkDetItr->pointECAL,
                                  a_mipR_,
                                  trkDetItr->directionECAL,
                                  eIds,
                                  eHit);
       double eEcal(0);
       for (unsigned int k = 0; k < eIds.size(); ++k) {
         if (eHit[k] > eThreshold(eIds[k], geo))
           eEcal += eHit[k];
       }
 #ifdef EDM_ML_DEBUG
       if (debug_)
         edm::LogVerbatim("HcalIsoTrack") << "eMIP before and after: " << t_eMipDR << ":" << eEcal;
 #endif
       t_eMipDR = eEcal;
       std::vector<DetId> eIds2;
       std::vector<double> eHit2;
       t_eMipDR2 = spr::eCone_ecal(geo,
                                   barrelRecHitsHandle,
                                   endcapRecHitsHandle,
                                   trkDetItr->pointHCAL,
                                   trkDetItr->pointECAL,
                                   a_mipR2_,
                                   trkDetItr->directionECAL,
                                   eIds2,
                                   eHit2);
       double eEcal2(0);
       for (unsigned int k = 0; k < eIds2.size(); ++k) {
         if (eHit2[k] > eThreshold(eIds2[k], geo))
           eEcal2 += eHit2[k];
       }
 #ifdef EDM_ML_DEBUG
       if (debug_)
         edm::LogVerbatim("HcalIsoTrack") << "eMIP before and after: " << t_eMipDR2 << ":" << eEcal2;
 #endif
       t_eMipDR2 = eEcal2;
       std::vector<DetId> eIds3;
       std::vector<double> eHit3;
       t_eMipDR3 = spr::eCone_ecal(geo,
                                   barrelRecHitsHandle,
                                   endcapRecHitsHandle,
                                   trkDetItr->pointHCAL,
                                   trkDetItr->pointECAL,
                                   a_mipR3_,
                                   trkDetItr->directionECAL,
                                   eIds3,
                                   eHit3);
       double eEcal3(0);
       for (unsigned int k = 0; k < eIds3.size(); ++k) {
         if (eHit3[k] > eThreshold(eIds3[k], geo))
           eEcal3 += eHit3[k];
       }
 #ifdef EDM_ML_DEBUG
       if (debug_)
         edm::LogVerbatim("HcalIsoTrack") << "eMIP before and after: " << t_eMipDR3 << ":" << eEcal3;
 #endif
       t_eMipDR3 = eEcal3;
       std::vector<DetId> eIds4;
       std::vector<double> eHit4;
       t_eMipDR4 = spr::eCone_ecal(geo,
                                   barrelRecHitsHandle,
                                   endcapRecHitsHandle,
                                   trkDetItr->pointHCAL,
                                   trkDetItr->pointECAL,
                                   a_mipR4_,
                                   trkDetItr->directionECAL,
                                   eIds4,
                                   eHit4);
       double eEcal4(0);
       for (unsigned int k = 0; k < eIds4.size(); ++k) {
         if (eHit4[k] > eThreshold(eIds4[k], geo))
           eEcal4 += eHit4[k];
       }
 #ifdef EDM_ML_DEBUG
       if (debug_)
         edm::LogVerbatim("HcalIsoTrack") << "eMIP before and after: " << t_eMipDR4 << ":" << eEcal4;
 #endif
       t_eMipDR4 = eEcal4;
       std::vector<DetId> eIds5;
       std::vector<double> eHit5;
       t_eMipDR5 = spr::eCone_ecal(geo,
                                   barrelRecHitsHandle,
                                   endcapRecHitsHandle,
                                   trkDetItr->pointHCAL,
                                   trkDetItr->pointECAL,
                                   a_mipR5_,
                                   trkDetItr->directionECAL,
                                   eIds5,
                                   eHit5);
       double eEcal5(0);
       for (unsigned int k = 0; k < eIds5.size(); ++k) {
         if (eHit5[k] > eThreshold(eIds5[k], geo))
           eEcal5 += eHit5[k];
       }
 #ifdef EDM_ML_DEBUG
       if (debug_)
         edm::LogVerbatim("HcalIsoTrack") << "eMIP before and after: " << t_eMipDR5 << ":" << eEcal5;
 #endif
       t_eMipDR5 = eEcal5;
 
       t_emaxNearP = spr::chargeIsolationEcal(nTracks, trkCaloDets, geo, caloTopology, 15, 15);
       const DetId cellE(trkDetItr->detIdECAL);
       std::pair<double, bool> e11x11P = spr::eECALmatrix(cellE,
                                                          barrelRecHitsHandle,
                                                          endcapRecHitsHandle,
                                                          *theEcalChStatus,
                                                          geo,
                                                          caloTopology,
                                                          theEcalSevlv,
                                                          5,
                                                          5,
                                                          -100.0,
                                                          -100.0,
                                                          -100.0,
                                                          100.0);
       std::pair<double, bool> e15x15P = spr::eECALmatrix(cellE,
                                                          barrelRecHitsHandle,
                                                          endcapRecHitsHandle,
                                                          *theEcalChStatus,
                                                          geo,
                                                          caloTopology,
                                                          theEcalSevlv,
                                                          7,
                                                          7,
                                                          -100.0,
                                                          -100.0,
                                                          -100.0,
                                                          100.0);
       if (e11x11P.second && e15x15P.second) {
         t_eAnnular = (e15x15P.first - e11x11P.first);
       } else {
         t_eAnnular = -(e15x15P.first - e11x11P.first);
       }
       t_hmaxNearP = spr::chargeIsolationCone(nTracks, trkCaloDirections, a_charIsoR_, nNearTRKs, false);
       const DetId cellH(trkDetItr->detIdHCAL);
       double h5x5 = spr::eHCALmatrix(
           theHBHETopology, cellH, hbhe, 2, 2, false, true, -100.0, -100.0, -100.0, -100.0, -100.0, 100.0);
       double h7x7 = spr::eHCALmatrix(
           theHBHETopology, cellH, hbhe, 3, 3, false, true, -100.0, -100.0, -100.0, -100.0, -100.0, 100.0);
       t_hAnnular = h7x7 - h5x5;
 #ifdef EDM_ML_DEBUG
       if (debug_)
         edm::LogVerbatim("HcalIsoTrack") << "max p Near (Ecal) " << t_emaxNearP << " (Hcal) " << t_hmaxNearP
                                          << " Annular E (Ecal) " << e11x11P.first << ":" << e15x15P.first << ":"
                                          << t_eAnnular << " (Hcal) " << h5x5 << ":" << h7x7 << ":" << t_hAnnular;
 #endif
       t_gentrackP = trackP(pTrack, genParticles);
       if (t_eMipDR < eEcalMax_ && t_hmaxNearP < eIsolation) {
         t_DetIds->clear();
         t_HitEnergies->clear();
         t_DetIds1->clear();
         t_HitEnergies1->clear();
         t_DetIds3->clear();
         t_HitEnergies3->clear();
         int nRecHits(-999), nRecHits1(-999), nRecHits3(-999);
         std::vector<DetId> ids, ids1, ids3;
         std::vector<double> edet0, edet1, edet3;
         t_eHcal = spr::eCone_hcal(geo,
                                   hbhe,
                                   trkDetItr->pointHCAL,
                                   trkDetItr->pointECAL,
                                   a_coneR_,
                                   trkDetItr->directionHCAL,
                                   nRecHits,
                                   ids,
                                   edet0,
                                   useRaw_);
         if (!oldID_.empty()) {
           for (unsigned k = 0; k < ids.size(); ++k)
             ids[k] = newId(ids[k]);
         }
         storeEnergy(0, respCorrs, ids, edet0, t_eHcal, t_DetIds, t_HitEnergies);
 
         //----- hcal energy in the extended cone 1 (a_coneR+10) --------------
         t_eHcal10 = spr::eCone_hcal(geo,
                                     hbhe,
                                     trkDetItr->pointHCAL,
                                     trkDetItr->pointECAL,
                                     a_coneR1_,
                                     trkDetItr->directionHCAL,
                                     nRecHits1,
                                     ids1,
                                     edet1,
                                     useRaw_);
         if (!oldID_.empty()) {
           for (unsigned k = 0; k < ids1.size(); ++k)
             ids1[k] = newId(ids1[k]);
         }
         storeEnergy(1, respCorrs, ids1, edet1, t_eHcal10, t_DetIds1, t_HitEnergies1);
 
         //----- hcal energy in the extended cone 3 (a_coneR+30) --------------
         t_eHcal30 = spr::eCone_hcal(geo,
                                     hbhe,
                                     trkDetItr->pointHCAL,
                                     trkDetItr->pointECAL,
                                     a_coneR2_,
                                     trkDetItr->directionHCAL,
                                     nRecHits3,
                                     ids3,
                                     edet3,
                                     useRaw_);
         if (!oldID_.empty()) {
           for (unsigned k = 0; k < ids3.size(); ++k)
             ids3[k] = newId(ids3[k]);
         }
         storeEnergy(3, respCorrs, ids3, edet3, t_eHcal30, t_DetIds3, t_HitEnergies3);
 
 #ifdef EDM_ML_DEBUG
         if (debug_) {
           edm::LogVerbatim("HcalIsoTrack")
               << "This track : " << nTracks << " (pt|eta|phi|p) : " << t_pt << "|" << pTrack->eta() << "|" << t_phi
               << "|" << t_p << " Generator Level p " << t_gentrackP;
           edm::LogVerbatim("HcalIsoTrack")
               << "e_MIP " << t_eMipDR << " Chg Isolation " << t_hmaxNearP << " eHcal" << t_eHcal << " ieta " << t_ieta
               << " Quality " << t_qltyMissFlag << ":" << t_qltyPVFlag << ":" << t_selectTk;
           for (unsigned int ll = 0; ll < t_DetIds->size(); ll++) {
             edm::LogVerbatim("HcalIsoTrack")
                 << "det id is = " << t_DetIds->at(ll) << "   hit enery is  = " << t_HitEnergies->at(ll);
           }
           for (unsigned int ll = 0; ll < t_DetIds1->size(); ll++) {
             edm::LogVerbatim("HcalIsoTrack")
                 << "det id is = " << t_DetIds1->at(ll) << "   hit enery is  = " << t_HitEnergies1->at(ll);
           }
           for (unsigned int ll = 0; ll < t_DetIds3->size(); ll++) {
             edm::LogVerbatim("HcalIsoTrack")
                 << "det id is = " << t_DetIds3->at(ll) << "   hit enery is  = " << t_HitEnergies3->at(ll);
           }
         }
 #endif
         bool accept(false);
         if (t_p > pTrackMin_) {
           if (t_p < pTrackLow_) {
             ++nLow_;
             if (prescaleLow_ <= 1)
               accept = true;
             else if (nLow_ % prescaleLow_ == 1)
               accept = true;
           } else if (t_p > pTrackHigh_) {
             ++nHigh_;
             if (prescaleHigh_ <= 1)
               accept = true;
             else if (nHigh_ % prescaleHigh_ == 1)
               accept = true;
           } else {
             accept = true;
           }
         }
         if (accept) {
           tree->Fill();
           edm::LogVerbatim("HcalIsoTrackX") << "Run " << t_RunNo << " Event " << t_EventNo;
           nSave++;
           int type(0);
           if (t_eMipDR < 1.0) {
             if (t_hmaxNearP < eIsolate2_) {
               ++nLoose;
               type = 1;
             }
             if (t_hmaxNearP < eIsolate1_) {
               ++nTight;
               type = 2;
             }
           }
           if (t_p > 40.0 && t_p <= 60.0 && t_selectTk) {
             t_ietaGood->emplace_back(t_ieta);
             t_trackType->emplace_back(type);
           }
 #ifdef EDM_ML_DEBUG
           if (debug_) {
             for (unsigned int k = 0; k < t_trgbits->size(); k++) {
               edm::LogVerbatim("HcalIsoTrack") << "trigger bit is  = " << t_trgbits->at(k);
             }
           }
 #endif
         }
       }
     }
   }
   std::array<int, 3> i3{{nSave, nLoose, nTight}};
   return i3;
 }

DetId HcalIsoTrkAnalyzer::newId ( const DetId & id )

private

Definition at line 1388 of file HcalIsoTrkAnalyzer.cc.

References HcalDetId::depth(), hep17_, gpuClustering::id, HcalDetId::ieta(), HcalDetId::ietaAbs(), HcalDetId::iphi(), isotrackApplyRegressor::k, newDepth_, oldDepth_, oldDet_, oldEta_, oldID_, HcalDetId::subdet(), DetId::subdetId(), and HcalDetId::zside().

Referenced by fillTree().

                                                {
   HcalDetId hid(id);
   if (hep17_ && ((hid.iphi() < 63) || (hid.iphi() > 66) || (hid.zside() < 0)))
     return id;
   for (unsigned int k = 0; k < oldID_.size(); ++k) {
     if ((hid.subdetId() == oldDet_[k]) && (hid.ietaAbs() == oldEta_[k]) && (hid.depth() == oldDepth_[k])) {
       return static_cast<DetId>(HcalDetId(hid.subdet(), hid.ieta(), hid.iphi(), newDepth_[k]));
     }
   }
   return id;
 }

bool HcalIsoTrkAnalyzer::notaMuon	(	const reco::Track *	pTrack0,
		const edm::Handle< reco::MuonCollection > &	muonh
	)

private

Definition at line 1461 of file HcalIsoTrkAnalyzer.cc.

References reco::deltaR(), dR(), reco::TrackBase::eta(), reco::TrackBase::phi(), and muon::segmentCompatibility().

Referenced by fillTree().

                                                                                                       {
   bool flag(true);
   for (reco::MuonCollection::const_iterator recMuon = muonh->begin(); recMuon != muonh->end(); ++recMuon) {
     if (recMuon->innerTrack().isNonnull()) {
       const reco::Track* pTrack = (recMuon->innerTrack()).get();
       bool mediumMuon = (((recMuon->isPFMuon()) && (recMuon->isGlobalMuon() || recMuon->isTrackerMuon())) &&
                          (recMuon->innerTrack()->validFraction() > 0.49));
       if (mediumMuon) {
         double chiGlobal = ((recMuon->globalTrack().isNonnull()) ? recMuon->globalTrack()->normalizedChi2() : 999);
         bool goodGlob = (recMuon->isGlobalMuon() && chiGlobal < 3 &&
                          recMuon->combinedQuality().chi2LocalPosition < 12 && recMuon->combinedQuality().trkKink < 20);
         mediumMuon = muon::segmentCompatibility(*recMuon) > (goodGlob ? 0.303 : 0.451);
       }
       if (mediumMuon) {
         double dR = reco::deltaR(pTrack->eta(), pTrack->phi(), pTrack0->eta(), pTrack0->phi());
         if (dR < 0.1) {
           flag = false;
           break;
         }
       }
     }
   }
   return flag;
 }

double HcalIsoTrkAnalyzer::rhoh ( const edm::Handle< CaloTowerCollection > & tower )

private

Definition at line 1342 of file HcalIsoTrkAnalyzer.cc.

References debug_, reco::deltaPhi(), PVValHelper::eta, etabins_, etadist_, etahalfdist_, phi, phibins_, phidist_, and phihalfdist_.

Referenced by fillTree().

                                                                            {
   std::vector<double> sumPFNallSMDQH2;
   sumPFNallSMDQH2.reserve(phibins_.size() * etabins_.size());
 
   for (auto eta : etabins_) {
     for (auto phi : phibins_) {
       double hadder = 0;
       for (const auto& pf_it : (*tower)) {
         if (fabs(eta - pf_it.eta()) > etahalfdist_)
           continue;
         if (fabs(reco::deltaPhi(phi, pf_it.phi())) > phihalfdist_)
           continue;
         hadder += pf_it.hadEt();
       }
       sumPFNallSMDQH2.emplace_back(hadder);
     }
   }
 
   double evt_smdq(0);
   std::sort(sumPFNallSMDQH2.begin(), sumPFNallSMDQH2.end());
   if (sumPFNallSMDQH2.size() % 2)
     evt_smdq = sumPFNallSMDQH2[(sumPFNallSMDQH2.size() - 1) / 2];
   else
     evt_smdq = (sumPFNallSMDQH2[sumPFNallSMDQH2.size() / 2] + sumPFNallSMDQH2[(sumPFNallSMDQH2.size() - 2) / 2]) / 2.;
   double rhoh = evt_smdq / (etadist_ * phidist_);
 #ifdef EDM_ML_DEBUG
   if (debug_)
     edm::LogVerbatim("HcalIsoTrack") << "Rho " << evt_smdq << ":" << rhoh;
 #endif
   return rhoh;
 }

void HcalIsoTrkAnalyzer::storeEnergy	(	int	indx,
		const HcalRespCorrs *	respCorrs,
		const std::vector< DetId > &	ids,
		std::vector< double > &	edet,
		double &	eHcal,
		std::vector< unsigned int > *	detIds,
		std::vector< double > *	hitEnergies
	)

private

Definition at line 1400 of file HcalIsoTrkAnalyzer.cc.

References funct::abs(), collapseDepth_, alignCSCRings::corr, debug_, reco::JetExtendedAssociation::getValue(), HcalCondObjectContainer< Item >::getValues(), hcalScale_, hdc_, gpuClustering::id, isotrackApplyRegressor::k, HcalDDDRecConstants::mergedDepthDetId(), and unCorrect_.

Referenced by fillTree().

                                                                      {
   double ehcal(0);
   if (unCorrect_) {
     for (unsigned int k = 0; k < ids.size(); ++k) {
       double corr = (respCorrs->getValues(ids[k]))->getValue();
       if (corr != 0)
         edet[k] /= corr;
       ehcal += edet[k];
     }
   } else {
     for (const auto& en : edet)
       ehcal += en;
   }
   if ((std::abs(ehcal - eHcal) > 0.001) && (!unCorrect_))
     edm::LogWarning("HcalIsoTrack") << "Check inconsistent energies: " << indx << " " << eHcal << ":" << ehcal
                                     << " from " << ids.size() << " cells";
   eHcal = hcalScale_ * ehcal;
 
   if (collapseDepth_) {
     std::map<HcalDetId, double> hitMap;
     for (unsigned int k = 0; k < ids.size(); ++k) {
       HcalDetId id = hdc_->mergedDepthDetId(HcalDetId(ids[k]));
       auto itr = hitMap.find(id);
       if (itr == hitMap.end()) {
         hitMap[id] = edet[k];
       } else {
         (itr->second) += edet[k];
       }
     }
     detIds->reserve(hitMap.size());
     hitEnergies->reserve(hitMap.size());
     for (const auto& hit : hitMap) {
       detIds->emplace_back(hit.first.rawId());
       hitEnergies->emplace_back(hit.second);
     }
   } else {
     detIds->reserve(ids.size());
     hitEnergies->reserve(ids.size());
     for (unsigned int k = 0; k < ids.size(); ++k) {
       detIds->emplace_back(ids[k].rawId());
       hitEnergies->emplace_back(edet[k]);
     }
   }
 #ifdef EDM_ML_DEBUG
   if (debug_) {
     edm::LogVerbatim("HcalIsoTrack") << "Input to storeEnergy with " << ids.size() << " cells";
     for (unsigned int k = 0; k < ids.size(); ++k)
       edm::LogVerbatim("HcalIsoTrack") << "Hit [" << k << "] " << HcalDetId(ids[k]) << " E " << edet[k];
     edm::LogVerbatim("HcalIsoTrack") << "Output of storeEnergy with " << detIds->size() << " cells and Etot " << eHcal;
     for (unsigned int k = 0; k < detIds->size(); ++k)
       edm::LogVerbatim("HcalIsoTrack") << "Hit [" << k << "] " << HcalDetId((*detIds)[k]) << " E " << (*hitEnergies)[k];
   }
 #endif
 }

double HcalIsoTrkAnalyzer::trackP	(	const reco::Track *	pTrack,
		const edm::Handle< reco::GenParticleCollection > &	genParticles
	)

private

Definition at line 1326 of file HcalIsoTrkAnalyzer.cc.

References reco::deltaR(), dR(), reco::TrackBase::eta(), edm::HandleBase::isValid(), AlCaHLTBitMon_ParallelJobs::p, reco::TrackBase::phi(), and isotrackTrainRegressor::pmom.

Referenced by fillTree().

                                                                                             {
   double pmom = -1.0;
   if (genParticles.isValid()) {
     double mindR(999.9);
     for (const auto& p : (*genParticles)) {
       double dR = reco::deltaR(pTrack->eta(), pTrack->phi(), p.momentum().Eta(), p.momentum().Phi());
       if (dR < mindR) {
         mindR = dR;
         pmom = p.momentum().R();
       }
     }
   }
   return pmom;
 }

Member Data Documentation

double HcalIsoTrkAnalyzer::a_charIsoR_

private

Definition at line 162 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

double HcalIsoTrkAnalyzer::a_coneR1_

private

Definition at line 162 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

double HcalIsoTrkAnalyzer::a_coneR2_

private

Definition at line 162 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

const double HcalIsoTrkAnalyzer::a_coneR_

private

Definition at line 143 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

const double HcalIsoTrkAnalyzer::a_mipR2_

private

Definition at line 143 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

const double HcalIsoTrkAnalyzer::a_mipR3_

private

Definition at line 143 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

const double HcalIsoTrkAnalyzer::a_mipR4_

private

Definition at line 144 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

const double HcalIsoTrkAnalyzer::a_mipR5_

private

Definition at line 144 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

const double HcalIsoTrkAnalyzer::a_mipR_

private

Definition at line 143 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

const bool HcalIsoTrkAnalyzer::collapseDepth_

private

Definition at line 151 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer(), and storeEnergy().

const int HcalIsoTrkAnalyzer::dataType_

private

Definition at line 149 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

const std::vector<int> HcalIsoTrkAnalyzer::debEvents_

private

Definition at line 160 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

bool HcalIsoTrkAnalyzer::debug_

private

Definition at line 212 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), fillTree(), rhoh(), and storeEnergy().

const double HcalIsoTrkAnalyzer::eEcalMax_

private

Definition at line 144 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

const double HcalIsoTrkAnalyzer::eIsolate1_

private

Definition at line 146 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

const double HcalIsoTrkAnalyzer::eIsolate2_

private

Definition at line 146 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

std::vector<double> HcalIsoTrkAnalyzer::etabins_

private

Definition at line 164 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer(), and rhoh().

double HcalIsoTrkAnalyzer::etadist_

private

Definition at line 166 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer(), and rhoh().

double HcalIsoTrkAnalyzer::etahalfdist_

private

Definition at line 166 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer(), and rhoh().

edm::Service<TFileService> HcalIsoTrkAnalyzer::fs

private

Definition at line 136 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob().

const double HcalIsoTrkAnalyzer::hcalScale_

private

Definition at line 146 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer(), and storeEnergy().

const HcalDDDRecConstants* HcalIsoTrkAnalyzer::hdc_

private

Definition at line 163 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginRun(), and storeEnergy().

const bool HcalIsoTrkAnalyzer::hep17_

private

Definition at line 159 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer(), and newId().

const double HcalIsoTrkAnalyzer::hitEthrEB_

private

Definition at line 152 of file HcalIsoTrkAnalyzer.cc.

Referenced by eThreshold(), and HcalIsoTrkAnalyzer().

const double HcalIsoTrkAnalyzer::hitEthrEE0_

private

Definition at line 152 of file HcalIsoTrkAnalyzer.cc.

Referenced by eThreshold(), and HcalIsoTrkAnalyzer().

const double HcalIsoTrkAnalyzer::hitEthrEE1_

private

Definition at line 152 of file HcalIsoTrkAnalyzer.cc.

Referenced by eThreshold(), and HcalIsoTrkAnalyzer().

const double HcalIsoTrkAnalyzer::hitEthrEE2_

private

Definition at line 153 of file HcalIsoTrkAnalyzer.cc.

Referenced by eThreshold(), and HcalIsoTrkAnalyzer().

const double HcalIsoTrkAnalyzer::hitEthrEE3_

private

Definition at line 153 of file HcalIsoTrkAnalyzer.cc.

Referenced by eThreshold(), and HcalIsoTrkAnalyzer().

const double HcalIsoTrkAnalyzer::hitEthrEEHi_

private

Definition at line 154 of file HcalIsoTrkAnalyzer.cc.

Referenced by eThreshold(), and HcalIsoTrkAnalyzer().

const double HcalIsoTrkAnalyzer::hitEthrEELo_

private

Definition at line 154 of file HcalIsoTrkAnalyzer.cc.

Referenced by eThreshold(), and HcalIsoTrkAnalyzer().

HLTConfigProvider HcalIsoTrkAnalyzer::hltConfig_

private

Definition at line 137 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginRun().

const bool HcalIsoTrkAnalyzer::ignoreTrigger_

private

Definition at line 150 of file HcalIsoTrkAnalyzer.cc.

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

const std::string HcalIsoTrkAnalyzer::l1Filter_

private

Definition at line 141 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

l1t::L1TGlobalUtil* HcalIsoTrkAnalyzer::l1GtUtils_

private

Definition at line 135 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

const std::string HcalIsoTrkAnalyzer::l1TrigName_

private

Definition at line 157 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

const std::string HcalIsoTrkAnalyzer::l2Filter_

private

Definition at line 142 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

const std::string HcalIsoTrkAnalyzer::l3Filter_

private

Definition at line 142 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

const std::string HcalIsoTrkAnalyzer::labelEB_

private

Definition at line 156 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

const std::string HcalIsoTrkAnalyzer::labelEE_

private

Definition at line 157 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

const std::string HcalIsoTrkAnalyzer::labelGenTrack_

private

Definition at line 156 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

const std::string HcalIsoTrkAnalyzer::labelHBHE_

private

Definition at line 157 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

const std::string HcalIsoTrkAnalyzer::labelRecVtx_

private

Definition at line 156 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer().

const std::string HcalIsoTrkAnalyzer::labelTower_

private

Definition at line 157 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer().

const double HcalIsoTrkAnalyzer::maxRestrictionP_

private

Definition at line 145 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

const int HcalIsoTrkAnalyzer::mode_

private

Definition at line 149 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer().

const std::vector<int> HcalIsoTrkAnalyzer::newDepth_

private

Definition at line 158 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer(), and newId().

unsigned int HcalIsoTrkAnalyzer::nHigh_

private

Definition at line 161 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree().

unsigned int HcalIsoTrkAnalyzer::nLow_

private

Definition at line 161 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree().

unsigned int HcalIsoTrkAnalyzer::nRun_

private

Definition at line 161 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginRun(), and endRun().

std::vector<int> HcalIsoTrkAnalyzer::oldDepth_

private

Definition at line 165 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer(), and newId().

std::vector<int> HcalIsoTrkAnalyzer::oldDet_

private

Definition at line 165 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer(), and newId().

std::vector<int> HcalIsoTrkAnalyzer::oldEta_

private

Definition at line 165 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer(), and newId().

const std::vector<int> HcalIsoTrkAnalyzer::oldID_

private

Definition at line 158 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), HcalIsoTrkAnalyzer(), and newId().

std::vector<double> HcalIsoTrkAnalyzer::phibins_

private

Definition at line 164 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer(), and rhoh().

double HcalIsoTrkAnalyzer::phidist_

private

Definition at line 166 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer(), and rhoh().

double HcalIsoTrkAnalyzer::phihalfdist_

private

Definition at line 166 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer(), and rhoh().

const int HcalIsoTrkAnalyzer::prescaleHigh_

private

Definition at line 148 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

const int HcalIsoTrkAnalyzer::prescaleLow_

private

Definition at line 148 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

const std::string HcalIsoTrkAnalyzer::processName_

private

Definition at line 141 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginRun(), and HcalIsoTrkAnalyzer().

const double HcalIsoTrkAnalyzer::pTrackHigh_

private

Definition at line 147 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

const double HcalIsoTrkAnalyzer::pTrackLow_

private

Definition at line 147 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

const double HcalIsoTrkAnalyzer::pTrackMin_

private

Definition at line 144 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

spr::trackSelectionParameters HcalIsoTrkAnalyzer::selectionParameter_

private

Definition at line 139 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

const double HcalIsoTrkAnalyzer::slopeRestrictionP_

private

Definition at line 145 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

int HcalIsoTrkAnalyzer::t_allvertex

private

Definition at line 209 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

int HcalIsoTrkAnalyzer::t_DataType

private

Definition at line 192 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

std::vector<unsigned int>* HcalIsoTrkAnalyzer::t_DetIds

private

Definition at line 204 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

std::vector<unsigned int> * HcalIsoTrkAnalyzer::t_DetIds1

private

Definition at line 204 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

std::vector<unsigned int> * HcalIsoTrkAnalyzer::t_DetIds3

private

Definition at line 204 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

double HcalIsoTrkAnalyzer::t_eAnnular

private

Definition at line 200 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

double HcalIsoTrkAnalyzer::t_eHcal

private

Definition at line 201 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

double HcalIsoTrkAnalyzer::t_eHcal10

private

Definition at line 201 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

double HcalIsoTrkAnalyzer::t_eHcal30

private

Definition at line 201 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

double HcalIsoTrkAnalyzer::t_emaxNearP

private

Definition at line 200 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

double HcalIsoTrkAnalyzer::t_eMipDR

private

Definition at line 198 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

double HcalIsoTrkAnalyzer::t_eMipDR2

private

Definition at line 198 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

double HcalIsoTrkAnalyzer::t_eMipDR3

private

Definition at line 198 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

double HcalIsoTrkAnalyzer::t_eMipDR4

private

Definition at line 198 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

double HcalIsoTrkAnalyzer::t_eMipDR5

private

Definition at line 199 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

int HcalIsoTrkAnalyzer::t_Event

private

Definition at line 192 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

unsigned int HcalIsoTrkAnalyzer::t_EventNo

private

Definition at line 191 of file HcalIsoTrkAnalyzer.cc.

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

double HcalIsoTrkAnalyzer::t_EventWeight

private

Definition at line 194 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

double HcalIsoTrkAnalyzer::t_gentrackP

private

Definition at line 199 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

int HcalIsoTrkAnalyzer::t_goodPV

private

Definition at line 193 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

double HcalIsoTrkAnalyzer::t_hAnnular

private

Definition at line 200 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

std::vector<double>* HcalIsoTrkAnalyzer::t_HitEnergies

private

Definition at line 205 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

std::vector<double> * HcalIsoTrkAnalyzer::t_HitEnergies1

private

Definition at line 205 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

std::vector<double> * HcalIsoTrkAnalyzer::t_HitEnergies3

private

Definition at line 205 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

std::vector<bool> * HcalIsoTrkAnalyzer::t_hltbits

private

Definition at line 206 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

double HcalIsoTrkAnalyzer::t_hmaxNearP

private

Definition at line 199 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

int HcalIsoTrkAnalyzer::t_ieta

private

Definition at line 192 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

std::vector<int>* HcalIsoTrkAnalyzer::t_ietaAll

private

Definition at line 210 of file HcalIsoTrkAnalyzer.cc.

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

std::vector<int> * HcalIsoTrkAnalyzer::t_ietaGood

private

Definition at line 210 of file HcalIsoTrkAnalyzer.cc.

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

int HcalIsoTrkAnalyzer::t_iphi

private

Definition at line 192 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

bool HcalIsoTrkAnalyzer::t_L1Bit

private

Definition at line 207 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

double HcalIsoTrkAnalyzer::t_l1eta

private

Definition at line 195 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

double HcalIsoTrkAnalyzer::t_l1phi

private

Definition at line 195 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

double HcalIsoTrkAnalyzer::t_l1pt

private

Definition at line 195 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

double HcalIsoTrkAnalyzer::t_l3eta

private

Definition at line 196 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

double HcalIsoTrkAnalyzer::t_l3phi

private

Definition at line 196 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

double HcalIsoTrkAnalyzer::t_l3pt

private

Definition at line 196 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

double HcalIsoTrkAnalyzer::t_mindR1

private

Definition at line 197 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

double HcalIsoTrkAnalyzer::t_mindR2

private

Definition at line 197 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

int HcalIsoTrkAnalyzer::t_nTrk

private

Definition at line 193 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

int HcalIsoTrkAnalyzer::t_nVtx

private

Definition at line 193 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

double HcalIsoTrkAnalyzer::t_p

private

Definition at line 194 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

double HcalIsoTrkAnalyzer::t_phi

private

Definition at line 194 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

double HcalIsoTrkAnalyzer::t_pt

private

Definition at line 194 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

bool HcalIsoTrkAnalyzer::t_qltyFlag

private

Definition at line 202 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

bool HcalIsoTrkAnalyzer::t_qltyMissFlag

private

Definition at line 202 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

bool HcalIsoTrkAnalyzer::t_qltyPVFlag

private

Definition at line 203 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

double HcalIsoTrkAnalyzer::t_rhoh

private

Definition at line 201 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

int HcalIsoTrkAnalyzer::t_Run

private

Definition at line 192 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

unsigned int HcalIsoTrkAnalyzer::t_RunNo

private

Definition at line 191 of file HcalIsoTrkAnalyzer.cc.

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

bool HcalIsoTrkAnalyzer::t_selectTk

private

Definition at line 202 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginJob(), and fillTree().

int HcalIsoTrkAnalyzer::t_Tracks

private

Definition at line 208 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

int HcalIsoTrkAnalyzer::t_TracksLoose

private

Definition at line 209 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

int HcalIsoTrkAnalyzer::t_TracksProp

private

Definition at line 208 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

int HcalIsoTrkAnalyzer::t_TracksSaved

private

Definition at line 208 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

int HcalIsoTrkAnalyzer::t_TracksTight

private

Definition at line 209 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

std::vector<int> * HcalIsoTrkAnalyzer::t_trackType

private

Definition at line 210 of file HcalIsoTrkAnalyzer.cc.

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

std::vector<bool>* HcalIsoTrkAnalyzer::t_trgbits

private

Definition at line 206 of file HcalIsoTrkAnalyzer.cc.

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

bool HcalIsoTrkAnalyzer::t_TrigPass

private

Definition at line 203 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

bool HcalIsoTrkAnalyzer::t_TrigPassSel

private

Definition at line 203 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

const std::string HcalIsoTrkAnalyzer::theTrackQuality_

private

Definition at line 140 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

const edm::InputTag HcalIsoTrkAnalyzer::theTriggerResultsLabel_

private

Definition at line 155 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer().

edm::EDGetTokenT<BXVector<GlobalAlgBlk> > HcalIsoTrkAnalyzer::tok_alg_

private

Definition at line 178 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

edm::ESGetToken<MagneticField, IdealMagneticFieldRecord> HcalIsoTrkAnalyzer::tok_bFieldH_

private

Definition at line 182 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

edm::EDGetTokenT<reco::BeamSpot> HcalIsoTrkAnalyzer::tok_bs_

private

Definition at line 172 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

edm::EDGetTokenT<CaloTowerCollection> HcalIsoTrkAnalyzer::tok_cala_

private

Definition at line 176 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

edm::ESGetToken<CaloTopology, CaloTopologyRecord> HcalIsoTrkAnalyzer::tok_caloTopology_

private

Definition at line 186 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

edm::ESGetToken<HcalDDDRecConstants, HcalRecNumberingRecord> HcalIsoTrkAnalyzer::tok_ddrec_

private

Definition at line 181 of file HcalIsoTrkAnalyzer.cc.

Referenced by beginRun(), and HcalIsoTrkAnalyzer().

edm::EDGetTokenT<EcalRecHitCollection> HcalIsoTrkAnalyzer::tok_EB_

private

Definition at line 173 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

edm::ESGetToken<EcalChannelStatus, EcalChannelStatusRcd> HcalIsoTrkAnalyzer::tok_ecalChStatus_

private

Definition at line 183 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

edm::EDGetTokenT<EcalRecHitCollection> HcalIsoTrkAnalyzer::tok_EE_

private

Definition at line 174 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

edm::EDGetTokenT<GenEventInfoProduct> HcalIsoTrkAnalyzer::tok_ew_

private

Definition at line 177 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

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

private

Definition at line 170 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

edm::ESGetToken<CaloGeometry, CaloGeometryRecord> HcalIsoTrkAnalyzer::tok_geom_

private

Definition at line 185 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

edm::EDGetTokenT<HBHERecHitCollection> HcalIsoTrkAnalyzer::tok_hbhe_

private

Definition at line 175 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

edm::ESGetToken<HcalTopology, HcalRecNumberingRecord> HcalIsoTrkAnalyzer::tok_htopo_

private

Definition at line 187 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

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

private

Definition at line 179 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

edm::EDGetTokenT<reco::GenParticleCollection> HcalIsoTrkAnalyzer::tok_parts_

private

Definition at line 169 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

edm::EDGetTokenT<reco::VertexCollection> HcalIsoTrkAnalyzer::tok_recVtx_

private

Definition at line 171 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

edm::ESGetToken<HcalRespCorrs, HcalRespCorrsRcd> HcalIsoTrkAnalyzer::tok_resp_

private

Definition at line 188 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

edm::ESGetToken<EcalSeverityLevelAlgo, EcalSeverityLevelAlgoRcd> HcalIsoTrkAnalyzer::tok_sevlv_

private

Definition at line 184 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

edm::EDGetTokenT<trigger::TriggerEvent> HcalIsoTrkAnalyzer::tok_trigEvt_

private

Definition at line 167 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

edm::EDGetTokenT<edm::TriggerResults> HcalIsoTrkAnalyzer::tok_trigRes_

private

Definition at line 168 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

TTree* HcalIsoTrkAnalyzer::tree

private

Definition at line 190 of file HcalIsoTrkAnalyzer.cc.

Referenced by eventstfile.Events::__iter__(), svgfig.SVG::__str__(), beginJob(), ZJetsTreeAnalyzer.ZJetsTreeAnalyzer::beginLoop(), MetTreeProducer.MetTreeProducer::declareVariables(), core.AutoFillTreeProducer.AutoFillTreeProducer::declareVariables(), fillTree(), core.AutoFillTreeProducer.AutoFillTreeProducer::fillTree(), ZJetsTreeAnalyzer.ZJetsTreeAnalyzer::process(), MetTreeProducer.MetTreeProducer::process(), python.cmstools.EventTree::SetAlias(), and eventstfile.Events::to().

TTree * HcalIsoTrkAnalyzer::tree2

private

Definition at line 190 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and beginJob().

const edm::InputTag HcalIsoTrkAnalyzer::triggerEvent_

private

Definition at line 155 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

const std::vector<std::string> HcalIsoTrkAnalyzer::trigNames_

private

Definition at line 138 of file HcalIsoTrkAnalyzer.cc.

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

const bool HcalIsoTrkAnalyzer::unCorrect_

private

Definition at line 151 of file HcalIsoTrkAnalyzer.cc.

Referenced by HcalIsoTrkAnalyzer(), and storeEnergy().

const bool HcalIsoTrkAnalyzer::useL1Trigger_

private

Definition at line 150 of file HcalIsoTrkAnalyzer.cc.

Referenced by analyze(), and HcalIsoTrkAnalyzer().

const int HcalIsoTrkAnalyzer::useRaw_

private

Definition at line 149 of file HcalIsoTrkAnalyzer.cc.

Referenced by fillTree(), and HcalIsoTrkAnalyzer().

Public Member Functions

Static Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation