OscProb
Functions
Utils.h File Reference
#include "colormod.h"
#include "../tutorial/SetNiceStyle.C"
#include "TMath.h"
#include "TFile.h"
#include "PMNS_OQS.h"
#include "PMNS_LIV.h"
#include "PMNS_Deco.h"
#include "PMNS_Iter.h"
#include "PMNS_NUNM.h"
#include "PMNS_SNSI.h"
#include "PMNS_Decay.h"
#include "PMNS_Sterile.h"
#include "PMNS_SiderealLIV.h"

Go to the source code of this file.

Functions

void SetNominalPars (OscProb::PMNS_Base *p)
 
OscProb::PMNS_Fast * GetFast (bool is_nominal)
 
OscProb::PMNS_Iter * GetIter (bool is_nominal)
 
OscProb::PMNS_Deco * GetDeco (bool is_nominal)
 
OscProb::PMNS_OQS * GetOQS (bool is_nominal)
 
OscProb::PMNS_Sterile * GetSterile (bool is_nominal)
 
OscProb::PMNS_Decay * GetDecay (bool is_nominal)
 
OscProb::PMNS_NSI * GetNSI (bool is_nominal)
 
OscProb::PMNS_SNSI * GetSNSI (bool is_nominal)
 
OscProb::PMNS_LIV * GetLIV (bool is_nominal)
 
OscProb::PMNS_SiderealLIV * GetSiderealLIV (bool is_nominal)
 
OscProb::PMNS_NUNM * GetNUNM (bool is_nominal)
 
OscProb::PMNS_Base * GetModel (string model, bool is_nominal=false)
 
vector< string > GetListOfModels ()
 
void SetTestPath (OscProb::PMNS_Base *p)
 
void SaveTestFile (OscProb::PMNS_Base *p, TString filename)
 
int CheckProb (OscProb::PMNS_Base *p, TString filename)
 

Function Documentation

◆ CheckProb()

int CheckProb ( OscProb::PMNS_Base *  p,
TString  filename 
)

Definition at line 281 of file Utils.h.

281 {
282
283 SetTestPath(p);
284
285 TFile* f = TFile::Open("data/"+filename, "read");
286
287 if(!f){
288 printf((Color::FAILED + " data/%s not found\n").c_str(), filename.Data());
289 return 1;
290 }
291
292 int ntests = 0;
293 int fails = 0;
294
295 TCanvas* c1 = 0;
296 TH1D* h0 = 0;
297 TH1D* h = 0;
298
299 for(int flvi=0; flvi<3; flvi++){
300 for(int flvf=0; flvf<3; flvf++){
301 for(int isnb=0; isnb<2; isnb++){
302 p->SetIsNuBar(isnb);
303 TString hname = TString::Format("h%d%d%d",flvi,flvf,isnb);
304 h0 = (TH1D*)f->Get(hname);
305 h = (TH1D*)h0->Clone();
306 bool plot = false;
307 for(int i=1; i<=h0->GetNbinsX(); i++){
308 double energy = h->GetBinCenter(i);
309 double dE = h->GetBinWidth(i);
310 double p0 = h0->GetBinContent(i);
311 double p1 = p->AvgProb(flvi, flvf, energy, dE);
312 ntests++;
313 if(abs(p0-p1)>1e-12){
314 plot = true;
315 fails++;
316 }
317 h->SetBinContent(i, p1);
318 }
319 if(plot){
320 c1 = new TCanvas();
321 c1->Divide(1,2);
322 c1->cd(1);
323 SetHist(h0, kBlue);
324 SetHist(h, kRed);
325 TString nu_lab = isnb ? "#bar{#nu}" : "#nu";
326 TString flv_lab[3] = {"e","#mu","#tau"};
327 TString ylab = "P(" + nu_lab + "_{" + flv_lab[flvi] +
328 "}#rightarrow" + nu_lab + "_{" + flv_lab[flvf] + "})";
329 h0->SetTitle(";Energy [GeV];"+ylab+";");
330 h->SetLineStyle(7);
331 double ymax = max(h->GetMaximum(),
332 h0->GetMaximum());
333 double ymin = min(h->GetMinimum(),
334 h0->GetMinimum());
335 h0->GetYaxis()->SetRangeUser(ymin, ymax);
336 h0->DrawCopy("hist");
337 h->DrawCopy("hist same");
338 SetTH1Margin();
339 gPad->SetLogx();
340 c1->cd(2);
341 h->SetTitle(";Energy [GeV];#Delta"+ylab+";");
342 h->Add(h0,-1);
343 h->SetLineStyle(1);
344 h->DrawCopy("hist");
345 SetTH1Margin();
346 gPad->SetLogx();
347 MiscText(0.55,0.8,0.1,filename,kGray,1);
348 c1->DrawClone();
349 TString pngfile = filename;
350 pngfile.ReplaceAll(".root",".png");
351 c1->SaveAs("plots/Failed_"+hname+"_"+pngfile);
352 delete c1;
353 }
354 if(h0) delete h0;
355 if(h) delete h;
356 }}}
357
358 if(fails>0){
359 printf((Color::FAILED + " Found %d differences in %d tests (%.3g%%) in %s\n").c_str(),
360 fails, ntests, 100.*fails/ntests, filename.Data());
361 }
362 else {
363 cout << Color::PASSED << " No differences found in " << filename << endl;
364 }
365
366 return fails;
367
368}
SetTestPath
void SetTestPath(OscProb::PMNS_Base *p)
Definition: Utils.h:241
OscProb::PMNS_Base::SetIsNuBar
virtual void SetIsNuBar(bool isNuBar)
Set the anti-neutrino flag.
Definition: PMNS_Base.cxx:243
OscProb::PMNS_Base::AvgProb
virtual double AvgProb(vectorC nu_in, int flvf, double E, double dE=0)
Compute the average probability over a bin of energy.
Definition: PMNS_Base.cxx:1568
Color::PASSED
static const string PASSED
Definition: colormod.h:23
Color::FAILED
static const string FAILED
Definition: colormod.h:22

References OscProb::PMNS_Base::AvgProb(), Color::FAILED, Color::PASSED, OscProb::PMNS_Base::SetIsNuBar(), and SetTestPath().

◆ GetDecay()

OscProb::PMNS_Decay * GetDecay ( bool  is_nominal)

Definition at line 96 of file Utils.h.

96 {
97
98 OscProb::PMNS_Decay* p = new OscProb::PMNS_Decay();
99 SetNominalPars(p);
100 if(!is_nominal){
101 p->SetAlpha3(1e-4);
102 }
103
104 return p;
105
106}
SetNominalPars
void SetNominalPars(OscProb::PMNS_Base *p)
Definition: Utils.h:20
OscProb::PMNS_Decay
Implementation of neutrino decay in a three-neutrino framework.
Definition: PMNS_Decay.h:30
OscProb::PMNS_Decay::SetAlpha3
virtual void SetAlpha3(double alpha3)
Definition: PMNS_Decay.cxx:62

References OscProb::PMNS_Decay::SetAlpha3(), and SetNominalPars().

Referenced by GetModel().

◆ GetDeco()

OscProb::PMNS_Deco * GetDeco ( bool  is_nominal)

Definition at line 51 of file Utils.h.

51 {
52
53 OscProb::PMNS_Deco* p = new OscProb::PMNS_Deco();
54 SetNominalPars(p);
55 if(!is_nominal){
56 p->SetGamma(2, 1e-23);
57 p->SetGamma(3, 1e-22);
58 }
59
60 return p;
61
62}
OscProb::PMNS_Deco
Implementation of oscillations of neutrinos in matter in a three-neutrino framework with decoherence.
Definition: PMNS_Deco.h:27
OscProb::PMNS_Deco::SetGamma
virtual void SetGamma(int j, double val)
Set any given decoherence parameter.
Definition: PMNS_Deco.cxx:49

References OscProb::PMNS_Deco::SetGamma(), and SetNominalPars().

Referenced by GetModel().

◆ GetFast()

OscProb::PMNS_Fast * GetFast ( bool  is_nominal)

Definition at line 33 of file Utils.h.

33 {
34
35 OscProb::PMNS_Fast* p = new OscProb::PMNS_Fast();
36 SetNominalPars(p);
37 return p;
38
39}
OscProb::PMNS_Fast
Implementation of oscillations of neutrinos in matter in a three-neutrino framework.
Definition: PMNS_Fast.h:40

References SetNominalPars().

Referenced by GetModel().

◆ GetIter()

OscProb::PMNS_Iter * GetIter ( bool  is_nominal)

Definition at line 42 of file Utils.h.

42 {
43
44 OscProb::PMNS_Iter* p = new OscProb::PMNS_Iter();
45 SetNominalPars(p);
46 return p;
47
48}
OscProb::PMNS_Iter
Implementation of oscillations of neutrinos in matter in a three-neutrino framework.
Definition: PMNS_Iter.h:35

References SetNominalPars().

Referenced by GetModel().

◆ GetListOfModels()

vector< string > GetListOfModels ( )

Definition at line 232 of file Utils.h.

232 {
233
234 return {"Fast", "Iter", "Sterile", "NSI",
235 "Deco", "Decay", "LIV", "SNSI",
236 "NUNM", "OQS", "SiderealLIV"};
237
238}

◆ GetLIV()

OscProb::PMNS_LIV * GetLIV ( bool  is_nominal)

Definition at line 145 of file Utils.h.

145 {
146
147 OscProb::PMNS_LIV* p = new OscProb::PMNS_LIV();
148 SetNominalPars(p);
149 if(!is_nominal){
150 p->SetaT(0,0, 0.1e-22, 0);
151 p->SetaT(0,1, 0.2e-22, 0);
152 p->SetaT(0,2, 0.3e-22, 0);
153 p->SetaT(1,1, 0.4e-22, 0);
154 p->SetaT(1,2, 0.5e-22, 0);
155 p->SetaT(2,2, 0.6e-22, 0);
156 p->SetcT(0,0, 0.1e-22, 0);
157 p->SetcT(0,1, 0.2e-22, 0);
158 p->SetcT(0,2, 0.3e-22, 0);
159 p->SetcT(1,1, 0.4e-22, 0);
160 p->SetcT(1,2, 0.5e-22, 0);
161 p->SetcT(2,2, 0.6e-22, 0);
162 }
163
164 return p;
165
166}
OscProb::PMNS_LIV
Implements oscillations with LIV as modelled by SME.
Definition: PMNS_LIV.h:28
OscProb::PMNS_LIV::SetaT
virtual void SetaT(int flvi, int flvj, int dim, double val, double phase)
Definition: PMNS_LIV.cxx:68
OscProb::PMNS_LIV::SetcT
virtual void SetcT(int flvi, int flvj, int dim, double val, double phase)
Definition: PMNS_LIV.cxx:120

References OscProb::PMNS_LIV::SetaT(), OscProb::PMNS_LIV::SetcT(), and SetNominalPars().

Referenced by GetModel().

◆ GetModel()

OscProb::PMNS_Base * GetModel ( string  model,
bool  is_nominal = false 
)

Definition at line 214 of file Utils.h.

214 {
215
216 if(model == "Iter") return GetIter(is_nominal);
217 if(model == "Deco") return GetDeco(is_nominal);
218 if(model == "Sterile") return GetSterile(is_nominal);
219 if(model == "Decay") return GetDecay(is_nominal);
220 if(model == "NSI") return GetNSI(is_nominal);
221 if(model == "LIV") return GetLIV(is_nominal);
222 if(model == "SNSI") return GetSNSI(is_nominal);
223 if(model == "NUNM") return GetNUNM(is_nominal);
224 if(model == "OQS") return GetOQS(is_nominal);
225 if(model == "SiderealLIV") return GetSiderealLIV(is_nominal);
226
227 return GetFast(is_nominal);
228
229}
GetNUNM
OscProb::PMNS_NUNM * GetNUNM(bool is_nominal)
Definition: Utils.h:196
GetLIV
OscProb::PMNS_LIV * GetLIV(bool is_nominal)
Definition: Utils.h:145
GetFast
OscProb::PMNS_Fast * GetFast(bool is_nominal)
Definition: Utils.h:33
GetDeco
OscProb::PMNS_Deco * GetDeco(bool is_nominal)
Definition: Utils.h:51
GetSterile
OscProb::PMNS_Sterile * GetSterile(bool is_nominal)
Definition: Utils.h:80
GetIter
OscProb::PMNS_Iter * GetIter(bool is_nominal)
Definition: Utils.h:42
GetNSI
OscProb::PMNS_NSI * GetNSI(bool is_nominal)
Definition: Utils.h:109
GetSiderealLIV
OscProb::PMNS_SiderealLIV * GetSiderealLIV(bool is_nominal)
Definition: Utils.h:169
GetDecay
OscProb::PMNS_Decay * GetDecay(bool is_nominal)
Definition: Utils.h:96
GetSNSI
OscProb::PMNS_SNSI * GetSNSI(bool is_nominal)
Definition: Utils.h:127
GetOQS
OscProb::PMNS_OQS * GetOQS(bool is_nominal)
Definition: Utils.h:65

References GetDecay(), GetDeco(), GetFast(), GetIter(), GetLIV(), GetNSI(), GetNUNM(), GetOQS(), GetSiderealLIV(), GetSNSI(), and GetSterile().

◆ GetNSI()

OscProb::PMNS_NSI * GetNSI ( bool  is_nominal)

Definition at line 109 of file Utils.h.

109 {
110
111 OscProb::PMNS_NSI* p = new OscProb::PMNS_NSI();
112 SetNominalPars(p);
113 if(!is_nominal){
114 p->SetEps(0,0, 0.1, 0);
115 p->SetEps(0,1, 0.2, 0);
116 p->SetEps(0,2, 0.3, 0);
117 p->SetEps(1,1, 0.4, 0);
118 p->SetEps(1,2, 0.5, 0);
119 p->SetEps(2,2, 0.6, 0);
120 }
121
122 return p;
123
124}
OscProb::PMNS_NSI
Implementation of oscillations of neutrinos in matter in a three-neutrino framework with NSI.
Definition: PMNS_NSI.h:28
OscProb::PMNS_NSI::SetEps
virtual void SetEps(int flvi, int flvj, double val, double phase)
Set any given NSI parameter.
Definition: PMNS_NSI.cxx:86

References OscProb::PMNS_NSI::SetEps(), and SetNominalPars().

Referenced by GetModel().

◆ GetNUNM()

OscProb::PMNS_NUNM * GetNUNM ( bool  is_nominal)

Definition at line 196 of file Utils.h.

196 {
197
198 OscProb::PMNS_NUNM* p = new OscProb::PMNS_NUNM(0);
199 SetNominalPars(p);
200 if(!is_nominal){
201 p->SetAlpha(0,0, 0.05, 0);
202 p->SetAlpha(1,0, 0.06, 0);
203 p->SetAlpha(2,0, 0.07, 0);
204 p->SetAlpha(1,1, 0.08, 0);
205 p->SetAlpha(2,1, 0.09, 0);
206 p->SetAlpha(2,2, 0.1, 0);
207 }
208
209 return p;
210
211}
OscProb::PMNS_NUNM
Implementation of oscillations of neutrinos in matter in a three-neutrino framework with Non unitary ...
Definition: PMNS_NUNM.h:30
OscProb::PMNS_NUNM::SetAlpha
virtual void SetAlpha(int i, int j, double val, double phase)
Set any given NUNM parameter.
Definition: PMNS_NUNM.cxx:95

References OscProb::PMNS_NUNM::SetAlpha(), and SetNominalPars().

Referenced by GetModel().

◆ GetOQS()

OscProb::PMNS_OQS * GetOQS ( bool  is_nominal)

Definition at line 65 of file Utils.h.

65 {
66
67 OscProb::PMNS_OQS* p = new OscProb::PMNS_OQS();
68 SetNominalPars(p);
69 if(!is_nominal){
70 p->SetDecoElement(3, sqrt(2e-23));
71 p->SetDecoElement(8, sqrt(4e-23));
72 p->SetDecoAngle(3,8, acos(0.5));
73 }
74
75 return p;
76
77}
OscProb::PMNS_OQS
Implements neutrino oscillations using an open quantum system approach.
Definition: PMNS_OQS.h:28
OscProb::PMNS_OQS::SetDecoAngle
virtual void SetDecoAngle(int i, int j, double th)
Set mixing angle between two decoherence parameters a_i, a_j.
Definition: PMNS_OQS.cxx:75
OscProb::PMNS_OQS::SetDecoElement
virtual void SetDecoElement(int i, double val)
Set value of the a_i decoherence element in Gell-Mann basis.
Definition: PMNS_OQS.cxx:58

References OscProb::PMNS_OQS::SetDecoAngle(), OscProb::PMNS_OQS::SetDecoElement(), and SetNominalPars().

Referenced by GetModel().

◆ GetSiderealLIV()

OscProb::PMNS_SiderealLIV * GetSiderealLIV ( bool  is_nominal)

Definition at line 169 of file Utils.h.

169 {
170
171 OscProb::PMNS_SiderealLIV* p = new OscProb::PMNS_SiderealLIV();
172 SetNominalPars(p);
173 if(!is_nominal){
174 p->SetA(0,0, 0, 0.1e-22);
175 p->SetA(0,1, 1, 0.2e-22);
176 p->SetA(0,2, 2, 0.3e-22);
177 p->SetA(1,1, 0, 0.4e-22);
178 p->SetA(1,2, 1, 0.5e-22);
179 p->SetA(2,2, 2, 0.6e-22);
180 p->SetC(0,0, 0,0, 0.1e-22);
181 p->SetC(0,1, 1,1, 0.2e-22);
182 p->SetC(0,2, 2,2, 0.3e-22);
183 p->SetC(1,1, 0,1, 0.4e-22);
184 p->SetC(1,2, 1,2, 0.5e-22);
185 p->SetC(2,2, 0,2, 0.6e-22);
186 p->SetColatitude(-89, -59, -24); // IceCube (South Pole)
187 p->SetNeutrinoDirection(57.3, 28.6); // ~1.0 rad, ~0.5 rad
188 p->SetTimeHours(6.0);
189 }
190
191 return p;
192
193}
OscProb::PMNS_SiderealLIV
Implements oscillations with Sidereal LIV as modelled by SME.
Definition: PMNS_SiderealLIV.h:26
OscProb::PMNS_SiderealLIV::SetA
virtual void SetA(int flvi, int flvj, int coord, double val)
Definition: PMNS_SiderealLIV.cxx:97
OscProb::PMNS_SiderealLIV::SetColatitude
virtual void SetColatitude(double chi)
Definition: PMNS_SiderealLIV.cxx:221
OscProb::PMNS_SiderealLIV::SetC
virtual void SetC(int flvi, int flvj, int coord1, int coord2, double val)
Definition: PMNS_SiderealLIV.cxx:120
OscProb::PMNS_SiderealLIV::SetNeutrinoDirection
virtual void SetNeutrinoDirection(double zenith, double azimuth)
Definition: PMNS_SiderealLIV.cxx:187
OscProb::PMNS_SiderealLIV::SetTimeHours
virtual void SetTimeHours(double hours)
Definition: PMNS_SiderealLIV.cxx:209

References OscProb::PMNS_SiderealLIV::SetA(), OscProb::PMNS_SiderealLIV::SetC(), OscProb::PMNS_SiderealLIV::SetColatitude(), OscProb::PMNS_SiderealLIV::SetNeutrinoDirection(), SetNominalPars(), and OscProb::PMNS_SiderealLIV::SetTimeHours().

Referenced by GetModel().

◆ GetSNSI()

OscProb::PMNS_SNSI * GetSNSI ( bool  is_nominal)

Definition at line 127 of file Utils.h.

127 {
128
129 OscProb::PMNS_SNSI* p = new OscProb::PMNS_SNSI();
130 SetNominalPars(p);
131 if(!is_nominal){
132 p->SetEps(0,0, 0.1, 0);
133 p->SetEps(0,1, 0.2, 0);
134 p->SetEps(0,2, 0.3, 0);
135 p->SetEps(1,1, 0.4, 0);
136 p->SetEps(1,2, 0.5, 0);
137 p->SetEps(2,2, 0.6, 0);
138 }
139
140 return p;
141
142}
OscProb::PMNS_SNSI
Implementation of oscillations of neutrinos in matter in a three-neutrino framework with scalar NSI.
Definition: PMNS_SNSI.h:26

References OscProb::PMNS_NSI::SetEps(), and SetNominalPars().

Referenced by GetModel().

◆ GetSterile()

OscProb::PMNS_Sterile * GetSterile ( bool  is_nominal)

Definition at line 80 of file Utils.h.

80 {
81
82 OscProb::PMNS_Sterile* p = new OscProb::PMNS_Sterile(4);
83 SetNominalPars(p);
84 if(!is_nominal){
85 p->SetDm(4, 0.1);
86 p->SetAngle(1,4, 0.1);
87 p->SetAngle(2,4, 0.1);
88 p->SetAngle(3,4, 0.1);
89 }
90
91 return p;
92
93}
OscProb::PMNS_Base::SetDm
virtual void SetDm(int j, double dm)
Set the mass-splitting dm_j1 in eV^2.
Definition: PMNS_Base.cxx:674
OscProb::PMNS_Base::SetAngle
virtual void SetAngle(int i, int j, double th)
Set the mixing angle theta_ij.
Definition: PMNS_Base.cxx:539
OscProb::PMNS_Sterile
Implementation of oscillations of neutrinos in matter in a N-neutrino framework.
Definition: PMNS_Sterile.h:34

References OscProb::PMNS_Base::SetAngle(), OscProb::PMNS_Base::SetDm(), and SetNominalPars().

Referenced by GetModel().

◆ SaveTestFile()

void SaveTestFile ( OscProb::PMNS_Base *  p,
TString  filename 
)

Definition at line 250 of file Utils.h.

250 {
251
252 SetTestPath(p);
253
254 int nbins = 100;
255 vector<double> xbins = GetLogAxis(nbins, 0.1, 10);
256 TH1D* h = 0;
257
258 TFile* f = new TFile("data/"+filename, "recreate");
259
260 for(int flvi=0; flvi<3; flvi++){
261 for(int flvf=0; flvf<3; flvf++){
262 for(int isnb=0; isnb<2; isnb++){
263 p->SetIsNuBar(isnb);
264 TString hname = TString::Format("h%d%d%d",flvi,flvf,isnb);
265 h = new TH1D(hname, "", nbins, &xbins[0]);
266 for(int i=1; i<=nbins; i++){
267 double energy = h->GetBinCenter(i);
268 double dE = h->GetBinWidth(i);
269 h->SetBinContent(i, p->AvgProb(flvi, flvf, energy, dE));
270 }
271 h->Write();
272 delete h;
273 }}}
274
275 f->Close();
276 cout << "Saved new test file: data/" + filename << endl;
277
278}

References OscProb::PMNS_Base::AvgProb(), OscProb::PMNS_Base::SetIsNuBar(), and SetTestPath().

◆ SetNominalPars()

void SetNominalPars ( OscProb::PMNS_Base *  p)

Definition at line 20 of file Utils.h.

20 {
21
22 // Set to NuFIT 5.2 values (NO w/ SK)
23 p->SetDm(2, 7.41e-5);
24 p->SetDm(3, 2.507e-3);
25 p->SetAngle(1,2, asin(sqrt(0.303)));
26 p->SetAngle(1,3, asin(sqrt(0.02225)));
27 p->SetAngle(2,3, asin(sqrt(0.451)));
28 p->SetDelta(1,3, 232 * TMath::DegToRad());
29
30}
OscProb::PMNS_Base::SetDelta
virtual void SetDelta(int i, int j, double delta)
Set the CP phase delta_ij.
Definition: PMNS_Base.cxx:602

References OscProb::PMNS_Base::SetAngle(), OscProb::PMNS_Base::SetDelta(), and OscProb::PMNS_Base::SetDm().

Referenced by GetDecay(), GetDeco(), GetFast(), GetIter(), GetLIV(), GetNSI(), GetNUNM(), GetOQS(), GetSiderealLIV(), GetSNSI(), and GetSterile().

◆ SetTestPath()

void SetTestPath ( OscProb::PMNS_Base *  p)

Definition at line 241 of file Utils.h.

241 {
242
243 p->SetPath(1000, 2);
244 p->AddPath(1000, 4);
245 p->AddPath(1000, 2);
246
247}
OscProb::PMNS_Base::AddPath
virtual void AddPath(NuPath p)
Add a path to the sequence.
Definition: PMNS_Base.cxx:307
OscProb::PMNS_Base::SetPath
virtual void SetPath(NuPath p)
Set a single path.
Definition: PMNS_Base.cxx:330

References OscProb::PMNS_Base::AddPath(), and OscProb::PMNS_Base::SetPath().

Referenced by CheckProb(), and SaveTestFile().