ESListOfFEDSProducer Class Reference

#include <ESListOfFEDSProducer.h>

Inheritance diagram for ESListOfFEDSProducer:

Public Member Functions
void	Egamma (edm::Event &e, const edm::EventSetup &es, std::vector< int > &done, std::vector< int > &FEDs)
	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) override
	~ESListOfFEDSProducer () override
Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default
bool	hasAbilityToProduceInLumis () const final
bool	hasAbilityToProduceInRuns () const final

Static Public Member Functions
static void	fillDescriptions (edm::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::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

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 ( )

override

Definition at line 86 of file ESListOfFEDSProducer.cc.

                                            {
 delete TheMapping;
}

Member Function Documentation

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

Definition at line 211 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::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 333 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 424 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 RectangularEtaPhiRegion 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 294 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  )

override

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 123 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 67 of file ESListOfFEDSProducer.h.

bool ESListOfFEDSProducer::debug_

private

Definition at line 78 of file ESListOfFEDSProducer.h.

bool ESListOfFEDSProducer::EGamma_

private

Definition at line 44 of file ESListOfFEDSProducer.h.

bool ESListOfFEDSProducer::EMdoIsolated_

private

Definition at line 47 of file ESListOfFEDSProducer.h.

bool ESListOfFEDSProducer::EMdoNonIsolated_

private

Definition at line 48 of file ESListOfFEDSProducer.h.

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

private

Definition at line 45 of file ESListOfFEDSProducer.h.

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

private

Definition at line 46 of file ESListOfFEDSProducer.h.

double ESListOfFEDSProducer::EMregionEtaMargin_

private

Definition at line 49 of file ESListOfFEDSProducer.h.

double ESListOfFEDSProducer::EMregionPhiMargin_

private

Definition at line 50 of file ESListOfFEDSProducer.h.

bool ESListOfFEDSProducer::first_

private

Definition at line 77 of file ESListOfFEDSProducer.h.

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

private

Definition at line 68 of file ESListOfFEDSProducer.h.

bool ESListOfFEDSProducer::Jets_

private

Definition at line 60 of file ESListOfFEDSProducer.h.

bool ESListOfFEDSProducer::JETSdoCentral_

private

Definition at line 61 of file ESListOfFEDSProducer.h.

bool ESListOfFEDSProducer::JETSdoForward_

private

Definition at line 62 of file ESListOfFEDSProducer.h.

bool ESListOfFEDSProducer::JETSdoTau_

private

Definition at line 63 of file ESListOfFEDSProducer.h.

double ESListOfFEDSProducer::JETSregionEtaMargin_

private

Definition at line 64 of file ESListOfFEDSProducer.h.

double ESListOfFEDSProducer::JETSregionPhiMargin_

private

Definition at line 65 of file ESListOfFEDSProducer.h.

bool ESListOfFEDSProducer::Muon_

private

Definition at line 54 of file ESListOfFEDSProducer.h.

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

private

Definition at line 58 of file ESListOfFEDSProducer.h.

double ESListOfFEDSProducer::MUregionEtaMargin_

private

Definition at line 55 of file ESListOfFEDSProducer.h.

double ESListOfFEDSProducer::MUregionPhiMargin_

private

Definition at line 56 of file ESListOfFEDSProducer.h.

std::string ESListOfFEDSProducer::OutputLabel_

private

Definition at line 71 of file ESListOfFEDSProducer.h.

double ESListOfFEDSProducer::Ptmin_iso_

private

Definition at line 51 of file ESListOfFEDSProducer.h.

double ESListOfFEDSProducer::Ptmin_jets_

private

Definition at line 66 of file ESListOfFEDSProducer.h.

double ESListOfFEDSProducer::Ptmin_muon_

private

Definition at line 57 of file ESListOfFEDSProducer.h.

double ESListOfFEDSProducer::Ptmin_noniso_

private

Definition at line 52 of file ESListOfFEDSProducer.h.

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

private

Definition at line 69 of file ESListOfFEDSProducer.h.

ESElectronicsMapper* ESListOfFEDSProducer::TheESMapping

private

ESElectronicsMapping* TheESMapping;.

Definition at line 75 of file ESListOfFEDSProducer.h.

EcalElectronicsMapping* ESListOfFEDSProducer::TheMapping

private

Definition at line 72 of file ESListOfFEDSProducer.h.