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

Inheritance diagram for HcalLaserHBHEHFFilter2012:

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

	~HcalLaserHBHEHFFilter2012 ()

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

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

Private Member Functions
virtual void	endJob () override

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

Private Attributes
std::vector< int >	CalibTS_

edm::InputTag	digiLabel_

bool	filterHBHE_

bool	filterHF_

bool	forceFilterTrue_

double	HBHEcalibThreshold_

int	minCalibChannelsHBHELaser_

int	minCalibChannelsHFLaser_

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
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 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 HcalLaserHBHEHFFilter2012.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 95 of file HcalLaserHBHEHFFilter2012.cc.

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

 {
    //now do what ever initialization is needed
   filterHBHE_                = ps.getParameter<bool>("filterHBHE");
   minCalibChannelsHBHELaser_ = ps.getParameter<int>("minCalibChannelsHBHELaser");
   minFracDiffHBHELaser_      = ps.getParameter<double>("minFracDiffHBHELaser");
   HBHEcalibThreshold_        = ps.getParameter<double>("HBHEcalibThreshold");
   CalibTS_                   = ps.getParameter<std::vector<int> >("CalibTS");
 
   filterHF_                  = ps.getParameter<bool>("filterHF");
   minCalibChannelsHFLaser_   = ps.getParameter<int>("minCalibChannelsHFLaser");
 
   digiLabel_                 = ps.getParameter<edm::InputTag>("digiLabel");
   tok_calib_ = consumes<HcalCalibDigiCollection>(digiLabel_);
   tok_hbhe_ = consumes<HBHEDigiCollection>(digiLabel_);
 
   verbose_                   = ps.getUntrackedParameter<bool>("verbose",false);
   prefix_                    = ps.getUntrackedParameter<std::string>("prefix","");
   WriteBadToFile_            = ps.getUntrackedParameter<bool>("WriteBadToFile",false);
   forceFilterTrue_           = ps.getUntrackedParameter<bool>("forceFilterTrue",false);
   if (WriteBadToFile_)
     outfile_.open("badHcalLaserList_hcalfilter.txt");
 } // HcalLaserHBHEHFFilter2012::HcalLaserHBHEHFFilter2012  constructor

HcalLaserHBHEHFFilter2012::~HcalLaserHBHEHFFilter2012 ( )

Definition at line 120 of file HcalLaserHBHEHFFilter2012.cc.

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

Member Function Documentation

void HcalLaserHBHEHFFilter2012::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDFilter.

Definition at line 253 of file HcalLaserHBHEHFFilter2012.cc.

References outfile_, and WriteBadToFile_.

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

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

static

Definition at line 260 of file HcalLaserHBHEHFFilter2012.cc.

References edm::ConfigurationDescriptions::addDefault(), 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 HcalLaserHBHEHFFilter2012::filter	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

overrideprivatevirtual

Implements edm::EDFilter.

Definition at line 135 of file HcalLaserHBHEHFFilter2012.cc.

References HcalCalibDetId::calibFlavor(), CalibTS_, gather_cfg::cout, digiLabel_, edm::EventID::event(), filterHBHE_, filterHF_, forceFilterTrue_, edm::Event::getByToken(), HBHEcalibThreshold_, HcalBarrel, HcalEndcap, HcalForward, HcalCalibDetId::hcalSubdet(), HcalCalibDetId::HOCrosstalk, i, HBHEDataFrame::id(), HcalCalibDataFrame::id(), edm::EventBase::id(), HcalDetId::ieta(), HcalDetId::iphi(), edm::HandleBase::isValid(), edm::EventBase::luminosityBlock(), minCalibChannelsHBHELaser_, minCalibChannelsHFLaser_, 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("HcalLaserFilter2012")<< digiLabel_<<" calib_digi not available";
       return true;
     }
   
   if (!(calib_digi.isValid()))
     {
       edm::LogWarning("HcalLaserFilter2012")<< digiLabel_<<" calib_digi is not valid";
       return true;
     }
 
   edm::Handle<HBHEDigiCollection> hbhe_digi;
   if (!(iEvent.getByToken(tok_hbhe_,hbhe_digi)))
     {
       edm::LogWarning("HcalLaserFilter2012")<< digiLabel_<<" hbhe_digi not available";
       return true;
     }
    
   if (!(hbhe_digi.isValid()))
     {
       edm::LogWarning("HcalLaserHBHEHFFilter2012")<< 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
   int ncalibHF=0;   // this will track number of HF 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 (filterHBHE_ && (myid.hcalSubdet()==HcalBarrel || myid.hcalSubdet()==HcalEndcap)) 
       {
     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;
     } // if  (ncalibHBHE>=minCalibChannelsHBHELaser_)
 
       }
     else if (filterHF_ && (myid.hcalSubdet()==HcalForward)) 
       {
     ++ncalibHF;
     if (ncalibHF>=minCalibChannelsHFLaser_) {
       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;
     }
       }
   }
 
   if (filterHBHE_) {
     // 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:  HBHEHF laser tests passed.  return true
   return true;
 
 }  // HcalLaserHBHEHFFilter2012::filter