#include <EcalSimpleProducer.h>

Inheritance diagram for EcalSimpleProducer:

Public Member Functions
void	beginJob ()

	EcalSimpleProducer (const edm::ParameterSet &pset)

virtual void	produce (edm::Event &evt, const edm::EventSetup &)

virtual	~EcalSimpleProducer ()

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

ModuleDescription const &	moduleDescription () const

virtual	~EDProducer ()

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

	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
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 (std::string const &iProcessName, std::string const &iModuleLabel, bool iPrint, std::vector< char const * > &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
int	cIndex2iEta (int iEta0) const

int	cIndex2iPhi (int iPhi0) const

int	cIndex2iTtEta (int iEta0) const

int	cIndex2iTtPhi (int iPhi0) const

void	replaceAll (std::string &s, const std::string &from, const std::string &to) const

Private Attributes
std::auto_ptr< TFormula >	formula_

std::auto_ptr< TFormula >	simHitFormula_

std::auto_ptr< TFormula >	tpFormula_

bool	verbose_

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

Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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::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

This edm producer generates Ecal Digis (data frames and TPGs) according to a given pattern. The pattern is defined as function of event id, crystal/TT, and time sample. Only barrel is currently supported for the crystal channel data.

Module parameters (in addition to standard source parameters):

string formula: formula of crystal channel time sample encoded ADC counts.
string tpFormula: formula of trigger primitives.
untracked bool verbose: verbosity switch

The crystal ADC formula is parametrized with the following variables:

ieta0: crystal eta index starting from 0 at eta- end of barrel
iphi0: crystal phi index starting at Phi=0deg. in std CMS coordinates
ievt0 event sequence number within the job run starting from 0
isample0 sample time position starting from 0

The trigger primitive formula is parametrized with the following variables:

ieta0: trigger tower eta index starting from 0 at eta- end of barrel
iphi0: trigger tower index starting at Phi=0deg. in std CMS coordinates
ievt0 event sequence number within the job run starting from 0
isample0 sample time position starting from 0

In both formulae 'itt0' shortcut can be used for the trigger tower index within the SM starting at 0 from lowest relative eta and lowest phi and increasing first with phi then with eta. The syntax for the formula is the syntax defined in ROOT TFormula

Definition at line 38 of file EcalSimpleProducer.h.

Constructor & Destructor Documentation

EcalSimpleProducer::EcalSimpleProducer ( const edm::ParameterSet & pset )

Constructs an EcalSimpleProducer

Parameters

pset	CMSSW configuration
sdesc	description of this input source

Definition at line 97 of file EcalSimpleProducer.cc.

References Exception, EcalElecEmulTccOnlyExample_cfg::formula, formula_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), replaceAll(), simHitFormula_, EcalElecEmulTccOnlyExample_cfg::tpFormula, tpFormula_, and verbose_.

                                                                  :
   EDProducer(){
   string formula = pset.getParameter<string>("formula");
   string tpFormula = pset.getParameter<string>("tpFormula");
   string simHitFormula = pset.getParameter<string>("simHitFormula");
   
   verbose_ = pset.getUntrackedParameter<bool>("verbose", false);
   //  replaceAll(formula, "itt0", "((((ieta0<85)*(84-ieta0)+(ieta0>=85)*(ieta0-85))/5-18)*4+((iphi0/5+2)%4))");
   replaceAll(formula, "ebm", "(ieta0<85)");
   replaceAll(formula, "ebp", "(ieta0>84)");
   replaceAll(formula, "ieta0", "x");
   replaceAll(formula, "iphi0", "y");
   replaceAll(formula, "ievt0", "z");
   replaceAll(formula, "isample0", "t");
   //  cout << "----------> " << formula << endl;
   
   replaceAll(tpFormula, "itt0", "((ieta0<28)*(27-ieta0)+(ieta0>=28)*(ieta0-28))*4+(iphi0+2)%4");
   replaceAll(tpFormula, "eb", "(ieta0>10 && ieta0<45)");
   replaceAll(tpFormula, "ebm", "(ieta0>10 && ieta0<28)");
   replaceAll(tpFormula, "ebp", "(ieta0>27 && ieta0<45)");
   replaceAll(tpFormula, "ee", "(ieta0<11 || ieta0>44)");
   replaceAll(tpFormula, "eem", "(ieta0<11)");
   replaceAll(tpFormula, "eep", "(ieta0>44)");
   replaceAll(tpFormula, "ieta0", "x");
   replaceAll(tpFormula, "iphi0", "y");
   replaceAll(tpFormula, "ievt0", "z");
   replaceAll(tpFormula, "isample0", "t");
   //cout << "----------> " << tpFormula << endl;
   
 
   //replaceAll(simHitormula, "itt0", "((((ieta0<85)*(84-ieta0)+(ieta0>=85)*(ieta0-85))/5-18)*4+((iphi0/5+2)%4))");
   replaceAll(simHitFormula, "ebm", "(ieta0<85)");
   replaceAll(simHitFormula, "ebp", "(ieta0>84)");
   replaceAll(simHitFormula, "ieta0", "x");
   replaceAll(simHitFormula, "iphi0", "y");
   replaceAll(simHitFormula, "ievt0", "z");
   
   if(formula.size()!=0){
     formula_ = auto_ptr<TFormula>(new TFormula("f", formula.c_str()));
     Int_t err = formula_->Compile();
     if(err!=0){
       throw cms::Exception("Error in EcalSimpleProducer 'formula' config.");
     }
     produces<EBDigiCollection>();
     produces<EEDigiCollection>();
   }
   if(tpFormula.size()!=0){
     tpFormula_ = auto_ptr<TFormula>(new TFormula("f", tpFormula.c_str()));
     Int_t err = tpFormula_->Compile();
     if(err!=0){
       throw cms::Exception("Error in EcalSimpleProducer 'tpFormula' config.");
     }
     produces<EcalTrigPrimDigiCollection>();
   }
   if(simHitFormula.size()!=0){
     simHitFormula_
       = auto_ptr<TFormula>(new TFormula("f", simHitFormula.c_str()));
     Int_t err = simHitFormula_->Compile();
     if(err!=0){
       throw cms::Exception("Error in EcalSimpleProducer "
                "'simHitFormula' config.");
     }
     produces<edm::PCaloHitContainer>("EcalHitsEB");
     produces<edm::PCaloHitContainer>("EcalHitsEE");
   }
 }

virtual EcalSimpleProducer::~EcalSimpleProducer ( )

inlinevirtual

Destructor

Definition at line 50 of file EcalSimpleProducer.h.

50 {};

Member Function Documentation

void EcalSimpleProducer::beginJob ( void )

inlinevirtual

Called at start of job.

Parameters

es	the event setup

Reimplemented from edm::EDProducer.

Definition at line 55 of file EcalSimpleProducer.h.

55 {};

int EcalSimpleProducer::cIndex2iEta ( int iEta0 ) const

inlineprivate

Converts c-array index (contiguous integer starting from 0) to std CMSSW ECAL crystal eta index.

Parameters

iEta0 c-array index. '0' postfix reminds the index starts from 0

Returns: std CMSSW ECAL crystal index.

Definition at line 79 of file EcalSimpleProducer.h.

                                   {
     return (iEta0<85)?iEta0-85:iEta0-84;
   }

int EcalSimpleProducer::cIndex2iPhi ( int iPhi0 ) const

inlineprivate

Converts c-array index (contiguous integer starting from 0) to std CMSSW ECAL crystal phi index.

Parameters

iPhi0 c-array index. '0' postfix reminds the index starts from 0

Returns: std CMSSW ECAL crystal index.

Definition at line 88 of file EcalSimpleProducer.h.

                                   {
     return (iPhi0+10)%360+1;
   }

int EcalSimpleProducer::cIndex2iTtEta ( int iEta0 ) const

inlineprivate

Converts c-array index (contiguous integer starting from 0) to std CMSSW ECAL trigger tower eta index.

Parameters

iEta0 c-array index. '0' postfix reminds the index starts from 0

Returns: std CMSSW ECAL trigger tower index.

Definition at line 97 of file EcalSimpleProducer.h.

                                     {
     return (iEta0<28)?iEta0-28:iEta0-27;
   }

int EcalSimpleProducer::cIndex2iTtPhi ( int iPhi0 ) const

inlineprivate

Converts c-array index (contiguous integer starting from 0) to std CMSSW ECAL trigger tower phi index.

Parameters

iPhi0 c-array index. '0' postfix reminds the index starts from 0

Returns: std CMSSW ECAL trigger tower index.

Definition at line 106 of file EcalSimpleProducer.h.

                                     {
     return iPhi0+1;
   }

void EcalSimpleProducer::produce	(	edm::Event &	evt,
		const edm::EventSetup &
	)

virtual

The main method. It produces the event.

Parameters

evt	[out] produced event.

Implements edm::EDProducer.

Definition at line 16 of file EcalSimpleProducer.cc.

References funct::abs(), gather_cfg::cout, EcalTriggerTower, edm::EventID::event(), edm::EventBase::id(), hgc_digi::nSamples, edm::Event::put(), EcalTriggerPrimitiveDigi::setSample(), EcalTriggerPrimitiveDigi::setSize(), lumiQTWidget::t, and ecaldqm::zside().

                                                                    {
   const int ievt = evt.id().event();
   if(formula_.get()!=0){
     auto_ptr<EBDigiCollection> digis(new EBDigiCollection);
     
     digis->reserve(170*360);
     
     const int nSamples = digis->stride();
     for(int iEta0=0; iEta0<170; ++iEta0){
       for(int iPhi0=0; iPhi0<360; ++iPhi0){
     int iEta1 = cIndex2iEta(iEta0);
     int iPhi = cIndex2iPhi(iPhi0);
     if(verbose_) cout << "(" << iEta0 << "," << iPhi0 << "): ";
     digis->push_back(EBDetId(iEta1,iPhi));
     DataFrame dframe(digis->back());
     
     for(int t = 0; t < nSamples; ++t){
       uint16_t encodedAdc = (uint16_t)formula_->Eval(iEta0, iPhi0, ievt-1, t);
       if(verbose_) cout << encodedAdc << ((t<(nSamples-1))?"\t":"\n");
       dframe[t]=encodedAdc;
     }
       }
     }
     evt.put(digis);
     //puts an empty digi collecion for endcap:
     evt.put(auto_ptr<EEDigiCollection>(new EEDigiCollection()));
   }
   if(tpFormula_.get()!=0){
     auto_ptr<EcalTrigPrimDigiCollection> tps
       = auto_ptr<EcalTrigPrimDigiCollection>(new EcalTrigPrimDigiCollection);
     tps->reserve(56*72);
     const int nSamples = 5;
     for(int iTtEta0=0; iTtEta0<56; ++iTtEta0){
       for(int iTtPhi0=0; iTtPhi0<72; ++iTtPhi0){
     int iTtEta1 = cIndex2iTtEta(iTtEta0);
     int iTtPhi = cIndex2iTtPhi(iTtPhi0);
     
     if(verbose_) cout << "(" << iTtEta0 << "," << iTtPhi0 << "): ";
     int zside = iTtEta1<0?-1:1;
     EcalTriggerPrimitiveDigi
       tpframe(EcalTrigTowerDetId(zside, EcalTriggerTower,
                      abs(iTtEta1), iTtPhi));
     
     tpframe.setSize(nSamples);
     
     if(verbose_) cout << "TP: ";
     for(int t = 0; t < nSamples; ++t){
       uint16_t encodedTp = (uint16_t)tpFormula_->Eval(iTtEta0, iTtPhi0, ievt-1, t);
       
       if(verbose_) cout << "TP(" << iTtEta0 << "," << iTtPhi0 << ") = "
                 << encodedTp << ((t<(nSamples-1))?"\t":"\n");
       tpframe.setSample(t, EcalTriggerPrimitiveSample(encodedTp));
     }
     tps->push_back(tpframe);
       }
     }
     evt.put(tps);
   }
   if(simHitFormula_.get()!=0){//generation of barrel sim hits
     auto_ptr<PCaloHitContainer> hits
       = auto_ptr<PCaloHitContainer>(new PCaloHitContainer);
     for(int iEta0=0; iEta0<170; ++iEta0){
       for(int iPhi0=0; iPhi0<360; ++iPhi0){
     int iEta1 = cIndex2iEta(iEta0);
     int iPhi = cIndex2iPhi(iPhi0);
     if(verbose_) cout << "(" << iEta0 << "," << iPhi0 << "): ";
     
     double em = simHitFormula_->Eval(iEta0, iPhi0, ievt-1);
         double eh = 0.;
     double t = 0.;
     const PCaloHit hit(EBDetId(iEta1,iPhi).rawId(), em, eh, t, 0);
     hits->push_back(hit);
       }
     }
     evt.put(hits, "EcalHitsEB");
     //puts an empty digi collecion for endcap:
     evt.put(auto_ptr<PCaloHitContainer>(new PCaloHitContainer()),
         "EcalHitsEE");
   }
 }

void EcalSimpleProducer::replaceAll	(	std::string &	s,
		const std::string &	from,
		const std::string &	to
	)		const

private

Help function to replace a pattern within a string. Every occurance of the pattern is replaced. An exact match is performed: no wild card.

Parameters

s	string to operate on
pattern	to replace.
string	to substitute to the pattern

Definition at line 164 of file EcalSimpleProducer.cc.

Referenced by EcalSimpleProducer().

                                             {
   string::size_type pos = 0;
   //  cout << "replaceAll(" << s << "," << from << "," << to << ")\n";
   while((pos=s.find(from, pos))!=string::npos){
     // cout << "replace(" << pos << "," << from.size() << "," << to << ")\n";
     s.replace(pos, from.size(), to);
     //   cout << " -> " << s << "\n"; 
   }
 }