EvtXsection.hh

Go to the documentation of this file.

//--------------------------------------------------------------------------
// Environment:
//      This software is part of models developed at BES collaboration
//      based on the EvtGen framework.  If you use all or part
//      of it, please give an appropriate acknowledgement.
//
// Copyright Information: See EvtGen/BesCopyright
//      Copyright (A) 2006  Ping Rong-Gang @IHEP, Wang Weiping @USTC
//
// Module: EvtXsection.hh
//
// Description: To define cross section for the continuum exclusive process
// Experimental cross section taken from PRD73,012005, PRD76,092006, also
// see a review: Rev. Mod. Phys. 83,1545
// Modification history:
//
//    Ping R.-G.       Nov., 2012       Module created
//
/*******************--- mode definition: also see EvtXsection.cc
mode==0: ppbar
mode==1: nnbar
mode==2: Lambda0 anti-Lambda0
mode==3: Sigma0 anti-Sigma0
mode==4: Lambda0 anti-Sigma0
mode==5: Sigma0 anti-Lambda0
mode==6: pi+ pi-
mode==7: pi+ pi- pi0
mode==8: K+K- pi0
mode==9: KsK+pi-
mode==10: KsK-pi+
mode==11: K+K-eta
mode==12: 2(pi+pi-)
mode==13: pi+pi-2pi0
mode==14: K+K-pi+pi-
mode==15: K+K-2pi0
mode==16: 2(K+K-)
mode==17: 2(pi+pi-)pi0
mode==18: 2(pi+pi-)eta
mode==19: K+K-pi+pi-pi0
mode==20: K+K-pi+pi-eta
mode==21: 3(pi+pi-)
mode==22: 2(pi+pi-pi0)
mode==23: phi eta
mode==24: phi pi0
mode==25: K+K*-
mode==26: K-K*+
mode==27: K_SK*0-bar
mode==28: K*0(892)K+pi-
mode==29: K*0(892)K-pi+
mode==30: K*+K-pi0
mode==31: K*-K+pi0
mode==32: K_2*(1430)0 K+pi-
mode==33: K_2*(1430)0 K-pi+
mode==34: K+K-rho
mode==35: phi pi+pi-
mode==36: phi f0(980)
mode==37: eta pi+pi-
mode==38: omega pi+ pi-
mode==39: omega f0(980)
mode==40: eta' pi+ pi-
mode==41: f_1(1285)pi+pi-
mode==42: omega K+K-
mode==43: omega pi+pi-pi0
mode==44: Sigma+ Sigma- (Sigma0 Sigma0-bar SU(3) extention )
mode==45: K+K-
mode==46: K_S K_L
mode==47: omega eta
mode==48: phi eta'
mode==49: phi K+K-
mode==50: ppbar pi0
mode==51: ppbar eta
mode==52: omega pi0
 
mode==70: D0 anti-D0
mode==71: D+D-
mode==72: D+D*-
mode==73: D-D*+
mode==74: D*+D*-
mode==68: D0 anti-D*0
mode==69: anti-D0 D*0
mode==67: D*0 anti-D*0
 
mode==59: D0 anti-D0 pi0
mode==63: D+D-pi0
mode==75: D0D-pi+
mode==76: anti-D0 D+pi-
 
mode==77: D0D*-pi+
mode==78: anti-D0 D*+pi-
mode==65: D-D*0pi+
mode==66: D+ anti-D*0pi-
mode==60: anti-D0 D*0pi0
mode==61: D0 anti-D*0pi0
mode==62: D+D*-pi0
mode==64: D-D*+pi0
 
mode==90: pi+pi- J/psi
mode==92: pi0pi0 J/psi
mode==91: pi+pi- psi(2S)
mode==79: pi0pi0 psi(2S)
mode==80: eta J/psi
mode==81: pi+pi- h_c
mode==82: pi0pi0 h_c
mode==83: K+K- J/psi
mode==84: KsKs J/psi
 
mode==93: D_s+ D_s-
mode==94: D_s^{*+}D_s^-
mode==95: D_s^{*-}D_s^+
mode==85: D_s^{*+}D_s^{*-}
 
mode==96: Lambda_c+ anti-Lamba_c-
//---
*************************************/
//
#ifndef EVTXSECTION_HH
#define EVTXSECTION_HH
 
#include "../EvtGenBase/EvtComplex.hh"
#include "../EvtGenBase/EvtId.hh"
#include <algorithm>
#include <fstream>
#include <iostream>
#include <math.h>
#include <ostream>
#include <stdlib.h>
#include <string.h>
#include <string>
#include <vector>
using namespace std;
 
class EvtXsection {
public:
  EvtXsection( std::vector<EvtId> evtdaugs ) {
    xx.clear();
    yy.clear();
    er.clear();
    _mode = getMode( evtdaugs );
    ini_data( _mode );
  }
  EvtXsection( int mode ) {
    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;
  }
 
  virtual ~EvtXsection();
 
  int    getMode( std::vector<EvtId> evtdaugs );
  void   ini_data0( int mode ); // exper. cross section data table
  void   ini_data( int mode );  // fitted results on the cross section of expe. xsection
  void   ini_data_diy();
  void   ini_data_multimode();
  double getXsection( double mx );
  double getErr( double mx );
  double Xsection_a( double mx ); // for bin given by up and lower edge
  double Xsection_b( double mx ); // for bin given by a center value
  double Xsection_c( double mx ); // for xsection calculated with breit-wigner
 
  double Err_a( double mx ); // for bin given by up and lower edge
  double Err_b( double mx ); // for bin given by a center value
 
  int         getXBin( double mx, std::vector<double> vy );
  int         getXBin_a( double mx ); // for bin given by up and lower edge
  int         getXBin_b( double mx ); // for bin given by a center value
  std::string getUnit() { return _unit; }
 
  std::vector<double> getXX() { return xx; }
  std::vector<double> getYY() { return yy; }
  std::vector<double> getEr() { return er; }
  double              getXup() { return xx[nbins - 1]; }
  double              getXlw() { return xx[0]; }
  std::string         getMsg() { return msg; }
  void                setBW( int pdg );
  void                setFile( std::string name ) { file = name; }
  double              getVP( double mx );
  void                ReadVP();
  void                setModes( std::vector<int> vmd );
 
private:
  // e+e- --> ppbar mode, xsection taken from PRD73,012005
  std::vector<double> xx, yy, er;
  int                 nbins, _mode;
  std::string         _unit, msg;
  // calculate cross section with Breit-Wigner
  int    pdgcode;
  double bree;
  // user provide xsection list
  std::string         file;
  std::vector<double> vpx, vpr, vpi;
  std::vector<int>    _vmd;
};
 
#endif