TwoFlavourFitter.cxx File Reference

#include "TwoFlavourFitter.h"
#include "TemplateReader.h"
#include "TemplateWriter.h"
#include "TemplateCalculator.h"
#include "TemplateGrid.h"
#include "ParametersUtil.h"
#include "Configuration.h"
#include "Likelihood.h"
#include "TString.h"
#include "TRandom.h"
#include "TMarker.h"
#include "TGraph.h"
#include "TH2.h"
#include "TH1.h"
#include "TFile.h"
#include "TTree.h"
#include "TDirectory.h"
#include <iostream>
#include <cmath>
#include <ctime>
#include <cassert>

Go to the source code of this file.

Functions

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

Variables

static TwoFlavourFitterfgTwoFlavourFitter = 0

Function Documentation

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

Definition at line 32 of file TwoFlavourFitter.cxx.

References OscFit::TwoFlavourFitter::fAtmosChgCV, OscFit::TwoFlavourFitter::fAtmosChgRock, OscFit::TwoFlavourFitter::fAtmosNCBkg, OscFit::TwoFlavourFitter::fAtmosNormCV, OscFit::TwoFlavourFitter::fAtmosNormRock, OscFit::TwoFlavourFitter::fAtmosNue, OscFit::TwoFlavourFitter::fAtmosShwEn, OscFit::TwoFlavourFitter::fAtmosSpecNuBarCV, OscFit::TwoFlavourFitter::fAtmosSpecNuBarRock, OscFit::TwoFlavourFitter::fAtmosSpecNuCV, OscFit::TwoFlavourFitter::fAtmosSpecNuRock, OscFit::TwoFlavourFitter::fAtmosTrkEn, OscFit::TwoFlavourFitter::fAtmosTrkEnExit, OscFit::TwoFlavourFitter::fAtmosZenith, OscFit::TwoFlavourFitter::fBeamNCBkg, OscFit::TwoFlavourFitter::fBeamNorm, OscFit::TwoFlavourFitter::fBeamShwEn, OscFit::TwoFlavourFitter::fBeamTrkEn, OscFit::TwoFlavourFitter::fBeamTrkEnExit, OscFit::TwoFlavourFitter::fFittingAtmosChgCV, OscFit::TwoFlavourFitter::fFittingAtmosChgRock, OscFit::TwoFlavourFitter::fFittingAtmosNCBkg, OscFit::TwoFlavourFitter::fFittingAtmosNormCV, OscFit::TwoFlavourFitter::fFittingAtmosNormRock, OscFit::TwoFlavourFitter::fFittingAtmosNue, OscFit::TwoFlavourFitter::fFittingAtmosShwEn, OscFit::TwoFlavourFitter::fFittingAtmosSpecNuBarCV, OscFit::TwoFlavourFitter::fFittingAtmosSpecNuBarRock, OscFit::TwoFlavourFitter::fFittingAtmosSpecNuCV, OscFit::TwoFlavourFitter::fFittingAtmosSpecNuRock, OscFit::TwoFlavourFitter::fFittingAtmosTrkEn, OscFit::TwoFlavourFitter::fFittingAtmosTrkEnExit, OscFit::TwoFlavourFitter::fFittingAtmosZenith, OscFit::TwoFlavourFitter::fFittingBeamNCBkg, OscFit::TwoFlavourFitter::fFittingBeamNorm, OscFit::TwoFlavourFitter::fFittingBeamShwEn, OscFit::TwoFlavourFitter::fFittingBeamTrkEn, OscFit::TwoFlavourFitter::fFittingBeamTrkEnExit, OscFit::TwoFlavourFitter::fFittingDmsq, OscFit::TwoFlavourFitter::fFittingSinsq, OscFit::TwoFlavourFitter::fOscDmsq, OscFit::TwoFlavourFitter::fOscSinsq, OscFit::TwoFlavourFitter::get_counter(), OscFit::TwoFlavourFitter::GetLnL(), and OscFit::TwoFlavourFitter::itr().

Referenced by OscFit::TwoFlavourFitter::RunFitWithMinuit().

00033 {
00034   Bool_t printDebug = 0;
00035 
00036   Double_t osc_dmsq              = fgTwoFlavourFitter->fOscDmsq;
00037   Double_t osc_sinsq             = fgTwoFlavourFitter->fOscSinsq;
00038   Double_t beam_norm             = fgTwoFlavourFitter->fBeamNorm; 
00039   Double_t beam_ncbkg            = fgTwoFlavourFitter->fBeamNCBkg; 
00040   Double_t beam_shwen            = fgTwoFlavourFitter->fBeamShwEn;
00041   Double_t beam_trken            = fgTwoFlavourFitter->fBeamTrkEn;
00042   Double_t beam_trken_exit       = fgTwoFlavourFitter->fBeamTrkEnExit;
00043   Double_t atmos_norm_cv         = fgTwoFlavourFitter->fAtmosNormCV;
00044   Double_t atmos_norm_rock       = fgTwoFlavourFitter->fAtmosNormRock; 
00045   Double_t atmos_chg_cv          = fgTwoFlavourFitter->fAtmosChgCV;
00046   Double_t atmos_chg_rock        = fgTwoFlavourFitter->fAtmosChgRock; 
00047   Double_t atmos_nue             = fgTwoFlavourFitter->fAtmosNue;     
00048   Double_t atmos_ncbkg           = fgTwoFlavourFitter->fAtmosNCBkg;
00049   Double_t atmos_spec_cv_nu      = fgTwoFlavourFitter->fAtmosSpecNuCV;  
00050   Double_t atmos_spec_cv_nubar   = fgTwoFlavourFitter->fAtmosSpecNuBarCV;
00051   Double_t atmos_spec_rock_nu    = fgTwoFlavourFitter->fAtmosSpecNuRock;  
00052   Double_t atmos_spec_rock_nubar = fgTwoFlavourFitter->fAtmosSpecNuBarRock;
00053   Double_t atmos_zenith          = fgTwoFlavourFitter->fAtmosZenith;
00054   Double_t atmos_trken           = fgTwoFlavourFitter->fAtmosTrkEn;
00055   Double_t atmos_trken_exit      = fgTwoFlavourFitter->fAtmosTrkEnExit;   
00056   Double_t atmos_shwen           = fgTwoFlavourFitter->fAtmosShwEn; 
00057 
00058   if( fgTwoFlavourFitter->fFittingDmsq )               osc_dmsq              = par[0];
00059   if( fgTwoFlavourFitter->fFittingSinsq )              osc_sinsq             = par[1];
00060   if( fgTwoFlavourFitter->fFittingBeamNorm )           beam_norm             = par[2];
00061   if( fgTwoFlavourFitter->fFittingBeamNCBkg )          beam_ncbkg            = par[3];
00062   if( fgTwoFlavourFitter->fFittingBeamShwEn )          beam_shwen            = par[4];
00063   if( fgTwoFlavourFitter->fFittingBeamTrkEn )          beam_trken            = par[5];
00064   if( fgTwoFlavourFitter->fFittingBeamTrkEnExit )      beam_trken_exit       = par[6];
00065   if( fgTwoFlavourFitter->fFittingAtmosNormCV )        atmos_norm_cv         = par[7];
00066   if( fgTwoFlavourFitter->fFittingAtmosNormRock )      atmos_norm_rock       = par[8];
00067   if( fgTwoFlavourFitter->fFittingAtmosChgCV )         atmos_chg_cv          = par[9];
00068   if( fgTwoFlavourFitter->fFittingAtmosChgRock )       atmos_chg_rock        = par[10];
00069   if( fgTwoFlavourFitter->fFittingAtmosNue )           atmos_nue             = par[11];
00070   if( fgTwoFlavourFitter->fFittingAtmosNCBkg )         atmos_ncbkg           = par[12];
00071   if( fgTwoFlavourFitter->fFittingAtmosSpecNuCV )      atmos_spec_cv_nu      = par[13];
00072   if( fgTwoFlavourFitter->fFittingAtmosSpecNuBarCV )   atmos_spec_cv_nubar   = par[14];
00073   if( fgTwoFlavourFitter->fFittingAtmosSpecNuRock )    atmos_spec_rock_nu    = par[15];
00074   if( fgTwoFlavourFitter->fFittingAtmosSpecNuBarRock ) atmos_spec_rock_nubar = par[16];
00075   if( fgTwoFlavourFitter->fFittingAtmosZenith )        atmos_zenith          = par[17];
00076   if( fgTwoFlavourFitter->fFittingAtmosTrkEn )         atmos_trken           = par[18];
00077   if( fgTwoFlavourFitter->fFittingAtmosTrkEnExit )     atmos_trken_exit      = par[19];
00078   if( fgTwoFlavourFitter->fFittingAtmosShwEn )         atmos_shwen           = par[20];  
00079 
00080   fgTwoFlavourFitter->GetLnL( osc_dmsq,           osc_sinsq,
00081                               beam_norm,          beam_ncbkg,
00082                               beam_shwen,         beam_trken,
00083                               beam_trken_exit,
00084                               atmos_norm_cv,      atmos_norm_rock, 
00085                               atmos_chg_cv,       atmos_chg_rock, 
00086                               atmos_nue,          atmos_ncbkg,
00087                               atmos_spec_cv_nu,   atmos_spec_cv_nubar,
00088                               atmos_spec_rock_nu, atmos_spec_rock_nubar,
00089                               atmos_zenith,
00090                               atmos_trken,        atmos_trken_exit,
00091                               atmos_shwen,
00092                               f );
00093 
00094   if( printDebug ){
00095     std::cout << " --- [" << fgTwoFlavourFitter->get_counter() << "]: " << std::endl;
00096     std::cout << "     osc_dmsq=" << osc_dmsq << ", osc_sinsq=" << osc_sinsq << std::endl;
00097     std::cout << "     beam_norm=" << beam_norm << ", beam_ncbkg=" << beam_ncbkg << ", beam_shwen=" << beam_shwen << ", beam_trken=" << beam_trken << ", beam_trken_exit=" << beam_trken_exit << std::endl;
00098     std::cout << "     atmos_norm_cv=" << atmos_norm_cv << ", atmos_norm_rock=" << atmos_norm_rock << ", atmos_chg_cv=" << atmos_chg_cv << ", atmos_chg_rock=" <<atmos_chg_rock << ", atmos_nue=" << atmos_nue << ", atmos_ncbkg=" << atmos_ncbkg << std::endl;
00099     std::cout << "     atmos_spec_cv_nu=" << atmos_spec_cv_nu << " atmos_spec_cv_nubar=" << atmos_spec_cv_nubar << ", atmos_spec_rock_nu=" << atmos_spec_rock_nu << ", atmos_spec_rock_nubar=" << atmos_spec_rock_nubar << ", atmos_zenith=" << atmos_zenith << std::endl;
00100     std::cout << "     atmos_trken=" << atmos_trken << ", atmos_trken_exit=" << atmos_trken_exit << ", atmos_shwen=" << atmos_shwen << std::endl;
00101     std::cout << "     f=" << f << std::endl;
00102   }
00103 
00104   fgTwoFlavourFitter->itr();
00105 
00106   return;
00107 }


Variable Documentation

Definition at line 30 of file TwoFlavourFitter.cxx.


Generated on 3 Oct 2018 for loon by  doxygen 1.6.1