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
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndRuns () const final

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
std::string	centralityLabel_
std::string	centralityMC_
edm::InputTag	centralityTag_
edm::EDGetTokenT< reco::Centrality >	centralityToken
std::string	centralityVariable_
HiEvtPlaneFlatten *	flat [hi::NumEPNames]
const int	FlatOrder_
edm::ESWatcher< HeavyIonRPRcd >	hirpWatcher
edm::ESWatcher< HeavyIonRcd >	hiWatcher
edm::InputTag	inputPlanesTag_
edm::EDGetTokenT< reco::EvtPlaneCollection >	inputPlanesToken
int	Noffmax_
int	Noffmin_
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<>
typedef CacheContexts< T... >	CacheTypes
typedef CacheTypes::GlobalCache	GlobalCache
typedef AbilityChecker< T... >	HasAbility
typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache
typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext
typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache
typedef CacheTypes::RunCache	RunCache
typedef RunContextT< RunCache, GlobalCache >	RunContext
typedef CacheTypes::RunSummaryCache	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 92 of file HiEvtPlaneFlatProducer.cc.

    : 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")),
      caloCentRef_(iConfig.getParameter<double>("caloCentRef")),
      caloCentRefWidth_(iConfig.getParameter<double>("caloCentRefWidth")),
      CentBinCompression_(iConfig.getParameter<int>("CentBinCompression")),
      Noffmin_(iConfig.getParameter<int>("Noffmin")),
      Noffmax_(iConfig.getParameter<int>("Noffmax")),
      useOffsetPsi_(iConfig.getParameter<bool>("useOffsetPsi")) {
  if (iConfig.exists("nonDefaultGlauberModel")) {
    centralityMC_ = iConfig.getParameter<std::string>("nonDefaultGlauberModel");
  }
  centralityLabel_ = centralityVariable_ + centralityMC_;
 
  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 < hi::NumEPNames; i++) {
    flat[i] = new HiEvtPlaneFlatten();
    flat[i]->init(FlatOrder_, NumFlatBins_, hi::EPNames[i], hi::EPOrder[i]);
  }
}

References centralityBinTag_, centralityBinToken, centralityLabel_, centralityMC_, centralityTag_, centralityToken, centralityVariable_, hi::EPNames, hi::EPOrder, edm::ParameterSet::exists(), flat, FlatOrder_, edm::ParameterSet::getParameter(), mps_fire::i, HiEvtPlaneFlatten::init(), inputPlanesTag_, inputPlanesToken, hi::NumEPNames, NumFlatBins_, AlCaHLTBitMon_QueryRunRegistry::string, trackTag_, trackToken, vertexTag_, and vertexToken.

~HiEvtPlaneFlatProducer()

HiEvtPlaneFlatProducer::~HiEvtPlaneFlatProducer ( )

override

Definition at line 131 of file HiEvtPlaneFlatProducer.cc.

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

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

Member Function Documentation

produce()

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

overrideprivate

Definition at line 140 of file HiEvtPlaneFlatProducer.cc.

                                                                                  {
  using namespace edm;
  using namespace std;
  using namespace reco;
  using namespace hi;
 
  //
  //Get Flattening Parameters
  //
  if (hiWatcher.check(iSetup)) {
    //
    //Get Size of Centrality Table
    //
    edm::ESHandle<CentralityTable> centDB_;
    iSetup.get<HeavyIonRcd>().get(centralityLabel_, centDB_);
    int nCentBins = centDB_->m_table.size();
    for (int i = 0; i < hi::NumEPNames; i++) {
      flat[i]->setCaloCentRefBins(-1, -1);
      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);
        }
      }
    }
  }
 
  if (hirpWatcher.check(iSetup)) {
    edm::ESHandle<RPFlatParams> flatparmsDB_;
    iSetup.get<HeavyIonRPRcd>().get(flatparmsDB_);
    LoadEPDB db(flatparmsDB_, flat);
    if (!db.IsSuccess())
      return;
  }
  //
  //Get Centrality
  //
 
  int bin = iEvent.get(centralityBinToken) / CentBinCompression_;
 
  if (Noffmin_ >= 0) {
    int nOff = iEvent.get(centralityToken).Ntracks();
    if ((nOff < Noffmin_) or (nOff >= Noffmax_)) {
      return;
    }
  }
  //
  //Get Vertex
  //
  int vs_sell;  // vertex collection size
  float vzr_sell;
  auto const& vertices3 = iEvent.get(vertexToken);
  vs_sell = vertices3.size();
  if (vs_sell > 0) {
    vzr_sell = vertices3.begin()->z();
  } else
    vzr_sell = -999.9;
 
  //
  //Get Event Planes
  //
 
  auto evtplaneOutput = std::make_unique<EvtPlaneCollection>();
  EvtPlane* ep[hi::NumEPNames];
  for (int i = 0; i < hi::NumEPNames; i++) {
    ep[i] = nullptr;
  }
  int indx = 0;
  for (auto const& rp : iEvent.get(inputPlanesToken)) {
    double psiOffset = rp.angle(0);
    double s = rp.sumSin(0);
    double c = rp.sumCos(0);
    uint m = rp.mult();
 
    double soff = s;
    double coff = c;
    if (useOffsetPsi_) {
      soff = flat[indx]->getSoffset(s, vzr_sell, bin);
      coff = flat[indx]->getCoffset(c, vzr_sell, bin);
      psiOffset = flat[indx]->getOffsetPsi(soff, coff);
    }
    double psiFlat = flat[indx]->getFlatPsi(psiOffset, vzr_sell, 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 < hi::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];
}

References newFWLiteAna::bin, HltBtagPostValidation_cff::c, caloCentRef_, caloCentRefWidth_, CentBinCompression_, centralityBinToken, centralityLabel_, centralityToken, edm::ESWatcher< T >::check(), dqmiodatasetharvest::db, SiStripBadComponentsDQMServiceTemplate_cfg::ep, hi::EPDet, flat, edm::EventSetup::get(), get, HiEvtPlaneFlatten::getCoffset(), HiEvtPlaneFlatten::getFlatPsi(), HiEvtPlaneFlatten::getOffsetPsi(), HiEvtPlaneFlatten::getSoffset(), CaloTPGTranscoder_cfi::HF, hirpWatcher, hiWatcher, mps_fire::i, iEvent, inputPlanesToken, visualization-live-secondInstance_cfg::m, CentralityTable::m_table, eostools::move(), Noffmax_, Noffmin_, hi::NumEPNames, or, alignCSCRings::s, HiEvtPlaneFlatten::setCaloCentRefBins(), parallelization::uint, useOffsetPsi_, and vertexToken.

Member Data Documentation

caloCentRef_

double HiEvtPlaneFlatProducer::caloCentRef_

private

Definition at line 83 of file HiEvtPlaneFlatProducer.cc.

Referenced by produce().

caloCentRefWidth_

double HiEvtPlaneFlatProducer::caloCentRefWidth_

private

Definition at line 84 of file HiEvtPlaneFlatProducer.cc.

Referenced by produce().

CentBinCompression_

int HiEvtPlaneFlatProducer::CentBinCompression_

private

Definition at line 85 of file HiEvtPlaneFlatProducer.cc.

Referenced by produce().

centralityBinTag_

edm::InputTag HiEvtPlaneFlatProducer::centralityBinTag_

private

Definition at line 62 of file HiEvtPlaneFlatProducer.cc.

Referenced by HiEvtPlaneFlatProducer().

centralityBinToken

edm::EDGetTokenT<int> HiEvtPlaneFlatProducer::centralityBinToken

private

Definition at line 63 of file HiEvtPlaneFlatProducer.cc.

Referenced by HiEvtPlaneFlatProducer(), and produce().

centralityLabel_

std::string HiEvtPlaneFlatProducer::centralityLabel_

private

Definition at line 59 of file HiEvtPlaneFlatProducer.cc.

Referenced by HiEvtPlaneFlatProducer(), and produce().

centralityMC_

std::string HiEvtPlaneFlatProducer::centralityMC_

private

Definition at line 60 of file HiEvtPlaneFlatProducer.cc.

Referenced by HiEvtPlaneFlatProducer().

centralityTag_

edm::InputTag HiEvtPlaneFlatProducer::centralityTag_

private

Definition at line 65 of file HiEvtPlaneFlatProducer.cc.

Referenced by HiEvtPlaneFlatProducer().

centralityToken

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

private

Definition at line 66 of file HiEvtPlaneFlatProducer.cc.

Referenced by HiEvtPlaneFlatProducer(), and produce().

centralityVariable_

std::string HiEvtPlaneFlatProducer::centralityVariable_

private

Definition at line 58 of file HiEvtPlaneFlatProducer.cc.

Referenced by HiEvtPlaneFlatProducer().

flat

HiEvtPlaneFlatten* HiEvtPlaneFlatProducer::flat[hi::NumEPNames]

private

Definition at line 88 of file HiEvtPlaneFlatProducer.cc.

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

FlatOrder_

const int HiEvtPlaneFlatProducer::FlatOrder_

private

Definition at line 81 of file HiEvtPlaneFlatProducer.cc.

Referenced by HiEvtPlaneFlatProducer().

hirpWatcher

edm::ESWatcher<HeavyIonRPRcd> HiEvtPlaneFlatProducer::hirpWatcher

private

Definition at line 79 of file HiEvtPlaneFlatProducer.cc.

Referenced by produce().

hiWatcher

edm::ESWatcher<HeavyIonRcd> HiEvtPlaneFlatProducer::hiWatcher

private

Definition at line 78 of file HiEvtPlaneFlatProducer.cc.

Referenced by produce().

inputPlanesTag_

edm::InputTag HiEvtPlaneFlatProducer::inputPlanesTag_

private

Definition at line 71 of file HiEvtPlaneFlatProducer.cc.

Referenced by HiEvtPlaneFlatProducer().

inputPlanesToken

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

private

Definition at line 72 of file HiEvtPlaneFlatProducer.cc.

Referenced by HiEvtPlaneFlatProducer(), and produce().

Noffmax_

int HiEvtPlaneFlatProducer::Noffmax_

private

Definition at line 87 of file HiEvtPlaneFlatProducer.cc.

Referenced by produce().

Noffmin_

int HiEvtPlaneFlatProducer::Noffmin_

private

Definition at line 86 of file HiEvtPlaneFlatProducer.cc.

Referenced by produce().

NumFlatBins_

int HiEvtPlaneFlatProducer::NumFlatBins_

private

Definition at line 82 of file HiEvtPlaneFlatProducer.cc.

Referenced by HiEvtPlaneFlatProducer().

trackCollection_

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

private

Definition at line 76 of file HiEvtPlaneFlatProducer.cc.

trackTag_

edm::InputTag HiEvtPlaneFlatProducer::trackTag_

private

Definition at line 74 of file HiEvtPlaneFlatProducer.cc.

Referenced by HiEvtPlaneFlatProducer().

trackToken

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

private

Definition at line 75 of file HiEvtPlaneFlatProducer.cc.

Referenced by HiEvtPlaneFlatProducer().

useOffsetPsi_

bool HiEvtPlaneFlatProducer::useOffsetPsi_

private

Definition at line 89 of file HiEvtPlaneFlatProducer.cc.

Referenced by produce().

vertexTag_

edm::InputTag HiEvtPlaneFlatProducer::vertexTag_

private

Definition at line 68 of file HiEvtPlaneFlatProducer.cc.

Referenced by HiEvtPlaneFlatProducer().

vertexToken

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

private

Definition at line 69 of file HiEvtPlaneFlatProducer.cc.

Referenced by HiEvtPlaneFlatProducer(), and produce().