#include <EcalFEtoDigi.h>

Inheritance diagram for EcalFEtoDigi:

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

	~EcalFEtoDigi ()

Public Member Functions inherited from edm::EDProducer
	EDProducer ()

virtual	~EDProducer ()

Public Member Functions inherited from edm::ProducerBase
	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

boost::function< void(const BranchDescription &)>	registrationCallback () const
	used by the fwk to register list of products More...

virtual	~ProducerBase ()

Private Member Functions
virtual void	beginJob ()
	method called once each job just before starting event loop More...

EcalTriggerPrimitiveSample	create_TPSample (TCCinput, const edm::EventSetup &)
	create EcalTriggerPrimitiveSample from input data (line) More...

EcalTriggerPrimitiveSample	create_TPSampleTcp (TCCinput, const edm::EventSetup &)
	create EcalTriggerPrimitiveSample in tcp format (uncomrpessed energy) More...

EcalTrigTowerDetId	create_TTDetId (TCCinput)
	create EcalTrigTowerDetId from input data (line) More...

virtual void	endJob ()
	method called once each job just after ending the event loop More...

void	getLUT (unsigned int *lut, const int towerId, const edm::EventSetup &) const
	return the LUT from eventSetup More...

virtual void	produce (edm::Event &, const edm::EventSetup &)
	method called to produce the data More...

void	readInput ()
	open and read in input (flat) data file More...

int	SMidToTCCid (const int) const
	translate input supermodule id into TCC id (barrel) More...

Private Attributes
std::string	basename_

bool	debug_

int	fileEventOffset_

TCCInputData	inputdata_ [N_SM]

std::ofstream	outfile

bool	singlefile

int	sm_

bool	useIdentityLUT_

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
template<class TProducer , class TMethod >
void	callWhenNewProductsRegistered (TProducer *iProd, TMethod iMethod)

Detailed Description

Definition at line 36 of file EcalFEtoDigi.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 7 of file EcalFEtoDigi.cc.

References basename_, debug_, fileEventOffset_, edm::ParameterSet::getUntrackedParameter(), singlefile, sm_, and useIdentityLUT_.

                                                          {
   basename_           = iConfig.getUntrackedParameter<std::string>("FlatBaseName","ecal_tcc_");
   sm_                 = iConfig.getUntrackedParameter<int>("SuperModuleId",-1);
   fileEventOffset_    = iConfig.getUntrackedParameter<int>("FileEventOffset",0);
   useIdentityLUT_     = iConfig.getUntrackedParameter<bool>("UseIdentityLUT",false);
   debug_              = iConfig.getUntrackedParameter<bool>("debugPrintFlag", false);
 
   singlefile = (sm_==-1)?false:true;
 
   produces<EcalTrigPrimDigiCollection>();
   produces<EcalTrigPrimDigiCollection >("formatTCP");
  
 }

EcalFEtoDigi::~EcalFEtoDigi ( )

inline

Definition at line 41 of file EcalFEtoDigi.h.

41 {}

Member Function Documentation

void EcalFEtoDigi::beginJob ( void )

privatevirtual

method called once each job just before starting event loop

check SM numbering convetion: 1-38 [or -1 flag to indicate all sm's are to be read in]

debug: open file for recreating input copy

Reimplemented from edm::EDProducer.

Definition at line 295 of file EcalFEtoDigi.cc.

References debug_, edm::hlt::Exception, outfile, readInput(), and sm_.

                       {
 
   if(sm_!=-1 && (sm_<1 || sm_>36)) 
     throw cms::Exception("EcalFEtoDigiInvalidDetId") 
       << "EcalFEtoDigi: Adapt SM numbering convention.\n";
 
   if(debug_)
     outfile.open("inputcopy.txt");
 
   readInput();
   
 }  

EcalTriggerPrimitiveSample EcalFEtoDigi::create_TPSample	(	TCCinput	data,
		const edm::EventSetup &	evtSetup
	)

private

create EcalTriggerPrimitiveSample from input data (line)

setup look up table

compress energy 10 -> 8 bit

crate sample

Definition at line 237 of file EcalFEtoDigi.cc.

References create_TTDetId(), debug_, TCCinput::get_energy(), TCCinput::get_fg(), getLUT(), i, TCCinput::input, DetId::rawId(), and useIdentityLUT_.

Referenced by produce().

                                                                           {
 
   int tower      = create_TTDetId(data).rawId() ; 
   int  Et        = data.get_energy();
   bool tt_fg     = data.get_fg();
   //unsigned input = data.input;
   //int  Et    = input & 0x3ff; //get bits 0-9
   //bool tt_fg = input & 0x400; //get bit number 10
 
   unsigned int lut_[1024] ;
   if(!useIdentityLUT_)
     getLUT(lut_, tower, evtSetup) ;
   else
     for(int i=0; i<1024; i++) lut_[i] = i ; //identity lut!
   
   int lut_out = lut_[Et];
   int ttFlag  = (lut_out & 0x700) >> 8;
   int cEt     = (lut_out & 0xff );
 
   if(debug_&&data.get_energy()!=0)
     printf("[EcalFEtoDigi] Creating sample; input:0x%X (Et:0x%x) cEt:0x%x fg:%d ttflag:0x%x \n",
        data.input, Et, cEt, tt_fg, ttFlag);
   
   EcalTriggerPrimitiveSample e_sample(cEt, tt_fg, ttFlag);
   
   return e_sample;
 }

EcalTriggerPrimitiveSample EcalFEtoDigi::create_TPSampleTcp	(	TCCinput	data,
		const edm::EventSetup &	evtSetup
	)

private

create EcalTriggerPrimitiveSample in tcp format (uncomrpessed energy)

setup look up table

Definition at line 270 of file EcalFEtoDigi.cc.

References create_TTDetId(), TCCinput::get_energy(), TCCinput::get_fg(), getLUT(), i, DetId::rawId(), and useIdentityLUT_.

Referenced by produce().

                                                                              {
 
   int tower      = create_TTDetId(data).rawId() ; 
   int  Et        = data.get_energy();
   bool tt_fg     = data.get_fg();
 
   unsigned int lut_[1024] ;
   if(!useIdentityLUT_)
     getLUT(lut_, tower, evtSetup) ;
   else
     for(int i=0; i<1024; i++) lut_[i] = i ; //identity lut!
   
   int lut_out = lut_[Et];
   int ttFlag  = (lut_out & 0x700) >> 8;
   int tcpdata = ((ttFlag&0x7)<<11) | ((tt_fg & 0x1)<<10) |  (Et & 0x3ff) ;
 
   EcalTriggerPrimitiveSample e_sample(tcpdata);
   
   return e_sample;
 }

EcalTrigTowerDetId EcalFEtoDigi::create_TTDetId ( TCCinput data )

private

create EcalTrigTowerDetId from input data (line)

construct the EcalTrigTowerDetId object

Definition at line 203 of file EcalFEtoDigi.cc.

References debug_, EcalBarrel, TCCinput::get_energy(), sm_, and TCCinput::tower.

Referenced by create_TPSample(), create_TPSampleTcp(), and produce().

                                           {
 
   // (EcalBarrel only)
   static const int kTowersInPhi = 4;
   
   int iTT   = data.tower;
   int zside = (sm_>18)?-1:+1;
   int SMid  = sm_;
 
   int jtower = iTT-1;
   int etaTT  = jtower / kTowersInPhi +1;
   int phiTT;
   if (zside < 0) 
     phiTT = (SMid-19) * kTowersInPhi + jtower % kTowersInPhi;
   else 
     phiTT = (SMid- 1) * kTowersInPhi + kTowersInPhi-(jtower % kTowersInPhi)-1;
 
   phiTT ++;
   //needed as phi=0 (iphi=1) is at middle of lower SMs (1 and 19), need shift by 2
   phiTT = phiTT -2; 
   if (phiTT <= 0) phiTT = 72+phiTT;
 
   if(debug_&&data.get_energy()!=0)
     printf("[EcalFEtoDigi] Creating EcalTrigTowerDetId (SMid,itt)=(%d,%d)->(eta,phi)=(%d,%d) \n", SMid, iTT, etaTT, phiTT);
   
   EcalTrigTowerDetId 
     e_id( zside , EcalBarrel, etaTT, phiTT, 0);
   
   return e_id;
 }

void EcalFEtoDigi::endJob ( void )

privatevirtual

method called once each job just after ending the event loop

Reimplemented from edm::EDProducer.

Definition at line 314 of file EcalFEtoDigi.cc.

References outfile.

                      {
   if(outfile.is_open())
     outfile.close();
 }

void EcalFEtoDigi::getLUT	(	unsigned int *	lut,
		const int	towerId,
		const edm::EventSetup &	evtSetup
	)		const

private

return the LUT from eventSetup

Definition at line 328 of file EcalFEtoDigi.cc.

References edm::EventSetup::get(), edm::eventsetup::EventSetupRecord::get(), i, and edm::ESHandle< class >::product().

Referenced by create_TPSample(), and create_TPSampleTcp().

 {
   edm::ESHandle<EcalTPGLutGroup> lutGrpHandle;
   evtSetup.get<EcalTPGLutGroupRcd>().get( lutGrpHandle );
   const EcalTPGGroups::EcalTPGGroupsMap & lutGrpMap = lutGrpHandle.product()->getMap() ;  
   EcalTPGGroups::EcalTPGGroupsMapItr itgrp = lutGrpMap.find(towerId) ;
   uint32_t lutGrp = 999 ;
   if (itgrp != lutGrpMap.end()) lutGrp = itgrp->second ;
   
   edm::ESHandle<EcalTPGLutIdMap> lutMapHandle;
   evtSetup.get<EcalTPGLutIdMapRcd>().get( lutMapHandle );
   const EcalTPGLutIdMap::EcalTPGLutMap & lutMap = lutMapHandle.product()->getMap() ;  
   EcalTPGLutIdMap::EcalTPGLutMapItr itLut = lutMap.find(lutGrp) ;
   if (itLut != lutMap.end()) {
     const unsigned int * theLut = (itLut->second).getLut() ;
     for (unsigned int i=0 ; i<1024 ; i++) lut[i] = theLut[i] ;
   }
 }

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

privatevirtual

method called to produce the data

event counter

re-read input (needed in case of event-by-event input production)

create EcalTrigTowerDetId

create EcalTriggerPrimitiveDigi

create EcalTriggerPrimitiveSample

set sample

add to EcalTrigPrimDigiCollection

set sample (uncompressed format)

add to EcalTrigPrimDigiCollection (uncompressed format)

print & debug

in case no info was found for the event:need to create something

Implements edm::EDProducer.

Definition at line 23 of file EcalFEtoDigi.cc.

References begin, EcalTriggerPrimitiveDigi::compressedEt(), gather_cfg::cout, create_TPSample(), create_TPSampleTcp(), create_TTDetId(), debug_, fileEventOffset_, EcalTriggerPrimitiveDigi::fineGrain(), i, EcalTriggerPrimitiveDigi::id(), EcalTrigTowerDetId::ieta(), inputdata_, EcalTrigTowerDetId::iphi(), N_SM, outfile, edm::Event::put(), EcalTriggerPrimitiveDigi::setSample(), EcalTriggerPrimitiveDigi::setSize(), singlefile, EcalTriggerPrimitiveDigi::size(), sm_, and EcalTriggerPrimitiveDigi::ttFlag().

 {
 
   static int current_bx = -1;
   current_bx++;
 
   //readInput();
 
   if(debug_)
     std::cout << "[EcalFEtoDigi::produce] producing event " << current_bx << std::endl;
   
   std::auto_ptr<EcalTrigPrimDigiCollection>  
     e_tpdigis (new EcalTrigPrimDigiCollection);
   std::auto_ptr<EcalTrigPrimDigiCollection>  
     e_tpdigisTcp (new EcalTrigPrimDigiCollection);
 
   std::vector<TCCinput>::const_iterator it;
 
   for(int i=0; i<N_SM; i++) {
 
     if(!singlefile)
       sm_=i+1;
 
     for(it = inputdata_[i].begin(); it != inputdata_[i].end(); it++) {
 
       if(!(*it).is_current(current_bx+fileEventOffset_)) 
     continue;
       else
         if(debug_ && (*it).input!=0 )
       std::cout << "[EcalFEtoDigi] " 
            << "\tsupermodule:" << sm_ 
            << "\tevent: "      << current_bx 
            << "\tbx: "         << (*it).bunchCrossing
            << "\tvalue:0x"     << std::setfill('0') << std::setw(4) 
            << std::hex << (*it).input << std::setfill(' ') << std::dec 
            << std::endl;
 
       
       const EcalTrigTowerDetId  e_id = create_TTDetId(*it);
       
       //EcalElectronicsMapping theMapping;
       //const EcalTrigTowerDetId  e_id 
       //= theMapping.getTrigTowerDetId(SMidToTCCid(sm_),(*it).tower);
       //EcalElectronicsMapping::getTrigTowerDetId(int TCCid, int iTT)
       
       EcalTriggerPrimitiveDigi *e_digi = new EcalTriggerPrimitiveDigi(e_id);
       EcalTriggerPrimitiveDigi *e_digiTcp = new EcalTriggerPrimitiveDigi(e_id);
       
       EcalTriggerPrimitiveSample e_sample = create_TPSample(*it, iSetup);
       EcalTriggerPrimitiveSample e_sampleTcp = create_TPSampleTcp(*it, iSetup);
 
       e_digi->setSize(1); //set sampleOfInterest to 0
       e_digi->setSample(0,e_sample);
       
       e_tpdigis->push_back(*e_digi);
     
       e_digiTcp->setSize(1); //set sampleOfInterest to 0
       e_digiTcp->setSample(0,e_sampleTcp);
       
       e_tpdigisTcp->push_back(*e_digiTcp);
 
       if(debug_) 
     outfile << (*it).tower << '\t' 
         << (*it).bunchCrossing << '\t'<< std::setfill('0') << std::hex
         << "0x" << std::setw(4) << (*it).input << '\t'
         << "0"  << std::dec << std::setfill(' ') 
         << std::endl;
 
       if(debug_ && (*it).input!=0 )
     std::cout << "[EcalFEtoDigi] debug id: " << e_digi->id() << "\n\t" 
           << std::dec 
           << "\tieta: "     << e_digi->id().ieta()
           << "\tiphi: "     << e_digi->id().iphi()
           << "\tsize: "     << e_digi->size()
           << "\tfg: "       <<(e_digi->fineGrain()?1:0)    
           << std::hex    
           << "\tEt: 0x"     << e_digi->compressedEt() 
           << " (0x"         << (*it).get_energy() << ")" 
           << "\tttflag: 0x" << e_digi->ttFlag()
           << std::dec
           << std::endl;
 
       delete e_digi;
       delete e_digiTcp;
 
     }
     
     if(singlefile)
       break;
   }
 
   if(e_tpdigis->size()==0) {
     std::cout << "[EcalFEtoDigi] creating empty collection for the event!\n";
     EcalTriggerPrimitiveDigi *e_digi = new EcalTriggerPrimitiveDigi();
     e_tpdigis->push_back(*e_digi);
   }
 
   iEvent.put(e_tpdigis);
   iEvent.put(e_tpdigisTcp,"formatTCP");
 
 }

void EcalFEtoDigi::readInput ( )

private

open and read in input (flat) data file

Definition at line 139 of file EcalFEtoDigi.cc.

References prof2calltree::back, basename_, hitfit::clear(), gather_cfg::cout, debug_, f, fileEventOffset_, i, inputdata_, N_SM, alignCSCRings::s, singlefile, sm_, SMidToTCCid(), and groupFilesInBlocks::tt.

Referenced by beginJob(), and parserPerfsuiteMetadata.parserPerfsuiteMetadata::readCmsScimarkTest().

                         {
 
   if(debug_)
     std::cout << "\n[EcalFEtoDigi::readInput] Reading input data\n";
   
   if(!singlefile)
     sm_=-1;
   for(int i=0; i<N_SM; i++)
     inputdata_[i].clear();
 
   std::stringstream s;
   int tcc;
 
   for (int i=0; i<N_SM; i++) {
 
     tcc = (sm_==-1)?SMidToTCCid(i+1):SMidToTCCid(sm_);
 
     s.str("");
     s << basename_ << tcc << ".txt"; 
 
     std::ifstream f(s.str().c_str());
     
     if(debug_) {
       std::cout << "  opening " << s.str().c_str() << "..." << std::endl;
       if(!f.good())
     std::cout << " skipped!"; 
       std::cout << std::endl;       
     }
     //if (!f.good() || f.eof()) 
     //  throw cms::Exception("BadInputFile") 
     //  << "EcalFEtoDigi: cannot open file " << s.str().c_str() << std::endl; 
     
     int n_bx=0;
     int tt; int bx; unsigned val; int dummy;
 
     while(f.good()) {
       if(f.eof()) break;
       tt=0; bx=-1; val=0x0; dummy=0;
       f >> tt >> bx >> std::hex >> val >> std::dec >> dummy;
       if(bx==-1 || bx < fileEventOffset_ ) continue;
       if( !n_bx || (bx!=(inputdata_[i].back()).bunchCrossing) )
     n_bx++;
       TCCinput ttdata(tt,bx,val);
       inputdata_[i].push_back(ttdata);
       
       if(debug_&&val!=0)
     printf("\treading tower:%d  bx:%d input:0x%x dummy:%2d\n", 
            tt, bx, val, dummy);
     } 
     
     f.close();
     
     if(sm_!=-1)
       break;    
   }
 
   if(debug_)
     std::cout << "[EcalFEtoDigi::readInput] Done reading." << std::endl;
 
   return; 
 }

int EcalFEtoDigi::SMidToTCCid ( const int smid ) const

private

translate input supermodule id into TCC id (barrel)

Definition at line 321 of file EcalFEtoDigi.cc.

Referenced by readInput().

                                                 {
 
   return (smid<=18) ? smid+55-1 : smid+37-19;
 
 }

Member Data Documentation

std::string EcalFEtoDigi::basename_

private

Definition at line 58 of file EcalFEtoDigi.h.

Referenced by EcalFEtoDigi(), and readInput().

bool EcalFEtoDigi::debug_

private

Definition at line 64 of file EcalFEtoDigi.h.

Referenced by beginJob(), create_TPSample(), create_TTDetId(), EcalFEtoDigi(), produce(), and readInput().

int EcalFEtoDigi::fileEventOffset_

private

Definition at line 63 of file EcalFEtoDigi.h.

Referenced by EcalFEtoDigi(), produce(), and readInput().

TCCInputData EcalFEtoDigi::inputdata_[N_SM]

private

Definition at line 56 of file EcalFEtoDigi.h.

Referenced by produce(), and readInput().

std::ofstream EcalFEtoDigi::outfile

private

Definition at line 65 of file EcalFEtoDigi.h.

Referenced by beginJob(), endJob(), and produce().

bool EcalFEtoDigi::singlefile

private

Definition at line 61 of file EcalFEtoDigi.h.

Referenced by EcalFEtoDigi(), produce(), and readInput().

int EcalFEtoDigi::sm_

private

Definition at line 60 of file EcalFEtoDigi.h.

Referenced by beginJob(), create_TTDetId(), EcalFEtoDigi(), produce(), and readInput().

bool EcalFEtoDigi::useIdentityLUT_

private

Definition at line 59 of file EcalFEtoDigi.h.

Referenced by create_TPSample(), create_TPSampleTcp(), and EcalFEtoDigi().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation