#include <UserCode/HcalLaserHBHEFilter2012/src/HcalLaserHBHEFilter2012.cc>

Inheritance diagram for HcalLaserHBHEFilter2012:

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

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

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

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

Private Member Functions
void	endJob () override

bool	filter (edm::Event &, const edm::EventSetup &) override

Private Attributes
std::vector< int >	CalibTS_

edm::InputTag	digiLabel_

bool	forceFilterTrue_

double	HBHEcalibThreshold_

int	minCalibChannelsHBHELaser_

double	minFracDiffHBHELaser_

std::ofstream	outfile_

std::string	prefix_

edm::EDGetTokenT< HcalCalibDigiCollection >	tok_calib_

edm::EDGetTokenT< HBHEDigiCollection >	tok_hbhe_

bool	verbose_

bool	WriteBadToFile_

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

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: [filters out HBHE laser events based on number of HBHE calib digis above threshold and the distribution of HBHE calib digis]

Implementation: [Notes on implementation]

Definition at line 51 of file HcalLaserHBHEFilter2012.cc.

Constructor & Destructor Documentation

HcalLaserHBHEFilter2012::HcalLaserHBHEFilter2012 ( const edm::ParameterSet & ps )

explicit

Definition at line 89 of file HcalLaserHBHEFilter2012.cc.

References CalibTS_, digiLabel_, forceFilterTrue_, edm::ParameterSet::getUntrackedParameter(), HBHEcalibThreshold_, mps_fire::i, minCalibChannelsHBHELaser_, minFracDiffHBHELaser_, outfile_, prefix_, AlCaHLTBitMon_QueryRunRegistry::string, tok_calib_, tok_hbhe_, verbose_, and WriteBadToFile_.

 {
    //now do what ever initialization is needed
   verbose_ = ps.getUntrackedParameter<bool>("verbose",false);
   prefix_  = ps.getUntrackedParameter<std::string>("prefix","");
   minCalibChannelsHBHELaser_ = ps.getUntrackedParameter<int>("minCalibChannelsHBHELaser",20);
   minFracDiffHBHELaser_      = ps.getUntrackedParameter<double>("minFracDiffHBHELaser",0.3);
   edm::InputTag digi_default("hcalDigis");
   digiLabel_                 = ps.getUntrackedParameter<edm::InputTag>("digiLabel",digi_default);
 
   tok_calib_ = consumes<HcalCalibDigiCollection>(digiLabel_);
   tok_hbhe_ = consumes<HBHEDigiCollection>(digiLabel_);
 
   HBHEcalibThreshold_        = ps.getUntrackedParameter<double>("HBHEcalibThreshold",15);
   std::vector <int> dummyTS;
   for (int i=3;i<=6;++i)
     dummyTS.push_back(i);
   CalibTS_                   = ps.getUntrackedParameter<std::vector<int> >("CalibTS",dummyTS);
   WriteBadToFile_=ps.getUntrackedParameter<bool>("WriteBadToFile",false);
   if (WriteBadToFile_)
     outfile_.open("badHcalLaserList_hbhefilter.txt");
   forceFilterTrue_=ps.getUntrackedParameter<bool>("forceFilterTrue",false);
 } // HcalLaserHBHEFilter2012::HcalLaserHBHEFilter2012  constructor

HcalLaserHBHEFilter2012::~HcalLaserHBHEFilter2012 ( )

override

Definition at line 114 of file HcalLaserHBHEFilter2012.cc.

 {
  
    // do anything here that needs to be done at destruction time
    // (e.g. close files, deallocate resources etc.)
 
 }

Member Function Documentation

void HcalLaserHBHEFilter2012::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDFilter.

Definition at line 235 of file HcalLaserHBHEFilter2012.cc.

References outfile_, and WriteBadToFile_.

                                 {
   if (WriteBadToFile_) outfile_.close();
 
 }

void HcalLaserHBHEFilter2012::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 242 of file HcalLaserHBHEFilter2012.cc.

References edm::ConfigurationDescriptions::addDefault(), DEFINE_FWK_MODULE, and edm::ParameterSetDescription::setUnknown().

                                                                                     {
   //The following says we do not know what parameters are allowed so do no validation
   // Please change this to state exactly what you do use, even if it is no parameters
   edm::ParameterSetDescription desc;
   desc.setUnknown();
   descriptions.addDefault(desc);
 }

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

overrideprivate

Definition at line 129 of file HcalLaserHBHEFilter2012.cc.

References edm::SortedCollection< T, SORT >::begin(), HcalCalibDetId::calibFlavor(), CalibTS_, gather_cfg::cout, digiLabel_, edm::SortedCollection< T, SORT >::end(), edm::EventID::event(), forceFilterTrue_, edm::Event::getByToken(), photonIsolationHIProducer_cfi::hbhe, HBHEcalibThreshold_, HcalBarrel, HcalEndcap, HcalCalibDetId::hcalSubdet(), HcalCalibDetId::HOCrosstalk, mps_fire::i, HcalCalibDataFrame::id(), HBHEDataFrame::id(), edm::EventBase::id(), HcalDetId::ieta(), HcalDetId::iphi(), edm::HandleBase::isValid(), edm::EventBase::luminosityBlock(), minCalibChannelsHBHELaser_, minFracDiffHBHELaser_, outfile_, prefix_, edm::EventID::run(), HcalDetId::subdet(), tok_calib_, tok_hbhe_, verbose_, WriteBadToFile_, and HcalCalibDataFrame::zsMarkAndPass().

Referenced by Vispa.Plugins.Browser.BrowserTabController.BrowserTabController::filter(), Vispa.Plugins.Browser.BrowserTabController.BrowserTabController::find(), Vispa.Plugins.Browser.BrowserTabController.BrowserTabController::setDataAccessor(), and Vispa.Plugins.Browser.BrowserTabController.BrowserTabController::switchCenterView().

 {
   // Step 1:: try to get calib digi and HBHE collections.
   // Return true if collection not found?  Or false?  What should default behavior be?
   edm::Handle<HcalCalibDigiCollection> calib_digi;
    if (!(iEvent.getByToken(tok_calib_,calib_digi)))
     {
       edm::LogWarning("HcalLaserHBHEFilter2012")<< digiLabel_<<" calib_digi not available";
       return true;
     }
    
    if (!(calib_digi.isValid()))
      {
        edm::LogWarning("HcalLaserHBHEFilter2012")<< digiLabel_<<" calib_digi is not valid";
        return true;
     }
 
    edm::Handle<HBHEDigiCollection> hbhe_digi;
    if (!(iEvent.getByToken(tok_hbhe_,hbhe_digi)))
     {
       edm::LogWarning("HcalLaserHBHEFilter2012")<< digiLabel_<<" hbhe_digi not available";
       return true;
     }
    
    if (!(hbhe_digi.isValid()))
      {
        edm::LogWarning("HcalLaserHBHEFilter2012")<< digiLabel_<<" hbhe_digi is not valid";
        return true;
     }
 
   // Step 2:  Count HBHE digi calib channels
   int ncalibHBHE=0; // this will track number of HBHE digi channels
 
   
   for (HcalCalibDigiCollection::const_iterator Calibiter = calib_digi->begin();
        Calibiter != calib_digi->end(); ++ Calibiter)
      {
        const HcalCalibDataFrame digi = (const HcalCalibDataFrame)(*Calibiter);
        if (digi.zsMarkAndPass()) continue; // skip digis labeled as "mark and pass" in NZS events
        HcalCalibDetId myid=(HcalCalibDetId)digi.id();
        if (myid.hcalSubdet()!=HcalBarrel && myid.hcalSubdet()!=HcalEndcap) continue;
        if ( myid.calibFlavor()==HcalCalibDetId::HOCrosstalk) continue;
        // Compute charge in current channel (for relevant TS only)
        // If total charge in channel exceeds threshold, increment count of calib channels
        double thischarge=0;
        for (unsigned int i=0;i<CalibTS_.size();++i)
      {
        thischarge+=digi[CalibTS_[i]].nominal_fC();
        if (thischarge>  HBHEcalibThreshold_ )
          {
            ++ncalibHBHE;
            break;
          }
      }
 
        if (ncalibHBHE>=minCalibChannelsHBHELaser_)
      {
        if (verbose_) 
          std::cout <<prefix_<<iEvent.id().run()<<":"<<iEvent.luminosityBlock()<<":"<<iEvent.id().event()<<std::endl;
        if (WriteBadToFile_)
          outfile_<<iEvent.id().run()<<":"<<iEvent.luminosityBlock()<<":"<<iEvent.id().event()<<std::endl;
        if (forceFilterTrue_) return true; // if special input boolean set, always return true, regardless of filter decision
        else return false;
      }
      }
 
   // Step 3:  Look at distribution of HBHE hits
   
   // Count digis in good, bad RBXes.  ('bad' RBXes see no laser signal)
   double badrbxfrac=0.;
   double goodrbxfrac=0.;
   int Nbad=72*3;  // 3 bad RBXes, 72 channels each
   int Ngood=2592*2-Nbad;  // remaining HBHE channels are 'good'
   for (HBHEDigiCollection::const_iterator hbhe = hbhe_digi->begin();
        hbhe != hbhe_digi->end(); ++ hbhe)
     {
       const HBHEDataFrame digi = (const HBHEDataFrame)(*hbhe);
       HcalDetId myid=(HcalDetId)digi.id();
       bool isbad=false; // assume channel is not bad
       if ( myid.subdet()==HcalBarrel && myid.ieta()<0 )
     {
       if (myid.iphi()>=15 && myid.iphi()<=18) isbad=true;
       else if (myid.iphi()>=27 && myid.iphi()<=34) isbad=true;
     }
       if (isbad==true)   badrbxfrac+=1.;
        else goodrbxfrac+=1.;
     }
   goodrbxfrac/=Ngood;
   badrbxfrac/=Nbad;
   if (goodrbxfrac-badrbxfrac>minFracDiffHBHELaser_)
     {
       if (verbose_) 
     std::cout <<prefix_<<iEvent.id().run()<<":"<<iEvent.luminosityBlock()<<":"<<iEvent.id().event()<<std::endl;
       if (WriteBadToFile_)
          outfile_<<iEvent.id().run()<<":"<<iEvent.luminosityBlock()<<":"<<iEvent.id().event()<<std::endl;
       
       if (forceFilterTrue_) return true; // if special input boolean set, always return true, regardless of filter decision
       else return false;
     }   
   // Step 4:  HBHE laser tests passed.  return true
   return true;
 
 }  // HcalLaserHBHEFilter2012::filter