RooFit v3.60 -- Developed by Wouter Verkerke and David Kirkby Copyright (C) 2000-2013 NIKHEF, University of California & Stanford University All rights reserved, please read http://roofit.sourceforge.net/license.txt Reading configuration from psychev3.cfg - no POT specified, setting to 3.01E20 - no nubar POT specified, setting to 3.01E20 - no livetime specified, setting to 1417.4 days (only matters for Super-K atmospheric analyses) - found fluxCovMatrix = inputs/flux_covariance_banff_13av2.root - didn't find fluxCovMatrix! Setting default = total_flux_cov - found xsecCovMatrix = inputs/xsec_covariance_2018a_NOMINAL_v11.root - didn't find xsecCovMatrix! Setting default = xsec_cov - found NDdetCovMatrix = inputs/ND280_MultiPi_1p1h_2to8_Cov_WithCorr_29March_convertedToMaCh3.root - didn't find ndDetCovMatrix! Setting default = nddet_cov - didn't find SKdetCovMatrix! Setting default = utils/skdet_total_error_matrix_2015/SKJointErrorMatrix2015_Total_v7_13thMarch2016.root - didn't find skDetCovMatrix! Setting default = SKJointError_Erec_Total - didn't find oscCovMatrix! Setting default = inputs/oscillation_covariance_6par_nondouble.root - didn't find oscCovMatrix! Setting default = osc_cov - didn't find STATONLY, setting default 0 - didn't find DOOSC, setting default 0 - didn't find USESKSYST, setting default 1 - didn't find CUSTOM_REWEIGHT, setting to false --------------------------------- General settings Output file: try_multi_pi_29march_withcovcorr.root Seed: 0 Debug: 0 GPU fit: 1 CPU MP fit: 1 N cores: 8 PlotByMode: 0 --------------------------------- MCMC settings N steps: 0 Flux step scale: 0.05 Xsec step scale: 0.02 ND det step scale:0.04 SK det stet scale:1 Osc step scale: 0.1 Anneal. temp: -1 Random start: 0 Process MCMC: 0 --------------------------------- T2K run settings ND runs: P6_S17 2a 2w 3b 3c 4a 4w 5 6b 6c 6d 6e 7 8a 8w SK nu POT: 3.01e+20 SK antinu POT: 3.01e+20 Livetime: 1417.4 Save nominal: 1 --------------------------------- Samples settings ND280: 1 FHC nu_mu CC0pi: 1 FHC nu_e CC0pi: 1 FHC nu_e CC1pi: 1 RHC nu_mu CC0pi: 1 RHC nu_e CC0pi: 1 SK SF->RFG tuning:1 SK RPA tuning: 1 SK CCCoh tuning: 1 --------------------------------- Samples at ND280 Samples binning ND280 Bin anu mom qe: Bin anu th qe: Bin anu mom nqe: Bin anu th nqe: Bin anu numu p qe: Bin anu numu th qe: Bin anu numu p nqe: Bin anu numu th nqe: Samples binning SK Numu bins Erec: 0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 0.55 0.6 0.65 0.7 0.75 0.8 0.85 0.9 0.95 1 1.05 1.1 1.15 1.2 1.25 1.3 1.35 1.4 1.45 1.5 1.55 1.6 1.65 1.7 1.75 1.8 1.85 1.9 1.95 2 2.05 2.1 2.15 2.2 2.25 2.3 2.35 2.4 2.45 2.5 2.55 2.6 2.65 2.7 2.75 2.8 2.85 2.9 2.95 3 3.25 3.5 3.75 4 4.5 5 5.5 6 7 8 9 10 30 Nue binning option: 0 Nue bins Erec: 0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 0.55 0.6 0.65 0.7 0.75 0.8 0.85 0.9 0.95 1 1.05 1.1 1.15 1.2 1.25 --------------------------------- Input settings Data file: --------------------------------- Real data settings RealDataFitFlag: 0 --------------------------------- Fake data settings FakeDataFitFlag: 1 FakeDataFile: --------------------------------- Asimov settings AsimovFitFlag: 0 NDAsimovFitFlag: 0 MakeAsimovOnly: 0 --------------------------------- Toy fit settings Post. pred. file: EMPTY ToyFitFlag: 0 Toy file: N toys: 20000 Full LLH: 1 --------------------------------- Systematics settings Flux syst: 1 Flux flat: Flux fixed: X-sec syst: 1 X-sec flat: 0 1 2 3 4 X-sec fixed: 12 ND280 detector syst: 1 ND280 detector flat: ND280 detector fixed: Sk detector syst: 1 Sk detector flat: Sk detector fixed: Oscill. flat: Oscill. fixed: Oscill. flat antinu: Oscill. fixed antinu: Reactor const.: 0 Stat only?: 0 Do oscillation?: 0 Use SK syst?: 1 Use BANFF tuning?:0 --------------------------------- Covariance matrices Flux file: inputs/flux_covariance_banff_13av2.root flux matrix: total_flux_cov Xsec file: inputs/xsec_covariance_2018a_NOMINAL_v11.root xsec matrix: xsec_cov BANFF file: BANFF matrix: BANFF vector: ND280 file: inputs/ND280_MultiPi_1p1h_2to8_Cov_WithCorr_29March_convertedToMaCh3.root ND280 matrix: nddet_cov SK file: utils/skdet_total_error_matrix_2015/SKJointErrorMatrix2015_Total_v7_13thMarch2016.root SK matrix: SKJointError_Erec_Total Osc file: inputs/oscillation_covariance_6par_nondouble.root osc matrix: osc_cov Osc file antinu: osc_cov osc matrix: Use BANFF covariance instead of ND data (and flux/xsec cov): 0 --------------------------------- Oscillation parameters Double angle : 0 Binned oscillation : 0 Oscillation binning: Starting parameters: 0.304 0.528 0.0217 7.53e-05 0.002509 -1.601 Starting parameters RHC: Input parameters (if different from starting parameters): Flat likelihood parameters: --------------------------------- Beta options: Use Beta : 0 Flip Beta : 0 Apply Beta to Nue : 0 Apply Diagnal Beta : 0 Succesfully read config psychev3.cfg, now proceeding Trying to load parameters file for parameter: psycheSteering.POT.Production WARNING: production has been overwritten by parameter psycheSteering.POT.Production == 6 WARNING: production has been overwritten by parameter psycheSteering.Bunching.Production == 6 WARNING: production has been overwritten by parameter psycheSteering.Corrections.Production == 6 WARNING: production has been overwritten by parameter psycheSteering.Systematics.Production == 6 AnalysisManager WARNING: The selection do not have to be exclusive to run. AnalysisManager WARNING: This is fine but you should be careful because the order in which you sel().AddSelection() now matters!! Trying to load parameters file for parameter: psycheSelections.Systematics.ApplyFineTuning SelectionBase::Validate(). Validating selection 'kTrackerNumuCC' of type 15numuCCSelection ..... OK SelectionBase::Validate(). Validating selection '' of type 15numuCCSelection ..... OK SelectionBase::AddSplit(). Add split with 3 branches to branch 'trunk' after step 'find_pions' SelectionBase::Validate(). Validating selection 'kTrackerNumuCCMultiPi' of type 22numuCCMultiPiSelection ..... OK SelectionBase::Validate(). Validating selection '' of type 15numuCCSelection ..... OK SelectionBase::Validate(). Validating selection 'kTrackerAntiNumuCC' of type 19antiNumuCCSelection ..... OK SelectionBase::Validate(). Validating selection '' of type 15numuCCSelection ..... OK SelectionBase::Validate(). Validating selection '' of type 19antiNumuCCSelection ..... OK SelectionBase::AddSplit(). Add split with 2 branches to branch 'trunk' after step 'fill summary antinu_CCMultiTrack' SelectionBase::Validate(). Validating selection 'kTrackerAntiNumuCCMultiTrack' of type 29antiNumuCCMultiTrackSelection ..... OK SelectionBase::Validate(). Validating selection '' of type 15numuCCSelection ..... OK SelectionBase::Validate(). Validating selection '' of type 19antiNumuCCSelection ..... OK SelectionBase::Validate(). Validating selection '' of type 15numuCCSelection ..... OK SelectionBase::AddSplit(). Add split with 3 branches to branch 'trunk' after step 'find_pions' SelectionBase::Validate(). Validating selection '' of type 22numuCCMultiPiSelection ..... OK SelectionBase::AddSplit(). Add split with 3 branches to branch 'trunk' after step 'find_pions' SelectionBase::Validate(). Validating selection 'kTrackerAntiNumuCCMultiPi' of type 26antiNumuCCMultiPiSelection ..... OK SelectionBase::Validate(). Validating selection '' of type 15numuCCSelection ..... OK SelectionBase::Validate(). Validating selection 'kTrackerNumuInAntiNuModeCC' of type 30numuBkgInAntiNuModeCCSelection ..... OK SelectionBase::Validate(). Validating selection '' of type 15numuCCSelection ..... OK SelectionBase::Validate(). Validating selection '' of type 30numuBkgInAntiNuModeCCSelection ..... OK SelectionBase::AddSplit(). Add split with 2 branches to branch 'trunk' after step 'fill summary numuBkg_CCMultiTrack' SelectionBase::Validate(). Validating selection 'kTrackerNumuInAntiNuModeCCMultiTrack' of type 40numuBkgInAntiNuModeCCMultiTrackSelection ..... OK SelectionBase::Validate(). Validating selection '' of type 15numuCCSelection ..... OK SelectionBase::Validate(). Validating selection '' of type 15numuCCSelection ..... OK SelectionBase::AddSplit(). Add split with 3 branches to branch 'trunk' after step 'find_pions' SelectionBase::Validate(). Validating selection '' of type 22numuCCMultiPiSelection ..... OK SelectionBase::Validate(). Validating selection '' of type 30numuBkgInAntiNuModeCCSelection ..... OK SelectionBase::AddSplit(). Add split with 3 branches to branch 'trunk' after step 'find_pions' SelectionBase::Validate(). Validating selection 'kTrackerNumuInAntiNuModeCCMultiPi' of type 37numuBkgInAntiNuModeCCMultiPiSelection ..... OK SelectionBase::Validate(). Validating selection '' of type 15numuCCSelection ..... OK SelectionBase::Validate(). Validating selection 'kTrackerNumuCCFGD2' of type 19numuCCFGD2Selection ..... OK SelectionBase::Validate(). Validating selection '' of type 15numuCCSelection ..... OK SelectionBase::Validate(). Validating selection '' of type 15numuCCSelection ..... OK SelectionBase::AddSplit(). Add split with 3 branches to branch 'trunk' after step 'find_pions' SelectionBase::Validate(). Validating selection '' of type 22numuCCMultiPiSelection ..... OK SelectionBase::Validate(). Validating selection '' of type 19numuCCFGD2Selection ..... OK SelectionBase::AddSplit(). Add split with 3 branches to branch 'trunk' after step 'find_pions' SelectionBase::Validate(). Validating selection 'kTrackerNumuCCMultiPiFGD2' of type 26numuCCMultiPiFGD2Selection ..... OK SelectionBase::Validate(). Validating selection '' of type 15numuCCSelection ..... OK SelectionBase::Validate(). Validating selection '' of type 19antiNumuCCSelection ..... OK SelectionBase::Validate(). Validating selection 'kTrackerAntiNumuCCFGD2' of type 23antiNumuCCFGD2Selection ..... OK SelectionBase::Validate(). Validating selection '' of type 15numuCCSelection ..... OK SelectionBase::Validate(). Validating selection '' of type 19antiNumuCCSelection ..... OK SelectionBase::Validate(). Validating selection '' of type 15numuCCSelection ..... OK SelectionBase::Validate(). Validating selection '' of type 19antiNumuCCSelection ..... OK SelectionBase::Validate(). Validating selection '' of type 23antiNumuCCFGD2Selection ..... OK SelectionBase::AddSplit(). Add split with 2 branches to branch 'trunk' after step 'fill summary antinu_CCMultiTrack' SelectionBase::Validate(). Validating selection '' of type 29antiNumuCCMultiTrackSelection ..... OK SelectionBase::AddSplit(). Add split with 2 branches to branch 'trunk' after step 'fill summary AntiNuMuCCMultiTrackFGD2' SelectionBase::Validate(). Validating selection 'kTrackerAntiNumuCCMultiTrackFGD2' of type 33antiNumuCCMultiTrackFGD2Selection ..... OK SelectionBase::Validate(). Validating selection '' of type 15numuCCSelection ..... OK SelectionBase::Validate(). Validating selection '' of type 19antiNumuCCSelection ..... OK SelectionBase::Validate(). Validating selection '' of type 15numuCCSelection ..... OK SelectionBase::AddSplit(). Add split with 3 branches to branch 'trunk' after step 'find_pions' SelectionBase::Validate(). Validating selection '' of type 22numuCCMultiPiSelection ..... OK SelectionBase::Validate(). Validating selection '' of type 15numuCCSelection ..... OK SelectionBase::Validate(). Validating selection '' of type 19antiNumuCCSelection ..... OK SelectionBase::Validate(). Validating selection '' of type 23antiNumuCCFGD2Selection ..... OK SelectionBase::AddSplit(). Add split with 3 branches to branch 'trunk' after step 'find_pions' SelectionBase::Validate(). Validating selection '' of type 26antiNumuCCMultiPiSelection ..... OK SelectionBase::AddSplit(). Add split with 3 branches to branch 'trunk' after step 'find_pions' SelectionBase::Validate(). Validating selection 'kTrackerAntiNumuCCMultiPiFGD2' of type 30antiNumuCCMultiPiFGD2Selection ..... OK SelectionBase::Validate(). Validating selection '' of type 15numuCCSelection ..... OK SelectionBase::Validate(). Validating selection '' of type 15numuCCSelection ..... OK SelectionBase::Validate(). Validating selection '' of type 19numuCCFGD2Selection ..... OK SelectionBase::Validate(). Validating selection '' of type 30numuBkgInAntiNuModeCCSelection ..... OK SelectionBase::Validate(). Validating selection 'kTrackerNumuInAntiNuModeCCFGD2' of type 34numuBkgInAntiNuModeCCFGD2Selection ..... OK SelectionBase::Validate(). Validating selection '' of type 15numuCCSelection ..... OK SelectionBase::Validate(). Validating selection '' of type 15numuCCSelection ..... OK SelectionBase::Validate(). Validating selection '' of type 19numuCCFGD2Selection ..... OK SelectionBase::Validate(). Validating selection '' of type 30numuBkgInAntiNuModeCCSelection ..... OK SelectionBase::Validate(). Validating selection '' of type 15numuCCSelection ..... OK SelectionBase::Validate(). Validating selection '' of type 30numuBkgInAntiNuModeCCSelection ..... OK SelectionBase::Validate(). Validating selection '' of type 34numuBkgInAntiNuModeCCFGD2Selection ..... OK SelectionBase::AddSplit(). Add split with 2 branches to branch 'trunk' after step 'fill summary numuBkg_CCMultiTrack' SelectionBase::Validate(). Validating selection '' of type 40numuBkgInAntiNuModeCCMultiTrackSelection ..... OK SelectionBase::AddSplit(). Add split with 2 branches to branch 'trunk' after step 'fill_summary numuBkgMultiTrackFGD2' SelectionBase::Validate(). Validating selection 'kTrackerNumuInAntiNuModeCCMultiTrackFGD2' of type 44numuBkgInAntiNuModeCCMultiTrackFGD2Selection ..... OK SelectionBase::Validate(). Validating selection '' of type 15numuCCSelection ..... OK SelectionBase::Validate(). Validating selection '' of type 15numuCCSelection ..... OK SelectionBase::AddSplit(). Add split with 3 branches to branch 'trunk' after step 'find_pions' SelectionBase::Validate(). Validating selection '' of type 22numuCCMultiPiSelection ..... OK SelectionBase::Validate(). Validating selection '' of type 15numuCCSelection ..... OK SelectionBase::Validate(). Validating selection '' of type 15numuCCSelection ..... OK SelectionBase::Validate(). Validating selection '' of type 19numuCCFGD2Selection ..... OK SelectionBase::Validate(). Validating selection '' of type 30numuBkgInAntiNuModeCCSelection ..... OK SelectionBase::Validate(). Validating selection '' of type 15numuCCSelection ..... OK SelectionBase::Validate(). Validating selection '' of type 15numuCCSelection ..... OK SelectionBase::AddSplit(). Add split with 3 branches to branch 'trunk' after step 'find_pions' SelectionBase::Validate(). Validating selection '' of type 22numuCCMultiPiSelection ..... OK SelectionBase::Validate(). Validating selection '' of type 19numuCCFGD2Selection ..... OK SelectionBase::AddSplit(). Add split with 3 branches to branch 'trunk' after step 'find_pions' SelectionBase::Validate(). Validating selection '' of type 26numuCCMultiPiFGD2Selection ..... OK SelectionBase::Validate(). Validating selection '' of type 34numuBkgInAntiNuModeCCFGD2Selection ..... OK SelectionBase::Validate(). Validating selection '' of type 30numuBkgInAntiNuModeCCSelection ..... OK SelectionBase::AddSplit(). Add split with 3 branches to branch 'trunk' after step 'find_pions' SelectionBase::Validate(). Validating selection '' of type 37numuBkgInAntiNuModeCCMultiPiSelection ..... OK SelectionBase::AddSplit(). Add split with 3 branches to branch 'trunk' after step 'find_pions' SelectionBase::Validate(). Validating selection 'kTrackerNumuInAntiNuModeCCMultiPiFGD2' of type 41numuBkgInAntiNuModeCCMultiPiFGD2Selection ..... OK SelectionBase::Validate(). Validating selection 'kTrackerNuECC' of type 14nueCCSelection ..... OK SelectionBase::Validate(). Validating selection '' of type 14nueCCSelection ..... OK SelectionBase::Validate(). Validating selection 'kTrackerAntiNuECC' of type 18antiNueCCSelection ..... OK SelectionBase::Validate(). Validating selection '' of type 14nueCCSelection ..... OK SelectionBase::Validate(). Validating selection 'kTrackerGamma' of type 14gammaSelection ..... OK SelectionBase::Validate(). Validating selection '' of type 14nueCCSelection ..... OK SelectionBase::Validate(). Validating selection 'kTrackerNuECCFGD2' of type 18nueCCFGD2Selection ..... OK SelectionBase::Validate(). Validating selection '' of type 14nueCCSelection ..... OK SelectionBase::Validate(). Validating selection '' of type 18nueCCFGD2Selection ..... OK SelectionBase::Validate(). Validating selection 'kTrackerAntiNuECCFGD2' of type 22antiNueCCFGD2Selection ..... OK SelectionBase::Validate(). Validating selection '' of type 14nueCCSelection ..... OK SelectionBase::Validate(). Validating selection '' of type 14gammaSelection ..... OK SelectionBase::Validate(). Validating selection 'kTrackerGammaFGD2' of type 18gammaFGD2Selection ..... OK SelectionManager::GetSelection(). Selection 'kTrackerNumuCC4pi' does not exist !!!! SelectionManager::GetSelection(). Selection 'kTrackerNumuCC4piMultiPi' does not exist !!!! SelectionManager::GetSelection(). Selection 'kTrackerNumuCC4piFGD2' does not exist !!!! SelectionManager::GetSelection(). Selection 'kTrackerNumuCC4piMultiPiFGD2' does not exist !!!! SelectionManager::GetSelection(). Selection 'kTrackerNumuCC4pi' does not exist !!!! SelectionManager::GetSelection(). Selection 'kTrackerNumuCC4piMultiPi' does not exist !!!! SelectionManager::GetSelection(). Selection 'kTrackerNumuCC4piFGD2' does not exist !!!! SelectionManager::GetSelection(). Selection 'kTrackerNumuCC4piMultiPiFGD2' does not exist !!!! Trying to load parameters file for parameter: psycheSystematics.Tracker.UseP0DFV Trying to load parameters file for parameter: psycheND280Utils.ToF.FGD1FGD2 Create FluxWeight in FluxWeightSystematics: ------------------------------------------------------ ------------------------------------------------------ Flux reweight is enabled Folder : $(PSYCHEND280UTILSROOT)/data/tuned13av2 File : $(PSYCHEND280UTILSROOT)/data/tuned13av2/runXXX/ndX_tuned13av2_13anom_runXXX.root Tuning : tuned13a ------------------------------------------------------ ------------------------------------------------------ constructing things! nVarSys = 1 Total parameters: 33 filling links Trying to load parameters file for parameter: psycheSystematics.gdPid.CovFile filling links Trying to load parameters file for parameter: psycheSystematics.luxWeight.CovFile Creating ND280 2018 instance Using TSpline3 reduced 0x1cb79c0 Warning in : no dictionary for class PackageVersion is available Warning in : no dictionary for class DocString is available Warning in : no dictionary for class CorrectionBase is available ****************************** No selections specified in config file, adding all! ****************************** ------------------------------ Adding selections at ND280 Adding selection: FGD1 numuCC 0pi Adding selection: FGD1 numuCC 1pi Adding selection: FGD1 numuCC other Adding selection: FGD2 numuCC 0pi Adding selection: FGD2 numuCC 1pi Adding selection: FGD2 numuCC other Adding selection: FGD1 anti-numuCC 0pi Adding selection: FGD1 anti-numuCC 1pi Adding selection: FGD1 anti-numuCC other Adding selection: FGD2 anti-numuCC 0pi Adding selection: FGD2 anti-numuCC 1pi Adding selection: FGD2 anti-numuCC other Adding selection: FGD1 NuMuBkg CC0pi in AntiNu Mode Adding selection: FGD1 NuMuBkg CC1pi in AntiNu Mode Adding selection: FGD1 NuMuBkg CCother in AntiNu Mode Adding selection: FGD2 NuMuBkg CC0pi in AntiNu Mode Adding selection: FGD2 NuMuBkg CC1pi in AntiNu Mode Adding selection: FGD2 NuMuBkg CCother in AntiNu Mode - Using ND280 GPU version ------------------------------ Created covariance matrix named: xsec_cov from file: inputs/xsec_covariance_2018a_NOMINAL_v11.root Scanned cross-section matrix, found xsecmodel = NIWG_2018a ################################################# Printing covarianceXsec2015: NIWG model: NIWG_2018a Number of parameters: 33 Global parameter map: # | Name | Nom. | Prior | Error | Lower | Upper 0 | MAQE | 1 | 1 | +/- 0.025 | 0 | 1e+04 1 | pF_C | 1 | 1 | +/- 0.0599 | 0.922 | 1.27 2 | pF_O | 1 | 1 | +/- 0.0578 | 0.889 | 1.22 3 | 2p2h_norm_nu | 1 | 1 | +/- 1 | 0 | 1e+04 4 | 2p2h_norm_nubar | 1 | 1 | +/- 1 | 0 | 1e+04 5 | 2p2h_normCtoO | 1 | 1 | +/- 0.2 | 0 | 1e+04 6 | 2p2h_shape_C | 0 | 0 | +/- 3 | -1 | 1 7 | 2p2h_shape_O | 0 | 0 | +/- 3 | -1 | 1 8 | BeRPA_A | 0.59 | 0.59 | +/- 0.118 | 0 | 1e+04 9 | BeRPA_B | 1.05 | 1.05 | +/- 0.21 | 0 | 1e+04 10 | BeRPA_D | 1.13 | 1.13 | +/- 0.17 | 0 | 1e+04 11 | BeRPA_E | 0.88 | 0.88 | +/- 0.352 | 0 | 1e+04 12 | BeRPA_U | 1.2 | 1.2 | +/- 0.1 | 0 | 1e+04 13 | CA5 | 1 | 0.95 | +/- 0.149 | 0 | 1e+04 14 | MARES | 1 | 1.13 | +/- 0.158 | 0 | 1e+04 15 | ISO_BKG | 1 | 0.738 | +/- 0.308 | 0 | 1e+04 16 | CC_norm_nu | 1 | 1 | +/- 0.02 | 0 | 1e+04 17 | CC_norm_nubar | 1 | 1 | +/- 0.01 | 0 | 1e+04 18 | nue_numu | 1 | 1 | +/- 0.0283 | 0 | 1e+04 19 | nuebar_numubar | 1 | 1 | +/- 0.0283 | 0 | 1e+04 20 | CC_DIS | 0 | 0 | +/- 0.4 | -1e+04 | 1e+04 21 | CC_Coh_C | 1 | 1 | +/- 0.3 | 0 | 1e+04 22 | CC_Coh_O | 1 | 1 | +/- 0.3 | 0 | 1e+04 23 | NC_Coh | 1 | 1 | +/- 0.3 | 0 | 1e+04 24 | NC_1gamma | 1 | 1 | +/- 1 | 0 | 1e+04 25 | NC_other_near | 1 | 1 | +/- 0.3 | 0 | 1e+04 26 | NC_other_far | 1 | 1 | +/- 0.3 | 0 | 1e+04 27 | FEFQE | 1 | 1.07 | +/- 0.313 | -1e+04 | 1e+04 28 | FEFQEH | 1 | 1.01 | +/- 0.477 | -1e+04 | 1e+04 29 | FEFINEL | 1 | 1 | +/- 1.1 | -1e+04 | 1e+04 30 | FEFABS | 1 | 1.28 | +/- 0.393 | -1e+04 | 1e+04 31 | FEFCX | 1 | 0.697 | +/- 0.305 | -1e+04 | 1e+04 32 | EB_dial | 0 | 0 | +/- 2 | -1e+04 | 1e+04 Normalisation parameters:12 # | Global # | Name | Int. mode | Target | Type 0 | 3 | 2p2h_norm_nu | 9 | 12 16 | 12 14 1 | 4 | 2p2h_norm_nubar | 9 | 12 16 | -12 -14 2 | 5 | 2p2h_normCtoO | 9 | 16 | 12 14 -12 -14 3 | 16 | CC_norm_nu | -3 | 12 16 | 12 14 4 | 17 | CC_norm_nubar | -3 | 12 16 | -12 -14 5 | 18 | nue_numu | -1 | 12 16 | 12 6 | 19 | nuebar_numubar | -1 | 12 16 | -12 7 | 21 | CC_Coh_C | 2 | 12 | 12 14 -12 -14 8 | 22 | CC_Coh_O | 2 | 16 | 12 14 -12 -14 9 | 23 | NC_Coh | 7 | | 12 14 -12 -14 10 | 24 | NC_1gamma | 10 | 12 16 | 12 14 -12 -14 11 | 25 | NC_other_near | 8 | | 12 14 -12 -14 SK Only normalisation parameters:1 # | Global # | Name | Int. mode | Target | Type 0 | 26 | NC_other_far | 8 | 16 | 12 SK Only spline parameter: 1 # | Global # | Name 0 | 32 | EB_dial Spline parameters: 14 # | Global # | Name 0 | 0 | MAQE 1 | 1 | pF_C 2 | 2 | pF_O 3 | 6 | 2p2h_shape_C 4 | 7 | 2p2h_shape_O 5 | 13 | CA5 6 | 14 | MARES 7 | 15 | ISO_BKG 8 | 20 | CC_DIS 9 | 27 | FEFQE 10 | 28 | FEFQEH 11 | 29 | FEFINEL 12 | 30 | FEFABS 13 | 31 | FEFCX Unique spline parameters: 14 # | Global # | Name 0 | 0 | MAQE 1 | 1 | pF_C 2 | 2 | pF_O 3 | 6 | 2p2h_shape_C 4 | 7 | 2p2h_shape_O 5 | 13 | CA5 6 | 14 | MARES 7 | 15 | ISO_BKG 8 | 20 | CC_DIS 9 | 27 | FEFQE 10 | 28 | FEFQEH 11 | 29 | FEFINEL 12 | 30 | FEFABS 13 | 31 | FEFCX What spline parameter uses what spline? Spline # Uses unique # | Param Name Uses Name 0 0 | MAQE MAQE 1 1 | pF_C pF_C 2 2 | pF_O pF_O 3 3 | 2p2h_shape_C 2p2h_shape_C 4 4 | 2p2h_shape_O 2p2h_shape_O 5 5 | CA5 CA5 6 6 | MARES MARES 7 7 | ISO_BKG ISO_BKG 8 8 | CC_DIS CC_DIS 9 9 | FEFQE FEFQE 10 10 | FEFQEH FEFQEH 11 11 | FEFINEL FEFINEL 12 12 | FEFABS FEFABS 13 13 | FEFCX FEFCX Functional parameters: 5 # | Global # | Name 0 | 8 | BeRPA_A 1 | 9 | BeRPA_B 2 | 10 | BeRPA_D 3 | 11 | BeRPA_E 4 | 12 | BeRPA_U Done printing covarianceXsec2015 ################################################# ------------------------------ Created covariance matrix named: total_flux_cov from file: inputs/flux_covariance_banff_13av2.root Flux covariance: Near detector size: 8 ND280_FHC_numu 0 - 0.4 - 0.5 - 0.6 - 0.7 - 1 - 1.5 - 2.5 - 3.5 - 5 - 7 - 30 ND280_FHC_numubar 0 - 0.7 - 1 - 1.5 - 2.5 - 30 ND280_FHC_nue 0 - 0.5 - 0.7 - 0.8 - 1.5 - 2.5 - 4 - 30 ND280_FHC_nuebar 0 - 2.5 - 30 ND280_RHC_numu 0 - 0.7 - 1 - 1.5 - 2.5 - 30 ND280_RHC_numubar 0 - 0.4 - 0.5 - 0.6 - 0.7 - 1 - 1.5 - 2.5 - 3.5 - 5 - 7 - 30 ND280_RHC_nue 0 - 2.5 - 30 ND280_RHC_nuebar 0 - 0.5 - 0.7 - 0.8 - 1.5 - 2.5 - 4 - 30 Far detector size: 8 SK_FHC_numu 0 - 0.4 - 0.5 - 0.6 - 0.7 - 1 - 1.5 - 2.5 - 3.5 - 5 - 7 - 30 SK_FHC_numubar 0 - 0.7 - 1 - 1.5 - 2.5 - 30 SK_FHC_nue 0 - 0.5 - 0.7 - 0.8 - 1.5 - 2.5 - 4 - 30 SK_FHC_nuebar 0 - 2.5 - 30 SK_RHC_numu 0 - 0.7 - 1 - 1.5 - 2.5 - 30 SK_RHC_numubar 0 - 0.4 - 0.5 - 0.6 - 0.7 - 1 - 1.5 - 2.5 - 3.5 - 5 - 7 - 30 SK_RHC_nue 0 - 2.5 - 30 SK_RHC_nuebar 0 - 0.5 - 0.7 - 0.8 - 1.5 - 2.5 - 4 - 30 ------------------------------ Created covariance matrix named: nddet_cov from file: inputs/ND280_MultiPi_1p1h_2to8_Cov_WithCorr_29March_convertedToMaCh3.root NDDET SIMPLE: 496 ------------------------------ ------------------------------ ------------------------------ Setting MAQE (parameter 0) to flat prior Setting pF_C (parameter 1) to flat prior Setting pF_O (parameter 2) to flat prior Setting 2p2h_norm_nu (parameter 3) to flat prior Setting 2p2h_norm_nubar (parameter 4) to flat prior Setting BeRPA_U (parameter 12) to fixed at 1.2 ------------------------------ Normalisation parameters: 0 |3 |2p2h_norm_nu |9 1 |4 |2p2h_norm_nubar |9 2 |5 |2p2h_normCtoO |9 3 |16 |CC_norm_nu |-3 4 |17 |CC_norm_nubar |-3 5 |18 |nue_numu |-1 6 |19 |nuebar_numubar |-1 7 |21 |CC_Coh_C |2 8 |22 |CC_Coh_O |2 9 |23 |NC_Coh |7 10 |24 |NC_1gamma |10 11 |25 |NC_other_near |8 Total spline parameters: 0 |0 |MAQE 1 |1 |pF_C 2 |2 |pF_O 3 |6 |2p2h_shape_C 4 |7 |2p2h_shape_O 5 |13 |CA5 6 |14 |MARES 7 |15 |ISO_BKG 8 |20 |CC_DIS 9 |27 |FEFQE 10 |28 |FEFQEH 11 |29 |FEFINEL 12 |30 |FEFABS 13 |31 |FEFCX Unique spline parameters: 0 |0 |MAQE 1 |1 |pF_C 2 |2 |pF_O 3 |6 |2p2h_shape_C 4 |7 |2p2h_shape_O 5 |13 |CA5 6 |14 |MARES 7 |15 |ISO_BKG 8 |20 |CC_DIS 9 |27 |FEFQE 10 |28 |FEFQEH 11 |29 |FEFINEL 12 |30 |FEFABS 13 |31 |FEFCX Shared spline parameters: Finished setXsecCov in samplePDFND2018, xsecmodel = NIWG_2018a ------------------------------ xsec_cov setStepScale() = 0.02 total_flux_cov setStepScale() = 0.05 nddet_cov setStepScale() = 0.04 ------------------------------ 121885 data events passed selection in psyche 0 data events didn't pass sample enabled test Looped 121885 events and 121885 in tree Sample Binned Unbinned FGD1 numuCC 0pi 33553 33615 FGD1 numuCC 1pi 7757 7774 FGD1 numuCC other 8068 8139 FGD2 numuCC 0pi 33462 33543 FGD2 numuCC 1pi 6133 6159 FGD2 numuCC other 7664 7734 FGD1 anti-numuCC 0pi 6368 6380 FGD1 anti-numuCC 1pi 535 536 FGD1 anti-numuCC other 1102 1110 FGD2 anti-numuCC 0pi 6451 6464 FGD2 anti-numuCC 1pi 465 465 FGD2 anti-numuCC other 1032 1046 FGD1 NuMuBkg CC0pi in AntiNu Mode 2707 2728 FGD1 NuMuBkg CC1pi in AntiNu Mode 847 852 FGD1 NuMuBkg CCother in AntiNu Mode1015 1030 FGD2 NuMuBkg CC0pi in AntiNu Mode 2648 2667 FGD2 NuMuBkg CC1pi in AntiNu Mode 693 697 FGD2 NuMuBkg CCother in AntiNu Mode932 946 Total: 121432 121885 Data fill took 27.3634s ------------------------------ Making TChain from: /vols/t2k/users/cvw09/software/T2KReWeight_OA2018/app/CCMultiPi_MAQEfix/MCRun2a_NewMAQE.root Making TChain from: /vols/t2k/users/cvw09/software/T2KReWeight_OA2018/app/CCMultiPi_MAQEfix/SANDRun2a_NewMAQE.root Making TChain from: /vols/t2k/users/cvw09/software/T2KReWeight_OA2018/app/CCMultiPi_MAQEfix/MCRun2w_NewMAQE.root Making TChain from: /vols/t2k/users/cvw09/software/T2KReWeight_OA2018/app/CCMultiPi_MAQEfix/SANDRun2w_NewMAQE.root Making TChain from: /vols/t2k/users/cvw09/software/T2KReWeight_OA2018/app/CCMultiPi_MAQEfix/MCRun3b_NewMAQE.root Making TChain from: /vols/t2k/users/cvw09/software/T2KReWeight_OA2018/app/CCMultiPi_MAQEfix/SANDRun3b_NewMAQE.root Making TChain from: /vols/t2k/users/cvw09/software/T2KReWeight_OA2018/app/CCMultiPi_MAQEfix/MCRun3c_NewMAQE.root Making TChain from: /vols/t2k/users/cvw09/software/T2KReWeight_OA2018/app/CCMultiPi_MAQEfix/SANDRun3c_NewMAQE.root Making TChain from: /vols/t2k/users/cvw09/software/T2KReWeight_OA2018/app/CCMultiPi_MAQEfix/MCRun4a_NewMAQE.root Making TChain from: /vols/t2k/users/cvw09/software/T2KReWeight_OA2018/app/CCMultiPi_MAQEfix/SANDRun4a_NewMAQE.root Making TChain from: /vols/t2k/users/cvw09/software/T2KReWeight_OA2018/app/CCMultiPi_MAQEfix/MCRun4w_NewMAQE.root Making TChain from: /vols/t2k/users/cvw09/software/T2KReWeight_OA2018/app/CCMultiPi_MAQEfix/SANDRun4w_NewMAQE.root Making TChain from: /vols/t2k/users/cvw09/software/T2KReWeight_OA2018/app/CCMultiPi_MAQEfix/MCRun5_NewMAQE.root Making TChain from: /vols/t2k/users/cvw09/software/T2KReWeight_OA2018/app/CCMultiPi_MAQEfix/SANDRun5_NewMAQE.root Making TChain from: /vols/t2k/users/cvw09/software/T2KReWeight_OA2018/app/CCMultiPi_MAQEfix/MCRun6b_NewMAQE.root Making TChain from: /vols/t2k/users/cvw09/software/T2KReWeight_OA2018/app/CCMultiPi_MAQEfix/SANDRun6b_NewMAQE.root Making TChain from: /vols/t2k/users/cvw09/software/T2KReWeight_OA2018/app/CCMultiPi_MAQEfix/MCRun6c_NewMAQE.root Making TChain from: /vols/t2k/users/cvw09/software/T2KReWeight_OA2018/app/CCMultiPi_MAQEfix/SANDRun6c_NewMAQE.root Making TChain from: /vols/t2k/users/cvw09/software/T2KReWeight_OA2018/app/CCMultiPi_MAQEfix/MCRun6d_NewMAQE.root Making TChain from: /vols/t2k/users/cvw09/software/T2KReWeight_OA2018/app/CCMultiPi_MAQEfix/SANDRun6d_NewMAQE.root Making TChain from: /vols/t2k/users/cvw09/software/T2KReWeight_OA2018/app/CCMultiPi_MAQEfix/MCRun6e_NewMAQE.root Making TChain from: /vols/t2k/users/cvw09/software/T2KReWeight_OA2018/app/CCMultiPi_MAQEfix/SANDRun6e_NewMAQE.root Making TChain from: /vols/t2k/users/cvw09/software/T2KReWeight_OA2018/app/CCMultiPi_MAQEfix/MCRun7_NewMAQE.root Making TChain from: /vols/t2k/users/cvw09/software/T2KReWeight_OA2018/app/CCMultiPi_MAQEfix/SANDRun7_NewMAQE.root Making TChain from: /vols/t2k/users/cvw09/software/T2KReWeight_OA2018/app/CCMultiPi_MAQEfix/MCRun8a_NewMAQE.root Making TChain from: /vols/t2k/users/cvw09/software/T2KReWeight_OA2018/app/CCMultiPi_MAQEfix/SANDRun8a_NewMAQE.root Making TChain from: /vols/t2k/users/cvw09/software/T2KReWeight_OA2018/app/CCMultiPi_MAQEfix/MCRun8w_NewMAQE.root Making TChain from: /vols/t2k/users/cvw09/software/T2KReWeight_OA2018/app/CCMultiPi_MAQEfix/SANDRun8w_NewMAQE.root Total entries in TChain: 1814809 Looping over MC events... Event 0/1814809 (0%) Running on file : /vols/t2k/users/cvw09/software/T2KReWeight_OA2018/app/CCMultiPi_MAQEfix/MCRun2a_NewMAQE.root Running on file : /vols/t2k/users/cvw09/software/T2KReWeight_OA2018/app/CCMultiPi_MAQEfix/SANDRun2a_NewMAQE.root Running on file : /vols/t2k/users/cvw09/software/T2KReWeight_OA2018/app/CCMultiPi_MAQEfix/MCRun2w_NewMAQE.root Event 90741/1814809 (5%) ***WARNING, detector weight from psyche in input file = 68.8171 --- Capping at 10! Selection: FGD1 numuCC 1pi Momentum: 16884.5 CosTheta: 0.997471 POTWeight: 0.0360446 PileUpWeight: 1 DetWeight: 1 DetWeight_a: 68.8171 FluxWeight: 1 FluxWeight_a: 0.677369 Running on file : /vols/t2k/users/cvw09/software/T2KReWeight_OA2018/app/CCMultiPi_MAQEfix/SANDRun2w_NewMAQE.root Running on file : /vols/t2k/users/cvw09/software/T2KReWeight_OA2018/app/CCMultiPi_MAQEfix/MCRun3b_NewMAQE.root Event 181482/1814809 (10%) Running on file : /vols/t2k/users/cvw09/software/T2KReWeight_OA2018/app/CCMultiPi_MAQEfix/SANDRun3b_NewMAQE.root Running on file : /vols/t2k/users/cvw09/software/T2KReWeight_OA2018/app/CCMultiPi_MAQEfix/MCRun3c_NewMAQE.root Event 272223/1814809 (15%) Event 362964/1814809 (20%) ***WARNING, detector weight from psyche in input file = -0.417013 --- Setting to 0! Selection: FGD1 numuCC other Momentum: 17823.8 CosTheta: 0.994803 POTWeight: 0.0518178 PileUpWeight: 1 DetWeight: 1 DetWeight_a: -0.417013 FluxWeight: 1 FluxWeight_a: 1.25482 Running on file : /vols/t2k/users/cvw09/software/T2KReWeight_OA2018/app/CCMultiPi_MAQEfix/SANDRun3c_NewMAQE.root Running on file : /vols/t2k/users/cvw09/software/T2KReWeight_OA2018/app/CCMultiPi_MAQEfix/MCRun4a_NewMAQE.root Event 453705/1814809 (25%) Event 544446/1814809 (30%) ***WARNING, detector weight from psyche in input file = -0.0137773 --- Setting to 0! Selection: FGD1 numuCC other Momentum: 814.214 CosTheta: 0.971249 POTWeight: 0.0509402 PileUpWeight: 1 DetWeight: 1 DetWeight_a: -0.0137773 FluxWeight: 1 FluxWeight_a: 1.21461 Event 635187/1814809 (35%) Event 725928/1814809 (40%) Running on file : /vols/t2k/users/cvw09/software/T2KReWeight_OA2018/app/CCMultiPi_MAQEfix/SANDRun4a_NewMAQE.root Running on file : /vols/t2k/users/cvw09/software/T2KReWeight_OA2018/app/CCMultiPi_MAQEfix/MCRun4w_NewMAQE.root ***WARNING, detector weight from psyche in input file = -0.728417 --- Setting to 0! Selection: FGD1 numuCC 0pi Momentum: 11.6877 CosTheta: 0.657659 POTWeight: 0.0726197 PileUpWeight: 1 DetWeight: 1 DetWeight_a: -0.728417 FluxWeight: 1 FluxWeight_a: 1.11217 Event 816669/1814809 (45%) ***WARNING, detector weight from psyche in input file = 72.247 --- Capping at 10! Selection: FGD1 numuCC other Momentum: 12594.5 CosTheta: 0.99788 POTWeight: 0.0726197 PileUpWeight: 1 DetWeight: 1 DetWeight_a: 72.247 FluxWeight: 1 FluxWeight_a: 0.923416 ***WARNING, detector weight from psyche in input file = -0.282073 --- Setting to 0! Selection: FGD1 numuCC other Momentum: 6070.47 CosTheta: 0.996679 POTWeight: 0.0726197 PileUpWeight: 1 DetWeight: 1 DetWeight_a: -0.282073 FluxWeight: 1 FluxWeight_a: 1.12631 Event 907410/1814809 (50%) Running on file : /vols/t2k/users/cvw09/software/T2KReWeight_OA2018/app/CCMultiPi_MAQEfix/SANDRun4w_NewMAQE.root Running on file : /vols/t2k/users/cvw09/software/T2KReWeight_OA2018/app/CCMultiPi_MAQEfix/MCRun5_NewMAQE.root ***WARNING, detector weight from psyche in input file = -0.746103 --- Setting to 0! Selection: FGD2 NuMuBkg CC0pi in AntiNu Mode Momentum: 2.84668 CosTheta: 0.305819 POTWeight: 0.0189299 PileUpWeight: 1 DetWeight: 1 DetWeight_a: -0.746103 FluxWeight: 1 FluxWeight_a: 1.04712 ***WARNING, detector weight from psyche in input file = -0.46981 --- Setting to 0! Selection: FGD1 NuMuBkg CC0pi in AntiNu Mode Momentum: 27823.5 CosTheta: 0.949169 POTWeight: 0.0189299 PileUpWeight: 1 DetWeight: 1 DetWeight_a: -0.46981 FluxWeight: 1 FluxWeight_a: 1.04476 Event 998151/1814809 (55%) ***WARNING, detector weight from psyche in input file = -0.0428515 --- Setting to 0! Selection: FGD1 anti-numuCC 0pi Momentum: 3242.28 CosTheta: 0.727517 POTWeight: 0.0189299 PileUpWeight: 1 DetWeight: 1 DetWeight_a: -0.0428515 FluxWeight: 1 FluxWeight_a: 0.969283 Running on file : /vols/t2k/users/cvw09/software/T2KReWeight_OA2018/app/CCMultiPi_MAQEfix/SANDRun5_NewMAQE.root Running on file : /vols/t2k/users/cvw09/software/T2KReWeight_OA2018/app/CCMultiPi_MAQEfix/MCRun6b_NewMAQE.root ***WARNING, detector weight from psyche in input file = 15.4898 --- Capping at 10! Selection: FGD1 anti-numuCC other Momentum: 20308.6 CosTheta: 0.997553 POTWeight: 0.0898129 PileUpWeight: 1 DetWeight: 1 DetWeight_a: 15.4898 FluxWeight: 1 FluxWeight_a: 1.76629 Running on file : /vols/t2k/users/cvw09/software/T2KReWeight_OA2018/app/CCMultiPi_MAQEfix/SANDRun6b_NewMAQE.root Running on file : /vols/t2k/users/cvw09/software/T2KReWeight_OA2018/app/CCMultiPi_MAQEfix/MCRun6c_NewMAQE.root Running on file : /vols/t2k/users/cvw09/software/T2KReWeight_OA2018/app/CCMultiPi_MAQEfix/SANDRun6c_NewMAQE.root Running on file : /vols/t2k/users/cvw09/software/T2KReWeight_OA2018/app/CCMultiPi_MAQEfix/MCRun6d_NewMAQE.root Event 1088892/1814809 (60%) ***WARNING, detector weight from psyche in input file = 39.4941 --- Capping at 10! Selection: FGD1 anti-numuCC 1pi Momentum: 55224.2 CosTheta: 0.999818 POTWeight: 0.11264 PileUpWeight: 1 DetWeight: 1 DetWeight_a: 39.4941 FluxWeight: 1 FluxWeight_a: 0.619521 Running on file : /vols/t2k/users/cvw09/software/T2KReWeight_OA2018/app/CCMultiPi_MAQEfix/SANDRun6d_NewMAQE.root Running on file : /vols/t2k/users/cvw09/software/T2KReWeight_OA2018/app/CCMultiPi_MAQEfix/MCRun6e_NewMAQE.root ***WARNING, detector weight from psyche in input file = -0.379341 --- Setting to 0! Selection: FGD1 NuMuBkg CCother in AntiNu Mode Momentum: 589.834 CosTheta: 0.329399 POTWeight: 0.0990679 PileUpWeight: 1 DetWeight: 1 DetWeight_a: -0.379341 FluxWeight: 1 FluxWeight_a: 1.0337 Running on file : /vols/t2k/users/cvw09/software/T2KReWeight_OA2018/app/CCMultiPi_MAQEfix/SANDRun6e_NewMAQE.root Running on file : /vols/t2k/users/cvw09/software/T2KReWeight_OA2018/app/CCMultiPi_MAQEfix/MCRun7_NewMAQE.root Event 1179633/1814809 (65%) ***WARNING, detector weight from psyche in input file = -0.295735 --- Setting to 0! Selection: FGD1 anti-numuCC 0pi Momentum: 125333 CosTheta: 0.990015 POTWeight: 0.0722969 PileUpWeight: 1 DetWeight: 1 DetWeight_a: -0.295735 FluxWeight: 1 FluxWeight_a: 0.9793 ***WARNING, detector weight from psyche in input file = -0.0878035 --- Setting to 0! Selection: FGD1 NuMuBkg CC0pi in AntiNu Mode Momentum: 12860.6 CosTheta: 0.771413 POTWeight: 0.0722969 PileUpWeight: 1 DetWeight: 1 DetWeight_a: -0.0878035 FluxWeight: 1 FluxWeight_a: 1.04862 Event 1270374/1814809 (70%) ***WARNING, detector weight from psyche in input file = -0.37712 --- Setting to 0! Selection: FGD1 NuMuBkg CCother in AntiNu Mode Momentum: 16205.1 CosTheta: 0.747908 POTWeight: 0.0722969 PileUpWeight: 1 DetWeight: 1 DetWeight_a: -0.37712 FluxWeight: 1 FluxWeight_a: 1.17336 ***WARNING, detector weight from psyche in input file = 10.6339 --- Capping at 10! Selection: FGD1 NuMuBkg CCother in AntiNu Mode Momentum: 11921 CosTheta: 0.999862 POTWeight: 0.0722969 PileUpWeight: 1 DetWeight: 1 DetWeight_a: 10.6339 FluxWeight: 1 FluxWeight_a: 0.856219 Running on file : /vols/t2k/users/cvw09/software/T2KReWeight_OA2018/app/CCMultiPi_MAQEfix/SANDRun7_NewMAQE.root Running on file : /vols/t2k/users/cvw09/software/T2KReWeight_OA2018/app/CCMultiPi_MAQEfix/MCRun8a_NewMAQE.root ***WARNING, detector weight from psyche in input file = -0.649987 --- Setting to 0! Selection: FGD1 numuCC other Momentum: 7340.66 CosTheta: 0.988366 POTWeight: 0.114283 PileUpWeight: 1 DetWeight: 1 DetWeight_a: -0.649987 FluxWeight: 1 FluxWeight_a: 1.24572 ***WARNING, detector weight from psyche in input file = 28.6031 --- Capping at 10! Selection: FGD1 numuCC other Momentum: 9135.96 CosTheta: 0.99871 POTWeight: 0.114283 PileUpWeight: 1 DetWeight: 1 DetWeight_a: 28.6031 FluxWeight: 1 FluxWeight_a: 1.67076 Event 1361115/1814809 (75%) Event 1451856/1814809 (80%) ***WARNING, detector weight from psyche in input file = 42.0261 --- Capping at 10! Selection: FGD1 numuCC other Momentum: 9420.76 CosTheta: 0.988459 POTWeight: 0.114283 PileUpWeight: 1 DetWeight: 1 DetWeight_a: 42.0261 FluxWeight: 1 FluxWeight_a: 1.30071 Event 1542597/1814809 (85%) Running on file : /vols/t2k/users/cvw09/software/T2KReWeight_OA2018/app/CCMultiPi_MAQEfix/SANDRun8a_NewMAQE.root Running on file : /vols/t2k/users/cvw09/software/T2KReWeight_OA2018/app/CCMultiPi_MAQEfix/MCRun8w_NewMAQE.root Event 1633338/1814809 (90%) Event 1724079/1814809 (95%) Running on file : /vols/t2k/users/cvw09/software/T2KReWeight_OA2018/app/CCMultiPi_MAQEfix/SANDRun8w_NewMAQE.root ------------------------------------------------------ Finished looping psyche events (took 12548.3s) There are 1695889 which passed cuts Allocated 67.8356 MB to the detector structs 67.8356 MB to the xsec structs (excluding new TSpline3 or TF1) Found 12 events with NEGATIVE detector weight from psyche! Found 7 events with >10 detector weight from psyche! Finished cross-section loading: 67.8356 MB + splines allocated Cross-section loading took 12548.3s MC fill took 12548.3s ------------------------------ -- GPUING WITH {X} and {Y,B,C,D} arrays and master spline containing TSpline3_red Found 1695889 events Found 7 knots at max Found 14 parameters Number of splines = 4784620 --- INITIALISED ND280 2018 {X}, {YBCD} ARRAYS --- 1695889 events with 4784620 splines On average 2.82131 splines per event (4784620/1695889) Size of x array = 133.969 MB Size of coefficient {y,b,c,d} array = 535.877 MB Size of parameter # array = 19.1385 MB Total size = 688.985 MB memory on CPU to move to GPU GPU weight array (GPU->CPU every step) = 19.1385 MB Parameter value array (CPU->GPU every step) = 5.6e-05 MB Allocated 4784620 entries for paramNo arrays, size = 19.138480 MB Allocated 33492340 entries for x coeff arrays, size = 133.969360 MB Allocated 133969360 entries for {ybcd} coeff arrays, size = 535.877440 MB Good GPU loading GPU spline fill took 11.6378s ------------------------------ Setup complete! Sample Data MC -LLH | FGD1 numuCC 0pi 33553 31529.2 222.492 | FGD1 numuCC 1pi 7757 7997.52 125.527 | FGD1 numuCC other 8068 6793.45 251.601 | FGD2 numuCC 0pi 33462 31732.1 205.517 | FGD2 numuCC 1pi 6133 6418.42 111.47 | FGD2 numuCC other 7664 6562.61 237.67 | FGD1 anti-numuCC 0pi 6368 6370.47 57.7938 | FGD1 anti-numuCC 1pi 535 533.342 21.2839 | FGD1 anti-numuCC other 1102 1023.3 22.334 | FGD2 anti-numuCC 0pi 6451 6283.55 63.1408 | FGD2 anti-numuCC 1pi 465 483.358 15.8615 | FGD2 anti-numuCC other 1032 944.017 26.1267 | FGD1 NuMuBkg CC0pi in AntiNu Mode 2707 2500.77 52.7746 | FGD1 NuMuBkg CC1pi in AntiNu Mode 847 845.922 24.0801 | FGD1 NuMuBkg CCother in AntiNu Mode1015 800.914 51.8522 | FGD2 NuMuBkg CC0pi in AntiNu Mode 2648 2554.46 56.7928 | FGD2 NuMuBkg CC1pi in AntiNu Mode 693 670.813 31.2077 | FGD2 NuMuBkg CCother in AntiNu Mode932 785.661 37.4582 | Total 121432 114830 1614.98 | Single reweight of MC took 0.0915s to complete ------------------------------------------ Have set Asimov fake-data to reweighted MC With cross-section settings: Printing parameters for xsec_cov Name Nominal Init Current Proposed MAQE 1 1 1 1 pF_C 1 1 1 1 pF_O 1 1 1 1 2p2h_norm_nu 1 1 1 1 2p2h_norm_nubar1 1 1 1 2p2h_normCtoO 1 1 1 1 2p2h_shape_C 0 0 0 0 2p2h_shape_O 0 0 0 0 BeRPA_A 0.59 0.59 0.59 0.59 BeRPA_B 1.05 1.05 1.05 1.05 BeRPA_D 1.13 1.13 1.13 1.13 BeRPA_E 0.88 0.88 0.88 0.88 BeRPA_U 1.2 1.2 1.2 1.2 CA5 0.950495 0.950495 0.950495 0.950495 MARES 1.12632 1.12632 1.12632 1.12632 ISO_BKG 0.738462 0.738462 0.738462 0.738462 CC_norm_nu 1 1 1 1 CC_norm_nubar 1 1 1 1 nue_numu 1 1 1 1 nuebar_numubar 1 1 1 1 CC_DIS 0 0 0 0 CC_Coh_C 1 1 1 1 CC_Coh_O 1 1 1 1 NC_Coh 1 1 1 1 NC_1gamma 1 1 1 1 NC_other_near 1 1 1 1 NC_other_far 1 1 1 1 FEFQE 1.069 1.069 1.069 1.069 FEFQEH 1.01333 1.01333 1.01333 1.01333 FEFINEL 1.002 1.002 1.002 1.002 FEFABS 1.27636 1.27636 1.27636 1.27636 FEFCX 0.697 0.697 0.697 0.697 EB_dial 0 0 0 0 ------------------------------------------ Sample Data MC -LLH | FGD1 numuCC 0pi 31529.2 31529.2 0 | FGD1 numuCC 1pi 7997.52 7997.52 0 | FGD1 numuCC other 6793.45 6793.45 0 | FGD2 numuCC 0pi 31732.1 31732.1 0 | FGD2 numuCC 1pi 6418.42 6418.42 0 | FGD2 numuCC other 6562.61 6562.61 0 | FGD1 anti-numuCC 0pi 6370.47 6370.47 0 | FGD1 anti-numuCC 1pi 533.342 533.342 0 | FGD1 anti-numuCC other 1023.3 1023.3 0 | FGD2 anti-numuCC 0pi 6283.55 6283.55 0 | FGD2 anti-numuCC 1pi 483.358 483.358 0 | FGD2 anti-numuCC other 944.017 944.017 0 | FGD1 NuMuBkg CC0pi in AntiNu Mode 2500.77 2500.77 0 | FGD1 NuMuBkg CC1pi in AntiNu Mode 845.922 845.922 0 | FGD1 NuMuBkg CCother in AntiNu Mode800.914 800.914 0 | FGD2 NuMuBkg CC0pi in AntiNu Mode 2554.46 2554.46 0 | FGD2 NuMuBkg CC1pi in AntiNu Mode 670.813 670.813 0 | FGD2 NuMuBkg CCother in AntiNu Mode785.661 785.661 0 | Total 114830 114830 0 | Saving data and MC distributions in try_multi_pi_29march_withcovcorr_ND280_nom.root ------------------------------ Adding samplePDF object Adding systematic object total_flux_cov Adding systematic object xsec_cov Adding systematic object nddet_cov ------------------------------ --------------------------------- General settings Output file: try_multi_pi_29march_withcovcorr.root Seed: 0 Debug: 0 GPU fit: 1 CPU MP fit: 1 N cores: 8 PlotByMode: 0 --------------------------------- MCMC settings N steps: 0 Flux step scale: 0.05 Xsec step scale: 0.02 ND det step scale:0.04 SK det stet scale:1 Osc step scale: 0.1 Anneal. temp: -1 Random start: 0 Process MCMC: 0 --------------------------------- T2K run settings ND runs: P6_S17 2a 2w 3b 3c 4a 4w 5 6b 6c 6d 6e 7 8a 8w SK nu POT: 3.01e+20 SK antinu POT: 3.01e+20 Livetime: 1417.4 Save nominal: 1 --------------------------------- Samples settings ND280: 1 FHC nu_mu CC0pi: 1 FHC nu_e CC0pi: 1 FHC nu_e CC1pi: 1 RHC nu_mu CC0pi: 1 RHC nu_e CC0pi: 1 SK SF->RFG tuning:1 SK RPA tuning: 1 SK CCCoh tuning: 1 --------------------------------- Samples at ND280 Samples binning ND280 Bin anu mom qe: Bin anu th qe: Bin anu mom nqe: Bin anu th nqe: Bin anu numu p qe: Bin anu numu th qe: Bin anu numu p nqe: Bin anu numu th nqe: Samples binning SK Numu bins Erec: 0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 0.55 0.6 0.65 0.7 0.75 0.8 0.85 0.9 0.95 1 1.05 1.1 1.15 1.2 1.25 1.3 1.35 1.4 1.45 1.5 1.55 1.6 1.65 1.7 1.75 1.8 1.85 1.9 1.95 2 2.05 2.1 2.15 2.2 2.25 2.3 2.35 2.4 2.45 2.5 2.55 2.6 2.65 2.7 2.75 2.8 2.85 2.9 2.95 3 3.25 3.5 3.75 4 4.5 5 5.5 6 7 8 9 10 30 Nue binning option: 0 Nue bins Erec: 0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 0.55 0.6 0.65 0.7 0.75 0.8 0.85 0.9 0.95 1 1.05 1.1 1.15 1.2 1.25 --------------------------------- Input settings Data file: --------------------------------- Real data settings RealDataFitFlag: 0 --------------------------------- Fake data settings FakeDataFitFlag: 1 FakeDataFile: --------------------------------- Asimov settings AsimovFitFlag: 0 NDAsimovFitFlag: 0 MakeAsimovOnly: 0 --------------------------------- Toy fit settings Post. pred. file: EMPTY ToyFitFlag: 0 Toy file: N toys: 20000 Full LLH: 1 --------------------------------- Systematics settings Flux syst: 1 Flux flat: Flux fixed: X-sec syst: 1 X-sec flat: 0 1 2 3 4 X-sec fixed: 12 ND280 detector syst: 1 ND280 detector flat: ND280 detector fixed: Sk detector syst: 1 Sk detector flat: Sk detector fixed: Oscill. flat: Oscill. fixed: Oscill. flat antinu: Oscill. fixed antinu: Reactor const.: 0 Stat only?: 0 Do oscillation?: 0 Use SK syst?: 1 Use BANFF tuning?:0 --------------------------------- Covariance matrices Flux file: inputs/flux_covariance_banff_13av2.root flux matrix: total_flux_cov Xsec file: inputs/xsec_covariance_2018a_NOMINAL_v11.root xsec matrix: xsec_cov BANFF file: BANFF matrix: BANFF vector: ND280 file: inputs/ND280_MultiPi_1p1h_2to8_Cov_WithCorr_29March_convertedToMaCh3.root ND280 matrix: nddet_cov SK file: utils/skdet_total_error_matrix_2015/SKJointErrorMatrix2015_Total_v7_13thMarch2016.root SK matrix: SKJointError_Erec_Total Osc file: inputs/oscillation_covariance_6par_nondouble.root osc matrix: osc_cov Osc file antinu: osc_cov osc matrix: Use BANFF covariance instead of ND data (and flux/xsec cov): 0 --------------------------------- Oscillation parameters Double angle : 0 Binned oscillation : 0 Oscillation binning: Starting parameters: 0.304 0.528 0.0217 7.53e-05 0.002509 -1.601 Starting parameters RHC: Input parameters (if different from starting parameters): Flat likelihood parameters: --------------------------------- Beta options: Use Beta : 0 Flip Beta : 0 Apply Beta to Nue : 0 Apply Diagnal Beta : 0 -------------------- Starting MCMC -------------------- 0 steps took 6.91414e-06 seconds to complete. (infs / step). 0 steps were accepted. ------------------------------ Sample Data MC -LLH | FGD1 numuCC 0pi 31529.2 32110.7 5.46497 | FGD1 numuCC 1pi 7997.52 8046.45 0.182899 | FGD1 numuCC other 6793.45 6807.02 0.0857593 | FGD2 numuCC 0pi 31732.1 32286.6 4.94668 | FGD2 numuCC 1pi 6418.42 6462.7 0.174218 | FGD2 numuCC other 6562.61 6590.09 0.0912648 | FGD1 anti-numuCC 0pi 6370.47 6485.39 1.12284 | FGD1 anti-numuCC 1pi 533.342 538.547 0.0292505 | FGD1 anti-numuCC other 1023.3 1029.97 0.0297489 | FGD2 anti-numuCC 0pi 6283.55 6396.17 1.06671 | FGD2 anti-numuCC 1pi 483.358 488.415 0.0309555 | FGD2 anti-numuCC other 944.017 951.6 0.0394359 | FGD1 NuMuBkg CC0pi in AntiNu Mode 2500.77 2539.36 0.315541 | FGD1 NuMuBkg CC1pi in AntiNu Mode 845.922 850.561 0.0260821 | FGD1 NuMuBkg CCother in AntiNu Mode800.914 802.533 0.00562613| FGD2 NuMuBkg CC0pi in AntiNu Mode 2554.46 2593.43 0.315413 | FGD2 NuMuBkg CC1pi in AntiNu Mode 670.813 675.02 0.0182555 | FGD2 NuMuBkg CCother in AntiNu Mode785.661 787.72 0.00682668| Total 114830 116442 13.9525 | ------------------------------