Non-standard physics and Non-standard physics and user-defined priors in user-defined priors in GLoBES GLoBES Workshop on physics and applications of Workshop on physics and applications of the GLoBES software the GLoBES software Max-Planck-Institut für Kernphysik Max-Planck-Institut für Kernphysik January 24, 2007 January 24, 2007 Walter Winter Walter Winter Universität Würzburg Universität Würzburg
21
Embed
Non-standard physics and user- defined priors in GLoBES Workshop on physics and applications of the GLoBES software Workshop on physics and applications.
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
Non-standard physics and user-Non-standard physics and user-defined priors in GLoBESdefined priors in GLoBES
Workshop on physics and applications of the Workshop on physics and applications of the GLoBES softwareGLoBES software
Max-Planck-Institut für KernphysikMax-Planck-Institut für Kernphysik
January 24, 2007January 24, 2007
Walter WinterWalter WinterUniversität WürzburgUniversität Würzburg
January 24, 2007 GLoBES Workshop - Walter Winter
2
OutlineOutline Introduction:Introduction:
Different levels in GLoBESDifferent levels in GLoBES Probability level and the simulation of Probability level and the simulation of
non-standard physicsnon-standard physics Systematics levelSystematics level Physics level and user-defined priorsPhysics level and user-defined priors SummarySummary
January 24, 2007 GLoBES Workshop - Walter Winter
3
Different levels in GLoBESDifferent levels in GLoBES
Probability level
Systematics level
Correlation anddegeneracy level
Additional “nuisance“ parameters. Example:
(Fogli et al, 2002)
Projection onto sub-space/marginalization:
Channel
Rule
Experi-ment(s)
AEDL:
January 24, 2007 GLoBES Workshop - Walter Winter
4
“Evolution“ operator in one layer, being diagonalized:
Probability levelProbability levelHamiltonian in constant matter density layer:
Probability calculation:
January 24, 2007 GLoBES Workshop - Walter Winter
5
Non-standard physics in GLoBESNon-standard physics in GLoBES1. Modify the Hamiltonian, probability calculation, etc. Example: Non-standard matter effect in e--sector:
2. GLoBES now carries k more oscillation parameters, which need to be maintained:
January 24, 2007 GLoBES Workshop - Walter Winter
6
Example: Decoherence at reactor exp.Example: Decoherence at reactor exp. Oscillation probabilities with damping effects:Oscillation probabilities with damping effects:
Damping factors for (wave packet) decoherence:Damping factors for (wave packet) decoherence:
E: Intrinsic wave packet widthCan be implemented analytically in quasi-vacuum (short L)
(Blennow, Ohlsson, Winter, hep-ph/0502147)
(J. Kopp for GLoBES 3.0)
January 24, 2007 GLoBES Workshop - Walter Winter
7
Modification of the access functionsModification of the access functions
const int const int GLB_SIGMA_E = 6GLB_SIGMA_E = 6;;double th12, th13,th23,deltacp,sdm,ldm,double th12, th13,th23,deltacp,sdm,ldm,sigma_Esigma_E;;
Define access functions for the oscillation parameters:Define access functions for the oscillation parameters:
Set oscillationparameters
from internalstructure
Get oscillationparameters
from internalstructure
January 24, 2007 GLoBES Workshop - Walter Winter
8
Definition of a probability engineDefinition of a probability engineint my_probability_matrix(double P[3][3], int cp_sign, double E, int psteps,int my_probability_matrix(double P[3][3], int cp_sign, double E, int psteps, const double *length, const double *density, const double *length, const double *density, double filter_sigma, void *user_data) double filter_sigma, void *user_data){{ int i, j; int i, j; double L, Delta21, Delta31, Delta32, D21, D31, D32, s13, c13, s12, c12, t; double L, Delta21, Delta31, Delta32, D21, D31, D32, s13, c13, s12, c12, t; /* Set all probabilities to zero initially *//* Set all probabilities to zero initially */ for (i=0; i < 3; i++) for (j=0; j < 3; j++) P[i][j] = 0.0; for (i=0; i < 3; i++) for (j=0; j < 3; j++) P[i][j] = 0.0;
/* Calculate total baseline from the input lists*//* Calculate total baseline from the input lists*/ L = 0.0; for (i=0; i < psteps; i++) L += length[i]; L = 0.0; for (i=0; i < psteps; i++) L += length[i]; L = GLB_KM_TO_EV(L) * 1.0e9; /* Convert to GeV^{-1} */ L = GLB_KM_TO_EV(L) * 1.0e9; /* Convert to GeV^{-1} */
/* Compute P_ee analytically with a piece of code*//* Compute P_ee analytically with a piece of code*/ s12 = sin(th12);c12 = cos(th12);s13 = sin(th13); s12 = sin(th12);c12 = cos(th12);s13 = sin(th13);c13 = cos(th13);c13 = cos(th13); t = L / (4.0 * E); t = L / (4.0 * E); Delta21 = sdm * t;Delta31 = ldm * t;Delta32 = Delta31 - Delta21; Delta21 = sdm * t;Delta31 = ldm * t;Delta32 = Delta31 - Delta21; t = M_SQRT2 * sigma_E / E; t = M_SQRT2 * sigma_E / E; D21 = exp(-square( Delta21 * t )); D31 = exp(-square( Delta31 * t )); D21 = exp(-square( Delta21 * t )); D31 = exp(-square( Delta31 * t )); D32 = exp(-square( Delta32 * t )); D32 = exp(-square( Delta32 * t )); P[0][0]P[0][0] = square(square(c13)) * ( 1 - 2.0*square(s12*c12)* (1–D21* = square(square(c13)) * ( 1 - 2.0*square(s12*c12)* (1–D21* cos(2.0*Delta21))) + 2.0*square(s13*c13) * ( D31*square(c12) * cos(2.0*Delta21))) + 2.0*square(s13*c13) * ( D31*square(c12) * cos(2.0*Delta31) + D32*square(s12) * cos(2.0*Delta32) ) + cos(2.0*Delta31) + D32*square(s12) * cos(2.0*Delta32) ) + square(square(s13)); square(square(s13));
return 0; return 0;}}
(GLoBES 3.0, example6.c)
January 24, 2007 GLoBES Workshop - Walter Winter
9
Using the new probability engineUsing the new probability engine Register the new probability engine and access functionsRegister the new probability engine and access functions
glbRegisterProbabilityEngine(glbRegisterProbabilityEngine(7, // Number of params7, // Number of params &my_probability_matrix, &my_probability_matrix, &my_set_oscillation_parameters, &my_set_oscillation_parameters, &my_get_oscillation_parameters, &my_get_oscillation_parameters, NULL); NULL);
Maintain the new oscillation parameter(s):Maintain the new oscillation parameter(s):– Use Use glbSetOscParamsglbSetOscParams and and glbGetOscParamsglbGetOscParams to to
access the non-standard parameter(s)access the non-standard parameter(s)
– Always use Always use glbChiNPglbChiNP (instead of (instead of glbChiTheta13glbChiTheta13 etc.) etc.)to define how the non-standard degrees of freedom behaveto define how the non-standard degrees of freedom behave
– Define your non-standard behavior in the projection with Define your non-standard behavior in the projection with glbSetProjectionFlagglbSetProjectionFlag
January 24, 2007 GLoBES Workshop - Walter Winter
10
Results for our exampleResults for our example glbDefineProjection(myproj, GLB_FREE, glbDefineProjection(myproj, GLB_FREE, GLB_FIXED, GLB_FIXED, GLB_FIXED, GLB_FIXED, GLB_FIXED, GLB_FIXED, GLB_FREE, GLB_FREE); GLB_FREE, GLB_FREE); glbSetDensityProjectionFlag(myproj, glbSetDensityProjectionFlag(myproj, GLB_FIXED, GLB_ALL); GLB_FIXED, GLB_ALL); glbSetProjectionFlag(myproj,GLB_FIXED,GLB_SIGMA_E);glbSetProjectionFlag(myproj,GLB_FIXED,GLB_SIGMA_E); glbSetProjection(myproj); glbSetProjection(myproj);
for(x=0; x < 0.05+0.001; x+=0.005) /* th13 */ for(x=0; x < 0.05+0.001; x+=0.005) /* th13 */ for(y=0.0; y < 0.010+0.001; y+=0.001) /* sigma_E */ for(y=0.0; y < 0.010+0.001; y+=0.001) /* sigma_E */ { { /* Set vector of test=fit values */ /* Set vector of test=fit values */ thetheta13=asin(sqrt(x))/2.0; thetheta13=asin(sqrt(x))/2.0; glbSetOscParams(test_values,thetheta13, glbSetOscParams(test_values,thetheta13, GLB_THETA_13); GLB_THETA_13); glbSetOscParams(test_values,y,GLB_SIGMA_E);glbSetOscParams(test_values,y,GLB_SIGMA_E); /* Compute Chi^2 with correlations */ /* Compute Chi^2 with correlations */ res= res=glbChiNPglbChiNP(test_values,NULL,GLB_ALL);(test_values,NULL,GLB_ALL); AddToOutput(x,y,res); AddToOutput(x,y,res); } } True values: 13=e=0
(Original figure from hep-ph/0502147)
Correlations
January 24, 2007 GLoBES Workshop - Walter Winter
11
More physics examples computed with More physics examples computed with GLoBESGLoBES
Simulation of Simulation of Hamiltonian-level effectsHamiltonian-level effects
These ratios are changed through averaged These ratios are changed through averaged neutrino oscillations:neutrino oscillations:Only CP-conserving effects remaining ~ cos Only CP-conserving effects remaining ~ cos CPCP
Measure muon track to shower ratio at neutrino Measure muon track to shower ratio at neutrino telescope: R = telescope: R = /(/(ee))(conservative, since in future also flavors!?)(conservative, since in future also flavors!?)
January 24, 2007 GLoBES Workshop - Walter Winter
16
Complementarity to beamsComplementarity to beams Use R to obtainUse R to obtain
information oninformation onosc. parameters?osc. parameters?Difficult, sinceDifficult, since– Low statisticsLow statistics– No spectral infoNo spectral info
But: ComplementaryBut: Complementarydependence on dependence on CPCPHere: Constant-rates/ Here: Constant-rates/ constant-R curvesconstant-R curves
Combine the Combine the information from two information from two “low-statistics”“low-statistics”sources?sources? (Winter, hep-ph/0604191)(Winter, hep-ph/0604191)
Best-fit
January 24, 2007 GLoBES Workshop - Walter Winter
17
Implementation in GLoBESImplementation in GLoBESdouble R_neutron_true; /* Simulated/true R */double R_neutron_true; /* Simulated/true R */double relerror = 0.2; /* Relative error */double relerror = 0.2; /* Relative error */
Double ChoozDouble Choozcould be thecould be thefirst experiment tofirst experiment toobserve observe CPCP
(Winter, hep-ph/0604191)(Winter, hep-ph/0604191)
(1(1, 90% CL; 1 d.o.f.), 90% CL; 1 d.o.f.)
January 24, 2007 GLoBES Workshop - Walter Winter
19
Mass hierarchyMass hierarchy Astrophysical source may Astrophysical source may
help mass hierarchy help mass hierarchy measurement at superbeam:measurement at superbeam:20% prec. good20% prec. good (Winter, hep-ph/0604191)(Winter, hep-ph/0604191)
No ext. info
20%
5%10%
January 24, 2007 GLoBES Workshop - Walter Winter
20
Further physics applicationsFurther physics applications Combination with other external measurements, Combination with other external measurements,
such as atmospheric neutrinossuch as atmospheric neutrinos(Huber, Maltoni, Schwetz, hep-ph/0501037; (Huber, Maltoni, Schwetz, hep-ph/0501037; Campagne, Maltoni, Mezzetto, Schwetz, hep-ph/0603172)Campagne, Maltoni, Mezzetto, Schwetz, hep-ph/0603172)
Penalties for degeneracy localization:Penalties for degeneracy localization:E.g. add penalty if in wrong octantE.g. add penalty if in wrong octant(Schwetz priors)(Schwetz priors)