#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 ()
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_

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, EcalElecEmulExample_cfg::formula, formula_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), replaceAll(), EcalElecEmulExample_cfg::simHitFormula, simHitFormula_, EcalElecEmulExample_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 ( ) [inline, virtual]

Destructor

Definition at line 50 of file EcalSimpleProducer.h.

{};

Member Function Documentation

void EcalSimpleProducer::beginJob ( void ) [inline, virtual]

Called at start of job.

Parameters:

es	the event setup

Reimplemented from edm::EDProducer.

Definition at line 55 of file EcalSimpleProducer.h.

{};

int EcalSimpleProducer::cIndex2iEta ( int iEta0 ) const [inline, private]

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 [inline, private]

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 [inline, private]

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 [inline, private]

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 abs, gather_cfg::cout, EcalTriggerTower, edm::EventID::event(), edm::EventBase::id(), edm::Event::put(), EcalTriggerPrimitiveDigi::setSample(), EcalTriggerPrimitiveDigi::setSize(), and lumiQTWidget::t.

                                                                   {
  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.

References pos.

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"; 
  }
}