#include <ReducedESRecHitCollectionProducer.h>

Inheritance diagram for ReducedESRecHitCollectionProducer:

Public Member Functions
void	beginRun (edm::Run &, const edm::EventSetup &)

void	collectIds (const ESDetId strip1, const ESDetId strip2, const int &row=0)

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

	ReducedESRecHitCollectionProducer (const edm::ParameterSet &pset)

virtual	~ReducedESRecHitCollectionProducer ()

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
std::set< DetId >	collectedIds_

const EcalPreshowerGeometry *	geometry_p

edm::InputTag	InputRecHitES_

edm::InputTag	InputSpuerClusterEE_

std::vector< edm::InputTag >	interestingDetIdCollections_

std::string	OutputLabelES_

double	scEtThresh_

CaloSubdetectorTopology *	topology_p

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

Definition at line 27 of file ReducedESRecHitCollectionProducer.h.

Constructor & Destructor Documentation

ReducedESRecHitCollectionProducer::ReducedESRecHitCollectionProducer ( const edm::ParameterSet & pset )

Definition at line 17 of file ReducedESRecHitCollectionProducer.cc.

References edm::ParameterSet::getParameter(), InputRecHitES_, InputSpuerClusterEE_, interestingDetIdCollections_, OutputLabelES_, and scEtThresh_.

                                                                                              :
   geometry_p(0),
   topology_p(0)
 {
 
  scEtThresh_          = ps.getParameter<double>("scEtThreshold");
 
  InputRecHitES_       = ps.getParameter<edm::InputTag>("EcalRecHitCollectionES");
  InputSpuerClusterEE_ = ps.getParameter<edm::InputTag>("EndcapSuperClusterCollection"); 
 
  OutputLabelES_       = ps.getParameter<std::string>("OutputLabel_ES");
  
  interestingDetIdCollections_         = ps.getParameter<std::vector< edm::InputTag> >("interestingDetIds");
  
  produces< EcalRecHitCollection > (OutputLabelES_);
  
 }

ReducedESRecHitCollectionProducer::~ReducedESRecHitCollectionProducer ( )

virtual

Definition at line 35 of file ReducedESRecHitCollectionProducer.cc.

References topology_p.

                                                                       {
   if (topology_p) delete topology_p;
 }

Member Function Documentation

void ReducedESRecHitCollectionProducer::beginRun	(	edm::Run &	,
		const edm::EventSetup &	iSetup
	)

virtual

Reimplemented from edm::EDProducer.

Definition at line 39 of file ReducedESRecHitCollectionProducer.cc.

References DetId::Ecal, EcalPreshower, geometry, geometry_p, edm::EventSetup::get(), and topology_p.

                                                                                      {
   ESHandle<CaloGeometry> geoHandle;
   iSetup.get<CaloGeometryRecord>().get(geoHandle);
   const CaloSubdetectorGeometry *geometry = geoHandle->getSubdetectorGeometry(DetId::Ecal, EcalPreshower);
   geometry_p = dynamic_cast<const EcalPreshowerGeometry *>(geometry);
   if (!geometry_p){
     edm::LogError("WrongGeometry")<<
       "could not cast the subdet geometry to preshower geometry";
   }
   
   if (geometry) topology_p = new EcalPreshowerTopology(geoHandle);
   
 }

void ReducedESRecHitCollectionProducer::collectIds	(	const ESDetId	strip1,
		const ESDetId	strip2,
		const int &	row = `0`
	)

Definition at line 121 of file ReducedESRecHitCollectionProducer.cc.

References collectedIds_, CaloNavigator< T >::east(), CaloNavigator< T >::home(), i, CaloNavigator< T >::north(), CaloNavigator< T >::setHome(), CaloNavigator< T >::south(), topology_p, and CaloNavigator< T >::west().

Referenced by produce().

                                                                                                                   {
 
   //cout<<row<<endl;
   
   map<DetId,const EcalRecHit*>::iterator it;
   map<DetId, int>::iterator itu;
   ESDetId next;
   ESDetId strip1;
   ESDetId strip2;
 
   strip1 = esDetId1;
   strip2 = esDetId2;
 
   EcalPreshowerNavigator theESNav1(strip1, topology_p);
   theESNav1.setHome(strip1);
   
   EcalPreshowerNavigator theESNav2(strip2, topology_p);
   theESNav2.setHome(strip2);
 
   if (row == 1) {
     if (strip1 != ESDetId(0)) strip1 = theESNav1.north();
     if (strip2 != ESDetId(0)) strip2 = theESNav2.east();
   } else if (row == -1) {
     if (strip1 != ESDetId(0)) strip1 = theESNav1.south();
     if (strip2 != ESDetId(0)) strip2 = theESNav2.west();
   }
 
   // Plane 1 
   if (strip1 == ESDetId(0)) {
   } else {
     collectedIds_.insert(strip1);
     //cout<<"center : "<<strip1<<endl;
     // east road 
     for (int i=0; i<15; ++i) {
       next = theESNav1.east();
       //cout<<"east : "<<i<<" "<<next<<endl;
       if (next != ESDetId(0)) {
     collectedIds_.insert(next);
       } else {
     break;
       }
     }
 
     // west road 
     theESNav1.setHome(strip1);
     theESNav1.home();
     for (int i=0; i<15; ++i) {
       next = theESNav1.west();
       //cout<<"west : "<<i<<" "<<next<<endl;
       if (next != ESDetId(0)) {
     collectedIds_.insert(next);
       } else {
     break;
       }
     }
 
   }
 
   if (strip2 == ESDetId(0)) {
   } else {
     collectedIds_.insert(strip2);
     //cout<<"center : "<<strip2<<endl;
     // north road 
     for (int i=0; i<15; ++i) {
       next = theESNav2.north();
       //cout<<"north : "<<i<<" "<<next<<endl;
       if (next != ESDetId(0)) {
     collectedIds_.insert(next);
       } else {
     break;
       } 
     }
 
     // south road 
     theESNav2.setHome(strip2);
     theESNav2.home();
     for (int i=0; i<15; ++i) {
       next = theESNav2.south();
       //cout<<"south : "<<i<<" "<<next<<endl;
       if (next != ESDetId(0)) {
     collectedIds_.insert(next);
       } else {
     break;
       }
     }
   }
 }

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

virtual

Implements edm::EDProducer.

Definition at line 53 of file ReducedESRecHitCollectionProducer.cc.

References collectedIds_, collectIds(), geometry_p, edm::Event::getByLabel(), EcalPreshowerGeometry::getClosestCellInPlane(), InputRecHitES_, InputSpuerClusterEE_, interestingDetIdCollections_, edm::HandleBase::isValid(), convertSQLitetoXML_cfg::output, OutputLabelES_, point, edm::Handle< T >::product(), edm::Event::put(), scEtThresh_, and matplotRender::t.

                                                                                         {
 
 
   edm::Handle<ESRecHitCollection> ESRecHits_;
   e.getByLabel(InputRecHitES_, ESRecHits_);
   
   std::auto_ptr<EcalRecHitCollection> output(new EcalRecHitCollection);
 
   edm::Handle<reco::SuperClusterCollection> pEndcapSuperClusters;
   e.getByLabel(InputSpuerClusterEE_, pEndcapSuperClusters);
   {
     const reco::SuperClusterCollection* eeSuperClusters = pEndcapSuperClusters.product();
     
     for (reco::SuperClusterCollection::const_iterator isc = eeSuperClusters->begin(); isc != eeSuperClusters->end(); ++isc) {
 
       if (isc->energy() < scEtThresh_) continue;
       if (fabs(isc->eta()) < 1.65 || fabs(isc->eta()) > 2.6) continue;
       //cout<<"SC energy : "<<isc->energy()<<" "<<isc->eta()<<endl;
 
       //Int_t nBC = 0;
       reco::CaloCluster_iterator ibc = isc->clustersBegin();
       for ( ; ibc != isc->clustersEnd(); ++ibc ) {  
 
     //cout<<"BC : "<<nBC<<endl;
 
     const GlobalPoint point((*ibc)->x(),(*ibc)->y(),(*ibc)->z());
 
     ESDetId esId1 = geometry_p->getClosestCellInPlane(point, 1);
     ESDetId esId2 = geometry_p->getClosestCellInPlane(point, 2);
     
     collectIds(esId1, esId2, 0);
     collectIds(esId1, esId2, 1);
     collectIds(esId1, esId2, -1);
 
     //nBC++;
       }
       
     }
     
   }
 
 
   edm::Handle<DetIdCollection > detId;
   for( unsigned int t = 0; t < interestingDetIdCollections_.size(); ++t )
     {
       e.getByLabel(interestingDetIdCollections_[t],detId);
       if (!detId.isValid()){
     edm::LogError("MissingInput")<<"the collection of interesting detIds:"<<interestingDetIdCollections_[t]<<" is not found.";
         continue;
       }
       collectedIds_.insert(detId->begin(),detId->end());
     }
 
 
   output->reserve( collectedIds_.size());
   EcalRecHitCollection::const_iterator it;
   for (it = ESRecHits_->begin(); it != ESRecHits_->end(); ++it) {
     if (it->recoFlag()==1 || it->recoFlag()==14 || (it->recoFlag()<=10 && it->recoFlag()>=5)) continue;
     if (collectedIds_.find(it->id())!=collectedIds_.end()){
       output->push_back(*it);
     }
   }
   collectedIds_.clear();
 
   e.put(output, OutputLabelES_);
 
 }