Inheritance diagram for UETableProducer:

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

	~UETableProducer ()

Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzer ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

virtual	~EDAnalyzer ()

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

	EDConsumerBase ()

ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Private Member Functions
virtual void	analyze (const edm::Event &, const edm::EventSetup &)

virtual void	beginRun (const edm::EventSetup &)

virtual void	endJob ()

Private Attributes
string	calibrationFile_

bool	debug_

unsigned int	index = 0

bool	jetCorrectorFormat_

unsigned int	ni0 [2]

unsigned int	ni1 [2]

unsigned int	ni2 [2]

unsigned int	np [5]

unsigned int	runnum_

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType

Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Definition at line 35 of file UETableProducer.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 76 of file UETableProducer.cc.

References calibrationFile_, debug_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), jetCorrectorFormat_, ni0, ni1, ni2, np, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                               :
   runnum_(0)
 {
   //now do what ever initialization is needed
   calibrationFile_ = iConfig.getParameter<std::string>("txtFile");
   //calibrationFile_ = "RecoHI/HiJetAlgos/data/ue_calibrations_pf_data.txt";
 
   debug_ = iConfig.getUntrackedParameter<bool>("debug",false);
   jetCorrectorFormat_ = iConfig.getUntrackedParameter<bool>("jetCorrectorFormat",false);
 
   np[0] = 3;// Number of reduced PF ID (track, ECAL, HCAL)
   np[1] = 15;// Number of pseudorapidity block
   np[2] = 5;// Fourier series order
   np[3] = 2;// Re or Im
   np[4] = 82;// Number of feature parameter
 
   ni0[0] = np[1];
   ni0[1] = 344;// Number of track binning (kept = ECAL)
 
   ni1[0] = np[1];
   ni1[1] = 344;// Number of ECAL binning
 
   ni2[0] = np[1];
   ni2[1] = 82;// Number of HCAL binning
 
 }

UETableProducer::~UETableProducer ( )

Definition at line 103 of file UETableProducer.cc.

 {
   // do anything here that needs to be done at desctruction time
   // (e.g. close files, deallocate resources etc.)
 }

Member Function Documentation

void UETableProducer::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

privatevirtual

Implements edm::EDAnalyzer.

Definition at line 115 of file UETableProducer.cc.

 {
   // nothing
 
 }

void UETableProducer::beginRun ( const edm::EventSetup & iSetup )

privatevirtual

Definition at line 123 of file UETableProducer.cc.

124 {

125 }

void UETableProducer::endJob ( void )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 129 of file UETableProducer.cc.

References cond::service::PoolDBOutputService::appendSinceTime(), cond::service::PoolDBOutputService::beginOfTime(), calibrationFile_, filterCSVwithJSON::copy, gather_cfg::cout, cond::service::PoolDBOutputService::createNewIOV(), cond::service::PoolDBOutputService::currentTime(), UETable::edgeEta, cond::service::PoolDBOutputService::endOfTime(), index, edm::Service< T >::isAvailable(), cond::service::PoolDBOutputService::isNewTagRequest(), jetCorrectorFormat_, JetCorrectorParametersCollection::L1Offset, geometryCSVtoXML::line, UETable::ni0, ni0, UETable::ni1, ni1, UETable::ni2, ni2, UETable::np, np, NULL, fetchall_from_DQM_v2::pool, JetCorrectorParametersCollection::push_back(), record, AlCaHLTBitMon_QueryRunRegistry::string, and UETable::values.

                         {
   std::string qpDataName = calibrationFile_.c_str();
   std::ifstream textTable_(qpDataName.c_str());
 
   std::vector<float> ue_vec;
   UETable *ue_predictor_pf = NULL;
 
   if (!jetCorrectorFormat_) {
     ue_predictor_pf = new UETable();
   }
   //  unsigned int Nnp_full = np[0] * np[1] * np[2] * np[3] * np[4];
   unsigned int Nnp = np[0] * np[1] * (1 + (np[2] - 1) * np[3]) * np[4];
   unsigned int Nni0 = ni0[0]*ni0[1];
   unsigned int Nni1 = ni1[0]*ni1[1];
   unsigned int Nni2 = ni2[0]*ni2[1];
 
   if (!jetCorrectorFormat_) {
   ue_predictor_pf->np.resize(5);
   ue_predictor_pf->ni0.resize(2);
   ue_predictor_pf->ni1.resize(2);
   ue_predictor_pf->ni2.resize(2);
 
   std::copy(np, np + 5, ue_predictor_pf->np.begin());
   std::copy(ni0, ni0 + 2, ue_predictor_pf->ni0.begin());
   std::copy(ni1, ni1 + 2, ue_predictor_pf->ni1.begin());
   std::copy(ni2, ni2 + 2, ue_predictor_pf->ni2.begin());
 
   static const float edge_pseudorapidity[16] = {
     -5.191, -2.650, -2.043, -1.740, -1.479, -1.131, -0.783, -0.522, 0.522, 0.783, 1.131, 1.479, 1.740, 2.043, 2.650, 5.191
   };
 
   ue_predictor_pf->edgeEta.resize(16);
 
   std::copy(edge_pseudorapidity, edge_pseudorapidity + 16, ue_predictor_pf->edgeEta.begin());
 
   }
 
   std::string line;
 
   while( std::getline( textTable_, line)){
     if(!line.size() || line[0]=='#') {
       std::cout<<" continue "<<std::endl;
       continue;
     }
     std::istringstream linestream(line);
     float val;
     int bin0, bin1, bin2, bin3, bin4;
     if(index < Nnp){
       linestream>>bin0>>bin1>>bin2>>bin3>>bin4>>val;
       ue_vec.push_back(val);
     }else if(index < Nnp + Nni0){
       linestream>>bin0>>bin1>>val;
       ue_vec.push_back(val);
     }else if(index < Nnp + Nni0 + Nni1){
       linestream>>bin0>>bin1>>val;
       ue_vec.push_back(val);
     }else if(index < Nnp + Nni0 + Nni1 + Nni2){
       linestream>>bin0>>bin1>>val;
       ue_vec.push_back(val);
     }
     ++index;
 
   }
 
   edm::Service<cond::service::PoolDBOutputService> pool;
 
   if( pool.isAvailable() ){
 
   if (jetCorrectorFormat_) {
     // A minimal dummy line that satisfies the JME # token >= 6 requirement, and has the correct key type
     JetCorrectorParameters::Definitions definition("1 0 0 0 Correction L1Offset");
     std::vector<JetCorrectorParameters::Record> record(1, JetCorrectorParameters::Record(ue_vec.size(), std::vector<float>(ue_vec.size(), 0), std::vector<float>(ue_vec.size(), 0), ue_vec));
     JetCorrectorParameters parameter(definition, record);
 
     JetCorrectorParametersCollection *jme_payload = new JetCorrectorParametersCollection();
 
     jme_payload->push_back(JetCorrectorParametersCollection::L1Offset, parameter);
 
     if( pool->isNewTagRequest( "HeavyIonUERcd" ) ){
       pool->createNewIOV<JetCorrectorParametersCollection>( jme_payload, pool->beginOfTime(), pool->endOfTime(), "HeavyIonUERcd" );
     }else{
       pool->appendSinceTime<JetCorrectorParametersCollection>( jme_payload, pool->currentTime(), "HeavyIonUERcd" );
     }
   }
   else {
     ue_predictor_pf->values = ue_vec;
 
     if( pool->isNewTagRequest( "HeavyIonUERcd" ) ){
       pool->createNewIOV<UETable>( ue_predictor_pf, pool->beginOfTime(), pool->endOfTime(), "HeavyIonUERcd" );
     }else{
       pool->appendSinceTime<UETable>( ue_predictor_pf, pool->currentTime(), "HeavyIonUERcd" );
     }
   }
   }
 }