HGCDigitizer Class Reference

#include <HGCDigitizer.h>

Public Member Functions
void	accumulate (edm::Event const &e, edm::EventSetup const &c, CLHEP::HepRandomEngine *)
	handle SimHit accumulation More...
void	accumulate (PileUpEventPrincipal const &e, edm::EventSetup const &c, CLHEP::HepRandomEngine *)
void	accumulate (edm::Handle< edm::PCaloHitContainer > const &hits, int bxCrossing, const edm::ESHandle< HGCalGeometry > &geom, CLHEP::HepRandomEngine *)
void	beginRun (const edm::EventSetup &es)
	actions at the start/end of run More...
std::string	digiCollection ()
void	endRun ()
void	finalizeEvent (edm::Event &e, edm::EventSetup const &c, CLHEP::HepRandomEngine *)
	HGCDigitizer (const edm::ParameterSet &ps, edm::ConsumesCollector &iC)
void	initializeEvent (edm::Event const &e, edm::EventSetup const &c)
	actions at the start/end of event More...
bool	producesEEDigis ()
bool	producesHEbackDigis ()
bool	producesHEfrontDigis ()
	~HGCDigitizer ()

Private Member Functions
void	resetSimHitDataAccumulator ()

Private Attributes
int	bxTime_
bool	checkValidDetIds_
std::string	digiCollection_
int	digitizationType_
std::string	hitCollection_
int	maxSimHitsAccTime_
ForwardSubdetector	mySubDet_
float	refSpeed_
std::unique_ptr < HGCSimHitDataAccumulator >	simHitAccumulator_
std::unique_ptr< HGCEEDigitizer >	theHGCEEDigitizer_
std::unique_ptr < HGCHEbackDigitizer >	theHGCHEbackDigitizer_
std::unique_ptr < HGCHEfrontDigitizer >	theHGCHEfrontDigitizer_
float	tofDelay_
bool	useAllChannels_
int	verbosity_

Detailed Description

Definition at line 34 of file HGCDigitizer.h.

Constructor & Destructor Documentation

HGCDigitizer::HGCDigitizer ( const edm::ParameterSet &  ps, edm::ConsumesCollector &  iC  )

explicit

Definition at line 20 of file HGCDigitizer.cc.

References bxTime_, edm::ConsumesCollector::consumes(), digiCollection_, digitizationType_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), HGCEE, HGCHEB, HGCHEF, hitCollection_, HLT_25ns14e33_v1_cff::InputTag, maxSimHitsAccTime_, mySubDet_, AlCaHLTBitMon_QueryRunRegistry::string, theHGCEEDigitizer_, theHGCHEbackDigitizer_, theHGCHEfrontDigitizer_, tofDelay_, useAllChannels_, and verbosity_.

                                                                               : 
  checkValidDetIds_(true),
  simHitAccumulator_( new HGCSimHitDataAccumulator ),
  mySubDet_(ForwardSubdetector::ForwardEmpty),
  refSpeed_(0.1*CLHEP::c_light) { //[CLHEP::c_light]=mm/ns convert to cm/ns
  
  //configure from cfg
  hitCollection_     = ps.getParameter< std::string >("hitCollection");
  digiCollection_    = ps.getParameter< std::string >("digiCollection");
  maxSimHitsAccTime_ = ps.getParameter< uint32_t >("maxSimHitsAccTime");
  bxTime_            = ps.getParameter< int32_t >("bxTime");
  digitizationType_  = ps.getParameter< uint32_t >("digitizationType");
  useAllChannels_    = ps.getParameter< bool >("useAllChannels");
  verbosity_         = ps.getUntrackedParameter< int32_t >("verbosity",0);
  tofDelay_          = ps.getParameter< double >("tofDelay");  

  iC.consumes<std::vector<PCaloHit> >(edm::InputTag("g4SimHits",hitCollection_));

  if(hitCollection_.find("HitsEE")!=std::string::npos) { 
    mySubDet_=ForwardSubdetector::HGCEE;  
    theHGCEEDigitizer_=std::unique_ptr<HGCEEDigitizer>(new HGCEEDigitizer(ps) ); 
  }
  if(hitCollection_.find("HitsHEfront")!=std::string::npos) { 
    mySubDet_=ForwardSubdetector::HGCHEF;
    theHGCHEfrontDigitizer_=std::unique_ptr<HGCHEfrontDigitizer>(new HGCHEfrontDigitizer(ps) );
  }
  if(hitCollection_.find("HitsHEback")!=std::string::npos) { 
    mySubDet_=ForwardSubdetector::HGCHEB;
    theHGCHEbackDigitizer_=std::unique_ptr<HGCHEbackDigitizer>(new HGCHEbackDigitizer(ps) );
  }
}

HGCDigitizer::~HGCDigitizer ( )

inline

Definition at line 39 of file HGCDigitizer.h.

{ }

Member Function Documentation

void HGCDigitizer::accumulate	(	edm::Event const &	e,
		edm::EventSetup const &	c,
		CLHEP::HepRandomEngine *	engine
	)

handle SimHit accumulation

Definition at line 80 of file HGCDigitizer.cc.

References relativeConstraints::geom, edm::EventSetup::get(), edm::Event::getByLabel(), hitCollection_, edm::HandleBase::isValid(), producesEEDigis(), producesHEbackDigis(), and producesHEfrontDigis().

Referenced by accumulate().

                                                                                                              {

  //get inputs
  edm::Handle<edm::PCaloHitContainer> hits;
  e.getByLabel(edm::InputTag("g4SimHits",hitCollection_),hits); 
  if( !hits.isValid() ){
    edm::LogError("HGCDigitizer") << " @ accumulate : can't find " << hitCollection_ << " collection of g4SimHits";
    return;
  }

  //get geometry
  edm::ESHandle<HGCalGeometry> geom;
  if( producesEEDigis() )      eventSetup.get<IdealGeometryRecord>().get("HGCalEESensitive"            , geom);
  if( producesHEfrontDigis() ) eventSetup.get<IdealGeometryRecord>().get("HGCalHESiliconSensitive"     , geom);
  if( producesHEbackDigis() )  eventSetup.get<IdealGeometryRecord>().get("HGCalHEScintillatorSensitive", geom);

  //accumulate in-time the main event
  accumulate(hits, 0, geom, engine);
}

void HGCDigitizer::accumulate	(	PileUpEventPrincipal const &	e,
		edm::EventSetup const &	c,
		CLHEP::HepRandomEngine *	engine
	)

Definition at line 101 of file HGCDigitizer.cc.

References accumulate(), PileUpEventPrincipal::bunchCrossing(), relativeConstraints::geom, edm::EventSetup::get(), PileUpEventPrincipal::getByLabel(), hitCollection_, edm::HandleBase::isValid(), producesEEDigis(), producesHEbackDigis(), and producesHEfrontDigis().

                                                                                                                          {

  //get inputs
  edm::Handle<edm::PCaloHitContainer> hits;
  e.getByLabel(edm::InputTag("g4SimHits",hitCollection_),hits); 
  if( !hits.isValid() ){
    edm::LogError("HGCDigitizer") << " @ accumulate : can't find " << hitCollection_ << " collection of g4SimHits";
    return;
  }

  //get geometry
  edm::ESHandle<HGCalGeometry> geom;
  if( producesEEDigis() )      eventSetup.get<IdealGeometryRecord>().get("HGCalEESensitive"            , geom);
  if( producesHEfrontDigis() ) eventSetup.get<IdealGeometryRecord>().get("HGCalHESiliconSensitive"     , geom);
  if( producesHEbackDigis() )  eventSetup.get<IdealGeometryRecord>().get("HGCalHEScintillatorSensitive", geom);

  //accumulate for the simulated bunch crossing
  accumulate(hits, e.bunchCrossing(), geom, engine);
}

void HGCDigitizer::accumulate	(	edm::Handle< edm::PCaloHitContainer > const &	hits,
		int	bxCrossing,
		const edm::ESHandle< HGCalGeometry > &	geom,
		CLHEP::HepRandomEngine *	engine
	)

Definition at line 122 of file HGCDigitizer.cc.

References bxTime_, checkValidDetIds_, gather_cfg::cout, HGCalTopology::dddConstants(), HGCalTopology::detectorType(), hitCollection_, edm::ESHandleBase::isValid(), mySubDet_, producesEEDigis(), refSpeed_, simHitAccumulator_, HGCalDDDConstants::simToReco(), tofDelay_, useAllChannels_, and verbosity_.

                                                                                                                                                              {
  if(!geom.isValid()) return;
  const HGCalTopology &topo=geom->topology();
  const HGCalDDDConstants &dddConst=topo.dddConstants();

  for(edm::PCaloHitContainer::const_iterator hit_it = hits->begin(); hit_it != hits->end(); ++hit_it) {
    HGCalDetId simId( hit_it->id() );

    //skip this hit if after ganging it is not valid
    int layer(simId.layer()), cell(simId.cell());
    std::pair<int,int> recoLayerCell=dddConst.simToReco(cell,layer,topo.detectorType());
    cell  = recoLayerCell.first;
    layer = recoLayerCell.second;
    if(layer<0 || cell<0) continue;
     
    //assign the RECO DetId
    DetId id( producesEEDigis() ?
          (uint32_t)HGCEEDetId(mySubDet_,simId.zside(),layer,simId.sector(),simId.subsector(),cell):
          (uint32_t)HGCHEDetId(mySubDet_,simId.zside(),layer,simId.sector(),simId.subsector(),cell) );

    if (verbosity_>0) {
      if (producesEEDigis())
    std::cout << "HGCDigitizer: i/p " << simId << " o/p " << HGCEEDetId(id) << std::endl;
      else
    std::cout << "HGCDigitizer: i/p " << simId << " o/p " << HGCHEDetId(id) << std::endl;
    }

    //distance to the center of the detector
    float dist2center( geom->getPosition(id).mag() );

    //hit time: [time()]=ns  [centerDist]=cm [refSpeed_]=cm/ns + delay by 1ns
    //accumulate in 6 buckets of 25ns (4 pre-samples, 1 in-time, 1 post-sample)
    float tof(hit_it->time()-dist2center/refSpeed_+tofDelay_);
    int itime=floor( tof/bxTime_ ) ;
      
    //no need to add bx crossing - tof comes already corrected from the mixing module
    //itime += bxCrossing;
    itime += 4;
      
    if(itime<0 || itime>5) continue; 
      
    //energy deposited 
    HGCSimEn_t ien( hit_it->energy() );
      
    //check if already existing (perhaps could remove this in the future - 2nd event should have all defined)
    HGCSimHitDataAccumulator::iterator simHitIt=simHitAccumulator_->find(id);
    if(simHitIt==simHitAccumulator_->end()) {
      HGCSimHitData baseData;
      baseData.fill(0.);
      simHitAccumulator_->insert( std::make_pair(id,baseData) );
      simHitIt=simHitAccumulator_->find(id);
    }
      
    //check if time index is ok and store energy
    if(itime >= (int)simHitIt->second.size() ) continue;
    (simHitIt->second)[itime] += ien;
  }
  
  //add base data for noise simulation
  if(!checkValidDetIds_) return;
  if(!geom.isValid()) return;
  HGCSimHitData baseData;
  baseData.fill(0.);
  const std::vector<DetId> &validIds=geom->getValidDetIds(); 
  int nadded(0);
  if (useAllChannels_) {
    for(std::vector<DetId>::const_iterator it=validIds.begin(); it!=validIds.end(); it++) {
      uint32_t id(it->rawId());
      if(simHitAccumulator_->find(id)!=simHitAccumulator_->end()) continue;
      simHitAccumulator_->insert( std::make_pair(id,baseData) );
      nadded++;
    }
  }
  if (verbosity_ > 0) 
    std::cout << "HGCDigitizer:Added " << nadded << ":" << validIds.size() 
          << " detIds without " << hitCollection_ 
          << " in first event processed" << std::endl;
  checkValidDetIds_=false;
}

void HGCDigitizer::beginRun

(

const edm::EventSetup &

es

)

actions at the start/end of run

Definition at line 203 of file HGCDigitizer.cc.

{ }

std::string HGCDigitizer::digiCollection ( )

inline

Definition at line 59 of file HGCDigitizer.h.

References digiCollection_.

Referenced by finalizeEvent().

{ return digiCollection_; }

void HGCDigitizer::endRun

(

void

)

Definition at line 206 of file HGCDigitizer.cc.

{ }

void HGCDigitizer::finalizeEvent	(	edm::Event &	e,
		edm::EventSetup const &	c,
		CLHEP::HepRandomEngine *	engine
	)

Definition at line 58 of file HGCDigitizer.cc.

References digiCollection(), digitizationType_, producesEEDigis(), producesHEbackDigis(), producesHEfrontDigis(), edm::Event::put(), simHitAccumulator_, theHGCEEDigitizer_, theHGCHEbackDigitizer_, and theHGCHEfrontDigitizer_.

                                                                                                   {
  if( producesEEDigis() ) {
    std::auto_ptr<HGCEEDigiCollection> digiResult(new HGCEEDigiCollection() );
    theHGCEEDigitizer_->run(digiResult,*simHitAccumulator_,digitizationType_,engine);
    edm::LogInfo("HGCDigitizer") << " @ finalize event - produced " << digiResult->size() <<  " EE hits";
    e.put(digiResult,digiCollection());
  }
  if( producesHEfrontDigis()) {
    std::auto_ptr<HGCHEDigiCollection> digiResult(new HGCHEDigiCollection() );
    theHGCHEfrontDigitizer_->run(digiResult,*simHitAccumulator_,digitizationType_,engine);
    edm::LogInfo("HGCDigitizer") << " @ finalize event - produced " << digiResult->size() <<  " HE front hits";
    e.put(digiResult,digiCollection());
  }
  if( producesHEbackDigis() )  {
    std::auto_ptr<HGCHEDigiCollection> digiResult(new HGCHEDigiCollection() );
    theHGCHEbackDigitizer_->run(digiResult,*simHitAccumulator_,digitizationType_,engine);
    edm::LogInfo("HGCDigitizer") << " @ finalize event - produced " << digiResult->size() <<  " HE back hits";
    e.put(digiResult,digiCollection());
  }
}

void HGCDigitizer::initializeEvent	(	edm::Event const &	e,
		edm::EventSetup const &	c
	)

actions at the start/end of event

Definition at line 53 of file HGCDigitizer.cc.

References resetSimHitDataAccumulator().

                                                                             {
  resetSimHitDataAccumulator(); 
}

bool HGCDigitizer::producesEEDigis ( )

inline

Definition at line 56 of file HGCDigitizer.h.

References HGCEE, and mySubDet_.

Referenced by accumulate(), and finalizeEvent().

{ return (mySubDet_==ForwardSubdetector::HGCEE);  }

bool HGCDigitizer::producesHEbackDigis ( )

inline

Definition at line 58 of file HGCDigitizer.h.

References HGCHEB, and mySubDet_.

Referenced by accumulate(), and finalizeEvent().

{ return (mySubDet_==ForwardSubdetector::HGCHEB); }

bool HGCDigitizer::producesHEfrontDigis ( )

inline

Definition at line 57 of file HGCDigitizer.h.

References HGCHEF, and mySubDet_.

Referenced by accumulate(), and finalizeEvent().

{ return (mySubDet_==ForwardSubdetector::HGCHEF); }

void HGCDigitizer::resetSimHitDataAccumulator ( )

private

Definition at line 209 of file HGCDigitizer.cc.

References simHitAccumulator_.

Referenced by initializeEvent().

                                              {
  for( HGCSimHitDataAccumulator::iterator it = simHitAccumulator_->begin(); it!=simHitAccumulator_->end(); it++)  it->second.fill(0.);
}

Member Data Documentation

int HGCDigitizer::bxTime_

private

Definition at line 80 of file HGCDigitizer.h.

Referenced by accumulate(), and HGCDigitizer().

bool HGCDigitizer::checkValidDetIds_

private

Definition at line 70 of file HGCDigitizer.h.

Referenced by accumulate().

std::string HGCDigitizer::digiCollection_

private

Definition at line 73 of file HGCDigitizer.h.

Referenced by digiCollection(), and HGCDigitizer().

int HGCDigitizer::digitizationType_

private

Definition at line 76 of file HGCDigitizer.h.

Referenced by finalizeEvent(), and HGCDigitizer().

std::string HGCDigitizer::hitCollection_

private

Definition at line 73 of file HGCDigitizer.h.

Referenced by accumulate(), and HGCDigitizer().

int HGCDigitizer::maxSimHitsAccTime_

private

Definition at line 79 of file HGCDigitizer.h.

Referenced by HGCDigitizer().

ForwardSubdetector HGCDigitizer::mySubDet_

private

Definition at line 90 of file HGCDigitizer.h.

Referenced by accumulate(), HGCDigitizer(), producesEEDigis(), producesHEbackDigis(), and producesHEfrontDigis().

float HGCDigitizer::refSpeed_

private

Definition at line 97 of file HGCDigitizer.h.

Referenced by accumulate().

std::unique_ptr<HGCSimHitDataAccumulator> HGCDigitizer::simHitAccumulator_

private

Definition at line 81 of file HGCDigitizer.h.

Referenced by accumulate(), finalizeEvent(), and resetSimHitDataAccumulator().

std::unique_ptr<HGCEEDigitizer> HGCDigitizer::theHGCEEDigitizer_

private

Definition at line 85 of file HGCDigitizer.h.

Referenced by finalizeEvent(), and HGCDigitizer().

std::unique_ptr<HGCHEbackDigitizer> HGCDigitizer::theHGCHEbackDigitizer_

private

Definition at line 86 of file HGCDigitizer.h.

Referenced by finalizeEvent(), and HGCDigitizer().

std::unique_ptr<HGCHEfrontDigitizer> HGCDigitizer::theHGCHEfrontDigitizer_

private

Definition at line 87 of file HGCDigitizer.h.

Referenced by finalizeEvent(), and HGCDigitizer().

float HGCDigitizer::tofDelay_

private

Definition at line 100 of file HGCDigitizer.h.

Referenced by accumulate(), and HGCDigitizer().

bool HGCDigitizer::useAllChannels_

private

Definition at line 93 of file HGCDigitizer.h.

Referenced by accumulate(), and HGCDigitizer().

int HGCDigitizer::verbosity_

private

Definition at line 94 of file HGCDigitizer.h.

Referenced by accumulate(), and HGCDigitizer().