#include <L1RCTLutWriter.h>

Inheritance diagram for L1RCTLutWriter:

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

	~L1RCTLutWriter ()

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

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

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

virtual	~EDConsumerBase () noexcept(false)

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

virtual void	endJob ()

void	writeEicLutFile (unsigned short card)

void	writeJscLutFile ()

void	writeRcLutFile (unsigned short card)

void	writeThresholdsFile (unsigned int eicThreshold, unsigned int jscThresholdBarrel, unsigned int jscThresholdEndcap)

Private Attributes
std::string	keyName_

L1RCTLookupTables *	lookupTable_

std::ofstream	lutFile_

const L1RCTParameters *	rctParameters_

bool	useDebugTpgScales_

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

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

L1Trigger/L1RCTLutWriter/src/L1RCTLutWriter.cc

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

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

Definition at line 48 of file L1RCTLutWriter.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 37 of file L1RCTLutWriter.cc.

                                                              :
   lookupTable_(new L1RCTLookupTables),
   keyName_(iConfig.getParameter<std::string>("key")),
   useDebugTpgScales_(iConfig.getParameter<bool>("useDebugTpgScales"))
 {
    //now do what ever initialization is needed
  
 }

L1RCTLutWriter::~L1RCTLutWriter ( )

Definition at line 47 of file L1RCTLutWriter.cc.

References lookupTable_.

 {
  
    // do anything here that needs to be done at destruction time
    // (e.g. close files, deallocate resources etc.)
   
   if (lookupTable_ != 0) delete lookupTable_;
 }

Member Function Documentation

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

privatevirtual

Definition at line 63 of file L1RCTLutWriter.cc.

References gather_cfg::cout, MillePedeFileConverter_cfg::e, EcalBarrel, EcalEndcap, L1RCTNoisyChannelMask::ecalMask, L1RCTChannelMask::ecalMask, L1RCTNoisyChannelMask::ecalThreshold, L1RCTParameters::eicIsolationThreshold(), edm::EventSetup::get(), EcalTPGScale::getTPGInGeV(), CaloTPGTranscoder::hcaletValue(), L1RCTChannelMask::hcalMask, L1RCTNoisyChannelMask::hcalMask, L1RCTNoisyChannelMask::hcalThreshold, L1RCTChannelMask::hfMask, L1RCTNoisyChannelMask::hfMask, L1RCTNoisyChannelMask::hfThreshold, mps_fire::i, L1RCTParameters::jscQuietThresholdBarrel(), L1RCTParameters::jscQuietThresholdEndcap(), gen::k, lookupTable_, funct::m, L1CaloEcalScale::nBinEta, L1CaloHcalScale::nBinEta, L1CaloEcalScale::nBinRank, L1CaloHcalScale::nBinRank, edm::ESHandle< T >::product(), rctParameters_, alignCSCRings::s, L1CaloEcalScale::setBin(), L1CaloHcalScale::setBin(), L1RCTLookupTables::setChannelMask(), L1RCTLookupTables::setEcalScale(), EcalTPGScale::setEventSetup(), L1RCTLookupTables::setHcalScale(), L1RCTLookupTables::setL1CaloEtScale(), L1RCTLookupTables::setNoisyChannelMask(), L1RCTLookupTables::setRCTParameters(), useDebugTpgScales_, writeEicLutFile(), writeJscLutFile(), writeRcLutFile(), and writeThresholdsFile().

 {
    
    // get all the configuration information from the event, set it
    // in the lookuptable
    edm::ESHandle<L1RCTParameters> rctParameters;
    iSetup.get<L1RCTParametersRcd>().get(rctParameters);
    rctParameters_ = rctParameters.product();
    edm::ESHandle<L1CaloEtScale> emScale;
    iSetup.get<L1EmEtScaleRcd>().get(emScale);
    const L1CaloEtScale* s = emScale.product();
 
    // make dummy channel mask -- we don't want to mask
    // any channels when writing LUTs, that comes afterwards
    L1RCTChannelMask* m = new L1RCTChannelMask;
    for (int i = 0; i < 18; i++)
      {
        for (int j = 0; j < 2; j++)
      {
        for (int k = 0; k < 28; k++)
          {
            m->ecalMask[i][j][k] = false;
            m->hcalMask[i][j][k] = false;
          }
        for (int k = 0; k < 4; k++)
          {
            m->hfMask[i][j][k] = false;
          }
      }
      }
 
 
    //Same for Noisy mask
    // make dummy channel mask -- we don't want to mask
    // any channels when writing LUTs, that comes afterwards
    L1RCTNoisyChannelMask* m2 = new L1RCTNoisyChannelMask;
    for (int i = 0; i < 18; i++)
      {
        for (int j = 0; j < 2; j++)
      {
        for (int k = 0; k < 28; k++)
          {
            m2->ecalMask[i][j][k] = false;
            m2->hcalMask[i][j][k] = false;
          }
        for (int k = 0; k < 4; k++)
          {
            m2->hfMask[i][j][k] = false;
          }
      }
      }
    m2->ecalThreshold = 0.0;
    m2->hcalThreshold = 0.0;
    m2->hfThreshold = 0.0;
 
 
    
   // use these dummies to get the delete right when using old-style
   // scales to create set of L1CaloXcalScales
   L1CaloEcalScale* dummyE(0);
   L1CaloHcalScale* dummyH(0);
 
   if (useDebugTpgScales_) // generate new-style scales from tpg scales
     {
 
       std::cout << "Using old-style TPG scales!" << std::endl;
 
       // old version of hcal energy scale to convert input
       edm::ESHandle<CaloTPGTranscoder> transcoder;
       iSetup.get<CaloTPGRecord>().get(transcoder);
       const CaloTPGTranscoder* h_tpg = transcoder.product();
 
       // old version of ecal energy scale to convert input
       EcalTPGScale* e_tpg = new EcalTPGScale();
       e_tpg->setEventSetup(iSetup);
 
       L1CaloEcalScale* ecalScale = new L1CaloEcalScale();
       L1CaloHcalScale* hcalScale = new L1CaloHcalScale();
 
       // generate L1CaloXcalScales from old-style scales (thanks, werner!)
 
       // ECAL
       for( unsigned short ieta = 1 ; ieta <= L1CaloEcalScale::nBinEta; ++ieta )
     {
       for( unsigned short irank = 0 ; irank < L1CaloEcalScale::nBinRank; ++irank )
         {
           EcalSubdetector subdet = ( ieta <= 17 ) ? EcalBarrel : EcalEndcap ;
           double etGeVPos =
         e_tpg->getTPGInGeV
         ( irank, EcalTrigTowerDetId(1, // +ve eta
                         subdet,
                         ieta,
                         1 )); // dummy phi value
           ecalScale->setBin( irank, ieta, 1, etGeVPos ) ;
         }
     }
       
       for( unsigned short ieta = 1 ; ieta <= L1CaloEcalScale::nBinEta; ++ieta )
     {
       for( unsigned short irank = 0 ; irank < L1CaloEcalScale::nBinRank; ++irank )
         {
           EcalSubdetector subdet = ( ieta <= 17 ) ? EcalBarrel : EcalEndcap ;
           
           double etGeVNeg =
         e_tpg->getTPGInGeV
         ( irank,
           EcalTrigTowerDetId(-1, // -ve eta
                      subdet,
                      ieta,
                      2 )); // dummy phi value
           ecalScale->setBin( irank, ieta, -1, etGeVNeg ) ;
         }
     }
 
       // HCAL
       for( unsigned short ieta = 1 ; ieta <= L1CaloHcalScale::nBinEta; ++ieta )
     {
       for( unsigned short irank = 0 ; irank < L1CaloHcalScale::nBinRank; ++irank )
         {
           double etGeV = h_tpg->hcaletValue( ieta, irank ) ;
 
           hcalScale->setBin( irank, ieta, 1, etGeV ) ;
           hcalScale->setBin( irank, ieta, -1, etGeV ) ;
         }
     }
 
       // set the input scales
       lookupTable_->setEcalScale(ecalScale);
       lookupTable_->setHcalScale(hcalScale);
 
       dummyE = ecalScale;
       dummyH = hcalScale;
 
       delete e_tpg;
 
     }
   else
     {
 
       // get energy scale to convert input from ECAL
       edm::ESHandle<L1CaloEcalScale> ecalScale;
       iSetup.get<L1CaloEcalScaleRcd>().get(ecalScale);
       const L1CaloEcalScale* e = ecalScale.product();
       
       // get energy scale to convert input from HCAL
       edm::ESHandle<L1CaloHcalScale> hcalScale;
       iSetup.get<L1CaloHcalScaleRcd>().get(hcalScale);
       const L1CaloHcalScale* h = hcalScale.product();
 
       // set scales
       lookupTable_->setEcalScale(e);
       lookupTable_->setHcalScale(h);
 
     }
 
    lookupTable_->setRCTParameters(rctParameters_);
    lookupTable_->setChannelMask(m);
    lookupTable_->setNoisyChannelMask(m2);
    //lookupTable_->setHcalScale(h);
    //lookupTable_->setEcalScale(e);
    lookupTable_->setL1CaloEtScale(s);
 
    for (unsigned short nCard = 0; nCard <= 6; nCard = nCard + 2)
      {
        writeRcLutFile(nCard);
        writeEicLutFile(nCard);
      }
    writeJscLutFile();
 
    unsigned int eicThreshold = rctParameters_->eicIsolationThreshold();
    unsigned int jscThresholdBarrel = rctParameters_->jscQuietThresholdBarrel();
    unsigned int jscThresholdEndcap = rctParameters_->jscQuietThresholdEndcap();
    writeThresholdsFile(eicThreshold, jscThresholdBarrel, jscThresholdEndcap);
 
   if (dummyE != 0) delete dummyE;
   if (dummyH != 0) delete dummyH;
 
 }

void L1RCTLutWriter::endJob ( void )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 247 of file L1RCTLutWriter.cc.

247 {

248 }

void L1RCTLutWriter::writeEicLutFile ( unsigned short card )

private

Definition at line 366 of file L1RCTLutWriter.cc.

References mps_alisetup::command, data, TauDecayModes::dec, L1RCTLookupTables::emRank(), corrVsCorr::filename, keyName_, lookupTable_, lutFile_, and estimatePileup_makeJSON::trunc.

Referenced by analyze().

 {
   // try timestamp
   char filename[256];
   char command[64];
   if (card != 6)
     {
       int card2 = card + 1;
       sprintf(filename,"EIC%i%i-%s.dat", card, card2, keyName_.c_str() );
     }
   else
     {
       sprintf(filename,"EIC6-%s.dat", keyName_.c_str() );
     }
   // open file for writing, delete any existing content
   lutFile_.open(filename, std::ios::trunc);
   lutFile_ << "Emulator-parameter generated EIC lut file, card " 
        << card << " key " << keyName_ << "   ";
   // close to append timestamp info
   lutFile_.close();
   sprintf(command, "date >> %s", filename);
   system(command);
 
   // reopen file for writing values
   lutFile_.open(filename, std::ios::app);
 
   unsigned long data = 0;
 
   // write all memory addresses in increasing order
   // address = (1<<22) + (etIn7Bits<<1)
 
   // 2^7 = 0x7f = 128
   for (int etIn7Bits = 0; etIn7Bits < 128; etIn7Bits++)
     {
       data = lookupTable_->emRank(etIn7Bits);
       if (data > 0x3f)
     {
       data = 0x3f;
     }
       lutFile_ << std::hex << data << std::dec << std::endl;
     }
   lutFile_.close();
   return;
 }

void L1RCTLutWriter::writeJscLutFile ( )

private

Definition at line 413 of file L1RCTLutWriter.cc.

References mps_alisetup::command, data, TauDecayModes::dec, corrVsCorr::filename, keyName_, L1RCTLookupTables::lookup(), lookupTable_, lutFile_, and estimatePileup_makeJSON::trunc.

Referenced by analyze().

 {
   char filename[256];
   char command[64];
   sprintf(filename, "JSC-%s.dat", keyName_.c_str() );
 
   // open file; if it already existed, delete existing content
   lutFile_.open(filename, std::ios::trunc);
   lutFile_ << "Emulator parameter-generated lut file, key " << keyName_ 
        << "   ";
   // close to append time-stamp
   lutFile_.close();
   sprintf(command, "date >> %s", filename);
   system(command);
   // reopen file for writing
   lutFile_.open(filename, std::ios::app);
 
   unsigned long data = 0;
   unsigned long data0 = 0;
   unsigned long data1 = 0;
 
   // write all memory addresses in increasing order
   // address = (1<<22) + (lutbits<<17) + (phi1et<<9) + (phi0et<<1);
 
   // ecl and U93/U225 lut id bits, identify eta segment of hf
   for (int lutbits = 0; lutbits < 4; lutbits++)
     {
       // 8-bit phi_1 et for each eta partition
       for (unsigned int phi1et = 0; phi1et < 256; phi1et++)
     {
       // 8-bit phi_0 et for each eta
       for (unsigned int phi0et = 0; phi0et < 256; phi0et++)
         {
           // lookup takes "(hf_et, crate, card, tower)"
           // "card" convention for hf is 999, tower is 0-7
           // but equivalent to 0-3 == lutbits
           // crate doesn't matter, take 0
           // only |ieta| matters
           data0 = lookupTable_->lookup(phi0et, 0, 999, lutbits);
           if (data0 > 0xff)
         {
           data0 = 0xff;   // 8-bit output energy for each phi region
         }
           data1 = lookupTable_->lookup(phi1et, 0, 999, lutbits);
           if (data1 > 0xff)
         {
           data1 = 0xff;   // 8-bit output energy for each phi region
         }
           data = (data1<<8) + (data0);
           lutFile_ << std::hex << data << std::dec << std::endl;
           /*
           if (phi0et < 10 && phi1et < 10)
         {
           std::cout << "Writer: jsc. lutbits=" << lutbits 
                 << " phi0et=" << phi0et << " data0=" << data0
                 << " phi1et=" << phi1et << " data1=" << data1
                 << std::endl;
         }
           */
         }
     }
     }
   
   lutFile_.close();
   return;
 }

void L1RCTLutWriter::writeRcLutFile ( unsigned short card )

private

Definition at line 253 of file L1RCTLutWriter.cc.

References L1RCTParameters::calcCrate(), L1RCTParameters::calcTower(), mps_alisetup::command, data, TauDecayModes::dec, corrVsCorr::filename, keyName_, lookupTable_, lutFile_, convertSQLitetoXML_cfg::output, rctParameters_, and estimatePileup_makeJSON::trunc.

Referenced by analyze().

 {
 
   // don't mess yet with name
   char filename[256];
   char command[64];
   if (card != 6)
     {
       int card2 = card + 1;
       sprintf(filename,"RC%i%i-%s.dat",card,card2,keyName_.c_str() );
       //sprintf(filename, "RC%i%i.dat", card, card2);
     }
   else
     {
       sprintf(filename,"RC6-%s.dat",keyName_.c_str() );
       //sprintf(filename, "RC6.dat");
     }
   // open file for writing, delete any existing content
   lutFile_.open(filename, std::ios::trunc);
   lutFile_ << "Emulator-parameter generated lut file, card " 
        << card << " key " << keyName_ << "   ";
 
   // close to append timestamp info
   lutFile_.close();
   sprintf(command, "date >> %s", filename);
   system(command);
 
   // reopen file for writing values
   lutFile_.open(filename, std::ios::app);
 
   unsigned long data = 0;
 
   // write all memory addresses in increasing order
   // address = (1<<22)+(nLUT<<19)+(eG?<18)+(hcalEt<<10)+(ecalfg<<9)+(ecalEt<<1)
 
   // loop through the physical LUTs on the card, 0-7
   for (unsigned short nLUT = 0; nLUT < 8; nLUT++)
     {
       // determine ieta, iphi, etc. everything
       unsigned short iAbsEta = 0;
       if (card != 6)
     {
       iAbsEta = (card/2)*8 + nLUT + 1;
     }
       else
     {
       if (nLUT < 4)
         {
           iAbsEta = (card/2)*8 + nLUT + 1;
         }
       else
         {
           iAbsEta = (card/2)*8 + (3 - (nLUT%4) ) + 1;         
         }
       //std::cout << "LUT is " << nLUT << " iAbsEta is " << iAbsEta << std::endl;
     }
 
 
       // All RCT stuff uniform in phi, symmetric wrt eta = 0
 
       // below line always gives crate in +eta; makes no difference to us
       unsigned short crate = rctParameters_->calcCrate(1, iAbsEta); 
       unsigned short tower = rctParameters_->calcTower(1, iAbsEta);
 
       // first do region sums half of LUTs, bit 18 of address is 0
       // loop through 8 bits of hcal energy, 2^8 is 256
       for (unsigned int hcalEt = 0; hcalEt < 256; hcalEt++)
     {
       // loop through 1 bit of ecal fine grain
       for (unsigned short ecalfg = 0; ecalfg < 2; ecalfg++)
         {
           // loop through 8 bits of ecal energy
           for (unsigned int ecalEt = 0; ecalEt < 256; ecalEt++)
         {
           // assign 10-bit (9et,1mip) sums data here!
           unsigned long output = lookupTable_->
             lookup(ecalEt, hcalEt, ecalfg, crate, card, tower);
           unsigned short etIn9Bits = (output>>8)&511;
           unsigned short tauActivityBit = (output>>17)&1;
           data = (tauActivityBit<<9)+etIn9Bits;
           lutFile_ << std::hex << data << std::dec << std::endl;
         }
         }
     }
       // second do egamma half of LUTs, bit 18 of address is 1
       // loop through 8 bits of hcal energy
       for (unsigned int hcalEt = 0; hcalEt < 256; hcalEt++)
         {
           // loop through 1 bit of ecal fine grain
           for (unsigned short ecalfg = 0; ecalfg < 2; ecalfg++)
             {
               // loop through 8 bits of ecal energy
               for (unsigned int ecalEt = 0; ecalEt < 256; ecalEt++)
                 {
                   // assign 8-bit (7et,1veto) egamma data here!
           unsigned long output = lookupTable_->
             lookup(ecalEt, hcalEt, ecalfg, crate, card, tower);
           unsigned short etIn7Bits = output&127;
           unsigned short heFgVetoBit = (output>>7)&1;
           data = (heFgVetoBit<<7)+etIn7Bits;
           lutFile_ << std::hex << data << std::dec << std::endl;
                 }
             }
         }
     }
 
   lutFile_.close();
   return;
 
 }

void L1RCTLutWriter::writeThresholdsFile	(	unsigned int	eicThreshold,
		unsigned int	jscThresholdBarrel,
		unsigned int	jscThresholdEndcap
	)

private

Definition at line 482 of file L1RCTLutWriter.cc.

References corrVsCorr::filename, keyName_, and estimatePileup_makeJSON::trunc.

Referenced by analyze().

 {
   //
   std::ofstream thresholdsFile;
   char filename[256];
   sprintf(filename, "Thresholds-%s.dat", keyName_.c_str() );
   thresholdsFile.open(filename, std::ios::trunc);
 
   thresholdsFile << "key is " << keyName_ << std::endl << std::endl;
   thresholdsFile << "eicIsolationThreshold " << eicThreshold << std::endl;
   thresholdsFile << "jscQuietThresholdBarrel " << jscThresholdBarrel << std::endl;
   thresholdsFile << "jscQuietThresholdEndcap " << jscThresholdEndcap << std::endl;
 
   thresholdsFile.close();
 }