#include "includes.h"
#include "numu_tools.h"
#include "3FlavorAna/Plotting/NuePlotStyle.h"

Functions
void	HornSettings (int quantile, double &pot, double &livetime, Sign::Sign_t &wrongs, double &maxy)

void	saveFDMCHists (std::string anaName="fhc", bool withExtrap=true)

Variables
const double	pot_fhc = kAna2018FHCPOT

const double	pot_rhc = kAna2018RHCPOT

const double	livetime_fhc = kAna2018FHCLivetime

const double	livetime_rhc = kAna2018RHCLivetime

double	pot = 0

double	pot_one = pot_rhc

double	pot_two = pot_fhc

double	livetime = 0

double	livetime_one = livetime_rhc

double	livetime_two = livetime_fhc

double	maxy = 0

double	maxy_one = 1.999

double	maxy_two = 3.999

Sign::Sign_t	wrongs = Sign::kNu

Sign::Sign_t	ws_one = Sign::kNu

Sign::Sign_t	ws_two = Sign::kAntiNu

std::string	extrapName = "extrap"

std::string	horn_one = "rhc"

std::string	horn_two = "fhc"

std::string	period_one = "full"

std::string	period_two = "full"

bool	twobeams = true

int	totquant = 8

bool	DianaIsLazy = true

bool	errorBand = true

Color_t	color_band = kViolet-9

Color_t	color_pred = kViolet+1

Color_t	color_ws = kBlue-7

Color_t	color_nc = kBlue-7

Color_t	color_data = kBlack

Color_t	color_cos = kGreen+1

Color_t	color_bkg = kGray+1

Style_t	line_ws = kSolid

Style_t	line_nc = kSolid

Style_t	line_data = kSolid

Style_t	line_pred = kSolid

Style_t	line_cos = kSolid

Style_t	line_bkg = kSolid

Style_t	marker_ws = kFullCircle

Style_t	marker_nc = kFullSquare

Style_t	marker_data = kFullCircle

Style_t	marker_pred = kFullCircle

Style_t	marker_cos = kFullTriangleUp

Style_t	marker_bkg = kFullCircle

Function Documentation

void HornSettings	(	int	quantile,
		double &	pot,
		double &	livetime,
		Sign::Sign_t &	wrongs,
		double &	maxy
	)

Definition at line 69 of file saveFDMCHists.C.

References livetime_one, livetime_two, maxy_one, maxy_two, pot_one, pot_two, ws_one, and ws_two.

                                                                                                 {
   if(quantile<4){
     pot = pot_one;
     livetime = livetime_one;
     wrongs = ws_one;
     maxy = maxy_one;
   }
   else{
     pot = pot_two;
     livetime = livetime_two;
     wrongs = ws_two;
     maxy = maxy_two;
   }
 
 }// end HornSettings

void saveFDMCHists	(	std::string	anaName = `"fhc"`,
		bool	withExtrap = `true`
	)

Definition at line 87 of file saveFDMCHists.C.

References std::asin(), calc, Close(), om::cout, ana::DefaultOscCalc(), allTimeWatchdog::endl, plotROC::fileName, makeBrightnessMap::fOut, horn_one, ana::Flavors::kAll, ana::Flavors::kAllNuE, ana::Flavors::kAllNuMu, ana::Flavors::kAllNuTau, ana::Sign::kAntiNu, ana::Sign::kBoth, ana::Current::kCC, ana::Current::kNC, ana::Sign::kNu, ana::Flavors::kNuEToNuE, ana::Flavors::kNuEToNuMu, ana::Flavors::kNuEToNuTau, ana::kNuMu, ana::kNumuNoExtrap, ana::Flavors::kNuMuToNuE, ana::Flavors::kNuMuToNuMu, ana::Flavors::kNuMuToNuTau, livetime_fhc, livetime_one, maxy_one, period_one, pot_fhc, pot_one, preliminary, ana::Spectrum::SaveTo(), osc::_IOscCalcAdjustable< T >::SetDmsq32(), osc::_IOscCalcAdjustable< T >::SetTh23(), std::sqrt(), string, totquant, twobeams, and ws_one.

 {
 
   bool preliminary = false;
 
   if(anaName=="rhc"){
     totquant = 4;
     twobeams = false;
   }
   // if fhc, change settings so the first 4 quantiles are fhc instead of rhc
   // dont do anything if rhc as the first 4 quantiles are rhc already
   if(anaName=="fhc"){
     horn_one = "fhc";
     period_one = "full";
     pot_one = pot_fhc;
     livetime_one = livetime_fhc;
     ws_one = Sign::kAntiNu;
     maxy_one = 3.999;
     twobeams = false;
     totquant = 4;
   }
 
   // calculator                                                                                                                                                            
   osc::IOscCalcAdjustable* calc = DefaultOscCalc();
   calc->SetTh23(asin(sqrt(0.558)));
   calc->SetDmsq32(0.00244);
 
   //inputs                                                                                                                             
   std::string dpred_name = "/nova/ana/nu_mu_ana/Ana2018/Predictions/";
 
   std::cout << "\nloading predictions" << std::endl;
   const std::string dpred_name_fhc= dpred_name+"pred_numuconcat_fhc__numu2018.root";
   const std::string dpred_name_rhc= dpred_name+"pred_numuconcat_rhc__numu2018.root";
   std::string dpred_name_one = dpred_name_rhc;
   if(anaName=="fhc") dpred_name_one = dpred_name_fhc;
 
   std::vector<PredictionInterp*> predictions;
   ENumu2018ExtrapType extrap = kNuMu;
   ENumu2018ExtrapType noextrap = kNumuNoExtrap;
   for(int quant = 1; quant < totquant+1; quant++){
     if(withExtrap) predictions.push_back(new PredictionSystNumu2018(extrap, calc, quant, {}, dpred_name_one));
     if(!withExtrap) predictions.push_back(new PredictionSystNumu2018(noextrap, calc, quant, {}, dpred_name_one));
   }
 
 
   std::map< TString, Flavors::Flavors_t > mapFlavor;
   mapFlavor["All"] = Flavors::kAll;
   mapFlavor["AllNuMu"] = Flavors::kAllNuMu;
   mapFlavor["AllNuE"] = Flavors::kAllNuE;
   mapFlavor["AllNuTau"] = Flavors::kAllNuTau;
   mapFlavor["NuEToNuE"] = Flavors::kNuEToNuE;
   mapFlavor["NuEToNuMu"] = Flavors::kNuEToNuMu;
   mapFlavor["NuEToNuTau"] = Flavors::kNuEToNuTau;
   mapFlavor["NuMuToNuE"] = Flavors::kNuMuToNuE;
   mapFlavor["NuMuToNuMu"] = Flavors::kNuMuToNuMu;
   mapFlavor["NuMuToNuTau"] = Flavors::kNuMuToNuTau;
 
   std::map< TString, Sign::Sign_t > mapSign;
   mapSign["Both"] = Sign::kBoth;
   mapSign["Nu"] = Sign::kNu;
   mapSign["AntiNu"] = Sign::kAntiNu;
 
 
   TString fileName = "FDMCHist";
   if(withExtrap) fileName += "_extrap_";
   if(!withExtrap) fileName += "_noExtrap_";
   fileName += anaName;
   TFile *fOut = new TFile(fileName + ".root","RECREATE");
 
   for(int quant = 0; quant < totquant; quant++){
     Spectrum hTot = predictions[quant] -> Predict(calc);
     Spectrum hNC = predictions[quant] -> PredictComponent(calc, Flavors::kAll, Current::kNC, Sign::kBoth);
     hTot.SaveTo(fOut, Form("Q%d_tot",quant+1));
     hNC.SaveTo(fOut, Form("Q%d_NC",quant+1));
     for(auto const& thisSign : mapSign){
       for(auto const& thisFlav : mapFlavor){
         Spectrum hComp = predictions[quant] -> PredictComponent(calc, thisFlav.second, Current::kCC, thisSign.second);
         hComp.SaveTo(fOut, Form("Q%d_",quant+1) + thisSign.first + "_" + thisFlav.first);
       }
     }
   }
 
   fOut -> Close();
   delete fOut;
 
   
 } // End of function