ESListOfFEDSProducer Class Reference

#include <ESListOfFEDSProducer.h>

Inheritance diagram for ESListOfFEDSProducer:

Public Member Functions
void	beginJob (void)
void	Egamma (edm::Event &e, const edm::EventSetup &es, std::vector< int > &done, std::vector< int > &FEDs)
void	endJob (void)
	ESListOfFEDSProducer (const edm::ParameterSet &pset)
void	Jets (edm::Event &e, const edm::EventSetup &es, std::vector< int > &done, std::vector< int > &FEDs)
void	Muon (edm::Event &e, const edm::EventSetup &es, std::vector< int > &done, std::vector< int > &FEDs)
void	produce (edm::Event &e, const edm::EventSetup &c)
virtual	~ESListOfFEDSProducer ()
Public Member Functions inherited from edm::EDProducer
	EDProducer ()
ModuleDescription const &	moduleDescription () const
virtual	~EDProducer ()
Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const
	ProducerBase ()
void	registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...
void	resolvePutIndicies (BranchType iBranchType, std::unordered_multimap< std::string, edm::ProductResolverIndex > const &iIndicies, std::string const &moduleLabel)
virtual	~ProducerBase () noexcept(false)
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...
	EDConsumerBase ()
	EDConsumerBase (EDConsumerBase const &)=delete
	EDConsumerBase (EDConsumerBase &&)=default
ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
EDConsumerBase &	operator= (EDConsumerBase &&)=default
bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
virtual	~EDConsumerBase () noexcept(false)

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)
Static Public Member Functions inherited from edm::EDProducer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)

Private Member Functions
std::vector< int >	ListOfFEDS (double etaLow, double etaHigh, double phiLow, double phiHigh, double etamargin, double phimargin)

Private Attributes
edm::EDGetTokenT< l1extra::L1JetParticleCollection >	CentralSource_
bool	debug_
bool	EGamma_
bool	EMdoIsolated_
bool	EMdoNonIsolated_
edm::EDGetTokenT< l1extra::L1EmParticleCollection >	EMl1TagIsolated_
edm::EDGetTokenT< l1extra::L1EmParticleCollection >	EMl1TagNonIsolated_
double	EMregionEtaMargin_
double	EMregionPhiMargin_
bool	first_
edm::EDGetTokenT< l1extra::L1JetParticleCollection >	ForwardSource_
bool	Jets_
bool	JETSdoCentral_
bool	JETSdoForward_
bool	JETSdoTau_
double	JETSregionEtaMargin_
double	JETSregionPhiMargin_
bool	Muon_
edm::EDGetTokenT< l1extra::L1MuonParticleCollection >	MuonSource_
double	MUregionEtaMargin_
double	MUregionPhiMargin_
std::string	OutputLabel_
double	Ptmin_iso_
double	Ptmin_jets_
double	Ptmin_muon_
double	Ptmin_noniso_
edm::EDGetTokenT< l1extra::L1JetParticleCollection >	TauSource_
ESElectronicsMapper *	TheESMapping
	ESElectronicsMapping* TheESMapping;. More...
EcalElectronicsMapping *	TheMapping

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType
Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
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 31 of file ESListOfFEDSProducer.h.

Constructor & Destructor Documentation

ESListOfFEDSProducer::ESListOfFEDSProducer

(

const edm::ParameterSet &

pset

)

Definition at line 31 of file ESListOfFEDSProducer.cc.

References Exception, edm::ParameterSet::getUntrackedParameter(), and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                      {
  
  debug_ = pset.getUntrackedParameter<bool>("debug");
  
  // Pi0ListToIgnore_ = pset.getParameter<edm::InputTag>("Pi0ListToIgnore");
  EGamma_ = pset.getUntrackedParameter<bool>("EGamma",false);
  Muon_ = pset.getUntrackedParameter<bool>("Muon",false);
  Jets_ = pset.getUntrackedParameter<bool>("Jets",false);
  
  if (EGamma_ && Muon_) {
    throw cms::Exception("ESListOfFEDSProducer") << 
      " Wrong configuration : EGamma and Muon should not be true at the same time." ;
  }
    
  if (EGamma_) {
    EMl1TagIsolated_    = consumes<l1extra::L1EmParticleCollection>(pset.getUntrackedParameter<edm::InputTag>("EM_l1TagIsolated"));
    EMl1TagNonIsolated_ = consumes<l1extra::L1EmParticleCollection>(pset.getUntrackedParameter<edm::InputTag>("EM_l1TagNonIsolated"));
    EMdoIsolated_       = pset.getUntrackedParameter<bool>("EM_doIsolated",true);
    EMdoNonIsolated_    = pset.getUntrackedParameter<bool>("EM_doNonIsolated",true);
    EMregionEtaMargin_  = pset.getUntrackedParameter<double>("EM_regionEtaMargin",0.25);
    EMregionPhiMargin_  = pset.getUntrackedParameter<double>("EM_regionPhiMargin",0.40);
    Ptmin_iso_          = pset.getUntrackedParameter<double>("Ptmin_iso",0.);
    Ptmin_noniso_       = pset.getUntrackedParameter<double>("Ptmin_noniso",0.);
  }

  if (Muon_) {
    MUregionEtaMargin_ = pset.getUntrackedParameter<double>("MU_regionEtaMargin",1.0);
    MUregionPhiMargin_ = pset.getUntrackedParameter<double>("MU_regionPhiMargin",1.0);
    Ptmin_muon_        = pset.getUntrackedParameter<double>("Ptmin_muon",0.);
    MuonSource_        = consumes<l1extra::L1MuonParticleCollection>(pset.getUntrackedParameter<edm::InputTag>("MuonSource"));
  }
  
  if (Jets_) {
    JETSregionEtaMargin_ = pset.getUntrackedParameter<double>("JETS_regionEtaMargin",1.0);
    JETSregionPhiMargin_ = pset.getUntrackedParameter<double>("JETS_regionPhiMargin",1.0);
    Ptmin_jets_          = pset.getUntrackedParameter<double>("Ptmin_jets",0.);
    CentralSource_       = consumes<l1extra::L1JetParticleCollection>(pset.getUntrackedParameter<edm::InputTag>("CentralSource"));
    ForwardSource_       = consumes<l1extra::L1JetParticleCollection>(pset.getUntrackedParameter<edm::InputTag>("ForwardSource"));
    TauSource_           = consumes<l1extra::L1JetParticleCollection>(pset.getUntrackedParameter<edm::InputTag>("TauSource"));
    JETSdoCentral_       = pset.getUntrackedParameter<bool>("JETS_doCentral",true);
    JETSdoForward_       = pset.getUntrackedParameter<bool>("JETS_doForward",true);
    JETSdoTau_           = pset.getUntrackedParameter<bool>("JETS_doTau",true);
  }

  OutputLabel_ = pset.getUntrackedParameter<std::string>("OutputLabel");
  
  TheMapping = new EcalElectronicsMapping();
  first_ = true;
  
  TheESMapping = new ESElectronicsMapper(pset);
  
  consumesMany<ESListOfFEDS>();
  produces<ESListOfFEDS>(OutputLabel_);
}

ESListOfFEDSProducer::~ESListOfFEDSProducer ( )

virtual

Definition at line 86 of file ESListOfFEDSProducer.cc.

                                            {
 delete TheMapping;
}

Member Function Documentation

void ESListOfFEDSProducer::beginJob ( void  )

virtual

Reimplemented from edm::EDProducer.

Definition at line 123 of file ESListOfFEDSProducer.cc.

{}

void ESListOfFEDSProducer::Egamma	(	edm::Event &	e,
		const edm::EventSetup &	es,
		std::vector< int > &	done,
		std::vector< int > &	FEDs
	)

Definition at line 217 of file ESListOfFEDSProducer.cc.

References gather_cfg::cout, L1CaloGeometry::emJetPhiBinHighEdge(), L1CaloGeometry::emJetPhiBinLowEdge(), L1CaloGeometry::etaBinHighEdge(), L1CaloGeometry::etaBinLowEdge(), spr::find(), edm::EventSetup::get(), edm::Event::getByToken(), mps_fire::i, createfilelist::int, and EnergyCorrector::pt.

                                                                                                                   {
  
  // std::vector<int> FEDs;
  
  if (debug_) std::cout << std::endl << std::endl << " enter in EcalListOfFEDSProducer::Egamma" << std::endl;
  
  //Get the L1 EM Particle Collection
  edm::Handle<l1extra::L1EmParticleCollection > emIsolColl ;
  if(EMdoIsolated_)
    e.getByToken(EMl1TagIsolated_, emIsolColl);
  //Get the L1 EM Particle Collection
  edm::Handle<l1extra::L1EmParticleCollection > emNonIsolColl ;
  if (EMdoNonIsolated_)
    e.getByToken(EMl1TagNonIsolated_, emNonIsolColl);
  
  // Get the CaloGeometry
  edm::ESHandle<L1CaloGeometry> l1CaloGeom ;
  es.get<L1CaloGeometryRecord>().get(l1CaloGeom) ;
  
  if(EMdoIsolated_) {
    
    for( l1extra::L1EmParticleCollection::const_iterator emItr = emIsolColl->begin();
     emItr != emIsolColl->end() ;++emItr ){
      
      float pt = emItr -> pt();
      if (pt < Ptmin_iso_ ) continue;
      if (debug_) std::cout << " Here is an L1 isoEM candidate of pt " << pt << std::endl;
      // Access the GCT hardware object corresponding to the L1Extra EM object.
      int etaIndex = emItr->gctEmCand()->etaIndex() ;
      int phiIndex = emItr->gctEmCand()->phiIndex() ;
      // Use the L1CaloGeometry to find the eta, phi bin boundaries.
      double etaLow  = l1CaloGeom->etaBinLowEdge( etaIndex ) ;
      double etaHigh = l1CaloGeom->etaBinHighEdge( etaIndex ) ;
      double phiLow  = l1CaloGeom->emJetPhiBinLowEdge( phiIndex ) ;
      double phiHigh = l1CaloGeom->emJetPhiBinHighEdge( phiIndex ) ;
      
      std::vector<int> feds = ListOfFEDS(etaLow, etaHigh, phiLow, phiHigh, EMregionEtaMargin_, EMregionPhiMargin_);
      for (int i=0; i < (int)feds.size(); i++) {
    if ( std::find(FEDs.begin(), FEDs.end(), feds[i]) == FEDs.end() &&
         std::find(done.begin(), done.end(), feds[i]) == done.end() ) FEDs.push_back(feds[i]);
      }
    } // end loop on L1EmParticleCollection
  }  // endif doIsolated_

  if (EMdoNonIsolated_) {
    
    for( l1extra::L1EmParticleCollection::const_iterator emItr = emNonIsolColl->begin();
     emItr != emNonIsolColl->end() ;++emItr ){
      
      float pt = emItr -> pt();
      if (debug_) std::cout << " Here is an L1 nonisoEM candidate of pt " << pt << std::endl;
      if (pt < Ptmin_noniso_ ) continue;
      // Access the GCT hardware object corresponding to the L1Extra EM object.
      int etaIndex = emItr->gctEmCand()->etaIndex() ;
      int phiIndex = emItr->gctEmCand()->phiIndex() ;
      // std::cout << " etaIndex phiIndex " << etaIndex << " " << phiIndex << std::endl;
      // Use the L1CaloGeometry to find the eta, phi bin boundaries.
      double etaLow  = l1CaloGeom->etaBinLowEdge( etaIndex ) ;
      double etaHigh = l1CaloGeom->etaBinHighEdge( etaIndex ) ;
      double phiLow  = l1CaloGeom->emJetPhiBinLowEdge( phiIndex ) ;
      double phiHigh = l1CaloGeom->emJetPhiBinHighEdge( phiIndex ) ;
      
      std::vector<int> feds = ListOfFEDS(etaLow, etaHigh, phiLow, phiHigh, EMregionEtaMargin_, EMregionPhiMargin_);
      for (int i=0; i < (int)feds.size(); i++) {
    if ( std::find(FEDs.begin(), FEDs.end(), feds[i]) == FEDs.end() &&
         std::find(done.begin(), done.end(), feds[i]) == done.end() ) FEDs.push_back(feds[i]);
    
      }
    } // end loop on L1EmParticleCollection
  }
  
  if (debug_) {
    std::cout << std::endl;
    for (int i=0; i < (int)FEDs.size(); i++) {
      std::cout << "Egamma: unpack FED " << FEDs[i] << std::endl;
    }
    std::cout << "Number of FEDS is " << FEDs.size() << std::endl;
  }
  
  // return FEDs;
  
}

void ESListOfFEDSProducer::endJob ( void  )

virtual

Reimplemented from edm::EDProducer.

Definition at line 126 of file ESListOfFEDSProducer.cc.

{}

void ESListOfFEDSProducer::fillDescriptions ( edm::ConfigurationDescriptions &  descriptions )

static

Definition at line 90 of file ESListOfFEDSProducer.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                      {

  edm::ParameterSetDescription desc;  
  desc.add<bool>("debug", false);
  desc.add<bool>("EGamma", true);
  desc.add<bool>("Muon", false);
  desc.add<bool>("Jets", false);
  desc.add<edm::InputTag>("EM_l1TagIsolated", edm::InputTag("l1extraParticles","Isolated"));
  desc.add<edm::InputTag>("EM_l1TagNonIsolated", edm::InputTag("l1extraParticles","NonIsolated"));
  desc.add<bool>("EM_doIsolated", true);
  desc.add<bool>("EM_doNonIsolated", true);
  desc.add<double>("EM_regionEtaMargin", 0.25);
  desc.add<double>("EM_regionPhiMargin", 0.40);
  desc.add<double>("Ptmin_iso", 5.0);
  desc.add<double>("Ptmin_noniso", 5.0);
  desc.add<double>("MU_regionEtaMargin", 1.0);
  desc.add<double>("MU_regionPhiMargin", 1.0);
  desc.add<double>("Ptmin_muon", 0.);
  desc.add<edm::InputTag>("MuonSource", edm::InputTag("l1extraParticles"));
  desc.add<double>("JETS_regionEtaMargin", 1.0);
  desc.add<double>("JETS_regionPhiMargin", 1.0);
  desc.add<double>("Ptmin_jets", 50.);
  desc.add<edm::InputTag>("CentralSource", edm::InputTag("l1extraParticles","Central"));
  desc.add<edm::InputTag>("ForwardSource", edm::InputTag("l1extraParticles","Forward"));
  desc.add<edm::InputTag>("TauSource", edm::InputTag("l1extraParticles","Tau"));
  desc.add<bool>("JETS_doCentral", true);
  desc.add<bool>("JETS_doForward", true);
  desc.add<bool>("JETS_doTau", true);
  desc.add<std::string>("OutputLabel", "");
  descriptions.add(("hltESListOfFEDSProducer"), desc);  
}

void ESListOfFEDSProducer::Jets	(	edm::Event &	e,
		const edm::EventSetup &	es,
		std::vector< int > &	done,
		std::vector< int > &	FEDs
	)

Definition at line 339 of file ESListOfFEDSProducer.cc.

References gather_cfg::cout, geometryDiff::epsilon, PVValHelper::eta, spr::find(), edm::Event::getByToken(), mps_fire::i, createfilelist::int, and EnergyCorrector::pt.

                                                                                                                {
  
  // std::vector<int> FEDs;
  
  if (debug_) std::cout << std::endl << std::endl << " enter in EcalListOfFEDSProducer::Jets" << std::endl;
  double epsilon = 0.01;

  if (JETSdoCentral_) {
    
    edm::Handle<L1JetParticleCollection> jetColl;
    e.getByToken(CentralSource_,jetColl);
    
    for (L1JetParticleCollection::const_iterator it=jetColl->begin(); it != jetColl->end(); it++) {
      
      double pt    =  it -> pt();
      double eta   =  it -> eta();
      double phi   =  it -> phi();
      
      if (debug_) std::cout << " here is a L1 CentralJet Seed  with (eta,phi) = " <<
            eta << " " << phi << " and pt " << pt << std::endl;
      if (pt < Ptmin_jets_ ) continue;
      
      std::vector<int> feds = ListOfFEDS(eta, eta, phi-epsilon, phi+epsilon, JETSregionEtaMargin_, JETSregionPhiMargin_);
      
      for (int i=0; i < (int)feds.size(); i++) {
    if ( std::find(FEDs.begin(), FEDs.end(), feds[i]) == FEDs.end() &&
         std::find(done.begin(), done.end(), feds[i]) == done.end() ) FEDs.push_back(feds[i]);
      }
    }
  }
  
  if (JETSdoForward_) {
    
    edm::Handle<L1JetParticleCollection> jetColl;
    e.getByToken(ForwardSource_,jetColl);

    for (L1JetParticleCollection::const_iterator it=jetColl->begin(); it != jetColl->end(); it++) {
      
      double pt    =  it -> pt();
      double eta   =  it -> eta();
      double phi   =  it -> phi();
      
      if (debug_) std::cout << " here is a L1 ForwardJet Seed  with (eta,phi) = " <<
            eta << " " << phi << " and pt " << pt << std::endl;
      if (pt < Ptmin_jets_ ) continue;
      
      std::vector<int> feds = ListOfFEDS(eta, eta, phi-epsilon, phi+epsilon, JETSregionEtaMargin_, JETSregionPhiMargin_);
      
      for (int i=0; i < (int)feds.size(); i++) {
    if ( std::find(FEDs.begin(), FEDs.end(), feds[i]) == FEDs.end() &&
         std::find(done.begin(), done.end(), feds[i]) == done.end() ) FEDs.push_back(feds[i]);
      }
    }
  }
  
  if (JETSdoTau_) {
    
    edm::Handle<L1JetParticleCollection> jetColl;
    e.getByToken(TauSource_,jetColl);
    
    for (L1JetParticleCollection::const_iterator it=jetColl->begin(); it != jetColl->end(); it++) {
      
      double pt    =  it -> pt();
      double eta   =  it -> eta();
      double phi   =  it -> phi();
      
      if (debug_) std::cout << " here is a L1 TauJet Seed  with (eta,phi) = " <<
            eta << " " << phi << " and pt " << pt << std::endl;
      if (pt < Ptmin_jets_ ) continue;
      
      std::vector<int> feds = ListOfFEDS(eta, eta, phi-epsilon, phi+epsilon, JETSregionEtaMargin_, JETSregionPhiMargin_);
      
      for (int i=0; i < (int)feds.size(); i++) {
    if ( std::find(FEDs.begin(), FEDs.end(), feds[i]) == FEDs.end() &&
         std::find(done.begin(), done.end(), feds[i]) == done.end() ) FEDs.push_back(feds[i]);
      }
    }
  }
    
  if (debug_) {
    std::cout << std::endl;
    for (int i=0; i < (int)FEDs.size(); i++) {
      std::cout << "Jets: unpack FED " << FEDs[i] << std::endl;
    }
    std::cout << "Number of FEDS is " << FEDs.size() << std::endl;
  }
  
  // return FEDs;
}

std::vector< int > ESListOfFEDSProducer::ListOfFEDS ( double  etaLow, double  etaHigh, double  phiLow, double  phiHigh, double  etamargin, double  phimargin  )

private

Definition at line 430 of file ESListOfFEDSProducer.cc.

References Vispa.Plugins.EdmBrowser.EdmDataAccessor::all(), gather_cfg::cout, createTree::dd, and Geom::pi().

{
  
  std::vector<int> FEDs;
  
  if (phimargin > Geom::pi()) phimargin =  Geom::pi() ;
  
  
  if (debug_) std::cout << " etaLow etaHigh phiLow phiHigh " << etaLow << " " << 
        etaHigh << " " << phiLow << " " << phiHigh << std::endl;
  
  etaLow -= etamargin;
  etaHigh += etamargin;
  double phiMinus = phiLow - phimargin;
  double phiPlus = phiHigh + phimargin;
  
  bool all = false;
  double dd = fabs(phiPlus-phiMinus);
  if (debug_) std::cout << " dd = " << dd << std::endl;
  if (dd > 2.*Geom::pi() ) all = true;
  
  while (phiPlus > Geom::pi()) { phiPlus -= 2.*Geom::pi() ; }
  while (phiMinus < 0) { phiMinus += 2.*Geom::pi() ; }
  if ( phiMinus > Geom::pi()) phiMinus -= 2.*Geom::pi() ;
  
  double dphi = phiPlus - phiMinus;
  if (dphi < 0) dphi += 2.*Geom::pi() ;
  if (debug_) std::cout << "dphi = " << dphi << std::endl;
  if (dphi > Geom::pi()) {
    int fed_low1 = TheMapping -> GetFED(etaLow,phiMinus*180./Geom::pi());
    int fed_low2 = TheMapping -> GetFED(etaLow,phiPlus*180./Geom::pi());
    if (debug_) std::cout << "fed_low1 fed_low2 " << fed_low1 << " " << fed_low2 << std::endl;
    if (fed_low1 == fed_low2) all = true;
    int fed_hi1 = TheMapping -> GetFED(etaHigh,phiMinus*180./Geom::pi());
    int fed_hi2 = TheMapping -> GetFED(etaHigh,phiPlus*180./Geom::pi());
    if (debug_) std::cout << "fed_hi1 fed_hi2 " << fed_hi1 << " " << fed_hi2 << std::endl;
    if (fed_hi1 == fed_hi2) all = true;
  }
  
  
  if (all) {
    if (debug_) std::cout << " unpack everything in phi ! " << std::endl;
    phiMinus = -20 * Geom::pi() / 180.;  // -20 deg
    phiPlus = -40 * Geom::pi() / 180.;  // -20 deg
  }
  
  if (debug_) std::cout << " with margins : " << etaLow << " " << etaHigh << " " << 
        phiMinus << " " << phiPlus << std::endl;
  
  
  const EcalEtaPhiRegion ecalregion(etaLow,etaHigh,phiMinus,phiPlus);
  
  FEDs = TheMapping -> GetListofFEDs(ecalregion);
  
  return FEDs;
}

void ESListOfFEDSProducer::Muon	(	edm::Event &	e,
		const edm::EventSetup &	es,
		std::vector< int > &	done,
		std::vector< int > &	FEDs
	)

Definition at line 300 of file ESListOfFEDSProducer.cc.

References gather_cfg::cout, geometryDiff::epsilon, PVValHelper::eta, spr::find(), edm::Event::getByToken(), mps_fire::i, createfilelist::int, and EnergyCorrector::pt.

                                                                                                                {
  
  if (debug_) std::cout << std::endl << std::endl << " enter in EcalListOfFEDSProducer::Muon" << std::endl;
  edm::Handle<L1MuonParticleCollection> muColl;
  e.getByToken(MuonSource_, muColl);

  double epsilon = 0.01;
  
  for (L1MuonParticleCollection::const_iterator it=muColl->begin(); it != muColl->end(); it++) {
    
    const L1MuGMTExtendedCand muonCand = (*it).gmtMuonCand();
    double pt    =  (*it).pt();
    double eta   =  (*it).eta();
    double phi   =  (*it).phi();
    
    if (debug_) std::cout << " here is a L1 muon Seed  with (eta,phi) = " << 
          eta << " " << phi << " and pt " << pt << std::endl;
    if (pt < Ptmin_muon_ ) continue;
    
    std::vector<int> feds = ListOfFEDS(eta, eta, phi-epsilon, phi+epsilon, MUregionEtaMargin_, MUregionPhiMargin_);
    
    for (int i=0; i < (int)feds.size(); i++) {
      if ( std::find(FEDs.begin(), FEDs.end(), feds[i]) == FEDs.end() &&
       std::find(done.begin(), done.end(), feds[i]) == done.end() ) FEDs.push_back(feds[i]);
    }
  }
  
  if (debug_) {
    std::cout << std::endl;
    for (int i=0; i < (int)FEDs.size(); i++) {
      std::cout << "Muons: unpack FED " << FEDs[i] << std::endl;
    }
    std::cout << "Number of FEDS is " << FEDs.size() << std::endl;
  }
  
  // return FEDs;
  
}

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

commented out at this moment, for muon and jets.

get ES List of FEDs ,given the ecal FED list

now push list of ES FEDs into event.

Definition at line 129 of file ESListOfFEDSProducer.cc.

References gather_cfg::cout, edm::EventID::event(), spr::find(), edm::EventSetup::get(), edm::Event::getManyByType(), mps_fire::i, edm::EventBase::id(), createfilelist::int, FEDNumbering::MINECALFEDID, eostools::move(), gen::n, edm::ESHandle< T >::product(), edm::Event::put(), and edm::EventID::run().

                                                                           {

  if (first_) {
    edm::ESHandle< EcalElectronicsMapping > ecalmapping;
    iSetup.get< EcalMappingRcd >().get(ecalmapping);
    const EcalElectronicsMapping* TheMapping_ = ecalmapping.product();
    *TheMapping = *TheMapping_;
    first_ = false;
  }                                                                                              
  
  auto productAddress = std::make_unique<ESListOfFEDS>();
  std::vector<int> feds;        // the list of Ecal FEDS produced 
  
  std::vector<int> Done;
  if (EGamma_) {
    // feds = Egamma(e, iSetup, Done);
    Egamma(e, iSetup, Done, feds);
  }
  
  if (Muon_) {
    // feds = Muon(e, iSetup, Done);
    // Muon(e, iSetup, Done, feds);
  }
  
  if (Jets_) {
    // feds = Jets(e, iSetup, Done);
    //Jets(e, iSetup, Done, feds);
  }
  
  if ( !EGamma_ && !Muon_ && ! Jets_)  {
    for (int i=1; i <= 54; i++) {
      if ( std::find(Done.begin(), Done.end(), i) == Done.end())
    feds.push_back(i);
    }
  }
  
  
  int nf = (int)feds.size();
  for (int i=0; i <nf; i++) {
    feds[i] += FEDNumbering::MINECALFEDID;
    if (debug_) std::cout << "Will unpack FED " << feds[i] << std::endl;
  }
  
  if (debug_ && nf < 1 ) 
    std::cout << " Warning : no ECAL FED to unpack for Run " << e.id().run() << "  Event " << e.id().event() << std::endl;
  
  std::vector<int> es_feds =  TheESMapping->GetListofFEDs(feds);
  
  //check ES List if some of them already done
  std::vector< edm::Handle<ESListOfFEDS> > esFEDs_Done;
  e.getManyByType(esFEDs_Done);
  unsigned int esnDone = esFEDs_Done.size();
  
  if(debug_){
    std::cout<<int(es_feds.size())<<" feds requested, and ESunpacking have already been done "<<esnDone<<" times."<<std::endl;
  }
  
  for(int i=0; i<int(esFEDs_Done.size());i++){
    std::vector<int> done = esFEDs_Done[i]-> GetList();
    for(int n=0; n<int(done.size());n++){
      for(int j=0; j<int(es_feds.size());j++){
    std::vector<int>::iterator it = find(es_feds.begin(),es_feds.end(),done[n]);
    if( it != es_feds.end()){
      es_feds.erase(it);
      j--; 
    }
      }
    }
  }
  
  if (debug_ ){
    std::cout<<"ESListofFEDSProducer: new ES_list will be put into event: total "<<int(es_feds.size())<<std::endl;
    for(int j=0; j<int(es_feds.size());j++){
      std::cout<<es_feds[j]<<" ";
    }
    std::cout<<std::endl;
  }
  
  
  productAddress.get() -> SetList(es_feds);
  e.put(std::move(productAddress),OutputLabel_);
  
}

Member Data Documentation

edm::EDGetTokenT<l1extra::L1JetParticleCollection> ESListOfFEDSProducer::CentralSource_

private

Definition at line 69 of file ESListOfFEDSProducer.h.

bool ESListOfFEDSProducer::debug_

private

Definition at line 80 of file ESListOfFEDSProducer.h.

bool ESListOfFEDSProducer::EGamma_

private

Definition at line 46 of file ESListOfFEDSProducer.h.

bool ESListOfFEDSProducer::EMdoIsolated_

private

Definition at line 49 of file ESListOfFEDSProducer.h.

bool ESListOfFEDSProducer::EMdoNonIsolated_

private

Definition at line 50 of file ESListOfFEDSProducer.h.

edm::EDGetTokenT<l1extra::L1EmParticleCollection> ESListOfFEDSProducer::EMl1TagIsolated_

private

Definition at line 47 of file ESListOfFEDSProducer.h.

edm::EDGetTokenT<l1extra::L1EmParticleCollection> ESListOfFEDSProducer::EMl1TagNonIsolated_

private

Definition at line 48 of file ESListOfFEDSProducer.h.

double ESListOfFEDSProducer::EMregionEtaMargin_

private

Definition at line 51 of file ESListOfFEDSProducer.h.

double ESListOfFEDSProducer::EMregionPhiMargin_

private

Definition at line 52 of file ESListOfFEDSProducer.h.

bool ESListOfFEDSProducer::first_

private

Definition at line 79 of file ESListOfFEDSProducer.h.

edm::EDGetTokenT<l1extra::L1JetParticleCollection> ESListOfFEDSProducer::ForwardSource_

private

Definition at line 70 of file ESListOfFEDSProducer.h.

bool ESListOfFEDSProducer::Jets_

private

Definition at line 62 of file ESListOfFEDSProducer.h.

bool ESListOfFEDSProducer::JETSdoCentral_

private

Definition at line 63 of file ESListOfFEDSProducer.h.

bool ESListOfFEDSProducer::JETSdoForward_

private

Definition at line 64 of file ESListOfFEDSProducer.h.

bool ESListOfFEDSProducer::JETSdoTau_

private

Definition at line 65 of file ESListOfFEDSProducer.h.

double ESListOfFEDSProducer::JETSregionEtaMargin_

private

Definition at line 66 of file ESListOfFEDSProducer.h.

double ESListOfFEDSProducer::JETSregionPhiMargin_

private

Definition at line 67 of file ESListOfFEDSProducer.h.

bool ESListOfFEDSProducer::Muon_

private

Definition at line 56 of file ESListOfFEDSProducer.h.

edm::EDGetTokenT<l1extra::L1MuonParticleCollection> ESListOfFEDSProducer::MuonSource_

private

Definition at line 60 of file ESListOfFEDSProducer.h.

double ESListOfFEDSProducer::MUregionEtaMargin_

private

Definition at line 57 of file ESListOfFEDSProducer.h.

double ESListOfFEDSProducer::MUregionPhiMargin_

private

Definition at line 58 of file ESListOfFEDSProducer.h.

std::string ESListOfFEDSProducer::OutputLabel_

private

Definition at line 73 of file ESListOfFEDSProducer.h.

double ESListOfFEDSProducer::Ptmin_iso_

private

Definition at line 53 of file ESListOfFEDSProducer.h.

double ESListOfFEDSProducer::Ptmin_jets_

private

Definition at line 68 of file ESListOfFEDSProducer.h.

double ESListOfFEDSProducer::Ptmin_muon_

private

Definition at line 59 of file ESListOfFEDSProducer.h.

double ESListOfFEDSProducer::Ptmin_noniso_

private

Definition at line 54 of file ESListOfFEDSProducer.h.

edm::EDGetTokenT<l1extra::L1JetParticleCollection> ESListOfFEDSProducer::TauSource_

private

Definition at line 71 of file ESListOfFEDSProducer.h.

ESElectronicsMapper* ESListOfFEDSProducer::TheESMapping

private

ESElectronicsMapping* TheESMapping;.

Definition at line 77 of file ESListOfFEDSProducer.h.

EcalElectronicsMapping* ESListOfFEDSProducer::TheMapping

private

Definition at line 74 of file ESListOfFEDSProducer.h.