Params.cxx

Go to the documentation of this file.

#include <iostream>
#include <fstream>
#include <string>
#include "TSystem.h"
#include "MessageService/MsgService.h"
#include "Params.h"

using std::cout;
using std::endl;
using std::ifstream;
using std::ofstream;

// ************************* STATIC DATA MEMBERS ************************8

CVSID("$Id: Params.cxx,v 1.3 2007/03/01 17:07:27 rhatcher Exp $");

// ************************** CONSTRUCTOR & DESTRUCTOR *************************
Params::Params(){

  //1st entry is neutrinos and second is antineutrinos
  IsFid[0]=1;          IsFid[1]=1;
  PittTrkQual[0]=1;    PittTrkQual[1]=1;
  FitSig[0]=0.1;       FitSig[1]=0.1;
  DavidPID[0]=0.1;     DavidPID[1]=-2;
  NuBarPID[0]=-0.5;    NuBarPID[1]=0.8;

}

Params::Params(const char* name){

  //Default values
  //1st entry is neutrinos and second is antineutrinos
  IsFid[0]=1;          IsFid[1]=1;
  PittTrkQual[0]=1;    PittTrkQual[1]=1;
  FitSig[0]=0.1;       FitSig[1]=0.1;
  DavidPID[0]=0.1;     DavidPID[1]=-2;
  NuBarPID[0]=-0.5;    NuBarPID[1]=0.8;
  
  ifstream file_sel (name);
  string line;
  
  if(file_sel.is_open()){
    while( !file_sel.eof() ){
      getline(file_sel,line);
      cout << line << endl;
    }
    file_sel.close();
  } else {
    MSG("NuBarPID",Msg::kWarning) << "Could not open file" << endl;
  }

}

Params::~Params(){}

// ************************* PUBLIC METHODS ************************
// -----------------------------------------------------------------
Int_t Params::WantIsFid(Int_t nunubar){

  if(nunubar==0 || nunubar==1){
    return IsFid[nunubar];
  }else {
    MSG("NuBarPID",Msg::kWarning) << "Entry must be 0 (neutrino selection) or 1 (antineutrino selection)" << endl;
    return -1;
  }
}

// -----------------------------------------------------------------
Int_t Params::WantPittTrkQual(Int_t nunubar){

  if(nunubar==0 || nunubar==1){
    return PittTrkQual[nunubar];
  }else {
    MSG("NuBarPID",Msg::kWarning) << "Entry must be 0 (neutrino selection) or 1 (antineutrino selection)" << endl;
    return -1;
  }
}

// -----------------------------------------------------------------
Float_t Params::FitSigCut(Int_t nunubar){

  if(nunubar==0 || nunubar==1){
    return FitSig[nunubar];
  }else {
    MSG("NuBarPID",Msg::kWarning) << "Entry must be 0 (neutrino selection) or 1 (antineutrino selection)" << endl;
    return -100;
  }
}

// -----------------------------------------------------------------
Float_t Params::DavidPIDCut(Int_t nunubar){

  if(nunubar==0 || nunubar==1){
    return DavidPID[nunubar];
  }else {
    MSG("NuBarPID",Msg::kWarning) << "Entry must be 0 (neutrino selection) or 1 (antineutrino selection)" << endl;
    return -100;
  }
}

// -----------------------------------------------------------------
Float_t Params::NuBarPIDCut(Int_t nunubar){

  if(nunubar==0 || nunubar==1){
    return NuBarPID[nunubar];
  }else {
    MSG("NuBarPID",Msg::kWarning) << "Entry must be 0 (neutrino selection) or 1 (antineutrino selection)" << endl;
    return -100;
  }
}

---