L1CaloInputScalesGenerator Class Reference

#include <L1TriggerConfig/L1ScalesProducers/src/L1CaloInputScalesGenerator.cc>

Inheritance diagram for L1CaloInputScalesGenerator:

Public Member Functions
	L1CaloInputScalesGenerator (const edm::ParameterSet &)
	~L1CaloInputScalesGenerator () override
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
	~EDAnalyzer () override
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...
	EDConsumerBase ()
	EDConsumerBase (EDConsumerBase const &)=delete
	EDConsumerBase (EDConsumerBase &&)=default
ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
EDConsumerBase &	operator= (EDConsumerBase &&)=default
bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override
void	beginJob () override
void	endJob () override

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &)
static bool	wantsGlobalLuminosityBlocks ()
static bool	wantsGlobalRuns ()
static bool	wantsStreamLuminosityBlocks ()
static bool	wantsStreamRuns ()
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

Description: <one line="" class="" summary>="">

Implementation: <Notes on="" implementation>="">

Definition at line 27 of file L1CaloInputScalesGenerator.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 40 of file L1CaloInputScalesGenerator.cc.

{
   //now do what ever initialization is needed

}

L1CaloInputScalesGenerator::~L1CaloInputScalesGenerator ( )

override

Definition at line 48 of file L1CaloInputScalesGenerator.cc.

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

}

Member Function Documentation

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

overrideprivate

Definition at line 63 of file L1CaloInputScalesGenerator.cc.

References EcalBarrel, EcalEndcap, edm::EventSetup::get(), EcalTPGScale::getTPGInGeV(), CaloTPGTranscoder::hcaletValue(), input, convertSQLitetoXML_cfg::output, and EcalTPGScale::setEventSetup().

{
   using namespace edm;


   ESHandle<CaloTPGTranscoder> caloTPGTranscoder;
   iSetup.get<CaloTPGRecord>().get(caloTPGTranscoder);
   
   EcalTPGScale* ecalTPGScale = new EcalTPGScale();
   ecalTPGScale->setEventSetup(iSetup);

   double output; 
   ofstream scalesFile("L1CaloInputScales_cfi.py");


   // Write the ecal scales, positive eta

   scalesFile << "import FWCore.ParameterSet.Config as cms\n" <<endl;

   scalesFile << "L1CaloInputScalesProducer =cms.ESProducer(\"L1CaloInputScalesProducer\"," << endl;
   scalesFile << "L1EcalEtThresholdsPositiveEta = cms.vdouble(" << endl; 


   //Python does not support arrays over 255 entries so we neeed ton accomodate it by creating new array after 255 entries
   int nEntries = 0;



   // loop over ietas, barrel
   for (unsigned short absIeta = 1; absIeta <= 28; absIeta++)
     {
       EcalSubdetector subdet = ( absIeta <= 17 ) ? EcalBarrel : EcalEndcap ;
       // 8 bits of input energy
       for (unsigned short input = 0; input <= 0xFF; input++)
     {
       output = ecalTPGScale->getTPGInGeV( (unsigned int) input, 
                          EcalTrigTowerDetId(1, subdet,
                                 absIeta, 1));
       scalesFile << setprecision (8) << output;
       nEntries++;

       if (absIeta == 28 && input == 0xFF)
         {
           scalesFile << "),";
         }
       else if(nEntries>254)
         {
           scalesFile <<")+cms.vdouble(";
           nEntries=0;
         }
       else
         {
           scalesFile << ", ";
         }
     }
       scalesFile << endl;
     }

   // Write the ecal scales, negative eta

   scalesFile << endl << "\tL1EcalEtThresholdsNegativeEta = cms.vdouble(" << endl;

   nEntries=0;
   // loop over ietas, barrel first
   for (unsigned short absIeta = 1; absIeta <= 28; absIeta++)
     {
       EcalSubdetector subdet = ( absIeta <= 17 ) ? EcalBarrel : EcalEndcap ;
       // 8 bits of input energy
       for (unsigned short input = 0; input <= 0xFF; input++)
     {
       // negative eta
       output = ecalTPGScale->
         getTPGInGeV( (unsigned int) input, EcalTrigTowerDetId(-1, subdet,
                               absIeta, 2));
       scalesFile << setprecision (8) << output;
       nEntries++;



       if (absIeta == 28 && input == 0xFF)
         {
           scalesFile << "),";
         }
       else if(nEntries>254)
         {
           scalesFile <<")+cms.vdouble(";
           nEntries=0;
         }
       else
         {
           scalesFile << ", ";
         }
     }
       scalesFile << endl;
     }

   // Write the hcal scales (Positive Eta)

     scalesFile << endl << "\tL1HcalEtThresholdsPositiveEta = cms.vdouble(" << endl;

   // loop over ietas

     nEntries=0;
   for (unsigned short absIeta = 1; absIeta <= 32; absIeta++)
     {
       for (unsigned short input = 0; input <= 0xFF; input++)
     {
       output = caloTPGTranscoder->hcaletValue(absIeta, input); 
       scalesFile << setprecision (8) << output ;
       nEntries++;


       if (absIeta == 32 && input == 0xFF)
         {
           scalesFile << "),";
         }
       else if(nEntries>254)
         {
           scalesFile <<")+cms.vdouble(";
           nEntries=0;
         }
       else
         {
           scalesFile << ", ";
         }
     }
       scalesFile << endl;
     }

   // Write the hcal scales (Negative Eta)

     scalesFile << endl << "\tL1HcalEtThresholdsNegativeEta = cms.vdouble(" << endl;

     nEntries=0;  
   // loop over ietas
   for (unsigned short absIeta = 1; absIeta <= 32; absIeta++)
     {
       for (unsigned short input = 0; input <= 0xFF; input++)
     {
       output = caloTPGTranscoder->hcaletValue(-absIeta, input); 
       scalesFile << setprecision (8) << output ;
       nEntries++;


       if (absIeta == 32 && input == 0xFF)
         {
           scalesFile << ")";
         }
       else if(nEntries>254)
         {
           scalesFile <<")+cms.vdouble(";
           nEntries=0;
         }
       else 
         {
           scalesFile << ", ";
         }
     }
       scalesFile << endl;
     }


   scalesFile << ")" << endl;

   scalesFile.close();
}

void L1CaloInputScalesGenerator::beginJob ( void  )

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 233 of file L1CaloInputScalesGenerator.cc.

{
}

void L1CaloInputScalesGenerator::endJob ( void  )

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 239 of file L1CaloInputScalesGenerator.cc.

                                   {
}