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 ( 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 ( )

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

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

overrideprivate

Definition at line 153 of file HiEvtPlaneFlatProducer.cc.

References reco::EvtPlane::addLevel(), newFWLiteAna::bin, c, caloCentRef_, caloCentRefWidth_, CentBinCompression_, centralityBinToken_, centralityESToken_, edm::ESWatcher< T >::check(), HiEvtPlaneFlatten::coffset(), EcalCondDB::db, hi::EPDet, flat, flatparmsToken_, edm::Event::get(), edm::EventSetup::getData(), HiEvtPlaneFlatten::getFlatPsi(), DigiToRawDM_cff::HF, hirpWatcher, hiWatcher, mps_fire::i, inputPlanesToken_, visualization-live-secondInstance_cfg::m, eostools::move(), nCentBins_, hi::NumEPNames, HiEvtPlaneFlatten::offsetPsi(), edm::Event::put(), dt_dqm_sourceclient_common_cff::reco, alignCSCRings::s, HiEvtPlaneFlatten::setCaloCentRefBins(), HiEvtPlaneFlatten::soffset(), parallelization::uint(), useOffsetPsi_, and vertexToken_.

                                                                                   {
   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];
 }