#include <EvtXsection.hh>

Public Member Functions
	EvtXsection (std::vector< EvtId > evtdaugs)

	EvtXsection (int mode)

virtual	~EvtXsection ()

int	getMode (std::vector< EvtId > evtdaugs)

void	ini_data0 (int mode)

void	ini_data (int mode)

void	ini_data_diy ()

void	ini_data_multimode ()

double	getXsection (double mx)

double	getErr (double mx)

double	Xsection_a (double mx)

double	Xsection_b (double mx)

double	Xsection_c (double mx)

double	Err_a (double mx)

double	Err_b (double mx)

int	getXBin (double mx, std::vector< double > vy)

int	getXBin_a (double mx)

int	getXBin_b (double mx)

std::string	getUnit ()

std::vector< double >	getXX ()

std::vector< double >	getYY ()

std::vector< double >	getEr ()

double	getXup ()

double	getXlw ()

std::string	getMsg ()

void	setBW (int pdg)

void	setFile (std::string name)

double	getVP (double mx)

void	ReadVP ()

void	setModes (std::vector< int > vmd)

Detailed Description

Definition at line 132 of file EvtXsection.hh.

Constructor & Destructor Documentation

◆ EvtXsection() [1/2]

EvtXsection::EvtXsection ( std::vector< EvtId > evtdaugs )

inline

Definition at line 134 of file EvtXsection.hh.

                                        {
    xx.clear();yy.clear();er.clear();
    _mode = getMode(evtdaugs);
    ini_data(_mode);
  
  }

◆ EvtXsection() [2/2]

EvtXsection::EvtXsection ( int mode )

inline

Definition at line 140 of file EvtXsection.hh.

                       {
    xx.clear();yy.clear();er.clear();
    _mode = mode;
    if(mode>=0){
    ini_data(_mode);
    }else if(mode==-1){//calculate the correction factor with Breit-wigner
      _unit="nb";
    }else if(mode==-2){//users provide the cross section list
      //std::cout<<"user diy"<<std::endl;
       ini_data_diy();
    }else if(mode==-100){//multi-exclusive modes
      _unit="nb";
    }
    //std::cout<<"The mode is : "<<_mode<<", "<<xx[0]<<", "<<yy[0]<<std::endl;
  }

◆ ~EvtXsection()

EvtXsection::~EvtXsection ( )

virtual

Definition at line 124 of file EvtXsection.cc.

124{}

Member Function Documentation

◆ Err_a()

double EvtXsection::Err_a ( double mx )

Definition at line 291 of file EvtXsection.cc.

                                  {
  if(_mode == -1) return 0.;
  if (mx <= xx[0]) return er[0];
  if (mx > xx[nbins]) return 0.;
  int thebin = getXBin_a(mx);
  //-- debug
  //std::cout<<"Mode= "<<_mode<<", thebin"<<thebin<<std::endl;
  return er[thebin]; 
}

◆ Err_b()

double EvtXsection::Err_b ( double mx )

Definition at line 301 of file EvtXsection.cc.

                                  {// interpolation of the cross section error between two bins
  if (mx <= xx[0]) return er[0];
  if (mx > xx[nbins-1]) return 0.;
  int thebin = getXBin_b(mx);
  double ylow = er[thebin];
  double yup  = er[thebin+1];
 
  double xlow = xx[thebin];
  double xup  = xx[thebin+1]; 
 
  double yi=ylow + (yup-ylow)/(xup-xlow)*(mx-xlow);
  return yi;
}

Referenced by getErr().

◆ getEr()

std::vector< double > EvtXsection::getEr ( )

inline

Definition at line 179 of file EvtXsection.hh.

179{return er;}

◆ getErr()

double EvtXsection::getErr ( double mx )

Definition at line 260 of file EvtXsection.cc.

                                   {
  //if(_mode <= 5 && _mode !=1 ){ return Err_a(mx);}
  //else {return  Err_b(mx);}
  return  Err_b(mx);
}

Referenced by ini_data_multimode().

◆ getMode()

int EvtXsection::getMode ( std::vector< EvtId > evtdaugs )

Definition at line 218 of file EvtXsection.cc.

                                                 {
  /*******************--- mode definition:
  0: ppbar
  1: nnbar
  2: Lambda0 anti-Lambda0
  3: Sigma0 anti-Sigma0
  4: Lambda0 anti-Sigma0
  5: Sigma0 anti-Lambda0
  6: pi+ pi-
  7: pi+ pi- pi0
  8: K+K- pi0
  9: KsK+pi- 
  10: KsK-pi+
  11: K+K-eta
  *************************************/
  int pp[2]={-2212,2212};
  int nn[2]={-2112,2112};
  int pipi[2]={-211,211};
  int pi3[3]={-211,111,211};
  int kkpi0[3]={-321,111,321}; //K+K-pi0
  int kskpi[3]={-211,310,321};//KsK+pi-
 
  std::vector<int> vson;vson.clear();
  for(int i=0;i<evtdaugs.size();i++){vson.push_back( EvtPDL::getStdHep(evtdaugs[i]) );}
  sort(vson.begin(),vson.end(),mysort);
 
  if(vson.size() ==2 ){
    if(vson[0]==pp[0] && vson[1]==pp[1] ) {return 0;}
    if(vson[0]==nn[0] && vson[1]==nn[1] ) {return 1;}
    if(vson[0]==pipi[0] && vson[1]==pipi[1] ) {return 2;}
  }else if(vson.size()==3){
 
  }
 
}

Referenced by EvtXsection().

◆ getMsg()

std::string EvtXsection::getMsg ( )

inline

Definition at line 182 of file EvtXsection.hh.

182{return msg;}

◆ getUnit()

std::string EvtXsection::getUnit ( )

inline

Definition at line 175 of file EvtXsection.hh.

175{return _unit;}

Referenced by ini_data_multimode().

◆ getVP()

double EvtXsection::getVP ( double mx )

◆ getXBin()

int EvtXsection::getXBin	(	double	mx,
		std::vector< double >	vy )

Definition at line 331 of file EvtXsection.cc.

                                                      {
  double nbins = vy.size();
  if(vy[0]<mx || mx<vy[nbins-1]) {std::cout<<"Outside vacuum pol. value"<<std::endl;abort();}
  for(int i=0;i<nbins-1;i++){
    if(mx <= vy[i+1] && mx > vy[i]) return i; //mx at i to i+1 bin
  }
}

◆ getXBin_a()

int EvtXsection::getXBin_a ( double mx )

Definition at line 315 of file EvtXsection.cc.

                                   {
  if(mx <= xx[0]) return 0;
  if(mx > xx[nbins]) return nbins;
  for(int i=0;i<nbins;i++){
    if(mx <= xx[i+1] && mx > xx[i]) return i; //mx at i-th bin
  }
}

Referenced by Err_a(), and Xsection_a().

◆ getXBin_b()

int EvtXsection::getXBin_b ( double mx )

Definition at line 323 of file EvtXsection.cc.

                                   {
  if(mx <= xx[0]) return 0;
  if(mx > xx[nbins-1]) return nbins;
  for(int i=0;i<nbins-1;i++){
    if(mx <= xx[i+1] && mx > xx[i]) return i; //mx at i to i+1 bin
  }
}

Referenced by Err_b(), and Xsection_b().

◆ getXlw()

double EvtXsection::getXlw ( )

inline

Definition at line 181 of file EvtXsection.hh.

181{return xx[0];}

Referenced by ini_data_multimode().

◆ getXsection()

double EvtXsection::getXsection ( double mx )

Definition at line 254 of file EvtXsection.cc.

                                        {
  if(_mode == -1){return Xsection_c(mx);} 
  //else if(_mode == 0 || _mode == 2 || _mode == 3 || _mode == 4 ||_mode == 5){ return Xsection_a(mx);}
  else {return  Xsection_b(mx);}
}

Referenced by ini_data_multimode().

◆ getXup()

double EvtXsection::getXup ( )

inline

Definition at line 180 of file EvtXsection.hh.

180{return xx[nbins-1];}

Referenced by ini_data_multimode().

◆ getXX()

std::vector< double > EvtXsection::getXX ( )

inline

Definition at line 177 of file EvtXsection.hh.

177{return xx;}

◆ getYY()

std::vector< double > EvtXsection::getYY ( )

inline

Definition at line 178 of file EvtXsection.hh.

178{return yy;}

◆ ini_data()

void EvtXsection::ini_data ( int mode )

Referenced by EvtXsection(), EvtXsection(), and ini_data_multimode().

◆ ini_data0()

void EvtXsection::ini_data0 ( int mode )

Definition at line 126 of file EvtXsection.cc.

                                   {
 
  xx.clear();yy.clear();er.clear();
  nbins=yy.size();
 
} 

◆ ini_data_diy()

void EvtXsection::ini_data_diy ( )

Definition at line 133 of file EvtXsection.cc.

                              {//user provide xs list
 
  xx.clear();yy.clear();er.clear();
 
  ifstream file("xs_user.txt");
 
  if (!file){
    cout << "EvtXsection.cc: The input file not found. The default file name should be xs_user.txt!"<<std::endl;
    exit(0);
  }
 
  double xm,xs,xs_er;
 
  while(!file.eof()){
    file>>xm>>xs>>xs_er;
    //std::cout<<"read XS: "<<xm<<" "<<xs<<" "<<xs_er<<std::endl;
    xx.push_back(xm);
    yy.push_back(xs);
    er.push_back(xs_er);
  }
 
  xx.pop_back();
  yy.pop_back();
  er.pop_back();
  nbins=yy.size();
  file.close();
  _unit="";
  msg="";
 
}

Referenced by EvtXsection().

◆ ini_data_multimode()

void EvtXsection::ini_data_multimode ( )

Definition at line 164 of file EvtXsection.cc.

                                    {//multi-exclusive modes
  xx.clear();yy.clear();er.clear();
 
  if(_vmd.size()==0){std::cout<<"EvtXsection::ini_data_multimode: No mode indexes are available"<<std::endl; abort();}
 
  double xL=10.0;
  double xU=-1.0;
 
  for(int i=0;i<_vmd.size();i++){
    ini_data(_vmd[i]);
    double x0=getXlw();
    double x1=getXup();
    if(x0<xL) xL=x0; 
    if(x1>xU) xU=x1;
  }
 
  std::cout<<"The low and up end of multimodes: "<<xL<<" ~ "<<xU<<std::endl; 
  double stp=0.0005; //5 MeV for bin width
  double xm;
  double xs=0; //xsection in nb;
  double xs_er=0;
  double mypb=1;
  std::vector<double>uxx,uyy,uer;
  uxx.clear();uyy.clear();uer.clear();
  for(double xxm=xL;xxm<xU;xxm+=stp){
    xm=xxm;
    xs   =0;
    xs_er=0;
    for(int  i=0;i<_vmd.size();i++){
      ini_data(_vmd[i]);
      std::string myunit=getUnit();
      if(myunit=="pb"){ mypb=0.001;}else{mypb=1;}      
      xs   += mypb*getXsection(xxm);
      xs_er+= mypb*getErr(xxm);
      //std::cout<<_vmd[i]<< ": "<<xm<<" "<<xs<<" "<<xs_er<<std::endl;
    }//loop over mode
    uxx.push_back(xm);
    uyy.push_back(xs);
    uer.push_back(xs_er);
  } //loop over mass
 
  xx.clear();yy.clear();er.clear();
  for(int i=0;i<uxx.size();i++){
    xx.push_back(uxx[i]);yy.push_back(uyy[i]);er.push_back(uer[i]);
  }
  //debugging
  //for(int i=0;i<yy.size();i++){
  //  std::cout<<xx[i]<<" "<<yy[i]<<std::endl;
  //}
  _unit="nb";
  nbins=yy.size();
  msg="multi-exclusive mode";
}

◆ ReadVP()

void EvtXsection::ReadVP ( )

◆ setBW()

void EvtXsection::setBW ( int pdg )

Definition at line 358 of file EvtXsection.cc.

                              {
  pdgcode = pdg;
  EvtId pid = EvtPDL::evtIdFromStdHep(pdgcode );
  double maxM= EvtPDL::getMaxMass(pid);
  double minM= EvtPDL::getMinMass(pid);
  double meanM = EvtPDL::getMeanMass(pid);
  double width = EvtPDL::getWidth(pid);
  // std::cout<<minM<<" <=M<= "<<maxM<<std::endl;
  double xstp=(maxM-minM)/100.;
  xx.clear();yy.clear();er.clear();
  for(int i=0;i<100;i++){
    double m=i*xstp;
    EvtComplex im(0.,1.);
    EvtComplex bw=1./(m*m-meanM*meanM+im*m*width);
    double amps = abs2(bw);
    xx.push_back(m); yy.push_back(amps); er.push_back(0.);
  }
  nbins=yy.size();
}

◆ setFile()

void EvtXsection::setFile ( std::string name )

inline

Definition at line 184 of file EvtXsection.hh.

184{file = name;}

◆ setModes()

void EvtXsection::setModes ( std::vector< int > vmd )

Definition at line 378 of file EvtXsection.cc.

                                            {
  _vmd.clear();
  for(int i=0;i<vmd.size();i++){
    _vmd.push_back(vmd[i]);
  }
}

◆ Xsection_a()

double EvtXsection::Xsection_a ( double mx )

Definition at line 267 of file EvtXsection.cc.

                                       {
  if (mx <= xx[0]) return 0.;
  if (mx > xx[nbins]) return 0.;
  int thebin = getXBin_a(mx);
  //-- debug
  //std::cout<<"Mode= "<<_mode<<", thebin"<<thebin<<std::endl;
  return yy[thebin]; 
}

◆ Xsection_b()

double EvtXsection::Xsection_b ( double mx )

Definition at line 276 of file EvtXsection.cc.

                                       {// interpolation of the cross section between two bins
  if (mx <= xx[0]) return 0.;
  if (mx > xx[nbins-1]) return 0.;
  int thebin = getXBin_b(mx);
  double ylow = yy[thebin];
  double yup  = yy[thebin+1];
 
  double xlow = xx[thebin];
  double xup  = xx[thebin+1]; 
 
  double yi=ylow + (yup-ylow)/(xup-xlow)*(mx-xlow);
  return yi;
}

Referenced by getXsection().

◆ Xsection_c()

double EvtXsection::Xsection_c ( double mx )

Definition at line 339 of file EvtXsection.cc.

                                       {// in unit nb
   double pi=3.1415926;
   double s= mx*mx;
   EvtId pid = EvtPDL::evtIdFromStdHep(pdgcode );
   double mass = EvtPDL::getMeanMass(pid);
   double width = EvtPDL::getWidth(pid);
   double mass2 = mass*mass;
   double width2 = width*width;
   double br = (s - mass2)*(s - mass2) + mass2 * width2;
   bree = 1; //this value is canceled when calculation the correction factor;
   double sigma = 12*pi*bree*width2/br;
   double nbar = 4E05;  // ! GeV-2 = 4*10^5 nbar
   double thexs = sigma*nbar;
   // std::cout<<"EvtXsection::Xsection_c: "<<mass<<" "<<width<<" "<<thexs<<std::endl;
   // msg = "Calculate the correction factor for " + EvtPDL::name(pid);
   
   return thexs;
 }

Referenced by getXsection().

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

Public Member Functions

Detailed Description

Constructor & Destructor Documentation

◆ EvtXsection() [1/2]

◆ EvtXsection() [2/2]

◆ ~EvtXsection()

Member Function Documentation

◆ Err_a()

◆ Err_b()

◆ getEr()

◆ getErr()

◆ getMode()

◆ getMsg()

◆ getUnit()

◆ getVP()

◆ getXBin()

◆ getXBin_a()

◆ getXBin_b()

◆ getXlw()

◆ getXsection()

◆ getXup()

◆ getXX()

◆ getYY()

◆ ini_data()

◆ ini_data0()

◆ ini_data_diy()

◆ ini_data_multimode()

◆ ReadVP()

◆ setBW()

◆ setFile()

◆ setModes()

◆ Xsection_a()

◆ Xsection_b()

◆ Xsection_c()