Optimized version of OscCalcPMNS. More...

#include "/cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-04-18/CAFAna/Core/StanTypedefs.h"

Inheritance diagram for osc::_OscCalcPMNSOpt< T >:

Classes
struct	Val_t

Public Member Functions
	_OscCalcPMNSOpt ()

virtual	~_OscCalcPMNSOpt ()

_IOscCalcAdjustable< T > *	Copy () const override

T	P (int flavBefore, int flavAfter, double E) override
	E in GeV; flavors as PDG codes (so, neg==>antinu) More...

void	SetL (double L) override

void	SetRho (double rho) override

void	SetDmsq21 (const T &dmsq21) override

void	SetDmsq32 (const T &dmsq32) override

void	SetTh13 (const T &th13) override

void	SetTh12 (const T &th12) override

void	SetTh23 (const T &th23) override

void	SetdCP (const T &dCP) override

TMD5 *	GetParamsHash () const override
	Use to check two calculators are in the same state. More...

virtual double	GetL () const

virtual double	GetRho () const

virtual T	GetDmsq21 () const

virtual T	GetDmsq32 () const

virtual T	GetTh12 () const

virtual T	GetTh13 () const

virtual T	GetTh23 () const

virtual T	GetdCP () const

virtual void	InvalidateCache ()
	Invalidate any caching used internally by the calculator. More...

virtual Eigen::Array< T, Eigen::Dynamic, 1 >	P (int flavBefore, int flavAfter, const std::vector< double > &E)

virtual Eigen::Array< T, Eigen::Dynamic, 1 >	P (int flavBefore, int flavAfter, const Eigen::ArrayXd &E)

Protected Member Functions
TMD5 *	GetParamsHashDefault (const std::string &txt) const
	This is only a safe implementation if your calculator only depends on these eight parameters. More...

Protected Attributes
long	fMixIdx

long	fDmIdx

long	fLRIdx

std::unordered_map< double, Val_t >	fPMNSOpt [2]

double	fRho

double	fL

T	fDmsq21

T	fDmsq32

T	fTh12

T	fTh13

T	fTh23

T	fdCP

Detailed Description

template<typename T>
class osc::_OscCalcPMNSOpt< T >

Optimized version of OscCalcPMNS.

Adapt the PMNSOpt calculator to standard interface

Definition at line 32 of file StanTypedefs.h.

Constructor & Destructor Documentation

template<typename T >

osc::_OscCalcPMNSOpt< T >::_OscCalcPMNSOpt ( )

Definition at line 14 of file OscCalcPMNSOpt.cxx.

       : fMixIdx(0), fDmIdx(0), fLRIdx(0)
   {
   }

template<typename T >

osc::_OscCalcPMNSOpt< T >::~_OscCalcPMNSOpt ( )

virtual

Definition at line 21 of file OscCalcPMNSOpt.cxx.

References osc::_OscCalcPMNSOpt< T >::fPMNSOpt, MECModelEnuComparisons::i, and it.

   {
     for(int i = 0; i < 2; ++i)
       for(auto it: fPMNSOpt[i])
         delete it.second.pmns;
   }

Member Function Documentation

template<typename T >

_IOscCalcAdjustable< T > * osc::_OscCalcPMNSOpt< T >::Copy ( ) const

overridevirtual

Implements osc::_IOscCalcAdjustable< T >.

Definition at line 30 of file OscCalcPMNSOpt.cxx.

References osc::_OscCalcPMNSOpt< T >::fPMNSOpt, and runNovaSAM::ret.

   {
     _OscCalcPMNSOpt<T>* ret = new _OscCalcPMNSOpt<T>(*this);
     // Raw pointers were blindly copied, so we'd be in trouble when the
     // destructors are called. More importantly, having two calculators sharing
     // one PMNS object will lead to both of them getting confused as to whether
     // they're up-to-date or not. Just clear out the cache entirely in the copy
     // and let it be repopulated.
     ret->fPMNSOpt[0].clear();
     ret->fPMNSOpt[1].clear();
     return ret;
   }

template<typename T>

virtual T osc::_IOscCalcAdjustable< T >::GetdCP ( ) const

inlinevirtualinherited

Reimplemented in osc::OscCalcCPT.

Definition at line 95 of file IOscCalc.h.

Referenced by osc::CopyParams(), osc::OscCalcCPT::GetdCP(), getHists_FNEX(), joint_fit_2017_make_fc_slice(), joint_fit_2017_make_fc_surf(), jointsensitivity(), ana::ResetOscCalcToTemplate(), ana::ResetSterileCalcToDefault(), ana::SaveTo(), ana::FitSinSq2Theta13::SetValue(), test_ana(), and test_nue2018_fitter().

95 { return fdCP ; }

osc::_IOscCalcAdjustable::fdCP

T fdCP

Definition: IOscCalc.h:119

template<typename T>

virtual T osc::_IOscCalcAdjustable< T >::GetDmsq21 ( ) const

inlinevirtualinherited

Reimplemented in osc::OscCalcCPT.

Definition at line 90 of file IOscCalc.h.

Referenced by osc::CopyParams(), ana::SolarConstraints::Dmsq21(), osc::OscCalcCPT::GetDmsq21(), ana::FitDmSq21::GetValue(), ana::ResetOscCalcToTemplate(), ana::ResetSterileCalcToDefault(), and ana::SaveTo().

90 { return fDmsq21 ; }

osc::_IOscCalcAdjustable::fDmsq21

T fDmsq21

Definition: IOscCalc.h:114

template<typename T>

virtual T osc::_IOscCalcAdjustable< T >::GetDmsq32 ( ) const

inlinevirtualinherited

Reimplemented in osc::OscCalcCPT.

Definition at line 91 of file IOscCalc.h.

Referenced by CAF_makeCAFSensitivities_for_FNEX(), cc(), ana::Dmsq32Constraint::ChiSq(), ana::AtmConstraint::ChiSq(), osc::CopyParams(), demoFitContours(), demoFitSlices(), getBestFit(), osc::OscCalcCPT::GetDmsq32(), getHists_FNEX(), ana::FitDmSq32::GetValue(), ana::FitDmSq32Scaled::GetValue(), ana::FitDmSq32NormalHierarchy::GetValue(), ana::FitDmSq32ScaledNH::GetValue(), ana::FitDmSq32ScaledIH::GetValue(), ana::FitDmSq32InvertedHierarchy::GetValue(), joint_fit_2017_contours(), joint_fit_2017_slices(), joint_fit_2018_contours(), joint_fit_2018_slices(), joint_fit_2019_contours(), joint_fit_2019_slices(), joint_fit_future_contour_univ(), jointsensitivity(), MakeCAFSensitivities_for_FNEX(), plot_3NDvsFD(), plot_NDvsFD_REW(), plot_NDvsFD_weights(), plot_nueFD_Signal_REWvsNOM(), ana::ResetOscCalcToTemplate(), ana::ResetSterileCalcToDefault(), run_joint_fit_2020_contours(), run_joint_fit_2020_slices(), ana::SaveTo(), sensitivity2018(), sensitivity2020(), test_nue2018_fitter(), test_nueextrapsyst(), Tutorial2019FitContours(), and Tutorial2019FitSlices().

91 { return fDmsq32 ; }

osc::_IOscCalcAdjustable::fDmsq32

T fDmsq32

Definition: IOscCalc.h:115

template<typename T>

virtual double osc::_IOscCalcAdjustable< T >::GetL ( ) const

inlinevirtualinherited

Reimplemented in osc::OscCalcCPT, and osc::OscCalcSterile.

Definition at line 88 of file IOscCalc.h.

Referenced by ana::T2KToyExperiment::ChiSq(), osc::CopyParams(), osc::OscCalcCPT::GetL(), ana::ResetSterileCalcToDefault(), ana::SaveTo(), and ana::T2KToyExperiment::SetTruthParams().

88 { return fL ; }

osc::_IOscCalcAdjustable::fL

double fL

Definition: IOscCalc.h:113

template<typename T>

TMD5* osc::_OscCalcPMNSOpt< T >::GetParamsHash ( ) const

inlineoverridevirtual

Use to check two calculators are in the same state.

Returns: Null means not implemented for this calculator

Reimplemented from osc::_IOscCalc< T >.

Definition at line 45 of file OscCalcPMNSOpt.h.

References osc::_IOscCalcAdjustable< T >::GetParamsHashDefault().

       {
         return _IOscCalcAdjustable<T>::GetParamsHashDefault("PMNSOpt");
       }

template<typename T >

TMD5 * osc::_IOscCalcAdjustable< T >::GetParamsHashDefault ( const std::string & txt ) const

protectedinherited

This is only a safe implementation if your calculator only depends on these eight parameters.

Parameters

txt	A string to uniquely identify your calculator class

Definition at line 47 of file IOscCalc.cxx.

References runNovaSAM::ret, and T.

Referenced by osc::_OscCalcPMNS< T >::GetParamsHash(), osc::_OscCalcPMNSOpt< T >::GetParamsHash(), and osc::analytic::_OscCalc< T >::GetParamsHash().

   {
     TMD5* ret = new TMD5;
     ret->Update((unsigned char*)txt.c_str(), txt.size());
     const int kNumParams = 8;
     T buf[kNumParams] = {fRho, fL, fDmsq21, fDmsq32,
                          fTh12, fTh13, fTh23, fdCP};
     ret->Update((unsigned char*)buf, sizeof(T)*kNumParams);
     ret->Final();
     return ret;
   }

template<typename T>

virtual double osc::_IOscCalcAdjustable< T >::GetRho ( ) const

inlinevirtualinherited

Reimplemented in osc::OscCalcCPT, and osc::OscCalcSterile.

Definition at line 89 of file IOscCalc.h.

Referenced by osc::CopyParams(), osc::OscCalcCPT::GetRho(), ana::ResetSterileCalcToDefault(), and ana::SaveTo().

89 { return fRho ; }

osc::_IOscCalcAdjustable::fRho

double fRho

Definition: IOscCalc.h:112

template<typename T>

virtual T osc::_IOscCalcAdjustable< T >::GetTh12 ( ) const

inlinevirtualinherited

Reimplemented in osc::OscCalcCPT.

Definition at line 92 of file IOscCalc.h.

Referenced by osc::CopyParams(), osc::OscCalcCPT::GetTh12(), ana::FitTanSqTheta12::GetValue(), ana::FitSinSq2Theta12::GetValue(), ana::ResetOscCalcToTemplate(), ana::ResetSterileCalcToDefault(), ana::SaveTo(), and ana::SolarConstraints::Th12().

92 { return fTh12 ; }

osc::_IOscCalcAdjustable::fTh12

T fTh12

Definition: IOscCalc.h:116

template<typename T>

virtual T osc::_IOscCalcAdjustable< T >::GetTh13 ( ) const

inlinevirtualinherited

Reimplemented in osc::OscCalcCPT.

Definition at line 93 of file IOscCalc.h.

Referenced by osc::CopyParams(), demoFitContours(), osc::OscCalcCPT::GetTh13(), ana::FitTheta13::GetValue(), ana::FitSinSq2Theta13::GetValue(), joint_fit_2018_contours(), joint_fit_future_bestfit_univ(), monoprob(), ana::ResetOscCalcToTemplate(), ana::ResetSterileCalcToDefault(), ana::SaveTo(), sensitivity2020(), test_ana(), and Tutorial2019FitContours().

93 { return fTh13 ; }

osc::_IOscCalcAdjustable::fTh13

T fTh13

Definition: IOscCalc.h:117

template<typename T>

virtual T osc::_IOscCalcAdjustable< T >::GetTh23 ( ) const

inlinevirtualinherited

template<typename T>

virtual void osc::_IOscCalcAdjustable< T >::InvalidateCache ( )

inlinevirtualinherited

Invalidate any caching used internally by the calculator.

Some calculators use a cache that can become stale in ways that the calculator may not know about (e.g., Stan var clearing). Default implementation does nothing.

Reimplemented in osc::_OscCalcDMP< T >.

Definition at line 102 of file IOscCalc.h.

102 {};

template<class T >

Eigen::Array< T, Eigen::Dynamic, 1 > osc::_IOscCalc< T >::P	(	int	flavBefore,
		int	flavAfter,
		const std::vector< double > &	E
	)

virtualinherited

Default implementation forwards to non-vector version using a simple loop. Override if your calculator has a more efficient implementation.

Reimplemented in osc::analytic::_OscCalc< T >, osc::OscCalcPMNSOptEigen, and osc::_OscCalcDMP< T >.

Definition at line 25 of file IOscCalc.cxx.

References MECModelEnuComparisons::i, osc::_IOscCalc< T >::P(), P, runNovaSAM::ret, and submit_hadd::u.

   {
     Eigen::Matrix<T,Eigen::Dynamic,1> ret(E.size());
     for(auto i = 0u; i < E.size(); i++) {
       ret(i) = this->P(flavBefore, flavAfter, E[i]);
     }
     return ret.array().isNaN().select(0, ret);
   }

template<class T >

Eigen::Array< T, Eigen::Dynamic, 1 > osc::_IOscCalc< T >::P	(	int	flavBefore,
		int	flavAfter,
		const Eigen::ArrayXd &	E
	)

virtualinherited

Default implementation forawrds to vector<double> version. Override if your calculator has a more efficient implementation.

Reimplemented in osc::analytic::_OscCalc< T >.

Definition at line 35 of file IOscCalc.cxx.

References P.

   {
     return P(flavBefore, flavAfter, std::vector<double>(&E[0], &E[0]+E.size()));
   }

template<typename T >

T osc::_OscCalcPMNSOpt< T >::P	(	int	flavBefore,
		int	flavAfter,
		double	E
	)

overridevirtual

E in GeV; flavors as PDG codes (so, neg==>antinu)

Implements osc::_IOscCalc< T >.

Definition at line 45 of file OscCalcPMNSOpt.cxx.

Referenced by Toy_analyses().

   {
     // Normal usage of a calculator in a fit is to configure one set of
     // oscillation parameters and then calculate probabilities for all flavours
     // and for a range of energies. The underlying PMNSOpt object has to redo
     // its expensive diagonalization every time the energy changes, but not
     // when simply switching to a different flavour. Unfortunately callers in
     // CAFAna wind up making the calls in the wrong order to take advantage of
     // this. So, we maintain a map of calculators for different energies (also
     // for neutrinos and antineutrinos) and select the right one if available,
     // on the hope that it's already configured to do what we're asked, since
     // it might already have done a different flavour. This works because
     // callers often pass the same set of energies (bin centers) over and over
     // again.
 
     // If the caches get too large clear them. The caller is probably giving
     // the precise energy of each event?
     for(int i = 0; i < 2; ++i){
       if(fPMNSOpt[i].size() > 10000){
         for(auto it: fPMNSOpt[i]) delete it.second.pmns;
         fPMNSOpt[i].clear();
       }
     }
 
     const int anti = (flavBefore > 0) ? +1 : -1;
     assert(flavAfter/anti > 0);
 
     Val_t& calc = fPMNSOpt[(1+anti)/2][E];
     if(!calc.pmns) calc.pmns = new _PMNSOpt<T>;
 
     int i = -1, j = -1;
     if(abs(flavBefore) == 12) i = 0;
     if(abs(flavBefore) == 14) i = 1;
     if(abs(flavBefore) == 16) i = 2;
     if(abs(flavAfter) == 12) j = 0;
     if(abs(flavAfter) == 14) j = 1;
     if(abs(flavAfter) == 16) j = 2;
     assert(i >= 0 && j >= 0);
 
     const bool dirty = (fMixIdx > calc.mixIdx ||
                         fDmIdx > calc.dmIdx ||
                         fLRIdx > calc.lrIdx);
 
     // If the calculator has out of date parameters update them
     if(fMixIdx > calc.mixIdx){
       calc.pmns->SetMix(this->fTh12, this->fTh23, this->fTh13, this->fdCP);
       calc.mixIdx = fMixIdx;
     }
     if(fDmIdx > calc.dmIdx){
       calc.pmns->SetDeltaMsqrs(this->fDmsq21, this->fDmsq32);
       calc.dmIdx = fDmIdx;
     }
 
     if(dirty){
       // Cache results for all nine flavour combinations
       for(int ii = 0; ii < 3; ++ii){
         calc.pmns->ResetToFlavour(ii);
         // Assume Z/A=0.5
         const double Ne = this->fRho/2;
         calc.pmns->PropMatter(this->fL, E, Ne, anti);
         for(int jj = 0; jj < 3; ++jj){
           calc.P[ii][jj] = calc.pmns->P(jj);
         }
       }
 
       calc.lrIdx = fLRIdx;
     }
 
     return calc.P[i][j];
   }