NCRunUtil.cxx

Go to the documentation of this file.

#include "NCRunUtil.h"
using namespace NC::RunUtil;

#include "AnalysisNtuples/ANtpHeaderInfo.h"
#include "AnalysisNtuples/ANtpRecoInfo.h"
#include "Conventions/Detector.h"
#include "Conventions/SimFlag.h"
#include "MessageService/MsgService.h"

#include "TRegexp.h"

#include <cassert>
#include <map>
#include <vector>
using namespace std;

CVSID("$Id: NCRunUtil.cxx,v 1.7 2010/01/19 17:02:22 rodriges Exp $");


//-----------------------------------------------------------------------
ERunType NC::RunUtil::FindRunType(ANtpHeaderInfo* header, ANtpRecoInfo* reco)
{
  assert(header);

  if(header->dataType == SimFlag::kData){
    if(header->detector == Detector::kFar){
      if     (header->run <= 35723) return kRunI;
      // According to Robert P., the last run in the run period is
      // 38449, while the last LE run is 38420. Be more inclusive here
      // in case we ever run on non-LE data
      else if(header->run > 35723 && header->run <= 38449) return kRunII;
      // TODO: Just used the end of run II as the start of run III,
      // but this includes shutdown. Probably not a problem
      //
      // The final run number comes from Justin
      else if(header->run > 38449 && header->run <= 43639) return kRunIII;
      assert(0 && "Unknown run");
    }
    else{
      // Near detector
      if     (header->run <= 10159) return kRunI;
      // According to Robert P., the last run in the run period is
      // 12623, while the last LE run is 12578. Be more inclusive here
      // in case we ever run on non-LE data
      else if(header->run > 10159 && header->run <= 12623) return kRunII;
      // TODO: Just used the end of run II as the start of run III,
      // but this includes shutdown. Probably not a problem
      //
      // The final run number comes from Andy's docdb 6572
      else if(header->run > 12623 && header->run <= 16502) return kRunIII;
      assert(0 && "Unknown run");
    } // end if near det
  } // end if data
  else{
    // MC
    assert(reco);
    return ERunType(reco->runPeriod);
  }
}

//-----------------------------------------------------------------------
//any data before march 2006 is run I, anything between may 2006 and august 2007
//is run II
ERunType NC::RunUtil::FindRunType(TString fileName)
{

  // Yay for string processing in C and its evil spawn. In perl these
  // five lines would be:
  // /(20\d\d)-(\d\d)/ and $epochMonth=100*$1+$2;
  // or something like that

  // Match something that looks like a month in yyyy-mm format between 2000 and 2100
  TRegexp re("20[0-9][0-9]-[0-9][0-9]");
  // TString::operator()(TRegexp) returns the TSubString that matches the regexp
  TString date=TString(fileName(re));
  if(date.Length()==7){
    // The filename had a month number in it - must be a data file
    int year=TString(date(0, 4)).Atoi();
    // Brief worry that this would try octal if the month has a
    // leading zero, but man 3 atoi assures me that it always uses base 10
    int month=TString(date(5, 2)).Atoi();
    
    int epochMonth=100*year+month;
    
    // May 2005 - Feb 2006: Run I
    // May 2006 - Aug 2007: Run II
    // Nov 2007 - Jun 2009: Run III
    if (epochMonth >= 2005*100 + 5  && epochMonth <= 2006*100 + 2) return kRunI;
    if (epochMonth >= 2006*100 + 5  && epochMonth <= 2007*100 + 8) return kRunII;
    if (epochMonth >= 2007*100 + 11 && epochMonth <= 2006*100 + 6) return kRunIII;

    assert(0 && "fileName with date outside known range");
  }
  else{
    // Must be an MC file. No harm in checking
    assert(fileName.Contains("mc") ||
           fileName.Contains("tau") ||
           fileName.Contains("electron") ||
           fileName.Contains("mock"));

    TRegexp runRE("run[0-9]");
    TString runPeriodStr=TString(fileName(runRE));
    assert(runPeriodStr.Length()==4);
    int runPeriod=TString(runPeriodStr(3,1)).Atoi();
    return ERunType(runPeriod);
  }
}

---