genie::flux::GFlavorMap Class Reference

GENIE interface for flavor modification. More...

#include <GFlavorMap.h>

Inheritance diagram for genie::flux::GFlavorMap:

Collaboration diagram for genie::flux::GFlavorMap:

Public Member Functions
	GFlavorMap ()
	~GFlavorMap ()
void	Config (std::string config)
double	Probability (int pdg_initial, int pdg_final, double energy, double dist)
void	PrintConfig (bool verbose=true)
	provide a means of printing the configuration More...
Public Member Functions inherited from genie::flux::GFlavorMixerI
	GFlavorMixerI ()
virtual	~GFlavorMixerI ()

Private Member Functions
void	ParseMapString (std::string config)
void	ParseFixedfracString (std::string config)
int	PDG2Indx (int pdg)
int	Indx2PDG (int indx)
const char *	IndxName (int indx)
const char *	NuName (int pdg)

Private Attributes
double	fProb [7][7]

Detailed Description

GENIE interface for flavor modification.

     Concrete instance of GFlavorMixerI that maps from
     one flavor to another independent of energy or distance.
     Users specify the transition probability from one flavor
     to any of the PDG codes { 0, 12, 14, 16, -12, -14, -16 }
     where 0 represents the complete disappearance (decay, sterile, ...).

     Probability is expected to be normalized (that is, the sum
     of all possible outcomes, including 0, must be 1).

     Supported config string formats:
     1)  " swap pdg1:pdg2  pdg3:pdg4 "
         Map all neutrinos of flavor "pdg1" to "pdg2",  "pdg3" to "pdg4"
         Use PDG values { 0, 12, 14, 16, -12, -14, -16 }
           for { sterile, nu_e, nu_mu, nu_tau, nu_e_bar, ...}
         Unnamed initial flavors are left unmodified.
         Use numeric values with spaces only between pairs.
         Must start with the literal "swap"
         (FMWK note:  param must be surrounded by explicit "'s)
     2)  " fixedfrac {pdg1:f0,f12,f14,f16,f-12,f-14,f-16} ..."
         For each group delineated by {}'s the map the "pdg1"
         by each pdg by the fraction given [0-1, sum=1].
         So {12:0.5,0.5,0,0,0,0,0} means nu_e => 50/50% nu_e/nu_mu.
         Each list *must* have an int + 7 fractions.

Author

Robert Hatcher <rhatcher fnal.gov> Fermi National Accelerator Laboratory

Created:

2010-10-31

License:

Copyright (c) 2003-2024, The GENIE Collaboration for the full text of the license visit http://copyright.genie-mc.org

Definition at line 51 of file GFlavorMap.h.

Constructor & Destructor Documentation

genie::flux::GFlavorMap::GFlavorMap

(

)

Definition at line 31 of file GFlavorMap.cxx.

References fProb.

{
  // Initialize with identity matrix
  size_t jin, jout;
  for (jin=0; jin<7; ++jin) {
    for (jout=0; jout<7; ++jout ) {
      fProb[jin][jout] = ( (jin==jout) ? 1. : 0. );
    }
  }
}

genie::flux::GFlavorMap::~GFlavorMap

(

)

Definition at line 42 of file GFlavorMap.cxx.

{ ; }

Member Function Documentation

void genie::flux::GFlavorMap::Config ( std::string  config )

virtual

each schema must take a string that configures it it is up to the individual model to parse said string and extract parameters (e.g. sin2th23, deltam12, etc)

Implements genie::flux::GFlavorMixerI.

Definition at line 45 of file GFlavorMap.cxx.

References LOG_BEGIN, LOG_END, ParseFixedfracString(), ParseMapString(), pINFO, pWARN, and genie::utils::str::TrimSpaces().

{
  std::string config = genie::utils::str::TrimSpaces(configIn);
  LOG_BEGIN("FluxBlender", pINFO)
    << "GFlavorMap::Config \"" << config << "\"" << LOG_END;

  if ( config.find("swap") == 0 ||
       config.find("map")  == 0 ||
       config.find("genie::flux::GFlavorMap") == 0 ) {
    ParseMapString(config);
  } else if ( config.find("fixedfrac") == 0 ) {
    ParseFixedfracString(config);
  } else {
    LOG_BEGIN("FluxBlender", pWARN)
      << "GFlavorMap::Config don't know how to parse \""
      << config << "\"" << LOG_END;
    LOG_BEGIN("FluxBlender", pWARN)
      << " ... will attempt \"map\" strategy" << LOG_END;

  }

}

int genie::flux::GFlavorMap::Indx2PDG ( int  indx )

inlineprivate

Definition at line 117 of file GFlavorMap.h.

Referenced by PrintConfig().

{
  switch ( indx ) {
  case  1: return  12; break;
  case  2: return  14; break;
  case  3: return  16; break;
  case  4: return -12; break;
  case  5: return -14; break;
  case  6: return -16; break;
  default: return   0; break;
  }
  return 0;
}

const char * genie::flux::GFlavorMap::IndxName ( int  indx )

private

Definition at line 180 of file GFlavorMap.cxx.

Referenced by NuName(), and PrintConfig().

{
  static const char* name[] = { "sterile",
                                "nu_e", "nu_mu", "nu_tau",
                                "nu_e_bar", "nu_mu_bar", "nu_tau_bar" };
  return name[indx];

}

const char* genie::flux::GFlavorMap::NuName ( int  pdg )

inlineprivate

Definition at line 85 of file GFlavorMap.h.

References IndxName(), and PDG2Indx().

{ return IndxName(PDG2Indx(pdg)); }

void genie::flux::GFlavorMap::ParseFixedfracString ( std::string  config )

private

Definition at line 99 of file GFlavorMap.cxx.

References fProb, LOG_BEGIN, LOG_END, PDG2Indx(), pINFO, pWARN, genie::utils::str::Split(), and genie::utils::str::TrimSpaces().

Referenced by Config().

{
  LOG_BEGIN("FluxBlender", pINFO)
    << "GFlavorMap::ParseFixedFracString \"" << config << "\"" << LOG_END;
  vector<string> tokens = genie::utils::str::Split(config,"{}");
  for (unsigned int jtok = 0; jtok< tokens.size(); ++jtok ) {
    string tok1 = genie::utils::str::TrimSpaces(tokens[jtok]);
    if ( tok1 == "" ) continue;
    if ( tok1 == "fixedfrac" ) continue;
    // should have the form pdg:f0,f12,f14,f16,f-12,f-14,f-16
    vector<string> pair = genie::utils::str::Split(tok1,":");
    if ( pair.size() != 2 ) {
      LOG_BEGIN("FluxBlender", pWARN)
        << "could not parse \"" << tok1 << "\" split size=" << pair.size()
        << LOG_END;
      continue;
    }
    int pdg_in   = strtol(pair[0].c_str(),NULL,0);
    int indx_in  = PDG2Indx(pdg_in);
    vector<string> fracs = genie::utils::str::Split(pair[1],",");
    if ( fracs.size() != 7 ) {
      LOG_BEGIN("FluxBlender", pWARN)
        << "could not parse frac list \"" << pair[1] << "\" split size=" << fracs.size()
        << LOG_END;
      continue;
    }
    // set each value
    double psum = 0;
    for (int indx_out = 0; indx_out < 7; ++indx_out ) {
      double p = strtod(fracs[indx_out].c_str(),NULL);
      psum += p;
      fProb[indx_in][indx_out] = p;
    }
    if ( psum > 0 ) {
      // normalize to 1.0
      for (int indx_out = 0; indx_out < 7; ++indx_out )
        fProb[indx_in][indx_out] /= psum;
    }
  }
}

void genie::flux::GFlavorMap::ParseMapString ( std::string  config )

private

Definition at line 69 of file GFlavorMap.cxx.

References fProb, LOG_BEGIN, LOG_END, PDG2Indx(), pINFO, pWARN, and genie::utils::str::Split().

Referenced by Config().

{
  LOG_BEGIN("FluxBlender", pINFO)
    << "GFlavorMap::ParseMapString \"" << config << "\"" << LOG_END;
  vector<string> tokens = genie::utils::str::Split(config," ");
  for (unsigned int jtok = 0; jtok < tokens.size(); ++jtok ) {
    string tok1 = tokens[jtok];
    if ( tok1 == "" ) continue;
    if ( tok1 == "swap" || tok1 == "map" ) continue;
    if ( tok1 == "genie::flux::GFlavorMap" ) continue;
    // should have the form  <int>:<int>
    vector<string> pair = genie::utils::str::Split(tok1,":");
    if ( pair.size() != 2 ) {
      LOG_BEGIN("FluxBlender", pWARN)
        << "could not parse " << tok1 << " split size=" << pair.size()
        << LOG_END;
      continue;
    }
    int pdg_in   = strtol(pair[0].c_str(),NULL,0);
    int indx_in  = PDG2Indx(pdg_in);
    int pdg_out  = strtol(pair[1].c_str(),NULL,0);
    int indx_out = PDG2Indx(pdg_out);
    for (int jout = 0; jout < 7; ++jout ) {
      fProb[indx_in][jout] = ( ( jout == indx_out ) ? 1 : 0 );
    }

  }
}

int genie::flux::GFlavorMap::PDG2Indx ( int  pdg )

inlineprivate

Definition at line 104 of file GFlavorMap.h.

Referenced by NuName(), ParseFixedfracString(), ParseMapString(), and Probability().

{
  switch ( pdg ) {
  case  12: return 1; break;
  case  14: return 2; break;
  case  16: return 3; break;
  case -12: return 4; break;
  case -14: return 5; break;
  case -16: return 6; break;
  default:  return 0; break;
  }
  return 0;
}

void genie::flux::GFlavorMap::PrintConfig ( bool  verbose = true )

virtual

provide a means of printing the configuration

Implements genie::flux::GFlavorMixerI.

Definition at line 154 of file GFlavorMap.cxx.

References fProb, Indx2PDG(), IndxName(), LOG_BEGIN, LOG_END, and pINFO.

{
  size_t jin, jout;
  LOG_BEGIN("FluxBlender", pINFO)
    << "GFlavorMap::PrintConfig():" << LOG_END;

  //            1234567890[xxx]:
  std::cout << "       in      \\ out  ";
  for (jout=0; jout<7; ++jout )
    std::cout << "  " << std::setw(3) << Indx2PDG(jout) << "    ";
  std::cout << std::endl;

  std::cout << "----------------+";
  for (jout=0; jout<7; ++jout ) std::cout << "----------";
  std::cout << std::endl;

  for (jin=0; jin<7; ++jin) {
    std::cout << std::setw(10) << IndxName(jin)
              << "[" << std::setw(3) << Indx2PDG(jin) << "] |  ";
    for (jout=0; jout<7; ++jout )
      std::cout <<  std::setw(8) << fProb[jin][jout] << " ";
    std::cout << std::endl;
  }

}

double genie::flux::GFlavorMap::Probability ( int  pdg_initial, int  pdg_final, double  energy, double  dist  )

virtual

for any pair of PDG codes the model must calculate the transition probability. This can also depend on neutrino energy (in GeV) and distance (in meters) from the neutrino origin.

Implements genie::flux::GFlavorMixerI.

Definition at line 141 of file GFlavorMap.cxx.

References fProb, LOG_BEGIN, LOG_END, PDG2Indx(), and pINFO.

{
  double prob = fProb[PDG2Indx(pdg_initial)][PDG2Indx(pdg_final)];
  if ( false ) {
    LOG_BEGIN("FluxBlender", pINFO)
      << "Probability " << pdg_initial << "=>" << pdg_final
      << " = " << prob << LOG_END;
  }
  return prob;
}

Member Data Documentation

double genie::flux::GFlavorMap::fProb[7][7]

private

Definition at line 87 of file GFlavorMap.h.

Referenced by GFlavorMap(), ParseFixedfracString(), ParseMapString(), PrintConfig(), and Probability().

---

The documentation for this class was generated from the following files:

• GFlavorMap.h
• GFlavorMap.cxx