make_nue_thirdana_pot.C File Reference

#include <bitset>
#include <iostream>
#include "CAFAna/Analysis/Exposures.h"
#include "CAFAna/Cuts/Cuts.h"
#include "CAFAna/Cuts/TimingCuts.h"
#include "CAFAna/Cuts/NueCutsFirstAna.h"
#include "CAFAna/Cuts/SpillCuts.h"
#include "CAFAna/Core/Loaders.h"
#include "CAFAna/Core/Spectrum.h"
#include "CAFAna/Core/SpectrumLoader.h"
#include "CAFAna/Prediction/IPrediction.h"
#include "CAFAna/Prediction/PredictionGenerator.h"
#include "CAFAna/Vars/HistAxes.h"
#include "CAFAna/Core/EventList.h"
#include "StandardRecord/Proxy/SRProxy.h"
#include "TFile.h"
#include "TH1D.h"

Go to the source code of this file.

Functions
void	make_nue_thirdana_pot ()

Variables
const SpillVar	kSpillTime ([](const caf::SRSpillProxy *spill){return spill->spilltimesec;})
const SpillVar	kRun ([](const caf::SRSpillProxy *spill){return spill->run;})
const SpillVar	fPOT ([](const caf::SRSpillProxy *spill){return spill->spillpot;})
const SpillVar	kLiveMass ([](const caf::SRSpillProxy *spill){std::string binary=std::bitset< 14 >(spill->dibmask).to_string();double live=12e-6;if(spill->run< 17946) live=24e-6;double sum=0;int dbcount=1;int mydb=0;for(int i=13;i >-1;--i){int test=binary[i]- '0';if(dbcount >=spill->dibfirst && dbcount<=spill->diblast && test==1) mydb++;if(dbcount >=spill->dibfirst && dbcount<=spill->diblast && test==0){if(mydb >=4) sum+=live *kAna2017FiducialMass[mydb];mydb=0;}dbcount++;}if(mydb >=4) sum+=live *kAna2017FiducialMass[mydb];return sum;})
const SpillVar	kLive ([](const caf::SRSpillProxy *spill){if(spill->run< 17946) return 24e-6;else return 12e-6;})
const SpillVar	fPOTMass ([](const caf::SRSpillProxy *spill){std::string binary=std::bitset< 14 >(spill->dibmask).to_string();double pot=0;int dbcount=1;int mydb=0;for(int i=13;i >-1;--i){int test=binary[i]- '0';if(dbcount >=spill->dibfirst && dbcount<=spill->diblast && test==1) mydb++;if(dbcount >=spill->dibfirst && dbcount<=spill->diblast && test==0){if(mydb >=4) pot+=spill->spillpot *kAna2017FiducialMass[mydb];mydb=0;}dbcount++;}if(mydb >=4) pot+=spill->spillpot *kAna2017FiducialMass[mydb];return pot;})
const Var	kSliceTime ([](const caf::SRProxy *sr){return sr->hdr.unixtime;})

Function Documentation

void make_nue_thirdana_pot

(

)

Definition at line 100 of file make_nue_thirdana_pot.C.

References ana::SpectrumLoaderBase::AddSpillHistogram(), fPOT, fPOTMass, ana::SpectrumLoader::Go(), ana::kCosmic, kLive, kLiveMass, ana::kRun, kSpillTime, ana::kStandardSpillCuts, outFile, and ana::SpectrumLoaderBase::SetSpillCut().

{
  TH1::SetDefaultSumw2();

  //binned in unix time encompassing first analysis period
  TH1D* hLiveMass = new TH1D("hLiveMass","livetime exposure;run;exposure (s #times kT)",12473, 12939.5, 25412.5);
  TH1D* hLive = new TH1D("hLive","livetime;run;livetime (s)",12473, 12939.5, 25412.5);
  TH1D* hPOTMass = new TH1D("hPOTMass","exposure;run;exposure (POT #times kT)",12473, 12939.5, 25412.5);
  TH1D* hPOT = new TH1D("hPOT","POT;run;POT",12473, 12939.5, 25412.5);
  TH1D* hLiveMassT = new TH1D("hLiveMassT","livetime exposure;run;exposure (s #times kT)",9900,1391000000,1490000000);
  TH1D* hLiveT = new TH1D("hLiveT","livetime;run;livetime (s)",9900,1391000000,1490000000);
  TH1D* hPOTMassT = new TH1D("hPOTMassT","exposure;run;exposure (POT #times kT)",9900,1391000000,1490000000);
  TH1D* hPOTT = new TH1D("hPOTT","POT;run;POT",9900,1391000000,1490000000);
  

  //load data files
  SpectrumLoader loaderFDData("/pnfs/nova/persistent/production/concat/R17-03-01-prod3reco.k/full/prod_decaf_R17-03-01-prod3reco.k_fd_numi_fhc_full_v1_goodruns_nue2017/*.root", kCosmic);

  loaderFDData.SetSpillCut(kStandardSpillCuts);
  loaderFDData.AddSpillHistogram(hLiveMass,kRun,kStandardSpillCuts, kLiveMass);
  loaderFDData.AddSpillHistogram(hLive,kRun,kStandardSpillCuts, kLive);
  loaderFDData.AddSpillHistogram(hPOTMass,kRun,kStandardSpillCuts, fPOTMass);
  loaderFDData.AddSpillHistogram(hPOT,kRun,kStandardSpillCuts, fPOT);
  loaderFDData.AddSpillHistogram(hLiveMassT,kSpillTime,kStandardSpillCuts, kLiveMass);
  loaderFDData.AddSpillHistogram(hLiveT,kSpillTime,kStandardSpillCuts, kLive);
  loaderFDData.AddSpillHistogram(hPOTMassT,kSpillTime,kStandardSpillCuts, fPOTMass);
  loaderFDData.AddSpillHistogram(hPOTT,kSpillTime,kStandardSpillCuts, fPOT);
  loaderFDData.Go();

  TFile* outFile = new TFile("nue_third_ana_pot.root", "RECREATE");
  hLiveMass->Write();
  hLive->Write();
  hPOTMass->Write();
  hPOT->Write();
  hLiveMassT->Write();
  hLiveT->Write();
  hPOTMassT->Write();
  hPOTT->Write();
  outFile->Close();

}

Variable Documentation

const SpillVar fPOT([](const caf::SRSpillProxy *spill){return spill->spillpot;})

Referenced by ana::_preview::NewSpectrum::HandleExposure(), and make_nue_thirdana_pot().

const SpillVar fPOTMass([](const caf::SRSpillProxy *spill){std::string binary=std::bitset< 14 >(spill->dibmask).to_string();double pot=0;int dbcount=1;int mydb=0;for(int i=13;i >-1;--i){int test=binary[i]- '0';if(dbcount >=spill->dibfirst && dbcount<=spill->diblast && test==1) mydb++;if(dbcount >=spill->dibfirst && dbcount<=spill->diblast && test==0){if(mydb >=4) pot+=spill->spillpot *kAna2017FiducialMass[mydb];mydb=0;}dbcount++;}if(mydb >=4) pot+=spill->spillpot *kAna2017FiducialMass[mydb];return pot;})

Referenced by make_nue_thirdana_pot().

const SpillVar kLive([](const caf::SRSpillProxy *spill){if(spill->run< 17946) return 24e-6;else return 12e-6;})

Referenced by make_nue_thirdana_pot().

const SpillVar kLiveMass([](const caf::SRSpillProxy *spill){std::string binary=std::bitset< 14 >(spill->dibmask).to_string();double live=12e-6;if(spill->run< 17946) live=24e-6;double sum=0;int dbcount=1;int mydb=0;for(int i=13;i >-1;--i){int test=binary[i]- '0';if(dbcount >=spill->dibfirst && dbcount<=spill->diblast && test==1) mydb++;if(dbcount >=spill->dibfirst && dbcount<=spill->diblast && test==0){if(mydb >=4) sum+=live *kAna2017FiducialMass[mydb];mydb=0;}dbcount++;}if(mydb >=4) sum+=live *kAna2017FiducialMass[mydb];return sum;})

Referenced by make_nue_thirdana_pot().

const SpillVar kRun([](const caf::SRSpillProxy *spill){return spill->run;})

Referenced by resolution_event_tree().

const Var kSliceTime([](const caf::SRProxy *sr){return sr->hdr.unixtime;})

const SpillVar kSpillTime([](const caf::SRSpillProxy *spill){return spill->spilltimesec;})

Referenced by drawLongTerm(), drawLongTerm_2019(), getStabilitySpectra(), and make_nue_thirdana_pot().