#include <EcalDeadCellDeltaRFilter.cc>

Inheritance diagram for EcalDeadCellDeltaRFilter:

Public Member Functions
	EcalDeadCellDeltaRFilter (const edm::ParameterSet &)

	~EcalDeadCellDeltaRFilter ()

Public Member Functions inherited from edm::EDFilter
	EDFilter ()

ModuleDescription const &	moduleDescription () const

virtual	~EDFilter ()

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 (const std::string &iProcessName, std::vector< const char * > &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 Member Functions
virtual void	beginJob () override

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

virtual void	endJob () override

virtual void	endRun (const edm::Run &, const edm::EventSetup &) override

virtual void	envSet (const edm::EventSetup &)

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

virtual 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_

EcalTPGScale	ecalScale_

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

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

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

Definition at line 97 of file EcalDeadCellDeltaRFilter.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 222 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 ( )

Definition at line 257 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 342 of file EcalDeadCellDeltaRFilter.cc.

References gather_cfg::cout, and debug_.

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

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

overrideprivatevirtual

Reimplemented from edm::EDFilter.

Definition at line 352 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 390 of file EcalDeadCellDeltaRFilter.cc.

References funct::abs(), SiPixelRawToDigiRegional_cfi::deltaPhi, reco::deltaPhi(), cuy::ii, metsig::jet, HLT_25ns14e33_v1_cff::minDphi, 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 419 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 347 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
	)

overrideprivatevirtual

Reimplemented from edm::EDFilter.

Definition at line 363 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 268 of file EcalDeadCellDeltaRFilter.cc.

References gather_cfg::cout, debug_, ecalScale_, ecalStatus, geometry, edm::EventSetup::get(), hcalStatus, edm::ESHandleBase::isValid(), EcalTPGScale::setEventSetup(), and ttMap_.

Referenced by beginRun().

                                                                  {
 
   if (debug_) std::cout << "***envSet***" << std::endl;
 
   ecalScale_.setEventSetup( iSetup );
   iSetup.get<IdealGeometryRecord>().get(ttMap_);
 
   iSetup.get<EcalChannelStatusRcd> ().get(ecalStatus);
   iSetup.get<HcalChannelQualityRcd>().get("withTopo",hcalStatus);
   iSetup.get<CaloGeometryRecord>   ().get(geometry);
 
   if( !ecalStatus.isValid() )  throw "Failed to get ECAL channel status!";
   if( !hcalStatus.isValid() )  throw "Failed to get HCAL channel status!";
   if( !geometry.isValid()   )  throw "Failed to get the geometry!";
 
 }

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

private

Definition at line 369 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
	)

overrideprivatevirtual

Implements edm::EDFilter.

Definition at line 286 of file EcalDeadCellDeltaRFilter.cc.

References funct::abs(), edm::View< T >::begin(), chnStatusToBeEvaluated_, debug_, doCracks_, dPhiToMETfunc(), dRtoMaskedChnsEvtFilterFunc(), EcalDeadCellDeltaRFilterInput_, edm::View< T >::end(), etaToBoundary(), event, 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);
   }
 
   std::auto_ptr<int> deadCellStatusPtr ( new int(deadCellStatus) );
   std::auto_ptr<int> boundaryStatusPtr ( new int(boundaryStatus) );
 
   iEvent.put( deadCellStatusPtr, "deadCellStatus");
   iEvent.put( boundaryStatusPtr, "boundaryStatus");
 
   if( deadCellStatus || (doCracks_ && boundaryStatus) ) pass = false;
 
   iEvent.put( std::auto_ptr<bool>(new bool(pass)) );
 
   return taggingMode_ || pass;
 }

int EcalDeadCellDeltaRFilter::getChannelStatusMaps ( )

private

Definition at line 475 of file EcalDeadCellDeltaRFilter.cc.

References cond::rpcobgas::detid, EcalAllDeadChannelsBitMap, EcalAllDeadChannelsTTMap, EcalAllDeadChannelsValMap, ecalStatus, eta, PV3DBase< T, PVType, FrameType >::eta(), EBDetId::ETAPHIMODE, geometry, CaloSubdetectorGeometry::getGeometry(), CaloCellGeometry::getPosition(), maskedEcalChannelStatusThreshold_, phi, PV3DBase< T, PVType, FrameType >::phi(), ntuplemaker::status, PV3DBase< T, PVType, FrameType >::theta(), theta(), 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);
         const CaloCellGeometry*        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);
            const CaloCellGeometry*        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 438 of file EcalDeadCellDeltaRFilter.cc.

References funct::abs(), gather_cfg::cout, reco::deltaR(), EcalAllDeadChannelsBitMap, EcalAllDeadChannelsValMap, eta, reco::LeafCandidate::eta(), reco::btau::jetEta, reco::btau::jetPhi, phi, reco::LeafCandidate::phi(), and ntuplemaker::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 193 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 209 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 187 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 133 of file EcalDeadCellDeltaRFilter.cc.

double EcalDeadCellDeltaRFilter::calometPhi

private

Definition at line 133 of file EcalDeadCellDeltaRFilter.cc.

const int EcalDeadCellDeltaRFilter::chnStatusToBeEvaluated_

private

Definition at line 145 of file EcalDeadCellDeltaRFilter.cc.

Referenced by filter().

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

private

Definition at line 172 of file EcalDeadCellDeltaRFilter.cc.

Referenced by etaToBoundary().

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

private

Definition at line 172 of file EcalDeadCellDeltaRFilter.cc.

Referenced by etaToBoundary().

const bool EcalDeadCellDeltaRFilter::debug_

private

Definition at line 120 of file EcalDeadCellDeltaRFilter.cc.

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

const bool EcalDeadCellDeltaRFilter::doCracks_

private

Definition at line 170 of file EcalDeadCellDeltaRFilter.cc.

Referenced by filter().

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

private

Definition at line 152 of file EcalDeadCellDeltaRFilter.cc.

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

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

private

Definition at line 155 of file EcalDeadCellDeltaRFilter.cc.

Referenced by getChannelStatusMaps().

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

private

Definition at line 149 of file EcalDeadCellDeltaRFilter.cc.

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

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

private

Definition at line 175 of file EcalDeadCellDeltaRFilter.cc.

Referenced by filter().

EcalTPGScale EcalDeadCellDeltaRFilter::ecalScale_

private

Definition at line 142 of file EcalDeadCellDeltaRFilter.cc.

Referenced by envSet().

edm::ESHandle<EcalChannelStatus> EcalDeadCellDeltaRFilter::ecalStatus

private

Definition at line 136 of file EcalDeadCellDeltaRFilter.cc.

Referenced by envSet(), and getChannelStatusMaps().

edm::EventNumber_t EcalDeadCellDeltaRFilter::event

private

Definition at line 129 of file EcalDeadCellDeltaRFilter.cc.

Referenced by Types.EventID::cppID(), and filter().

int EcalDeadCellDeltaRFilter::evtProcessedCnt

private

Definition at line 159 of file EcalDeadCellDeltaRFilter.cc.

edm::ESHandle<CaloGeometry> EcalDeadCellDeltaRFilter::geometry

private

Definition at line 138 of file EcalDeadCellDeltaRFilter.cc.

Referenced by envSet(), getChannelStatusMaps(), and Vispa.Gui.ConnectableWidget.ConnectableWidget::leaveEvent().

TH1F* EcalDeadCellDeltaRFilter::h1_dummy

private

Definition at line 165 of file EcalDeadCellDeltaRFilter.cc.

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

edm::ESHandle<HcalChannelQuality> EcalDeadCellDeltaRFilter::hcalStatus

private

Definition at line 137 of file EcalDeadCellDeltaRFilter.cc.

Referenced by envSet().

bool EcalDeadCellDeltaRFilter::isdata

private

Definition at line 131 of file EcalDeadCellDeltaRFilter.cc.

Referenced by loadEventInfo().

bool EcalDeadCellDeltaRFilter::isPrintedOnce

private

Definition at line 122 of file EcalDeadCellDeltaRFilter.cc.

Referenced by loadEventInfo().

const bool EcalDeadCellDeltaRFilter::isProd_

private

Definition at line 167 of file EcalDeadCellDeltaRFilter.cc.

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

private

Definition at line 112 of file EcalDeadCellDeltaRFilter.cc.

Referenced by filter(), and loadJets().

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

private

Definition at line 115 of file EcalDeadCellDeltaRFilter.cc.

Referenced by filter().

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

private

Definition at line 111 of file EcalDeadCellDeltaRFilter.cc.

Referenced by loadJets().

edm::LuminosityBlockNumber_t EcalDeadCellDeltaRFilter::ls

private

Definition at line 130 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 162 of file EcalDeadCellDeltaRFilter.cc.

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

const int EcalDeadCellDeltaRFilter::maskedEcalChannelStatusThreshold_

private

Definition at line 144 of file EcalDeadCellDeltaRFilter.cc.

Referenced by getChannelStatusMaps().

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

private

Definition at line 118 of file EcalDeadCellDeltaRFilter.cc.

Referenced by filter(), and loadMET().

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

private

Definition at line 117 of file EcalDeadCellDeltaRFilter.cc.

Referenced by loadMET().

double EcalDeadCellDeltaRFilter::pfmet

private

Definition at line 133 of file EcalDeadCellDeltaRFilter.cc.

double EcalDeadCellDeltaRFilter::pfmetPhi

private

Definition at line 133 of file EcalDeadCellDeltaRFilter.cc.

const bool EcalDeadCellDeltaRFilter::printSkimInfo_

private

Definition at line 120 of file EcalDeadCellDeltaRFilter.cc.

TFile* EcalDeadCellDeltaRFilter::profFile

private

Definition at line 164 of file EcalDeadCellDeltaRFilter.cc.

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

const std::string EcalDeadCellDeltaRFilter::profileRootName_

private

Definition at line 163 of file EcalDeadCellDeltaRFilter.cc.

Referenced by EcalDeadCellDeltaRFilter().

edm::RunNumber_t EcalDeadCellDeltaRFilter::run

private

Definition at line 128 of file EcalDeadCellDeltaRFilter.cc.

Referenced by Types.EventID::cppID(), Types.LuminosityBlockID::cppID(), filter(), loadEventInfo(), and dqm_interfaces.DirWalkerFile::walk().

const bool EcalDeadCellDeltaRFilter::taggingMode_

private

Definition at line 184 of file EcalDeadCellDeltaRFilter.cc.

Referenced by filter().

double EcalDeadCellDeltaRFilter::tcmet

private

Definition at line 133 of file EcalDeadCellDeltaRFilter.cc.

double EcalDeadCellDeltaRFilter::tcmetPhi

private

Definition at line 133 of file EcalDeadCellDeltaRFilter.cc.

int EcalDeadCellDeltaRFilter::totTPFilteredCnt

private

Definition at line 159 of file EcalDeadCellDeltaRFilter.cc.

edm::ESHandle<EcalTrigTowerConstituentsMap> EcalDeadCellDeltaRFilter::ttMap_

private

Definition at line 140 of file EcalDeadCellDeltaRFilter.cc.

Referenced by envSet(), and getChannelStatusMaps().

const int EcalDeadCellDeltaRFilter::verbose_

private

Definition at line 168 of file EcalDeadCellDeltaRFilter.cc.

double EcalDeadCellDeltaRFilter::wtdEvtProcessed

private

Definition at line 160 of file EcalDeadCellDeltaRFilter.cc.

double EcalDeadCellDeltaRFilter::wtdTPFiltered

private

Definition at line 160 of file EcalDeadCellDeltaRFilter.cc.

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation