HiEvtPlaneFlatProducer Class Reference

Inheritance diagram for HiEvtPlaneFlatProducer:

Public Member Functions
	HiEvtPlaneFlatProducer (const edm::ParameterSet &)
	~HiEvtPlaneFlatProducer () 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	produce (edm::Event &, const edm::EventSetup &) override

Private Attributes
double	caloCentRef_
double	caloCentRefWidth_
int	CentBinCompression_
edm::InputTag	centralityBinTag_
edm::EDGetTokenT< int >	centralityBinToken_
edm::ESGetToken< CentralityTable, HeavyIonRcd >	centralityESToken_
std::string	centralityMC_
edm::InputTag	centralityTag_
edm::EDGetTokenT< reco::Centrality >	centralityToken_
std::string	centralityVariable_
HiEvtPlaneFlatten *	flat [NumEPNames]
double	flatdelvtx_
double	flatminvtx_
int	flatnvtxbins_
const int	FlatOrder_
edm::ESGetToken< RPFlatParams, HeavyIonRPRcd >	flatparmsToken_
edm::ESWatcher< HeavyIonRPRcd >	hirpWatcher
edm::ESWatcher< HeavyIonRcd >	hiWatcher
edm::InputTag	inputPlanesTag_
edm::EDGetTokenT< reco::EvtPlaneCollection >	inputPlanesToken_
double	nCentBins_
int	NumFlatBins_
edm::Handle< reco::TrackCollection >	trackCollection_
edm::InputTag	trackTag_
edm::EDGetTokenT< reco::TrackCollection >	trackToken_
bool	useOffsetPsi_
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

Definition at line 48 of file HiEvtPlaneFlatProducer.cc.

Constructor & Destructor Documentation

HiEvtPlaneFlatProducer()

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

explicit

Definition at line 97 of file HiEvtPlaneFlatProducer.cc.

References centralityBinTag_, centralityBinToken_, centralityESToken_, centralityMC_, centralityTag_, centralityToken_, centralityVariable_, hi::EPNames, hi::EPOrder, deDxTools::esConsumes(), edm::ParameterSet::exists(), flat, flatdelvtx_, flatminvtx_, flatnvtxbins_, FlatOrder_, flatparmsToken_, edm::ParameterSet::getParameter(), mps_fire::i, HiEvtPlaneFlatten::init(), inputPlanesTag_, inputPlanesToken_, nCentBins_, hi::NumEPNames, NumFlatBins_, AlCaHLTBitMon_QueryRunRegistry::string, trackTag_, trackToken_, vertexTag_, and vertexToken_.

    : centralityVariable_(iConfig.getParameter<std::string>("centralityVariable")),
      centralityBinTag_(iConfig.getParameter<edm::InputTag>("centralityBinTag")),
      centralityTag_(iConfig.getParameter<edm::InputTag>("centralityTag")),
      vertexTag_(iConfig.getParameter<edm::InputTag>("vertexTag")),
      inputPlanesTag_(iConfig.getParameter<edm::InputTag>("inputPlanesTag")),
      trackTag_(iConfig.getParameter<edm::InputTag>("trackTag")),
      FlatOrder_(iConfig.getParameter<int>("FlatOrder")),
      NumFlatBins_(iConfig.getParameter<int>("NumFlatBins")),
      flatnvtxbins_(iConfig.getParameter<int>("flatnvtxbins")),
      flatminvtx_(iConfig.getParameter<double>("flatminvtx")),
      flatdelvtx_(iConfig.getParameter<double>("flatdelvtx")),
      caloCentRef_(iConfig.getParameter<double>("caloCentRef")),
      caloCentRefWidth_(iConfig.getParameter<double>("caloCentRefWidth")),
      CentBinCompression_(iConfig.getParameter<int>("CentBinCompression")),
      useOffsetPsi_(iConfig.getParameter<bool>("useOffsetPsi")) {
  nCentBins_ = 200.;

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

  centralityBinToken_ = consumes<int>(centralityBinTag_);

  centralityToken_ = consumes<reco::Centrality>(centralityTag_);

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

  trackToken_ = consumes<reco::TrackCollection>(trackTag_);

  inputPlanesToken_ = consumes<reco::EvtPlaneCollection>(inputPlanesTag_);

  //register your products
  produces<reco::EvtPlaneCollection>();
  //now do what ever other initialization is needed
  for (int i = 0; i < NumEPNames; i++) {
    flat[i] = new HiEvtPlaneFlatten();
    flat[i]->init(FlatOrder_, NumFlatBins_, flatnvtxbins_, flatminvtx_, flatdelvtx_, EPNames[i], EPOrder[i]);
  }
}

~HiEvtPlaneFlatProducer()

HiEvtPlaneFlatProducer::~HiEvtPlaneFlatProducer ( )

override

Definition at line 140 of file HiEvtPlaneFlatProducer.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

produce()

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

overrideprivate

Definition at line 153 of file HiEvtPlaneFlatProducer.cc.

References newFWLiteAna::bin, DummyCfis::c, caloCentRef_, caloCentRefWidth_, CentBinCompression_, centralityBinToken_, centralityESToken_, edm::ESWatcher< T >::check(), HiEvtPlaneFlatten::coffset(), dqmiodatasetharvest::db, SiStripBadComponentsDQMServiceTemplate_cfg::ep, hi::EPDet, flat, flatparmsToken_, edm::EventSetup::getData(), HiEvtPlaneFlatten::getFlatPsi(), tpScales_cff::HF, hirpWatcher, hiWatcher, mps_fire::i, iEvent, inputPlanesToken_, visualization-live-secondInstance_cfg::m, eostools::move(), nCentBins_, hi::NumEPNames, HiEvtPlaneFlatten::offsetPsi(), alignCSCRings::s, HiEvtPlaneFlatten::setCaloCentRefBins(), HiEvtPlaneFlatten::soffset(), parallelization::uint, useOffsetPsi_, vertexToken_, and L1BJetProducer_cff::vtx.

                                                                                  {
  using namespace edm;
  using namespace std;
  using namespace reco;

  if (hiWatcher.check(iSetup)) {
    //
    //Get Size of Centrality Table
    //
    auto const& centDB = iSetup.getData(centralityESToken_);
    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 (hirpWatcher.check(iSetup)) {
    LoadEPDB db(iSetup.getData(flatparmsToken_), flat);
  }  //rp record change

  //
  //Get Centrality
  //
  int bin = 0;
  int cbin = 0;
  cbin = iEvent.get(centralityBinToken_);
  bin = cbin / CentBinCompression_;
  //
  //Get Vertex
  //

  //best vertex
  double bestvz = -999.9;
  const reco::Vertex& vtx = iEvent.get(vertexToken_)[0];
  bestvz = vtx.z();

  //
  //Get Event Planes
  //

  auto const& evtPlanes = iEvent.get(inputPlanesToken_);

  auto evtplaneOutput = std::make_unique<EvtPlaneCollection>();
  EvtPlane* ep[NumEPNames];
  for (int i = 0; i < NumEPNames; i++) {
    ep[i] = nullptr;
  }
  int indx = 0;
  for (auto&& rp : (evtPlanes)) {
    double s = rp.sumSin(0);
    double c = rp.sumCos(0);
    uint m = rp.mult();
    double soff = s;
    double coff = c;
    double psiOffset = -10;
    double psiFlat = -10;
    if (rp.angle(0) > -5) {
      if (useOffsetPsi_) {
        soff = flat[indx]->soffset(s, bestvz, bin);
        coff = flat[indx]->coffset(c, bestvz, bin);
        psiOffset = flat[indx]->offsetPsi(soff, coff);
      }
      psiFlat = flat[indx]->getFlatPsi(psiOffset, bestvz, bin);
    }
    ep[indx] = new EvtPlane(indx, 2, psiFlat, soff, coff, rp.sumw(), rp.sumw2(), rp.sumPtOrEt(), rp.sumPtOrEt2(), m);
    ep[indx]->addLevel(0, rp.angle(0), s, c);
    ep[indx]->addLevel(3, 0., rp.sumSin(3), rp.sumCos(3));
    if (useOffsetPsi_)
      ep[indx]->addLevel(1, psiOffset, soff, coff);
    ++indx;
  }

  for (int i = 0; i < NumEPNames; i++) {
    if (ep[i] != nullptr)
      evtplaneOutput->push_back(*ep[i]);
  }
  iEvent.put(std::move(evtplaneOutput));
  for (int i = 0; i < indx; i++)
    delete ep[i];
}

Member Data Documentation

caloCentRef_

double HiEvtPlaneFlatProducer::caloCentRef_

private

Definition at line 87 of file HiEvtPlaneFlatProducer.cc.

Referenced by produce().

caloCentRefWidth_

double HiEvtPlaneFlatProducer::caloCentRefWidth_

private

Definition at line 88 of file HiEvtPlaneFlatProducer.cc.

Referenced by produce().

CentBinCompression_

int HiEvtPlaneFlatProducer::CentBinCompression_

private

Definition at line 89 of file HiEvtPlaneFlatProducer.cc.

Referenced by produce().

centralityBinTag_

edm::InputTag HiEvtPlaneFlatProducer::centralityBinTag_

private

Definition at line 61 of file HiEvtPlaneFlatProducer.cc.

Referenced by HiEvtPlaneFlatProducer().

centralityBinToken_

edm::EDGetTokenT<int> HiEvtPlaneFlatProducer::centralityBinToken_

private

Definition at line 62 of file HiEvtPlaneFlatProducer.cc.

Referenced by HiEvtPlaneFlatProducer(), and produce().

centralityESToken_

edm::ESGetToken<CentralityTable, HeavyIonRcd> HiEvtPlaneFlatProducer::centralityESToken_

private

Definition at line 79 of file HiEvtPlaneFlatProducer.cc.

Referenced by HiEvtPlaneFlatProducer(), and produce().

centralityMC_

std::string HiEvtPlaneFlatProducer::centralityMC_

private

Definition at line 59 of file HiEvtPlaneFlatProducer.cc.

Referenced by HiEvtPlaneFlatProducer().

centralityTag_

edm::InputTag HiEvtPlaneFlatProducer::centralityTag_

private

Definition at line 64 of file HiEvtPlaneFlatProducer.cc.

Referenced by HiEvtPlaneFlatProducer().

centralityToken_

edm::EDGetTokenT<reco::Centrality> HiEvtPlaneFlatProducer::centralityToken_

private

Definition at line 65 of file HiEvtPlaneFlatProducer.cc.

Referenced by HiEvtPlaneFlatProducer().

centralityVariable_

std::string HiEvtPlaneFlatProducer::centralityVariable_

private

Definition at line 58 of file HiEvtPlaneFlatProducer.cc.

Referenced by HiEvtPlaneFlatProducer().

flat

HiEvtPlaneFlatten* HiEvtPlaneFlatProducer::flat[NumEPNames]

private

Definition at line 90 of file HiEvtPlaneFlatProducer.cc.

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

flatdelvtx_

double HiEvtPlaneFlatProducer::flatdelvtx_

private

Definition at line 86 of file HiEvtPlaneFlatProducer.cc.

Referenced by HiEvtPlaneFlatProducer().

flatminvtx_

double HiEvtPlaneFlatProducer::flatminvtx_

private

Definition at line 85 of file HiEvtPlaneFlatProducer.cc.

Referenced by HiEvtPlaneFlatProducer().

flatnvtxbins_

int HiEvtPlaneFlatProducer::flatnvtxbins_

private

Definition at line 84 of file HiEvtPlaneFlatProducer.cc.

Referenced by HiEvtPlaneFlatProducer().

FlatOrder_

const int HiEvtPlaneFlatProducer::FlatOrder_

private

Definition at line 82 of file HiEvtPlaneFlatProducer.cc.

Referenced by HiEvtPlaneFlatProducer().

flatparmsToken_

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

private

Definition at line 80 of file HiEvtPlaneFlatProducer.cc.

Referenced by HiEvtPlaneFlatProducer(), and produce().

hirpWatcher

edm::ESWatcher<HeavyIonRPRcd> HiEvtPlaneFlatProducer::hirpWatcher

private

Definition at line 78 of file HiEvtPlaneFlatProducer.cc.

Referenced by produce().

hiWatcher

edm::ESWatcher<HeavyIonRcd> HiEvtPlaneFlatProducer::hiWatcher

private

Definition at line 77 of file HiEvtPlaneFlatProducer.cc.

Referenced by produce().

inputPlanesTag_

edm::InputTag HiEvtPlaneFlatProducer::inputPlanesTag_

private

Definition at line 70 of file HiEvtPlaneFlatProducer.cc.

Referenced by HiEvtPlaneFlatProducer().

inputPlanesToken_

edm::EDGetTokenT<reco::EvtPlaneCollection> HiEvtPlaneFlatProducer::inputPlanesToken_

private

Definition at line 71 of file HiEvtPlaneFlatProducer.cc.

Referenced by HiEvtPlaneFlatProducer(), and produce().

nCentBins_

double HiEvtPlaneFlatProducer::nCentBins_

private

Definition at line 92 of file HiEvtPlaneFlatProducer.cc.

Referenced by HiEvtPlaneFlatProducer(), and produce().

NumFlatBins_

int HiEvtPlaneFlatProducer::NumFlatBins_

private

Definition at line 83 of file HiEvtPlaneFlatProducer.cc.

Referenced by HiEvtPlaneFlatProducer().

trackCollection_

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

private

Definition at line 75 of file HiEvtPlaneFlatProducer.cc.

trackTag_

edm::InputTag HiEvtPlaneFlatProducer::trackTag_

private

Definition at line 73 of file HiEvtPlaneFlatProducer.cc.

Referenced by HiEvtPlaneFlatProducer().

trackToken_

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

private

Definition at line 74 of file HiEvtPlaneFlatProducer.cc.

Referenced by HiEvtPlaneFlatProducer().

useOffsetPsi_

bool HiEvtPlaneFlatProducer::useOffsetPsi_

private

Definition at line 91 of file HiEvtPlaneFlatProducer.cc.

Referenced by produce().

vertexTag_

edm::InputTag HiEvtPlaneFlatProducer::vertexTag_

private

Definition at line 67 of file HiEvtPlaneFlatProducer.cc.

Referenced by HiEvtPlaneFlatProducer().

vertexToken_

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

private

Definition at line 68 of file HiEvtPlaneFlatProducer.cc.

Referenced by HiEvtPlaneFlatProducer(), and produce().