IsoTrackCalib Class Reference

Inheritance diagram for IsoTrackCalib:

Public Member Functions
	IsoTrackCalib (const edm::ParameterSet &)
	~IsoTrackCalib () override
Public Member Functions inherited from edm::one::EDAnalyzer< edm::one::WatchRuns, edm::one::SharedResources >
	EDAnalyzer ()=default
bool	wantsGlobalLuminosityBlocks () const final
bool	wantsGlobalRuns () const final
Public Member Functions inherited from edm::one::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzerBase ()
ModuleDescription const &	moduleDescription () const
bool	wantsStreamLuminosityBlocks () const
bool	wantsStreamRuns () const
	~EDAnalyzerBase () override
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...
	EDConsumerBase ()
	EDConsumerBase (EDConsumerBase const &)=delete
	EDConsumerBase (EDConsumerBase &&)=default
ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
EDConsumerBase &	operator= (EDConsumerBase &&)=default
bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
virtual	~EDConsumerBase () noexcept(false)

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

Private Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override
void	beginJob () override
void	beginRun (edm::Run const &, edm::EventSetup const &) override
double	deltaR (double eta1, double eta2, double phi1, double phi2)
double	dEta (math::XYZTLorentzVector &, math::XYZTLorentzVector &)
double	dPhi (math::XYZTLorentzVector &, math::XYZTLorentzVector &)
double	dR (math::XYZTLorentzVector &, math::XYZTLorentzVector &)
void	endJob () override
void	endRun (edm::Run const &, edm::EventSetup const &) override

Private Attributes
double	a_charIsoR_
double	a_coneR_
double	a_mipR_
edm::Service< TFileService >	fs_
TH1I *	h_iEta
TH1F *	h_jetp
TH1F *	h_jetpt [4]
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::map< std::pair< unsigned int, std::string >, int >	l1AlgoMap_
std::vector< std::string >	l1Names_
const L1GtTriggerMenu *	m_l1GtMenu
L1GtUtils	m_l1GtUtils
std::vector< double >	pbin
spr::trackSelectionParameters	selectionParameters_
std::vector< unsigned int > *	t_DetIds
double	t_eHcal
double	t_eMipDR
int	t_Event
double	t_EventWeight
std::vector< double > *	t_HitEnergies
double	t_hmaxNearP
int	t_ieta
std::vector< bool > *	t_l1bits
double	t_l1eta
double	t_l1phi
double	t_l1pt
double	t_l3eta
double	t_l3phi
double	t_l3pt
double	t_mindR1
double	t_mindR2
double	t_p
bool	t_qltyFlag
bool	t_qltyMissFlag
bool	t_qltyPVFlag
int	t_Run
bool	t_selectTk
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< reco::PFJetCollection >	tok_pfjets_
edm::EDGetTokenT< reco::VertexCollection >	tok_recVtx_
edm::EDGetTokenT< trigger::TriggerEvent >	tok_trigEvt
edm::EDGetTokenT< edm::TriggerResults >	tok_trigRes
TTree *	tree
edm::InputTag	triggerEvent_
int	verbosity_

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

Detailed Description

Definition at line 82 of file IsoTrackCalib.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 149 of file IsoTrackCalib.cc.

References a_charIsoR_, a_coneR_, a_mipR_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), gen::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_, and verbosity_.

                                                           :
  m_l1GtUtils(iConfig, consumesCollector(), useL1GtTriggerMenuLite, *this){

  usesResource("TFileService");

 //now do whatever initialization is needed
  verbosity_                          = iConfig.getUntrackedParameter<int>("Verbosity",0);
  l1Names_                = iConfig.getUntrackedParameter<std::vector<std::string> >("L1Seed"); 
  theTrackQuality_                    = iConfig.getUntrackedParameter<std::string>("TrackQuality","highPurity");
  reco::TrackBase::TrackQuality trackQuality_=reco::TrackBase::qualityByName(theTrackQuality_);
  selectionParameters_.minPt          = iConfig.getUntrackedParameter<double>("MinTrackPt", 10.0);
  selectionParameters_.minQuality     = trackQuality_;
  selectionParameters_.maxDxyPV       = iConfig.getUntrackedParameter<double>("MaxDxyPV", 0.2);
  selectionParameters_.maxDzPV        = iConfig.getUntrackedParameter<double>("MaxDzPV",  5.0);
  selectionParameters_.maxChi2        = iConfig.getUntrackedParameter<double>("MaxChi2",  5.0);
  selectionParameters_.maxDpOverP     = iConfig.getUntrackedParameter<double>("MaxDpOverP",  0.1);
  selectionParameters_.minOuterHit    = iConfig.getUntrackedParameter<int>("MinOuterHit", 4);
  selectionParameters_.minLayerCrossed= iConfig.getUntrackedParameter<int>("MinLayerCrossed", 8);
  selectionParameters_.maxInMiss      = iConfig.getUntrackedParameter<int>("MaxInMiss", 0);
  selectionParameters_.maxOutMiss     = iConfig.getUntrackedParameter<int>("MaxOutMiss", 0);
  a_coneR_                            = iConfig.getUntrackedParameter<double>("ConeRadius",34.98);
  a_charIsoR_                         = a_coneR_ + 28.9;
  a_mipR_                             = iConfig.getUntrackedParameter<double>("ConeRadiusMIP",14.0);
  bool isItAOD                        = iConfig.getUntrackedParameter<bool>("IsItAOD", false);
  triggerEvent_                       = edm::InputTag("hltTriggerSummaryAOD","","HLT");
  theTriggerResultsLabel_             = edm::InputTag("TriggerResults","","HLT");
  edm::InputTag L1extraTauJetSource_  = iConfig.getParameter<edm::InputTag>  ("L1extraTauJetSource"   );
  edm::InputTag L1extraCenJetSource_  = iConfig.getParameter<edm::InputTag>  ("L1extraCenJetSource"   );
  edm::InputTag L1extraFwdJetSource_  = iConfig.getParameter<edm::InputTag>  ("L1extraFwdJetSource"   );
  
  // define tokens for access
  tok_trigEvt   = consumes<trigger::TriggerEvent>(triggerEvent_);
  tok_trigRes   = consumes<edm::TriggerResults>(theTriggerResultsLabel_);
  tok_genTrack_ = consumes<reco::TrackCollection>(edm::InputTag("generalTracks"));
  tok_recVtx_   = consumes<reco::VertexCollection>(edm::InputTag("offlinePrimaryVertices"));
  tok_bs_       = consumes<reco::BeamSpot>(edm::InputTag("offlineBeamSpot"));
  tok_ew_       = consumes<GenEventInfoProduct>(edm::InputTag("generatorSmeared")); 
  tok_L1extTauJet_  = consumes<l1extra::L1JetParticleCollection>(L1extraTauJetSource_);
  tok_L1extCenJet_  = consumes<l1extra::L1JetParticleCollection>(L1extraCenJetSource_);
  tok_L1extFwdJet_  = consumes<l1extra::L1JetParticleCollection>(L1extraFwdJetSource_);
  if (isItAOD) {
    tok_EB_     = consumes<EcalRecHitCollection>(edm::InputTag("reducedEcalRecHitsEB"));
    tok_EE_     = consumes<EcalRecHitCollection>(edm::InputTag("reducedEcalRecHitsEE"));
    tok_hbhe_   = consumes<HBHERecHitCollection>(edm::InputTag("reducedHcalRecHits", "hbhereco"));
  } else {
    tok_EB_     = consumes<EcalRecHitCollection>(edm::InputTag("ecalRecHit","EcalRecHitsEB"));
    tok_EE_     = consumes<EcalRecHitCollection>(edm::InputTag("ecalRecHit","EcalRecHitsEE"));
    tok_hbhe_   = consumes<HBHERecHitCollection>(edm::InputTag("hbhereco"));
  }
  tok_jets_     = consumes<reco::GenJetCollection>(iConfig.getParameter<edm::InputTag>("JetSource"));
  tok_pfjets_   = consumes<reco::PFJetCollection>(edm::InputTag("ak5PFJets"));
  if (verbosity_>=0) {
    edm::LogInfo("IsoTrack") 
      <<"Parameters read from config file \n" 
      <<"\t minPt "           << selectionParameters_.minPt   
      <<"\t theTrackQuality " << theTrackQuality_
      <<"\t minQuality "      << selectionParameters_.minQuality
      <<"\t maxDxyPV "        << selectionParameters_.maxDxyPV          
      <<"\t maxDzPV "         << selectionParameters_.maxDzPV          
      <<"\t maxChi2 "         << selectionParameters_.maxChi2          
      <<"\t maxDpOverP "      << selectionParameters_.maxDpOverP
      <<"\t minOuterHit "     << selectionParameters_.minOuterHit
      <<"\t minLayerCrossed " << selectionParameters_.minLayerCrossed
      <<"\t maxInMiss "       << selectionParameters_.maxInMiss
      <<"\t maxOutMiss "      << selectionParameters_.maxOutMiss
      <<"\t a_coneR "         << a_coneR_          
      <<"\t a_charIsoR "      << a_charIsoR_          
      <<"\t a_mipR "          << a_mipR_ 
      <<"\t isItAOD "         << isItAOD;
    edm::LogInfo("IsoTrack") << l1Names_.size() << " triggers to be studied";
    for (unsigned int k=0; k<l1Names_.size(); ++k)
      edm::LogInfo("IsoTrack") << "[" << k << "]: " << l1Names_[k];
  }
 
}

IsoTrackCalib::~IsoTrackCalib ( )

override

Definition at line 225 of file IsoTrackCalib.cc.

                              {
  // do anything here that needs to be done at desctruction time
  // (e.g. close files, deallocate resources etc.)

}

Member Function Documentation

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

overrideprivate

Definition at line 256 of file IsoTrackCalib.cc.

References a_charIsoR_, a_coneR_, a_mipR_, L1GtUtils::availableL1Configuration(), edm::SortedCollection< T, SORT >::begin(), ecalTB2006H4_GenSimDigiReco_cfg::bField, edm::EventBase::bunchCrossing(), spr::chargeIsolationCone(), L1GtUtils::decision(), deltaR(), spr::eCone_ecal(), spr::eCone_hcal(), edm::SortedCollection< T, SORT >::end(), reco::LeafCandidate::eta(), reco::TrackBase::eta(), edm::EventID::event(), particlelevel_cff::genEventInfo, ttbarCategorization_cff::genJets, 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, photonIsolationHIProducer_cfi::hbhe, mps_fire::i, edm::EventBase::id(), electrons_cff::ids, HcalDetId::ieta(), edm::HandleBase::isValid(), gen::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, pfJetBenchmark_cfi::pfJets, reco::TrackBase::phi(), reco::BeamSpot::position(), edm::ESHandle< T >::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_, useL1EventSetup(), useL1GtTriggerMenuLite(), verbosity_, GenEventInfoProduct::weight(), x, y, and z.

                                             {

  t_Run   = iEvent.id().run();
  t_Event = iEvent.id().event();
  if (verbosity_%10 > 0) 
    edm::LogInfo("IsoTrack") 
      << "Run " << t_Run << " Event " << t_Event << " Luminosity " 
      << iEvent.luminosityBlock() << " Bunch " << iEvent.bunchCrossing()
      << " starts ==========";

  //Get magnetic field and ECAL channel status
  edm::ESHandle<MagneticField> bFieldH;
  iSetup.get<IdealMagneticFieldRecord>().get(bFieldH);
  const MagneticField *bField = bFieldH.product();
  
  edm::ESHandle<EcalSeverityLevelAlgo> sevlv;
  iSetup.get<EcalSeverityLevelAlgoRcd>().get(sevlv);

  // get handles to calogeometry and calotopology
  edm::ESHandle<CaloGeometry> pG;
  iSetup.get<CaloGeometryRecord>().get(pG);
  const CaloGeometry* geo = pG.product();
  
  //Get track collection
  edm::Handle<reco::TrackCollection> trkCollection;
  iEvent.getByToken(tok_genTrack_, trkCollection);
  reco::TrackCollection::const_iterator trkItr;
 
  //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();  

  // genJet information
  edm::Handle<reco::GenJetCollection> genJets;
  iEvent.getByToken(tok_jets_, genJets);
  if (genJets.isValid()) {
    for (unsigned iGenJet = 0; iGenJet < genJets->size(); ++iGenJet) {
      const reco::GenJet& genJet = (*genJets) [iGenJet];
      double genJetPt  = genJet.pt();
      double genJetEta = genJet.eta();
      h_jetpt[0]->Fill(genJetPt);
      h_jetpt[1]->Fill(genJetPt,t_EventWeight);
      if (genJetEta>-2.5 && genJetEta<2.5) {
    h_jetpt[2]->Fill(genJetPt);
    h_jetpt[3]->Fill(genJetPt,t_EventWeight);
      }
      break;
    }
  }

  //pf jets 
  edm::Handle<reco::PFJetCollection> pfJets;
  iEvent.getByToken(tok_pfjets_, pfJets);
  reco::PFJetCollection::const_iterator pfItr;  

  //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);
  if (!recVtxs->empty() && !((*recVtxs)[0].isFake())) {
    leadPV = math::XYZPoint( (*recVtxs)[0].x(),(*recVtxs)[0].y(), (*recVtxs)[0].z() );
  } else if (beamSpotH.isValid()) {
    leadPV = beamSpotH->position();
  }
  if (verbosity_>10) {
    if ((verbosity_%100)/10>2)
      edm::LogInfo("IsoTrack") << "Primary Vertex " << leadPV;
    if (beamSpotH.isValid()) edm::LogInfo("IsoTrack") << " Beam Spot " 
                              << beamSpotH->position();
  }
  
  // RecHits
  edm::Handle<EcalRecHitCollection> barrelRecHitsHandle;
  edm::Handle<EcalRecHitCollection> endcapRecHitsHandle;
  iEvent.getByToken(tok_EB_, barrelRecHitsHandle);
  iEvent.getByToken(tok_EE_, endcapRecHitsHandle);
  edm::Handle<HBHERecHitCollection> hbhe;
  iEvent.getByToken(tok_hbhe_, hbhe);
  HBHERecHitCollection::const_iterator rhitItr;

  for (rhitItr=hbhe->begin();rhitItr!=hbhe->end();rhitItr++) {
    double rec_energy = rhitItr->energy();
    int    rec_ieta   = rhitItr->id().ieta();
    int    rec_depth  = rhitItr->id().depth();
    int    rec_zside  = rhitItr->id().zside();
    double num1_1     = rec_zside*(rec_ieta+0.2*(rec_depth-1));
    if (verbosity_%10>0) 
      edm::LogInfo("IsoTrack") << "detid/rechit/ieta/zside/depth/num " << " = "
                   << rhitItr->id() << "/" << rec_energy << "/" 
                   << rec_ieta << "/" << rec_zside << "/" 
                   << rec_depth << "/" << num1_1;
    h_iEta->Fill(rec_ieta);
    h_Rechit_E->Fill(rec_energy);
    h_RecHit_iEta->Fill(rec_ieta,rec_energy);
    h_RecHit_num->Fill(num1_1,rec_energy);
  }
      
  //Propagate tracks to calorimeter surface)
  std::vector<spr::propagatedTrackDirection> trkCaloDirections;
  spr::propagateCALO(trkCollection, geo, bField, theTrackQuality_,
             trkCaloDirections, ((verbosity_/100)%10>2));
  std::vector<spr::propagatedTrackDirection>::const_iterator trkDetItr;
  for (trkDetItr = trkCaloDirections.begin(); 
       trkDetItr != trkCaloDirections.end(); trkDetItr++) {
    if (trkDetItr->okHCAL) {
      HcalDetId detId = (HcalDetId)(trkDetItr->detIdHCAL);
      int tk_ieta = detId.ieta();
      const reco::Track* pTrack = &(*(trkDetItr->trkItr));
      double tk_p = pTrack->p();
      h_tketa0[0]->Fill(tk_ieta);
      for (unsigned int k=1; k<pbin.size(); ++k) {
    if (tk_p >= pbin[k-1] && tk_p < pbin[k]) {
      h_tketa0[k]->Fill(tk_ieta);
      break;
    }
      }
    }
  }
  t_l1bits->clear();
  for (unsigned int i=0; i<l1Names_.size(); ++i) t_l1bits->push_back(false);
  bool useL1EventSetup = true;
  bool useL1GtTriggerMenuLite = true;

  m_l1GtUtils.getL1GtRunCache(iEvent, iSetup, useL1EventSetup, useL1GtTriggerMenuLite);
  int iErrorCode = -1;
  int l1ConfCode = -1;
  const bool l1Conf = m_l1GtUtils.availableL1Configuration(iErrorCode, l1ConfCode);
  if( !l1Conf) {
    edm::LogInfo("IsoTrack") << "\nL1 configuration code:" << l1ConfCode
                 << "\nNo valid L1 trigger configuration available."
                 << "\nSee text above for error code interpretation"
                 << "\nNo return here, in order to test each method"
                 << ", protected against configuration error.";
  }

  const AlgorithmMap& algorithmMap = m_l1GtMenu->gtAlgorithmMap();
  const std::string&  menuName     = m_l1GtMenu->gtTriggerMenuName();
  if (verbosity_%10>0)  
    edm::LogInfo("IsoTrack") << "menuName " << menuName << std::endl;

  std::vector<int> algbits;
  for (CItAlgo itAlgo  = algorithmMap.begin(); itAlgo != algorithmMap.end(); 
       itAlgo++) {
    std::string algName      = itAlgo->first;
    int algBitNumber         = ( itAlgo->second ).algoBitNumber();
    bool decision            = m_l1GtUtils.decision(iEvent, itAlgo->first, iErrorCode);

    bool l1ok(false);
    if (verbosity_%10>0)  
      edm::LogInfo("IsoTrack") << algName <<"  "<< algBitNumber <<  "  " 
                   << decision;
    for (unsigned int i=0; i<l1Names_.size(); ++i) {  
      if (algName.find(l1Names_[i])!=std::string::npos){
    if (verbosity_%10>0)  
      edm::LogInfo("IsoTrack")  << "match found" << " " << algName << "  " 
                    << decision;
    t_l1bits->at(i) = (decision>0);                             
    if (decision > 0) l1ok = true;
      }  
    }
    if (verbosity_%10>0) edm::LogInfo("IsoTrack")  << "l1 ok =" << l1ok;
    
    if(l1ok){
      edm::Handle<l1extra::L1JetParticleCollection> l1TauHandle;
      iEvent.getByToken(tok_L1extTauJet_,l1TauHandle);
      l1extra::L1JetParticleCollection::const_iterator itr;
      double ptTriggered  = -10;
      double etaTriggered = -100;
      double phiTriggered = -100;

      for(itr = l1TauHandle->begin(); itr != l1TauHandle->end(); ++itr) {
        if(itr->pt()>ptTriggered){  
      ptTriggered = itr->pt();
      etaTriggered = itr->eta();
      phiTriggered = itr->phi() ;
    }
    if (verbosity_%10>0) 
      edm::LogInfo("IsoTrack") << "tauJ pt " << itr->pt() << "  eta/phi " 
                   << itr->eta() << " " << itr->phi();
      }
      edm::Handle<l1extra::L1JetParticleCollection> l1CenJetHandle;
      iEvent.getByToken(tok_L1extCenJet_,l1CenJetHandle);
      for (itr = l1CenJetHandle->begin();  itr != l1CenJetHandle->end(); ++itr){
    if (itr->pt()>ptTriggered) {
      ptTriggered = itr->pt();
      etaTriggered = itr->eta();
      phiTriggered = itr->phi() ;
        }
    if (verbosity_%10>0) 
      edm::LogInfo("IsoTrack") << "cenJ pt     " << itr->pt() 
                   << "  eta/phi " << itr->eta() << " " 
                   << itr->phi();
      } 
      edm::Handle<l1extra::L1JetParticleCollection> l1FwdJetHandle;
      iEvent.getByToken(tok_L1extFwdJet_,l1FwdJetHandle);
      for (itr = l1FwdJetHandle->begin(); itr != l1FwdJetHandle->end(); ++itr) {
    if (itr->pt()>ptTriggered) {
      ptTriggered = itr->pt();
      etaTriggered = itr->eta();
      phiTriggered = itr->phi() ;
    }
    if (verbosity_%10>0) 
      edm::LogInfo("IsoTrack") << "forJ pt     " << itr->pt() << " eta/phi "
                   << itr->eta() << " " << itr->phi();
     }
      if (verbosity_%10>0) 
    edm::LogInfo("IsoTrack") << "jets pt/eta/phi = " << ptTriggered << "/" 
                 << etaTriggered << "/" <<phiTriggered;
      unsigned int nTracks(0),nselTracks(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_mindR1 = deltaR(etaTriggered,v4.eta(),phiTriggered,v4.phi());
    t_mindR2 =  -999;
    if (verbosity_%10>0)  
      edm::LogInfo("IsoTrack") << "This track : " << nTracks << " (pt/eta/phi/p) :"
                   << pTrack->pt() << "/" << pTrack->eta() << "/"
                   << pTrack->phi() << "/" << pTrack->p();

    if (verbosity_%10>0) edm::LogInfo("IsoTrack") << "dr values are = " << t_mindR1;  
      
    t_l1pt  = ptTriggered;
    t_l1eta = etaTriggered;
    t_l1phi = phiTriggered;
    t_l3pt  = -999;
    t_l3eta = -999;
    t_l3phi = -999;
 
    //Selection of good track
    t_selectTk = spr::goodTrack(pTrack,leadPV,selectionParameters_,((verbosity_/100)%10>2));
    spr::trackSelectionParameters oneCutParameters = selectionParameters_;
    oneCutParameters.maxDxyPV  = 10;
    oneCutParameters.maxDzPV   = 100;
    oneCutParameters.maxInMiss = 2;
    oneCutParameters.maxOutMiss= 2;
    bool qltyFlag  = spr::goodTrack(pTrack,leadPV,oneCutParameters,((verbosity_/100)%10>2));
    oneCutParameters           = selectionParameters_;
    oneCutParameters.maxDxyPV  = 10;
    oneCutParameters.maxDzPV   = 100;
    t_qltyMissFlag = spr::goodTrack(pTrack,leadPV,oneCutParameters,((verbosity_/100)%10>2));
    oneCutParameters           = selectionParameters_;
    oneCutParameters.maxInMiss = 2;
    oneCutParameters.maxOutMiss= 2;
    t_qltyPVFlag   = spr::goodTrack(pTrack,leadPV,oneCutParameters,((verbosity_/100)%10>2));
    t_ieta = 0;
    if (trkDetItr->okHCAL) {
      HcalDetId detId = (HcalDetId)(trkDetItr->detIdHCAL);
      t_ieta = detId.ieta();
    }
    if (verbosity_%10 > 0) 
      edm::LogInfo("IsoTrack") << "qltyFlag|okECAL|okHCAL : " << qltyFlag << "|" 
                   << trkDetItr->okECAL << "/" << trkDetItr->okHCAL;
    t_qltyFlag = (qltyFlag && trkDetItr->okECAL && trkDetItr->okHCAL);
    t_p        = pTrack->p();
    h_tketa1[0]->Fill(t_ieta);
    for (unsigned int k=1; k<pbin.size(); ++k) {
      if (t_p >= pbin[k-1] && t_p < pbin[k]) {
        h_tketa1[k]->Fill(t_ieta);
        break;
      }
    }
    if (t_qltyFlag) {
      nselTracks++;
      h_tketa2[0]->Fill(t_ieta);
      for (unsigned int k=1; k<pbin.size(); ++k) {
        if (t_p >= pbin[k-1] && t_p < pbin[k]) {
          h_tketa2[k]->Fill(t_ieta);
          break;
        }
      }
      int nRH_eMipDR(0), nNearTRKs(0), nRecHits(-999);
      t_eMipDR = spr::eCone_ecal(geo, barrelRecHitsHandle, 
                     endcapRecHitsHandle,
                     trkDetItr->pointHCAL,
                     trkDetItr->pointECAL,
                     a_mipR_, trkDetItr->directionECAL, 
                     nRH_eMipDR);
      t_DetIds->clear(); t_HitEnergies->clear();
      std::vector<DetId> ids;
      t_eHcal = spr::eCone_hcal(geo, hbhe, trkDetItr->pointHCAL, 
                    trkDetItr->pointECAL, a_coneR_, 
                    trkDetItr->directionHCAL,nRecHits, 
                    ids, *t_HitEnergies);
      for (unsigned int k=0; k<ids.size(); ++k) {
        t_DetIds->push_back(ids[k].rawId());
      }
      t_hmaxNearP = spr::chargeIsolationCone(nTracks,trkCaloDirections,
                         a_charIsoR_, nNearTRKs, 
                         ((verbosity_/100)%10>2));
      if (t_hmaxNearP < 2) {
        h_tketa3[0]->Fill(t_ieta);
        for (unsigned int k=1; k<pbin.size(); ++k) {
          if (t_p >= pbin[k-1] && t_p < pbin[k]) {
        h_tketa3[k]->Fill(t_ieta);
        break;
          }
        }
        if (t_eMipDR < 1) {
          h_tketa4[0]->Fill(t_ieta);
          for (unsigned int k=1; k<pbin.size(); ++k) {
        if (t_p >= pbin[k-1] && t_p < pbin[k]) {
          h_tketa4[k]->Fill(t_ieta);
          break;
        }
          }
          if (t_mindR1 > 1) {
        h_tketa5[0]->Fill(t_ieta);
        for (unsigned int k=1; k<pbin.size(); ++k) {
          if (t_p >= pbin[k-1] && t_p < pbin[k]) {
            h_tketa5[k]->Fill(t_ieta);
            break;
          }
        }
          }
        }
      }
      if (verbosity_%10 > 0) {
        edm::LogInfo("IsoTrack") << "This track : " << nTracks << " (pt/eta/phi/p) :" 
                     << pTrack->pt() << "/" << pTrack->eta() << "/" 
                     << pTrack->phi() << "/" << t_p;
        edm::LogInfo("IsoTrack") << "e_MIP " << t_eMipDR << " Chg Isolation "
                     << t_hmaxNearP << " eHcal" << t_eHcal << " ieta " 
                     << t_ieta << " Quality " << t_qltyMissFlag
                     << ":" << t_qltyPVFlag << ":" << t_selectTk;
        for (unsigned int lll=0;lll<t_DetIds->size();lll++) {
          edm::LogInfo("IsoTrack") << "det id is = " << t_DetIds->at(lll)  << "  " 
                       << " hit enery is  = "  << t_HitEnergies->at(lll);
        }
      }
      if (t_p>20.0 && t_eMipDR<2.0 && t_hmaxNearP<10.0) {
        tree->Fill();
      }
    }
      }
    }
  } 
}

void IsoTrackCalib::beginJob ( void  )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 604 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, gen::k, TFileService::make(), dataset::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, and fftjetcommon_cfi::title.

                             {
  h_RecHit_iEta = fs_->make<TProfile>("rechit_ieta","Rec hit vs. ieta",60,-30,30,0,1000);
  h_RecHit_num = fs_->make<TProfile>("rechit_num","Rec hit vs. num",100,0,20,0,1000);
  h_iEta       = fs_->make<TH1I>("iEta","iEta",60,-30,30);
  h_Rechit_E   = fs_->make<TH1F>("Rechit_E","Rechit_E",100,0,1000);

  double prange[5] = {20,30,40,60,100};
  for (int k=0; k<5; ++k) pbin.push_back(prange[k]);
  std::string type[6] = {"All", "Trigger OK", "Tree Selected", 
             "Charge Isolation", "MIP Cut", "L1 Cut"};
  for (unsigned int k=0; k<pbin.size(); ++k) {
    char name[20], namp[20], title[100];
    if (k == 0) sprintf (namp, "all momentum");
    else        sprintf (namp, "p = %4.0f:%4.0f GeV", pbin[k-1], pbin[k]);
    sprintf (name, "TrackEta0%d", k);
    sprintf (title, "Track #eta for tracks with %s (%s)",namp,type[0].c_str());
    h_tketa0[k] = fs_->make<TH1I>(name, title, 60, -30, 30);
    sprintf (name, "TrackEta1%d", k);
    sprintf (title, "Track #eta for tracks with %s (%s)",namp,type[1].c_str());
    h_tketa1[k] = fs_->make<TH1I>(name, title, 60, -30, 30);
    sprintf (name, "TrackEta2%d", k);
    sprintf (title, "Track #eta for tracks with %s (%s)",namp,type[2].c_str());
    h_tketa2[k] = fs_->make<TH1I>(name, title, 60, -30, 30);
    sprintf (name, "TrackEta3%d", k);
    sprintf (title, "Track #eta for tracks with %s (%s)",namp,type[3].c_str());
    h_tketa3[k] = fs_->make<TH1I>(name, title, 60, -30, 30);
    sprintf (name, "TrackEta4%d", k);
    sprintf (title, "Track #eta for tracks with %s (%s)",namp,type[4].c_str());
    h_tketa4[k] = fs_->make<TH1I>(name, title, 60, -30, 30);
    sprintf (name, "TrackEta5%d", k);
    sprintf (title, "Track #eta for tracks with %s (%s)",namp,type[5].c_str());
    h_tketa5[k] = fs_->make<TH1I>(name, title, 60, -30, 30);
  }
  h_jetpt[0] = fs_->make<TH1F>("Jetpt0","Jet p_T (All)", 500,0.,2500.);
  h_jetpt[1] = fs_->make<TH1F>("Jetpt1","Jet p_T (All Weighted)", 500,0.,2500.);
  h_jetpt[2] = fs_->make<TH1F>("Jetpt2","Jet p_T (|#eta| < 2.5)", 500,0.,2500.);
  h_jetpt[3] = fs_->make<TH1F>("Jetpt3","Jet p_T (|#eta| < 2.5 Weighted)", 500,0.,2500.);

  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_ieta",        &t_ieta,        "t_ieta/I");
  tree->Branch("t_EventWeight", &t_EventWeight, "t_EventWeight/D");
  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_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_eHcal",       &t_eHcal,       "t_eHcal/D");
  tree->Branch("t_hmaxNearP",   &t_hmaxNearP,   "t_hmaxNearP/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)");

  t_DetIds      = new std::vector<unsigned int>();
  t_HitEnergies = new std::vector<double>();
  t_l1bits      = 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_l1bits","std::vector<bool>", &t_l1bits);
}

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

overrideprivate

Definition at line 674 of file IsoTrackCalib.cc.

References plotBeamSpotDB::first, L1GtTriggerMenu::gtAlgorithmMap(), L1GtTriggerMenu::gtTriggerMenuName(), hltConfig_, HLTConfigProvider::init(), l1AlgoMap_, m_l1GtMenu, m_l1GtUtils, convertSQLiteXML::ok, L1GtUtils::ptrL1TriggerMenuEventSetup(), edm::RunBase::run(), AlCaHLTBitMon_QueryRunRegistry::string, and verbosity_.

Referenced by endJob().

                                                                            {
  bool changed = false;
  bool ok = hltConfig_.init(iRun,iSetup,"HLT",changed);
  edm::LogInfo("IsoTrack")  << "Run " << iRun.run() << " hltconfig.init " << ok;

  int iErrorCode = -1;
  m_l1GtMenu   = m_l1GtUtils.ptrL1TriggerMenuEventSetup(iErrorCode);
  const AlgorithmMap& algorithmMap = m_l1GtMenu->gtAlgorithmMap();
  const std::string&  menuName     = m_l1GtMenu->gtTriggerMenuName();

  if (verbosity_%10>0) edm::LogInfo("IsoTrack") << "menuName " << menuName;
  for (CItAlgo itAlgo = algorithmMap.begin(); itAlgo != algorithmMap.end(); 
       itAlgo++) {
    std::string algName = itAlgo->first;
    int algBitNumber = ( itAlgo->second ).algoBitNumber();
    l1AlgoMap_.insert( std::pair<std::pair<unsigned int,std::string>,int>( std::pair<unsigned int,std::string>(algBitNumber, algName) , 0) ) ;
  }
  std::map< std::pair<unsigned int,std::string>, int>::iterator itr;
  for (itr=l1AlgoMap_.begin(); itr!=l1AlgoMap_.end(); itr++) {
    if (verbosity_%10>0)  
      edm::LogInfo("IsoTrack") << " ********** " << (itr->first).first << " "
                   << (itr->first).second <<" "<<itr->second;
  }
}

double IsoTrackCalib::deltaR ( double  eta1, double  eta2, double  phi1, double  phi2  )

private

Definition at line 726 of file IsoTrackCalib.cc.

References DEFINE_FWK_MODULE, M_PI, and mathSSE::sqrt().

Referenced by analyze(), and endJob().

                                                                           {
  double deta = eta1-eta2;
  if (phi1 < 0) phi1 += 2.0*M_PI;
  if (phi2 < 0) phi2 += 2.0*M_PI;
  double dphi = phi1-phi2;
  if (dphi > M_PI)       dphi -= 2.*M_PI;
  else if (dphi < -M_PI) dphi += 2.*M_PI;
  return std::sqrt(deta*deta + dphi*dphi); 
}

double IsoTrackCalib::dEta ( math::XYZTLorentzVector &  vec1, math::XYZTLorentzVector &  vec2  )

private

Definition at line 704 of file IsoTrackCalib.cc.

Referenced by dR(), and endJob().

                                                                                   {
  return (vec1.eta()-vec2.eta());
}

double IsoTrackCalib::dPhi ( math::XYZTLorentzVector &  vec1, math::XYZTLorentzVector &  vec2  )

private

Definition at line 708 of file IsoTrackCalib.cc.

References M_PI.

Referenced by dR(), and endJob().

                                                                                   {

  double phi1 = vec1.phi();
  if (phi1 < 0) phi1 += 2.0*M_PI;
  double phi2 = vec2.phi();
  if (phi2 < 0) phi2 += 2.0*M_PI;
  double dphi = phi1-phi2;
  if (dphi > M_PI)       dphi -= 2.*M_PI;
  else if (dphi < -M_PI) dphi += 2.*M_PI;
  return dphi;
}

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

private

Definition at line 720 of file IsoTrackCalib.cc.

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

Referenced by endJob().

                                                                                 {
  double deta = dEta(vec1,vec2);
  double dphi = dPhi(vec1,vec2);
  return std::sqrt(deta*deta + dphi*dphi);
}

void IsoTrackCalib::endJob ( void  )

inlineoverrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 93 of file IsoTrackCalib.cc.

References beginRun(), deltaR(), dEta(), dPhi(), dR(), and endRun().

Referenced by o2olib.O2ORunMgr::executeJob().

{}

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

overrideprivate

Definition at line 700 of file IsoTrackCalib.cc.

References edm::RunBase::run().

Referenced by endJob().

                                                                   {
  edm::LogInfo("IsoTrack") << "endRun " << iRun.run() << std::endl;
}

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

static

Definition at line 232 of file IsoTrackCalib.cc.

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

                                                                               {

  std::vector<std::string> seeds = {"L1_SingleJet36","L1_SingleJet52",
                    "L1_SingleJet68","L1_SingleJet92",
                    "L1_SingleJet128"};
  edm::ParameterSetDescription desc;
  desc.addUntracked<int>("Verbosity",       0);
  desc.addUntracked<std::vector<std::string> >("L1Seed", seeds);
  desc.addUntracked<std::string>("TrackQuality", "highPurity");
  desc.addUntracked<double>("MinTrackPt",   10.0);
  desc.addUntracked<double>("MaxDxyPV",     0.02);
  desc.addUntracked<double>("MaxDzPV",      0.02);
  desc.addUntracked<double>("MaxChi2",      5.0);
  desc.addUntracked<double>("MaxDpOverP",   0.1);
  desc.addUntracked<int>("MinOuterHit",     4);
  desc.addUntracked<int>("MinLayerCrossed", 8);
  desc.addUntracked<int>("MaxInMiss",       0);
  desc.addUntracked<int>("MaxOutMiss",      0);
  desc.addUntracked<double>("ConeRadius",   34.98);
  desc.addUntracked<double>("ConeRadiusMIP",14.0);
  desc.addUntracked<bool>("IsItAOD",        false);
  descriptions.add("isoTrackCalib",desc);
}

Member Data Documentation

double IsoTrackCalib::a_charIsoR_

private

Definition at line 109 of file IsoTrackCalib.cc.

Referenced by analyze(), and IsoTrackCalib().

double IsoTrackCalib::a_coneR_

private

Definition at line 109 of file IsoTrackCalib.cc.

Referenced by analyze(), and IsoTrackCalib().

double IsoTrackCalib::a_mipR_

private

Definition at line 109 of file IsoTrackCalib.cc.

Referenced by analyze(), and IsoTrackCalib().

edm::Service<TFileService> IsoTrackCalib::fs_

private

Definition at line 102 of file IsoTrackCalib.cc.

Referenced by beginJob().

TH1I* IsoTrackCalib::h_iEta

private

Definition at line 139 of file IsoTrackCalib.cc.

Referenced by analyze(), and beginJob().

TH1F * IsoTrackCalib::h_jetp

private

Definition at line 141 of file IsoTrackCalib.cc.

TH1F* IsoTrackCalib::h_jetpt[4]

private

Definition at line 142 of file IsoTrackCalib.cc.

Referenced by analyze(), and beginJob().

TH1F* IsoTrackCalib::h_Rechit_E

private

Definition at line 141 of file IsoTrackCalib.cc.

Referenced by analyze(), and beginJob().

TProfile* IsoTrackCalib::h_RecHit_iEta

private

Definition at line 138 of file IsoTrackCalib.cc.

Referenced by analyze(), and beginJob().

TProfile * IsoTrackCalib::h_RecHit_num

private

Definition at line 138 of file IsoTrackCalib.cc.

Referenced by analyze(), and beginJob().

TH1I * IsoTrackCalib::h_tkEta0[5]

private

Definition at line 139 of file IsoTrackCalib.cc.

TH1I* IsoTrackCalib::h_tketa0[6]

private

Definition at line 143 of file IsoTrackCalib.cc.

Referenced by analyze(), and beginJob().

TH1I * IsoTrackCalib::h_tkEta1[5]

private

Definition at line 139 of file IsoTrackCalib.cc.

TH1I * IsoTrackCalib::h_tketa1[6]

private

Definition at line 143 of file IsoTrackCalib.cc.

Referenced by analyze(), and beginJob().

TH1I * IsoTrackCalib::h_tkEta2[5]

private

Definition at line 139 of file IsoTrackCalib.cc.

TH1I * IsoTrackCalib::h_tketa2[6]

private

Definition at line 143 of file IsoTrackCalib.cc.

Referenced by analyze(), and beginJob().

TH1I* IsoTrackCalib::h_tkEta3[5]

private

Definition at line 140 of file IsoTrackCalib.cc.

TH1I* IsoTrackCalib::h_tketa3[6]

private

Definition at line 144 of file IsoTrackCalib.cc.

Referenced by analyze(), and beginJob().

TH1I * IsoTrackCalib::h_tkEta4[5]

private

Definition at line 140 of file IsoTrackCalib.cc.

TH1I * IsoTrackCalib::h_tketa4[6]

private

Definition at line 144 of file IsoTrackCalib.cc.

Referenced by analyze(), and beginJob().

TH1I * IsoTrackCalib::h_tkEta5[5]

private

Definition at line 140 of file IsoTrackCalib.cc.

TH1I * IsoTrackCalib::h_tketa5[6]

private

Definition at line 144 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::map< std::pair<unsigned int,std::string>, int> IsoTrackCalib::l1AlgoMap_

private

Definition at line 145 of file IsoTrackCalib.cc.

Referenced by beginRun().

std::vector<std::string> IsoTrackCalib::l1Names_

private

Definition at line 105 of file IsoTrackCalib.cc.

Referenced by analyze(), and IsoTrackCalib().

const L1GtTriggerMenu* IsoTrackCalib::m_l1GtMenu

private

Definition at line 110 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().

std::vector<double> IsoTrackCalib::pbin

private

Definition at line 137 of file IsoTrackCalib.cc.

Referenced by analyze(), and beginJob().

spr::trackSelectionParameters IsoTrackCalib::selectionParameters_

private

Definition at line 107 of file IsoTrackCalib.cc.

Referenced by analyze(), and IsoTrackCalib().

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

private

Definition at line 136 of file IsoTrackCalib.cc.

Referenced by analyze(), and beginJob().

double IsoTrackCalib::t_eHcal

private

Definition at line 133 of file IsoTrackCalib.cc.

Referenced by analyze(), and beginJob().

double IsoTrackCalib::t_eMipDR

private

Definition at line 133 of file IsoTrackCalib.cc.

Referenced by analyze(), and beginJob().

int IsoTrackCalib::t_Event

private

Definition at line 130 of file IsoTrackCalib.cc.

Referenced by analyze(), and beginJob().

double IsoTrackCalib::t_EventWeight

private

Definition at line 131 of file IsoTrackCalib.cc.

Referenced by analyze(), and beginJob().

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

private

Definition at line 137 of file IsoTrackCalib.cc.

Referenced by analyze(), and beginJob().

double IsoTrackCalib::t_hmaxNearP

private

Definition at line 133 of file IsoTrackCalib.cc.

Referenced by analyze(), and beginJob().

int IsoTrackCalib::t_ieta

private

Definition at line 130 of file IsoTrackCalib.cc.

Referenced by analyze(), and beginJob().

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

private

Definition at line 111 of file IsoTrackCalib.cc.

Referenced by analyze(), and beginJob().

double IsoTrackCalib::t_l1eta

private

Definition at line 131 of file IsoTrackCalib.cc.

Referenced by analyze(), and beginJob().

double IsoTrackCalib::t_l1phi

private

Definition at line 131 of file IsoTrackCalib.cc.

Referenced by analyze(), and beginJob().

double IsoTrackCalib::t_l1pt

private

Definition at line 131 of file IsoTrackCalib.cc.

Referenced by analyze(), and beginJob().

double IsoTrackCalib::t_l3eta

private

Definition at line 132 of file IsoTrackCalib.cc.

Referenced by analyze(), and beginJob().

double IsoTrackCalib::t_l3phi

private

Definition at line 132 of file IsoTrackCalib.cc.

Referenced by analyze(), and beginJob().

double IsoTrackCalib::t_l3pt

private

Definition at line 132 of file IsoTrackCalib.cc.

Referenced by analyze(), and beginJob().

double IsoTrackCalib::t_mindR1

private

Definition at line 132 of file IsoTrackCalib.cc.

Referenced by analyze(), and beginJob().

double IsoTrackCalib::t_mindR2

private

Definition at line 133 of file IsoTrackCalib.cc.

Referenced by analyze(), and beginJob().

double IsoTrackCalib::t_p

private

Definition at line 132 of file IsoTrackCalib.cc.

Referenced by analyze(), and beginJob().

bool IsoTrackCalib::t_qltyFlag

private

Definition at line 135 of file IsoTrackCalib.cc.

Referenced by analyze(), and beginJob().

bool IsoTrackCalib::t_qltyMissFlag

private

Definition at line 135 of file IsoTrackCalib.cc.

Referenced by analyze(), and beginJob().

bool IsoTrackCalib::t_qltyPVFlag

private

Definition at line 135 of file IsoTrackCalib.cc.

Referenced by analyze(), and beginJob().

int IsoTrackCalib::t_Run

private

Definition at line 130 of file IsoTrackCalib.cc.

Referenced by analyze(), and beginJob().

bool IsoTrackCalib::t_selectTk

private

Definition at line 135 of file IsoTrackCalib.cc.

Referenced by analyze(), and beginJob().

std::string IsoTrackCalib::theTrackQuality_

private

Definition at line 108 of file IsoTrackCalib.cc.

Referenced by analyze(), and IsoTrackCalib().

edm::InputTag IsoTrackCalib::theTriggerResultsLabel_

private

Definition at line 112 of file IsoTrackCalib.cc.

Referenced by IsoTrackCalib().

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

private

Definition at line 118 of file IsoTrackCalib.cc.

Referenced by analyze(), and IsoTrackCalib().

edm::EDGetTokenT<EcalRecHitCollection> IsoTrackCalib::tok_EB_

private

Definition at line 119 of file IsoTrackCalib.cc.

Referenced by analyze(), and IsoTrackCalib().

edm::EDGetTokenT<EcalRecHitCollection> IsoTrackCalib::tok_EE_

private

Definition at line 120 of file IsoTrackCalib.cc.

Referenced by analyze(), and IsoTrackCalib().

edm::EDGetTokenT<GenEventInfoProduct> IsoTrackCalib::tok_ew_

private

Definition at line 122 of file IsoTrackCalib.cc.

Referenced by analyze(), and IsoTrackCalib().

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

private

Definition at line 116 of file IsoTrackCalib.cc.

Referenced by analyze(), and IsoTrackCalib().

edm::EDGetTokenT<HBHERecHitCollection> IsoTrackCalib::tok_hbhe_

private

Definition at line 121 of file IsoTrackCalib.cc.

Referenced by analyze(), and IsoTrackCalib().

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

private

Definition at line 123 of file IsoTrackCalib.cc.

Referenced by analyze(), and IsoTrackCalib().

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

private

Definition at line 126 of file IsoTrackCalib.cc.

Referenced by analyze(), and IsoTrackCalib().

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

private

Definition at line 127 of file IsoTrackCalib.cc.

Referenced by analyze(), and IsoTrackCalib().

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

private

Definition at line 125 of file IsoTrackCalib.cc.

Referenced by analyze(), and IsoTrackCalib().

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

private

Definition at line 124 of file IsoTrackCalib.cc.

Referenced by analyze(), and IsoTrackCalib().

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

private

Definition at line 117 of file IsoTrackCalib.cc.

Referenced by analyze(), and IsoTrackCalib().

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

private

Definition at line 113 of file IsoTrackCalib.cc.

Referenced by IsoTrackCalib().

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

private

Definition at line 114 of file IsoTrackCalib.cc.

Referenced by IsoTrackCalib().

TTree* IsoTrackCalib::tree

private

Definition at line 129 of file IsoTrackCalib.cc.

Referenced by eventstfile.Events::__iter__(), svgfig.SVG::__str__(), 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().

edm::InputTag IsoTrackCalib::triggerEvent_

private

Definition at line 112 of file IsoTrackCalib.cc.

Referenced by IsoTrackCalib().

int IsoTrackCalib::verbosity_

private

Definition at line 106 of file IsoTrackCalib.cc.

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