#include <EventSelector.h>

Classes
struct	BitInfo

Public Types
typedef std::vector< std::string >	Strings

Public Member Functions
bool	acceptEvent (TriggerResults const &)

bool	acceptEvent (unsigned char const *, int) const

	EventSelector (Strings const &pathspecs, Strings const &names)

	EventSelector (Strings const &pathspecs)

	EventSelector (edm::ParameterSet const &pset, Strings const &triggernames)

boost::shared_ptr< TriggerResults >	maskTriggerResults (TriggerResults const &inputResults)

bool	wantAll () const

Static Public Member Functions
static void	fillDescription (ParameterSetDescription &desc)

static std::vector< std::string >	getEventSelectionVString (edm::ParameterSet const &pset)

static bool	selectionIsValid (Strings const &pathspec, Strings const &fullTriggerList)

static evtSel::OverlapResult	testSelectionOverlap (Strings const &pathspec1, Strings const &pathspec2, Strings const &fullTriggerList)

Private Types
typedef std::vector< BitInfo >	Bits

Private Member Functions
bool	acceptAllBits (Bits const &b, HLTGlobalStatus const &tr) const

bool	acceptOneBit (Bits const &b, HLTGlobalStatus const &tr, hlt::HLTState const &s=hlt::Ready) const

bool	acceptTriggerPath (HLTPathStatus const &, BitInfo const &) const

bool	containsExceptions (HLTGlobalStatus const &tr) const

void	init (Strings const &paths, Strings const &triggernames)

bool	selectionDecision (HLTGlobalStatus const &tr) const

Static Private Member Functions
static std::vector< bool >	combine (std::vector< bool > const &a, std::vector< bool > const &b)

static std::vector< bool >	expandDecisionList (Bits const &b, bool PassOrFail, unsigned int n)

static std::string	glob2reg (std::string const &s)

static bool	identical (std::vector< bool > const &a, std::vector< bool > const &b)

static bool	identical (EventSelector const &a, EventSelector const &b, unsigned int N)

static std::vector < Strings::const_iterator >	matching_triggers (Strings const &trigs, std::string const &s)

static bool	overlapping (std::vector< bool > const &a, std::vector< bool > const &b)

static bool	subset (std::vector< bool > const &a, std::vector< bool > const &b)

Private Attributes
Bits	absolute_acceptors_

bool	accept_all_

std::vector< Bits >	all_must_fail_

std::vector< Bits >	all_must_fail_noex_

Bits	conditional_acceptors_

Bits	exception_acceptors_

bool	notStarPresent_

int	nTriggerNames_

Strings	paths_

ParameterSetID	psetID_

bool	psetID_initialized_

bool	results_from_current_process_

Detailed Description

Definition at line 40 of file EventSelector.h.

Member Typedef Documentation

typedef std::vector<BitInfo> edm::EventSelector::Bits

private

Definition at line 88 of file EventSelector.h.

typedef std::vector<std::string> edm::EventSelector::Strings

Definition at line 44 of file EventSelector.h.

Constructor & Destructor Documentation

EventSelector::EventSelector	(	Strings const &	pathspecs,
		Strings const &	names
	)

Definition at line 49 of file EventSelector.cc.

References init().

                                         :
     accept_all_(false),
     absolute_acceptors_(),
     conditional_acceptors_(),
     exception_acceptors_(),
     all_must_fail_(),
     all_must_fail_noex_(),
     results_from_current_process_(true),
     psetID_initialized_(false),
     psetID_(),
     paths_(),
     nTriggerNames_(0),
     notStarPresent_(false)
   {
     init(pathspecs, names);
   }

EventSelector::EventSelector ( Strings const & pathspecs )

explicit

Definition at line 67 of file EventSelector.cc.

                                                       :
     accept_all_(false),
     absolute_acceptors_(),
     conditional_acceptors_(),
     exception_acceptors_(),
     all_must_fail_(),
     all_must_fail_noex_(),
     results_from_current_process_(false),
     psetID_initialized_(false),
     psetID_(),
     paths_(pathspecs),
     nTriggerNames_(0),
     notStarPresent_(false)
   {
   }

EventSelector::EventSelector	(	edm::ParameterSet const &	pset,
		Strings const &	triggernames
	)

Definition at line 83 of file EventSelector.cc.

References edm::ParameterSet::empty(), edm::ParameterSet::getParameter(), init(), and mergeAndRegister::paths.

                                                :
     accept_all_(false),
     absolute_acceptors_(),
     conditional_acceptors_(),
     exception_acceptors_(),
     all_must_fail_(),
     all_must_fail_noex_(),
     results_from_current_process_(true),
     psetID_initialized_(false),
     psetID_(),
     paths_(),
     nTriggerNames_(0),
     notStarPresent_(false)
   {
     Strings paths; // default is empty...
 
     if (!config.empty())
       paths = config.getParameter<Strings>("SelectEvents");
 
     init(paths, triggernames);
   }

Member Function Documentation

bool EventSelector::acceptAllBits	(	Bits const &	b,
		HLTGlobalStatus const &	tr
	)		const

private

Definition at line 433 of file EventSelector.cc.

References alignCSCRings::e, edm::hlt::Fail, i, edm::hlt::Pass, and edm::HLTGlobalStatus::state().

Referenced by selectionDecision().

   {
     Bits::const_iterator i(b.begin());
     Bits::const_iterator e(b.end());
     for(;i!=e;++i) {
       hlt::HLTState bstate = i->accept_state_ ? hlt::Pass : hlt::Fail;
       if (tr[i->pos_].state() != bstate) return false;
     }
     return true;    
   } // acceptAllBits                   

bool EventSelector::acceptEvent ( TriggerResults const & tr )

Definition at line 270 of file EventSelector.cc.

References accept_all_, edm::hlt::Exception, init(), edm::TriggerResults::parameterSetID(), paths_, psetID_, psetID_initialized_, results_from_current_process_, selectionDecision(), and edm::errors::Unknown.

Referenced by edm::detail::NamedEventSelector::match(), and selectionIsValid().

   {
     if (accept_all_) return true;
     
     // For the current process we already initialized in the constructor,
     // The trigger names will not change so we can skip initialization.
     if (!results_from_current_process_) {
   
       // For previous processes we need to get the trigger names that
       // correspond to the bits in TriggerResults from the ParameterSet
       // set registry, which is stored once per file.  The ParameterSetID
       // stored in TriggerResults is the key used to find the info in the
       // registry.  We optimize using the fact the ID is unique. If the ID
       // has not changed since the last time we initialized with new triggernames,
       // then the names have not changed and we can skip this initialization.
       if (!(psetID_initialized_ && psetID_ == tr.parameterSetID())) {
 
         Strings triggernames;
         bool fromPSetRegistry;
 
         Service<service::TriggerNamesService> tns;
         if (tns->getTrigPaths(tr, triggernames, fromPSetRegistry)) {
 
           init(paths_, triggernames);
 
           if (fromPSetRegistry) {
             psetID_ = tr.parameterSetID();
             psetID_initialized_ = true;
           }
           else {
             psetID_initialized_ = false;
           }
         }
         // This should never happen
         else {
           throw edm::Exception(errors::Unknown)
             << "EventSelector::acceptEvent cannot find the trigger names for\n"
                "a process for which the configuration has requested that the\n"
                "OutputModule use TriggerResults to select events from.  This should\n"
                "be impossible, please send information to reproduce this problem to\n"
                "the edm developers.\n"; 
     }
       }
     }
 
     // Now make the decision, based on the supplied TriggerResults tr,
     // which of course can be treated as an HLTGlobalStatus by inheritance
     
     return selectionDecision(tr);
     
   } // acceptEvent(TriggerResults const& tr)

bool EventSelector::acceptEvent	(	unsigned char const *	array_of_trigger_results,
		int	number_of_trigger_paths
	)		const

Definition at line 323 of file EventSelector.cc.

References accept_all_, edm::errors::Configuration, edm::hlt::Exception, results_from_current_process_, and selectionDecision().

   {
 
     // This should never occur unless someone uses this function in
     // an incorrect way ...
     if (!results_from_current_process_) {
       throw edm::Exception(errors::Configuration)
         << "\nEventSelector.cc::acceptEvent, you are attempting to\n"
         << "use a bit array for trigger results instead of the\n"
         << "TriggerResults object for a previous process.  This\n"
         << "will not work and ought to be impossible\n";
     }
 
     if (accept_all_) return true;
 
     // Form HLTGlobalStatus object to represent the array_of_trigger_results
     HLTGlobalStatus tr(number_of_trigger_paths);
     int byteIndex = 0;
     int subIndex  = 0;
     for (int pathIndex = 0; pathIndex < number_of_trigger_paths; ++pathIndex)
     {
       int state = array_of_trigger_results[byteIndex] >> (subIndex * 2);
       state &= 0x3;
       HLTPathStatus pathStatus(static_cast<hlt::HLTState>(state));
       tr[pathIndex] = pathStatus;
       ++subIndex;
       if (subIndex == 4)
       { ++byteIndex;
         subIndex = 0;
       }
     }    
 
     // Now make the decision, based on the HLTGlobalStatus tr,
     // which we have created from the supplied array of results
     
     return selectionDecision(tr);
 
   } // acceptEvent(array_of_trigger_results, number_of_trigger_paths)

bool EventSelector::acceptOneBit	(	Bits const &	b,
		HLTGlobalStatus const &	tr,
		hlt::HLTState const &	s = `hlt::Ready`
	)		const

private

Definition at line 413 of file EventSelector.cc.

References alignCSCRings::e, edm::hlt::Exception, edm::hlt::Fail, i, edm::hlt::Pass, and edm::HLTGlobalStatus::state().

Referenced by selectionDecision().

   {
     bool lookForException = (s == hlt::Exception);
     Bits::const_iterator i(b.begin());
     Bits::const_iterator e(b.end());
     for(;i!=e;++i) {
       hlt::HLTState bstate = 
           lookForException ? hlt::Exception
                    : i->accept_state_ ? hlt::Pass
                               : hlt::Fail;
       if (tr[i->pos_].state() == bstate) return true;
     }
     return false;    
   } // acceptOneBit                

bool EventSelector::acceptTriggerPath	(	HLTPathStatus const &	pathStatus,
		BitInfo const &	pathInfo
	)		const

private

Definition at line 401 of file EventSelector.cc.

References edm::EventSelector::BitInfo::accept_state_, edm::hlt::Exception, edm::hlt::Fail, edm::hlt::Pass, and edm::HLTPathStatus::state().

   {
     return (((pathStatus.state()==hlt::Pass) && (pathInfo.accept_state_)) ||
             ((pathStatus.state()==hlt::Fail) && !(pathInfo.accept_state_)) ||
             ((pathStatus.state()==hlt::Exception)));
   }

std::vector< bool > EventSelector::combine	(	std::vector< bool > const &	a,
		std::vector< bool > const &	b
	)

staticprivate

Definition at line 1116 of file EventSelector.cc.

References i, and x.

Referenced by testSelectionOverlap().

   {
     assert(a.size() == b.size());
     std::vector<bool> x(a.size());
     for (unsigned int i = 0; i != a.size(); ++i) {
       x[i] = a[i] || b[i];
     } // a really sharp compiler will optimize the hell out of this, 
       // exploiting word-size OR operations.
     return x;
   } // combine                            

bool EventSelector::containsExceptions ( HLTGlobalStatus const & tr ) const

private

Definition at line 963 of file EventSelector.cc.

References alignCSCRings::e, edm::hlt::Exception, i, and edm::HLTGlobalStatus::size().

Referenced by selectionDecision().

   {
     unsigned int e = tr.size();
     for (unsigned int i = 0; i < e; ++i) {
       if (tr[i].state() == hlt::Exception) return true;
     }
     return false;
   }

std::vector< bool > EventSelector::expandDecisionList	(	Bits const &	b,
		bool	PassOrFail,
		unsigned int	n
	)

staticprivate

Definition at line 1053 of file EventSelector.cc.

References i, and x.

Referenced by identical(), maskTriggerResults(), and testSelectionOverlap().

   {
     std::vector<bool> x(n, false);
     for (unsigned int i = 0; i != b.size(); ++i) {
       if (b[i].accept_state_ == PassOrFail) x[b[i].pos_] = true;
     }
     return x;
   } // expandDecisionList   

void EventSelector::fillDescription ( ParameterSetDescription & desc )

static

Definition at line 1129 of file EventSelector.cc.

References edm::ParameterSetDescription::addOptional(), and edm::ParameterSetDescription::addUntracked().

Referenced by edm::OutputModule::fillDescription(), and edm::one::OutputModuleBase::fillDescription().

                                                               {
     ParameterSetDescription selector;
     selector.addOptional<std::vector<std::string> >("SelectEvents");
     desc.addUntracked<ParameterSetDescription>("SelectEvents", selector);
   }

std::vector< std::string > EventSelector::getEventSelectionVString ( edm::ParameterSet const & pset )

static

Returns the list of strings that correspond to the trigger selection request in the specified parameter set (the list of strings contained in the "SelectEvents" parameter).

Parameters

pset	The ParameterSet that contains the trigger selection.

Returns: the trigger selection list (vector of string).

Definition at line 939 of file EventSelector.cc.

References edm::ParameterSet::empty(), edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), and corrVsCorr::selection.

Referenced by dqmservices::TriggerSelector::getEventSelectionVString(), and edm::StreamerOutputModuleBase::StreamerOutputModuleBase().

   {
     // default the selection to everything (wildcard)
     Strings selection;
     selection.push_back("*");
     selection.push_back("!*");
     selection.push_back("exception@*");
 
     // the SelectEvents parameter is a ParameterSet within
     // a ParameterSet, so we have to pull it out twice
     ParameterSet selectEventsParamSet =
       pset.getUntrackedParameter("SelectEvents", ParameterSet());
     if (!selectEventsParamSet.empty()) {
       Strings path_specs = 
         selectEventsParamSet.getParameter<Strings>("SelectEvents");
       if (!path_specs.empty()) {
         selection = path_specs;
       }
     }
 
     // return the result
     return selection;
   }

static std::string edm::EventSelector::glob2reg ( std::string const & s )

staticprivate

bool EventSelector::identical	(	std::vector< bool > const &	a,
		std::vector< bool > const &	b
	)

staticprivate

Definition at line 975 of file EventSelector.cc.

References i, and n.

Referenced by identical(), and testSelectionOverlap().

                                          {
      unsigned int n = a.size();
      if (n != b.size()) return false;
      for (unsigned int i=0; i!=n; ++i) {
        if (a[i] != b[i]) return false;
      }
      return true;
   }

bool EventSelector::identical	(	EventSelector const &	a,
		EventSelector const &	b,
		unsigned int	N
	)

staticprivate

Definition at line 986 of file EventSelector.cc.

References absolute_acceptors_, all_must_fail_, all_must_fail_noex_, conditional_acceptors_, exception_acceptors_, expandDecisionList(), identical(), gen::k, m, and edm::match().

   {
         // create the expanded masks for the various decision lists in a and b
     if (!identical(expandDecisionList(a.absolute_acceptors_,true,N),
                    expandDecisionList(b.absolute_acceptors_,true,N))) 
            return false;
     if (!identical(expandDecisionList(a.conditional_acceptors_,true,N),
                    expandDecisionList(b.conditional_acceptors_,true,N))) 
            return false;
     if (!identical(expandDecisionList(a.absolute_acceptors_,false,N),
                    expandDecisionList(b.absolute_acceptors_,false,N)))
            return false;
     if (!identical(expandDecisionList(a.conditional_acceptors_,false,N),
                    expandDecisionList(b.conditional_acceptors_,false,N))) 
            return false;
     if (!identical(expandDecisionList(a.exception_acceptors_,true,N),
                    expandDecisionList(b.exception_acceptors_,true,N)))
            return false;
     if (a.all_must_fail_.size() != b.all_must_fail_.size()) return false;
     
     std::vector< std::vector<bool> > aMustFail;
     for (unsigned int m = 0; m != a.all_must_fail_.size(); ++m) {
       aMustFail.push_back(expandDecisionList(a.all_must_fail_[m],false,N));
     }
     std::vector< std::vector<bool> > aMustFailNoex;
     for (unsigned int m = 0; m != a.all_must_fail_noex_.size(); ++m) {
       aMustFailNoex.push_back 
           (expandDecisionList(a.all_must_fail_noex_[m],false,N));
     }
     std::vector< std::vector<bool> > bMustFail;
     for (unsigned int m = 0; m != b.all_must_fail_.size(); ++m) {
       bMustFail.push_back(expandDecisionList(b.all_must_fail_[m],false,N));
     }
     std::vector< std::vector<bool> > bMustFailNoex;
     for (unsigned int m = 0; m != b.all_must_fail_noex_.size(); ++m) {
       bMustFailNoex.push_back 
           (expandDecisionList(b.all_must_fail_noex_[m],false,N));
     }
     
     for (unsigned int m = 0; m != aMustFail.size(); ++m) {
       bool match = false;
       for (unsigned int k = 0; k != bMustFail.size(); ++k) {
         if (identical(aMustFail[m],bMustFail[k])) {
           match = true;
       break;
     }
       }
       if (!match) return false;
     }
     for (unsigned int m = 0; m != aMustFailNoex.size(); ++m) {
       bool match = false;
       for (unsigned int k = 0; k != bMustFailNoex.size(); ++k) {
          if (identical(aMustFailNoex[m],bMustFailNoex[k])) {
           match = true;
       break;
     }
       }
       if (!match) return false;
     }
 
     return true;
     
   } // identical (EventSelector, EventSelector, N);

void EventSelector::init	(	Strings const &	paths,
		Strings const &	triggernames
	)

private

Definition at line 107 of file EventSelector.cc.

References absolute_acceptors_, accept_all_, all_must_fail_, all_must_fail_noex_, conditional_acceptors_, edm::errors::Configuration, end, edm::hlt::Exception, exception_acceptors_, i, edm::is_glob(), notStarPresent_, nTriggerNames_, edm::regexMatch(), AlCaHLTBitMon_QueryRunRegistry::string, and edmStreamStallGrapher::t.

Referenced by acceptEvent(), and EventSelector().

   {
     // std::cerr << "### init entered\n";
     accept_all_ = false;
     absolute_acceptors_.clear(),
     conditional_acceptors_.clear(),
     exception_acceptors_.clear(),
     all_must_fail_.clear();
     all_must_fail_noex_.clear();
     nTriggerNames_ = triggernames.size();
     notStarPresent_ = false;
 
     if (paths.empty())
       {
     accept_all_ = true;
     return;
       }
 
     // The following are for the purpose of establishing accept_all_ by 
     // virtue of an inclusive set of paths:
     bool unrestricted_star = false;
     bool negated_star      = false;
     bool exception_star    = false;
     
     for (Strings::const_iterator i(paths.begin()), end(paths.end()); 
      i!=end; ++i)
     {
       std::string pathSpecifier(*i);
       boost::erase_all(pathSpecifier, " \t"); // whitespace eliminated
       if (pathSpecifier == "*")           unrestricted_star = true;
       if (pathSpecifier == "!*")          negated_star = true;
       if (pathSpecifier == "exception@*") exception_star = true;
 
       std::string basePathSpec(pathSpecifier);
       bool noex_demanded = false;
       std::string::size_type 
           and_noexception = pathSpecifier.find("&noexception");
       if (and_noexception != std::string::npos) {
     basePathSpec = pathSpecifier.substr(0,and_noexception);
         noex_demanded = true;
       }
       std::string::size_type and_noex = pathSpecifier.find("&noex");
       if (and_noex != std::string::npos) {
     basePathSpec = pathSpecifier.substr(0,and_noexception);
         noex_demanded = true;
       }
       and_noexception = basePathSpec.find("&noexception");
       and_noex = basePathSpec.find("&noex");    
       if (and_noexception != std::string::npos ||
        and_noex != std::string::npos)
           throw edm::Exception(errors::Configuration)
             << "EventSelector::init, An OutputModule is using SelectEvents\n"
                "to request a trigger name, but specifying &noexceptions twice\n"
             << "The improper trigger name is: " << pathSpecifier << "\n";  
 
       std::string realname(basePathSpec);
       bool negative_criterion = false;
       if (basePathSpec[0] == '!') {
     negative_criterion = true;
     realname = basePathSpec.substr(1,std::string::npos);
       }
       bool exception_spec = false;
       if (realname.find("exception@") == 0) {
     exception_spec = true;  
     realname = realname.substr(10, std::string::npos); 
     // strip off 10 chars, which is length of "exception@" 
       } 
       if (negative_criterion &&  exception_spec)
           throw edm::Exception(errors::Configuration)
             << "EventSelector::init, An OutputModule is using SelectEvents\n"
                "to request a trigger name starting with !exception@.\n"
            "This is not supported.\n"
             << "The improper trigger name is: " << pathSpecifier << "\n";  
       if (noex_demanded &&  exception_spec)
           throw edm::Exception(errors::Configuration)
             << "EventSelector::init, An OutputModule is using SelectEvents\n"
                "to request a trigger name starting with exception@ "
            "and also demanding no &exceptions.\n"
             << "The improper trigger name is: " << pathSpecifier << "\n";  
 
 
       // instead of "see if the name can be found in the full list of paths"
       // we want to find all paths that match this name.    
       std::vector<Strings::const_iterator> matches =
           regexMatch(triggernames, realname);
 
       if (matches.empty() && !is_glob(realname)) 
       {
           throw edm::Exception(errors::Configuration)
             << "EventSelector::init, An OutputModule is using SelectEvents\n"
                "to request a trigger name that does not exist\n"
             << "The unknown trigger name is: " << realname << "\n";  
       }
       if (matches.empty() && is_glob(realname)) 
       {
           LogWarning("Configuration")
             << "EventSelector::init, An OutputModule is using SelectEvents\n"
                "to request a wildcarded trigger name that does not match any trigger \n"
             << "The wildcarded trigger name is: " << realname << "\n";  
       }
 
       if (!negative_criterion && !noex_demanded && !exception_spec) {
     for (unsigned int t = 0; t != matches.size(); ++t) {
       BitInfo bi(distance(triggernames.begin(),matches[t]), true);
       absolute_acceptors_.push_back(bi);
     }
       } else if (!negative_criterion && noex_demanded) {
     for (unsigned int t = 0; t != matches.size(); ++t) {
       BitInfo bi(distance(triggernames.begin(),matches[t]), true);
       conditional_acceptors_.push_back(bi);
     }
       } else if (exception_spec) {
     for (unsigned int t = 0; t != matches.size(); ++t) {
       BitInfo bi(distance(triggernames.begin(),matches[t]), true);
       exception_acceptors_.push_back(bi);
     }
       } else if (negative_criterion && !noex_demanded) {
     if (matches.empty()) {
             throw edm::Exception(errors::Configuration)
             << "EventSelector::init, An OutputModule is using SelectEvents\n"
                "to request all fails on a set of trigger names that do not exist\n"
             << "The problematic name is: " << pathSpecifier << "\n";  
 
     } else if (matches.size() == 1) {
       BitInfo bi(distance(triggernames.begin(),matches[0]), false);
       absolute_acceptors_.push_back(bi);
     } else {
       Bits mustfail;
       for (unsigned int t = 0; t != matches.size(); ++t) {
         BitInfo bi(distance(triggernames.begin(),matches[t]), false);
         // We set this to false because that will demand bits are Fail. 
         mustfail.push_back(bi);
       }
       all_must_fail_.push_back(mustfail);
     }   
       } else if (negative_criterion && noex_demanded) {
     if (matches.empty()) {
             throw edm::Exception(errors::Configuration)
             << "EventSelector::init, An OutputModule is using SelectEvents\n"
                "to request all fails on a set of trigger names that do not exist\n"
             << "The problematic name is: " << pathSpecifier << "\n";  
 
     } else if (matches.size() == 1) {
       BitInfo bi(distance(triggernames.begin(),matches[0]), false);
       conditional_acceptors_.push_back(bi);
     } else {
       Bits mustfail;
       for (unsigned int t = 0; t != matches.size(); ++t) {
         BitInfo bi(distance(triggernames.begin(),matches[t]), false);
         mustfail.push_back(bi);
       }
       all_must_fail_noex_.push_back(mustfail);
     }
       } 
     } // end of the for loop on i(paths.begin()), end(paths.end())
 
     if (unrestricted_star && negated_star && exception_star) accept_all_ = true;
 
     // std::cerr << "### init exited\n";
 
   } // EventSelector::init

boost::shared_ptr< TriggerResults > EventSelector::maskTriggerResults ( TriggerResults const & inputResults )

Applies a trigger selection mask to a specified trigger result object. Within the trigger result object, each path status is left unchanged if it satisfies the trigger selection (path specs) or cleared if it does not satisfy the trigger selection. In this way, the resulting trigger result object contains only path status values that "pass" the selection criteria.

Parameters

inputResults The raw trigger results object that will be masked.

Returns: a copy of the input trigger results object with only the path status results that match the trigger selection.

Exceptions

edm::Exception if the number of paths in the TriggerResults object does not match the specified full trigger list, or if the trigger selection is invalid in the context of the full trigger list.

Definition at line 772 of file EventSelector.cc.

References absolute_acceptors_, all_must_fail_, all_must_fail_noex_, conditional_acceptors_, edm::errors::EventCorruption, edm::hlt::Exception, exception_acceptors_, expandDecisionList(), f, edm::hlt::Fail, m, N, nTriggerNames_, edm::TriggerResults::parameterSetID(), edm::hlt::Pass, alignCSCRings::s, edm::HLTGlobalStatus::size(), edm::HLTGlobalStatus::state(), and old-fu_pass::TriggerResults.

   {
     // fetch and validate the total number of paths
     unsigned int fullTriggerCount = nTriggerNames_;
     unsigned int N = fullTriggerCount;
     if (fullTriggerCount != inputResults.size())
     {
       throw edm::Exception(errors::EventCorruption)
         << "EventSelector::maskTriggerResults, the TriggerResults\n"
         << "size (" << inputResults.size()
         << ") does not match the number of paths in the\n"
         << "full trigger list (" << fullTriggerCount << ").\n";
     }
 
     // create a suitable global status object to work with, all in Ready state
     HLTGlobalStatus mask(fullTriggerCount);
     
     // Deal with must_fail acceptors that would cause selection
     for (unsigned int m = 0; m < this->all_must_fail_.size(); ++m) {
       std::vector<bool>  
         f = expandDecisionList(this->all_must_fail_[m],false,N);
       bool all_fail = true;
       for (unsigned int ipath = 0; ipath < N; ++ipath) {        
     if  ((f[ipath]) && (inputResults [ipath].state() != hlt::Fail)) { 
       all_fail = false;
       break;
     }
       }
       if (all_fail) {
     for (unsigned int ipath = 0; ipath < N; ++ipath) {
           if  (f[ipath]) { 
         mask[ipath] = hlt::Fail;
       }
     }
       }
     }
     for (unsigned int m = 0; m < this->all_must_fail_noex_.size(); ++m) {
       std::vector<bool>  
         f = expandDecisionList(this->all_must_fail_noex_[m],false,N);
       bool all_fail = true;
       for (unsigned int ipath = 0; ipath < N; ++ipath) {        
     if ((f[ipath]) && (inputResults [ipath].state() != hlt::Fail)) { 
       all_fail = false;
       break;
     }
       }
       if (all_fail) {
     for (unsigned int ipath = 0; ipath < N; ++ipath) {
           if  (f[ipath]) { 
         mask[ipath] = hlt::Fail;
       }
     }
       }
     } // factoring opportunity - work done for fail_noex_ is same as for fail_
     
     // Deal with normal acceptors that would cause selection
     std::vector<bool> 
       aPassAbs = expandDecisionList(this->absolute_acceptors_,true,N);
     std::vector<bool> 
       aPassCon = expandDecisionList(this->conditional_acceptors_,true,N);
     std::vector<bool> 
       aFailAbs = expandDecisionList(this->absolute_acceptors_,false,N);
     std::vector<bool> 
       aFailCon = expandDecisionList(this->conditional_acceptors_,false,N);
     std::vector<bool> 
       aExc = expandDecisionList(this->exception_acceptors_,true,N);
     for (unsigned int ipath = 0; ipath < N; ++ipath) {
       hlt::HLTState s = inputResults [ipath].state();  
       if (((aPassAbs[ipath]) && (s == hlt::Pass))
             ||
       ((aPassCon[ipath]) && (s == hlt::Pass))       
             ||
       ((aFailAbs[ipath]) && (s == hlt::Fail))       
             ||
       ((aFailCon[ipath]) && (s == hlt::Fail))
             ||
       ((aExc[ipath]) && (s == hlt::Exception)))
       {
         mask[ipath] = s;
       }     
     }
  
     // Based on the global status for the mask, create and return a 
     // TriggerResults
     boost::shared_ptr<TriggerResults>
       maskedResults(new TriggerResults(mask, inputResults.parameterSetID()));
     return maskedResults;
   }  // maskTriggerResults

static std::vector< Strings::const_iterator > edm::EventSelector::matching_triggers	(	Strings const &	trigs,
		std::string const &	s
	)

staticprivate

bool EventSelector::overlapping	(	std::vector< bool > const &	a,
		std::vector< bool > const &	b
	)

staticprivate

Definition at line 1065 of file EventSelector.cc.

References i.

Referenced by testSelectionOverlap().

   {
     if (a.size() != b.size()) return false;
     for (unsigned int i = 0; i != a.size(); ++i) {
       if (a[i] && b[i]) return true;
     }
     return false;
   } // overlapping

bool EventSelector::selectionDecision ( HLTGlobalStatus const & tr ) const

private

Definition at line 364 of file EventSelector.cc.

References absolute_acceptors_, accept_all_, acceptAllBits(), acceptOneBit(), all_must_fail_, all_must_fail_noex_, conditional_acceptors_, containsExceptions(), edm::hlt::Exception, exception_acceptors_, and f.

Referenced by acceptEvent().

   {
     if (accept_all_) return true;
 
     bool exceptionPresent = false;
     bool exceptionsLookedFor = false;
     
     if (acceptOneBit(absolute_acceptors_, tr)) return true;
     if (acceptOneBit(conditional_acceptors_, tr)) {
       exceptionPresent = containsExceptions(tr);
       if (!exceptionPresent) return true;
       exceptionsLookedFor = true;
     }
     if (acceptOneBit(exception_acceptors_, tr, hlt::Exception)) return true;
 
     for (std::vector<Bits>::const_iterator f =  all_must_fail_.begin();
                            f != all_must_fail_.end(); ++f)
     {
       if (acceptAllBits(*f, tr)) return true;
     }
     for (std::vector<Bits>::const_iterator fn =  all_must_fail_noex_.begin();
                            fn != all_must_fail_noex_.end(); ++fn)
     {
       if (acceptAllBits(*fn, tr)) {
         if (!exceptionsLookedFor) exceptionPresent = containsExceptions(tr);
         return (!exceptionPresent);
       }
     }
     
     // If we have not accepted based on any of the acceptors, nor on any one of
     // the all_must_fail_ collections, then we reject this event.
     
     return false;
   
   }  // selectionDecision()

bool EventSelector::selectionIsValid	(	Strings const &	pathspec,
		Strings const &	fullTriggerList
	)

static

Tests if the specified trigger selection list (path spec) is valid in the context of the specified full trigger list. Each element in the selection list is tested to see if it possible for some combination of trigger results to satisfy the selection. If all selection elements can be satisfied one way or another, then this method returns true. If one or more selection elements could never be satisfied given the input full trigger list, then this method returns false. At some level, this method tests whether the selection list is a "subset" of the full trigger list.

Parameters

pathspec	The trigger selection list (vector of string).
fullTriggerList	The full list of trigger names (vector of string).

Returns: true if the selection list is valid, false otherwise.

Definition at line 460 of file EventSelector.cc.

References acceptEvent(), edm::hlt::Exception, edm::hlt::Fail, customizeTrackingMonitorSeedNumber::idx, edm::hlt::Pass, edm::HLTGlobalStatus::reset(), and wantAll().

Referenced by testSelectionOverlap().

   {
     // an empty selection list is not valid
     // (we default an empty "SelectEvents" parameter to {"*","!*"} in
     // the getEventSelectionVString method below to help avoid this)
     if (pathspec.size() == 0)
     {
       return false;
     }
 
     // loop over each element in the selection list
     for (unsigned int idx = 0; idx < pathspec.size(); idx++)
     {
       Strings workingList;
       workingList.push_back(pathspec[idx]);
 
       // catch exceptions from the EventSelector constructor
       // (and anywhere else) and mark those as failures.
       // The EventSelector constructor seems to do the work of
       // checking if the selection is outside the full trigger list.
       try
       {
         // create an EventSelector instance for this selection
         EventSelector evtSelector(workingList, fullTriggerList);
 
         // create the TriggerResults instance that we'll use for testing
         unsigned int fullTriggerCount = fullTriggerList.size();
         HLTGlobalStatus hltGS(fullTriggerCount);
         TriggerResults sampleResults(hltGS, fullTriggerList);
 
         // loop over each path
         bool oneResultMatched = false;
         for (unsigned int iPath = 0; iPath < fullTriggerCount; iPath++)
         {
           // loop over the possible values for the path status
           for (int iState = static_cast<int>(hlt::Pass);
                iState <= static_cast<int>(hlt::Exception);
                iState++)
           {
             sampleResults[iPath] = HLTPathStatus(static_cast<hlt::HLTState>(iState), 0);
             if (evtSelector.wantAll() || evtSelector.acceptEvent(sampleResults))
             {
               oneResultMatched = true;
               break;
             }
 
             sampleResults.reset(iPath);
           }
 
           if (oneResultMatched) break;
         }
 
     // Finally, check in case the selection element was a wildcarded 
     // negative such as "!*":
     
         if (!oneResultMatched)  {
       for (unsigned int iPath = 0; iPath < fullTriggerCount; iPath++) {
             sampleResults[iPath] = HLTPathStatus(hlt::Fail, 0);
           }
       if (evtSelector.acceptEvent(sampleResults)) {
               oneResultMatched = true;
           }
         }
     
         // if none of the possible trigger results matched the
         // selection element, then we declare the whole selection
         // list invalid
         if (!oneResultMatched)
         {
           return false;
         }
       }
       catch (edm::Exception const& excpt)
       {
         return false;
       }
     }
 
     // if we made it to this point, then it must have been possible
     // to satisfy every selection element one way or another
     return true;
   }

bool EventSelector::subset	(	std::vector< bool > const &	a,
		std::vector< bool > const &	b
	)

staticprivate

Definition at line 1077 of file EventSelector.cc.

References i.

Referenced by testSelectionOverlap().

   {
     if (a.size() != b.size()) return false;
     // First test whether a is a non-empty subset of b 
     bool aPresent = false;
     bool aSubset = true;
     for (unsigned int i = 0; i != a.size(); ++i) {
       if (a[i]) {
         aPresent = true;
         if (!b[i]) {
       aSubset = false;
       break; 
     }
       }
     }   
     if (!aPresent) return false;
     if (aSubset) return true;
     
     // Now test whether b is a non-empty subset of a 
     bool bPresent = false;
     bool bSubset = true;
     for (unsigned int i = 0; i != b.size(); ++i) {
       if (b[i]) {
         bPresent = true;
         if (!a[i]) {
       bSubset = false;
       break; 
     }
       }
     }   
     if (!bPresent) return false;
     if (bSubset) return true;
  
     return false;                        
   } // subset

evtSel::OverlapResult EventSelector::testSelectionOverlap	(	Strings const &	pathspec1,
		Strings const &	pathspec2,
		Strings const &	fullTriggerList
	)

static

Tests if the specified trigger selection lists (path specs) overlap, where "overlap" means that a valid trigger result (given the full trigger list) could satisfy both selections.

Parameters

pathspec1	The first trigger selection list (vector of string).
pathspec2	The second trigger selection list (vector of string).
fullTriggerList	The full list of trigger names (vector of string).

Returns: OverlapResult which indicates the degree of overlap.

Definition at line 556 of file EventSelector.cc.

References a, absolute_acceptors_, all_must_fail_, all_must_fail_noex_, b, combine(), conditional_acceptors_, edm::evtSel::ExactMatch, exception_acceptors_, expandDecisionList(), f, g, identical(), edm::evtSel::InvalidSelection, m, N, edm::evtSel::NoOverlap, muon::overlap(), overlapping(), edm::evtSel::PartialOverlap, selectionIsValid(), and subset().

   {
     bool overlap = false;
     
     // first, test that the selection lists are valid
     if (!selectionIsValid(pathspec1, fullTriggerList) ||
         !selectionIsValid(pathspec2, fullTriggerList))
     {
       return evtSel::InvalidSelection;
     }
  
     // catch exceptions from the EventSelector constructor
     // (and anywhere else) and mark those as failures
     try
     {
       // create an EventSelector instance for each selection list
       EventSelector a(pathspec1, fullTriggerList);
       EventSelector b(pathspec2, fullTriggerList);
 
       unsigned int N = fullTriggerList.size();
 
       // create the expanded masks for the various decision lists in a and b
       std::vector<bool> 
         aPassAbs = expandDecisionList(a.absolute_acceptors_,true,N);
       std::vector<bool> 
         aPassCon = expandDecisionList(a.conditional_acceptors_,true,N);
       std::vector<bool> 
         aFailAbs = expandDecisionList(a.absolute_acceptors_,false,N);
       std::vector<bool> 
         aFailCon = expandDecisionList(a.conditional_acceptors_,false,N);
       std::vector<bool> 
         aExc = expandDecisionList(a.exception_acceptors_,true,N);
       std::vector< std::vector<bool> > aMustFail;
       for (unsigned int m = 0; m != a.all_must_fail_.size(); ++m) {
         aMustFail.push_back(expandDecisionList(a.all_must_fail_[m],false,N));
       }
       std::vector< std::vector<bool> > aMustFailNoex;
       for (unsigned int m = 0; m != a.all_must_fail_noex_.size(); ++m) {
         aMustFailNoex.push_back 
         (expandDecisionList(a.all_must_fail_noex_[m],false,N));
       }
 
       std::vector<bool> 
         bPassAbs = expandDecisionList(b.absolute_acceptors_,true,N);
       std::vector<bool> 
         bPassCon = expandDecisionList(b.conditional_acceptors_,true,N);
       std::vector<bool> 
         bFailAbs = expandDecisionList(b.absolute_acceptors_,false,N);
       std::vector<bool> 
         bFailCon = expandDecisionList(b.conditional_acceptors_,false,N);
       std::vector<bool> 
         bExc = expandDecisionList(b.exception_acceptors_,true,N);
       std::vector< std::vector<bool> > bMustFail;
       for (unsigned int m = 0; m != b.all_must_fail_.size(); ++m) {
         bMustFail.push_back(expandDecisionList(b.all_must_fail_[m],false,N));
       }
       std::vector< std::vector<bool> > bMustFailNoex;
       for (unsigned int m = 0; m != b.all_must_fail_noex_.size(); ++m) {
         bMustFailNoex.push_back 
         (expandDecisionList(b.all_must_fail_noex_[m],false,N));
       }
 
       std::vector<bool> aPass = combine(aPassAbs, aPassCon);
       std::vector<bool> bPass = combine(bPassAbs, bPassCon);
       std::vector<bool> aFail = combine(aFailAbs, aFailCon);
       std::vector<bool> bFail = combine(bFailAbs, bFailCon);
 
       // Check for overlap in the primary masks
       overlap = overlapping(aPass, bPass) || 
             overlapping(aFail, bFail) || 
             overlapping(aExc, bExc);
       if (overlap) return identical(a,b,N) ? evtSel::ExactMatch 
                              : evtSel::PartialOverlap;
 
       // Check for overlap of a primary fail mask with a must fail mask
       for (unsigned int f = 0; f != aMustFail.size(); ++f) {
         overlap = overlapping(aMustFail[f], bFail);
     if (overlap) return evtSel::PartialOverlap;
     for (unsigned int g = 0; g != bMustFail.size(); ++g) {
           overlap = subset(aMustFail[f], bMustFail[g]);
       if (overlap) return evtSel::PartialOverlap;
     } 
     for (unsigned int g = 0; g != bMustFailNoex.size(); ++g) {
           overlap = subset(aMustFail[f], bMustFailNoex[g]);
       if (overlap) return evtSel::PartialOverlap;
     }
       }
       for (unsigned int f = 0; f != aMustFailNoex.size(); ++f) {
         overlap = overlapping(aMustFailNoex[f], bFail);
     if (overlap) return evtSel::PartialOverlap;
     for (unsigned int g = 0; g != bMustFail.size(); ++g) {
           overlap = subset(aMustFailNoex[f], bMustFail[g]);
       if (overlap) return evtSel::PartialOverlap;
     } 
     for (unsigned int g = 0; g != bMustFailNoex.size(); ++g) {
           overlap = subset(aMustFailNoex[f], bMustFailNoex[g]);
       if (overlap) return evtSel::PartialOverlap;
     }
       }
       for (unsigned int g = 0; g != bMustFail.size(); ++g) {
         overlap = overlapping(bMustFail[g], aFail);
     if (overlap) return evtSel::PartialOverlap;
       }
       for (unsigned int g = 0; g != bMustFailNoex.size(); ++g) {
         overlap = overlapping(bMustFail[g], aFail);
     if (overlap) return evtSel::PartialOverlap;
       }
 
     }
     catch (edm::Exception const& excpt)
     {
       return evtSel::InvalidSelection;
     }
 
     // If we get to here without overlap becoming true, there is no overlap
 
     return evtSel::NoOverlap;
 
   } // testSelectionOverlap