BMPTConfig.cxx

Go to the documentation of this file.

#ifndef bmptconfig_cxx
#define bmptconfig_cxx
#include "MessageService/MsgService.h"
#include "MCReweight/BMPTConfig.h"

CVSID("$Id: BMPTConfig.cxx,v 1.2 2005/02/28 14:08:24 cbs Exp $");

//*********************************************
BMPTConfig::BMPTConfig() :
  CPA_PI(0),CPA_K(0),CPA_P(0),CPA_PBAR(0),
  CPB_PI(0),CPB_K(0),CPB_P(0),CPB_PBAR(0),
  ALPHA_PI(0),ALPHA_K(0),ALPHA_P(0),ALPHA_PBAR(0),
  BETA_PI(0),BETA_K(0),BETA_P(0),BETA_PBAR(0),
  GAMMA_PI(0),GAMMA_K(0),GAMMA_P(0),GAMMA_PBAR(0),
  DELTA_PI(0),DELTA_K(0),DELTA_P(0),DELTA_PBAR(0), 
  SMA_PI(0),SMA_K(0),SMA_P(0),SMA_PBAR(0),
  SMB_PI(0),SMB_K(0),SMB_P(0),SMB_PBAR(0),
  R0_PI(0),R0_K(0),R0_P(0),R0_PBAR(0),
  R1_PI(0),R1_K(0),R1_P(0),R1_PBAR(0),
  BEAM_P(0),A_TARGET(0),A_BE(0),
  TARGET_RHO(0),LAMBDA_P(0),LAMBDA_S(0)
{

  MSG("BMPT",Msg::kDebug) << "In BMPTConfig Constructor" << endl;  
  SetDefaults();

}

//*********************************************
BMPTConfig::~BMPTConfig()
{  
}

//*********************************************
void BMPTConfig::SetDefaults()
{

  BEAM_P   = 120.0;
  A_TARGET = 12.0;
  A_BE     = 4.0;

  CPA_PI   = 62.3;
  CPA_K    = 7.74;
  CPA_P    = 8.69;
  CPA_PBAR = 5.20;

  CPB_PI   = 1.57;
  CPB_K    = 0.0;
  CPB_P    = 12.3;
  CPB_PBAR = 0.0;

  ALPHA_PI   = 3.45;
  ALPHA_K    = 2.45;
  ALPHA_P    = 0.0;
  ALPHA_PBAR = 7.56;

  BETA_PI   = 0.517;
  BETA_K    = 0.444;
  BETA_P    = 0.0;
  BETA_PBAR = 0.362;

  SMA_PI   = 6.10;
  SMA_K    = 5.04;
  SMA_P    = 5.77;
  SMA_PBAR = 5.77;

  SMB_PI   = 0.0;
  SMB_K    = 0.0;
  SMB_P    = 1.47;
  SMB_PBAR = 0.0;

  GAMMA_PI   = 0.153;
  GAMMA_K    = 0.121;
  GAMMA_P    = 0.0;
  GAMMA_PBAR = 0.0;

  DELTA_PI   = 0.478;
  DELTA_K    = 2.0*GAMMA_K;
  DELTA_P    = 0.0;
  DELTA_PBAR = 0.0;

  R0_PI   = 1.05;
  R0_K    = 1.15;
  R0_P    = 0.0;
  R0_PBAR = 0.0;

  R1_PI   =  2.65;
  R1_K    = -3.17;
  R1_P    =  0.0;
  R1_PBAR =  0.0;

  TARGET_RHO = 1.754; // g/cm^3
  LAMBDA_P   = 45.0;   // Proton interaction length in cm
  LAMBDA_S   = 45.0;   // Secondary interaction length in cm

}

//*********************************************
void BMPTConfig::Print(ostream & stream)
{
  stream << "=============================" << endl;
  stream << "BMPTConfig parameters:" << endl;
  stream << "----------------------" << endl;
  stream << "CPA_PI   = " << CPA_PI << endl;
  stream << "CPA_K    = " << CPA_K << endl;
  stream << "CPA_P    = " << CPA_P << endl;
  stream << "CPA_PBAR = " << CPA_PBAR << endl;
  
  stream << "CPB_PI   = " << CPB_PI << endl;
  stream << "CPB_K    = " << CPB_K << endl;
  stream << "CPB_P    = " << CPB_P << endl;
  stream << "CPB_PBAR = " << CPB_PBAR << endl;
  
  stream << "ALPHA_PI   = " << ALPHA_PI << endl;
  stream << "ALPHA_K    = " << ALPHA_K << endl;
  stream << "ALPHA_P    = " << ALPHA_P << endl;
  stream << "ALPHA_PBAR = " << ALPHA_PBAR << endl;
  
  stream << "BETA_PI   = " << BETA_PI << endl;
  stream << "BETA_K    = " << BETA_K << endl;
  stream << "BETA_P    = " << BETA_P << endl;
  stream << "BETA_PBAR = " << BETA_PBAR << endl;
  
  stream << "GAMMA_PI   = " << GAMMA_PI << endl;
  stream << "GAMMA_K    = " << GAMMA_K << endl;
  stream << "GAMMA_P    = " << GAMMA_P << endl;
  stream << "GAMMA_PBAR = " << GAMMA_PBAR << endl;
  
  stream << "DELTA_PI   = " << DELTA_PI << endl;
  stream << "DELTA_K    = " << DELTA_K << endl;
  stream << "DELTA_P    = " << DELTA_P << endl;
  stream << "DELTA_PBAR = " << DELTA_PBAR << endl;
  
  stream << "SMA_PI   = " << SMA_PI << endl;
  stream << "SMA_K    = " << SMA_K << endl;
  stream << "SMA_P    = " << SMA_P << endl;
  stream << "SMA_PBAR = " << SMA_PBAR << endl;
  
  stream << "SMB_PI   = " << SMB_PI << endl;
  stream << "SMB_K    = " << SMB_K << endl;
  stream << "SMB_P    = " << SMB_P << endl;
  stream << "SMB_PBAR = " << SMB_PBAR << endl;
  
  stream << "R0_PI   = " << R0_PI << endl;
  stream << "R0_K    = " << R0_K << endl;
  stream << "R0_P    = " << R0_P << endl;
  stream << "R0_PBAR = " << R0_PBAR << endl;
  
  stream << "R1_PI   = " << R1_PI << endl;
  stream << "R1_K    = " << R1_K << endl;
  stream << "R1_P    = " << R1_P << endl;
  stream << "R1_PBAR = " << R1_PBAR << endl;
  
  stream << "BEAM_P   = " << BEAM_P << endl;
  stream << "A_TARGET = " << A_TARGET << endl;
  stream << "A_BE     = " << A_BE << endl;
  
  stream << "TARGET_RHO = " << TARGET_RHO << endl;
  stream << "LAMBDA_P   = " << LAMBDA_P << endl;
  stream << "LAMBDA_S   = " << LAMBDA_S << endl;
  stream << "=============================" << endl;

}

//**********************************8
BMPTConfig BMPTConfig::operator=(BMPTConfig source)
{

  CPA_PI   = source.GetCPA_PI();
  CPA_K    = source.GetCPA_K();
  CPA_P    = source.GetCPA_P();
  CPA_PBAR = source.GetCPA_PBAR();
  
  CPB_PI   = source.GetCPB_PI();
  CPB_K    = source.GetCPB_K();
  CPB_P    = source.GetCPB_P();
  CPB_PBAR = source.GetCPB_PBAR();
  
  ALPHA_PI   = source.GetALPHA_PI();
  ALPHA_K    = source.GetALPHA_K();
  ALPHA_P    = source.GetALPHA_P();
  ALPHA_PBAR = source.GetALPHA_PBAR();
  
  BETA_PI   = source.GetBETA_PI();
  BETA_K    = source.GetBETA_K();
  BETA_P    = source.GetBETA_P();
  BETA_PBAR = source.GetBETA_PBAR();
  
  GAMMA_PI   = source.GetGAMMA_PI();
  GAMMA_K    = source.GetGAMMA_K();
  GAMMA_P    = source.GetGAMMA_P();
  GAMMA_PBAR = source.GetGAMMA_PBAR();
  
  DELTA_PI   = source.GetDELTA_PI();
  DELTA_K    = source.GetDELTA_K();
  DELTA_P    = source.GetDELTA_P();
  DELTA_PBAR = source.GetDELTA_PBAR();
  
  SMA_PI   = source.GetSMA_PI();
  SMA_K    = source.GetSMA_K();
  SMA_P    = source.GetSMA_P();
  SMA_PBAR = source.GetSMA_PBAR();
  
  SMB_PI   = source.GetSMB_PI();
  SMB_K    = source.GetSMB_K();
  SMB_P    = source.GetSMB_P();
  SMB_PBAR = source.GetSMB_PBAR();
  
  R0_PI   = source.GetR0_PI();
  R0_K    = source.GetR0_K();
  R0_P    = source.GetR0_P();
  R0_PBAR = source.GetR0_PBAR();
  
  R1_PI   = source.GetR1_PI();
  R1_K    = source.GetR1_K();
  R1_P    = source.GetR1_P();
  R1_PBAR = source.GetR1_PBAR();
  
  BEAM_P   = source.GetBEAM_P();
  A_TARGET = source.GetA_TARGET();
  A_BE     = source.GetA_BE();
  
  TARGET_RHO = source.GetTARGET_RHO();
  LAMBDA_P   = source.GetLAMBDA_P();
  LAMBDA_S   = source.GetLAMBDA_S();

  return *this;
}


//**********************************8
const double BMPTConfig::GetCPA(int pid) const
{
  if(pid==8||pid==9){
    return CPA_PI;
  }
  else if(pid==10||pid==11||pid==12||pid==16){
    return CPA_K;
  }
  else if(pid==14){
    return CPA_P;
  }
  else if(pid==15){
    return CPA_PBAR;
  }
  return 0;
}

//**********************************8
const double BMPTConfig::GetCPB(int pid) const
{
  if(pid==8||pid==9){
    return CPB_PI;
  }
  else if(pid==10||pid==11||pid==12||pid==16){
    return CPB_K;
  }
  else if(pid==14){
    return CPB_P;
  }
  else if(pid==15){
    return CPB_PBAR;
  }
  return 0;
}

//**********************************8
const double BMPTConfig::GetALPHA(int pid) const
{
  if(pid==8||pid==9){
    return ALPHA_PI;
  }
  else if(pid==10||pid==11||pid==12||pid==16){
    return ALPHA_K;
  }
  else if(pid==14){
    return ALPHA_P;
  }
  else if(pid==15){
    return ALPHA_PBAR;
  }
  return 0;
}

//**********************************8
const double BMPTConfig::GetBETA(int pid) const
{
  if(pid==8||pid==9){
    return BETA_PI;
  }
  else if(pid==10||pid==11||pid==12||pid==16){
    return BETA_K;
  }
  else if(pid==14){
    return BETA_P;
  }
  else if(pid==15){
    return BETA_PBAR;
  }
  return 0;
}

//**********************************8
const double BMPTConfig::GetGAMMA(int pid) const
{
  if(pid==8||pid==9){
    return GAMMA_PI;
  }
  else if(pid==10||pid==11||pid==12||pid==16){
    return GAMMA_K;
  }
  else if(pid==14){
    return GAMMA_P;
  }
  else if(pid==15){
    return GAMMA_PBAR;
  }
  return 0;
}

//**********************************8
const double BMPTConfig::GetDELTA(int pid) const
{
  if(pid==8||pid==9){
    return DELTA_PI;
  }
  else if(pid==10||pid==11||pid==12||pid==16){
    return DELTA_K;
  }
  else if(pid==14){
    return DELTA_P;
  }
  else if(pid==15){
    return DELTA_PBAR;
  }
  return 0;
}

//**********************************8
const double BMPTConfig::GetSMA(int pid) const
{
  if(pid==8||pid==9){
    return SMA_PI;
  }
  else if(pid==10||pid==11||pid==12||pid==16){
    return SMA_K;
  }
  else if(pid==14){
    return SMA_P;
  }
  else if(pid==15){
    return SMA_PBAR;
  }
  return 0;
}

//**********************************8
const double BMPTConfig::GetSMB(int pid) const
{
  if(pid==8||pid==9){
    return SMB_PI;
  }
  else if(pid==10||pid==11||pid==12||pid==16){
    return SMB_K;
  }
  else if(pid==14){
    return SMB_P;
  }
  else if(pid==15){
    return SMB_PBAR;
  }
  return 0;
}

//**********************************8
const double BMPTConfig::GetR0(int pid) const
{
  if(pid==8||pid==9){
    return R0_PI;
  }
  else if(pid==10||pid==11||pid==12||pid==16){
    return R0_K;
  }
  else if(pid==14){
    return R0_P;
  }
  else if(pid==15){
    return R0_PBAR;
  }
  return 0;
}

//**********************************8
const double BMPTConfig::GetR1(int pid) const
{
  if(pid==8||pid==9){
    return R1_PI;
  }
  else if(pid==10||pid==11||pid==12||pid==16){
    return R1_K;
  }
  else if(pid==14){
    return R1_P;
  }
  else if(pid==15){
    return R1_PBAR;
  }
  return 0;
}

//**********************************8
void BMPTConfig::SetCPA(int pid,double val)
{
  if(pid==8||pid==9){
    CPA_PI = val;
  }
  else if(pid==10||pid==11||pid==12||pid==16){
    CPA_K = val;
  }
  else if(pid==14){
    CPA_P = val;
  }
  else if(pid==15){
    CPA_PBAR = val;
  }
}

//**********************************8
void BMPTConfig::SetCPB(int pid,double val)
{
  if(pid==8||pid==9){
    CPB_PI = val;
  }
  else if(pid==10||pid==11||pid==12||pid==16){
    CPB_K = val;
  }
  else if(pid==14){
    CPB_P = val;
  }
  else if(pid==15){
    CPB_PBAR = val;
  }
}

//**********************************8
void BMPTConfig::SetALPHA(int pid,double val)
{
  if(pid==8||pid==9){
    ALPHA_PI = val;
  }
  else if(pid==10||pid==11||pid==12||pid==16){
    ALPHA_K = val;
  }
  else if(pid==14){
    ALPHA_P = val;
  }
  else if(pid==15){
    ALPHA_PBAR = val;
  }
}

//**********************************8
void BMPTConfig::SetBETA(int pid,double val)
{
  if(pid==8||pid==9){
    BETA_PI = val;
  }
  else if(pid==10||pid==11||pid==12||pid==16){
    BETA_K = val;
  }
  else if(pid==14){
    BETA_P = val;
  }
  else if(pid==15){
    BETA_PBAR = val;
  }
}

//**********************************8
void BMPTConfig::SetGAMMA(int pid,double val)
{
  if(pid==8||pid==9){
    GAMMA_PI = val;
  }
  else if(pid==10||pid==11||pid==12||pid==16){
    GAMMA_K = val;
  }
  else if(pid==14){
    GAMMA_P = val;
  }
  else if(pid==15){
    GAMMA_PBAR = val;
  }
}

//**********************************8
void BMPTConfig::SetDELTA(int pid,double val)
{
  if(pid==8||pid==9){
    DELTA_PI = val;
  }
  else if(pid==10||pid==11||pid==12||pid==16){
    DELTA_K = val;
  }
  else if(pid==14){
    DELTA_P = val;
  }
  else if(pid==15){
    DELTA_PBAR = val;
  }
}

//**********************************8
void BMPTConfig::SetSMA(int pid,double val)
{
  if(pid==8||pid==9){
    SMA_PI = val;
  }
  else if(pid==10||pid==11||pid==12||pid==16){
    SMA_K = val;
  }
  else if(pid==14){
    SMA_P = val;
  }
  else if(pid==15){
    SMA_PBAR = val;
  }
}

//**********************************8
void BMPTConfig::SetSMB(int pid,double val)
{
  if(pid==8||pid==9){
    SMB_PI = val;
  }
  else if(pid==10||pid==11||pid==12||pid==16){
    SMB_K = val;
  }
  else if(pid==14){
    SMB_P = val;
  }
  else if(pid==15){
    SMB_PBAR = val;
  }
}

//**********************************8
void BMPTConfig::SetR0(int pid,double val)
{
  if(pid==8||pid==9){
    R0_PI = val;
  }
  else if(pid==10||pid==11||pid==12||pid==16){
    R0_K = val;
  }
  else if(pid==14){
    R0_P = val;
  }
  else if(pid==15){
    R0_PBAR = val;
  }
}

//**********************************8
void BMPTConfig::SetR1(int pid,double val)
{
  if(pid==8||pid==9){
    R1_PI = val;
  }
  else if(pid==10||pid==11||pid==12||pid==16){
    R1_K = val;
  }
  else if(pid==14){
    R1_P = val;
  }
  else if(pid==15){
    R1_PBAR = val;
  }
}

#endif

---