EEBadScFilter Class Reference

Inheritance diagram for EEBadScFilter:

Public Member Functions
	EEBadScFilter (const edm::ParameterSet &iConfig)
	~EEBadScFilter ()
Public Member Functions inherited from edm::EDFilter
	EDFilter ()
virtual	~EDFilter ()
Public Member Functions inherited from edm::ProducerBase
	ProducerBase ()
void	registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
boost::function< void(const BranchDescription &)>	registrationCallback () const
	used by the fwk to register list of products More...
virtual	~ProducerBase ()

Private Member Functions
virtual bool	filter (edm::Event &iEvent, const edm::EventSetup &iSetup)
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::InputTag	eeRHSrc_
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
typedef WorkerT< EDFilter >	WorkerType
Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList
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::EDFilter
CurrentProcessingContext const *	currentContext () const
Protected Member Functions inherited from edm::ProducerBase
template<class TProducer , class TMethod >
void	callWhenNewProductsRegistered (TProducer *iProd, TMethod iMethod)

Detailed Description

Definition at line 36 of file EEBadScFilter.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 71 of file EEBadScFilter.cc.

  : eeRHSrc_     (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 41 of file EEBadScFilter.cc.

{}

Member Function Documentation

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

privatevirtual

Implements edm::EDFilter.

Definition at line 135 of file EEBadScFilter.cc.

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

                                                                          {


  // load required collections


  // EE rechit collection
  edm::Handle<EcalRecHitCollection> eeRHs;
  iEvent.getByLabel(eeRHSrc_, 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 85 of file EEBadScFilter.cc.

References EcalRecHit::checkFlag(), Emin_, CaloRecHit::energy(), CaloNavigator< T >::home(), i, j, EcalRecHit::kGood, CaloNavigator< T >::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 63 of file EEBadScFilter.cc.

Referenced by filter().

const bool EEBadScFilter::debug_

private

Definition at line 66 of file EEBadScFilter.cc.

Referenced by filter().

const edm::InputTag EEBadScFilter::eeRHSrc_

private

Definition at line 56 of file EEBadScFilter.cc.

Referenced by filter().

const double EEBadScFilter::Emin_

private

Definition at line 59 of file EEBadScFilter.cc.

Referenced by scan5x5().

const double EEBadScFilter::EtminSC_

private

Definition at line 60 of file EEBadScFilter.cc.

Referenced by filter().

const int EEBadScFilter::nBadHitsSC_

private

Definition at line 62 of file EEBadScFilter.cc.

Referenced by filter().

const int EEBadScFilter::side_

private

Definition at line 61 of file EEBadScFilter.cc.

Referenced by scan5x5().

const bool EEBadScFilter::taggingMode_

private

Definition at line 65 of file EEBadScFilter.cc.

Referenced by filter().