EvtPlaneProducer Class Reference

#include <RecoHI/EvtPlaneProducer/src/EvtPlaneProducer.cc>

Inheritance diagram for EvtPlaneProducer:

Public Member Functions
	EvtPlaneProducer (const edm::ParameterSet &)
	~EvtPlaneProducer () override
Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default
	EDProducer (const EDProducer &)=delete
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () const final
bool	hasAbilityToProduceInEndRuns () const final
const EDProducer &	operator= (const EDProducer &)=delete

Private Member Functions
void	fillCastor (const TrackStructure &track, double vz, int bin)
void	fillHF (const TrackStructure &track, double vz, int bin)
void	fillTracker (const TrackStructure &track, double vz, int bin)
bool	passEta (float eta, int i)
void	produce (edm::Event &, const edm::EventSetup &) override

Private Attributes
bool	bScalo_particleFlow_
bool	bStrack_packedPFCandidates_
double	caloCentRef_
double	caloCentRefWidth_
edm::Handle< CaloTowerCollection >	caloCollection_
edm::InputTag	caloTag_
edm::EDGetTokenT< CaloTowerCollection >	caloToken_
edm::EDGetTokenT< reco::PFCandidateCollection >	caloTokenPF_
edm::Handle< std::vector< reco::CastorTower > >	castorCollection_
edm::InputTag	castorTag_
edm::EDGetTokenT< std::vector< reco::CastorTower > >	castorToken_
int	CentBinCompression_
edm::InputTag	centralityBinTag_
edm::EDGetTokenT< int >	centralityBinToken_
std::string	centralityMC_
edm::ESGetToken< CentralityTable, HeavyIonRcd >	centralityToken_
std::string	centralityVariable_
double	chi2_
edm::InputTag	chi2MapLostTag_
edm::EDGetTokenT< edm::ValueMap< float > >	chi2MapLostToken_
edm::InputTag	chi2MapTag_
edm::EDGetTokenT< edm::ValueMap< float > >	chi2MapToken_
double	chi2perlayer_
int	cutEra_
EPCuts	cuts_
double	d0d0error_
double	dzdzerror_
double	dzdzerror_pix_
double	dzerr_
HiEvtPlaneFlatten *	flat [NumEPNames]
double	flatdelvtx_
double	flatminvtx_
int	flatnvtxbins_
int	FlatOrder_
edm::ESGetToken< RPFlatParams, HeavyIonRPRcd >	flatparmsToken_
edm::ESWatcher< HeavyIonRPRcd >	hirpWatcher_
edm::ESWatcher< HeavyIonRcd >	hiWatcher_
bool	loadDB_
edm::EDGetTokenT< pat::PackedCandidateCollection >	lostToken_
edm::InputTag	losttrackTag_
double	maxet_
double	maxpt_
double	minet_
double	minpt_
double	nCentBins_
int	nhitsValid_
int	NumFlatBins_
edm::EDGetTokenT< pat::PackedCandidateCollection >	packedToken_
double	pterror_
GenPlane *	rp [NumEPNames]
TrackStructure	track_
edm::Handle< reco::TrackCollection >	trackCollection_
edm::InputTag	trackTag_
edm::EDGetTokenT< reco::TrackCollection >	trackToken_
edm::Handle< std::vector< reco::Vertex > >	vertex_
edm::InputTag	vertexTag_
edm::EDGetTokenT< std::vector< reco::Vertex > >	vertexToken_

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
using	CacheTypes = CacheContexts< T... >
using	GlobalCache = typename CacheTypes::GlobalCache
using	HasAbility = AbilityChecker< T... >
using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache
using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache
using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >
using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache
using	RunCache = typename CacheTypes::RunCache
using	RunContext = RunContextT< RunCache, GlobalCache >
using	RunSummaryCache = typename CacheTypes::RunSummaryCache

Detailed Description

Description: <one line="" class="" summary>="">

Implementation: <Notes on="" implementation>="">

Definition at line 167 of file EvtPlaneProducer.cc.

Constructor & Destructor Documentation

EvtPlaneProducer()

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

explicit

Definition at line 345 of file EvtPlaneProducer.cc.

References bScalo_particleFlow_, bStrack_packedPFCandidates_, caloTag_, caloToken_, caloTokenPF_, castorTag_, castorToken_, centralityBinTag_, centralityBinToken_, centralityMC_, centralityToken_, centralityVariable_, chi2_, chi2MapLostTag_, chi2MapLostToken_, chi2MapTag_, chi2MapToken_, chi2perlayer_, edm::errors::Configuration, cutEra_, cuts_, d0d0error_, dzdzerror_, dzdzerror_pix_, hi::EPEtaMax1, hi::EPEtaMax2, hi::EPEtaMin1, hi::EPEtaMin2, hi::EPNames, hi::EPOrder, deDxTools::esConsumes(), Exception, edm::ParameterSet::exists(), flat, flatdelvtx_, flatminvtx_, flatnvtxbins_, FlatOrder_, flatparmsToken_, edm::ParameterSet::getParameter(), mps_fire::i, HiEvtPlaneFlatten::init(), edm::InputTag::label(), loadDB_, lostToken_, losttrackTag_, nCentBins_, nhitsValid_, hi::NumEPNames, NumFlatBins_, packedToken_, pterror_, rp, AlCaHLTBitMon_QueryRunRegistry::string, trackTag_, trackToken_, vertexTag_, and vertexToken_.

    : centralityVariable_(iConfig.getParameter<std::string>("centralityVariable")),
      centralityBinTag_(iConfig.getParameter<edm::InputTag>("centralityBinTag")),
      vertexTag_(iConfig.getParameter<edm::InputTag>("vertexTag")),
      caloTag_(iConfig.getParameter<edm::InputTag>("caloTag")),
      castorTag_(iConfig.getParameter<edm::InputTag>("castorTag")),
      trackTag_(iConfig.getParameter<edm::InputTag>("trackTag")),
      losttrackTag_(iConfig.getParameter<edm::InputTag>("lostTag")),
      chi2MapTag_(iConfig.getParameter<edm::InputTag>("chi2MapTag")),
      chi2MapLostTag_(iConfig.getParameter<edm::InputTag>("chi2MapLostTag")),
      loadDB_(iConfig.getParameter<bool>("loadDB")),
      minet_(iConfig.getParameter<double>("minet")),
      maxet_(iConfig.getParameter<double>("maxet")),
      minpt_(iConfig.getParameter<double>("minpt")),
      maxpt_(iConfig.getParameter<double>("maxpt")),
      flatnvtxbins_(iConfig.getParameter<int>("flatnvtxbins")),
      flatminvtx_(iConfig.getParameter<double>("flatminvtx")),
      flatdelvtx_(iConfig.getParameter<double>("flatdelvtx")),
      dzdzerror_(iConfig.getParameter<double>("dzdzerror")),
      d0d0error_(iConfig.getParameter<double>("d0d0error")),
      pterror_(iConfig.getParameter<double>("pterror")),
      chi2perlayer_(iConfig.getParameter<double>("chi2perlayer")),
      dzdzerror_pix_(iConfig.getParameter<double>("dzdzerror_pix")),
      chi2_(iConfig.getParameter<double>("chi2")),
      nhitsValid_(iConfig.getParameter<int>("nhitsValid")),
      FlatOrder_(iConfig.getParameter<int>("FlatOrder")),
      NumFlatBins_(iConfig.getParameter<int>("NumFlatBins")),
      caloCentRef_(iConfig.getParameter<double>("caloCentRef")),
      caloCentRefWidth_(iConfig.getParameter<double>("caloCentRefWidth")),
      CentBinCompression_(iConfig.getParameter<int>("CentBinCompression")),
      cutEra_(iConfig.getParameter<int>("cutEra"))

{
  if (cutEra_ > 3)
    throw edm::Exception(edm::errors::Configuration) << "wrong range in cutEra parameter";
  cuts_ = EPCuts(
      static_cast<EP_ERA>(cutEra_), pterror_, dzdzerror_, d0d0error_, chi2perlayer_, dzdzerror_pix_, chi2_, nhitsValid_);
  nCentBins_ = 200.;

  if (iConfig.exists("nonDefaultGlauberModel")) {
    centralityMC_ = iConfig.getParameter<std::string>("nonDefaultGlauberModel");
  }
  centralityToken_ = esConsumes(edm::ESInputTag("", centralityVariable_ + centralityMC_));
  if (loadDB_) {
    flatparmsToken_ = esConsumes();
  }

  centralityBinToken_ = consumes<int>(centralityBinTag_);

  vertexToken_ = consumes<std::vector<reco::Vertex>>(vertexTag_);

  bStrack_packedPFCandidates_ = (trackTag_.label().find("packedPFCandidates") != std::string::npos);
  bScalo_particleFlow_ = (caloTag_.label().find("particleFlow") != std::string::npos);
  if (bStrack_packedPFCandidates_) {
    packedToken_ = consumes<pat::PackedCandidateCollection>(trackTag_);
    lostToken_ = consumes<pat::PackedCandidateCollection>(losttrackTag_);
    chi2MapToken_ = consumes<edm::ValueMap<float>>(chi2MapTag_);
    chi2MapLostToken_ = consumes<edm::ValueMap<float>>(chi2MapLostTag_);

  } else {
    if (bScalo_particleFlow_) {
      caloTokenPF_ = consumes<reco::PFCandidateCollection>(caloTag_);
    } else {
      caloToken_ = consumes<CaloTowerCollection>(caloTag_);
    }
    castorToken_ = consumes<std::vector<reco::CastorTower>>(castorTag_);
    trackToken_ = consumes<reco::TrackCollection>(trackTag_);
  }

  produces<reco::EvtPlaneCollection>();
  for (int i = 0; i < NumEPNames; i++) {
    rp[i] = new GenPlane(EPNames[i], EPEtaMin1[i], EPEtaMax1[i], EPEtaMin2[i], EPEtaMax2[i], EPOrder[i]);
  }
  for (int i = 0; i < NumEPNames; i++) {
    flat[i] = new HiEvtPlaneFlatten();
    flat[i]->init(FlatOrder_, NumFlatBins_, flatnvtxbins_, flatminvtx_, flatdelvtx_, EPNames[i], EPOrder[i]);
  }
}

~EvtPlaneProducer()

EvtPlaneProducer::~EvtPlaneProducer ( )

override

Definition at line 424 of file EvtPlaneProducer.cc.

References flat, mps_fire::i, and hi::NumEPNames.

                                    {
  // do anything here that needs to be done at desctruction time
  // (e.g. close files, deallocate resources etc.)
  for (int i = 0; i < NumEPNames; i++) {
    delete flat[i];
  }
}

Member Function Documentation

fillCastor()

void EvtPlaneProducer::fillCastor ( const TrackStructure &  track, double  vz, int  bin  )

inlineprivate

Definition at line 273 of file EvtPlaneProducer.cc.

References hi::GenPlane::addParticle(), newFWLiteAna::bin, funct::cos(), hi::EPDet, hi::EPOrder, HiEvtPlaneFlatten::getW(), mps_fire::i, HiEvtPlane_cfi::maxet, hi::maxTransverse, HiEvtPlane_cfi::minet, hi::minTransverse, hi::MomConsWeight, hi::NumEPNames, funct::sin(), HLT_2024v11_cff::track, btvMC_cff::vz, and w().

Referenced by produce().

                                                                   {
    double minet = minet_;
    double maxet = maxet_;
    for (int i = 0; i < NumEPNames; i++) {
      if (EPDet[i] == Castor) {
        if (minet_ < 0)
          minet = minTransverse[i];
        if (maxet_ < 0)
          maxet = maxTransverse[i];
        if (track.et < minet)
          continue;
        if (track.et > maxet)
          continue;
        if (not passEta(track.eta, i))
          continue;
        double w = track.et;
        if (EPOrder[i] == 1) {
          if (MomConsWeight[i][0] == 'y' && loadDB_) {
            w = flat[i]->getW(track.et, vz, bin);
          }
        }
        rp[i]->addParticle(w, track.et, sin(EPOrder[i] * track.phi), cos(EPOrder[i] * track.phi), track.eta);
      }
    }
  }

fillHF()

void EvtPlaneProducer::fillHF ( const TrackStructure &  track, double  vz, int  bin  )

inlineprivate

Definition at line 245 of file EvtPlaneProducer.cc.

References hi::GenPlane::addParticle(), newFWLiteAna::bin, funct::cos(), hi::EPDet, hi::EPOrder, HiEvtPlaneFlatten::etScale(), HiEvtPlaneFlatten::getW(), tpScales_cff::HF, mps_fire::i, HiEvtPlane_cfi::maxet, hi::maxTransverse, HiEvtPlane_cfi::minet, hi::minTransverse, hi::MomConsWeight, hi::NumEPNames, funct::sin(), HLT_2024v11_cff::track, btvMC_cff::vz, and w().

Referenced by produce().

                                                               {
    double minet = minet_;
    double maxet = maxet_;
    for (int i = 0; i < NumEPNames; i++) {
      if (EPDet[i] != HF)
        continue;
      if (minet_ < 0)
        minet = minTransverse[i];
      if (maxet_ < 0)
        maxet = maxTransverse[i];
      if (track.et < minet)
        continue;
      if (track.et > maxet)
        continue;
      if (not passEta(track.eta, i))
        continue;
      double w = track.et;
      if (loadDB_)
        w = track.et * flat[i]->etScale(vz, bin);
      if (EPOrder[i] == 1) {
        if (MomConsWeight[i][0] == 'y' && loadDB_) {
          w = flat[i]->getW(track.et, vz, bin);
        }
      }
      rp[i]->addParticle(w, track.et, sin(EPOrder[i] * track.phi), cos(EPOrder[i] * track.phi), track.eta);
    }
  };

fillTracker()

void EvtPlaneProducer::fillTracker ( const TrackStructure &  track, double  vz, int  bin  )

inlineprivate

Definition at line 316 of file EvtPlaneProducer.cc.

References hi::GenPlane::addParticle(), newFWLiteAna::bin, funct::cos(), hi::EPDet, hi::EPOrder, HiEvtPlaneFlatten::getW(), mps_fire::i, HiEvtPlane_cfi::maxpt, hi::maxTransverse, HiEvtPlane_cfi::minpt, hi::minTransverse, hi::MomConsWeight, hi::NumEPNames, funct::sin(), HLT_2024v11_cff::track, align::Tracker, btvMC_cff::vz, and w().

Referenced by produce().

                                                                    {
    double minpt = minpt_;
    double maxpt = maxpt_;
    for (int i = 0; i < NumEPNames; i++) {
      if (EPDet[i] == Tracker) {
        if (minpt_ < 0)
          minpt = minTransverse[i];
        if (maxpt_ < 0)
          maxpt = maxTransverse[i];
        if (track.pt < minpt)
          continue;
        if (track.pt > maxpt)
          continue;
        if (not passEta(track.eta, i))
          continue;
        double w = track.pt;
        if (w > 2.5)
          w = 2.0;  //v2 starts decreasing above ~2.5 GeV/c
        if (EPOrder[i] == 1) {
          if (MomConsWeight[i][0] == 'y' && loadDB_) {
            w = flat[i]->getW(track.pt, vz, bin);
          }
        }
        rp[i]->addParticle(w, track.pt, sin(EPOrder[i] * track.phi), cos(EPOrder[i] * track.phi), track.eta);
      }
    }
  };

passEta()

bool EvtPlaneProducer::passEta ( float  eta, int  i  )

inlineprivate

Definition at line 299 of file EvtPlaneProducer.cc.

References hi::EPEtaMax1, hi::EPEtaMax2, hi::EPEtaMin1, hi::EPEtaMin2, PVValHelper::eta, and mps_fire::i.

                                 {
    if (EPEtaMin2[i] == EPEtaMax2[i]) {
      if (eta < EPEtaMin1[i])
        return false;
      if (eta > EPEtaMax1[i])
        return false;
    } else {
      if (eta < EPEtaMin1[i])
        return false;
      if (eta > EPEtaMax2[i])
        return false;
      if (eta > EPEtaMax1[i] && eta < EPEtaMin2[i])
        return false;
    }
    return true;
  }

produce()

void EvtPlaneProducer::produce ( edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

overrideprivate

Definition at line 437 of file EvtPlaneProducer.cc.

References funct::abs(), reco::TrackBase::algo(), hi::TrackStructure::algos, newFWLiteAna::bin, bScalo_particleFlow_, bStrack_packedPFCandidates_, caloCentRef_, caloCentRefWidth_, caloCollection_, caloToken_, caloTokenPF_, HLT_2024v11_cff::cands, castorCollection_, castorToken_, hi::TrackStructure::centbin, CentBinCompression_, centralityBinToken_, centralityToken_, hi::TrackStructure::charge, reco::TrackBase::charge(), edm::ESWatcher< T >::check(), hi::TrackStructure::chi2layer, chi2MapLostToken_, chi2MapToken_, hi::TrackStructure::collection, cuts_, cuy::cv, dqmiodatasetharvest::db, hi::TrackStructure::dxy, reco::TrackBase::dxy(), hi::TrackStructure::dxyError, reco::TrackBase::dxyError(), hi::TrackStructure::dxySig, hi::TrackStructure::dz, reco::TrackBase::dz(), hi::TrackStructure::dzError, reco::TrackBase::dzError(), hi::TrackStructure::dzSig, hi::EPDet, hi::TrackStructure::et, hi::TrackStructure::eta, reco::TrackBase::eta(), fillCastor(), fillHF(), fillTracker(), flat, flatparmsToken_, hi::GenPlane::getAngle(), edm::EventSetup::getData(), tpScales_cff::HF, hi::TrackStructure::highPurity, reco::TrackBase::highPurity, hirpWatcher_, reco::TrackBase::hitPattern(), hiWatcher_, Matriplex::hypot(), mps_fire::i, heavyIonCSV_trainingSettings::idx, iEvent, hi::EPCuts::isGoodCastor(), hi::EPCuts::isGoodHF(), hi::EPCuts::isGoodTrack(), loadDB_, lostToken_, eostools::move(), dqmiodumpmetadata::n, nCentBins_, hi::TrackStructure::normalizedChi2, hi::TrackStructure::numberOfValidHits, reco::TrackBase::numberOfValidHits(), hi::NumEPNames, packedToken_, hi::TrackStructure::pdgid, packedPFCandidateRefMixer_cfi::pf, hi::TrackStructure::phi, reco::TrackBase::phi(), hi::TrackStructure::pt, reco::TrackBase::pt(), hi::TrackStructure::ptError, reco::TrackBase::ptError(), reset(), rp, HiEvtPlaneFlatten::setCaloCentRefBins(), pfDeepBoostedJetPreprocessParams_cfi::sv, l1tHGCalTowerProducer_cfi::tower, track_, trackCollection_, reco::HitPattern::trackerLayersWithMeasurement(), trackToken_, parallelization::uint, vertexToken_, L1BJetProducer_cff::vtx, and ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::wv.

                                                                            {
  using namespace edm;
  using namespace std;
  using namespace reco;
  if (hiWatcher_.check(iSetup)) {
    //
    //Get Size of Centrality Table
    //
    auto const &centDB = iSetup.getData(centralityToken_);
    nCentBins_ = centDB.m_table.size();
    for (int i = 0; i < NumEPNames; i++) {
      if (caloCentRef_ > 0) {
        int minbin = (caloCentRef_ - caloCentRefWidth_ / 2.) * nCentBins_ / 100.;
        int maxbin = (caloCentRef_ + caloCentRefWidth_ / 2.) * nCentBins_ / 100.;
        minbin /= CentBinCompression_;
        maxbin /= CentBinCompression_;
        if (minbin > 0 && maxbin >= minbin) {
          if (EPDet[i] == HF || EPDet[i] == Castor)
            flat[i]->setCaloCentRefBins(minbin, maxbin);
        }
      }
    }
  }
  //
  //Get flattening parameter file.
  //
  if (loadDB_ && hirpWatcher_.check(iSetup)) {
    LoadEPDB db(iSetup.getData(flatparmsToken_), flat);
    if (!db.IsSuccess()) {
      loadDB_ = kFALSE;
    }
  }
  //
  //Get Centrality
  //
  int bin = 0;
  int cbin = 0;
  if (loadDB_) {
    cbin = iEvent.get(centralityBinToken_);
    bin = cbin / CentBinCompression_;
  }
  //
  //Get Vertex
  //
  //best vertex
  const reco::Vertex &vtx = iEvent.get(vertexToken_)[0];
  double bestvz = vtx.z();
  double bestvx = vtx.x();
  double bestvy = vtx.y();
  double bestvzError = vtx.zError();
  math::XYZPoint bestvtx(bestvx, bestvy, bestvz);
  math::Error<3>::type vtx_cov = vtx.covariance();

  for (int i = 0; i < NumEPNames; i++)
    rp[i]->reset();
  edm::Handle<edm::ValueMap<float>> chi2Map;
  edm::Handle<pat::PackedCandidateCollection> cands;
  edm::Handle<reco::PFCandidateCollection> calocands;
  if (bStrack_packedPFCandidates_) {
    for (int idx = 1; idx < 3; idx++) {
      if (idx == 1) {
        iEvent.getByToken(packedToken_, cands);
        iEvent.getByToken(chi2MapToken_, chi2Map);
      }
      if (idx == 2) {
        iEvent.getByToken(lostToken_, cands);
        iEvent.getByToken(chi2MapLostToken_, chi2Map);
      }
      for (unsigned int i = 0, n = cands->size(); i < n; ++i) {
        track_ = {};
        track_.centbin = cbin;
        const pat::PackedCandidate &pf = (*cands)[i];
        track_.et = pf.et();
        track_.eta = pf.eta();
        track_.phi = pf.phi();
        track_.pdgid = pf.pdgId();
        if ((idx == 1) and cuts_.isGoodHF(track_)) {
          fillHF(track_, bestvz, bin);
        }
        if (!pf.hasTrackDetails())
          continue;
        const reco::Track &trk = pf.pseudoTrack();
        track_.highPurity = pf.trackHighPurity();
        track_.charge = trk.charge();
        if (!track_.highPurity || track_.charge == 0)
          continue;
        track_.collection = idx;
        track_.eta = trk.eta();
        track_.phi = trk.phi();
        track_.pt = trk.pt();
        track_.ptError = trk.ptError();
        track_.numberOfValidHits = trk.numberOfValidHits();
        track_.algos = trk.algo();
        track_.dz = std::abs(trk.dz(bestvtx));
        track_.dxy = std::abs(trk.dxy(bestvtx));
        track_.dzError = std::hypot(trk.dzError(), bestvzError);
        track_.dxyError = trk.dxyError(bestvtx, vtx_cov);
        track_.dzSig = track_.dz / track_.dzError;
        track_.dxySig = track_.dxy / track_.dxyError;
        const reco::HitPattern &hit_pattern = trk.hitPattern();
        track_.normalizedChi2 = (*chi2Map)[pat::PackedCandidateRef(cands, i)];
        track_.chi2layer = (*chi2Map)[pat::PackedCandidateRef(cands, i)] / hit_pattern.trackerLayersWithMeasurement();
        if (cuts_.isGoodTrack(track_)) {
          fillTracker(track_, bestvz, bin);
        }
      }
    }
  } else {
    //calorimetry part
    if (bScalo_particleFlow_) {
      iEvent.getByToken(caloTokenPF_, calocands);
      for (unsigned int i = 0, n = calocands->size(); i < n; ++i) {
        track_ = {};
        track_.centbin = cbin;
        const reco::PFCandidate &pf = (*calocands)[i];
        track_.et = pf.et();
        track_.eta = pf.eta();
        track_.phi = pf.phi();
        track_.pdgid = pf.pdgId();
        if (cuts_.isGoodHF(track_)) {
          fillHF(track_, bestvz, bin);
        }
      }
    } else {
      iEvent.getByToken(caloToken_, caloCollection_);
      for (const auto &tower : *caloCollection_) {
        track_.eta = tower.eta();
        track_.phi = tower.phi();
        track_.et = tower.emEt() + tower.hadEt();
        track_.pdgid = 1;
        if (cuts_.isGoodHF(track_))
          fillHF(track_, bestvz, bin);
      }
    }

    //Castor part
    iEvent.getByToken(castorToken_, castorCollection_);
    for (const auto &tower : *castorCollection_) {
      track_.eta = tower.eta();
      track_.phi = tower.phi();
      track_.et = tower.et();
      track_.pdgid = 1;
      if (cuts_.isGoodCastor(track_))
        fillCastor(track_, bestvz, bin);
    }
    //Tracking part
    iEvent.getByToken(trackToken_, trackCollection_);
    for (const auto &trk : *trackCollection_) {
      track_.highPurity = trk.quality(reco::TrackBase::highPurity);
      track_.charge = trk.charge();
      if (!track_.highPurity || track_.charge == 0)
        continue;
      track_.centbin = cbin;
      track_.collection = 0;
      track_.eta = trk.eta();
      track_.phi = trk.phi();
      track_.pt = trk.pt();
      track_.ptError = trk.ptError();
      track_.numberOfValidHits = trk.numberOfValidHits();
      track_.algos = trk.algo();
      track_.dz = std::abs(trk.dz(bestvtx));
      track_.dxy = std::abs(trk.dxy(bestvtx));
      track_.dzError = std::hypot(trk.dzError(), bestvzError);
      track_.dxyError = trk.dxyError(bestvtx, vtx_cov);
      track_.dzSig = track_.dz / track_.dzError;
      track_.dxySig = track_.dxy / track_.dxyError;
      track_.normalizedChi2 = trk.normalizedChi2();
      track_.chi2layer = track_.normalizedChi2 / trk.hitPattern().trackerLayersWithMeasurement();
      if (cuts_.isGoodTrack(track_))
        fillTracker(track_, bestvz, bin);
    }
  }

  auto evtplaneOutput = std::make_unique<EvtPlaneCollection>();

  double ang = -10;
  double sv = 0;
  double cv = 0;
  double svNoWgt = 0;
  double cvNoWgt = 0;

  double wv = 0;
  double wv2 = 0;
  double pe = 0;
  double pe2 = 0;
  uint epmult = 0;

  for (int i = 0; i < NumEPNames; i++) {
    rp[i]->getAngle(ang, sv, cv, svNoWgt, cvNoWgt, wv, wv2, pe, pe2, epmult);
    evtplaneOutput->push_back(EvtPlane(i, 0, ang, sv, cv, wv, wv2, pe, pe2, epmult));
    evtplaneOutput->back().addLevel(3, 0., svNoWgt, cvNoWgt);
  }

  iEvent.put(std::move(evtplaneOutput));
}

Member Data Documentation

bScalo_particleFlow_

bool EvtPlaneProducer::bScalo_particleFlow_

private

Definition at line 204 of file EvtPlaneProducer.cc.

Referenced by EvtPlaneProducer(), and produce().

bStrack_packedPFCandidates_

bool EvtPlaneProducer::bStrack_packedPFCandidates_

private

Definition at line 203 of file EvtPlaneProducer.cc.

Referenced by EvtPlaneProducer(), and produce().

caloCentRef_

double EvtPlaneProducer::caloCentRef_

private

Definition at line 235 of file EvtPlaneProducer.cc.

Referenced by produce().

caloCentRefWidth_

double EvtPlaneProducer::caloCentRefWidth_

private

Definition at line 236 of file EvtPlaneProducer.cc.

Referenced by produce().

caloCollection_

edm::Handle<CaloTowerCollection> EvtPlaneProducer::caloCollection_

private

Definition at line 192 of file EvtPlaneProducer.cc.

Referenced by produce().

caloTag_

edm::InputTag EvtPlaneProducer::caloTag_

private

Definition at line 190 of file EvtPlaneProducer.cc.

Referenced by EvtPlaneProducer().

caloToken_

edm::EDGetTokenT<CaloTowerCollection> EvtPlaneProducer::caloToken_

private

Definition at line 191 of file EvtPlaneProducer.cc.

Referenced by EvtPlaneProducer(), and produce().

caloTokenPF_

edm::EDGetTokenT<reco::PFCandidateCollection> EvtPlaneProducer::caloTokenPF_

private

Definition at line 193 of file EvtPlaneProducer.cc.

Referenced by EvtPlaneProducer(), and produce().

castorCollection_

edm::Handle<std::vector<reco::CastorTower> > EvtPlaneProducer::castorCollection_

private

Definition at line 197 of file EvtPlaneProducer.cc.

Referenced by produce().

castorTag_

edm::InputTag EvtPlaneProducer::castorTag_

private

Definition at line 195 of file EvtPlaneProducer.cc.

Referenced by EvtPlaneProducer().

castorToken_

edm::EDGetTokenT<std::vector<reco::CastorTower> > EvtPlaneProducer::castorToken_

private

Definition at line 196 of file EvtPlaneProducer.cc.

Referenced by EvtPlaneProducer(), and produce().

CentBinCompression_

int EvtPlaneProducer::CentBinCompression_

private

Definition at line 237 of file EvtPlaneProducer.cc.

Referenced by produce().

centralityBinTag_

edm::InputTag EvtPlaneProducer::centralityBinTag_

private

Definition at line 183 of file EvtPlaneProducer.cc.

Referenced by EvtPlaneProducer().

centralityBinToken_

edm::EDGetTokenT<int> EvtPlaneProducer::centralityBinToken_

private

Definition at line 184 of file EvtPlaneProducer.cc.

Referenced by EvtPlaneProducer(), and produce().

centralityMC_

std::string EvtPlaneProducer::centralityMC_

private

Definition at line 181 of file EvtPlaneProducer.cc.

Referenced by EvtPlaneProducer().

centralityToken_

edm::ESGetToken<CentralityTable, HeavyIonRcd> EvtPlaneProducer::centralityToken_

private

Definition at line 213 of file EvtPlaneProducer.cc.

Referenced by EvtPlaneProducer(), and produce().

centralityVariable_

std::string EvtPlaneProducer::centralityVariable_

private

Definition at line 180 of file EvtPlaneProducer.cc.

Referenced by EvtPlaneProducer().

chi2_

double EvtPlaneProducer::chi2_

private

Definition at line 230 of file EvtPlaneProducer.cc.

Referenced by EvtPlaneProducer().

chi2MapLostTag_

edm::InputTag EvtPlaneProducer::chi2MapLostTag_

private

Definition at line 210 of file EvtPlaneProducer.cc.

Referenced by EvtPlaneProducer().

chi2MapLostToken_

edm::EDGetTokenT<edm::ValueMap<float> > EvtPlaneProducer::chi2MapLostToken_

private

Definition at line 211 of file EvtPlaneProducer.cc.

Referenced by EvtPlaneProducer(), and produce().

chi2MapTag_

edm::InputTag EvtPlaneProducer::chi2MapTag_

private

Definition at line 208 of file EvtPlaneProducer.cc.

Referenced by EvtPlaneProducer().

chi2MapToken_

edm::EDGetTokenT<edm::ValueMap<float> > EvtPlaneProducer::chi2MapToken_

private

Definition at line 209 of file EvtPlaneProducer.cc.

Referenced by EvtPlaneProducer(), and produce().

chi2perlayer_

double EvtPlaneProducer::chi2perlayer_

private

Definition at line 227 of file EvtPlaneProducer.cc.

Referenced by EvtPlaneProducer().

cutEra_

int EvtPlaneProducer::cutEra_

private

Definition at line 238 of file EvtPlaneProducer.cc.

Referenced by EvtPlaneProducer().

cuts_

EPCuts EvtPlaneProducer::cuts_

private

Definition at line 178 of file EvtPlaneProducer.cc.

Referenced by EvtPlaneProducer(), and produce().

d0d0error_

double EvtPlaneProducer::d0d0error_

private

Definition at line 225 of file EvtPlaneProducer.cc.

Referenced by EvtPlaneProducer().

dzdzerror_

double EvtPlaneProducer::dzdzerror_

private

Definition at line 224 of file EvtPlaneProducer.cc.

Referenced by EvtPlaneProducer().

dzdzerror_pix_

double EvtPlaneProducer::dzdzerror_pix_

private

Definition at line 229 of file EvtPlaneProducer.cc.

Referenced by EvtPlaneProducer().

dzerr_

double EvtPlaneProducer::dzerr_

private

Definition at line 228 of file EvtPlaneProducer.cc.

flat

HiEvtPlaneFlatten* EvtPlaneProducer::flat[NumEPNames]

private

Definition at line 239 of file EvtPlaneProducer.cc.

Referenced by EvtPlaneProducer(), produce(), and ~EvtPlaneProducer().

flatdelvtx_

double EvtPlaneProducer::flatdelvtx_

private

Definition at line 223 of file EvtPlaneProducer.cc.

Referenced by EvtPlaneProducer().

flatminvtx_

double EvtPlaneProducer::flatminvtx_

private

Definition at line 222 of file EvtPlaneProducer.cc.

Referenced by EvtPlaneProducer().

flatnvtxbins_

int EvtPlaneProducer::flatnvtxbins_

private

Definition at line 221 of file EvtPlaneProducer.cc.

Referenced by EvtPlaneProducer().

FlatOrder_

int EvtPlaneProducer::FlatOrder_

private

Definition at line 232 of file EvtPlaneProducer.cc.

Referenced by EvtPlaneProducer().

flatparmsToken_

edm::ESGetToken<RPFlatParams, HeavyIonRPRcd> EvtPlaneProducer::flatparmsToken_

private

Definition at line 214 of file EvtPlaneProducer.cc.

Referenced by EvtPlaneProducer(), and produce().

hirpWatcher_

edm::ESWatcher<HeavyIonRPRcd> EvtPlaneProducer::hirpWatcher_

private

Definition at line 243 of file EvtPlaneProducer.cc.

Referenced by produce().

hiWatcher_

edm::ESWatcher<HeavyIonRcd> EvtPlaneProducer::hiWatcher_

private

Definition at line 242 of file EvtPlaneProducer.cc.

Referenced by produce().

loadDB_

bool EvtPlaneProducer::loadDB_

private

Definition at line 216 of file EvtPlaneProducer.cc.

Referenced by EvtPlaneProducer(), and produce().

lostToken_

edm::EDGetTokenT<pat::PackedCandidateCollection> EvtPlaneProducer::lostToken_

private

Definition at line 206 of file EvtPlaneProducer.cc.

Referenced by EvtPlaneProducer(), and produce().

losttrackTag_

edm::InputTag EvtPlaneProducer::losttrackTag_

private

Definition at line 201 of file EvtPlaneProducer.cc.

Referenced by EvtPlaneProducer().

maxet_

double EvtPlaneProducer::maxet_

private

Definition at line 218 of file EvtPlaneProducer.cc.

maxpt_

double EvtPlaneProducer::maxpt_

private

Definition at line 220 of file EvtPlaneProducer.cc.

minet_

double EvtPlaneProducer::minet_

private

Definition at line 217 of file EvtPlaneProducer.cc.

minpt_

double EvtPlaneProducer::minpt_

private

Definition at line 219 of file EvtPlaneProducer.cc.

nCentBins_

double EvtPlaneProducer::nCentBins_

private

Definition at line 234 of file EvtPlaneProducer.cc.

Referenced by EvtPlaneProducer(), and produce().

nhitsValid_

int EvtPlaneProducer::nhitsValid_

private

Definition at line 231 of file EvtPlaneProducer.cc.

Referenced by EvtPlaneProducer().

NumFlatBins_

int EvtPlaneProducer::NumFlatBins_

private

Definition at line 233 of file EvtPlaneProducer.cc.

Referenced by EvtPlaneProducer().

packedToken_

edm::EDGetTokenT<pat::PackedCandidateCollection> EvtPlaneProducer::packedToken_

private

Definition at line 205 of file EvtPlaneProducer.cc.

Referenced by EvtPlaneProducer(), and produce().

pterror_

double EvtPlaneProducer::pterror_

private

Definition at line 226 of file EvtPlaneProducer.cc.

Referenced by EvtPlaneProducer().

rp

GenPlane* EvtPlaneProducer::rp[NumEPNames]

private

Definition at line 173 of file EvtPlaneProducer.cc.

Referenced by EvtPlaneProducer(), and produce().

track_

TrackStructure EvtPlaneProducer::track_

private

Definition at line 240 of file EvtPlaneProducer.cc.

Referenced by produce().

trackCollection_

edm::Handle<reco::TrackCollection> EvtPlaneProducer::trackCollection_

private

Definition at line 202 of file EvtPlaneProducer.cc.

Referenced by produce().

trackTag_

edm::InputTag EvtPlaneProducer::trackTag_

private

Definition at line 199 of file EvtPlaneProducer.cc.

Referenced by EvtPlaneProducer().

trackToken_

edm::EDGetTokenT<reco::TrackCollection> EvtPlaneProducer::trackToken_

private

Definition at line 200 of file EvtPlaneProducer.cc.

Referenced by EvtPlaneProducer(), and produce().

vertex_

edm::Handle<std::vector<reco::Vertex> > EvtPlaneProducer::vertex_

private

Definition at line 188 of file EvtPlaneProducer.cc.

vertexTag_

edm::InputTag EvtPlaneProducer::vertexTag_

private

Definition at line 186 of file EvtPlaneProducer.cc.

Referenced by EvtPlaneProducer().

vertexToken_

edm::EDGetTokenT<std::vector<reco::Vertex> > EvtPlaneProducer::vertexToken_

private

Definition at line 187 of file EvtPlaneProducer.cc.

Referenced by EvtPlaneProducer(), and produce().