EcalDeadCellDeltaRFilter Class Reference

#include <EcalDeadCellDeltaRFilter.cc>

Inheritance diagram for EcalDeadCellDeltaRFilter:

Public Member Functions
	EcalDeadCellDeltaRFilter (const edm::ParameterSet &)
	~EcalDeadCellDeltaRFilter () override
Public Member Functions inherited from edm::EDFilter
	EDFilter ()
SerialTaskQueue *	globalLuminosityBlocksQueue ()
SerialTaskQueue *	globalRunsQueue ()
ModuleDescription const &	moduleDescription () const
	~EDFilter () override
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 ()
std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const
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, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)
	~ProducerBase () noexcept(false) override
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
ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const
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
ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)
virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
void	beginJob () override
void	beginLuminosityBlock (const edm::LuminosityBlock &, const edm::EventSetup &) override
void	beginRun (const edm::Run &, const edm::EventSetup &) override
int	dPhiToMETfunc (const std::vector< reco::Jet > &jetTVec, const double &dPhiCutVal, std::vector< reco::Jet > &closeToMETjetsVec)
int	dRtoMaskedChnsEvtFilterFunc (const std::vector< reco::Jet > &jetTVec, const int &chnStatus, const double &dRCutVal)
void	endJob () override
void	endRun (const edm::Run &, const edm::EventSetup &) override
virtual void	envSet (const edm::EventSetup &)
int	etaToBoundary (const std::vector< reco::Jet > &jetTVec)
bool	filter (edm::Event &, const edm::EventSetup &) override
int	getChannelStatusMaps ()
int	isCloseToBadEcalChannel (const reco::Jet &jet, const double &deltaRCut, const int &chnStatus, std::map< double, DetId > &deltaRdetIdMap)
void	loadEventInfo (const edm::Event &iEvent, const edm::EventSetup &iSetup)
void	loadJets (const edm::Event &iEvent, const edm::EventSetup &iSetup)
void	loadMET (const edm::Event &iEvent, const edm::EventSetup &iSetup)

Private Attributes
double	calomet
double	calometPhi
const int	chnStatusToBeEvaluated_
const std::vector< double >	cracksHBHEdef_
const std::vector< double >	cracksHEHFdef_
const bool	debug_
const bool	doCracks_
std::map< DetId, std::vector< int > >	EcalAllDeadChannelsBitMap
std::map< DetId, EcalTrigTowerDetId >	EcalAllDeadChannelsTTMap
std::map< DetId, std::vector< double > >	EcalAllDeadChannelsValMap
const std::vector< double >	EcalDeadCellDeltaRFilterInput_
edm::ESHandle< EcalChannelStatus >	ecalStatus
edm::EventNumber_t	event
int	evtProcessedCnt
edm::ESHandle< CaloGeometry >	geometry
TH1F *	h1_dummy
edm::ESHandle< HcalChannelQuality >	hcalStatus
bool	isdata
bool	isPrintedOnce
const bool	isProd_
edm::Handle< edm::View< reco::Jet > >	jets
const std::vector< double >	jetSelCuts_
edm::EDGetTokenT< edm::View< reco::Jet > >	jetToken_
edm::LuminosityBlockNumber_t	ls
const bool	makeProfileRoot_
const int	maskedEcalChannelStatusThreshold_
edm::Handle< edm::View< reco::MET > >	met
edm::EDGetTokenT< edm::View< reco::MET > >	metToken_
double	pfmet
double	pfmetPhi
const bool	printSkimInfo_
TFile *	profFile
const std::string	profileRootName_
edm::RunNumber_t	run
const bool	taggingMode_
double	tcmet
double	tcmetPhi
int	totTPFilteredCnt
edm::ESHandle< EcalTrigTowerConstituentsMap >	ttMap_
const int	verbose_
double	wtdEvtProcessed
double	wtdTPFiltered

Additional Inherited Members
Public Types inherited from edm::EDFilter
typedef EDFilter	ModuleType
Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const *, const char *, edm::ProductResolverIndex >>
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::EDFilter
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &)
static bool	wantsGlobalLuminosityBlocks ()
static bool	wantsGlobalRuns ()
static bool	wantsStreamLuminosityBlocks ()
static bool	wantsStreamRuns ()
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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)
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

Description: <one line="" class="" summary>=""> Event filtering for RA2 analysis (filtering status is stored in the event)

Definition at line 96 of file EcalDeadCellDeltaRFilter.cc.

Constructor & Destructor Documentation

EcalDeadCellDeltaRFilter::EcalDeadCellDeltaRFilter ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 221 of file EcalDeadCellDeltaRFilter.cc.

References h1_dummy, makeProfileRoot_, profFile, and profileRootName_.

  : jetToken_ (consumes<edm::View<reco::Jet> >(iConfig.getParameter<edm::InputTag>("jetInputTag")))
  , jetSelCuts_ (iConfig.getParameter<std::vector<double> >("jetSelCuts"))

  , metToken_ (consumes<edm::View<reco::MET> >(iConfig.getParameter<edm::InputTag>("metInputTag")))

  , debug_ (iConfig.getUntrackedParameter<bool>("debug",false))
  , printSkimInfo_ (iConfig.getUntrackedParameter<bool>("printSkimInfo",false))

  , maskedEcalChannelStatusThreshold_ (iConfig.getParameter<int>("maskedEcalChannelStatusThreshold"))
  , chnStatusToBeEvaluated_ (iConfig.getParameter<int>("chnStatusToBeEvaluated"))

  , makeProfileRoot_ (iConfig.getUntrackedParameter<bool>("makeProfileRoot", true))
  , profileRootName_ (iConfig.getUntrackedParameter<std::string>("profileRootName", "EcalDeadCellDeltaRFilter.root"))

  , isProd_ (iConfig.getUntrackedParameter<bool>("isProd"))
  , verbose_ (iConfig.getParameter<int>("verbose"))

  , doCracks_ (iConfig.getUntrackedParameter<bool>("doCracks"))
  , cracksHBHEdef_ (iConfig.getParameter<std::vector<double> > ("cracksHBHEdef"))
  , cracksHEHFdef_ (iConfig.getParameter<std::vector<double> > ("cracksHEHFdef"))

  , EcalDeadCellDeltaRFilterInput_ (iConfig.getParameter<std::vector<double> >("EcalDeadCellDeltaRFilterInput"))

  , taggingMode_ (iConfig.getParameter<bool>("taggingMode"))
{
  produces<int> ("deadCellStatus"); produces<int> ("boundaryStatus");
  produces<bool>();

  if( makeProfileRoot_ ){
     profFile = new TFile(profileRootName_.c_str(), "RECREATE");
     h1_dummy = new TH1F("dummy", "dummy", 500, 0, 500);
  }
}

EcalDeadCellDeltaRFilter::~EcalDeadCellDeltaRFilter ( )

override

Definition at line 256 of file EcalDeadCellDeltaRFilter.cc.

References h1_dummy, makeProfileRoot_, and profFile.

                                                    {
  if( makeProfileRoot_ ){
     profFile->cd();

     h1_dummy->Write();

     profFile->Close();
     delete profFile;
  }
}

Member Function Documentation

void EcalDeadCellDeltaRFilter::beginJob ( void  )

overrideprivatevirtual

Reimplemented from edm::EDFilter.

Definition at line 335 of file EcalDeadCellDeltaRFilter.cc.

References gather_cfg::cout, and debug_.

                                        {
  if (debug_) std::cout << "beginJob" << std::endl;
}

void EcalDeadCellDeltaRFilter::beginLuminosityBlock ( const edm::LuminosityBlock &  iLSblock, const edm::EventSetup &  iSetup  )

overrideprivate

Definition at line 362 of file EcalDeadCellDeltaRFilter.cc.

References edm::EventSetup::get(), hcalStatus, and edm::ESHandleBase::isValid().

                                                                                                                   {

  // needs per-LS access, if used at all
  iSetup.get<HcalChannelQualityRcd>().get("withTopo",hcalStatus);
  if( !hcalStatus.isValid() )  throw "Failed to get HCAL channel status!";
  return; 
}

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

overrideprivate

Definition at line 345 of file EcalDeadCellDeltaRFilter.cc.

References gather_cfg::cout, debug_, EcalAllDeadChannelsBitMap, EcalAllDeadChannelsValMap, envSet(), and getChannelStatusMaps().

                                                                                      {
  if (debug_) std::cout << "beginRun" << std::endl;
// Channel status might change for each run (data)
// Event setup
  envSet(iSetup);
  getChannelStatusMaps();
  if( debug_) std::cout<< "EcalAllDeadChannelsValMap.size() : "<<EcalAllDeadChannelsValMap.size()<<"  EcalAllDeadChannelsBitMap.size() : "<<EcalAllDeadChannelsBitMap.size()<<std::endl;
  return;
}

int EcalDeadCellDeltaRFilter::dPhiToMETfunc ( const std::vector< reco::Jet > &  jetTVec, const double &  dPhiCutVal, std::vector< reco::Jet > &  closeToMETjetsVec  )

private

Definition at line 391 of file EcalDeadCellDeltaRFilter.cc.

References funct::abs(), reco::deltaPhi(), hiPixelPairStep_cff::deltaPhi, cuy::ii, metsig::jet, and reco::LeafCandidate::phi().

Referenced by filter().

                                                                                                                                             {

  closeToMETjetsVec.clear();

  double minDphi = 999.0;
  int minIdx = -1;
  for(unsigned int ii=0; ii<jetTVec.size(); ii++){

     const reco::Jet& jet = jetTVec[ii];

     double deltaPhi = std::abs(reco::deltaPhi( jet.phi(), (*met)[0].phi() ) );
     if( deltaPhi > dPhiCutVal ) continue;

     closeToMETjetsVec.push_back(jetTVec[ii]);

     if( deltaPhi < minDphi ){
        minDphi = deltaPhi;
        minIdx = ii;
     }
  }

  if( minIdx == -1 ){} // removing a stupid compiling WARNING that minIdx NOT used.
//  if( minIdx == -1 ) return 0;
//  closeToMETjetsVec.push_back(jetTVec[minIdx]);

  return (int)closeToMETjetsVec.size();
}

int EcalDeadCellDeltaRFilter::dRtoMaskedChnsEvtFilterFunc ( const std::vector< reco::Jet > &  jetTVec, const int &  chnStatus, const double &  dRCutVal  )

private

Definition at line 420 of file EcalDeadCellDeltaRFilter.cc.

References cuy::ii, isCloseToBadEcalChannel(), and metsig::jet.

Referenced by filter().

                                                                                                                                        {

  int isClose = 0;

  for(unsigned int ii=0; ii<jetTVec.size(); ii++){

     const reco::Jet& jet = jetTVec[ii];

     std::map<double, DetId> dummy;
     int isPerJetClose = isCloseToBadEcalChannel(jet, dRCutVal, chnStatus, dummy);
//     if( isPerJetClose ){ isClose = 1; break; }
     if( isPerJetClose ){ isClose ++; }
  }

  return isClose;

}

void EcalDeadCellDeltaRFilter::endJob ( void  )

overrideprivatevirtual

Reimplemented from edm::EDFilter.

Definition at line 340 of file EcalDeadCellDeltaRFilter.cc.

References gather_cfg::cout, and debug_.

                                      {
  if (debug_) std::cout << "endJob" << std::endl;
}

void EcalDeadCellDeltaRFilter::endRun ( const edm::Run &  run, const edm::EventSetup &  iSetup  )

overrideprivate

Definition at line 356 of file EcalDeadCellDeltaRFilter.cc.

References gather_cfg::cout, and debug_.

                                                                                    {
  if (debug_) std::cout << "endRun" << std::endl;
  return;
}

void EcalDeadCellDeltaRFilter::envSet ( const edm::EventSetup &  iSetup )

privatevirtual

Definition at line 267 of file EcalDeadCellDeltaRFilter.cc.

References gather_cfg::cout, debug_, ecalStatus, edm::EventSetup::get(), and ttMap_.

Referenced by beginRun().

                                                                 {

  if (debug_) std::cout << "***envSet***" << std::endl;

  iSetup.get<IdealGeometryRecord>().get(ttMap_);

  iSetup.get<EcalChannelStatusRcd> ().get(ecalStatus);
  iSetup.get<CaloGeometryRecord>   ().get(geometry);

  if( !ecalStatus.isValid() )  throw "Failed to get ECAL channel status!";
  if( !geometry.isValid()   )  throw "Failed to get the geometry!";

}

int EcalDeadCellDeltaRFilter::etaToBoundary ( const std::vector< reco::Jet > &  jetTVec )

private

Definition at line 370 of file EcalDeadCellDeltaRFilter.cc.

References funct::abs(), cracksHBHEdef_, cracksHEHFdef_, and funct::pow().

Referenced by filter().

                                                                            {

  int isClose = 0;

  int cntOrder10 = 0;
  for(unsigned int ij=0; ij<jetTVec.size(); ij++){

     double recoJetEta = jetTVec[ij].eta();

     if( std::abs(recoJetEta)>cracksHBHEdef_[0] && std::abs(recoJetEta)<cracksHBHEdef_[1] ) isClose += (cntOrder10*10 + 1);
     if( std::abs(recoJetEta)>cracksHEHFdef_[0] && std::abs(recoJetEta)<cracksHEHFdef_[1] ) isClose += (cntOrder10*10 + 2);

     if( isClose/pow(10, cntOrder10) >=3 ) cntOrder10 = isClose/10 + 1;
  }

  return isClose;

}

bool EcalDeadCellDeltaRFilter::filter ( edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

overrideprivate

Definition at line 282 of file EcalDeadCellDeltaRFilter.cc.

References funct::abs(), chnStatusToBeEvaluated_, debug_, doCracks_, dPhiToMETfunc(), dRtoMaskedChnsEvtFilterFunc(), EcalDeadCellDeltaRFilterInput_, etaToBoundary(), jets, jetSelCuts_, loadEventInfo(), loadJets(), loadMET(), ls, makeProfileRoot_, met, phi, EnergyCorrector::pt, edm::Event::put(), run, and taggingMode_.

                                                                                   {

  loadEventInfo(iEvent, iSetup);
  loadJets(iEvent, iSetup);
  loadMET(iEvent, iSetup);

// XXX: In the following, never assign pass to true again
// Currently, always true
  bool pass = true;

  using namespace edm;

  std::vector<reco::Jet> seledJets;

  for( edm::View<reco::Jet>::const_iterator ij = jets->begin(); ij != jets->end(); ij++){
     if( ij->pt() > jetSelCuts_[0] && std::abs(ij->eta()) < jetSelCuts_[1] ){
        seledJets.push_back(reco::Jet(*ij));
     }
  }

  if( seledJets.empty() ) return pass;

  double dPhiToMET = EcalDeadCellDeltaRFilterInput_[0], dRtoDeadCell = EcalDeadCellDeltaRFilterInput_[1];

  std::vector<reco::Jet> closeToMETjetsVec;

  int dPhiToMETstatus = dPhiToMETfunc(seledJets, dPhiToMET, closeToMETjetsVec);

// Get event filter for simple dR cut
  int deadCellStatus = dRtoMaskedChnsEvtFilterFunc(closeToMETjetsVec, chnStatusToBeEvaluated_, dRtoDeadCell);

  int boundaryStatus = etaToBoundary(closeToMETjetsVec);

  if(debug_ ){
     printf("\nrun : %8u  event : %12llu  ls : %8u  dPhiToMETstatus : %d  deadCellStatus : %d  boundaryStatus : %d\n", run, event, ls, dPhiToMETstatus, deadCellStatus, boundaryStatus);
     printf("met : %6.2f  metphi : % 6.3f  dPhiToMET : %5.3f  dRtoDeadCell : %5.3f\n", (*met)[0].pt(), (*met)[0].phi(), dPhiToMET, dRtoDeadCell);
  }

  if( makeProfileRoot_ ){
//     h1_dummy->Fill(xxx);
  }

  iEvent.put(std::make_unique<int>(deadCellStatus), "deadCellStatus");
  iEvent.put(std::make_unique<int>(boundaryStatus), "boundaryStatus");

  if( deadCellStatus || (doCracks_ && boundaryStatus) ) pass = false;

  iEvent.put(std::make_unique<bool>(pass));

  return taggingMode_ || pass;
}

int EcalDeadCellDeltaRFilter::getChannelStatusMaps ( )

private

Definition at line 476 of file EcalDeadCellDeltaRFilter.cc.

References DEFINE_FWK_MODULE, EcalAllDeadChannelsBitMap, EcalAllDeadChannelsTTMap, EcalAllDeadChannelsValMap, ecalStatus, PVValHelper::eta, EBDetId::ETAPHIMODE, CaloSubdetectorGeometry::getGeometry(), maskedEcalChannelStatusThreshold_, phi, mps_update::status, theta(), EcalTrigTowerConstituentsMap::towerOf(), ttMap_, EBDetId::validDetId(), EEDetId::validDetId(), and EEDetId::XYMODE.

Referenced by beginRun().

                                                  {

  EcalAllDeadChannelsValMap.clear(); EcalAllDeadChannelsBitMap.clear();

// Loop over EB ...
  for( int ieta=-85; ieta<=85; ieta++ ){
     for( int iphi=0; iphi<=360; iphi++ ){
        if(! EBDetId::validDetId( ieta, iphi ) )  continue;

        const EBDetId detid = EBDetId( ieta, iphi, EBDetId::ETAPHIMODE );
        EcalChannelStatus::const_iterator chit = ecalStatus->find( detid );
// refer https://twiki.cern.ch/twiki/bin/viewauth/CMS/EcalChannelStatus
        int status = ( chit != ecalStatus->end() ) ? chit->getStatusCode() & 0x1F : -1;

    const CaloSubdetectorGeometry* subGeom = geometry->getSubdetectorGeometry (detid);
        auto cellGeom = subGeom->getGeometry (detid);
        double eta = cellGeom->getPosition ().eta ();
        double phi = cellGeom->getPosition ().phi ();
        double theta = cellGeom->getPosition().theta();

        if(status >= maskedEcalChannelStatusThreshold_){
           std::vector<double> valVec; std::vector<int> bitVec;
           valVec.push_back(eta); valVec.push_back(phi); valVec.push_back(theta);
           bitVec.push_back(1); bitVec.push_back(ieta); bitVec.push_back(iphi); bitVec.push_back(status);
           EcalAllDeadChannelsValMap.insert( std::make_pair(detid, valVec) );
           EcalAllDeadChannelsBitMap.insert( std::make_pair(detid, bitVec) );
        }
     } // end loop iphi
  } // end loop ieta

// Loop over EE detid
  for( int ix=0; ix<=100; ix++ ){
     for( int iy=0; iy<=100; iy++ ){
        for( int iz=-1; iz<=1; iz++ ){
           if(iz==0)  continue;
           if(! EEDetId::validDetId( ix, iy, iz ) )  continue;

           const EEDetId detid = EEDetId( ix, iy, iz, EEDetId::XYMODE );
           EcalChannelStatus::const_iterator chit = ecalStatus->find( detid );
           int status = ( chit != ecalStatus->end() ) ? chit->getStatusCode() & 0x1F : -1;

           const CaloSubdetectorGeometry* subGeom = geometry->getSubdetectorGeometry (detid);
           auto cellGeom = subGeom->getGeometry (detid);
           double eta = cellGeom->getPosition ().eta () ;
           double phi = cellGeom->getPosition ().phi () ;
           double theta = cellGeom->getPosition().theta();

           if(status >= maskedEcalChannelStatusThreshold_){
              std::vector<double> valVec; std::vector<int> bitVec;
              valVec.push_back(eta); valVec.push_back(phi); valVec.push_back(theta);
              bitVec.push_back(2); bitVec.push_back(ix); bitVec.push_back(iy); bitVec.push_back(iz); bitVec.push_back(status);
              EcalAllDeadChannelsValMap.insert( std::make_pair(detid, valVec) );
              EcalAllDeadChannelsBitMap.insert( std::make_pair(detid, bitVec) );
           }
        } // end loop iz
     } // end loop iy
  } // end loop ix

  EcalAllDeadChannelsTTMap.clear();
  std::map<DetId, std::vector<int> >::iterator bitItor;
  for(bitItor = EcalAllDeadChannelsBitMap.begin(); bitItor != EcalAllDeadChannelsBitMap.end(); bitItor++){
     const DetId id = bitItor->first;
     EcalTrigTowerDetId ttDetId = ttMap_->towerOf(id);
     EcalAllDeadChannelsTTMap.insert(std::make_pair(id, ttDetId) );
  }

  return 1;
}

int EcalDeadCellDeltaRFilter::isCloseToBadEcalChannel ( const reco::Jet &  jet, const double &  deltaRCut, const int &  chnStatus, std::map< double, DetId > &  deltaRdetIdMap  )

private

Definition at line 439 of file EcalDeadCellDeltaRFilter.cc.

References funct::abs(), gather_cfg::cout, reco::deltaR(), EcalAllDeadChannelsBitMap, EcalAllDeadChannelsValMap, PVValHelper::eta, reco::LeafCandidate::eta(), reco::btau::jetEta, reco::btau::jetPhi, phi, reco::LeafCandidate::phi(), and mps_update::status.

Referenced by dRtoMaskedChnsEvtFilterFunc().

                                                                                                                                                             {

   double jetEta = jet.eta(), jetPhi = jet.phi();

   deltaRdetIdMap.clear();

   double min_dist = 999;
   DetId min_detId;

   std::map<DetId, std::vector<int> >::iterator bitItor;
   for(bitItor = EcalAllDeadChannelsBitMap.begin(); bitItor != EcalAllDeadChannelsBitMap.end(); bitItor++){

      DetId maskedDetId = bitItor->first;
//      int subdet = bitItor->second.front();
      int status = bitItor->second.back();

      if( chnStatus >0 && status != chnStatus ) continue;
      if( chnStatus <0 && status < abs(chnStatus) ) continue;

      std::map<DetId, std::vector<double> >::iterator valItor = EcalAllDeadChannelsValMap.find(maskedDetId);
      if( valItor == EcalAllDeadChannelsValMap.end() ){ std::cout<<"Error cannot find maskedDetId in EcalAllDeadChannelsValMap ?!"<<std::endl; continue; }

      double eta = (valItor->second)[0], phi = (valItor->second)[1];

      double dist = reco::deltaR(eta, phi, jetEta, jetPhi);

      if( min_dist > dist ){ min_dist = dist; min_detId = maskedDetId; }
   }

   if( min_dist > deltaRCut && deltaRCut >0 ) return 0;

   deltaRdetIdMap.insert( std::make_pair(min_dist, min_detId) );

   return 1;
}

void EcalDeadCellDeltaRFilter::loadEventInfo ( const edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

private

Definition at line 192 of file EcalDeadCellDeltaRFilter.cc.

References gather_cfg::cout, debug_, edm::EventID::event(), edm::EventBase::id(), isdata, isPrintedOnce, edm::EventBase::isRealData(), ls, edm::EventBase::luminosityBlock(), edm::EventID::run(), and run.

Referenced by filter().

                                                                                               {
   run = iEvent.id().run();
   event = iEvent.id().event();
   ls = iEvent.luminosityBlock();
   isdata = iEvent.isRealData();

   if( !isPrintedOnce ){
      if( debug_ ){
         if( isdata ) std::cout<<"\nInput dataset is DATA"<<std::endl<<std::endl;
         else std::cout<<"\nInput dataset is MC"<<std::endl<<std::endl;
      }
      isPrintedOnce = true;
   }

}

void EcalDeadCellDeltaRFilter::loadJets ( const edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

private

Definition at line 208 of file EcalDeadCellDeltaRFilter.cc.

References edm::Event::getByToken(), jets, and jetToken_.

Referenced by filter().

                                                                                           {

  iEvent.getByToken(jetToken_, jets);

}

void EcalDeadCellDeltaRFilter::loadMET ( const edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

private

Definition at line 186 of file EcalDeadCellDeltaRFilter.cc.

References edm::Event::getByToken(), met, and metToken_.

Referenced by filter().

                                                                                         {

  iEvent.getByToken(metToken_, met);

}

Member Data Documentation

double EcalDeadCellDeltaRFilter::calomet

private

Definition at line 134 of file EcalDeadCellDeltaRFilter.cc.

double EcalDeadCellDeltaRFilter::calometPhi

private

Definition at line 134 of file EcalDeadCellDeltaRFilter.cc.

const int EcalDeadCellDeltaRFilter::chnStatusToBeEvaluated_

private

Definition at line 144 of file EcalDeadCellDeltaRFilter.cc.

Referenced by filter().

const std::vector<double> EcalDeadCellDeltaRFilter::cracksHBHEdef_

private

Definition at line 171 of file EcalDeadCellDeltaRFilter.cc.

Referenced by etaToBoundary().

const std::vector<double> EcalDeadCellDeltaRFilter::cracksHEHFdef_

private

Definition at line 171 of file EcalDeadCellDeltaRFilter.cc.

Referenced by etaToBoundary().

const bool EcalDeadCellDeltaRFilter::debug_

private

Definition at line 121 of file EcalDeadCellDeltaRFilter.cc.

Referenced by beginJob(), beginRun(), endJob(), endRun(), envSet(), filter(), and loadEventInfo().

const bool EcalDeadCellDeltaRFilter::doCracks_

private

Definition at line 169 of file EcalDeadCellDeltaRFilter.cc.

Referenced by filter().

std::map<DetId, std::vector<int> > EcalDeadCellDeltaRFilter::EcalAllDeadChannelsBitMap

private

Definition at line 151 of file EcalDeadCellDeltaRFilter.cc.

Referenced by beginRun(), getChannelStatusMaps(), and isCloseToBadEcalChannel().

std::map<DetId, EcalTrigTowerDetId> EcalDeadCellDeltaRFilter::EcalAllDeadChannelsTTMap

private

Definition at line 154 of file EcalDeadCellDeltaRFilter.cc.

Referenced by getChannelStatusMaps().

std::map<DetId, std::vector<double> > EcalDeadCellDeltaRFilter::EcalAllDeadChannelsValMap

private

Definition at line 148 of file EcalDeadCellDeltaRFilter.cc.

Referenced by beginRun(), getChannelStatusMaps(), and isCloseToBadEcalChannel().

const std::vector<double> EcalDeadCellDeltaRFilter::EcalDeadCellDeltaRFilterInput_

private

Definition at line 174 of file EcalDeadCellDeltaRFilter.cc.

Referenced by filter().

edm::ESHandle<EcalChannelStatus> EcalDeadCellDeltaRFilter::ecalStatus

private

Definition at line 137 of file EcalDeadCellDeltaRFilter.cc.

Referenced by envSet(), and getChannelStatusMaps().

edm::EventNumber_t EcalDeadCellDeltaRFilter::event

private

Definition at line 130 of file EcalDeadCellDeltaRFilter.cc.

Referenced by Types.EventID::cppID().

int EcalDeadCellDeltaRFilter::evtProcessedCnt

private

Definition at line 158 of file EcalDeadCellDeltaRFilter.cc.

edm::ESHandle<CaloGeometry> EcalDeadCellDeltaRFilter::geometry

private

Definition at line 139 of file EcalDeadCellDeltaRFilter.cc.

Referenced by Vispa.Gui.ConnectableWidget.ConnectableWidget::leaveEvent().

TH1F* EcalDeadCellDeltaRFilter::h1_dummy

private

Definition at line 164 of file EcalDeadCellDeltaRFilter.cc.

Referenced by EcalDeadCellDeltaRFilter(), and ~EcalDeadCellDeltaRFilter().

edm::ESHandle<HcalChannelQuality> EcalDeadCellDeltaRFilter::hcalStatus

private

Definition at line 138 of file EcalDeadCellDeltaRFilter.cc.

Referenced by beginLuminosityBlock().

bool EcalDeadCellDeltaRFilter::isdata

private

Definition at line 132 of file EcalDeadCellDeltaRFilter.cc.

Referenced by loadEventInfo().

bool EcalDeadCellDeltaRFilter::isPrintedOnce

private

Definition at line 123 of file EcalDeadCellDeltaRFilter.cc.

Referenced by loadEventInfo().

const bool EcalDeadCellDeltaRFilter::isProd_

private

Definition at line 166 of file EcalDeadCellDeltaRFilter.cc.

edm::Handle<edm::View<reco::Jet> > EcalDeadCellDeltaRFilter::jets

private

Definition at line 113 of file EcalDeadCellDeltaRFilter.cc.

Referenced by filter(), and loadJets().

const std::vector<double> EcalDeadCellDeltaRFilter::jetSelCuts_

private

Definition at line 116 of file EcalDeadCellDeltaRFilter.cc.

Referenced by filter().

edm::EDGetTokenT<edm::View<reco::Jet> > EcalDeadCellDeltaRFilter::jetToken_

private

Definition at line 112 of file EcalDeadCellDeltaRFilter.cc.

Referenced by loadJets().

edm::LuminosityBlockNumber_t EcalDeadCellDeltaRFilter::ls

private

Definition at line 131 of file EcalDeadCellDeltaRFilter.cc.

Referenced by filter(), dqm_interfaces.DQMcommunicator::get_root_objects(), dqm_interfaces.DQMcommunicator::get_root_objects_list(), dqm_interfaces.DQMcommunicator::get_root_objects_list_recursive(), dqm_interfaces.DQMcommunicator::get_root_objects_names_list_recursive(), dqm_interfaces.DQMcommunicator::get_root_objects_recursive(), loadEventInfo(), and dqm_interfaces.DirWalkerFile::walk().

const bool EcalDeadCellDeltaRFilter::makeProfileRoot_

private

Definition at line 161 of file EcalDeadCellDeltaRFilter.cc.

Referenced by EcalDeadCellDeltaRFilter(), filter(), and ~EcalDeadCellDeltaRFilter().

const int EcalDeadCellDeltaRFilter::maskedEcalChannelStatusThreshold_

private

Definition at line 143 of file EcalDeadCellDeltaRFilter.cc.

Referenced by getChannelStatusMaps().

edm::Handle<edm::View<reco::MET> > EcalDeadCellDeltaRFilter::met

private

Definition at line 119 of file EcalDeadCellDeltaRFilter.cc.

Referenced by filter(), and loadMET().

edm::EDGetTokenT<edm::View<reco::MET> > EcalDeadCellDeltaRFilter::metToken_

private

Definition at line 118 of file EcalDeadCellDeltaRFilter.cc.

Referenced by loadMET().

double EcalDeadCellDeltaRFilter::pfmet

private

Definition at line 134 of file EcalDeadCellDeltaRFilter.cc.

double EcalDeadCellDeltaRFilter::pfmetPhi

private

Definition at line 134 of file EcalDeadCellDeltaRFilter.cc.

const bool EcalDeadCellDeltaRFilter::printSkimInfo_

private

Definition at line 121 of file EcalDeadCellDeltaRFilter.cc.

TFile* EcalDeadCellDeltaRFilter::profFile

private

Definition at line 163 of file EcalDeadCellDeltaRFilter.cc.

Referenced by EcalDeadCellDeltaRFilter(), and ~EcalDeadCellDeltaRFilter().

const std::string EcalDeadCellDeltaRFilter::profileRootName_

private

Definition at line 162 of file EcalDeadCellDeltaRFilter.cc.

Referenced by EcalDeadCellDeltaRFilter().

edm::RunNumber_t EcalDeadCellDeltaRFilter::run

private

Definition at line 129 of file EcalDeadCellDeltaRFilter.cc.

Referenced by DTWorkflow.DTWorkflow::all(), Types.EventID::cppID(), Types.LuminosityBlockID::cppID(), o2olib.O2OTool::execute(), filter(), loadEventInfo(), and dqm_interfaces.DirWalkerFile::walk().

const bool EcalDeadCellDeltaRFilter::taggingMode_

private

Definition at line 183 of file EcalDeadCellDeltaRFilter.cc.

Referenced by filter().

double EcalDeadCellDeltaRFilter::tcmet

private

Definition at line 134 of file EcalDeadCellDeltaRFilter.cc.

double EcalDeadCellDeltaRFilter::tcmetPhi

private

Definition at line 134 of file EcalDeadCellDeltaRFilter.cc.

int EcalDeadCellDeltaRFilter::totTPFilteredCnt

private

Definition at line 158 of file EcalDeadCellDeltaRFilter.cc.

edm::ESHandle<EcalTrigTowerConstituentsMap> EcalDeadCellDeltaRFilter::ttMap_

private

Definition at line 141 of file EcalDeadCellDeltaRFilter.cc.

Referenced by envSet(), and getChannelStatusMaps().

const int EcalDeadCellDeltaRFilter::verbose_

private

Definition at line 167 of file EcalDeadCellDeltaRFilter.cc.

double EcalDeadCellDeltaRFilter::wtdEvtProcessed

private

Definition at line 159 of file EcalDeadCellDeltaRFilter.cc.

double EcalDeadCellDeltaRFilter::wtdTPFiltered

private

Definition at line 159 of file EcalDeadCellDeltaRFilter.cc.