IsoTrackCalibration Class Reference

Inheritance diagram for IsoTrackCalibration:

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

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

Private Member Functions
void	analyze (edm::Event const &, edm::EventSetup const &) override
void	beginJob () override
virtual void	beginLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &)
void	beginRun (edm::Run const &, edm::EventSetup const &) override
virtual void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &)
void	endRun (edm::Run const &, edm::EventSetup const &) override

Private Attributes
double	a_charIsoR_
double	a_coneR_
double	a_mipR_
double	constTrackPt_
edm::Service< TFileService >	fs_
TH1I *	h_iEta
TH1F *	h_jetp
TH1F *	h_jetpt [4]
TH1F *	h_nTrk
TH1F *	h_nVtx
TH1F *	h_Rechit_E
TProfile *	h_RecHit_iEta
TProfile *	h_RecHit_num
TH1I *	h_tketa0 [5]
TH1I *	h_tketa1 [5]
TH1I *	h_tketa2 [5]
TH1I *	h_tketa3 [5]
TH1I *	h_tketa4 [5]
TH1I *	h_tketa5 [5]
TH1I *	h_tketav1 [5][6]
TH1I *	h_tketav2 [5][6]
HLTConfigProvider	hltConfig_
std::vector< std::string >	HLTNames_
bool	isAOD_
bool	isMC_
bool	isQCD_
double	maxEcalEnr_
double	maxNeighborTrackEnr_
int	nRun_
spr::trackSelectionParameters	selectionParameters_
double	slopeTrackPt_
std::vector< unsigned int > *	t_DetIds
std::vector< unsigned int > *	t_DetIds1
std::vector< unsigned int > *	t_DetIds3
double	t_eHcal
double	t_eHcal10
double	t_eHcal30
double	t_eMipDR
int	t_Event
double	t_EventWeight
std::vector< double > *	t_HitEnergies
std::vector< double > *	t_HitEnergies1
std::vector< double > *	t_HitEnergies3
double	t_hmaxNearP
int	t_ieta
int	t_nTrk
int	t_nVtx
double	t_p
double	t_phi
double	t_pt
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::VertexCollection >	tok_recVtx_
edm::EDGetTokenT< trigger::TriggerEvent >	tok_trigEvt_
edm::EDGetTokenT< edm::TriggerResults >	tok_trigRes_
TTree *	tree
edm::InputTag	triggerEvent_
std::vector< std::string >	trigNames_
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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Definition at line 74 of file IsoTrackCalibration.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 131 of file IsoTrackCalibration.cc.

References a_charIsoR_, a_coneR_, a_mipR_, constTrackPt_, gather_cfg::cout, edm::ParameterSet::getUntrackedParameter(), HLT_2018_cff::InputTag, isAOD_, isMC_, isQCD_, dqmdumpme::k, spr::trackSelectionParameters::maxChi2, spr::trackSelectionParameters::maxDpOverP, spr::trackSelectionParameters::maxDxyPV, spr::trackSelectionParameters::maxDzPV, maxEcalEnr_, spr::trackSelectionParameters::maxInMiss, maxNeighborTrackEnr_, spr::trackSelectionParameters::maxOutMiss, spr::trackSelectionParameters::minLayerCrossed, spr::trackSelectionParameters::minOuterHit, spr::trackSelectionParameters::minPt, spr::trackSelectionParameters::minQuality, reco::TrackBase::qualityByName(), selectionParameters_, slopeTrackPt_, AlCaHLTBitMon_QueryRunRegistry::string, theTrackQuality_, theTriggerResultsLabel_, tok_bs_, tok_EB_, tok_EE_, tok_ew_, tok_genTrack_, tok_hbhe_, tok_recVtx_, tok_trigEvt_, tok_trigRes_, triggerEvent_, trigNames_, and verbosity_.

                                                                       : nRun_(0) {
  usesResource("TFileService");
  //now do whatever initialization is needed
  const double isolationRadius(28.9);
  verbosity_ = iConfig.getUntrackedParameter<int>("Verbosity", 0);
  trigNames_ = iConfig.getUntrackedParameter<std::vector<std::string> >("Triggers");
  theTrackQuality_ = iConfig.getUntrackedParameter<std::string>("TrackQuality", "highPurity");
  reco::TrackBase::TrackQuality trackQuality_ = reco::TrackBase::qualityByName(theTrackQuality_);
  constTrackPt_ = iConfig.getUntrackedParameter<double>("ConstTrackPt", 10.0);
  slopeTrackPt_ = iConfig.getUntrackedParameter<double>("SlopeTrackPt", 0.16);
  selectionParameters_.minPt = constTrackPt_;
  selectionParameters_.minQuality = trackQuality_;
  selectionParameters_.maxDxyPV = iConfig.getUntrackedParameter<double>("MaxDxyPV", 0.02);
  selectionParameters_.maxDzPV = iConfig.getUntrackedParameter<double>("MaxDzPV", 0.02);
  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_ + isolationRadius;
  a_mipR_ = iConfig.getUntrackedParameter<double>("ConeRadiusMIP", 14.0);
  maxEcalEnr_ = iConfig.getUntrackedParameter<double>("MaxEcalEnergyInCone", 2.5);
  maxNeighborTrackEnr_ = iConfig.getUntrackedParameter<double>("MaxNeighborTrackEnergy", 40.0);
  isMC_ = iConfig.getUntrackedParameter<bool>("IsMC", false);
  isQCD_ = iConfig.getUntrackedParameter<bool>("IsQCD", false);
  isAOD_ = iConfig.getUntrackedParameter<bool>("IsAOD", true);
  triggerEvent_ = edm::InputTag("hltTriggerSummaryAOD", "", "HLT");
  theTriggerResultsLabel_ = edm::InputTag("TriggerResults", "", "HLT");

  // 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("generator"));

  if (isAOD_) {
    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"));
  }

#ifdef DebugLog
  if (verbosity_ >= 0) {
    std::cout << "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 isMC " << isMC_ << "\t isQCD " << isQCD_ << "\t isAOD " << isAOD_
              << std::endl;
    std::cout << trigNames_.size() << " triggers to be studied:";
    for (unsigned int k = 0; k < trigNames_.size(); ++k)
      std::cout << " " << trigNames_[k];
    std::cout << std::endl;
  }
#endif
}

IsoTrackCalibration::~IsoTrackCalibration ( )

override

Definition at line 200 of file IsoTrackCalibration.cc.

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

Member Function Documentation

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

overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 205 of file IsoTrackCalibration.cc.

References a_charIsoR_, a_coneR_, a_mipR_, funct::abs(), edm::HLTGlobalStatus::accept(), Calorimetry_cff::bField, edm::EventBase::bunchCrossing(), spr::chargeIsolationCone(), constTrackPt_, gather_cfg::cout, spr::eCone_ecal(), spr::eCone_hcal(), reco::TrackBase::eta(), edm::EventID::event(), spr::find(), genParticles2HepMC_cfi::genEventInfo, edm::EventSetup::get(), edm::Event::getByToken(), spr::goodTrack(), h_nTrk, h_nVtx, photonIsolationHIProducer_cfi::hbhe, ValidationMatrix::hlt, mps_fire::i, edm::EventBase::id(), photons_cff::ids, HcalDetId::ieta(), isMC_, isQCD_, edm::HandleBase::isValid(), dqmdumpme::k, edm::InputTag::label(), edm::EventBase::luminosityBlock(), spr::trackSelectionParameters::maxDxyPV, spr::trackSelectionParameters::maxDzPV, maxEcalEnr_, spr::trackSelectionParameters::maxInMiss, maxNeighborTrackEnr_, spr::trackSelectionParameters::maxOutMiss, spr::trackSelectionParameters::minPt, nTracks(), reco::TrackBase::p(), reco::TrackBase::phi(), reco::BeamSpot::position(), edm::ESHandle< T >::product(), spr::propagateCALO(), reco::TrackBase::pt(), edm::EventID::run(), selectionParameters_, edm::HLTGlobalStatus::size(), slopeTrackPt_, t_DetIds, t_DetIds1, t_DetIds3, t_eHcal, t_eHcal10, t_eHcal30, t_eMipDR, t_Event, t_EventWeight, t_HitEnergies, t_HitEnergies1, t_HitEnergies3, t_hmaxNearP, t_ieta, t_nTrk, t_nVtx, t_p, t_phi, t_pt, t_qltyMissFlag, t_qltyPVFlag, t_Run, t_selectTk, theTrackQuality_, tok_bs_, tok_EB_, tok_EE_, tok_ew_, tok_genTrack_, tok_hbhe_, tok_recVtx_, tok_trigEvt_, tok_trigRes_, PDWG_DiPhoton_SD_cff::triggerEvent, triggerEvent_, L1TEGammaOffline_cfi::triggerNames, edm::TriggerNames::triggerNames(), edm::Event::triggerNames(), triggerResults, trigNames_, verbosity_, GenEventInfoProduct::weight(), x, y, and z.

                                                                                     {
  t_Run = iEvent.id().run();
  t_Event = iEvent.id().event();
#ifdef DebugLog
  if (verbosity_ % 10 > 0)
    std::cout << "Run " << t_Run << " Event " << t_Event << " Luminosity " << iEvent.luminosityBlock() << " Bunch "
              << iEvent.bunchCrossing() << " starts ==========" << std::endl;
#endif
  //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);

  //event weight for FLAT sample and PU information
  t_EventWeight = 1.0;
  edm::Handle<GenEventInfoProduct> genEventInfo;
  iEvent.getByToken(tok_ew_, genEventInfo);
  if (genEventInfo.isValid())
    t_EventWeight = genEventInfo->weight();

  //Define the best vertex and the beamspot
  edm::Handle<reco::VertexCollection> recVtxs;
  iEvent.getByToken(tok_recVtx_, recVtxs);
  edm::Handle<reco::BeamSpot> beamSpotH;
  iEvent.getByToken(tok_bs_, beamSpotH);
  math::XYZPoint leadPV(0, 0, 0);

  t_nVtx = recVtxs->size();
  h_nVtx->Fill(t_nVtx);

  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();
  }
#ifdef DebugLog
  if ((verbosity_ / 100) % 10 > 0) {
    std::cout << "Primary Vertex " << leadPV;
    if (beamSpotH.isValid())
      std::cout << " Beam Spot " << beamSpotH->position();
    std::cout << std::endl;
  }
#endif

  // 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);

  //Trigger
  bool triggerOK = false;
  if (isMC_ && !isQCD_) {
    triggerOK = true;  // ignore HLT for single pion MC
  } else {
    trigger::TriggerEvent triggerEvent;
    edm::Handle<trigger::TriggerEvent> triggerEventHandle;
    iEvent.getByToken(tok_trigEvt_, triggerEventHandle);
    if (!triggerEventHandle.isValid()) {
#ifdef DebugLog
      std::cout << "Error! Can't get the product " << triggerEvent_.label() << std::endl;
#endif
    } else {
      edm::Handle<edm::TriggerResults> triggerResults;
      iEvent.getByToken(tok_trigRes_, triggerResults);
      if (triggerResults.isValid()) {
        const edm::TriggerNames &triggerNames = iEvent.triggerNames(*triggerResults);
        const std::vector<std::string> &triggerNames_ = triggerNames.triggerNames();
        for (unsigned int iHLT = 0; iHLT < triggerResults->size(); iHLT++) {
          int hlt = triggerResults->accept(iHLT);
          if (hlt > 0) {
            for (unsigned int i = 0; i < trigNames_.size(); ++i) {
              if (triggerNames_[iHLT].find(trigNames_[i]) != std::string::npos) {
                triggerOK = true;
#ifdef DebugLog
                if (verbosity_ % 10 > 0)
                  std::cout << "This is the trigger we are looking for " << triggerNames_[iHLT] << " Flag " << hlt
                            << std::endl;
#endif
              }
            }
          }
        }
      }
    }
  }

  if (triggerOK) {
    //Propagate tracks to calorimeter surface)
    std::vector<spr::propagatedTrackDirection> trkCaloDirections;
    spr::propagateCALO(trkCollection, geo, bField, theTrackQuality_, trkCaloDirections, ((verbosity_ / 100) % 10 > 2));
    //Loop over tracks
    std::vector<spr::propagatedTrackDirection>::const_iterator trkDetItr;
    unsigned int nTracks(0), nselTracks(0);
    t_nTrk = trkCaloDirections.size();

    for (trkDetItr = trkCaloDirections.begin(), nTracks = 0; trkDetItr != trkCaloDirections.end();
         trkDetItr++, nTracks++) {
      const reco::Track *pTrack = &(*(trkDetItr->trkItr));
#ifdef DebugLog
      if (verbosity_ % 10 > 0)
        std::cout << "This track : " << nTracks << " (pt/eta/phi/p) :" << pTrack->pt() << "/" << pTrack->eta() << "/"
                  << pTrack->phi() << "/" << pTrack->p() << std::endl;
#endif

      t_ieta = 0;
      if (trkDetItr->okHCAL) {
        HcalDetId detId = (HcalDetId)(trkDetItr->detIdHCAL);
        t_ieta = detId.ieta();
      }
      // ---------- eta-dependent restriction on Pt ----------------------------
      selectionParameters_.minPt = constTrackPt_ - std::abs((double)t_ieta) * slopeTrackPt_;
      // -----------------------------------------------------------------------

      //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));

#ifdef DebugLog
      if (verbosity_ % 10 > 0)
        std::cout << "qltyFlag|okECAL|okHCAL : " << qltyFlag << "|" << trkDetItr->okECAL << "/" << trkDetItr->okHCAL
                  << std::endl;
#endif
      if (qltyFlag && trkDetItr->okECAL && trkDetItr->okHCAL) {
        nselTracks++;
        int nRH_eMipDR(0), nNearTRKs(0);
        //------ ecal energy around track -------------------------------
        t_eMipDR = spr::eCone_ecal(geo,
                                   barrelRecHitsHandle,
                                   endcapRecHitsHandle,
                                   trkDetItr->pointHCAL,
                                   trkDetItr->pointECAL,
                                   a_mipR_,
                                   trkDetItr->directionECAL,
                                   nRH_eMipDR);
        //---- isolation criteria ----------------------------------------------
        t_hmaxNearP =
            spr::chargeIsolationCone(nTracks, trkCaloDirections, a_charIsoR_, nNearTRKs, ((verbosity_ / 100) % 10 > 2));

        if (t_eMipDR < maxEcalEnr_ && t_hmaxNearP < maxNeighborTrackEnr_) {
          //------------ HCAL --------------------------------------------------
          //------ initialize arrays of DetID and hit energies -----------------
          t_DetIds->clear();
          t_DetIds1->clear();
          t_DetIds3->clear();
          t_HitEnergies->clear();
          t_HitEnergies1->clear();
          t_HitEnergies3->clear();
          int nRecHits(-999), nRecHits1(-999), nRecHits3(-999);
          std::vector<DetId> ids, ids1, ids3;
          //------ hcal energy in the main cone -------------------------------
          t_eHcal = spr::eCone_hcal(geo,
                                    hbhe,
                                    trkDetItr->pointHCAL,
                                    trkDetItr->pointECAL,
                                    a_coneR_,
                                    trkDetItr->directionHCAL,
                                    nRecHits,
                                    ids,
                                    *t_HitEnergies);
          t_DetIds->reserve(ids.size());
          for (unsigned int k = 0; k < ids.size(); ++k) {
            t_DetIds->push_back(ids[k].rawId());
          }
          //----- hcal energy in the extended cone 1 (a_coneR+10) --------------
          t_eHcal10 = spr::eCone_hcal(geo,
                                      hbhe,
                                      trkDetItr->pointHCAL,
                                      trkDetItr->pointECAL,
                                      a_coneR_ + 10,
                                      trkDetItr->directionHCAL,
                                      nRecHits1,
                                      ids1,
                                      *t_HitEnergies1);
          t_DetIds1->reserve(ids1.size());
          for (unsigned int k = 0; k < ids1.size(); ++k) {
            t_DetIds1->push_back(ids1[k].rawId());
          }
          //----- hcal energy in the extended cone 3 (a_coneR+30) --------------
          t_eHcal30 = spr::eCone_hcal(geo,
                                      hbhe,
                                      trkDetItr->pointHCAL,
                                      trkDetItr->pointECAL,
                                      a_coneR_ + 30,
                                      trkDetItr->directionHCAL,
                                      nRecHits3,
                                      ids3,
                                      *t_HitEnergies3);
          t_DetIds3->reserve(ids3.size());
          for (unsigned int k = 0; k < ids3.size(); ++k) {
            t_DetIds3->push_back(ids3[k].rawId());
          }

          t_p = pTrack->p();
          t_pt = pTrack->pt();
          t_phi = pTrack->phi();
#ifdef DebugLog
          if (verbosity_ % 10 > 0) {
            std::cout << "This track : " << nTracks << " (pt/eta/phi/p) :" << pTrack->pt() << "/" << pTrack->eta()
                      << "/" << pTrack->phi() << "/" << t_p << std::endl;
            std::cout << "e_MIP " << t_eMipDR << " Chg Isolation " << t_hmaxNearP << " eHcal" << t_eHcal << " ieta "
                      << t_ieta << " Quality " << t_qltyMissFlag << ":" << t_qltyPVFlag << ":" << t_selectTk
                      << std::endl;
            for (unsigned int lll = 0; lll < t_DetIds->size(); lll++) {
              std::cout << "det id is = " << t_DetIds->at(lll) << "  "
                        << " hit enery is  = " << t_HitEnergies->at(lll) << std::endl;
            }
          }
#endif
          tree->Fill();
        }  // end of conditions on t_eMipDR and t_hmaxNearP
      }    // end of loose check of track quality
    }      // end of loop over tracks

    h_nTrk->Fill(nTracks);
  }  // end of triggerOK
}

void IsoTrackCalibration::beginJob ( void  )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 451 of file IsoTrackCalibration.cc.

References fs_, h_nTrk, h_nVtx, TFileService::make(), t_DetIds, t_DetIds1, t_DetIds3, t_eHcal, t_eHcal10, t_eHcal30, t_eMipDR, t_Event, t_EventWeight, t_HitEnergies, t_HitEnergies1, t_HitEnergies3, t_hmaxNearP, t_ieta, t_nTrk, t_nVtx, t_p, t_phi, t_pt, t_qltyMissFlag, t_qltyPVFlag, t_Run, and t_selectTk.

                                   {
  h_nVtx = fs_->make<TH1F>("h_nVtx", "h_nVtx", 100, 0, 100);
  h_nTrk = fs_->make<TH1F>("h_nTrk", "h_nTrk", 100, 0, 2000);

  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_nVtx", &t_nVtx, "t_nVtx/I");
  tree->Branch("t_nTrk", &t_nTrk, "t_nTrk/I");
  tree->Branch("t_EventWeight", &t_EventWeight, "t_EventWeight/D");
  tree->Branch("t_p", &t_p, "t_p/D");
  tree->Branch("t_pt", &t_pt, "t_pt/D");
  tree->Branch("t_ieta", &t_ieta, "t_ieta/I");
  tree->Branch("t_phi", &t_phi, "t_phi/D");
  tree->Branch("t_eMipDR", &t_eMipDR, "t_eMipDR/D");
  tree->Branch("t_eHcal", &t_eHcal, "t_eHcal/D");
  tree->Branch("t_eHcal10", &t_eHcal10, "t_eHcal10/D");
  tree->Branch("t_eHcal30", &t_eHcal30, "t_eHcal30/D");
  tree->Branch("t_hmaxNearP", &t_hmaxNearP, "t_hmaxNearP/D");
  tree->Branch("t_selectTk", &t_selectTk, "t_selectTk/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_DetIds1 = new std::vector<unsigned int>();
  t_DetIds3 = new std::vector<unsigned int>();
  t_HitEnergies = new std::vector<double>();
  t_HitEnergies1 = new std::vector<double>();
  t_HitEnergies3 = new std::vector<double>();

  tree->Branch("t_DetIds", "std::vector<unsigned int>", &t_DetIds);
  //tree->Branch("t_DetIds1",      "std::vector<unsigned int>", &t_DetIds1);
  //tree->Branch("t_DetIds3",      "std::vector<unsigned int>", &t_DetIds3);
  tree->Branch("t_HitEnergies", "std::vector<double>", &t_HitEnergies);
  //tree->Branch("t_HitEnergies1", "std::vector<double>",       &t_HitEnergies1);
  //tree->Branch("t_HitEnergies3", "std::vector<double>",       &t_HitEnergies3);
}

void IsoTrackCalibration::beginLuminosityBlock ( edm::LuminosityBlock const &  , edm::EventSetup const &    )

privatevirtual

Definition at line 521 of file IsoTrackCalibration.cc.

{}

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

overrideprivate

Definition at line 491 of file IsoTrackCalibration.cc.

References RemoveAddSevLevel::flag, hltConfig_, HLTConfigProvider::init(), dqmiodumpmetadata::n, nRun_, edm::RunBase::run(), HLTConfigProvider::size(), HLTConfigProvider::triggerIndex(), and trigNames_.

                                                                                  {
  bool changed_(true);
  bool flag = hltConfig_.init(iRun, iSetup, "HLT", changed_);
  edm::LogInfo("HcalIsoTrack") << "Run[" << nRun_ << "] " << iRun.run() << " process HLT init flag " << flag
                               << " change flag " << changed_;

  // check if trigger names in (new) config
  if (changed_) {
    changed_ = false;
    edm::LogInfo("HcalIsoTrack") << "New trigger menu found !!!";
    const unsigned int n(hltConfig_.size());
    for (unsigned itrig = 0; itrig < trigNames_.size(); itrig++) {
      unsigned int triggerindx = hltConfig_.triggerIndex(trigNames_[itrig]);
      if (triggerindx >= n) {
        edm::LogWarning("HcalIsoTrack") << trigNames_[itrig] << " " << triggerindx << " does not exist in "
                                        << "the current menu";
      } else {
        edm::LogInfo("HcalIsoTrack") << trigNames_[itrig] << " " << triggerindx << " exists";
      }
    }
  }
}

void IsoTrackCalibration::endLuminosityBlock ( edm::LuminosityBlock const &  , edm::EventSetup const &    )

privatevirtual

Definition at line 523 of file IsoTrackCalibration.cc.

{}

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

overrideprivate

Definition at line 515 of file IsoTrackCalibration.cc.

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

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

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

static

Definition at line 525 of file IsoTrackCalibration.cc.

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

                                                                                     {
  //The following says we do not know what parameters are allowed so do no validation
  // Please change this to state exactly what you do use, even if it is no parameters
  edm::ParameterSetDescription desc;
  desc.setUnknown();
  descriptions.addDefault(desc);
}

Member Data Documentation

double IsoTrackCalibration::a_charIsoR_

private

Definition at line 95 of file IsoTrackCalibration.cc.

Referenced by analyze(), and IsoTrackCalibration().

double IsoTrackCalibration::a_coneR_

private

Definition at line 95 of file IsoTrackCalibration.cc.

Referenced by analyze(), and IsoTrackCalibration().

double IsoTrackCalibration::a_mipR_

private

Definition at line 95 of file IsoTrackCalibration.cc.

Referenced by analyze(), and IsoTrackCalibration().

double IsoTrackCalibration::constTrackPt_

private

Definition at line 97 of file IsoTrackCalibration.cc.

Referenced by analyze(), and IsoTrackCalibration().

edm::Service<TFileService> IsoTrackCalibration::fs_

private

Definition at line 89 of file IsoTrackCalibration.cc.

Referenced by beginJob().

TH1I* IsoTrackCalibration::h_iEta

private

Definition at line 124 of file IsoTrackCalibration.cc.

TH1F * IsoTrackCalibration::h_jetp

private

Definition at line 126 of file IsoTrackCalibration.cc.

TH1F* IsoTrackCalibration::h_jetpt[4]

private

Definition at line 127 of file IsoTrackCalibration.cc.

TH1F* IsoTrackCalibration::h_nTrk

private

Definition at line 122 of file IsoTrackCalibration.cc.

Referenced by analyze(), and beginJob().

TH1F * IsoTrackCalibration::h_nVtx

private

Definition at line 122 of file IsoTrackCalibration.cc.

Referenced by analyze(), and beginJob().

TH1F* IsoTrackCalibration::h_Rechit_E

private

Definition at line 126 of file IsoTrackCalibration.cc.

TProfile* IsoTrackCalibration::h_RecHit_iEta

private

Definition at line 123 of file IsoTrackCalibration.cc.

TProfile * IsoTrackCalibration::h_RecHit_num

private

Definition at line 123 of file IsoTrackCalibration.cc.

TH1I * IsoTrackCalibration::h_tketa0[5]

private

Definition at line 124 of file IsoTrackCalibration.cc.

TH1I * IsoTrackCalibration::h_tketa1[5]

private

Definition at line 124 of file IsoTrackCalibration.cc.

TH1I * IsoTrackCalibration::h_tketa2[5]

private

Definition at line 124 of file IsoTrackCalibration.cc.

TH1I* IsoTrackCalibration::h_tketa3[5]

private

Definition at line 125 of file IsoTrackCalibration.cc.

TH1I * IsoTrackCalibration::h_tketa4[5]

private

Definition at line 125 of file IsoTrackCalibration.cc.

TH1I * IsoTrackCalibration::h_tketa5[5]

private

Definition at line 125 of file IsoTrackCalibration.cc.

TH1I* IsoTrackCalibration::h_tketav1[5][6]

private

Definition at line 128 of file IsoTrackCalibration.cc.

TH1I * IsoTrackCalibration::h_tketav2[5][6]

private

Definition at line 128 of file IsoTrackCalibration.cc.

HLTConfigProvider IsoTrackCalibration::hltConfig_

private

Definition at line 90 of file IsoTrackCalibration.cc.

Referenced by beginRun().

std::vector<std::string> IsoTrackCalibration::HLTNames_

private

Definition at line 91 of file IsoTrackCalibration.cc.

bool IsoTrackCalibration::isAOD_

private

Definition at line 96 of file IsoTrackCalibration.cc.

Referenced by IsoTrackCalibration().

bool IsoTrackCalibration::isMC_

private

Definition at line 96 of file IsoTrackCalibration.cc.

Referenced by analyze(), and IsoTrackCalibration().

bool IsoTrackCalibration::isQCD_

private

Definition at line 96 of file IsoTrackCalibration.cc.

Referenced by analyze(), and IsoTrackCalibration().

double IsoTrackCalibration::maxEcalEnr_

private

Definition at line 98 of file IsoTrackCalibration.cc.

Referenced by analyze(), and IsoTrackCalibration().

double IsoTrackCalibration::maxNeighborTrackEnr_

private

Definition at line 99 of file IsoTrackCalibration.cc.

Referenced by analyze(), and IsoTrackCalibration().

int IsoTrackCalibration::nRun_

private

Definition at line 100 of file IsoTrackCalibration.cc.

Referenced by beginRun(), and endRun().

spr::trackSelectionParameters IsoTrackCalibration::selectionParameters_

private

Definition at line 93 of file IsoTrackCalibration.cc.

Referenced by analyze(), and IsoTrackCalibration().

double IsoTrackCalibration::slopeTrackPt_

private

Definition at line 97 of file IsoTrackCalibration.cc.

Referenced by analyze(), and IsoTrackCalibration().

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

private

Definition at line 119 of file IsoTrackCalibration.cc.

Referenced by analyze(), and beginJob().

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

private

Definition at line 119 of file IsoTrackCalibration.cc.

Referenced by analyze(), and beginJob().

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

private

Definition at line 119 of file IsoTrackCalibration.cc.

Referenced by analyze(), and beginJob().

double IsoTrackCalibration::t_eHcal

private

Definition at line 116 of file IsoTrackCalibration.cc.

Referenced by analyze(), and beginJob().

double IsoTrackCalibration::t_eHcal10

private

Definition at line 116 of file IsoTrackCalibration.cc.

Referenced by analyze(), and beginJob().

double IsoTrackCalibration::t_eHcal30

private

Definition at line 116 of file IsoTrackCalibration.cc.

Referenced by analyze(), and beginJob().

double IsoTrackCalibration::t_eMipDR

private

Definition at line 117 of file IsoTrackCalibration.cc.

Referenced by analyze(), and beginJob().

int IsoTrackCalibration::t_Event

private

Definition at line 114 of file IsoTrackCalibration.cc.

Referenced by analyze(), and beginJob().

double IsoTrackCalibration::t_EventWeight

private

Definition at line 115 of file IsoTrackCalibration.cc.

Referenced by analyze(), and beginJob().

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

private

Definition at line 120 of file IsoTrackCalibration.cc.

Referenced by analyze(), and beginJob().

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

private

Definition at line 120 of file IsoTrackCalibration.cc.

Referenced by analyze(), and beginJob().

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

private

Definition at line 120 of file IsoTrackCalibration.cc.

Referenced by analyze(), and beginJob().

double IsoTrackCalibration::t_hmaxNearP

private

Definition at line 117 of file IsoTrackCalibration.cc.

Referenced by analyze(), and beginJob().

int IsoTrackCalibration::t_ieta

private

Definition at line 114 of file IsoTrackCalibration.cc.

Referenced by analyze(), and beginJob().

int IsoTrackCalibration::t_nTrk

private

Definition at line 114 of file IsoTrackCalibration.cc.

Referenced by analyze(), and beginJob().

int IsoTrackCalibration::t_nVtx

private

Definition at line 114 of file IsoTrackCalibration.cc.

Referenced by analyze(), and beginJob().

double IsoTrackCalibration::t_p

private

Definition at line 115 of file IsoTrackCalibration.cc.

Referenced by analyze(), and beginJob().

double IsoTrackCalibration::t_phi

private

Definition at line 115 of file IsoTrackCalibration.cc.

Referenced by analyze(), and beginJob().

double IsoTrackCalibration::t_pt

private

Definition at line 115 of file IsoTrackCalibration.cc.

Referenced by analyze(), and beginJob().

bool IsoTrackCalibration::t_qltyMissFlag

private

Definition at line 118 of file IsoTrackCalibration.cc.

Referenced by analyze(), and beginJob().

bool IsoTrackCalibration::t_qltyPVFlag

private

Definition at line 118 of file IsoTrackCalibration.cc.

Referenced by analyze(), and beginJob().

int IsoTrackCalibration::t_Run

private

Definition at line 114 of file IsoTrackCalibration.cc.

Referenced by analyze(), and beginJob().

bool IsoTrackCalibration::t_selectTk

private

Definition at line 118 of file IsoTrackCalibration.cc.

Referenced by analyze(), and beginJob().

std::string IsoTrackCalibration::theTrackQuality_

private

Definition at line 94 of file IsoTrackCalibration.cc.

Referenced by analyze(), and IsoTrackCalibration().

edm::InputTag IsoTrackCalibration::theTriggerResultsLabel_

private

Definition at line 101 of file IsoTrackCalibration.cc.

Referenced by IsoTrackCalibration().

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

private

Definition at line 107 of file IsoTrackCalibration.cc.

Referenced by analyze(), and IsoTrackCalibration().

edm::EDGetTokenT<EcalRecHitCollection> IsoTrackCalibration::tok_EB_

private

Definition at line 108 of file IsoTrackCalibration.cc.

Referenced by analyze(), and IsoTrackCalibration().

edm::EDGetTokenT<EcalRecHitCollection> IsoTrackCalibration::tok_EE_

private

Definition at line 109 of file IsoTrackCalibration.cc.

Referenced by analyze(), and IsoTrackCalibration().

edm::EDGetTokenT<GenEventInfoProduct> IsoTrackCalibration::tok_ew_

private

Definition at line 111 of file IsoTrackCalibration.cc.

Referenced by analyze(), and IsoTrackCalibration().

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

private

Definition at line 105 of file IsoTrackCalibration.cc.

Referenced by analyze(), and IsoTrackCalibration().

edm::EDGetTokenT<HBHERecHitCollection> IsoTrackCalibration::tok_hbhe_

private

Definition at line 110 of file IsoTrackCalibration.cc.

Referenced by analyze(), and IsoTrackCalibration().

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

private

Definition at line 106 of file IsoTrackCalibration.cc.

Referenced by analyze(), and IsoTrackCalibration().

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

private

Definition at line 102 of file IsoTrackCalibration.cc.

Referenced by analyze(), and IsoTrackCalibration().

edm::EDGetTokenT<edm::TriggerResults> IsoTrackCalibration::tok_trigRes_

private

Definition at line 103 of file IsoTrackCalibration.cc.

Referenced by analyze(), and IsoTrackCalibration().

TTree* IsoTrackCalibration::tree

private

Definition at line 113 of file IsoTrackCalibration.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 IsoTrackCalibration::triggerEvent_

private

Definition at line 101 of file IsoTrackCalibration.cc.

Referenced by analyze(), and IsoTrackCalibration().

std::vector<std::string> IsoTrackCalibration::trigNames_

private

Definition at line 91 of file IsoTrackCalibration.cc.

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

int IsoTrackCalibration::verbosity_

private

Definition at line 92 of file IsoTrackCalibration.cc.

Referenced by analyze(), and IsoTrackCalibration().