NCType.cxx

Go to the documentation of this file.

// $Id: NCType.cxx,v 1.39 2009/06/18 13:32:49 bckhouse Exp $
//
// NCType
//
// Author:  B. Rebel 10/2005
//
#include "NCUtils/NCType.h"

#include "AnalysisNtuples/ANtpHeaderInfo.h"

#include "NCUtils/NCUtility.h"

#include <cassert>

#include "MessageService/MsgService.h"
CVSID("$Id: NCType.cxx,v 1.39 2009/06/18 13:32:49 bckhouse Exp $");


//-----------------------------------------------------------------------
NCType::EFileType NCType::FindFileType(ANtpHeaderInfo* header)
{
  assert(header);
  //check if the event is from a tau file.  the mc run numbers have the convention
  //<ikey><itgt><iflv><S><RRRR>
  //<ikey> = 1 near, 2 far
  //<itgt> = 1 in detector only, 2 in rock only, 3 detector + rock, 4 in fiducial only
  //<iflv> = 0 beam, 1 nue, 3 nutau, 4 numu, nue inverted, 5 far oscillated mock data
  //<S> = split set
  //<RRRR> = run number

  if(header->run > 21300000 &&
     header->run < 21400000) return NCType::kTauFile;

  //dont know what these runs are just now
  if(header->run > 21400000 &&
     header->run < 21510000) return NCType::kElectronFile;

  if(header->run > 21910000 &&
     header->run < 21910004) return NCType::kMockFile;

  return NCType::kBeamFile;
}

TString NCType::FileTypeAsString(EFileType f)
{
  switch(f){
  case kUnknownFile:  return "Unknown filetype";
  case kBeamFile:     return "Beam";
  case kTauFile:      return "Tau";
  case kElectronFile: return "Electron";
  case kMockFile:     return "Mock";
  default:
    assert(0 && "This really shouldn't happen");
  }
}

vector<BeamType::BeamType_t> NCType::BeamListFromString(const char* str)
{
  const vector<TString> sl = NC::Utility::ParseStringList(str);

  using BeamType::BeamType_t;

  vector<BeamType_t> ret;

  for(unsigned int n = 0; n < sl.size(); ++n){
    // BeamType doesn't have a FromString function
    for(BeamType_t beam = BeamType_t(0);
        beam < BeamType::kEndOfList;
        beam = BeamType_t(int(beam)+1)){
      if(sl[n] == BeamType::AsString(beam)) ret.push_back(beam);
    } // end for beam
  } // end for n

  return ret;
}

---