Inheritance diagram for HiEvtPlaneFlatProducer:

Public Member Functions
	HiEvtPlaneFlatProducer (const edm::ParameterSet &)

	~HiEvtPlaneFlatProducer () override

Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () 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 [NumEPNames]

const int	FlatOrder_

edm::ESWatcher< HeavyIonRPRcd >	hirpWatcher

edm::ESWatcher< HeavyIonRcd >	hiWatcher

edm::InputTag	inputPlanesTag_

edm::EDGetTokenT< reco::EvtPlaneCollection >	inputPlanesToken

double	nCentBins_

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 63 of file HiEvtPlaneFlatProducer.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 122 of file HiEvtPlaneFlatProducer.cc.

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, 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")),
   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"))
 {
 //  UseEtHF = kFALSE;
   nCentBins_ = 200.;
 
   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<NumEPNames; i++) {
     flat[i] = new HiEvtPlaneFlatten();
     flat[i]->init(FlatOrder_,NumFlatBins_,EPNames[i],EPOrder[i]);
   }
 
 }

HiEvtPlaneFlatProducer::~HiEvtPlaneFlatProducer ( )

override

Definition at line 167 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 185 of file HiEvtPlaneFlatProducer.cc.

References reco::EvtPlane::addLevel(), reco::EvtPlane::angle(), stringResolutionProvider_cfi::bin, EnergyCorrector::c, caloCentRef_, caloCentRefWidth_, CentBinCompression_, centralityBinToken, centralityLabel_, centralityToken, edm::ESWatcher< T >::check(), DEFINE_FWK_MODULE, hi::EPDet, flat, edm::EventSetup::get(), edm::Event::getByToken(), HiEvtPlaneFlatten::getCoffset(), HiEvtPlaneFlatten::getFlatPsi(), HiEvtPlaneFlatten::getOffsetPsi(), HiEvtPlaneFlatten::getSoffset(), CaloTPGTranscoder_cfi::HF, hirpWatcher, hiWatcher, mps_fire::i, inputPlanesToken, LoadEPDB::IsSuccess(), edm::HandleBase::isValid(), funct::m, CentralityTable::m_table, eostools::move(), nCentBins_, Noffmax_, Noffmin_, reco::Centrality::Ntracks(), hi::NumEPNames, or, edm::Handle< T >::product(), edm::Event::put(), mathSSE::return(), alignCSCRings::s, HiEvtPlaneFlatten::setCaloCentRefBins(), parallelization::uint(), useOffsetPsi_, and vertexToken.

 {
   using namespace edm;
   using namespace std;
   using namespace reco;
 
   //
   //Get Flattening Parameters
   //
   if( hiWatcher.check(iSetup) ) {
 
     //
     //Get Size of Centrality Table
     //
     edm::ESHandle<CentralityTable> centDB_;
     iSetup.get<HeavyIonRcd>().get(centralityLabel_,centDB_);
     nCentBins_ = centDB_->m_table.size();
     for(int i = 0; i<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 = 0;
   edm::Handle<int> cbin_;
   iEvent.getByToken(centralityBinToken, cbin_);
   int cbin = *cbin_;
   bin = cbin/CentBinCompression_;
 
   if(Noffmin_>=0) {
     edm::Handle<reco::Centrality> centrality_;
     iEvent.getByToken(centralityToken, centrality_);
     int Noff = centrality_->Ntracks();
     if ( (Noff < Noffmin_) or (Noff >= Noffmax_) ) {
       return;
     }
   }
   //
   //Get Vertex
   //
   int vs_sell;   // vertex collection size
   float vzr_sell;
   edm::Handle<std::vector<reco::Vertex>> vertex_;
   iEvent.getByToken(vertexToken,vertex_);
   const reco::VertexCollection * vertices3 = vertex_.product();
   vs_sell = vertices3->size();
   if(vs_sell>0) {
     vzr_sell = vertices3->begin()->z();
   } else
     vzr_sell = -999.9;
 
   //
   //Get Event Planes
   //
   
   edm::Handle<reco::EvtPlaneCollection> evtPlanes_;
   iEvent.getByToken(inputPlanesToken,evtPlanes_);
   
   if(!evtPlanes_.isValid()){
     //    cout << "Error! Can't get hiEvtPlane product!" << endl;
     return ;
   }
 
   auto evtplaneOutput = std::make_unique<EvtPlaneCollection>();
   EvtPlane * ep[NumEPNames];
   for(int i = 0; i<NumEPNames; i++) {
     ep[i]=nullptr;
   }
   int indx = 0;
   for (EvtPlaneCollection::const_iterator rp = evtPlanes_->begin();rp !=evtPlanes_->end(); rp++) {
     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< 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];
 }