00001
00002
00003
00004
00005
00006
00007
00008
00009 #include "Riostream.h"
00010
00011 #include "NuMuToNuE_Oscillation.h"
00012 #include "RooAbsReal.h"
00013 #include "RooAbsCategory.h"
00014 #include <math.h>
00015 #include "TMath.h"
00016
00017
00018
00019 NuMuToNuE_Oscillation::NuMuToNuE_Oscillation(const char *name, const char *title,
00020 RooAbsReal& _L,
00021 RooAbsReal& _E,
00022 RooAbsReal& _deltaMSq) :
00023 RooAbsPdf(name,title),
00024 L("L","L",this,_L),
00025 E("E","E",this,_E),
00026 deltaMSq("deltaMSq","deltaMSq",this,_deltaMSq)
00027 {
00028 }
00029
00030
00031 NuMuToNuE_Oscillation::NuMuToNuE_Oscillation(const NuMuToNuE_Oscillation& other, const char* name) :
00032 RooAbsPdf(other,name),
00033 L("L",this,other.L),
00034 E("E",this,other.E),
00035 deltaMSq("deltaMSq",this,other.deltaMSq)
00036 {
00037 }
00038
00039
00040
00041 Double_t NuMuToNuE_Oscillation::evaluate() const
00042 {
00043
00044 return pow(sin(1.27*deltaMSq*L/E),2) ;
00045 }
00046