~c-tunnell1/+junk/vlenf_scripts

%!GLoBES

$version="3.0.0"

/* Copyright 2011 Christopher Tunnell

 */

nuflux(#two)<

@flux_file = "flux_file_plus.dat"

@time = 1       /* years */

@stored_muons = 1

@norm =  0.958453635e-10

>

nuflux(#three)<

@flux_file = "flux_file_bkg.flx"

@time = 1       /* years */

@stored_muons = 1

@norm =  1

>

$baseline =             1.0

$target_mass = 	        1.3	/* kt (fiducial mass) */

/* ####### Energy window ####### */

$emin           =       0.0     /* GeV */

$emax           =       4.0     /* GeV */

$binsize        =       {0.25, 0.25,  0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25}

$sampling_min           =       0.0     /* GeV */

$sampling_max           =       4.0     /* GeV */

$sampling_stepsize      =       {0.25, 0.25,  0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25}

/* ######## Baseline setting ######## */

$profiletype = 1

/* ######## Technical information ######## */

$filter_state = 	1

$filter_value = 	0.1

/* ####### Energy resolution ####### */

energy(#Mu_CC_app_eff)<

        @energy = include "mm/MuApp_2cm.txt"

>

energy(#Mu_Bar_CC_disp_eff)<

	@energy = include "mm/MuBarCCBack_2cm.txt"

>

energy(#E_CC_disp_eff)<

	@energy = include "mm/eCCBack_2cm.txt"

>

energy(#NC_disp_eff)<

	@energy = include "mm/MuBarNCBack_2cm.txt"

>

/* ####### Cross sections ####### */

cross(#CC)<

	@cross_file = "XCC.dat"

>

cross(#NC)<

	@cross_file = "XNC.dat"

>

/* ####### Energy dependent efficiencies ###### */

%posteffs={1.0,1.0}

%energ={0.00, 7.0}

%bc=bincenter()

%inter=interpolation(%energ,%posteffs,1,%bc)

/* ####### Channel definitions ####### */

channel(#nu_mu_appearance)<

	@channel = 	#two:	+:	e:	m:	#CC:    #Mu_CC_app_eff

	@post_smearing_efficiencies = copy(%inter)

>

channel(#nu_mu_bar_disappearance)<

	@channel = 	#two:	-:	m:	m:	#CC:    #Mu_Bar_CC_disp_eff

	@post_smearing_efficiencies = copy(%inter)

>

channel(#nu_e_disappearance)<

        @channel =      #two:   +:      e:      e:      #CC:    #E_CC_disp_eff

        @post_smearing_efficiencies = copy(%inter)

>

channel(#nu_e_NC)<

        @channel =      #two:   +:      e:      e:      #NC:    #NC_disp_eff

        @post_smearing_efficiencies = copy(%inter)

>

channel(#nu_mu_bar_NC)<

	@channel = 	#two:	-:	m:	m:	#NC:    #NC_disp_eff

>

/* ####### Setting the rules ####### */

/* ORDER CHANGED BECAUSE OF DEGFINDER! */

rule(#Nu_Mu_Appearance)<

	@signal =	      1.0@#nu_mu_appearance

	@signalerror =	     0.02  :	    0.0001

	@background =	     1.0@#nu_mu_bar_NC : 1.0@#nu_mu_bar_disappearance : 1.0@#nu_e_disappearance : 1.0@#nu_e_NC

	@backgrounderror = 	     0.35  :	    0.0001

	@sys_on_function = 	"chiSpectrumTilt"

	@sys_off_function = 	"chiNoSysSpectrum"

>