HcalSimpleReconstructor Class Reference

#include <HcalSimpleReconstructor.h>

Inheritance diagram for HcalSimpleReconstructor:

Public Member Functions
virtual void	beginRun (edm::Run &r, edm::EventSetup const &es)
virtual void	endRun (edm::Run &r, edm::EventSetup const &es)
	HcalSimpleReconstructor (const edm::ParameterSet &ps)
virtual void	produce (edm::Event &e, const edm::EventSetup &c)
virtual	~HcalSimpleReconstructor ()
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 Attributes
DetId::Detector	det_
bool	dropZSmarkedPassed_
int	firstSample_
edm::InputTag	inputLabel_
HcalRecoParams *	paramTS
HcalSimpleRecAlgo	reco_
int	samplesToAdd_
int	subdet_
HcalOtherSubdetector	subdetOther_
bool	tsFromDB_

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

Date:

2011/02/22 20:44:52

Revision:

1.4

Author

J. Mans - Minnesota

Date:

2010/01/21 14:37:40

Revision:

1.1

Author

E. Garcia - CSU Based on HcalSimpleReconstructor.h by J. Mans

Definition at line 22 of file HcalSimpleReconstructor.h.

Constructor & Destructor Documentation

HcalSimpleReconstructor::HcalSimpleReconstructor ( const edm::ParameterSet &  ps )

explicit

Definition at line 16 of file HcalSimpleReconstructor.cc.

References gather_cfg::cout, edm::ParameterSet::getParameter(), HcalBarrel, HcalForward, HcalOuter, and subdet_.

                                                                           :
  reco_(conf.getParameter<bool>("correctForTimeslew"),
    conf.getParameter<bool>("correctForPhaseContainment"),conf.getParameter<double>("correctionPhaseNS")),
  det_(DetId::Hcal),
  inputLabel_(conf.getParameter<edm::InputTag>("digiLabel")),
  dropZSmarkedPassed_(conf.getParameter<bool>("dropZSmarkedPassed")),
  firstSample_(conf.getParameter<int>("firstSample")),
  samplesToAdd_(conf.getParameter<int>("samplesToAdd")),
  tsFromDB_(conf.getParameter<bool>("tsFromDB"))
{

  std::string subd=conf.getParameter<std::string>("Subdetector");
  if (!strcasecmp(subd.c_str(),"HBHE")) {
    subdet_=HcalBarrel;
    produces<HBHERecHitCollection>();
  } else if (!strcasecmp(subd.c_str(),"HO")) {
    subdet_=HcalOuter;
    produces<HORecHitCollection>();
  } else if (!strcasecmp(subd.c_str(),"HF")) {
    subdet_=HcalForward;
    produces<HFRecHitCollection>();
  } else {
    std::cout << "HcalSimpleReconstructor is not associated with a specific subdetector!" << std::endl;
  }       
  
}

HcalSimpleReconstructor::~HcalSimpleReconstructor ( )

virtual

Definition at line 43 of file HcalSimpleReconstructor.cc.

{ }

Member Function Documentation

void HcalSimpleReconstructor::beginRun ( edm::Run &  r, edm::EventSetup const &  es  )

virtual

Reimplemented from edm::EDProducer.

Definition at line 45 of file HcalSimpleReconstructor.cc.

References edm::EventSetup::get(), L1TEmulatorMonitor_cff::p, paramTS, and edm::ESHandle< class >::product().

                                                                        {

    edm::ESHandle<HcalRecoParams> p;
    es.get<HcalRecoParamsRcd>().get(p);
    paramTS = new HcalRecoParams(*p.product());

}

void HcalSimpleReconstructor::endRun ( edm::Run &  r, edm::EventSetup const &  es  )

virtual

Reimplemented from edm::EDProducer.

Definition at line 53 of file HcalSimpleReconstructor.cc.

References paramTS.

                                                                      {
  if (paramTS) delete paramTS;
}

void HcalSimpleReconstructor::produce ( edm::Event &  e, const edm::EventSetup &  c  )

virtual

Implements edm::EDProducer.

Definition at line 60 of file HcalSimpleReconstructor.cc.

References det_, dropZSmarkedPassed_, first, HcalRecoParam::firstSample(), firstSample_, edm::EventSetup::get(), edm::Event::getByLabel(), HcalCondObjectContainer< Item >::getValues(), DetId::Hcal, HcalBarrel, HcalCalibration, HcalEndcap, HcalForward, HcalOther, HcalOuter, i, HcalSimpleRecAlgo::initPulseCorr(), inputLabel_, edm::EventBase::isRealData(), paramTS, edm::Event::put(), DetId::rawId(), reco_, HcalSimpleRecAlgo::reconstruct(), HcalRecoParam::samplesToAdd(), samplesToAdd_, HcalSimpleRecAlgo::setForData(), subdet_, subdetOther_, and tsFromDB_.

{
  // get conditions
  edm::ESHandle<HcalDbService> conditions;
  eventSetup.get<HcalDbRecord>().get(conditions);
  const HcalQIEShape* shape = conditions->getHcalShape (); // this one is generic

  // HACK related to HB- corrections
  if(e.isRealData()) reco_.setForData();
 
  
  if (det_==DetId::Hcal) {
    if (subdet_==HcalBarrel || subdet_==HcalEndcap) {
      edm::Handle<HBHEDigiCollection> digi;
      
      e.getByLabel(inputLabel_,digi);
      
      // create empty output
      std::auto_ptr<HBHERecHitCollection> rec(new HBHERecHitCollection);
      rec->reserve(digi->size());
      // run the algorithm
      int toaddMem = 0;
      int first = firstSample_;
      int toadd = samplesToAdd_;
      HBHEDigiCollection::const_iterator i;
      for (i=digi->begin(); i!=digi->end(); i++) {
    HcalDetId cell = i->id();
    DetId detcell=(DetId)cell;
    // rof 27.03.09: drop ZS marked and passed digis:
    if (dropZSmarkedPassed_)
      if (i->zsMarkAndPass()) continue;

    const HcalCalibrations& calibrations=conditions->getHcalCalibrations(cell);
    const HcalQIECoder* channelCoder = conditions->getHcalCoder (cell);
    HcalCoderDb coder (*channelCoder, *shape);

    //>>> firstSample & samplesToAdd
        if(tsFromDB_) {
      const HcalRecoParam* param_ts = paramTS->getValues(detcell.rawId());
      first = param_ts->firstSample();    
      toadd = param_ts->samplesToAdd();
    }    
        if(toaddMem != toadd) {
      reco_.initPulseCorr(toadd);
          toaddMem = toadd;
    }

    rec->push_back(reco_.reconstruct(*i,first,toadd,coder,calibrations));

      }
      // return result
      e.put(rec);
    } else if (subdet_==HcalOuter) {
      edm::Handle<HODigiCollection> digi;
      e.getByLabel(inputLabel_,digi);
      
      // create empty output
      std::auto_ptr<HORecHitCollection> rec(new HORecHitCollection);
      rec->reserve(digi->size());
      // run the algorithm
      int toaddMem = 0;
      int first = firstSample_;
      int toadd = samplesToAdd_;
      HODigiCollection::const_iterator i;
      for (i=digi->begin(); i!=digi->end(); i++) {
    HcalDetId cell = i->id();
    DetId detcell=(DetId)cell;
    // rof 27.03.09: drop ZS marked and passed digis:
    if (dropZSmarkedPassed_)
      if (i->zsMarkAndPass()) continue;

    const HcalCalibrations& calibrations=conditions->getHcalCalibrations(cell);
    const HcalQIECoder* channelCoder = conditions->getHcalCoder (cell);
    HcalCoderDb coder (*channelCoder, *shape);
       

    //>>> firstSample & samplesToAdd
        if(tsFromDB_) {
      const HcalRecoParam* param_ts = paramTS->getValues(detcell.rawId());
      first = param_ts->firstSample();    
      toadd = param_ts->samplesToAdd();    
    }
        if(toaddMem != toadd) {
      reco_.initPulseCorr(toadd);
          toaddMem = toadd;
    }

    rec->push_back(reco_.reconstruct(*i,first,toadd,coder,calibrations));
      }
      // return result
      e.put(rec);    
    } else if (subdet_==HcalForward) {
      edm::Handle<HFDigiCollection> digi;
      e.getByLabel(inputLabel_,digi);
      
      // create empty output
      std::auto_ptr<HFRecHitCollection> rec(new HFRecHitCollection);
      rec->reserve(digi->size());
      // run the algorithm
      int toaddMem = 0;
      int first = firstSample_;
      int toadd = samplesToAdd_;
      HFDigiCollection::const_iterator i;
      for (i=digi->begin(); i!=digi->end(); i++) {
    HcalDetId cell = i->id();    
    DetId detcell=(DetId)cell;
    // rof 27.03.09: drop ZS marked and passed digis:
    if (dropZSmarkedPassed_)
      if (i->zsMarkAndPass()) continue;
 
    const HcalCalibrations& calibrations=conditions->getHcalCalibrations(cell);
    const HcalQIECoder* channelCoder = conditions->getHcalCoder (cell);
    HcalCoderDb coder (*channelCoder, *shape);

    //>>> firstSample & samplesToAdd
        if(tsFromDB_) {
      const HcalRecoParam* param_ts = paramTS->getValues(detcell.rawId());
      first = param_ts->firstSample();    
      toadd = param_ts->samplesToAdd();    
    }
        if(toaddMem != toadd) {
      reco_.initPulseCorr(toadd);
          toaddMem = toadd;
    }

    rec->push_back(reco_.reconstruct(*i,first,toadd,coder,calibrations));
      }
      // return result
      e.put(rec);     
    } else if (subdet_==HcalOther && subdetOther_==HcalCalibration) {
      edm::Handle<HcalCalibDigiCollection> digi;
      e.getByLabel(inputLabel_,digi);
      
      // create empty output
      std::auto_ptr<HcalCalibRecHitCollection> rec(new HcalCalibRecHitCollection);
      rec->reserve(digi->size());
      // run the algorithm
      int toaddMem = 0;
      int first = firstSample_;
      int toadd = samplesToAdd_;
      HcalCalibDigiCollection::const_iterator i;
      for (i=digi->begin(); i!=digi->end(); i++) {
    HcalCalibDetId cell = i->id();    
    DetId detcell=(DetId)cell;
    // rof 27.03.09: drop ZS marked and passed digis:
    if (dropZSmarkedPassed_)
      if (i->zsMarkAndPass()) continue;

    const HcalCalibrations& calibrations=conditions->getHcalCalibrations(cell);
    const HcalQIECoder* channelCoder = conditions->getHcalCoder (cell);
    HcalCoderDb coder (*channelCoder, *shape);

    //>>> firstSample & samplesToAdd
        if(tsFromDB_) {
      const HcalRecoParam* param_ts = paramTS->getValues(detcell.rawId());
      first = param_ts->firstSample();    
      toadd = param_ts->samplesToAdd();    
    }
        if(toaddMem != toadd) {
      reco_.initPulseCorr(toadd);
          toaddMem = toadd;
    }

    rec->push_back(reco_.reconstruct(*i,first,toadd,coder,calibrations));
      }
      // return result
      e.put(rec);     
    }
  } 
}

Member Data Documentation

DetId::Detector HcalSimpleReconstructor::det_

private

Definition at line 31 of file HcalSimpleReconstructor.h.

Referenced by produce().

bool HcalSimpleReconstructor::dropZSmarkedPassed_

private

Definition at line 36 of file HcalSimpleReconstructor.h.

Referenced by produce().

int HcalSimpleReconstructor::firstSample_

private

Definition at line 40 of file HcalSimpleReconstructor.h.

Referenced by produce().

edm::InputTag HcalSimpleReconstructor::inputLabel_

private

Definition at line 34 of file HcalSimpleReconstructor.h.

Referenced by produce().

HcalRecoParams* HcalSimpleReconstructor::paramTS

private

Definition at line 44 of file HcalSimpleReconstructor.h.

Referenced by beginRun(), endRun(), and produce().

HcalSimpleRecAlgo HcalSimpleReconstructor::reco_

private

Definition at line 30 of file HcalSimpleReconstructor.h.

Referenced by produce().

int HcalSimpleReconstructor::samplesToAdd_

private

Definition at line 41 of file HcalSimpleReconstructor.h.

Referenced by produce().

int HcalSimpleReconstructor::subdet_

private

Definition at line 32 of file HcalSimpleReconstructor.h.

Referenced by HcalSimpleReconstructor(), and produce().

HcalOtherSubdetector HcalSimpleReconstructor::subdetOther_

private

Definition at line 33 of file HcalSimpleReconstructor.h.

Referenced by produce().

bool HcalSimpleReconstructor::tsFromDB_

private

Definition at line 42 of file HcalSimpleReconstructor.h.

Referenced by produce().