EEBadScFilter Class Reference

Inheritance diagram for EEBadScFilter:

Public Member Functions
	EEBadScFilter (const edm::ParameterSet &iConfig)
	~EEBadScFilter ()
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 bool	filter (edm::Event &iEvent, const edm::EventSetup &iSetup) override
virtual void	scan5x5 (const DetId &det, const edm::Handle< EcalRecHitCollection > &hits, const edm::ESHandle< CaloTopology > &caloTopo, const edm::ESHandle< CaloGeometry > &geometry, int &nHits, float &totEt)

Private Attributes
const std::vector< int >	badsc_
const bool	debug_
const edm::EDGetTokenT < EcalRecHitCollection >	eeRHSrcToken_
const double	Emin_
const double	EtminSC_
const int	nBadHitsSC_
const int	side_
const bool	taggingMode_

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

Definition at line 37 of file EEBadScFilter.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 72 of file EEBadScFilter.cc.

  : eeRHSrcToken_     (consumes<EcalRecHitCollection>(iConfig.getParameter<edm::InputTag>("EERecHitSource")))
  , Emin_        (iConfig.getParameter<double>("EminHit"))
  , EtminSC_     (iConfig.getParameter<double>("EtminSC"))
  , side_        (iConfig.getParameter<int>("SCsize"))
  , nBadHitsSC_  (iConfig.getParameter<int>("nBadHitsSC"))
  , badsc_       (iConfig.getParameter<std::vector<int> >("badscEE"))
  , taggingMode_ (iConfig.getParameter<bool>("taggingMode"))
  , debug_       (iConfig.getParameter<bool>("debug"))
{
  produces<bool>();
}

EEBadScFilter::~EEBadScFilter ( )

inline

Definition at line 42 of file EEBadScFilter.cc.

{}

Member Function Documentation

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

overrideprivatevirtual

Implements edm::EDFilter.

Definition at line 136 of file EEBadScFilter.cc.

References badsc_, gather_cfg::cout, debug_, eeRHSrcToken_, EtminSC_, edm::EventSetup::get(), edm::Event::getByToken(), nBadHitsSC_, nhits, edm::Event::put(), scan5x5(), and taggingMode_.

                                                                          {


  // load required collections


  // EE rechit collection
  edm::Handle<EcalRecHitCollection> eeRHs;
  iEvent.getByToken(eeRHSrcToken_, eeRHs);

  // Calo Geometry - needed for computing E_t
  edm::ESHandle<CaloGeometry> pG;
  iSetup.get<CaloGeometryRecord>().get(pG);

  // Calo Toplology - needed for navigating the 5x5 xtal array around the centre of a SC
  edm::ESHandle<CaloTopology> pTopology;
  iSetup.get<CaloTopologyRecord>().get(pTopology);

  // by default the event is OK
  bool pass = true;


  // set discriminating variables to zero
  int nhits=0;
  float totEt=0.0;


  // define detid variables and ix,iy,iz coordinates

  EEDetId scdet;
  DetId det;

  int ix,iy,iz;

  ix=0,iy=0,iz=0;


  // loop over the list of bad SCs (defined in the python file)

  for (std::vector<int>::const_iterator scit = badsc_.begin(); scit != badsc_.end(); ++ scit) {



    // unpack the SC coordinates from the python file into ix,iy,iz

    iz=int(*scit/1000000);
    iy=*scit%100*iz;
    ix=int((*scit-iy-1000000*iz)/1000)*iz;

    // make the DetId from these coordinates
    scdet=EEDetId(ix,iy,iz);
    det=scdet;

    // loop over the 5x5 SC centered on this DetId and fill discriminating variables

    scan5x5(det,eeRHs,pTopology,pG,nhits,totEt);

    // print some debug info

    if (debug_) {
      std::cout << "EEBadScFilter.cc:  SCID=" <<  *scit << std::endl;
      std::cout << "EEBadScFilter.cc:  ix=" << ix << " iy=" << iy << " iz=" << iz << std::endl;
      std::cout << "EEBadScFilter.cc:  Et(5x5)=" << totEt << " nbadhits=" << nhits << std::endl;
    }


    // if 5x5 transverse energy is above threshold and number of bad hits above threshold
    // event is bad

    if (totEt>EtminSC_ && nhits>=nBadHitsSC_) pass=false;


  }

  // print the decision if event is bad
  if (pass==false && debug_) std::cout << "EEBadScFilter.cc:  REJECT EVENT!!!" << std::endl;


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

  // return the decision

  return taggingMode_ || pass;
}

void EEBadScFilter::scan5x5 ( const DetId &  det, const edm::Handle< EcalRecHitCollection > &  hits, const edm::ESHandle< CaloTopology > &  caloTopo, const edm::ESHandle< CaloGeometry > &  geometry, int &  nHits, float &  totEt  )

privatevirtual

Definition at line 86 of file EEBadScFilter.cc.

References EcalRecHit::checkFlag(), Emin_, EcalRecHit::energy(), CaloNavigator< T, TOPO >::home(), i, j, EcalRecHit::kGood, CaloNavigator< T, TOPO >::offsetBy(), AlCaHLTBitMon_ParallelJobs::p, and side_.

Referenced by filter().

{

  // function to compute:  total transverse energy in a given supercrystal (totEt)
  //                       number of hits with E>Emin_ and rechit flag != kGood (nHits)
  // bad events have large totEt and many high energy hits with rechit flags !kGood

  nHits = 0;
  totEt = 0;

  // navigator to define a 5x5 region around the input DetId

  CaloNavigator<DetId> cursor = CaloNavigator<DetId>(det,caloTopo->getSubdetectorTopology(det));

  // loop over a 5x5 array centered on the input DetId

  for(int j=side_/2; j>=-side_/2; --j)
    {
      for(int i=-side_/2; i<=side_/2; ++i)
    {
      cursor.home();
      cursor.offsetBy(i,j);
      if(hits->find(*cursor)!=hits->end()) // if hit exists in the rechit collection
        {
          EcalRecHit tmpHit = *hits->find(*cursor); // get rechit with detID at cursor


          const GlobalPoint p ( geometry->getPosition(*cursor) ) ; // calculate Et of the rechit
          TVector3 hitPos(p.x(),p.y(),p.z());
          hitPos *= 1.0/hitPos.Mag();
          hitPos *= tmpHit.energy();
          float rechitEt =  hitPos.Pt();

          //--- add rechit E_t to the total for this supercrystal
          totEt += rechitEt;

          // increment nHits if E>Emin and rechit flag is not kGood
          if(tmpHit.energy()>Emin_ && !tmpHit.checkFlag(EcalRecHit::kGood))nHits++;

        }
    }
    }


}

Member Data Documentation

const std::vector<int> EEBadScFilter::badsc_

private

Definition at line 64 of file EEBadScFilter.cc.

Referenced by filter().

const bool EEBadScFilter::debug_

private

Definition at line 67 of file EEBadScFilter.cc.

Referenced by filter().

const edm::EDGetTokenT<EcalRecHitCollection> EEBadScFilter::eeRHSrcToken_

private

Definition at line 57 of file EEBadScFilter.cc.

Referenced by filter().

const double EEBadScFilter::Emin_

private

Definition at line 60 of file EEBadScFilter.cc.

Referenced by scan5x5().

const double EEBadScFilter::EtminSC_

private

Definition at line 61 of file EEBadScFilter.cc.

Referenced by filter().

const int EEBadScFilter::nBadHitsSC_

private

Definition at line 63 of file EEBadScFilter.cc.

Referenced by filter().

const int EEBadScFilter::side_

private

Definition at line 62 of file EEBadScFilter.cc.

Referenced by scan5x5().

const bool EEBadScFilter::taggingMode_

private

Definition at line 66 of file EEBadScFilter.cc.

Referenced by filter().