AtNuOscFit.cxx File Reference

#include "AtNuOscFit.h"
#include <iostream>
#include <cmath>
#include <cassert>

Go to the source code of this file.

Functions

static void fcn (Int_t &, Double_t *, Double_t &f, Double_t *par, Int_t)

Variables

static AtNuOscFitfAtNuOscFit = 0

Function Documentation

static void fcn ( Int_t &  ,
Double_t *  ,
Double_t &  f,
Double_t *  par,
Int_t   
) [static]

Definition at line 11 of file AtNuOscFit.cxx.

Referenced by FitNdNonlinQuad(), AtNuOscFit::GetLikelihoodSyst(), PANAnalysis::SetOscillationFunction(), and MadAnalysis::SetOscillationFunction().

00012 {
00013   // input parameters
00014   Double_t norm      = fAtNuOscFit->fFitNorm;       
00015   Double_t normupmu  = fAtNuOscFit->fFitNormUpmu;
00016   Double_t elmu      = fAtNuOscFit->fFitElMu;
00017   Double_t ccnc      = fAtNuOscFit->fFitCCNC;      
00018   Double_t updn      = fAtNuOscFit->fFitUpDn;      
00019   Double_t chg       = fAtNuOscFit->fFitChg;       
00020   Double_t chgupmu   = fAtNuOscFit->fFitChgUpmu;    
00021   Double_t specnu    = fAtNuOscFit->fFitSpecNu;    
00022   Double_t specnubar = fAtNuOscFit->fFitSpecNuBar; 
00023   Double_t emurange  = fAtNuOscFit->fFitEmuRange;  
00024   Double_t emucurve  = fAtNuOscFit->fFitEmuCurve;  
00025   Double_t eshw      = fAtNuOscFit->fFitEshw;
00026   Double_t fsnu      = fAtNuOscFit->fFitSinNu;
00027   Double_t fsnubar   = fAtNuOscFit->fFitSinNuBar;
00028   Double_t dmnu      = fAtNuOscFit->fFitDmsqNu;
00029   Double_t dmnubar   = fAtNuOscFit->fFitDmsqNuBar;
00030 
00031   if(fAtNuOscFit->fUsingNorm)      norm      = par[0];    
00032   if(fAtNuOscFit->fUsingNormUpmu)  normupmu  = par[1];
00033   if(fAtNuOscFit->fUsingElMu)      elmu      = par[2];   
00034   if(fAtNuOscFit->fUsingCCNC)      ccnc      = par[3];
00035   if(fAtNuOscFit->fUsingUpDn)      updn      = par[4];    
00036   if(fAtNuOscFit->fUsingChg)       chg       = par[5];     
00037   if(fAtNuOscFit->fUsingChgUpmu)   chgupmu   = par[6];   
00038   if(fAtNuOscFit->fUsingSpecNu)    specnu    = par[7];  
00039   if(fAtNuOscFit->fUsingSpecNuBar) specnubar = par[8];  
00040   if(fAtNuOscFit->fUsingEmuRange)  emurange  = par[9];     
00041   if(fAtNuOscFit->fUsingEmuCurve)  emucurve  = par[10];
00042   if(fAtNuOscFit->fUsingEshw)      eshw      = par[11];    
00043   if(fAtNuOscFit->fUsingSinNu)     fsnu      = par[12];
00044   if(fAtNuOscFit->fUsingSinNuBar)  fsnubar   = par[13];
00045   if(fAtNuOscFit->fUsingDmsqNu)    dmnu      = par[14];
00046   if(fAtNuOscFit->fUsingDmsqNuBar) dmnubar   = par[15];
00047 
00048   if( fAtNuOscFit->fDebugMinuit ){
00049     std::cout << "  [0]: "  << fAtNuOscFit->fUsingNorm      << " " << fAtNuOscFit->fFitNorm      << " " << norm      << std::endl
00050               << "  [1]: "  << fAtNuOscFit->fUsingNormUpmu  << " " << fAtNuOscFit->fFitNormUpmu  << " " << normupmu  << std::endl
00051               << "  [2]: "  << fAtNuOscFit->fUsingElMu      << " " << fAtNuOscFit->fFitElMu      << " " << elmu      << std::endl
00052               << "  [3]: "  << fAtNuOscFit->fUsingCCNC      << " " << fAtNuOscFit->fFitCCNC      << " " << ccnc      << std::endl
00053               << "  [4]: "  << fAtNuOscFit->fUsingUpDn      << " " << fAtNuOscFit->fFitUpDn      << " " << updn      << std::endl
00054               << "  [5]: "  << fAtNuOscFit->fUsingChg       << " " << fAtNuOscFit->fFitChg       << " " << chg       << std::endl
00055               << "  [6]: "  << fAtNuOscFit->fUsingChgUpmu   << " " << fAtNuOscFit->fFitChgUpmu   << " " << chgupmu   << std::endl
00056               << "  [7]: "  << fAtNuOscFit->fUsingSpecNu    << " " << fAtNuOscFit->fFitSpecNu    << " " << specnu    << std::endl
00057               << "  [8]: "  << fAtNuOscFit->fUsingSpecNuBar << " " << fAtNuOscFit->fFitSpecNuBar << " " << specnubar << std::endl
00058               << "  [9]: "  << fAtNuOscFit->fUsingEmuRange  << " " << fAtNuOscFit->fFitEmuRange  << " " << emurange  << std::endl
00059               << "  [10]: " << fAtNuOscFit->fUsingEmuCurve  << " " << fAtNuOscFit->fFitEmuCurve  << " " << emucurve  << std::endl
00060               << "  [11]: " << fAtNuOscFit->fUsingEshw      << " " << fAtNuOscFit->fFitEshw      << " " << eshw      << std::endl
00061               << "  [12]: " << fAtNuOscFit->fUsingSinNu     << " " << fAtNuOscFit->fFitSinNu     << " " << fsnu      << std::endl
00062               << "  [13]: " << fAtNuOscFit->fUsingSinNuBar  << " " << fAtNuOscFit->fFitSinNuBar  << " " << fsnubar   << std::endl
00063               << "  [14]: " << fAtNuOscFit->fUsingDmsqNu    << " " << fAtNuOscFit->fFitDmsqNu    << " " << dmnu      << std::endl
00064               << "  [15]: " << fAtNuOscFit->fUsingDmsqNuBar << " " << fAtNuOscFit->fFitDmsqNuBar << " " << dmnubar   << std::endl;
00065   }
00066 
00067   // calculate likelihood
00068   fAtNuOscFit->GetLikelihood(norm,normupmu,
00069                              elmu,ccnc,updn,chg,chgupmu,
00070                              specnu,specnubar,
00071                              emurange,emucurve,eshw,
00072                              fsnu,fsnubar, 
00073                              dmnu,dmnubar, f);
00074 
00075   return;
00076 }


Variable Documentation

AtNuOscFit* fAtNuOscFit = 0 [static]

Definition at line 9 of file AtNuOscFit.cxx.

Referenced by AtNuOscFit::Instance().


Generated on 19 Jan 2018 for loon by  doxygen 1.6.1