CMS 3D CMS Logo

EcalPreshowerRecHitsMaker Class Reference

#include <FastSimulation/CaloRecHitsProducer/interface/EcalPreshowerRecHitsMaker.h>

List of all members.

Public Member Functions

 EcalPreshowerRecHitsMaker (edm::ParameterSet const &p, const RandomEngine *random)
void init (const edm::EventSetup &es)
void loadEcalPreshowerRecHits (edm::Event &iEvent, ESRecHitCollection &esRecHits)
 ~EcalPreshowerRecHitsMaker ()

Private Member Functions

void clean ()
unsigned createVectorsOfCells (const edm::EventSetup &es)
void Fill (uint32_t id, float energy, std::map< uint32_t, std::pair< float, bool > > &myHits, bool signal=true)
void loadPCaloHits (const edm::Event &iEvent)
void noisify ()
void noisifySignal (std::map< uint32_t, std::pair< float, bool > > &theMap)
void noisifySubdet (std::map< uint32_t, std::pair< float, bool > > &theMap, const std::vector< uint32_t > &thecells, unsigned ncells)

Private Attributes

std::map< uint32_t, std::pair
< float, bool > > 		ecalsRecHits_
std::vector< uint32_t > escells_
bool initialized_
edm::InputTag inputCol_
const GaussianTail * myGaussianTailGenerator_
unsigned ncells_
double noise_
double preshowerHotFraction_
const RandomEngine * random_
double threshold_

Detailed Description

Definition at line 21 of file EcalPreshowerRecHitsMaker.h.

Constructor & Destructor Documentation

EcalPreshowerRecHitsMaker::EcalPreshowerRecHitsMaker ( edm::ParameterSet const &  p,
const RandomEngine *  random 
)

Definition at line 19 of file EcalPreshowerRecHitsMaker.cc.

References edm::ParameterSet::getParameter(), initialized_, inputCol_, myGaussianTailGenerator_, noise_, preshowerHotFraction_, random_, funct::sqrt(), and threshold_.

00022   :
00023   random_(myrandom),
00024   myGaussianTailGenerator_(0)
00025 {
00026   edm::ParameterSet RecHitsParameters 
00027     = p.getParameter<edm::ParameterSet>("ECALPreshower");
00028   noise_ = RecHitsParameters.getParameter<double>("Noise");
00029   threshold_ = RecHitsParameters.getParameter<double>("Threshold");
00030   inputCol_=RecHitsParameters.getParameter<edm::InputTag>("MixedSimHits");
00031 
00032   initialized_=false;
00033 
00034   Genfun::Erf myErf;
00035   if(  noise_>0. ) {
00036     preshowerHotFraction_ = 0.5-0.5*myErf(threshold_/noise_/sqrt(2.));
00037     myGaussianTailGenerator_ = new GaussianTail(random_, noise_, threshold_);
00038   } else {
00039     preshowerHotFraction_ =0.;
00040   }
00041 }

EcalPreshowerRecHitsMaker::~EcalPreshowerRecHitsMaker (  ) 

Definition at line 43 of file EcalPreshowerRecHitsMaker.cc.

References initialized_, and myGaussianTailGenerator_.

00044 {
00045   initialized_=false;
00046   delete myGaussianTailGenerator_;
00047 }

Member Function Documentation

void EcalPreshowerRecHitsMaker::clean (  )  [private]

Definition at line 49 of file EcalPreshowerRecHitsMaker.cc.

References ecalsRecHits_.

Referenced by loadPCaloHits().

00050 {
00051 
00052   ecalsRecHits_.clear();
00053 }

unsigned EcalPreshowerRecHitsMaker::createVectorsOfCells ( const edm::EventSetup &  es  )  [private]

Definition at line 106 of file EcalPreshowerRecHitsMaker.cc.

References DetId::Ecal, EcalPreshower, escells_, edm::EventSetup::get(), CaloSubdetectorGeometry::getValidDetIds(), and i.

Referenced by init().

00107 {
00108     edm::ESHandle<CaloGeometry> pG;
00109     es.get<CaloGeometryRecord>().get(pG);     
00110     unsigned total=0;
00111 
00112     escells_.reserve(137728);
00113 
00114     const CaloSubdetectorGeometry* geom=pG->getSubdetectorGeometry(DetId::Ecal,EcalPreshower);  
00115     if(geom==0) return 0;
00116     const std::vector<DetId>& ids(geom->getValidDetIds(DetId::Ecal,EcalPreshower));  
00117     for (std::vector<DetId>::const_iterator i=ids.begin(); i!=ids.end(); i++) 
00118       {
00119         escells_.push_back(i->rawId());
00120         ++total;
00121       }
00122     return escells_.size();
00123 }

void EcalPreshowerRecHitsMaker::Fill ( uint32_t  id,
float  energy,
std::map< uint32_t, std::pair< float, bool > > &  myHits,
bool  signal = true 
) [private]

Definition at line 168 of file EcalPreshowerRecHitsMaker.cc.

References RandomEngine::gaussShoot(), noise_, and random_.

Referenced by loadPCaloHits().

00169 {
00170   // The signal hits are singletons (no need to look into a map)
00171   if(signal)
00172     {
00173       // a new hit is created
00174       // we can give a hint for the insert
00175       // Add the noise at this point. We are sure that it won't be added several times
00176       energy += random_->gaussShoot(0.,noise_);
00177       std::pair<float,bool> hit(energy,false); 
00178       // if it is signal, it is already ordered, so we can give a hint for the 
00179       // insert
00180       if(signal)
00181         myHits.insert(myHits.end(),std::pair<uint32_t,std::pair<float,bool> >(id,hit));
00182     }
00183   else       // In this case,there is a risk of duplication. Need to look into the map
00184     {
00185       std::map<uint32_t,std::pair<float,bool> >::iterator itcheck=myHits.find(id);
00186       if(itcheck==myHits.end())
00187         {
00188           energy += random_->gaussShoot(0.,noise_);
00189           std::pair<float,bool> hit(energy,false);        
00190           myHits.insert(std::pair<uint32_t,std::pair<float,bool> >(id,hit));
00191         }
00192       else
00193         {
00194           itcheck->second.first += energy;
00195         }
00196     }  
00197 }

void EcalPreshowerRecHitsMaker::init ( const edm::EventSetup &  es  ) 

Definition at line 98 of file EcalPreshowerRecHitsMaker.cc.

References createVectorsOfCells(), lat::endl(), initialized_, and ncells_.

Referenced by CaloRecHitsProducer::beginRun().

00099 {
00100   if(initialized_) return;
00101   ncells_=createVectorsOfCells(es);
00102   edm::LogInfo("CaloRecHitsProducer") << "Total number of cells in Preshower " << ncells_ << std::endl;
00103   initialized_=true;
00104 }

void EcalPreshowerRecHitsMaker::loadEcalPreshowerRecHits ( edm::Event &  iEvent,
ESRecHitCollection &  esRecHits 
)

Definition at line 56 of file EcalPreshowerRecHitsMaker.cc.

References ecalsRecHits_, it, loadPCaloHits(), myGaussianTailGenerator_, ncells_, noisify(), edm::SortedCollection< T, SORT >::push_back(), and threshold_.

Referenced by CaloRecHitsProducer::produce().

00057 {
00058   // if no preshower, do nothing
00059 
00060   if(ncells_==0) return;
00061   loadPCaloHits(iEvent);
00062   if( myGaussianTailGenerator_ ) noisify();
00063 
00064   std::map<uint32_t,std::pair<float,bool> >::const_iterator 
00065     it=ecalsRecHits_.begin();
00066   std::map<uint32_t,std::pair<float,bool> >::const_iterator 
00067     itend=ecalsRecHits_.end();
00068 
00069   for(;it!=itend;++it)
00070     {
00071       // check if the hit has been killed 
00072       if(it->second.second) continue;
00073       // check if it is above the threshold
00074       if(it->second.first<threshold_) continue;
00075       ESDetId detid(it->first);
00076       ecalHits.push_back(EcalRecHit(detid,it->second.first,0.)); 
00077     }
00078 }

void EcalPreshowerRecHitsMaker::loadPCaloHits ( const edm::Event &  iEvent  )  [private]

Definition at line 80 of file EcalPreshowerRecHitsMaker.cc.

References clean(), ecalsRecHits_, Fill(), edm::Event::getByLabel(), inputCol_, it, and edm::Handle< T >::product().

Referenced by loadEcalPreshowerRecHits().

00081 {
00082 
00083   clean();
00084 
00085   edm::Handle<CrossingFrame<PCaloHit> > cf;
00086   iEvent.getByLabel(inputCol_,cf);
00087   std::auto_ptr<MixCollection<PCaloHit> > colcalo(new MixCollection<PCaloHit>(cf.product(),std::pair<int,int>(0,0) ));
00088 
00089   MixCollection<PCaloHit>::iterator it=colcalo->begin();
00090   MixCollection<PCaloHit>::iterator itend=colcalo->end();
00091   for(;it!=itend;++it)
00092     {
00093       Fill(it->id(),it->energy(),ecalsRecHits_,it.getTrigger());
00094       //      Fill(it->id(),it->energy(),ecalsRecHits_);
00095     }
00096 }

void EcalPreshowerRecHitsMaker::noisify (  )  [private]

Definition at line 125 of file EcalPreshowerRecHitsMaker.cc.

References ecalsRecHits_, lat::endl(), escells_, ncells_, and noisifySubdet().

Referenced by loadEcalPreshowerRecHits().

00126 {
00127   if(ecalsRecHits_.size()<ncells_) 
00128     {
00129       // Not needed anymore, the noise is added when loading the PCaloHits
00130       // noisifySignal(ecalsRecHits_);
00131       noisifySubdet(ecalsRecHits_,escells_,ncells_);
00132     }
00133   else
00134     edm::LogWarning("CaloRecHitsProducer") << "All HCAL(HB-HE) cells on ! " << std::endl;
00135 }

void EcalPreshowerRecHitsMaker::noisifySignal ( std::map< uint32_t, std::pair< float, bool > > &  theMap  )  [private]

void EcalPreshowerRecHitsMaker::noisifySubdet ( std::map< uint32_t, std::pair< float, bool > > &  theMap,
const std::vector< uint32_t > &  thecells,
unsigned  ncells 
) [private]

Definition at line 138 of file EcalPreshowerRecHitsMaker.cc.

References RandomEngine::flatShoot(), mean(), myGaussianTailGenerator_, RandomEngine::poissonShoot(), preshowerHotFraction_, random_, and GaussianTail::shoot().

Referenced by noisify().

00139 {
00140   // noise won't be injected in cells that contain signal
00141   double mean = (double)(ncells-theMap.size())*preshowerHotFraction_;
00142   unsigned nps = random_->poissonShoot(mean);
00143   
00144   unsigned ncell=0;
00145   unsigned cellindex=0;
00146   uint32_t cellnumber=0;
00147   std::map<uint32_t,std::pair<float,bool> >::const_iterator itcheck;
00148 
00149   while(ncell < nps)
00150     {
00151       cellindex = (unsigned)(random_->flatShoot()*ncells);
00152       cellnumber = thecells[cellindex];
00153       itcheck=theMap.find(cellnumber);
00154       if(itcheck==theMap.end()) // inject only in empty cells
00155         {
00156           std::pair <float,bool> noisehit(myGaussianTailGenerator_->shoot(),
00157                                           false);
00158           theMap.insert(std::pair<uint32_t,std::pair<float,bool> >
00159                         (cellnumber,noisehit));
00160           ++ncell;
00161         }
00162     }
00163   //  edm::LogInfo("CaloRecHitsProducer") << "CaloRecHitsProducer : added noise in "<<  ncell << " HCAL cells "  << std::endl;
00164 }

Member Data Documentation

std::map<uint32_t,std::pair<float,bool> > EcalPreshowerRecHitsMaker::ecalsRecHits_ [private]

Definition at line 54 of file EcalPreshowerRecHitsMaker.h.

Referenced by clean(), loadEcalPreshowerRecHits(), loadPCaloHits(), and noisify().

std::vector<uint32_t> EcalPreshowerRecHitsMaker::escells_ [private]

Definition at line 55 of file EcalPreshowerRecHitsMaker.h.

Referenced by createVectorsOfCells(), and noisify().

bool EcalPreshowerRecHitsMaker::initialized_ [private]

Definition at line 52 of file EcalPreshowerRecHitsMaker.h.

Referenced by EcalPreshowerRecHitsMaker(), init(), and ~EcalPreshowerRecHitsMaker().

edm::InputTag EcalPreshowerRecHitsMaker::inputCol_ [private]

Definition at line 48 of file EcalPreshowerRecHitsMaker.h.

Referenced by EcalPreshowerRecHitsMaker(), and loadPCaloHits().

const GaussianTail* EcalPreshowerRecHitsMaker::myGaussianTailGenerator_ [private]

Definition at line 57 of file EcalPreshowerRecHitsMaker.h.

Referenced by EcalPreshowerRecHitsMaker(), loadEcalPreshowerRecHits(), noisifySubdet(), and ~EcalPreshowerRecHitsMaker().

unsigned EcalPreshowerRecHitsMaker::ncells_ [private]

Definition at line 53 of file EcalPreshowerRecHitsMaker.h.

Referenced by init(), loadEcalPreshowerRecHits(), and noisify().

double EcalPreshowerRecHitsMaker::noise_ [private]

Definition at line 50 of file EcalPreshowerRecHitsMaker.h.

Referenced by EcalPreshowerRecHitsMaker(), and Fill().

double EcalPreshowerRecHitsMaker::preshowerHotFraction_ [private]

Definition at line 51 of file EcalPreshowerRecHitsMaker.h.

Referenced by EcalPreshowerRecHitsMaker(), and noisifySubdet().

const RandomEngine* EcalPreshowerRecHitsMaker::random_ [private]

Definition at line 56 of file EcalPreshowerRecHitsMaker.h.

Referenced by EcalPreshowerRecHitsMaker(), Fill(), and noisifySubdet().

double EcalPreshowerRecHitsMaker::threshold_ [private]

Definition at line 49 of file EcalPreshowerRecHitsMaker.h.

Referenced by EcalPreshowerRecHitsMaker(), and loadEcalPreshowerRecHits().

The documentation for this class was generated from the following files:
Generated on Tue Jun 9 18:19:47 2009 for CMSSW by  doxygen 1.5.4