1 #include <boost/regex.hpp>
4 #include "TTreeCache.h"
15 #define private public
24 m_name(name), m_file(0), m_eventTree(0), m_event(0),
25 m_needUpdate(
true), m_globalEventList(0)
33 delete (*i)->m_eventList;
40 gErrorIgnoreLevel = 3000;
41 TFile *newFile = TFile::Open(
m_name.c_str());
42 if (newFile == 0 || newFile->IsZombie() || !newFile->Get(
"Events")) {
45 throw std::runtime_error(
"Invalid file. Ignored.");
47 gErrorIgnoreLevel = -1;
53 typedef std::vector<edm::ProcessConfiguration> provList;
55 TTree *metaData =
dynamic_cast<TTree*
>(
m_file->Get(
"MetaData"));
56 TBranch *
b = metaData->GetBranch(
"ProcessConfiguration");
63 int currentVersionArr[] = {0, 0, 0};
64 for (provList::iterator
i = x->begin();
i != x->end(); ++
i)
67 TString dcv =
i->releaseVersion();
69 int nvv = currentVersionArr[0]*100 + currentVersionArr[1]*10 + currentVersionArr[2];
70 if (nvv > latestVersion) {
84 TString
msg = Form(
"incompatible data: Process version does not mactch major data formats version. File produced with %s. Data formats version \"CMSSW_%d_%d_%d\".\n",
86 msg +=
"Use --no-version-check option if you still want to view the file.\n";
87 throw std::runtime_error(msg.Data());
97 throw std::runtime_error(
"Cannot find TTree 'Events' in the data file");
111 printf(
"Reading %lld bytes in %d transactions.\n",
159 Long64_t val = tree_entry;
174 Long64_t val = tree_entry;
176 if (list[idx] == val)
189 if ((*it)->m_selector->m_enabled)
209 if ((*it)->m_selector->m_enabled && (*it)->m_needsUpdate)
215 if ((*it)->m_selector->m_enabled)
217 if ((*it)->hasSelectedEvents())
266 boost::regex re(
std::string(
"\\$") + (*i)->name());
272 interpretedSelection = boost::regex_replace(interpretedSelection, re,
273 fullBranchName +
".obj");
287 TObjArray* branches =
m_eventTree->GetListOfBranches();
288 std::vector<void*> previousBranchAddresses;
289 previousBranchAddresses.reserve(branches->GetEntriesFast());
291 std::auto_ptr<TIterator> pIt( branches->MakeIterator());
292 while(TObject* branchObj = pIt->Next()) {
293 TBranch*
b =
dynamic_cast<TBranch*
> (branchObj);
295 const char *
name = b->GetName();
296 unsigned int length = strlen(name);
297 if(length > 1 && name[length-1]!=
'.') {
299 previousBranchAddresses.push_back(0);
303 if(0!=b->GetAddress()) {
306 previousBranchAddresses.push_back(b->GetAddress());
308 previousBranchAddresses.push_back(0);
317 Int_t
result =
m_eventTree->Draw(
">>fworks_filter", interpretedSelection.c_str());
333 std::auto_ptr<TIterator> pIt( branches->MakeIterator());
334 std::vector<void*>::const_iterator itAddress = previousBranchAddresses.begin();
335 while(TObject* branchObj = pIt->Next()) {
336 TBranch*
b =
dynamic_cast<TBranch*
> (branchObj);
337 if(0!=b && 0!=*itAddress) {
338 b->SetAddress(*itAddress);
354 boost::regex re_spaces(
"\\s+");
355 selection = boost::regex_replace(selection,re_spaces,
"");
356 if (selection.find(
"&&") != std::string::npos &&
357 selection.find(
"||") != std::string::npos )
381 bool junction_mode =
true;
382 if (selection.find(
"||")!=std::string::npos)
383 junction_mode =
false;
385 boost::regex re(
"\\&\\&|\\|\\|");
387 boost::sregex_token_iterator
i(selection.begin(), selection.end(), re, -1);
388 boost::sregex_token_iterator
j;
391 std::vector<std::pair<unsigned int,bool> >
filters;
397 if (filter[0] ==
'!')
400 filter.erase(filter.begin());
403 if (index == triggerNames->
size())
408 filters.push_back(std::make_pair(index, flag));
426 std::vector<std::pair<unsigned int,bool> >::const_iterator
filter = filters.begin();
427 bool passed = hTriggerResults->
accept(filter->first) == filter->second;
428 while (++filter != filters.end())
431 passed &= hTriggerResults->
accept(filter->first) == filter->second;
433 passed |= hTriggerResults->
accept(filter->first) == filter->second;
FWEventSelector * m_selector
virtual void Enter(Long64_t entry)
virtual edm::TriggerNames const & triggerNames(edm::TriggerResults const &triggerResults) const
bool accept() const
Has at least one path accepted the event?
FWTEventList * m_eventList
void getPrimaryData() const
std::string m_triggerProcess
FWFileEntry(const std::string &name, bool checkVersion)
void getDecomposedVersion(const TString &s, int *out)
Strings::size_type size() const
const std::string & processName() const
virtual std::string const getBranchNameFor(std::type_info const &, char const *iModuleLabel, char const *iProductInstanceLabel, char const *iProcessName) const
Return the branch name in the TFile which contains the data.
bool isEventSelected(int event)
std::string const & processName() const
void getByLabel(const P &iP, const char *iModuleLabel, const char *iProductInstanceLabel=0, const char *iProcessLabel=0)
Event const & toBegin()
Go to the very first Event.
const std::string & productInstanceLabel() const
std::list< Filter * > m_filterEntries
int * supportedDataFormatsVersion()
unsigned int triggerIndex(std::string const &name) const
virtual void Add(const TEventList *list)
const TClass * type() const
bool filterEventsWithCustomParser(Filter *filter)
bool to(Long64_t iIndex)
Go to the event at index iIndex.
std::list< Filter * > & filters()
virtual bool atEnd() const
const_iterator begin() const
NOTE: iterator is allowed to return a null object for items that have been removed.
tuple idx
DEBUGGING if hasattr(process,"trackMonIterativeTracking2012"): print "trackMonIterativeTracking2012 D...
ReleaseVersion const & releaseVersion() const
void runFilter(Filter *fe, const FWEventItemsManager *eiMng)
int nextSelectedEvent(int event)
const edm::EventBase * m_event
bool acceptDataFormatsVersion(TString &n)
void updateFilters(const FWEventItemsManager *eiMng, bool isOR)
FWTEventList * m_globalEventList
int previousSelectedEvent(int event)
std::vector< FWEventItem * >::const_iterator const_iterator
const std::string & moduleLabel() const
How EventSelector::AcceptEvent() decides whether to accept an event for output otherwise it is excluding the probing of A single or multiple positive and the trigger will pass if any such matching triggers are PASS or EXCEPTION[A criterion thatmatches no triggers at all is detected and causes a throw.] A single negative with an expectation of appropriate bit checking in the decision and the trigger will pass if any such matching triggers are FAIL or EXCEPTION A wildcarded negative criterion that matches more than one trigger in the trigger list("!*","!HLTx*"if it matches 2 triggers or more) will accept the event if all the matching triggers are FAIL.It will reject the event if any of the triggers are PASS or EXCEPTION(this matches the behavior of"!*"before the partial wildcard feature was incorporated).Triggers which are in the READY state are completely ignored.(READY should never be returned since the trigger paths have been run
const_iterator end() const