#include <EcalDigiSelector.h>

Inheritance diagram for EcalDigiSelector:

Public Member Functions
	EcalDigiSelector (const edm::ParameterSet &ps)

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

Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default

Public Member Functions inherited from edm::stream::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

virtual	~EDProducerBase ()

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

virtual	~ProducerBase ()

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

	EDConsumerBase ()

ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (std::string const &iProcessName, std::string const &iModuleLabel, bool iPrint, std::vector< char const * > &oModuleLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Private Attributes
edm::EDGetTokenT < reco::SuperClusterCollection >	barrelSuperClusterProducer_

double	cluster_pt_thresh_

edm::EDGetTokenT < EBDigiCollection >	EcalEBDigiToken_

edm::EDGetTokenT < EcalRecHitCollection >	EcalEBRecHitToken_

edm::EDGetTokenT < EEDigiCollection >	EcalEEDigiToken_

edm::EDGetTokenT < EcalRecHitCollection >	EcalEERecHitToken_

edm::EDGetTokenT < reco::SuperClusterCollection >	endcapSuperClusterProducer_

int	nclus_sel_

std::string	selectedEcalEBDigiCollection_

std::string	selectedEcalEEDigiCollection_

double	single_cluster_thresh_

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
typedef CacheContexts< T...>	CacheTypes

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T...>	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT < LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	RunSummaryCache

Public Types inherited from edm::stream::EDProducerBase
typedef EDProducerAdaptorBase	ModuleType

Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::stream::EDProducerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Definition at line 20 of file EcalDigiSelector.h.

Constructor & Destructor Documentation

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

Definition at line 23 of file EcalDigiSelector.cc.

References edm::ParameterSet::getParameter(), and AlCaHLTBitMon_QueryRunRegistry::string.

 {
   selectedEcalEBDigiCollection_ = ps.getParameter<std::string>("selectedEcalEBDigiCollection");
   selectedEcalEEDigiCollection_ = ps.getParameter<std::string>("selectedEcalEEDigiCollection");
 
   barrelSuperClusterProducer_ = 
       consumes<reco::SuperClusterCollection>(ps.getParameter<edm::InputTag>("barrelSuperClusterProducer"));
   endcapSuperClusterProducer_ = 
       consumes<reco::SuperClusterCollection>(ps.getParameter<edm::InputTag>("endcapSuperClusterProducer"));
 
   
   EcalEBRecHitToken_ = 
       consumes<EcalRecHitCollection>(ps.getParameter<edm::InputTag>("EcalEBRecHitTag"));
   EcalEERecHitToken_ = 
       consumes<EcalRecHitCollection>(ps.getParameter<edm::InputTag>("EcalEERecHitTag"));
 
   EcalEBDigiToken_ = 
       consumes<EBDigiCollection>(ps.getParameter<edm::InputTag>("EcalEBDigiTag"));
   EcalEEDigiToken_ = 
       consumes<EEDigiCollection>(ps.getParameter<edm::InputTag>("EcalEEDigiTag"));
   
   cluster_pt_thresh_ = ps.getParameter<double>("cluster_pt_thresh");
   single_cluster_thresh_ = ps.getParameter<double>("single_cluster_thresh");
 
   nclus_sel_ = ps.getParameter<int>("nclus_sel");
   produces<EBDigiCollection>(selectedEcalEBDigiCollection_);
   produces<EEDigiCollection>(selectedEcalEEDigiCollection_);
 
 }

Member Function Documentation

void EcalDigiSelector::produce	(	edm::Event &	evt,
		const edm::EventSetup &	es
	)

virtual

Implements edm::stream::EDProducerBase.

Definition at line 53 of file EcalDigiSelector.cc.

References edm::DataFrameContainer::begin(), edm::DataFrameContainer::end(), reco::CaloCluster::energy(), reco::CaloCluster::eta(), create_public_lumi_plots::exp, edm::EventSetup::get(), edm::Ptr< T >::get(), edm::Event::getByToken(), cmsHarvester::loop, edm::Handle< T >::product(), edm::ESHandle< class >::product(), edm::Event::put(), EcalMGPASample::raw(), HI_PhotonSkim_cff::rechits, EcalDataFrame::sample(), reco::SuperCluster::seed(), EcalDataFrame::setSample(), funct::sin(), EcalDataFrame::size(), and ecaldqm::topology().

 {
   
   //Get BarrelSuperClusters to start.
   edm::Handle<reco::SuperClusterCollection> pBarrelSuperClusters;
   
   evt.getByToken(barrelSuperClusterProducer_, pBarrelSuperClusters);
 
   const reco::SuperClusterCollection & BarrelSuperClusters = *pBarrelSuperClusters;
   //Got BarrelSuperClusters
 
   //Get BarrelSuperClusters to start.
   edm::Handle<reco::SuperClusterCollection> pEndcapSuperClusters;
   
   evt.getByToken(endcapSuperClusterProducer_, pEndcapSuperClusters);
 
   const reco::SuperClusterCollection & EndcapSuperClusters = *pEndcapSuperClusters;
   //Got EndcapSuperClusters
 
   reco::SuperClusterCollection saveBarrelSuperClusters;
   reco::SuperClusterCollection saveEndcapSuperClusters;
   bool meet_single_thresh = false;
   //Loop over barrel superclusters, and apply threshold
   for (int loop=0;loop<int(BarrelSuperClusters.size());loop++){
     SuperCluster clus1 = BarrelSuperClusters[loop];
     float eta1 = clus1.eta();
     float energy1 = clus1.energy();
     float theta1 = 2*atan(exp(-1.*eta1));
     float cluspt1 = energy1 * sin(theta1);
     if (cluspt1 > cluster_pt_thresh_){
       saveBarrelSuperClusters.push_back(clus1);
       if (cluspt1 > single_cluster_thresh_)
     meet_single_thresh = true;
       
     }
   }
 
   //Loop over endcap superclusters, and apply threshold
   for (int loop=0;loop<int(EndcapSuperClusters.size());loop++){
     SuperCluster clus1 = EndcapSuperClusters[loop];
     float eta1 = clus1.eta();
     float energy1 = clus1.energy();
     float theta1 = 2*atan(exp(-1.*eta1));
     float cluspt1 = energy1 * sin(theta1);
     if (cluspt1 > cluster_pt_thresh_){
       saveEndcapSuperClusters.push_back(clus1);
       if (cluspt1 > single_cluster_thresh_)
     meet_single_thresh = true;
     }
   }
   
   std::auto_ptr<EBDigiCollection> SEBDigiCol(new EBDigiCollection);
   std::auto_ptr<EEDigiCollection> SEEDigiCol(new EEDigiCollection);
   int TotClus = saveBarrelSuperClusters.size() + saveEndcapSuperClusters.size();
 
   if (TotClus >= nclus_sel_ || meet_single_thresh){
     
     if (saveBarrelSuperClusters.size() > 0){
       
         
       edm::ESHandle<CaloTopology> pTopology;
       es.get<CaloTopologyRecord>().get(pTopology);
       const CaloTopology *topology = pTopology.product();   
 
       //get barrel digi collection
       edm::Handle<EBDigiCollection> pdigis;
       const EBDigiCollection* digis=0;
       evt.getByToken(EcalEBDigiToken_,pdigis);
       digis = pdigis.product(); // get a ptr to the product
 
       edm::Handle<EcalRecHitCollection> prechits;
       const EcalRecHitCollection* rechits=0;
       evt.getByToken(EcalEBRecHitToken_,prechits);
       rechits = prechits.product(); // get a ptr to the product
 
 
       if ( digis ) {
          std::vector<DetId> saveTheseDetIds;
          //pick out the detids for the 3x3 in each of the selected superclusters
          for (int loop = 0;loop < int(saveBarrelSuperClusters.size());loop++){
            SuperCluster clus1 = saveBarrelSuperClusters[loop];
            CaloClusterPtr bcref = clus1.seed();
            const BasicCluster *bc = bcref.get();
            //Get the maximum detid
            std::pair<DetId, float> EDetty = 
                EcalClusterTools::getMaximum(*bc,rechits);
            //get the 3x3 array centered on maximum detid.
            std::vector<DetId> detvec = 
                EcalClusterTools::matrixDetId(topology,EDetty.first, -1, 1, -1, 1);
            //Loop over the 3x3
            for (int ik = 0;ik<int(detvec.size());++ik)
              saveTheseDetIds.push_back(detvec[ik]);
            
          }
          for (int detloop=0; detloop < int(saveTheseDetIds.size());++detloop){
             EBDetId detL = EBDetId(saveTheseDetIds[detloop]);
            
            for (EBDigiCollection::const_iterator blah = digis->begin();
                 blah!=digis->end();blah++){
              
              if (detL == blah->id()){
                EBDataFrame myDigi = (*blah);
                SEBDigiCol->push_back(detL);
            
                EBDataFrame df( SEBDigiCol->back());
                for (int iq =0;iq<myDigi.size();++iq){
                  df.setSample(iq, myDigi.sample(iq).raw());
                }
                //ebcounter++;
              } 
            }
            //if (ebcounter >= int(saveTheseDetIds.size())) break;
          }//loop over dets
     
       }
 
     }//If barrel superclusters need saving.
     
     
     if (saveEndcapSuperClusters.size() > 0){
  
       edm::ESHandle<CaloTopology> pTopology;
       es.get<CaloTopologyRecord>().get(pTopology);
       const CaloTopology *topology = pTopology.product();   
 
       //Get endcap rec hit collection
       //get endcap digi collection
       edm::Handle<EEDigiCollection> pdigis;
       const EEDigiCollection* digis=0;
       evt.getByToken(EcalEEDigiToken_,pdigis);
       digis = pdigis.product(); // get a ptr to the product
   
       edm::Handle<EcalRecHitCollection> prechits;
       const EcalRecHitCollection* rechits=0;
       evt.getByToken(EcalEERecHitToken_,prechits);
       rechits = prechits.product(); // get a ptr to the product
 
       if ( digis ) {
          //std::vector<DetId> saveTheseDetIds;
          std::set<DetId> saveTheseDetIds;
          //pick out the digis for the 3x3 in each of the selected superclusters
          for (int loop = 0;loop < int(saveEndcapSuperClusters.size());loop++){
            SuperCluster clus1 = saveEndcapSuperClusters[loop];
            CaloClusterPtr bcref = clus1.seed();
            const BasicCluster *bc = bcref.get();
            //Get the maximum detid
            std::pair<DetId, float> EDetty = EcalClusterTools::getMaximum(*bc,rechits);
            //get the 3x3 array centered on maximum detid.
            std::vector<DetId> detvec = 
                EcalClusterTools::matrixDetId(topology,EDetty.first, -1, 1, -1, 1);
            //Loop over the 3x3
            for (int ik = 0;ik<int(detvec.size());++ik)
              saveTheseDetIds.insert(detvec[ik]);
          }
          int eecounter=0;
          for (EEDigiCollection::const_iterator blah = digis->begin();
               blah!=digis->end();blah++){
            std::set<DetId>::const_iterator finder = saveTheseDetIds.find(blah->id());
            if (finder!=saveTheseDetIds.end()){
              EEDetId detL = EEDetId(*finder);
 
              if (detL == blah->id()){
                EEDataFrame myDigi = (*blah);
                SEEDigiCol->push_back(detL);
                EEDataFrame df( SEEDigiCol->back());
                for (int iq =0;iq<myDigi.size();++iq){
                  df.setSample(iq, myDigi.sample(iq).raw());       
                }
                eecounter++;
                
              }
            }
            if (eecounter >= int(saveTheseDetIds.size())) break;
          }//loop over digis
       }
     }//If endcap superclusters need saving.
     
   }//If we're actually saving stuff
   
   //Okay, either my collections have been filled with the requisite Digis, or they haven't.
   
 
   //Empty collection, or full, still put in event.
   SEBDigiCol->sort();
   SEEDigiCol->sort();
   evt.put(SEBDigiCol, selectedEcalEBDigiCollection_);
   evt.put(SEEDigiCol, selectedEcalEEDigiCollection_);
   
 }

Member Data Documentation

edm::EDGetTokenT<reco::SuperClusterCollection> EcalDigiSelector::barrelSuperClusterProducer_

private

Definition at line 35 of file EcalDigiSelector.h.

double EcalDigiSelector::cluster_pt_thresh_

private

Definition at line 47 of file EcalDigiSelector.h.

edm::EDGetTokenT<EBDigiCollection> EcalDigiSelector::EcalEBDigiToken_

private

Definition at line 41 of file EcalDigiSelector.h.

edm::EDGetTokenT<EcalRecHitCollection> EcalDigiSelector::EcalEBRecHitToken_

private

Definition at line 39 of file EcalDigiSelector.h.

edm::EDGetTokenT<EEDigiCollection> EcalDigiSelector::EcalEEDigiToken_

private

Definition at line 42 of file EcalDigiSelector.h.

edm::EDGetTokenT<EcalRecHitCollection> EcalDigiSelector::EcalEERecHitToken_

private

Definition at line 40 of file EcalDigiSelector.h.

edm::EDGetTokenT<reco::SuperClusterCollection> EcalDigiSelector::endcapSuperClusterProducer_

private

Definition at line 36 of file EcalDigiSelector.h.

int EcalDigiSelector::nclus_sel_

private

Definition at line 49 of file EcalDigiSelector.h.

std::string EcalDigiSelector::selectedEcalEBDigiCollection_

private

Definition at line 32 of file EcalDigiSelector.h.

std::string EcalDigiSelector::selectedEcalEEDigiCollection_

private

Definition at line 33 of file EcalDigiSelector.h.

double EcalDigiSelector::single_cluster_thresh_

private

Definition at line 48 of file EcalDigiSelector.h.

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation