Event.cpp

Go to the documentation of this file.

#include "cafe/Event.hpp"

namespace cafe {
    
  Event::Event()
  {
  }

  bool Event::isMC() const
  //
  // Purpose: Tell if event is mc.  The algorithm is similar to 
  //          HistorySelector::is_monte_carlo, namely, it checks a) history,
  //          b) TMBMCevtInfo::nchunks() > 0, or c) run number == 1 
  //          (else data).
  //
  {
    bool is_mc = false;

    // First check run number in TMBGlobal.  Assume run 1 is always mc.

    const TMBGlobal* global = getGlobal();
    is_mc = global != 0 && global->runno() == 1;
      
    // Next look for program "d0gstar" in TMBHistory.

    if(!is_mc) {
      const Collection<TMBHistory> history = getHistory();
      for(Collection<TMBHistory>::const_iterator i = history.begin();
          i != history.end(); ++i) {
        const TMBHistory& hist = *i;
        if(hist.program() == std::string("d0gstar")) {
          is_mc = true;
          break;
        }
      }
    }

    // Finally check if TMBMCevtInfo::nchunks() > 0.

    if(!is_mc) {
      const TMBMCevtInfo* mcinfo = getMCEventInfo();
      is_mc = mcinfo != 0 && mcinfo->nchunks() > 0;
    }

    // Done

    return is_mc;
  }

  bool Event::isRun2b() const
  {
    bool is_2b = false;

    // Data or mc?

    if(isMC()) {

      // For mc, the run 2a vs. run 2b determination is based on the 
      // d0gstar program version from TMBHistory.  Any version later than
      // t05.09.00 or p20.00.00 is assumed to be run 2b.  It is a fatal
      // error if a d0gstar TMBHistory object isn't found.

      const TMBHistory* phist = 0;
      const Collection<TMBHistory> history = getHistory();
      for(Collection<TMBHistory>::const_iterator i = history.begin();
          i != history.end(); ++i) {
        phist = &*i;
        assert(phist != 0);
        if(phist->program() == std::string("d0gstar"))
          break;
      }
      if(phist == 0) {
        assert(phist != 0);
        abort();
      }
          
      std::string version_tag = phist->version();

      // We only know how to definitively interpret standard 9-character
      // test and production version tags: txx.yy.zz or pxx.yy.zz.  In
      // these cases, do a straight lexical comparison with the cutoff
      // version.  Assume any version tag that doesn't conform to these
      // standard forms is run2a (since there were a few early
      // production releases with weird names), but print an error
      // message.  Don't think anyone should be using them any more, but
      // you never know.

      if(version_tag.size() == 9 && 
         version_tag.substr(0,1) == std::string("t")) {

        // Check test release.

        is_2b = version_tag >= std::string("t05.09.00");
      }
      else if(version_tag.size() == 9 && 
              version_tag.substr(0,1) == std::string("p")) {

        // Check production release.

        is_2b = version_tag >= std::string("p20.00.00");
      }
      else {

        // Nonstandard version tag.

        std::cerr << "cafe::Event::isRun2b - Nonstandard d0gstar version "
                  << version_tag << ", assuming run 2a." << std::endl;
      }
    }
    else {

      // For data, the determination whether this is run 2a or run 2b is
      // based on run number, with the cutoff being run 219000.

      const TMBGlobal* global = getGlobal();
      if(global == 0) {
        assert(global != 0);
        abort();
      }
      is_2b = global->runno() >= 219000;
    }

    return is_2b;
  }
}

ClassImp(cafe::Event);

ClassImp(cafe::Collection<TMBHistory>);
ClassImp(cafe::Collection<TMBTrigger>);

ClassImp(cafe::Collection<TMBLorentzVector>);
ClassImp(cafe::Collection<TMBMuon>);
ClassImp(cafe::Collection<TMBEMCluster>);
ClassImp(cafe::Collection<TMBJet>);
ClassImp(cafe::Collection<TMBParticleJet>);
ClassImp(cafe::Collection<TMBTrack>);
ClassImp(cafe::Collection<TMBIsoTrack>);
ClassImp(cafe::Collection<TMBVertex>);
ClassImp(cafe::Collection<TMBPrimaryVertex>);
ClassImp(cafe::Collection<TMBTrackCal>);
ClassImp(cafe::Collection<TMBTrackCalJet>);
ClassImp(cafe::Collection<TMBTau>);

ClassImp(cafe::Collection<TMBCps>);
ClassImp(cafe::Collection<TMBCpsDigi>);
ClassImp(cafe::Collection<TMBFps>);
ClassImp(cafe::Collection<TMBFpsData>);
ClassImp(cafe::Collection<TMBFpd>);

ClassImp(cafe::Collection<TMBMCvtx>);
ClassImp(cafe::Collection<TMBMCpart>);

ClassImp(cafe::Collection<TMBL1CalTile>);
ClassImp(cafe::Collection<TMBL1CalTower>);
ClassImp(cafe::Collection<TMBL1Track>);
ClassImp(cafe::Collection<TMBL1toL2CTT>);
ClassImp(cafe::Collection<TMBL2Base>);
ClassImp(cafe::Collection<TMBL2CPS>);
ClassImp(cafe::Collection<TMBL2EM>);
ClassImp(cafe::Collection<TMBL2Muon>);
ClassImp(cafe::Collection<TMBL2Jet>);
ClassImp(cafe::Collection<TMBL2Track>);
ClassImp(cafe::Collection<TMBL2GblJet>);
ClassImp(cafe::Collection<TMBL2GblMuon>);
ClassImp(cafe::Collection<TMBL2GblEM>);
ClassImp(cafe::Collection<TMBL2GblTrack>);
ClassImp(cafe::Collection<TMBL2GblInvMass>);
ClassImp(cafe::Collection<TMBL2GblMJt>);
ClassImp(cafe::Collection<TMBL2GblHt>);

ClassImp(cafe::Collection<TMBL3Base>);
ClassImp(cafe::Collection<TMBL3BTagIP>);
ClassImp(cafe::Collection<TMBL3CFTVtx>);
ClassImp(cafe::Collection<TMBL3Ele>);
ClassImp(cafe::Collection<TMBL3IPTrack>);
ClassImp(cafe::Collection<TMBL3Isolation>);
ClassImp(cafe::Collection<TMBL3Jet>);
ClassImp(cafe::Collection<TMBL3Muon>);
ClassImp(cafe::Collection<TMBL3Photon>);
ClassImp(cafe::Collection<TMBL3Tau>);
ClassImp(cafe::Collection<TMBL3Track>);
ClassImp(cafe::Collection<TMBL3MEt>);

ClassImp(cafe::Collection<TMBLum>);
ClassImp(cafe::Collection<TMBLumV>);
ClassImp(cafe::Collection<TMBSimCAEP>);

ClassImp(cafe::Collection<EventWeight>);

ClassImp(cafe::Collection<TMBBTag>);
ClassImp(cafe::Collection<TMBBTagSLT>);
ClassImp(cafe::Collection<TMBBTagJLIP>);
ClassImp(cafe::Collection<TMBBTagSVT>);
ClassImp(cafe::Collection<TMBBTagNN>);

ClassImp(cafe::Collection<TMBLeBob>);
ClassImp(cafe::Collection<TMBHiPt>);

ClassImp(cafe::Collection<TMBTdc>);
ClassImp(cafe::Collection<TMBFPDTrack>);

---