test
CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
List of all members | Public Member Functions | Static Public Member Functions | Private Member Functions | Private Attributes
IsoTrackCalib Class Reference

#include <IsoTrackCalib.h>

Inheritance diagram for IsoTrackCalib:
edm::EDAnalyzer edm::EDAnalyzer edm::EDConsumerBase edm::EDConsumerBase

Public Member Functions

double dR (double eta1, double eta2, double phi1, double phi2)
 
double dR (double eta1, double eta2, double phi1, double phi2)
 
 IsoTrackCalib (const edm::ParameterSet &)
 
 IsoTrackCalib (const edm::ParameterSet &)
 
 ~IsoTrackCalib ()
 
 ~IsoTrackCalib ()
 
- Public Member Functions inherited from edm::EDAnalyzer
void callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
 
 EDAnalyzer ()
 
ModuleDescription const & moduleDescription () const
 
std::string workerType () const
 
virtual ~EDAnalyzer ()
 
- Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo > consumesInfo () const
 
 EDConsumerBase ()
 
 EDConsumerBase (EDConsumerBase const &)=delete
 
 EDConsumerBase (EDConsumerBase &&)=default
 
ProductResolverIndexAndSkipBit indexFrom (EDGetToken, BranchType, TypeID const &) const
 
void itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
 
void itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
 
std::vector
< ProductResolverIndexAndSkipBit >
const & 		itemsToGetFromEvent () const
 
void labelsForToken (EDGetToken iToken, Labels &oLabels) const
 
void modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
 
EDConsumerBase const & operator= (EDConsumerBase const &)=delete
 
EDConsumerBase & operator= (EDConsumerBase &&)=default
 
bool registeredToConsume (ProductResolverIndex, bool, BranchType) const
 
bool registeredToConsumeMany (TypeID const &, BranchType) const
 
void updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
 
virtual ~EDConsumerBase () noexcept(false)
 

Static Public Member Functions

static void fillDescriptions (edm::ConfigurationDescriptions &descriptions)
 
static void fillDescriptions (edm::ConfigurationDescriptions &descriptions)
 
- Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string & baseType ()
 
static void fillDescriptions (ConfigurationDescriptions &descriptions)
 
static void prevalidate (ConfigurationDescriptions &)
 

Private Member Functions

virtual void analyze (const edm::Event &, const edm::EventSetup &)
 
virtual void analyze (const edm::Event &, const edm::EventSetup &)
 
virtual void beginJob ()
 
virtual void beginJob ()
 
virtual void beginLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &)
 
virtual void beginLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &)
 
virtual void beginRun (edm::Run const &, edm::EventSetup const &)
 
virtual void beginRun (edm::Run const &, edm::EventSetup const &)
 
void clearTreeVectors ()
 
double deltaR (double eta1, double eta2, double phi1, double phi2)
 
double dEta (math::XYZTLorentzVector &, math::XYZTLorentzVector &)
 
double dEta (math::XYZTLorentzVector &, math::XYZTLorentzVector &)
 
double dinvPt (math::XYZTLorentzVector &, math::XYZTLorentzVector &)
 
double dP (math::XYZTLorentzVector &, math::XYZTLorentzVector &)
 
double dPhi (math::XYZTLorentzVector &, math::XYZTLorentzVector &)
 
double dPhi (math::XYZTLorentzVector &, math::XYZTLorentzVector &)
 
double dPt (math::XYZTLorentzVector &, math::XYZTLorentzVector &)
 
double dR (math::XYZTLorentzVector &, math::XYZTLorentzVector &)
 
double dR (math::XYZTLorentzVector &, math::XYZTLorentzVector &)
 
virtual void endJob ()
 
virtual void endJob ()
 
virtual void endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &)
 
virtual void endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &)
 
virtual void endRun (edm::Run const &, edm::EventSetup const &)
 
virtual void endRun (edm::Run const &, edm::EventSetup const &)
 

Private Attributes

double a_charIsoR
 
double a_coneR
 
double a_mipR
 
double a_neutIsoR
 
double a_neutR1
 
double a_neutR2
 
bool changed
 
double cutCharge
 
double cutMip
 
double cutNeutral
 
double dr_L1
 
std::vector< double > drCuts
 
int Event
 
bool firstEvent
 
edm::Service< TFileService > fs
 
TH1I * g_Accepts
 
TH1I * g_Pre
 
TH1I * g_PreHLT
 
TH1I * g_PreL1
 
TH1I * h_HLTAccept
 
std::vector< TH1I * > h_HLTAccepts
 
TH1I * h_iEta
 
TH1F * h_jetp
 
TH1F * h_jetpt [4]
 
TH1I * h_nHLT
 
TH1F * h_Rechit_E
 
TProfile * h_RecHit_iEta
 
TProfile * h_RecHit_num
 
TH1I * h_tkEta0 [5]
 
TH1I * h_tketa0 [6]
 
TH1I * h_tkEta1 [5]
 
TH1I * h_tketa1 [6]
 
TH1I * h_tkEta2 [5]
 
TH1I * h_tketa2 [6]
 
TH1I * h_tkEta3 [5]
 
TH1I * h_tketa3 [6]
 
TH1I * h_tkEta4 [5]
 
TH1I * h_tketa4 [6]
 
TH1I * h_tkEta5 [5]
 
TH1I * h_tketa5 [6]
 
HLTConfigProvider hltConfig_
 
std::vector< std::string > HLTNames
 
HLTPrescaleProvider hltPrescaleProvider_
 
bool initL1
 
std::map< std::pair< unsigned
int, std::string >, int > 		l1AlgoMap
 
std::vector< std::string > l1Names
 
const L1GtTriggerMenu * m_l1GtMenu
 
L1GtUtils m_l1GtUtils
 
int maxRunNo
 
int minRunNo
 
int nRun
 
std::vector< double > pbin
 
int Run
 
spr::trackSelectionParameters selectionParameters
 
std::vector< double > * t_charge_iso
 
std::vector< unsigned int > * t_DetIds
 
std::vector< double > * t_disthotcell
 
std::vector< double > * t_ehcal
 
double t_eHcal
 
std::vector< double > * t_emip
 
double t_eMipDR
 
int t_Event
 
double t_EventWeight
 
std::vector< double > * t_eventweight
 
std::vector< double > * t_HitEnergies
 
double t_hmaxNearP
 
int t_ieta
 
std::vector< int > * t_ieta
 
std::vector< double > * t_ietahotcell
 
std::vector< bool > * t_l1bits
 
double t_l1eta
 
std::vector< double > * t_l1eta
 
double t_l1phi
 
std::vector< double > * t_l1phi
 
double t_l1pt
 
std::vector< double > * t_l1pt
 
double t_l3eta
 
std::vector< double > * t_l3eta
 
double t_l3phi
 
std::vector< double > * t_l3phi
 
double t_l3pt
 
std::vector< double > * t_l3pt
 
std::vector< double > * t_leadingeta
 
std::vector< double > * t_leadingphi
 
std::vector< double > * t_leadingpt
 
double t_mindR1
 
double t_mindR2
 
std::vector< double > * t_neu_iso
 
double t_p
 
bool t_qltyFlag
 
bool t_qltyMissFlag
 
bool t_qltyPVFlag
 
int t_Run
 
bool t_selectTk
 
std::vector< double > * t_trackEta
 
std::vector< double > * t_trackP
 
std::vector< double > * t_trackPhi
 
std::vector< double > * t_trackPt
 
std::vector< double > * t_trackPx
 
std::vector< double > * t_trackPy
 
std::vector< double > * t_trackPz
 
std::vector< double > * t_trkL3mindr
 
std::string theTrackQuality
 
edm::InputTag theTriggerResultsLabel
 
edm::EDGetTokenT< reco::BeamSpot > tok_bs_
 
edm::EDGetTokenT
< EcalRecHitCollection > 		tok_EB_
 
edm::EDGetTokenT
< EcalRecHitCollection > 		tok_EE_
 
edm::EDGetTokenT
< GenEventInfoProduct > 		tok_ew_
 
edm::EDGetTokenT
< reco::TrackCollection > 		tok_genTrack_
 
edm::EDGetTokenT
< HBHERecHitCollection > 		tok_hbhe_
 
edm::EDGetTokenT
< reco::GenJetCollection > 		tok_jets_
 
edm::EDGetTokenT
< l1extra::L1JetParticleCollection > 		tok_L1extCenJet_
 
edm::EDGetTokenT
< l1extra::L1JetParticleCollection > 		tok_L1extFwdJet_
 
edm::EDGetTokenT
< l1extra::L1JetParticleCollection > 		tok_L1extTauJet_
 
edm::EDGetTokenT< LumiDetails > tok_lumi
 
edm::EDGetTokenT
< reco::PFJetCollection > 		tok_pf_
 
edm::EDGetTokenT
< reco::PFJetCollection > 		tok_pfjets_
 
edm::EDGetTokenT
< reco::VertexCollection > 		tok_recVtx_
 
edm::EDGetTokenT
< trigger::TriggerEvent > 		tok_trigEvt
 
edm::EDGetTokenT
< edm::TriggerResults > 		tok_trigRes
 
TTree * tree
 
std::vector< std::string > trgnames
 
edm::InputTag triggerEvent_
 
std::map< std::pair< unsigned
int, std::string >, unsigned
int > 		TrigList
 
std::vector< std::string > trigNames
 
std::map< std::pair< unsigned
int, std::string >, const
std::pair< int, int > > 		TrigPreList
 
int verbosity
 

Additional Inherited Members

- Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer ModuleType
 
- Public Types inherited from edm::EDConsumerBase
typedef ProductLabels Labels
 
- Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > consumes (edm::InputTag const &tag)
 
EDGetToken consumes (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken consumes (TypeToGet const &id, edm::InputTag const &tag)
 
ConsumesCollector consumesCollector ()
 Use a ConsumesCollector to gather consumes information from helper functions. More...
 
template<typename ProductType , BranchType B = InEvent>
void consumesMany ()
 
void consumesMany (const TypeToGet &id)
 
template<BranchType B>
void consumesMany (const TypeToGet &id)
 
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > mayConsume (edm::InputTag const &tag)
 
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 

Detailed Description

Definition at line 79 of file IsoTrackCalib.cc.

Constructor & Destructor Documentation

IsoTrackCalib::IsoTrackCalib ( const edm::ParameterSet &  iConfig)
explicit

Definition at line 151 of file IsoTrackCalib.cc.

References a_charIsoR, a_coneR, a_mipR, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), HLT_25ns10e33_v2_cff::InputTag, relval_2017::k, l1Names, spr::trackSelectionParameters::maxChi2, spr::trackSelectionParameters::maxDpOverP, spr::trackSelectionParameters::maxDxyPV, spr::trackSelectionParameters::maxDzPV, spr::trackSelectionParameters::maxInMiss, spr::trackSelectionParameters::maxOutMiss, spr::trackSelectionParameters::minLayerCrossed, spr::trackSelectionParameters::minOuterHit, spr::trackSelectionParameters::minPt, spr::trackSelectionParameters::minQuality, reco::TrackBase::qualityByName(), selectionParameters, AlCaHLTBitMon_QueryRunRegistry::string, theTrackQuality, theTriggerResultsLabel, tok_bs_, tok_EB_, tok_EE_, tok_ew_, tok_genTrack_, tok_hbhe_, tok_jets_, tok_L1extCenJet_, tok_L1extFwdJet_, tok_L1extTauJet_, tok_pfjets_, tok_recVtx_, tok_trigEvt, tok_trigRes, triggerEvent_, trigNames, and verbosity.

151  :
152  m_l1GtUtils(iConfig, consumesCollector(), useL1GtTriggerMenuLite, *this){
153  //now do whatever initialization is needed
154  verbosity = iConfig.getUntrackedParameter<int>("Verbosity",0);
155  l1Names = iConfig.getUntrackedParameter<std::vector<std::string> >("L1Seed");
156  theTrackQuality = iConfig.getUntrackedParameter<std::string>("TrackQuality","highPurity");
157  reco::TrackBase::TrackQuality trackQuality_=reco::TrackBase::qualityByName(theTrackQuality);
158  selectionParameters.minPt = iConfig.getUntrackedParameter<double>("MinTrackPt", 10.0);
159  selectionParameters.minQuality = trackQuality_;
160  selectionParameters.maxDxyPV = iConfig.getUntrackedParameter<double>("MaxDxyPV", 0.2);
161  selectionParameters.maxDzPV = iConfig.getUntrackedParameter<double>("MaxDzPV", 5.0);
162  selectionParameters.maxChi2 = iConfig.getUntrackedParameter<double>("MaxChi2", 5.0);
163  selectionParameters.maxDpOverP = iConfig.getUntrackedParameter<double>("MaxDpOverP", 0.1);
164  selectionParameters.minOuterHit = iConfig.getUntrackedParameter<int>("MinOuterHit", 4);
165  selectionParameters.minLayerCrossed = iConfig.getUntrackedParameter<int>("MinLayerCrossed", 8);
166  selectionParameters.maxInMiss = iConfig.getUntrackedParameter<int>("MaxInMiss", 0);
167  selectionParameters.maxOutMiss = iConfig.getUntrackedParameter<int>("MaxOutMiss", 0);
168  a_coneR = iConfig.getUntrackedParameter<double>("ConeRadius",34.98);
169  a_charIsoR = a_coneR + 28.9;
170  a_mipR = iConfig.getUntrackedParameter<double>("ConeRadiusMIP",14.0);
171  bool isItAOD = iConfig.getUntrackedParameter<bool>("IsItAOD", false);
172  triggerEvent_ = edm::InputTag("hltTriggerSummaryAOD","","HLT");
173  theTriggerResultsLabel = edm::InputTag("TriggerResults","","HLT");
174  edm::InputTag L1extraTauJetSource_ = iConfig.getParameter<edm::InputTag> ("L1extraTauJetSource" );
175  edm::InputTag L1extraCenJetSource_ = iConfig.getParameter<edm::InputTag> ("L1extraCenJetSource" );
176  edm::InputTag L1extraFwdJetSource_ = iConfig.getParameter<edm::InputTag> ("L1extraFwdJetSource" );
177 
178  // define tokens for access
179  tok_trigEvt = consumes<trigger::TriggerEvent>(triggerEvent_);
180  tok_trigRes = consumes<edm::TriggerResults>(theTriggerResultsLabel);
181  tok_genTrack_ = consumes<reco::TrackCollection>(edm::InputTag("generalTracks"));
182  tok_recVtx_ = consumes<reco::VertexCollection>(edm::InputTag("offlinePrimaryVertices"));
183  tok_bs_ = consumes<reco::BeamSpot>(edm::InputTag("offlineBeamSpot"));
184  tok_ew_ = consumes<GenEventInfoProduct>(edm::InputTag("generatorSmeared"));
185  tok_L1extTauJet_ = consumes<l1extra::L1JetParticleCollection>(L1extraTauJetSource_);
186  tok_L1extCenJet_ = consumes<l1extra::L1JetParticleCollection>(L1extraCenJetSource_);
187  tok_L1extFwdJet_ = consumes<l1extra::L1JetParticleCollection>(L1extraFwdJetSource_);
188  if (isItAOD) {
189  tok_EB_ = consumes<EcalRecHitCollection>(edm::InputTag("reducedEcalRecHitsEB"));
190  tok_EE_ = consumes<EcalRecHitCollection>(edm::InputTag("reducedEcalRecHitsEE"));
191  tok_hbhe_ = consumes<HBHERecHitCollection>(edm::InputTag("reducedHcalRecHits", "hbhereco"));
192  } else {
193  tok_EB_ = consumes<EcalRecHitCollection>(edm::InputTag("ecalRecHit","EcalRecHitsEB"));
194  tok_EE_ = consumes<EcalRecHitCollection>(edm::InputTag("ecalRecHit","EcalRecHitsEE"));
195  tok_hbhe_ = consumes<HBHERecHitCollection>(edm::InputTag("hbhereco"));
196  }
197  tok_jets_ = consumes<reco::GenJetCollection>(iConfig.getParameter<edm::InputTag>("JetSource"));
198  tok_pfjets_ = consumes<reco::PFJetCollection>(edm::InputTag("ak5PFJets"));
199  if (verbosity>=0) {
200  edm::LogInfo("IsoTrack")
201  <<"Parameters read from config file \n"
202  <<"\t minPt " << selectionParameters.minPt
203  <<"\t theTrackQuality " << theTrackQuality
204  <<"\t minQuality " << selectionParameters.minQuality
205  <<"\t maxDxyPV " << selectionParameters.maxDxyPV
206  <<"\t maxDzPV " << selectionParameters.maxDzPV
207  <<"\t maxChi2 " << selectionParameters.maxChi2
208  <<"\t maxDpOverP " << selectionParameters.maxDpOverP
209  <<"\t minOuterHit " << selectionParameters.minOuterHit
210  <<"\t minLayerCrossed " << selectionParameters.minLayerCrossed
211  <<"\t maxInMiss " << selectionParameters.maxInMiss
212  <<"\t maxOutMiss " << selectionParameters.maxOutMiss
213  <<"\t a_coneR " << a_coneR
214  <<"\t a_charIsoR " << a_charIsoR
215  <<"\t a_mipR " << a_mipR
216  <<"\t isItAOD " << isItAOD;
217  edm::LogInfo("IsoTrack") << trigNames.size() << " triggers to be studied";
218  for (unsigned int k=0; k<trigNames.size(); ++k)
219  edm::LogInfo("IsoTrack") << "[" << k << "]: " << trigNames[k];
220  }
221 
222 }
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
IsoTrackCalib::tok_L1extCenJet_
edm::EDGetTokenT< l1extra::L1JetParticleCollection > tok_L1extCenJet_
Definition: IsoTrackCalib.cc:127
IsoTrackCalib::tok_jets_
edm::EDGetTokenT< reco::GenJetCollection > tok_jets_
Definition: IsoTrackCalib.cc:124
IsoTrackCalib::triggerEvent_
edm::InputTag triggerEvent_
Definition: IsoTrackCalib.cc:113
reco::TrackBase::TrackQuality
TrackQuality
track quality
Definition: TrackBase.h:151
IsoTrackCalib::l1Names
std::vector< std::string > l1Names
Definition: IsoTrackCalib.cc:106
relval_2017.k
list k
Definition: relval_2017.py:21
IsoTrackCalib::tok_pfjets_
edm::EDGetTokenT< reco::PFJetCollection > tok_pfjets_
Definition: IsoTrackCalib.cc:125
IsoTrackCalib::tok_L1extTauJet_
edm::EDGetTokenT< l1extra::L1JetParticleCollection > tok_L1extTauJet_
Definition: IsoTrackCalib.cc:126
IsoTrackCalib::tok_hbhe_
edm::EDGetTokenT< HBHERecHitCollection > tok_hbhe_
Definition: IsoTrackCalib.cc:122
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
IsoTrackCalib::tok_genTrack_
edm::EDGetTokenT< reco::TrackCollection > tok_genTrack_
Definition: IsoTrackCalib.cc:117
HLT_25ns10e33_v2_cff.InputTag
tuple InputTag
Definition: HLT_25ns10e33_v2_cff.py:58467
IsoTrackCalib::theTrackQuality
std::string theTrackQuality
Definition: IsoTrackCalib.cc:109
IsoTrackCalib::tok_trigRes
edm::EDGetTokenT< edm::TriggerResults > tok_trigRes
Definition: IsoTrackCalib.cc:115
useL1GtTriggerMenuLite
static const bool useL1GtTriggerMenuLite(true)
edm::EDConsumerBase::consumesCollector
ConsumesCollector consumesCollector()
Use a ConsumesCollector to gather consumes information from helper functions.
Definition: EDConsumerBase.cc:72
IsoTrackCalib::tok_L1extFwdJet_
edm::EDGetTokenT< l1extra::L1JetParticleCollection > tok_L1extFwdJet_
Definition: IsoTrackCalib.cc:128
IsoTrackCalib::tok_EE_
edm::EDGetTokenT< EcalRecHitCollection > tok_EE_
Definition: IsoTrackCalib.cc:121
reco::TrackBase::qualityByName
static TrackQuality qualityByName(const std::string &name)
Definition: TrackBase.cc:125
IsoTrackCalib::m_l1GtUtils
L1GtUtils m_l1GtUtils
Definition: IsoTrackCalib.cc:104
IsoTrackCalib::tok_ew_
edm::EDGetTokenT< GenEventInfoProduct > tok_ew_
Definition: IsoTrackCalib.cc:123
IsoTrackCalib::selectionParameters
spr::trackSelectionParameters selectionParameters
Definition: IsoTrackCalib.cc:108
spr::trackSelectionParameters::minQuality
reco::TrackBase::TrackQuality minQuality
Definition: TrackSelection.h:22
IsoTrackCalib::tok_EB_
edm::EDGetTokenT< EcalRecHitCollection > tok_EB_
Definition: IsoTrackCalib.cc:120
IsoTrackCalib::theTriggerResultsLabel
edm::InputTag theTriggerResultsLabel
Definition: IsoTrackCalib.cc:113
IsoTrackCalib::tok_bs_
edm::EDGetTokenT< reco::BeamSpot > tok_bs_
Definition: IsoTrackCalib.cc:119
IsoTrackCalib::tok_trigEvt
edm::EDGetTokenT< trigger::TriggerEvent > tok_trigEvt
Definition: IsoTrackCalib.cc:114
IsoTrackCalib::trigNames
std::vector< std::string > trigNames
Definition: IsoTrackCalib.cc:106
IsoTrackCalib::tok_recVtx_
edm::EDGetTokenT< reco::VertexCollection > tok_recVtx_
Definition: IsoTrackCalib.cc:118
IsoTrackCalib::~IsoTrackCalib ( )

Definition at line 224 of file IsoTrackCalib.cc.

224  {
225  // do anything here that needs to be done at desctruction time
226  // (e.g. close files, deallocate resources etc.)
227 
228 }
IsoTrackCalib::IsoTrackCalib ( const edm::ParameterSet &  )
explicit
IsoTrackCalib::~IsoTrackCalib ( )

Member Function Documentation

virtual void IsoTrackCalib::analyze ( const edm::Event &  ,
const edm::EventSetup &   
)
privatevirtual

Implements edm::EDAnalyzer.

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

Implements edm::EDAnalyzer.

Definition at line 230 of file IsoTrackCalib.cc.

References a_charIsoR, a_coneR, a_mipR, L1GtUtils::availableL1Configuration(), ecalTB2006H4_GenSimDigiReco_cfg::bField, edm::EventBase::bunchCrossing(), spr::chargeIsolationCone(), gather_cfg::cout, L1GtUtils::decision(), deltaR(), spr::eCone_ecal(), spr::eCone_hcal(), reco::LeafCandidate::eta(), reco::TrackBase::eta(), edm::EventID::event(), edm::EventSetup::get(), edm::Event::getByToken(), L1GtUtils::getL1GtRunCache(), spr::goodTrack(), L1GtTriggerMenu::gtAlgorithmMap(), L1GtTriggerMenu::gtTriggerMenuName(), h_iEta, h_jetpt, h_Rechit_E, h_RecHit_iEta, h_RecHit_num, h_tketa0, h_tketa1, h_tketa2, h_tketa3, h_tketa4, h_tketa5, i, edm::EventBase::id(), iErrorCode, HcalDetId::ieta(), edm::HandleBase::isValid(), relval_2017::k, l1Names, edm::EventBase::luminosityBlock(), m_l1GtMenu, m_l1GtUtils, spr::trackSelectionParameters::maxDxyPV, spr::trackSelectionParameters::maxDzPV, spr::trackSelectionParameters::maxInMiss, spr::trackSelectionParameters::maxOutMiss, nTracks(), reco::TrackBase::p(), pbin, pfJets_cff::pfJets, reco::TrackBase::phi(), edm::ESHandle< class >::product(), spr::propagateCALO(), reco::LeafCandidate::pt(), reco::TrackBase::pt(), reco::TrackBase::px(), reco::TrackBase::py(), reco::TrackBase::pz(), edm::EventID::run(), selectionParameters, AlCaHLTBitMon_QueryRunRegistry::string, t_DetIds, t_eHcal, t_eMipDR, t_Event, t_EventWeight, t_HitEnergies, t_hmaxNearP, t_ieta, t_l1bits, t_l1eta, t_l1phi, t_l1pt, t_l3eta, t_l3phi, t_l3pt, t_mindR1, t_mindR2, t_p, t_qltyFlag, t_qltyMissFlag, t_qltyPVFlag, t_Run, t_selectTk, theTrackQuality, tok_bs_, tok_EB_, tok_EE_, tok_ew_, tok_genTrack_, tok_hbhe_, tok_jets_, tok_L1extCenJet_, tok_L1extFwdJet_, tok_L1extTauJet_, tok_pfjets_, tok_recVtx_, tree, useL1EventSetup(), useL1GtTriggerMenuLite(), verbosity, x, y, and z.

231  {
232 
233  t_Run = iEvent.id().run();
234  t_Event = iEvent.id().event();
235  if (verbosity%10 > 0)
236  edm::LogInfo("IsoTrack")
237  << "Run " << t_Run << " Event " << t_Event << " Luminosity "
238  << iEvent.luminosityBlock() << " Bunch " << iEvent.bunchCrossing()
239  << " starts ==========";
240 
241  //Get magnetic field and ECAL channel status
242  edm::ESHandle<MagneticField> bFieldH;
243  iSetup.get<IdealMagneticFieldRecord>().get(bFieldH);
244  const MagneticField *bField = bFieldH.product();
245 
246  edm::ESHandle<EcalSeverityLevelAlgo> sevlv;
247  iSetup.get<EcalSeverityLevelAlgoRcd>().get(sevlv);
248 
249  // get handles to calogeometry and calotopology
250  edm::ESHandle<CaloGeometry> pG;
251  iSetup.get<CaloGeometryRecord>().get(pG);
252  const CaloGeometry* geo = pG.product();
253 
254  //Get track collection
255  edm::Handle<reco::TrackCollection> trkCollection;
256  iEvent.getByToken(tok_genTrack_, trkCollection);
257  reco::TrackCollection::const_iterator trkItr;
258 
259  //event weight for FLAT sample
260  t_EventWeight = 1.0;
261  edm::Handle<GenEventInfoProduct> genEventInfo;
262  iEvent.getByToken(tok_ew_, genEventInfo);
263  if (genEventInfo.isValid()) t_EventWeight = genEventInfo->weight();
264 
265  // genJet information
266  edm::Handle<reco::GenJetCollection> genJets;
267  iEvent.getByToken(tok_jets_, genJets);
268  if (genJets.isValid()) {
269  for (unsigned iGenJet = 0; iGenJet < genJets->size(); ++iGenJet) {
270  const reco::GenJet& genJet = (*genJets) [iGenJet];
271  double genJetPt = genJet.pt();
272  double genJetEta = genJet.eta();
273  h_jetpt[0]->Fill(genJetPt);
274  h_jetpt[1]->Fill(genJetPt,t_EventWeight);
275  if (genJetEta>-2.5 && genJetEta<2.5) {
276  h_jetpt[2]->Fill(genJetPt);
277  h_jetpt[3]->Fill(genJetPt,t_EventWeight);
278  }
279  break;
280  }
281  }
282 
283  //pf jets
284  edm::Handle<reco::PFJetCollection> pfJets;
285  iEvent.getByToken(tok_pfjets_, pfJets);
286  reco::PFJetCollection::const_iterator pfItr;
287 
288  //Define the best vertex and the beamspot
289  edm::Handle<reco::VertexCollection> recVtxs;
290  iEvent.getByToken(tok_recVtx_, recVtxs);
291  edm::Handle<reco::BeamSpot> beamSpotH;
292  iEvent.getByToken(tok_bs_, beamSpotH);
293  math::XYZPoint leadPV(0,0,0);
294  if (recVtxs->size()>0 && !((*recVtxs)[0].isFake())) {
295  leadPV = math::XYZPoint( (*recVtxs)[0].x(),(*recVtxs)[0].y(), (*recVtxs)[0].z() );
296  } else if (beamSpotH.isValid()) {
297  leadPV = beamSpotH->position();
298  }
299  if (verbosity>10) {
300  if ((verbosity%100)/10>2)
301  edm::LogInfo("IsoTrack") << "Primary Vertex " << leadPV;
302  if (beamSpotH.isValid()) edm::LogInfo("IsoTrack") << " Beam Spot "
303  << beamSpotH->position();
304  }
305 
306  // RecHits
307  edm::Handle<EcalRecHitCollection> barrelRecHitsHandle;
308  edm::Handle<EcalRecHitCollection> endcapRecHitsHandle;
309  iEvent.getByToken(tok_EB_, barrelRecHitsHandle);
310  iEvent.getByToken(tok_EE_, endcapRecHitsHandle);
311  edm::Handle<HBHERecHitCollection> hbhe;
312  iEvent.getByToken(tok_hbhe_, hbhe);
313  HBHERecHitCollection::const_iterator rhitItr;
314 
315  for (rhitItr=hbhe->begin();rhitItr!=hbhe->end();rhitItr++) {
316  double rec_energy = rhitItr->energy();
317  int rec_ieta = rhitItr->id().ieta();
318  int rec_depth = rhitItr->id().depth();
319  int rec_zside = rhitItr->id().zside();
320  double num1_1 = rec_zside*(rec_ieta+0.2*(rec_depth-1));
321  if (verbosity%10>0)
322  edm::LogInfo("IsoTrack") << "detid/rechit/ieta/zside/depth/num " << " = "
323  << rhitItr->id() << "/" << rec_energy << "/"
324  << rec_ieta << "/" << rec_zside << "/"
325  << rec_depth << "/" << num1_1;
326  h_iEta->Fill(rec_ieta);
327  h_Rechit_E->Fill(rec_energy);
328  h_RecHit_iEta->Fill(rec_ieta,rec_energy);
329  h_RecHit_num->Fill(num1_1,rec_energy);
330  }
331 
332  //Propagate tracks to calorimeter surface)
333  std::vector<spr::propagatedTrackDirection> trkCaloDirections;
334  spr::propagateCALO(trkCollection, geo, bField, theTrackQuality,
335  trkCaloDirections, ((verbosity/100)%10>2));
336  std::vector<spr::propagatedTrackDirection>::const_iterator trkDetItr;
337  for (trkDetItr = trkCaloDirections.begin();
338  trkDetItr != trkCaloDirections.end(); trkDetItr++) {
339  if (trkDetItr->okHCAL) {
340  HcalDetId detId = (HcalDetId)(trkDetItr->detIdHCAL);
341  int tk_ieta = detId.ieta();
342  const reco::Track* pTrack = &(*(trkDetItr->trkItr));
343  double tk_p = pTrack->p();
344  h_tketa0[0]->Fill(tk_ieta);
345  for (unsigned int k=1; k<pbin.size(); ++k) {
346  if (tk_p >= pbin[k-1] && tk_p < pbin[k]) {
347  h_tketa0[k]->Fill(tk_ieta);
348  break;
349  }
350  }
351  }
352  }
354  t_l1bits->clear();
355  for (unsigned int i=0; i<l1Names.size(); ++i) t_l1bits->push_back(false);
356  bool useL1EventSetup = true;
357  bool useL1GtTriggerMenuLite = true;
358 
359  m_l1GtUtils.getL1GtRunCache(iEvent, iSetup, useL1EventSetup, useL1GtTriggerMenuLite);
360  int iErrorCode = -1;
361  int l1ConfCode = -1;
362  const bool l1Conf = m_l1GtUtils.availableL1Configuration(iErrorCode, l1ConfCode);
363  if( !l1Conf) {
364  edm::LogInfo("IsoTrack") << "\nL1 configuration code:" << l1ConfCode
365  << "\nNo valid L1 trigger configuration available."
366  << "\nSee text above for error code interpretation"
367  << "\nNo return here, in order to test each method"
368  << ", protected against configuration error.";
369  }
370 
371  const AlgorithmMap& algorithmMap = m_l1GtMenu->gtAlgorithmMap();
372  const std::string& menuName = m_l1GtMenu->gtTriggerMenuName();
373  std::cout << "menuName " << menuName << std::endl;
374 
375  std::vector<int> algbits;
376  for (CItAlgo itAlgo = algorithmMap.begin(); itAlgo != algorithmMap.end();
377  itAlgo++) {
378  std::string algName = itAlgo->first;
379  int algBitNumber = ( itAlgo->second ).algoBitNumber();
380  bool decision = m_l1GtUtils.decision(iEvent, itAlgo->first, iErrorCode);
381 
382  bool l1ok(false);
383  if (verbosity%10>0)
384  edm::LogInfo("IsoTrack") << algName <<" "<< algBitNumber << " "
385  << decision;
386  for (unsigned int i=0; i<l1Names.size(); ++i) {
387  if (algName.find(l1Names[i])!=std::string::npos){
388  if (verbosity%10>0)
389  edm::LogInfo("IsoTrack") << "match found" << " " << algName << " "
390  << decision;
391  t_l1bits->at(i) = (decision>0);
392  if (decision > 0) l1ok = true;
393  }
394  }
395  if (verbosity%10>0) edm::LogInfo("IsoTrack") << "l1 ok =" << l1ok;
396 
397  if(l1ok){
398  edm::Handle<l1extra::L1JetParticleCollection> l1TauHandle;
399  iEvent.getByToken(tok_L1extTauJet_,l1TauHandle);
400  l1extra::L1JetParticleCollection::const_iterator itr;
401  double ptTriggered = -10;
402  double etaTriggered = -100;
403  double phiTriggered = -100;
404 
405  for(itr = l1TauHandle->begin(); itr != l1TauHandle->end(); ++itr) {
406  if(itr->pt()>ptTriggered){
407  ptTriggered = itr->pt();
408  etaTriggered = itr->eta();
409  phiTriggered = itr->phi() ;
410  }
411  if (verbosity%10>0)
412  edm::LogInfo("IsoTrack") << "tauJ pt " << itr->pt() << " eta/phi "
413  << itr->eta() << " " << itr->phi();
414  }
415  edm::Handle<l1extra::L1JetParticleCollection> l1CenJetHandle;
416  iEvent.getByToken(tok_L1extCenJet_,l1CenJetHandle);
417  for (itr = l1CenJetHandle->begin(); itr != l1CenJetHandle->end(); ++itr){
418  if (itr->pt()>ptTriggered) {
419  ptTriggered = itr->pt();
420  etaTriggered = itr->eta();
421  phiTriggered = itr->phi() ;
422  }
423  if (verbosity%10>0)
424  edm::LogInfo("IsoTrack") << "cenJ pt " << itr->pt()
425  << " eta/phi " << itr->eta() << " "
426  << itr->phi();
427  }
428  edm::Handle<l1extra::L1JetParticleCollection> l1FwdJetHandle;
429  iEvent.getByToken(tok_L1extFwdJet_,l1FwdJetHandle);
430  for (itr = l1FwdJetHandle->begin(); itr != l1FwdJetHandle->end(); ++itr) {
431  if (itr->pt()>ptTriggered) {
432  ptTriggered = itr->pt();
433  etaTriggered = itr->eta();
434  phiTriggered = itr->phi() ;
435  }
436  if (verbosity%10>0)
437  edm::LogInfo("IsoTrack") << "forJ pt " << itr->pt() << " eta/phi "
438  << itr->eta() << " " << itr->phi();
439  }
440  if (verbosity%10>0)
441  edm::LogInfo("IsoTrack") << "jets pt/eta/phi = " << ptTriggered << "/"
442  << etaTriggered << "/" <<phiTriggered;
444  unsigned int nTracks(0),nselTracks(0);
445  for (trkDetItr = trkCaloDirections.begin(),nTracks=0;
446  trkDetItr != trkCaloDirections.end(); trkDetItr++,nTracks++) {
447  const reco::Track* pTrack = &(*(trkDetItr->trkItr));
448  math::XYZTLorentzVector v4(pTrack->px(), pTrack->py(),
449  pTrack->pz(), pTrack->p());
450 
451  t_mindR1 = deltaR(etaTriggered,v4.eta(),phiTriggered,v4.phi());
452  t_mindR2 = -999;
453  if (verbosity%10>0)
454  edm::LogInfo("IsoTrack") << "This track : " << nTracks << " (pt/eta/phi/p) :"
455  << pTrack->pt() << "/" << pTrack->eta() << "/"
456  << pTrack->phi() << "/" << pTrack->p();
457 
458  if (verbosity%10>0) edm::LogInfo("IsoTrack") << "dr values are = " << t_mindR1;
459 
460  t_l1pt = ptTriggered;
461  t_l1eta = etaTriggered;
462  t_l1phi = phiTriggered;
463  t_l3pt = -999;
464  t_l3eta = -999;
465  t_l3phi = -999;
466 
467  //Selection of good track
468  t_selectTk = spr::goodTrack(pTrack,leadPV,selectionParameters,((verbosity/100)%10>2));
469  spr::trackSelectionParameters oneCutParameters = selectionParameters;
470  oneCutParameters.maxDxyPV = 10;
471  oneCutParameters.maxDzPV = 100;
472  oneCutParameters.maxInMiss = 2;
473  oneCutParameters.maxOutMiss= 2;
474  bool qltyFlag = spr::goodTrack(pTrack,leadPV,oneCutParameters,((verbosity/100)%10>2));
475  oneCutParameters = selectionParameters;
476  oneCutParameters.maxDxyPV = 10;
477  oneCutParameters.maxDzPV = 100;
478  t_qltyMissFlag = spr::goodTrack(pTrack,leadPV,oneCutParameters,((verbosity/100)%10>2));
479  oneCutParameters = selectionParameters;
480  oneCutParameters.maxInMiss = 2;
481  oneCutParameters.maxOutMiss= 2;
482  t_qltyPVFlag = spr::goodTrack(pTrack,leadPV,oneCutParameters,((verbosity/100)%10>2));
483  t_ieta = 0;
484  if (trkDetItr->okHCAL) {
485  HcalDetId detId = (HcalDetId)(trkDetItr->detIdHCAL);
486  t_ieta = detId.ieta();
487  }
488  if (verbosity%10 > 0)
489  edm::LogInfo("IsoTrack") << "qltyFlag|okECAL|okHCAL : " << qltyFlag << "|"
490  << trkDetItr->okECAL << "/" << trkDetItr->okHCAL;
491  t_qltyFlag = (qltyFlag && trkDetItr->okECAL && trkDetItr->okHCAL);
492  t_p = pTrack->p();
493  h_tketa1[0]->Fill(t_ieta);
494  for (unsigned int k=1; k<pbin.size(); ++k) {
495  if (t_p >= pbin[k-1] && t_p < pbin[k]) {
496  h_tketa1[k]->Fill(t_ieta);
497  break;
498  }
499  }
500  if (t_qltyFlag) {
501  nselTracks++;
502  h_tketa2[0]->Fill(t_ieta);
503  for (unsigned int k=1; k<pbin.size(); ++k) {
504  if (t_p >= pbin[k-1] && t_p < pbin[k]) {
505  h_tketa2[k]->Fill(t_ieta);
506  break;
507  }
508  }
509  int nRH_eMipDR(0), nNearTRKs(0), nRecHits(-999);
510  t_eMipDR = spr::eCone_ecal(geo, barrelRecHitsHandle,
511  endcapRecHitsHandle,
512  trkDetItr->pointHCAL,
513  trkDetItr->pointECAL,
514  a_mipR, trkDetItr->directionECAL,
515  nRH_eMipDR);
516  t_DetIds->clear(); t_HitEnergies->clear();
517  std::vector<DetId> ids;
518  t_eHcal = spr::eCone_hcal(geo, hbhe, trkDetItr->pointHCAL,
519  trkDetItr->pointECAL, a_coneR,
520  trkDetItr->directionHCAL,nRecHits,
521  ids, *t_HitEnergies);
522  for (unsigned int k=0; k<ids.size(); ++k) {
523  t_DetIds->push_back(ids[k].rawId());
524  }
525  t_hmaxNearP = spr::chargeIsolationCone(nTracks,trkCaloDirections,
526  a_charIsoR, nNearTRKs,
527  ((verbosity/100)%10>2));
528  if (t_hmaxNearP < 2) {
529  h_tketa3[0]->Fill(t_ieta);
530  for (unsigned int k=1; k<pbin.size(); ++k) {
531  if (t_p >= pbin[k-1] && t_p < pbin[k]) {
532  h_tketa3[k]->Fill(t_ieta);
533  break;
534  }
535  }
536  if (t_eMipDR < 1) {
537  h_tketa4[0]->Fill(t_ieta);
538  for (unsigned int k=1; k<pbin.size(); ++k) {
539  if (t_p >= pbin[k-1] && t_p < pbin[k]) {
540  h_tketa4[k]->Fill(t_ieta);
541  break;
542  }
543  }
544  if (t_mindR1 > 1) {
545  h_tketa5[0]->Fill(t_ieta);
546  for (unsigned int k=1; k<pbin.size(); ++k) {
547  if (t_p >= pbin[k-1] && t_p < pbin[k]) {
548  h_tketa5[k]->Fill(t_ieta);
549  break;
550  }
551  }
552  }
553  }
554  }
555  if (verbosity%10 > 0) {
556  edm::LogInfo("IsoTrack") << "This track : " << nTracks << " (pt/eta/phi/p) :"
557  << pTrack->pt() << "/" << pTrack->eta() << "/"
558  << pTrack->phi() << "/" << t_p;
559  edm::LogInfo("IsoTrack") << "e_MIP " << t_eMipDR << " Chg Isolation "
560  << t_hmaxNearP << " eHcal" << t_eHcal << " ieta "
561  << t_ieta << " Quality " << t_qltyMissFlag
562  << ":" << t_qltyPVFlag << ":" << t_selectTk;
563  for (unsigned int lll=0;lll<t_DetIds->size();lll++) {
564  edm::LogInfo("IsoTrack") << "det id is = " << t_DetIds->at(lll) << " "
565  << " hit enery is = " << t_HitEnergies->at(lll);
566  }
567  }
568  if (t_p>20.0 && t_eMipDR<2.0 && t_hmaxNearP<2.0) {
569  tree->Fill();
570  }
571  }
572  }
573  }
574  }
575 }
edm::EventID::run
RunNumber_t run() const
Definition: EventID.h:39
reco::TrackBase::p
double p() const
momentum vector magnitude
Definition: TrackBase.h:610
L1GtTriggerMenu::gtTriggerMenuName
const std::string & gtTriggerMenuName() const
Definition: L1GtTriggerMenu.h:99
edm::EventID::event
EventNumber_t event() const
Definition: EventID.h:41
IsoTrackCalib::t_HitEnergies
std::vector< double > * t_HitEnergies
Definition: IsoTrackCalib.cc:139
i
int i
Definition: DBlmapReader.cc:9
nTracks
const unsigned int nTracks(const reco::Vertex &sv)
Definition: TemplatedVertexArbitrator.h:48
IsoTrackCalib::m_l1GtMenu
const L1GtTriggerMenu * m_l1GtMenu
Definition: IsoTrackCalib.cc:111
IsoTrackCalib::t_DetIds
std::vector< unsigned int > * t_DetIds
Definition: IsoTrackCalib.cc:138
spr::propagateCALO
std::vector< spr::propagatedTrackID > propagateCALO(edm::Handle< reco::TrackCollection > &trkCollection, const CaloGeometry *geo, const MagneticField *bField, std::string &theTrackQuality, bool debug=false)
Definition: CaloPropagateTrack.cc:110
IsoTrackCalib::h_RecHit_num
TProfile * h_RecHit_num
Definition: IsoTrackCalib.cc:140
IsoTrackCalib::tok_L1extCenJet_
edm::EDGetTokenT< l1extra::L1JetParticleCollection > tok_L1extCenJet_
Definition: IsoTrackCalib.cc:127
IsoTrackCalib::tok_jets_
edm::EDGetTokenT< reco::GenJetCollection > tok_jets_
Definition: IsoTrackCalib.cc:124
edm::Event::getByToken
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:457
IsoTrackCalib::l1Names
std::vector< std::string > l1Names
Definition: IsoTrackCalib.cc:106
AlgorithmMap
std::map< std::string, L1GtAlgorithm > AlgorithmMap
map containing the algorithms
Definition: L1GtTriggerMenuFwd.h:31
relval_2017.k
list k
Definition: relval_2017.py:21
L1GtUtils::availableL1Configuration
const bool availableL1Configuration(int &errorCode, int &l1ConfCode) const
Definition: L1GtUtils.cc:1948
IsoTrackCalib::h_tketa0
TH1I * h_tketa0[6]
Definition: IsoTrackCalib.cc:145
IsoTrackCalib::h_tketa3
TH1I * h_tketa3[6]
Definition: IsoTrackCalib.cc:145
edm::SortedCollection< HBHERecHit >::const_iterator
std::vector< HBHERecHit >::const_iterator const_iterator
Definition: SortedCollection.h:82
edm::EventBase::bunchCrossing
int bunchCrossing() const
Definition: EventBase.h:64
edm::EventBase::luminosityBlock
edm::LuminosityBlockNumber_t luminosityBlock() const
Definition: EventBase.h:61
reco::TrackBase::phi
double phi() const
azimuthal angle of momentum vector
Definition: TrackBase.h:640
L1GtUtils::decision
const bool decision(const edm::Event &iEvent, const std::string &nameAlgoTechTrig, int &errorCode) const
Definition: L1GtUtils.cc:1208
IsoTrackCalib::tok_pfjets_
edm::EDGetTokenT< reco::PFJetCollection > tok_pfjets_
Definition: IsoTrackCalib.cc:125
spr::chargeIsolationCone
double chargeIsolationCone(unsigned int trkIndex, std::vector< spr::propagatedTrackDirection > &trkDirs, double dR, int &nNearTRKs, bool debug=false)
Definition: ChargeIsolation.cc:236
reco::TrackBase::px
double px() const
x coordinate of momentum vector
Definition: TrackBase.h:622
IsoTrackCalib::tok_L1extTauJet_
edm::EDGetTokenT< l1extra::L1JetParticleCollection > tok_L1extTauJet_
Definition: IsoTrackCalib.cc:126
IsoTrackCalib::tok_hbhe_
edm::EDGetTokenT< HBHERecHitCollection > tok_hbhe_
Definition: IsoTrackCalib.cc:122
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
IsoTrackCalib::h_tketa5
TH1I * h_tketa5[6]
Definition: IsoTrackCalib.cc:145
IsoTrackCalib::tok_genTrack_
edm::EDGetTokenT< reco::TrackCollection > tok_genTrack_
Definition: IsoTrackCalib.cc:117
IsoTrackCalib::h_tketa1
TH1I * h_tketa1[6]
Definition: IsoTrackCalib.cc:145
CItAlgo
AlgorithmMap::const_iterator CItAlgo
iterators through map containing the algorithms
Definition: L1GtTriggerMenuFwd.h:38
IsoTrackCalib::theTrackQuality
std::string theTrackQuality
Definition: IsoTrackCalib.cc:109
math::XYZTLorentzVector
XYZTLorentzVectorD XYZTLorentzVector
Lorentz vector with cylindrical internal representation using pseudorapidity.
Definition: LorentzVector.h:29
IsoTrackCalib::h_tketa4
TH1I * h_tketa4[6]
Definition: IsoTrackCalib.cc:145
reco::TrackBase::eta
double eta() const
pseudorapidity of momentum vector
Definition: TrackBase.h:646
IsoTrackCalib::pbin
std::vector< double > pbin
Definition: IsoTrackCalib.cc:139
useL1GtTriggerMenuLite
static const bool useL1GtTriggerMenuLite(true)
IsoTrackCalib::deltaR
double deltaR(double eta1, double eta2, double phi1, double phi2)
Definition: IsoTrackCalib.cc:717
spr::goodTrack
bool goodTrack(const reco::Track *pTrack, math::XYZPoint leadPV, trackSelectionParameters parameters, bool debug=false)
Definition: TrackSelection.cc:8
IsoTrackCalib::t_l1bits
std::vector< bool > * t_l1bits
Definition: IsoTrackCalib.cc:112
reco::TrackBase::pt
double pt() const
track transverse momentum
Definition: TrackBase.h:616
IsoTrackCalib::h_tketa2
TH1I * h_tketa2[6]
Definition: IsoTrackCalib.cc:145
HcalDetId::ieta
int ieta() const
get the cell ieta
Definition: HcalDetId.h:56
reco::GenJet
Jets made from MC generator particles.
Definition: GenJet.h:24
ecalTB2006H4_GenSimDigiReco_cfg.bField
tuple bField
Definition: ecalTB2006H4_GenSimDigiReco_cfg.py:235
spr::eCone_ecal
double eCone_ecal(const CaloGeometry *geo, edm::Handle< T > &barrelhits, edm::Handle< T > &endcaphits, const GlobalPoint &hpoint1, const GlobalPoint &point1, double dR, const GlobalVector &trackMom, int &nRecHits, double ebThr=-100, double eeThr=-100, double tMin=-500, double tMax=500, bool debug=false)
IsoTrackCalib::tok_L1extFwdJet_
edm::EDGetTokenT< l1extra::L1JetParticleCollection > tok_L1extFwdJet_
Definition: IsoTrackCalib.cc:128
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:75
reco::TrackBase::pz
double pz() const
z coordinate of momentum vector
Definition: TrackBase.h:634
IsoTrackCalib::tok_EE_
edm::EDGetTokenT< EcalRecHitCollection > tok_EE_
Definition: IsoTrackCalib.cc:121
IsoTrackCalib::t_EventWeight
double t_EventWeight
Definition: IsoTrackCalib.cc:133
IsoTrackCalib::m_l1GtUtils
L1GtUtils m_l1GtUtils
Definition: IsoTrackCalib.cc:104
IsoTrackCalib::h_RecHit_iEta
TProfile * h_RecHit_iEta
Definition: IsoTrackCalib.cc:140
math::XYZPoint
XYZPointD XYZPoint
point in space with cartesian internal representation
Definition: Point3D.h:12
edm::EventSetup::get
const T & get() const
Definition: EventSetup.h:56
edm::ESHandle::product
T const * product() const
Definition: ESHandle.h:86
IsoTrackCalib::tok_ew_
edm::EDGetTokenT< GenEventInfoProduct > tok_ew_
Definition: IsoTrackCalib.cc:123
edm::EventBase::id
edm::EventID id() const
Definition: EventBase.h:58
IsoTrackCalib::selectionParameters
spr::trackSelectionParameters selectionParameters
Definition: IsoTrackCalib.cc:108
iErrorCode
int iErrorCode
Definition: IsoTrackCalib.cc:650
L1GtUtils::getL1GtRunCache
void getL1GtRunCache(const edm::Run &, const edm::EventSetup &, const bool, const bool)
get all the run-constant quantities for L1 trigger and cache them
Definition: L1GtUtils.cc:319
IsoTrackCalib::tok_EB_
edm::EDGetTokenT< EcalRecHitCollection > tok_EB_
Definition: IsoTrackCalib.cc:120
L1GtTriggerMenu::gtAlgorithmMap
const AlgorithmMap & gtAlgorithmMap() const
get / set the algorithm map (by name)
Definition: L1GtTriggerMenu.h:241
IsoTrackCalib::tok_bs_
edm::EDGetTokenT< reco::BeamSpot > tok_bs_
Definition: IsoTrackCalib.cc:119
gather_cfg.cout
tuple cout
Definition: gather_cfg.py:145
reco::LeafCandidate::eta
virtual double eta() const final
momentum pseudorapidity
Definition: LeafCandidate.h:137
IsoTrackCalib::h_jetpt
TH1F * h_jetpt[4]
Definition: IsoTrackCalib.cc:144
pfJets_cff.pfJets
tuple pfJets
Definition: pfJets_cff.py:8
reco::TrackBase::py
double py() const
y coordinate of momentum vector
Definition: TrackBase.h:628
useL1EventSetup
static const bool useL1EventSetup(true)
spr::eCone_hcal
double eCone_hcal(const CaloGeometry *geo, edm::Handle< T > &hits, const GlobalPoint &hpoint1, const GlobalPoint &point1, double dR, const GlobalVector &trackMom, int &nRecHits, double hbThr=-100, double heThr=-100, double hfThr=-100, double hoThr=-100, double tMin=-500, double tMax=500, int detOnly=-1, bool useRaw=false, bool debug=false)
IsoTrackCalib::tok_recVtx_
edm::EDGetTokenT< reco::VertexCollection > tok_recVtx_
Definition: IsoTrackCalib.cc:118
reco::LeafCandidate::pt
virtual double pt() const final
transverse momentum
Definition: LeafCandidate.h:131
virtual void IsoTrackCalib::beginJob ( )
privatevirtual

Reimplemented from edm::EDAnalyzer.

void IsoTrackCalib::beginJob ( void  )
privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 577 of file IsoTrackCalib.cc.

References fs, h_iEta, h_jetpt, h_Rechit_E, h_RecHit_iEta, h_RecHit_num, h_tketa0, h_tketa1, h_tketa2, h_tketa3, h_tketa4, h_tketa5, relval_2017::k, TFileService::make(), mergeVDriftHistosByStation::name, pbin, AlCaHLTBitMon_QueryRunRegistry::string, t_DetIds, t_eHcal, t_eMipDR, t_Event, t_EventWeight, t_HitEnergies, t_hmaxNearP, t_ieta, t_l1bits, t_l1eta, t_l1phi, t_l1pt, t_l3eta, t_l3phi, t_l3pt, t_mindR1, t_mindR2, t_p, t_qltyFlag, t_qltyMissFlag, t_qltyPVFlag, t_Run, t_selectTk, SiPixelPhase1TrackClusters_cfi::title, and tree.

577  {
578  h_RecHit_iEta = fs->make<TProfile>("rechit_ieta","Rec hit vs. ieta",60,-30,30,0,1000);
579  h_RecHit_num = fs->make<TProfile>("rechit_num","Rec hit vs. num",100,0,20,0,1000);
580  h_iEta = fs->make<TH1I>("iEta","iEta",60,-30,30);
581  h_Rechit_E = fs->make<TH1F>("Rechit_E","Rechit_E",100,0,1000);
582 
583  double prange[5] = {20,30,40,60,100};
584  for (int k=0; k<5; ++k) pbin.push_back(prange[k]);
585  std::string type[6] = {"All", "Trigger OK", "Tree Selected",
586  "Charge Isolation", "MIP Cut", "L1 Cut"};
587  for (unsigned int k=0; k<pbin.size(); ++k) {
588  char name[20], namp[20], title[100];
589  if (k == 0) sprintf (namp, "all momentum");
590  else sprintf (namp, "p = %4.0f:%4.0f GeV", pbin[k-1], pbin[k]);
591  sprintf (name, "TrackEta0%d", k);
592  sprintf (title, "Track #eta for tracks with %s (%s)",namp,type[0].c_str());
593  h_tketa0[k] = fs->make<TH1I>(name, title, 60, -30, 30);
594  sprintf (name, "TrackEta1%d", k);
595  sprintf (title, "Track #eta for tracks with %s (%s)",namp,type[1].c_str());
596  h_tketa1[k] = fs->make<TH1I>(name, title, 60, -30, 30);
597  sprintf (name, "TrackEta2%d", k);
598  sprintf (title, "Track #eta for tracks with %s (%s)",namp,type[2].c_str());
599  h_tketa2[k] = fs->make<TH1I>(name, title, 60, -30, 30);
600  sprintf (name, "TrackEta3%d", k);
601  sprintf (title, "Track #eta for tracks with %s (%s)",namp,type[3].c_str());
602  h_tketa3[k] = fs->make<TH1I>(name, title, 60, -30, 30);
603  sprintf (name, "TrackEta4%d", k);
604  sprintf (title, "Track #eta for tracks with %s (%s)",namp,type[4].c_str());
605  h_tketa4[k] = fs->make<TH1I>(name, title, 60, -30, 30);
606  sprintf (name, "TrackEta5%d", k);
607  sprintf (title, "Track #eta for tracks with %s (%s)",namp,type[5].c_str());
608  h_tketa5[k] = fs->make<TH1I>(name, title, 60, -30, 30);
609  }
610  h_jetpt[0] = fs->make<TH1F>("Jetpt0","Jet p_T (All)", 500,0.,2500.);
611  h_jetpt[1] = fs->make<TH1F>("Jetpt1","Jet p_T (All Weighted)", 500,0.,2500.);
612  h_jetpt[2] = fs->make<TH1F>("Jetpt2","Jet p_T (|#eta| < 2.5)", 500,0.,2500.);
613  h_jetpt[3] = fs->make<TH1F>("Jetpt3","Jet p_T (|#eta| < 2.5 Weighted)", 500,0.,2500.);
614 
615  tree = fs->make<TTree>("CalibTree", "CalibTree");
616 
617  tree->Branch("t_Run", &t_Run, "t_Run/I");
618  tree->Branch("t_Event", &t_Event, "t_Event/I");
619  tree->Branch("t_ieta", &t_ieta, "t_ieta/I");
620  tree->Branch("t_EventWeight", &t_EventWeight, "t_EventWeight/D");
621  tree->Branch("t_l1pt", &t_l1pt, "t_l1pt/D");
622  tree->Branch("t_l1eta", &t_l1eta, "t_l1eta/D");
623  tree->Branch("t_l1phi", &t_l1phi, "t_l1phi/D");
624  tree->Branch("t_l3pt", &t_l3pt, "t_l3pt/D");
625  tree->Branch("t_l3eta", &t_l3eta, "t_l3eta/D");
626  tree->Branch("t_l3phi", &t_l3phi, "t_l3phi/D");
627  tree->Branch("t_p", &t_p, "t_p/D");
628  tree->Branch("t_mindR1", &t_mindR1, "t_mindR1/D");
629  tree->Branch("t_mindR2", &t_mindR2, "t_mindR2/D");
630  tree->Branch("t_eMipDR", &t_eMipDR, "t_eMipDR/D");
631  tree->Branch("t_eHcal", &t_eHcal, "t_eHcal/D");
632  tree->Branch("t_hmaxNearP", &t_hmaxNearP, "t_hmaxNearP/D");
633  tree->Branch("t_selectTk", &t_selectTk, "t_selectTk/O");
634  tree->Branch("t_qltyFlag", &t_qltyFlag, "t_qltyFlag/O");
635  tree->Branch("t_qltyMissFlag",&t_qltyMissFlag,"t_qltyMissFlag/O");
636  tree->Branch("t_qltyPVFlag", &t_qltyPVFlag, "t_qltyPVFlag/O)");
637 
638  t_DetIds = new std::vector<unsigned int>();
639  t_HitEnergies = new std::vector<double>();
640  t_l1bits = new std::vector<bool>();
641  tree->Branch("t_DetIds", "std::vector<unsigned int>", &t_DetIds);
642  tree->Branch("t_HitEnergies", "std::vector<double>", &t_HitEnergies);
643  tree->Branch("t_l1bits","std::vector<bool>", &t_l1bits);
644 }
type
type
Definition: HCALResponse.h:21
IsoTrackCalib::t_HitEnergies
std::vector< double > * t_HitEnergies
Definition: IsoTrackCalib.cc:139
IsoTrackCalib::t_DetIds
std::vector< unsigned int > * t_DetIds
Definition: IsoTrackCalib.cc:138
IsoTrackCalib::h_RecHit_num
TProfile * h_RecHit_num
Definition: IsoTrackCalib.cc:140
relval_2017.k
list k
Definition: relval_2017.py:21
IsoTrackCalib::h_tketa0
TH1I * h_tketa0[6]
Definition: IsoTrackCalib.cc:145
IsoTrackCalib::h_tketa3
TH1I * h_tketa3[6]
Definition: IsoTrackCalib.cc:145
TFileService::make
T * make(const Args &...args) const
make new ROOT object
Definition: TFileService.h:64
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
IsoTrackCalib::h_tketa5
TH1I * h_tketa5[6]
Definition: IsoTrackCalib.cc:145
IsoTrackCalib::h_tketa1
TH1I * h_tketa1[6]
Definition: IsoTrackCalib.cc:145
IsoTrackCalib::h_tketa4
TH1I * h_tketa4[6]
Definition: IsoTrackCalib.cc:145
IsoTrackCalib::pbin
std::vector< double > pbin
Definition: IsoTrackCalib.cc:139
IsoTrackCalib::fs
edm::Service< TFileService > fs
Definition: IsoTrackCalib.cc:102
IsoTrackCalib::t_l1bits
std::vector< bool > * t_l1bits
Definition: IsoTrackCalib.cc:112
IsoTrackCalib::h_tketa2
TH1I * h_tketa2[6]
Definition: IsoTrackCalib.cc:145
IsoTrackCalib::t_EventWeight
double t_EventWeight
Definition: IsoTrackCalib.cc:133
IsoTrackCalib::h_RecHit_iEta
TProfile * h_RecHit_iEta
Definition: IsoTrackCalib.cc:140
SiPixelPhase1TrackClusters_cfi.title
string title
Definition: SiPixelPhase1TrackClusters_cfi.py:75
IsoTrackCalib::h_jetpt
TH1F * h_jetpt[4]
Definition: IsoTrackCalib.cc:144
virtual void IsoTrackCalib::beginLuminosityBlock ( edm::LuminosityBlock const &  ,
edm::EventSetup const &   
)
privatevirtual

Reimplemented from edm::EDAnalyzer.

void IsoTrackCalib::beginLuminosityBlock ( edm::LuminosityBlock const &  ,
edm::EventSetup const &   
)
privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 683 of file IsoTrackCalib.cc.

683 {}
virtual void IsoTrackCalib::beginRun ( edm::Run const &  ,
edm::EventSetup const &   
)
privatevirtual

Reimplemented from edm::EDAnalyzer.

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

Reimplemented from edm::EDAnalyzer.

Definition at line 652 of file IsoTrackCalib.cc.

References changed, gather_cfg::cout, plotBeamSpotDB::first, L1GtTriggerMenu::gtAlgorithmMap(), L1GtTriggerMenu::gtTriggerMenuName(), hltConfig_, iErrorCode, HLTConfigProvider::init(), l1AlgoMap, m_l1GtMenu, m_l1GtUtils, convertSQLiteXML::ok, L1GtUtils::ptrL1TriggerMenuEventSetup(), edm::RunBase::run(), AlCaHLTBitMon_QueryRunRegistry::string, and verbosity.

652  {
653  bool changed = false;
654  bool ok = hltConfig_.init(iRun,iSetup,"HLT",changed);
655  edm::LogInfo("IsoTrack") << "Run " << iRun.run() << " hltconfig.init " << ok;
656 
657  m_l1GtMenu = m_l1GtUtils.ptrL1TriggerMenuEventSetup(iErrorCode);
658  const AlgorithmMap& algorithmMap = m_l1GtMenu->gtAlgorithmMap();
659  const std::string& menuName = m_l1GtMenu->gtTriggerMenuName();
660 
661  // if (verbosity%10>0) edm::LogInfo("IsoTrack") << "menuName " << menuName;
662  std::cout << "menuName " << menuName << std::endl;;
663  for (CItAlgo itAlgo = algorithmMap.begin(); itAlgo != algorithmMap.end();
664  itAlgo++) {
665  std::string algName = itAlgo->first;
666  int algBitNumber = ( itAlgo->second ).algoBitNumber();
667  l1AlgoMap.insert( std::pair<std::pair<unsigned int,std::string>,int>( std::pair<unsigned int,std::string>(algBitNumber, algName) , 0) ) ;
668  }
669  std::map< std::pair<unsigned int,std::string>, int>::iterator itr;
670  for (itr=l1AlgoMap.begin(); itr!=l1AlgoMap.end(); itr++) {
671  if (verbosity%10>0)
672  edm::LogInfo("IsoTrack") << " ********** " << (itr->first).first << " "
673  << (itr->first).second <<" "<<itr->second;
674  }
675 }
L1GtTriggerMenu::gtTriggerMenuName
const std::string & gtTriggerMenuName() const
Definition: L1GtTriggerMenu.h:99
IsoTrackCalib::m_l1GtMenu
const L1GtTriggerMenu * m_l1GtMenu
Definition: IsoTrackCalib.cc:111
IsoTrackCalib::hltConfig_
HLTConfigProvider hltConfig_
Definition: IsoTrackCalib.cc:103
AlgorithmMap
std::map< std::string, L1GtAlgorithm > AlgorithmMap
map containing the algorithms
Definition: L1GtTriggerMenuFwd.h:31
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
CItAlgo
AlgorithmMap::const_iterator CItAlgo
iterators through map containing the algorithms
Definition: L1GtTriggerMenuFwd.h:38
IsoTrackCalib::l1AlgoMap
std::map< std::pair< unsigned int, std::string >, int > l1AlgoMap
Definition: IsoTrackCalib.cc:147
IsoTrackCalib::m_l1GtUtils
L1GtUtils m_l1GtUtils
Definition: IsoTrackCalib.cc:104
HLTConfigProvider::init
bool init(const edm::Run &iRun, const edm::EventSetup &iSetup, const std::string &processName, bool &changed)
d&#39;tor
Definition: HLTConfigProvider.cc:38
iErrorCode
int iErrorCode
Definition: IsoTrackCalib.cc:650
L1GtTriggerMenu::gtAlgorithmMap
const AlgorithmMap & gtAlgorithmMap() const
get / set the algorithm map (by name)
Definition: L1GtTriggerMenu.h:241
gather_cfg.cout
tuple cout
Definition: gather_cfg.py:145
L1GtUtils::ptrL1TriggerMenuEventSetup
const L1GtTriggerMenu * ptrL1TriggerMenuEventSetup(int &errorCode)
return a pointer to the L1 trigger menu from event setup
Definition: L1GtUtils.cc:1883
void IsoTrackCalib::clearTreeVectors ( )
private
double IsoTrackCalib::deltaR ( double  eta1,
double  eta2,
double  phi1,
double  phi2 
)
private

Definition at line 717 of file IsoTrackCalib.cc.

References M_PI, and mathSSE::sqrt().

Referenced by analyze().

717  {
718  double deta = eta1-eta2;
719  if (phi1 < 0) phi1 += 2.0*M_PI;
720  if (phi2 < 0) phi2 += 2.0*M_PI;
721  double dphi = phi1-phi2;
722  if (dphi > M_PI) dphi -= 2.*M_PI;
723  else if (dphi < -M_PI) dphi += 2.*M_PI;
724  return std::sqrt(deta*deta + dphi*dphi);
725 }
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:18
M_PI
#define M_PI
Definition: BXVectorInputProducer.cc:51
double IsoTrackCalib::dEta ( math::XYZTLorentzVector &  ,
math::XYZTLorentzVector &   
)
private
double IsoTrackCalib::dEta ( math::XYZTLorentzVector &  vec1,
math::XYZTLorentzVector &  vec2 
)
private

Definition at line 695 of file IsoTrackCalib.cc.

Referenced by dR().

695  {
696  return (vec1.eta()-vec2.eta());
697 }
vec1
std::vector< double > vec1
Definition: HCALResponse.h:15
vec2
std::vector< vec1 > vec2
Definition: HCALResponse.h:16
double IsoTrackCalib::dinvPt ( math::XYZTLorentzVector &  ,
math::XYZTLorentzVector &   
)
private
double IsoTrackCalib::dP ( math::XYZTLorentzVector &  ,
math::XYZTLorentzVector &   
)
private
double IsoTrackCalib::dPhi ( math::XYZTLorentzVector &  ,
math::XYZTLorentzVector &   
)
private
double IsoTrackCalib::dPhi ( math::XYZTLorentzVector &  vec1,
math::XYZTLorentzVector &  vec2 
)
private

Definition at line 699 of file IsoTrackCalib.cc.

References M_PI.

Referenced by dR().

699  {
700 
701  double phi1 = vec1.phi();
702  if (phi1 < 0) phi1 += 2.0*M_PI;
703  double phi2 = vec2.phi();
704  if (phi2 < 0) phi2 += 2.0*M_PI;
705  double dphi = phi1-phi2;
706  if (dphi > M_PI) dphi -= 2.*M_PI;
707  else if (dphi < -M_PI) dphi += 2.*M_PI;
708  return dphi;
709 }
vec1
std::vector< double > vec1
Definition: HCALResponse.h:15
M_PI
#define M_PI
Definition: BXVectorInputProducer.cc:51
vec2
std::vector< vec1 > vec2
Definition: HCALResponse.h:16
double IsoTrackCalib::dPt ( math::XYZTLorentzVector &  ,
math::XYZTLorentzVector &   
)
private
double IsoTrackCalib::dR ( double  eta1,
double  eta2,
double  phi1,
double  phi2 
)
double IsoTrackCalib::dR ( double  eta1,
double  eta2,
double  phi1,
double  phi2 
)
double IsoTrackCalib::dR ( math::XYZTLorentzVector &  ,
math::XYZTLorentzVector &   
)
private
double IsoTrackCalib::dR ( math::XYZTLorentzVector &  vec1,
math::XYZTLorentzVector &  vec2 
)
private

Definition at line 711 of file IsoTrackCalib.cc.

References dEta(), dPhi(), and mathSSE::sqrt().

711  {
712  double deta = dEta(vec1,vec2);
713  double dphi = dPhi(vec1,vec2);
714  return std::sqrt(deta*deta + dphi*dphi);
715 }
vec1
std::vector< double > vec1
Definition: HCALResponse.h:15
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:18
IsoTrackCalib::dEta
double dEta(math::XYZTLorentzVector &, math::XYZTLorentzVector &)
Definition: IsoTrackCalib.cc:695
IsoTrackCalib::dPhi
double dPhi(math::XYZTLorentzVector &, math::XYZTLorentzVector &)
Definition: IsoTrackCalib.cc:699
vec2
std::vector< vec1 > vec2
Definition: HCALResponse.h:16
virtual void IsoTrackCalib::endJob ( )
privatevirtual

Reimplemented from edm::EDAnalyzer.

Referenced by o2o.O2ORunMgr::executeJob().

void IsoTrackCalib::endJob ( void  )
privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 647 of file IsoTrackCalib.cc.

Referenced by o2o.O2ORunMgr::executeJob().

647  {
648 }
virtual void IsoTrackCalib::endLuminosityBlock ( edm::LuminosityBlock const &  ,
edm::EventSetup const &   
)
privatevirtual

Reimplemented from edm::EDAnalyzer.

void IsoTrackCalib::endLuminosityBlock ( edm::LuminosityBlock const &  ,
edm::EventSetup const &   
)
privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 685 of file IsoTrackCalib.cc.

685 {}
virtual void IsoTrackCalib::endRun ( edm::Run const &  ,
edm::EventSetup const &   
)
privatevirtual

Reimplemented from edm::EDAnalyzer.

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

Reimplemented from edm::EDAnalyzer.

Definition at line 678 of file IsoTrackCalib.cc.

References edm::RunBase::run().

678  {
679  edm::LogInfo("IsoTrack") << "endRun " << iRun.run() << std::endl;
680 }
static void IsoTrackCalib::fillDescriptions ( edm::ConfigurationDescriptions &  descriptions)
static
void IsoTrackCalib::fillDescriptions ( edm::ConfigurationDescriptions &  descriptions)
static

Definition at line 687 of file IsoTrackCalib.cc.

References edm::ConfigurationDescriptions::addDefault(), and edm::ParameterSetDescription::setUnknown().

687  {
688  //The following says we do not know what parameters are allowed so do no validation
689  // Please change this to state exactly what you do use, even if it is no parameters
690  edm::ParameterSetDescription desc;
691  desc.setUnknown();
692  descriptions.addDefault(desc);
693 }
edm::ConfigurationDescriptions::addDefault
void addDefault(ParameterSetDescription const &psetDescription)
Definition: ConfigurationDescriptions.cc:99

Member Data Documentation

double IsoTrackCalib::a_charIsoR
private

Definition at line 110 of file IsoTrackCalib.cc.

Referenced by analyze(), and IsoTrackCalib().

double IsoTrackCalib::a_coneR
private

Definition at line 110 of file IsoTrackCalib.cc.

Referenced by analyze(), and IsoTrackCalib().

double IsoTrackCalib::a_mipR
private

Definition at line 110 of file IsoTrackCalib.cc.

Referenced by analyze(), and IsoTrackCalib().

double IsoTrackCalib::a_neutIsoR
private

Definition at line 103 of file IsoTrackCalib.h.

double IsoTrackCalib::a_neutR1
private

Definition at line 104 of file IsoTrackCalib.h.

double IsoTrackCalib::a_neutR2
private

Definition at line 104 of file IsoTrackCalib.h.

bool IsoTrackCalib::changed
private

Definition at line 96 of file IsoTrackCalib.h.

Referenced by beginRun().

double IsoTrackCalib::cutCharge
private

Definition at line 104 of file IsoTrackCalib.h.

double IsoTrackCalib::cutMip
private

Definition at line 104 of file IsoTrackCalib.h.

double IsoTrackCalib::cutNeutral
private

Definition at line 104 of file IsoTrackCalib.h.

double IsoTrackCalib::dr_L1
private

Definition at line 103 of file IsoTrackCalib.h.

std::vector<double> IsoTrackCalib::drCuts
private

Definition at line 106 of file IsoTrackCalib.h.

int IsoTrackCalib::Event
private

Definition at line 130 of file IsoTrackCalib.h.

bool IsoTrackCalib::firstEvent
private

Definition at line 122 of file IsoTrackCalib.h.

Referenced by looper.Looper::loop(), and Printer.Printer::process().

edm::Service< TFileService > IsoTrackCalib::fs
private

Definition at line 102 of file IsoTrackCalib.cc.

Referenced by beginJob().

TH1I * IsoTrackCalib::g_Accepts
private

Definition at line 127 of file IsoTrackCalib.h.

TH1I* IsoTrackCalib::g_Pre
private

Definition at line 127 of file IsoTrackCalib.h.

TH1I * IsoTrackCalib::g_PreHLT
private

Definition at line 127 of file IsoTrackCalib.h.

TH1I * IsoTrackCalib::g_PreL1
private

Definition at line 127 of file IsoTrackCalib.h.

TH1I * IsoTrackCalib::h_HLTAccept
private

Definition at line 125 of file IsoTrackCalib.h.

std::vector<TH1I*> IsoTrackCalib::h_HLTAccepts
private

Definition at line 126 of file IsoTrackCalib.h.

TH1I* IsoTrackCalib::h_iEta
private

Definition at line 141 of file IsoTrackCalib.cc.

Referenced by analyze(), and beginJob().

TH1F * IsoTrackCalib::h_jetp
private

Definition at line 143 of file IsoTrackCalib.cc.

TH1F* IsoTrackCalib::h_jetpt[4]
private

Definition at line 144 of file IsoTrackCalib.cc.

Referenced by analyze(), and beginJob().

TH1I* IsoTrackCalib::h_nHLT
private

Definition at line 125 of file IsoTrackCalib.h.

TH1F* IsoTrackCalib::h_Rechit_E
private

Definition at line 143 of file IsoTrackCalib.cc.

Referenced by analyze(), and beginJob().

TProfile* IsoTrackCalib::h_RecHit_iEta
private

Definition at line 140 of file IsoTrackCalib.cc.

Referenced by analyze(), and beginJob().

TProfile * IsoTrackCalib::h_RecHit_num
private

Definition at line 140 of file IsoTrackCalib.cc.

Referenced by analyze(), and beginJob().

TH1I * IsoTrackCalib::h_tkEta0[5]
private

Definition at line 141 of file IsoTrackCalib.cc.

TH1I* IsoTrackCalib::h_tketa0[6]
private

Definition at line 145 of file IsoTrackCalib.cc.

Referenced by analyze(), and beginJob().

TH1I * IsoTrackCalib::h_tkEta1[5]
private

Definition at line 141 of file IsoTrackCalib.cc.

TH1I * IsoTrackCalib::h_tketa1[6]
private

Definition at line 145 of file IsoTrackCalib.cc.

Referenced by analyze(), and beginJob().

TH1I * IsoTrackCalib::h_tkEta2[5]
private

Definition at line 141 of file IsoTrackCalib.cc.

TH1I * IsoTrackCalib::h_tketa2[6]
private

Definition at line 145 of file IsoTrackCalib.cc.

Referenced by analyze(), and beginJob().

TH1I* IsoTrackCalib::h_tkEta3[5]
private

Definition at line 142 of file IsoTrackCalib.cc.

TH1I * IsoTrackCalib::h_tketa3[6]
private

Definition at line 145 of file IsoTrackCalib.cc.

Referenced by analyze(), and beginJob().

TH1I * IsoTrackCalib::h_tkEta4[5]
private

Definition at line 142 of file IsoTrackCalib.cc.

TH1I * IsoTrackCalib::h_tketa4[6]
private

Definition at line 145 of file IsoTrackCalib.cc.

Referenced by analyze(), and beginJob().

TH1I * IsoTrackCalib::h_tkEta5[5]
private

Definition at line 142 of file IsoTrackCalib.cc.

TH1I * IsoTrackCalib::h_tketa5[6]
private

Definition at line 145 of file IsoTrackCalib.cc.

Referenced by analyze(), and beginJob().

HLTConfigProvider IsoTrackCalib::hltConfig_
private

Definition at line 103 of file IsoTrackCalib.cc.

Referenced by beginRun().

std::vector< std::string > IsoTrackCalib::HLTNames
private

Definition at line 106 of file IsoTrackCalib.cc.

HLTPrescaleProvider IsoTrackCalib::hltPrescaleProvider_
private

Definition at line 98 of file IsoTrackCalib.h.

bool IsoTrackCalib::initL1
private

Definition at line 105 of file IsoTrackCalib.cc.

std::map< std::pair<unsigned int,std::string>, int> IsoTrackCalib::l1AlgoMap
private

Definition at line 147 of file IsoTrackCalib.cc.

Referenced by beginRun().

std::vector<std::string> IsoTrackCalib::l1Names
private

Definition at line 106 of file IsoTrackCalib.cc.

Referenced by analyze(), and IsoTrackCalib().

const L1GtTriggerMenu* IsoTrackCalib::m_l1GtMenu
private

Definition at line 111 of file IsoTrackCalib.cc.

Referenced by analyze(), and beginRun().

L1GtUtils IsoTrackCalib::m_l1GtUtils
private

Definition at line 104 of file IsoTrackCalib.cc.

Referenced by analyze(), and beginRun().

int IsoTrackCalib::maxRunNo
private

Definition at line 105 of file IsoTrackCalib.h.

int IsoTrackCalib::minRunNo
private

Definition at line 105 of file IsoTrackCalib.h.

int IsoTrackCalib::nRun
private

Definition at line 105 of file IsoTrackCalib.h.

std::vector<double> IsoTrackCalib::pbin
private

Definition at line 139 of file IsoTrackCalib.cc.

Referenced by analyze(), and beginJob().

int IsoTrackCalib::Run
private

Definition at line 130 of file IsoTrackCalib.h.

Referenced by BeamSpotObj.BeamSpot::Reset().

spr::trackSelectionParameters IsoTrackCalib::selectionParameters
private

Definition at line 108 of file IsoTrackCalib.cc.

Referenced by analyze(), and IsoTrackCalib().

std::vector<double>* IsoTrackCalib::t_charge_iso
private

Definition at line 133 of file IsoTrackCalib.h.

std::vector<unsigned int>* IsoTrackCalib::t_DetIds
private

Definition at line 138 of file IsoTrackCalib.cc.

Referenced by analyze(), and beginJob().

std::vector<double>* IsoTrackCalib::t_disthotcell
private

Definition at line 135 of file IsoTrackCalib.h.

std::vector<double> * IsoTrackCalib::t_ehcal
private

Definition at line 133 of file IsoTrackCalib.h.

double IsoTrackCalib::t_eHcal
private

Definition at line 135 of file IsoTrackCalib.cc.

Referenced by analyze(), and beginJob().

std::vector<double> * IsoTrackCalib::t_emip
private

Definition at line 133 of file IsoTrackCalib.h.

double IsoTrackCalib::t_eMipDR
private

Definition at line 135 of file IsoTrackCalib.cc.

Referenced by analyze(), and beginJob().

int IsoTrackCalib::t_Event
private

Definition at line 132 of file IsoTrackCalib.cc.

Referenced by analyze(), and beginJob().

double IsoTrackCalib::t_EventWeight
private

Definition at line 133 of file IsoTrackCalib.cc.

Referenced by analyze(), and beginJob().

std::vector<double> * IsoTrackCalib::t_eventweight
private

Definition at line 135 of file IsoTrackCalib.h.

std::vector<double>* IsoTrackCalib::t_HitEnergies
private

Definition at line 139 of file IsoTrackCalib.cc.

Referenced by analyze(), and beginJob().

double IsoTrackCalib::t_hmaxNearP
private

Definition at line 135 of file IsoTrackCalib.cc.

Referenced by analyze(), and beginJob().

int IsoTrackCalib::t_ieta
private

Definition at line 132 of file IsoTrackCalib.cc.

Referenced by analyze(), and beginJob().

std::vector<int>* IsoTrackCalib::t_ieta
private

Definition at line 134 of file IsoTrackCalib.h.

std::vector<double> * IsoTrackCalib::t_ietahotcell
private

Definition at line 135 of file IsoTrackCalib.h.

std::vector<bool>* IsoTrackCalib::t_l1bits
private

Definition at line 112 of file IsoTrackCalib.cc.

Referenced by analyze(), and beginJob().

double IsoTrackCalib::t_l1eta
private

Definition at line 133 of file IsoTrackCalib.cc.

Referenced by analyze(), and beginJob().

std::vector<double> * IsoTrackCalib::t_l1eta
private

Definition at line 136 of file IsoTrackCalib.h.

double IsoTrackCalib::t_l1phi
private

Definition at line 133 of file IsoTrackCalib.cc.

Referenced by analyze(), and beginJob().

std::vector<double> * IsoTrackCalib::t_l1phi
private

Definition at line 136 of file IsoTrackCalib.h.

double IsoTrackCalib::t_l1pt
private

Definition at line 133 of file IsoTrackCalib.cc.

Referenced by analyze(), and beginJob().

std::vector<double>* IsoTrackCalib::t_l1pt
private

Definition at line 136 of file IsoTrackCalib.h.

double IsoTrackCalib::t_l3eta
private

Definition at line 134 of file IsoTrackCalib.cc.

Referenced by analyze(), and beginJob().

std::vector<double> * IsoTrackCalib::t_l3eta
private

Definition at line 137 of file IsoTrackCalib.h.

double IsoTrackCalib::t_l3phi
private

Definition at line 134 of file IsoTrackCalib.cc.

Referenced by analyze(), and beginJob().

std::vector<double> * IsoTrackCalib::t_l3phi
private

Definition at line 137 of file IsoTrackCalib.h.

double IsoTrackCalib::t_l3pt
private

Definition at line 134 of file IsoTrackCalib.cc.

Referenced by analyze(), and beginJob().

std::vector<double>* IsoTrackCalib::t_l3pt
private

Definition at line 137 of file IsoTrackCalib.h.

std::vector<double> * IsoTrackCalib::t_leadingeta
private

Definition at line 138 of file IsoTrackCalib.h.

std::vector<double> * IsoTrackCalib::t_leadingphi
private

Definition at line 138 of file IsoTrackCalib.h.

std::vector<double>* IsoTrackCalib::t_leadingpt
private

Definition at line 138 of file IsoTrackCalib.h.

double IsoTrackCalib::t_mindR1
private

Definition at line 134 of file IsoTrackCalib.cc.

Referenced by analyze(), and beginJob().

double IsoTrackCalib::t_mindR2
private

Definition at line 135 of file IsoTrackCalib.cc.

Referenced by analyze(), and beginJob().

std::vector<double> * IsoTrackCalib::t_neu_iso
private

Definition at line 132 of file IsoTrackCalib.h.

double IsoTrackCalib::t_p
private

Definition at line 134 of file IsoTrackCalib.cc.

Referenced by analyze(), and beginJob().

bool IsoTrackCalib::t_qltyFlag
private

Definition at line 137 of file IsoTrackCalib.cc.

Referenced by analyze(), and beginJob().

bool IsoTrackCalib::t_qltyMissFlag
private

Definition at line 137 of file IsoTrackCalib.cc.

Referenced by analyze(), and beginJob().

bool IsoTrackCalib::t_qltyPVFlag
private

Definition at line 137 of file IsoTrackCalib.cc.

Referenced by analyze(), and beginJob().

int IsoTrackCalib::t_Run
private

Definition at line 132 of file IsoTrackCalib.cc.

Referenced by analyze(), and beginJob().

bool IsoTrackCalib::t_selectTk
private

Definition at line 137 of file IsoTrackCalib.cc.

Referenced by analyze(), and beginJob().

std::vector<double>* IsoTrackCalib::t_trackEta
private

Definition at line 132 of file IsoTrackCalib.h.

std::vector<double>* IsoTrackCalib::t_trackP
private

Definition at line 131 of file IsoTrackCalib.h.

std::vector<double> * IsoTrackCalib::t_trackPhi
private

Definition at line 132 of file IsoTrackCalib.h.

std::vector<double> * IsoTrackCalib::t_trackPt
private

Definition at line 132 of file IsoTrackCalib.h.

std::vector<double> * IsoTrackCalib::t_trackPx
private

Definition at line 131 of file IsoTrackCalib.h.

std::vector<double> * IsoTrackCalib::t_trackPy
private

Definition at line 131 of file IsoTrackCalib.h.

std::vector<double> * IsoTrackCalib::t_trackPz
private

Definition at line 131 of file IsoTrackCalib.h.

std::vector<double> * IsoTrackCalib::t_trkL3mindr
private

Definition at line 133 of file IsoTrackCalib.h.

std::string IsoTrackCalib::theTrackQuality
private

Definition at line 109 of file IsoTrackCalib.cc.

Referenced by analyze(), and IsoTrackCalib().

edm::InputTag IsoTrackCalib::theTriggerResultsLabel
private

Definition at line 113 of file IsoTrackCalib.cc.

Referenced by IsoTrackCalib().

edm::EDGetTokenT< reco::BeamSpot > IsoTrackCalib::tok_bs_
private

Definition at line 119 of file IsoTrackCalib.cc.

Referenced by analyze(), and IsoTrackCalib().

edm::EDGetTokenT< EcalRecHitCollection > IsoTrackCalib::tok_EB_
private

Definition at line 120 of file IsoTrackCalib.cc.

Referenced by analyze(), and IsoTrackCalib().

edm::EDGetTokenT< EcalRecHitCollection > IsoTrackCalib::tok_EE_
private

Definition at line 121 of file IsoTrackCalib.cc.

Referenced by analyze(), and IsoTrackCalib().

edm::EDGetTokenT< GenEventInfoProduct > IsoTrackCalib::tok_ew_
private

Definition at line 123 of file IsoTrackCalib.cc.

Referenced by analyze(), and IsoTrackCalib().

edm::EDGetTokenT< reco::TrackCollection > IsoTrackCalib::tok_genTrack_
private

Definition at line 117 of file IsoTrackCalib.cc.

Referenced by analyze(), and IsoTrackCalib().

edm::EDGetTokenT< HBHERecHitCollection > IsoTrackCalib::tok_hbhe_
private

Definition at line 122 of file IsoTrackCalib.cc.

Referenced by analyze(), and IsoTrackCalib().

edm::EDGetTokenT<reco::GenJetCollection> IsoTrackCalib::tok_jets_
private

Definition at line 124 of file IsoTrackCalib.cc.

Referenced by analyze(), and IsoTrackCalib().

edm::EDGetTokenT<l1extra::L1JetParticleCollection> IsoTrackCalib::tok_L1extCenJet_
private

Definition at line 127 of file IsoTrackCalib.cc.

Referenced by analyze(), and IsoTrackCalib().

edm::EDGetTokenT<l1extra::L1JetParticleCollection> IsoTrackCalib::tok_L1extFwdJet_
private

Definition at line 128 of file IsoTrackCalib.cc.

Referenced by analyze(), and IsoTrackCalib().

edm::EDGetTokenT<l1extra::L1JetParticleCollection> IsoTrackCalib::tok_L1extTauJet_
private

Definition at line 126 of file IsoTrackCalib.cc.

Referenced by analyze(), and IsoTrackCalib().

edm::EDGetTokenT<LumiDetails> IsoTrackCalib::tok_lumi
private

Definition at line 109 of file IsoTrackCalib.h.

edm::EDGetTokenT<reco::PFJetCollection> IsoTrackCalib::tok_pf_
private

Definition at line 120 of file IsoTrackCalib.h.

edm::EDGetTokenT<reco::PFJetCollection> IsoTrackCalib::tok_pfjets_
private

Definition at line 125 of file IsoTrackCalib.cc.

Referenced by analyze(), and IsoTrackCalib().

edm::EDGetTokenT< reco::VertexCollection > IsoTrackCalib::tok_recVtx_
private

Definition at line 118 of file IsoTrackCalib.cc.

Referenced by analyze(), and IsoTrackCalib().

edm::EDGetTokenT< trigger::TriggerEvent > IsoTrackCalib::tok_trigEvt
private

Definition at line 114 of file IsoTrackCalib.cc.

Referenced by IsoTrackCalib().

edm::EDGetTokenT< edm::TriggerResults > IsoTrackCalib::tok_trigRes
private

Definition at line 115 of file IsoTrackCalib.cc.

Referenced by IsoTrackCalib().

TTree * IsoTrackCalib::tree
private

Definition at line 131 of file IsoTrackCalib.cc.

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

std::vector<std::string> IsoTrackCalib::trgnames
private

Definition at line 146 of file IsoTrackCalib.cc.

edm::InputTag IsoTrackCalib::triggerEvent_
private

Definition at line 113 of file IsoTrackCalib.cc.

Referenced by IsoTrackCalib().

std::map<std::pair<unsigned int, std::string>, unsigned int> IsoTrackCalib::TrigList
private

Definition at line 123 of file IsoTrackCalib.h.

std::vector< std::string > IsoTrackCalib::trigNames
private

Definition at line 106 of file IsoTrackCalib.cc.

Referenced by IsoTrackCalib().

std::map<std::pair<unsigned int, std::string>, const std::pair<int, int> > IsoTrackCalib::TrigPreList
private

Definition at line 124 of file IsoTrackCalib.h.

int IsoTrackCalib::verbosity
private

Definition at line 107 of file IsoTrackCalib.cc.

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

Generated for CMSSW Reference Manual by doxygen 1.8.5