RctInputTextToDigi Class Reference

#include <L1Trigger/TextToDigi/plugins/RctInputTextToDigi.h>

Inheritance diagram for RctInputTextToDigi:

Public Member Functions
	RctInputTextToDigi (const edm::ParameterSet &)
	~RctInputTextToDigi ()
Public Member Functions inherited from edm::EDProducer
	EDProducer ()
virtual	~EDProducer ()
Public Member Functions inherited from edm::ProducerBase
	ProducerBase ()
void	registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...
virtual	~ProducerBase ()
Public Member Functions inherited from edm::EDConsumerBase
	EDConsumerBase ()
ProductHolderIndex	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductHolderIndex > &) const
void	itemsToGet (BranchType, std::vector< ProductHolderIndex > &) const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
virtual	~EDConsumerBase ()

Private Member Functions
virtual void	beginJob ()
virtual void	endJob ()
virtual void	produce (edm::Event &, const edm::EventSetup &)

Private Attributes
edm::FileInPath	inputFile_
std::ifstream	inputStream_
L1RCTLookupTables *	lookupTables_
int	nEvent_
bool	oldVersion_

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType
typedef WorkerT< EDProducer >	WorkerType
Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList
Static Public Member Functions inherited from edm::EDProducer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
Protected Member Functions inherited from edm::EDProducer
CurrentProcessingContext const *	currentContext () const
Protected Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
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: Creates an EcalTrigPrimDigiCollection and an HcalTrigPrimDigiCollection from a text file (formatted as read out from saveRCTInput module), for input to the L1Trigger/RegionalCaloTrigger module.

Definition at line 50 of file RctInputTextToDigi.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 7 of file RctInputTextToDigi.cc.

References dtNoiseDBValidation_cfg::cerr, and inputStream_.

                                                                     :
  inputFile_(iConfig.getParameter<edm::FileInPath>("inputFile")),
  inputStream_(inputFile_.fullPath().c_str()),
  nEvent_(0),
  oldVersion_(false),
  lookupTables_(new L1RCTLookupTables)
{
  //register your products
  /* Examples
     produces<ExampleData2>();
     
     //if do put with a label
     produces<ExampleData2>("label");
  */
  
  produces<EcalTrigPrimDigiCollection>();
  produces<HcalTrigPrimDigiCollection>();

  //now do what ever other initialization is needed

  if ((!inputStream_.is_open())||(!inputStream_))
    {
      // not good!!
      std::cerr << "Input file didn't open!!" << std::endl;
    }
  //if (inputStream_.eof()) {std::cout << "Real zeroth eof! " << std::endl;}

}

RctInputTextToDigi::~RctInputTextToDigi

(

)

Definition at line 37 of file RctInputTextToDigi.cc.

References inputStream_.

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

  inputStream_.close();

}

Member Function Documentation

void RctInputTextToDigi::beginJob ( void  )

privatevirtual

Reimplemented from edm::EDProducer.

Definition at line 221 of file RctInputTextToDigi.cc.

{
  // open input file to read all events
  //inputStream_.open(inputFile_.fullPath().c_str());
  //std::cout << "beginJob entered" << std::endl;
}

void RctInputTextToDigi::endJob ( void  )

privatevirtual

Reimplemented from edm::EDProducer.

Definition at line 230 of file RctInputTextToDigi.cc.

{
  // close input file
  //inputStream_.close();
}

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

privatevirtual

Implements edm::EDProducer.

Definition at line 54 of file RctInputTextToDigi.cc.

References abs, L1RCTParameters::calcIEta(), L1RCTParameters::calcIPhi(), gather_cfg::cout, EcalTriggerTower, edm::EventSetup::get(), i, inputStream_, j, lookupTables_, nEvent_, oldVersion_, edm::Event::put(), alignCSCRings::r, L1RCTLookupTables::rctParameters(), L1RCTLookupTables::setRCTParameters(), EcalTriggerPrimitiveDigi::setSample(), HcalTriggerPrimitiveDigi::setSample(), EcalTriggerPrimitiveDigi::setSize(), HcalTriggerPrimitiveDigi::setSize(), and AlCaHLTBitMon_QueryRunRegistry::string.

{
  using namespace edm;

  //std::cout << std::endl << std::endl << "Event number " << nEvent_ << std::endl;

  // This next section taken directly from 
  // L1Trigger/RegionalCaloTrigger/plugins/L1RCTProducer.cc rev. 1.6
  // Refresh configuration information every event
  // Hopefully doesn't take too much time
  ESHandle<L1RCTParameters> rctParameters;
  iSetup.get<L1RCTParametersRcd>().get(rctParameters);
  const L1RCTParameters* r = rctParameters.product();
  lookupTables_->setRCTParameters(r);

  std::auto_ptr<EcalTrigPrimDigiCollection>
    ecalTPs(new EcalTrigPrimDigiCollection());
  std::auto_ptr<HcalTrigPrimDigiCollection> 
    hcalTPs(new HcalTrigPrimDigiCollection());
  ecalTPs->reserve(56*72); 
  hcalTPs->reserve(56*72+18*8);  // includes HF
  const int nEcalSamples = 1;   // we only use 1 sample for each
  const int nHcalSamples = 1;
  
  int fileEventNumber; 
  
  // check to see if need to skip file header and do so before
  // looping through entire event
  
  std::string junk;
  //bool old_version = false;
  if (nEvent_ == 0)
    {
      //std::string junk;
      unsigned short junk_counter = 0;
      //bool old_version = false;
      do
    {
      if(inputStream_ >> junk) {/*std::cout << "Good: ";*/}
      //std::cout << "header junk was input: \"" << junk << "\"."
      //             << std::endl;
      // for oldest version, which is same as newest version
      //      if((junk_counter == 11) && (junk.compare("0-32") == 0))
      //        {
      //          oldVersion_ = true;
      //        }
      if((junk_counter == 11) && (junk.compare("1-32") == 0))
        {
          oldVersion_ = true;
        }
      junk_counter++;
    }
      while (junk.compare("LUTOut") != 0);
      std::cout << "Skipped file header" << std::endl;
      if (oldVersion_) {std::cout << "oldVersion_ TRUE (tower 1-32)" << std::endl;}
      else {std::cout << "oldVersion_ FALSE (tower 0-31)" << std::endl;}
    }
  
  // can't actually read in phi and eta, file has crate card tower instead
  // do a while loop for event number instead??  dunno
  for (int i = 0; i < 72; i++)
    {
      // negative eta, iEta -28 to -1
      for (int j = 0; j < 56; j++)
    {
      // calc ieta, iphi coords of tower 
      // ieta -28 to -1 or 1 to 28, iphi 1 to 72
      // methods in CondFormats/L1TObjects/src/L1RCTParameters.cc
      
      unsigned short crate;
      unsigned short card;
      unsigned short tower;
      unsigned eAddr;
      unsigned hAddr;
      
      inputStream_ >> std::hex >> fileEventNumber >> crate >> card  
               >> tower
               >> eAddr >> hAddr >> junk >> std::dec;
      
      if (oldVersion_)
        {
          tower = tower - 1;
        }
      int encodedEtEcal = (int) (eAddr>>1);
      bool fineGrainEcal = (bool) (eAddr&1);
      int encodedEtHcal = (int) (hAddr>>1);
      bool fineGrainHcal = (bool) (hAddr&1);  // mip bit
      
      //std::cout << "Eventnumber " << fileEventNumber << "\tCrate " 
      //    << crate << "\tCard " << card << "\tTower " 
      //    << tower << " \teAddr " << eAddr <<"\thAddr "
      //    << hAddr << "\tjunk " << junk << std::endl;
      
      int iEta = lookupTables_->rctParameters()->calcIEta(crate,card,tower);
      int iPhi = lookupTables_->rctParameters()->calcIPhi(crate,card,tower);
      // transform rct iphi coords into global coords
      iPhi = ((72 + 18 - iPhi) % 72);
      if (iPhi == 0) {iPhi = 72;}
      unsigned absIeta = abs(iEta);
      int zSide = (iEta/absIeta);
      
      /*std::cout << "iEta " << iEta << "\tabsiEta " << absIeta 
        << "\tiPhi " << iPhi << "\tzSide " 
        << zSide << std::endl;
      */
      
      // args to detid are zside, type of tower, absieta, iphi
      // absieta and iphi must be between 1 and 127 inclusive

      EcalTriggerPrimitiveDigi 
        ecalDigi(EcalTrigTowerDetId(zSide, EcalTriggerTower, absIeta,
                    iPhi));
      ecalDigi.setSize(nEcalSamples);
      
      // last arg is 3-bit trigger tower flag, which we don't use
      // we only use 8-bit encoded et and 1-bit fg
      ecalDigi.setSample(0, EcalTriggerPrimitiveSample(encodedEtEcal, 
                               fineGrainEcal, 0));
      //std::cout << ecalDigi << std::endl;
      ecalTPs->push_back(ecalDigi);
      
      HcalTriggerPrimitiveDigi
        hcalDigi(HcalTrigTowerDetId(iEta, iPhi));
      
      hcalDigi.setSize(nHcalSamples);
      
      // last two arg's are slb and slb channel, which we don't need
      hcalDigi.setSample(0, HcalTriggerPrimitiveSample(encodedEtHcal,
                               fineGrainHcal,
                               0, 0));
      //std::cout << hcalDigi << std::endl;
      hcalTPs->push_back(hcalDigi);
    }
      
      // also need to push_back HF digis!  
      // file input doesn't include HF, so need empty digis
      for (int i = 0; i < 18; i++)
    {
      for (int j = 0; j < 8; j++)
        {
          // HF ieta: +- 29 through 32.  HF iphi: 1,5,9,13,etc.
          int hfIEta = (j%4)+29;
          if (i < 9)
        {
          hfIEta = hfIEta*(-1);
        }
          // iphi shift not implemented, but not necessary here --
          // everything's filled with zeros so it's symmetric anyhow
          int hfIPhi = (i%9)*8 + (j/4)*4 + 1;
          
          HcalTriggerPrimitiveDigi
        hfDigi(HcalTrigTowerDetId(hfIEta, hfIPhi));
          hfDigi.setSize(1);
          hfDigi.setSample(0, HcalTriggerPrimitiveSample(0,0,0,0));
          hcalTPs->push_back(hfDigi);
        }
    }
    }
  iEvent.put(ecalTPs);
  iEvent.put(hcalTPs);
  
  nEvent_++;
  //std::cout << "Produce done" << std::endl;
}

Member Data Documentation

edm::FileInPath RctInputTextToDigi::inputFile_

private

Definition at line 62 of file RctInputTextToDigi.h.

std::ifstream RctInputTextToDigi::inputStream_

private

Definition at line 63 of file RctInputTextToDigi.h.

Referenced by produce(), RctInputTextToDigi(), and ~RctInputTextToDigi().

L1RCTLookupTables* RctInputTextToDigi::lookupTables_

private

Definition at line 66 of file RctInputTextToDigi.h.

Referenced by produce().

int RctInputTextToDigi::nEvent_

private

Definition at line 64 of file RctInputTextToDigi.h.

Referenced by produce().

bool RctInputTextToDigi::oldVersion_

private

Definition at line 65 of file RctInputTextToDigi.h.

Referenced by produce().