Main Page | Modules | Namespace List | Class Hierarchy | Alphabetical List | Class List | Directories | File List | Namespace Members | Class Members | File Members

NueRWHelpers Namespace Reference


Functions

float Oscillate (NtpMCTruth *mcth, float L, float dm2, float theta23, float UE32)
float Oscillate (ANtpTruthInfoBeam *ib, float L, float dm2, float theta23, float UE32)
float Oscillate (int nuFlavor, int nonOscNuFlavor, float Energy, float L, float dm2, float theta23, float U)
float OscillateMatter (int nuFlavor, int nonOscNuFlavor, float Energy, float L, float dm2, float theta23, float UE32, float delta=0, Int_t hierarchy=1)
float OscillateMatter (NtpMCTruth *mcth, float L, float dm2, float theta23, float UE32, float delta=0, int hierarchy=1)
float OscillateMatter (ANtpTruthInfoBeam *ib, float L, float dm2, float theta23, float UE32, float delta=0, int hierarchy=1)

Function Documentation

float NueRWHelpers::Oscillate (  int  nuFlavor,
int  nonOscNuFlavor,
float  Energy,
float  L,
float  dm2,
float  theta23,
float  U
) 
 

Definition at line 22 of file NueRWHelpers.cxx.

References abs().

Referenced by MadCluAnalysis::FOM(), and MadCluAnalysis::PIDPlot().

00024 {
00025    float oscterm = TMath::Sin(1.27*dm2*L/Energy);      
00026 
00027    float pmt=pow((1-UE32)*oscterm*TMath::Sin(2*theta23),2);
00028    float pme=pow(TMath::Sin(theta23),2)*4.*UE32*(1-UE32)*pow(oscterm,2);
00029    float pmm=1.-pmt-pme;
00030    
00031    float pet=4*(1-UE32)*UE32*pow(TMath::Cos(theta23)*oscterm,2);
00032    float pem=pow(TMath::Sin(theta23),2)*4.*UE32*(1-UE32)*pow(oscterm,2);
00033    float pee=1.-pet-pem;
00034 
00035 
00036    if(abs(nonOscNuFlavor)==14){
00037       if(abs(nuFlavor)==12){
00038          return pme;
00039       }
00040       else if(abs(nuFlavor)==14){
00041          return pmm;
00042       }
00043       else if(abs(nuFlavor)==16){
00044          return pmt;
00045       }
00046    }
00047    else if(abs(nonOscNuFlavor)==12){
00048       if(abs(nuFlavor)==12){
00049          return pee;
00050       }
00051       else if(abs(nuFlavor)==14){
00052          return pem;
00053       }
00054       else if(abs(nuFlavor)==16){
00055          return pet;
00056       }
00057    }
00058    else{
00059      std::cout<<"I don't know what to do with "<<nonOscNuFlavor
00060          <<" "<<nuFlavor<<" "<<pee<<std::endl;
00061    }
00062    return 0.;
00063 }

float NueRWHelpers::Oscillate (  ANtpTruthInfoBeam *  ib,
float  L,
float  dm2,
float  theta23,
float  UE32
) 
 

Definition at line 15 of file NueRWHelpers.cxx.

References ANtpTruthInfoBeam::nonOscNuFlavor, ANtpTruthInfo::nuEnergy, ANtpTruthInfo::nuFlavor, and Oscillate().

00017 {
00018   return NueRWHelpers::Oscillate(ib->nuFlavor, ib->nonOscNuFlavor,ib->nuEnergy,
00019                                    L, dm2, theta23, UE32);
00020 }

float NueRWHelpers::Oscillate (  NtpMCTruth *  mcth,
float  L,
float  dm2,
float  theta23,
float  UE32
) 
 

Definition at line 8 of file NueRWHelpers.cxx.

References NtpMCTruth::inu, NtpMCTruth::inunoosc, and NtpMCTruth::p4neu.

Referenced by ANtpTruthInfoBeamAna::GetOscProb(), Oscillate(), NuePrint::PrintValues(), and NueReweight::Reco().

00010 {
00011     return NueRWHelpers::Oscillate(mcth->inu, mcth->inunoosc, mcth->p4neu[3],
00012                                    L, dm2, theta23, UE32);
00013 }

float NueRWHelpers::OscillateMatter (  ANtpTruthInfoBeam *  ib,
float  L,
float  dm2,
float  theta23,
float  UE32,
float  delta = 0,
int  hierarchy = 1
) 
 

Definition at line 117 of file NueRWHelpers.cxx.

References ANtpTruthInfoBeam::nonOscNuFlavor, ANtpTruthInfo::nuEnergy, ANtpTruthInfo::nuFlavor, and OscillateMatter().

00120 {
00121   return NueRWHelpers::OscillateMatter(ib->nuFlavor, ib->nonOscNuFlavor,
00122                                        ib->nuEnergy,
00123                                        L, dm2, theta23, UE32,
00124                                        delta,hierarchy);
00125 }

float NueRWHelpers::OscillateMatter (  NtpMCTruth *  mcth,
float  L,
float  dm2,
float  theta23,
float  UE32,
float  delta = 0,
int  hierarchy = 1
) 
 

Definition at line 107 of file NueRWHelpers.cxx.

References NtpMCTruth::inu, NtpMCTruth::inunoosc, OscillateMatter(), and NtpMCTruth::p4neu.

00110 {
00111   return NueRWHelpers::OscillateMatter(mcth->inu, mcth->inunoosc, 
00112                                        mcth->p4neu[3],
00113                                        L, dm2, theta23, UE32,
00114                                        delta,hierarchy);
00115 }

float OscillateMatter (  int  nuFlavor,
int  nonOscNuFlavor,
float  Energy,
float  L,
float  dm2,
float  theta23,
float  UE32,
float  delta = 0,
Int_t  hierarchy = 1
) 
 

Referenced by OscillateMatter().

Generated on Thu Nov 1 15:58:51 2007 for loon by  doxygen 1.3.9.1