FindDuplicateEvents.cpp

Go to the documentation of this file.

#include "caf_util/FindDuplicateEvents.hpp"
#include "cafe/Config.hpp"

using namespace std ;
using namespace cafe ;
 
namespace caf_util {

  FindDuplicateEvents::FindDuplicateEvents(const char *name) 
    : cafe::Processor(name), _vars("_runno", "_evtno")
 {
    // get parameters from the configuration file
    cafe::Config config(name);

    _nscreams   = config.get("Nscreams", 1);
    _ignoreMC   = config.get("IgnoreMC", true);
    _ndups = 0;
  };
  
  bool FindDuplicateEvents::processEvent(cafe::Event &event)
  {
    const TMBGlobal* global = event.getGlobal (_vars) ;
    assert (global);
    const int run = global->runno();
    const int evtno = global->evtno();
    
    // in lieu of an official way to check for MC, count MC chunks:
    const TMBMCevtInfo* pMC = event.getMCEventInfo ();
    int   nchunks = pMC? pMC->nchunks(): -1;
    if (debug() > 9) out()<<"run: "<<run<<" evt: "<<evtno<<" pMC: "<<(int)pMC<<" nc: "<<nchunks<<endl;
    if (_ignoreMC && nchunks > 0) return true;

    Key key (run, evtno);

    // look for it
    KeyIter pos = _keys.find (key);
    if (debug() > 1) {
      out()<<"FindDuplicatedEvents processor: size: "
           <<_keys.size()<<" ndups: "<<_ndups<<". run: "<<run<<" evt: "
           << evtno ;
      if (pos == _keys.end())
        out() << " New event." ;
      else 
        out() << " Duplicated event." ;
      out() << endl ;
    }
    // didn't find it in the set --> first appearance --> pass it through
    if (pos == _keys.end()) {
      _keys.insert (key);
      return true;
    }

    if (++_ndups <= _nscreams) {
      err()<<endl<<_line<<endl;
      err()<<"!  Warning: caf_util::FindDuplicateEvents["<<name()
           <<"] removed a duplicate event, run# "<<run<<", event# "<<evtno<<endl;
      err()<<_line<<endl<<endl;
    }
    return false;
  };

  void FindDuplicateEvents::finish() {
    if (_ndups > 0) {
      if (_nscreams > 0) out()<<endl<<_line<<endl;
      out()<<"WARNING: caf_util::FindDuplicateEvents["<<name()
           <<"] removed "<<_ndups<<" duplicate events."<<endl;
      if (_nscreams > 0) out()<<_line<<endl<<endl;
    }
  }

  const TString FindDuplicateEvents::_line
                ("!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@");
}
ClassImp(caf_util::FindDuplicateEvents);
  

---