---

HcalDummyHitProducer Class Reference

#include <SimCalorimetry/HcalTrigPrimProducers/src/HcalDummyHitProducer.h>

Inheritance diagram for HcalDummyHitProducer:

List of all members.

Public Types
typedef std::multimap < HcalTrigTowerDetId, HcalDetId >	tid2cid
Public Member Functions
	HcalDummyHitProducer (const edm::ParameterSet &ps)
virtual void	produce (edm::Event &e, const edm::EventSetup &c)
	Produces the EDM products,.
virtual	~HcalDummyHitProducer ()
Private Attributes
double	energyEM
double	energyHad
double	step_size
HcalTrigTowerGeometry	theTrigTowerGeometry
tid2cid	Tower_map

---

Detailed Description

Definition at line 17 of file HcalDummyHitProducer.h.

---

Member Typedef Documentation

typedef std::multimap<HcalTrigTowerDetId, HcalDetId> HcalDummyHitProducer::tid2cid

Definition at line 21 of file HcalDummyHitProducer.h.

---

Constructor & Destructor Documentation

HcalDummyHitProducer::HcalDummyHitProducer

(

const edm::ParameterSet &

ps

)

[explicit]

Definition at line 17 of file HcalDummyHitProducer.cc.

References energyEM, energyHad, edm::ParameterSet::getParameter(), and step_size.

{
  energyEM = 0;
  produces<edm::PCaloHitContainer>("HcalHits");
  step_size = ps.getParameter<double>("step_size");
  energyHad = ps.getParameter<double>("energy_start");
}

HcalDummyHitProducer::~HcalDummyHitProducer

(

)

[virtual]

Definition at line 25 of file HcalDummyHitProducer.cc.

{}

---

Member Function Documentation

void HcalDummyHitProducer::produce	(	edm::Event &	e,
		const edm::EventSetup &	c
	)			[virtual]

Produces the EDM products,.

Implements edm::EDProducer.

Definition at line 78 of file HcalDummyHitProducer.cc.

References calib, HcalDetId::depth(), energyHad, edm::EventSetup::get(), DetId::Hcal, HcalBarrel, HcalEndcap, HcalForward, int, n, edm::Event::put(), DetId::rawId(), HLT_VtxMuL3::result, sample_factor(), step_size, theTrigTowerGeometry, Tower_map, and HcalTrigTowerGeometry::towerIds().

                                                                               {
  // get the appropriate gains, noises, & widths for this event
  //edm::ESHandle<HcalDbService> conditions;
  //eventSetup.get<HcalDbRecord>().get(conditions);
  //theParameters->setDbService(conditions.product());

  std::auto_ptr<edm::PCaloHitContainer> result(new edm::PCaloHitContainer);
  // get the correct geometry
  edm::ESHandle<CaloGeometry> geometry;
  eventSetup.get<CaloGeometryRecord>().get(geometry);
  vector<DetId> hbCells =  geometry->getValidDetIds(DetId::Hcal, HcalBarrel);
  vector<DetId> heCells =  geometry->getValidDetIds(DetId::Hcal, HcalEndcap);
  vector<DetId> hfCells =  geometry->getValidDetIds(DetId::Hcal, HcalForward);

  vector<HcalTrigTowerDetId> towerids;
  Tower_map.clear();
  for(vector<DetId>::const_iterator itr = hbCells.begin(); itr != hbCells.end(); ++itr)
    {
      HcalDetId barrelId(*itr);
      towerids = theTrigTowerGeometry.towerIds(barrelId);
      for(unsigned int n = 0; n < towerids.size(); ++n)
        {
          Tower_map.insert(tid2cid::value_type(towerids[n],barrelId));
        }
    }
  for(vector<DetId>::const_iterator itr = heCells.begin(); itr != heCells.end(); ++itr)
    {      
      HcalDetId endcapId(*itr);
      towerids = theTrigTowerGeometry.towerIds(endcapId);
      for(unsigned int n = 0; n < towerids.size(); ++n)
      {
        Tower_map.insert(tid2cid::value_type(towerids[n],endcapId));
      }
    }
  for(vector<DetId>::const_iterator itr = hfCells.begin(); itr != hfCells.end(); ++itr)
    {
      HcalDetId forwardId(*itr);
      towerids = theTrigTowerGeometry.towerIds(forwardId);
      for(unsigned int n = 0; n < towerids.size(); ++n)
      {
        Tower_map.insert(tid2cid::value_type(towerids[n],forwardId));
      }
    }

  double time;
  double calib;
  int ncells;
  int ntowers;
  double cell_energyEM, cell_energyHad;

  for(vector<DetId>::const_iterator itr = hbCells.begin(); itr != hbCells.end(); ++itr)
    {
      HcalDetId barrelId(*itr);
      //const CaloSimParameters & parameters = theParameterMap->simParameters(*itr);
      //calib = parameters.simHitToPhotoelectrons()
      //    * parameters.photoelectronsToAnalog()
      //    * theParameters->fCtoGeV(*itr);
      //time = 8.4 - parameters.timePhase(); 
      calib = sample_factor(barrelId);
      time = 8.4;
      towerids = theTrigTowerGeometry.towerIds(barrelId);
      ncells = Tower_map.count(towerids[0]);
      ntowers = towerids.size();
      cell_energyEM = 0;
      //cell_energyHad= (energyHad*ntowers)/(ncells*calib);
      if(barrelId.depth()==1)
        {
          cell_energyHad= (energyHad*ntowers)/(calib);
        }
      else
        {
          cell_energyHad = 0;
        }
      PCaloHit barrelHit(barrelId.rawId(), cell_energyEM, cell_energyHad, time, 0);
      result->push_back(barrelHit);
    }
  for(vector<DetId>::const_iterator itr = heCells.begin(); itr != heCells.end(); ++itr)
    {
      HcalDetId endcapId(*itr);
      //const CaloSimParameters & parameters = theParameterMap->simParameters(*itr);
      //calib = parameters.simHitToPhotoelectrons()
      //* parameters.photoelectronsToAnalog()
      //* theParameters->fCtoGeV(*itr);
      //time = 13.0 - parameters.timePhase();
      calib = sample_factor(endcapId);
      time = 13.0;
      towerids = theTrigTowerGeometry.towerIds(endcapId);
      ncells = Tower_map.count(towerids[0]);
      ntowers = towerids.size();
      cell_energyEM =  0;
      //cell_energyHad= (energyHad*ntowers)/(ncells*calib);
      if(endcapId.depth()==1)
        {
          cell_energyHad= (energyHad*ntowers)/(calib);
        }
      else
        {
          cell_energyHad = 0;
        }
      PCaloHit endcapHit(endcapId.rawId(), cell_energyEM, cell_energyHad, time, 0);
      result->push_back(endcapHit);
    }
  for(vector<DetId>::const_iterator itr = hfCells.begin(); itr != hfCells.end(); ++itr)
    {
      HcalDetId forwardId(*itr);
      //const CaloSimParameters & parameters = theParameterMap->simParameters(*itr);
      //calib = parameters.simHitToPhotoelectrons()
      //* parameters.photoelectronsToAnalog()
      //* theParameters->fCtoGeV(*itr);
      //time = 37.0 - parameters.timePhase();
      calib = sample_factor(forwardId);
      time = 37.0;
      towerids = theTrigTowerGeometry.towerIds(forwardId);
      ncells = Tower_map.count(towerids[0]);
      ntowers = towerids.size();
      cell_energyEM = 0;
      //cell_energyHad= (energyHad*ntowers)/(ncells*calib);
      if(forwardId.depth()==1)
        {
          cell_energyHad= (int)((energyHad)/(calib));
        }
      else
        {
          cell_energyHad = 0;
        }

      PCaloHit forwardHit(forwardId.rawId(), cell_energyEM, cell_energyHad, time, 0);
      result->push_back(forwardHit);
    }
  e.put(result,"HcalHits");
  energyHad += step_size;
}

---

Member Data Documentation

double HcalDummyHitProducer::energyEM [private]

Definition at line 31 of file HcalDummyHitProducer.h.

Referenced by HcalDummyHitProducer().

double HcalDummyHitProducer::energyHad [private]

Definition at line 32 of file HcalDummyHitProducer.h.

Referenced by HcalDummyHitProducer(), and produce().

double HcalDummyHitProducer::step_size [private]

Definition at line 33 of file HcalDummyHitProducer.h.

Referenced by HcalDummyHitProducer(), and produce().

HcalTrigTowerGeometry HcalDummyHitProducer::theTrigTowerGeometry [private]

Definition at line 30 of file HcalDummyHitProducer.h.

Referenced by produce().

tid2cid HcalDummyHitProducer::Tower_map [private]

Definition at line 29 of file HcalDummyHitProducer.h.

Referenced by produce().

---

The documentation for this class was generated from the following files:

• SimCalorimetry/HcalTrigPrimProducers/src/HcalDummyHitProducer.h
• SimCalorimetry/HcalTrigPrimProducers/src/HcalDummyHitProducer.cc

---