NuCutter.cxx

Go to the documentation of this file.

#include "NtupleUtils/NuCutter.h"
#include "NtupleUtils/NuCutImps.h"

#include "MessageService/MsgService.h"

#include "TString.h"


CVSID("$Id: NuCutter.cxx,v 1.9 2010/02/01 15:24:36 sjc Exp $");

// Function to make new cut instances

NuCut *NuCutter::GetNewCutInstance(TString ana, const NuPlots *plots) const
{
  // If we are passed in a blank string, return nothing
  if (ana == "") return 0;
  
  // Only compare upper case
  ana.ToUpper();
  
  if (ana == "CC0325STD") {
    return new NuCutImps::CC0325Std(plots);
  } else if (ana == "CC0720TEST") {
    return new NuCutImps::CC0720Test(plots);
  } else if (ana == "CCAPRESEL") {
    return new NuCutImps::CCAPresel(plots);
  } else if (ana == "CCA" || ana == "CCA_NUBAR" || ana == "CCA_0" || ana == "CCA_1" ||
             ana == "CCA_2" || ana == "CCA_3" || ana == "CCA_4") {
    return new NuCutImps::CCA(plots,ana);
  } else if (ana == "CCA_NC") {
    return new NuCutImps::CCA_NC(plots);
  } else if (ana == "NMBBRAVO" || ana == "BRAVO") {
    return new NuCutImps::Bravo(plots);
  } else if (ana == "NMBCHARLIE" || ana == "CHARLIE") {
    return new NuCutImps::Charlie(plots);
  } else if (ana == "NMBDELTA" || ana == "DELTA") {
    return new NuCutImps::Delta(plots);
  } else if (ana == "NMBBRAVOPRIME" || ana == "BRAVOPRIME") {
    return new NuCutImps::BravoPrime(plots);
  } else if (ana == "CHAIRSOUND") {
    return new NuCutImps::ChairSound(plots);
  } else if (ana == "PASS" || ana == "PASSTHROUGH") {
    return new NuCutImps::PassThrough(plots);
  } else if (ana == "MSROCK_NOV09") {
    return new NuCutImps::MSRock_Nov09(plots);
  } else {
    // MSG("NuCutter",Msg::kError) << "Cut named " << ana << " is not recognised!" << endl;
    return 0;
  }
}


// 
// NuCutter::NuCutter(const TString &anaVersion) : fCut(0)
// {
//   // If anaVersion is all numbers, use the backwards compatability cut
//   if (anaVersion.IsDigit()) {
//     MSG("NuCutter",Msg::kInfo) << "Falling back to NuCuts-style Cut for anaVersion " << anaVersion << endl;
//     fCut = new NuCutImps::NuCutsSelection(anaVersion.Atoi(), 0);
//     
//   } else {
//     MSG("NuCutter",Msg::kFatal) << "AnaVersion " << anaVersion << " can not be initialised through NuCuts" << endl;
//   }
// }
//

NuCutter::NuCutter() : fCut(0)
{
  
}

NuCutter::NuCutter(NuCut *cut) : fCut(cut)
{
  
}


NuCutter::NuCutter(const TString &anaVersion, const NuPlots *plots) : fCut(0)
{
  // If anaVersion is all numbers, use the backwards compatability cut
  if (anaVersion.IsDigit()) {
    MSG("NuCutter",Msg::kInfo) << "Falling back to NuCuts-style Cut for anaVersion " << anaVersion << endl;
    fCut = new NuCutImps::NuCutsSelection(anaVersion.Atoi(), plots);
    
  } else {
    fCut = GetNewCutInstance(anaVersion, plots);
    if (fCut == 0)
      MSG("NuCutter",Msg::kFatal) << "AnaVersion " << anaVersion << " can not be initialised through NuCuts" << endl;
  }
}


void NuCutter::SetCut(NuCut *cut)
{
  if (fCut) delete fCut;
  fCut = cut;
}


void NuCutter::SetCut(const TString &anaVersion, const NuPlots *plots)
{
  SetCut(GetNewCutInstance(anaVersion, plots));
  // This function will not let you set an invalid cut - this prevents
  // errors from not validating the input
  if (fCut == 0)
    MSG("NuCutter",Msg::kFatal) << "AnaVersion " << anaVersion << " can not be initialised through NuCuts" << endl;
}


Bool_t NuCutter::MakeCuts(const NuEvent &nu)
{
  if (!fCut) return false;
  
  return fCut->MakeCuts(nu);
}


Bool_t NuCutter::MakePreselectionCuts(const NuEvent &nu)
{
  if (!fCut) return false;
  
  return fCut->MakePreselectionCuts(nu);
}


Bool_t NuCutter::MakeSelectionCuts(const NuEvent &nu)
{
  if (!fCut) return false;
  
  return fCut->MakeSelectionCuts(nu);
}


  

---