b2bii modules

Collaboration diagram for b2bii modules:

Classes
class	B2BIIMCParticlesMonitorModule
	Declaration of class B2BIIMCParticlesMonitor. More...
class	B2BIIConvertBeamParamsModule
	Convert the IpProfile and Benergy to BeamParameters, BeamSpot, CollisionInvariantMass and CollisionBoostVector. More...
class	B2BIIConvertMdstModule
	Module converts Belle MDST objects (Panther records) to Belle II MDST objects. More...
class	B2BIIFixMdstModule
	Declaration of class B2BIIFixMdst. More...
class	MuidProb
	Class computes probability density for Muid calculation. More...
class	B2BIIMdstInputModule
	Module reads Belle MDST files and converts Belle_event record to Belle II EventMetaData StoreObjectPointer. More...
struct	cal_scale_error_func_set_t
	Structure type cal_scale_error_func_set_t. More...
class	BelleMCOutputModule
	KLM digitization module. More...

Macros
#define	NOMINAL_ENERGY   5.290024915
	Nomial beam energy.

Typedefs
typedef void(*	cal_scale_error_func_t) (double scale[5], const double pt, const double tanl)
	Function pointer type.

Functions
	REG_MODULE (B2BIIConvertBeamParams)
	Register the module.
static double	ecl_adhoc_corr (int Exp, int Run, int iflag05th, double Energy, double)
	The function giving correction factor. More...
static double	ecl_adhoc_corr_45 (int exp, int, int cid)
	The function giving correction factor in Exp.45.
static double	ecl_mcx3_corr (int, int, double energy, double)
	Correct energy scale (MC) to make pi0 peak nominal.
static double	mpi0pdg (double Energy)
	Make MC mass peak to PDG value.
static bool	operator== (const cal_scale_error_func_set_t &lhs, const cal_scale_error_func_set_t &rhs)
	Operator definition for cal_scale_error_func_set_t.
static void	null_scale (double[5], double, double)
	Dummy function.
static void	get_event_id (int *no_exp, int *no_run, int *no_evt, int *no_frm, int *expmc)
	Get event ID.
static bool	is_already_scaled (void)
	Check if event is already scaled.
double	vfunc (const double x, const double x1, const double yc, const double a1, const double a2)
	vfunc
double	cupfunc (const double x, const double x1, const double x2, const double yc, const double a1, const double a2)
	cupfunc
double	rootfunc (const double x, const double x1, const double x2, const double yc, const double a1, const double a2)
	rootfunc
double	lambdafunc (const double x, const double x1, const double x2, const double yc, const double a1, const double a2, const double a3)
	lambdafunc
static void	cal_scale_error_EXP0723_cosmic_mc (double scale[5], const double pt, const double)
	Scale error for Exp.7-23 Cosmic MC.
static void	cal_scale_error_EXP2527_cosmic_mc (double scale[5], const double pt, const double)
	Scale error for Exp.25-27 Cosmic MC.
static void	cal_scale_error_EXP31_cosmic_mc (double scale[5], double pt, double)
	Scale error for Exp.31 Cosmic MC.
static void	cal_scale_error_EXP33_cosmic_mc (double scale[5], double pt, double)
	Scale error for Exp.33 Cosmic MC.
static void	cal_scale_error_EXP35_cosmic_mc (double scale[5], double pt, double tanl)
	Scale error for Exp.35 Cosmic MC.
static void	cal_scale_error_EXP37_cosmic_mc (double scale[5], double pt, double)
	Scale error for Exp.37 Cosmic MC.
static void	cal_scale_error_EXP3941_cosmic_mc (double scale[5], double pt, double)
	Scale error for Exp.39,41 Cosmic MC.
static void	cal_scale_error_EXP43_cosmic_mc (double scale[5], double pt, double)
	Scale error for Exp.43 Cosmic MC.
static void	cal_scale_error_EXP45_cosmic_mc (double scale[5], double pt, double)
	Scale error for Exp.45 Cosmic MC.
static void	cal_scale_error_EXP47_cosmic_mc (double scale[5], double pt, double)
	Scale error for Exp.47 Cosmic MC.
static void	cal_scale_error_EXP49_cosmic_mc (double scale[5], double pt, double)
	Scale error for Exp.49 Cosmic MC.
static void	cal_scale_error_EXP51_cosmic_mc (double scale[5], double pt, double)
	Scale error for Exp.51 Cosmic MC.
static void	cal_scale_error_EXP53_cosmic_mc (double scale[5], double pt, double)
	Scale error for Exp.53 Cosmic MC.
static void	cal_scale_error_EXP55_cosmic_mc (double scale[5], double pt, double)
	Scale error for Exp.55 Cosmic MC.
static void	cal_scale_error_EXP6165_cosmic_mc (double scale[5], double pt, double tanl)
	Scale error for Exp.61-65 Cosmic MC.
static void	cal_scale_error_EXP67_cosmic_mc (double scale[5], double pt, double tanl)
	Scale error for Exp.67 Cosmic MC.
static void	cal_scale_error_EXP6971_cosmic_mc (double scale[5], double pt, double tanl)
	Scale error for Exp.69-71 Cosmic MC.
static void	cal_scale_error_EXP0723_cosmic_data (double scale[5], const double pt, const double)
	Scale error for Exp.7-23 Cosmic data.
static void	cal_scale_error_EXP2527_cosmic_data (double scale[5], const double pt, const double)
	Scale error for Exp.25-27 Cosmic data.
static void	cal_scale_error_EXP31_cosmic_data (double scale[5], double pt, double)
	Scale error for Exp.31 Cosmic data.
static void	cal_scale_error_EXP33_cosmic_data (double scale[5], double pt, double)
	Scale error for Exp.33 Cosmic data.
static void	cal_scale_error_EXP35_cosmic_data (double scale[5], double pt, double tanl)
	Scale error for Exp.35 Cosmic data.
static void	cal_scale_error_EXP37_cosmic_data (double scale[5], double pt, double)
	Scale error for Exp.37 Cosmic data.
static void	cal_scale_error_EXP3941_cosmic_data (double scale[5], double pt, double)
	Scale error for Exp.39-41 Cosmic data.
static void	cal_scale_error_EXP43_cosmic_data (double scale[5], double pt, double)
	Scale error for Exp.43 Cosmic data.
static void	cal_scale_error_EXP4547_cosmic_data (double scale[5], double pt, double)
	Scale error for Exp.45,47 Cosmic data.
static void	cal_scale_error_EXP49_cosmic_data (double scale[5], double pt, double)
	Scale error for Exp.49 Cosmic data.
static void	cal_scale_error_EXP51_cosmic_data (double scale[5], double pt, double tanl)
	Scale error for Exp.51 Cosmic data.
static void	cal_scale_error_EXP53_cosmic_data (double scale[5], double pt, double)
	Scale error for Exp.53 Cosmic data.
static void	cal_scale_error_EXP55_cosmic_data (double scale[5], double pt, double tanl)
	Scale error for Exp.55 Cosmic data.
static void	cal_scale_error_EXP6165_cosmic_data (double scale[5], double pt, double tanl)
	Scale error for Exp.61-65 Cosmic data.
static void	cal_scale_error_EXP67_cosmic_data (double scale[5], double pt, double tanl)
	Scale error for Exp.67 Cosmic data.
static void	cal_scale_error_EXP6971_cosmic_data (double scale[5], double pt, double tanl)
	Scale error for Exp.69,71 Cosmic data.
static void	cal_scale_error_EXP0723_hadronic_mc (double scale[5], const double pt, const double tanl)
	Scale error for Exp.7-23 Hadron MC.
static void	cal_scale_error_EXP2527_hadronic_mc (double scale[5], const double pt, const double tanl)
	Scale error for Exp.25-27 Hadron MC.
static void	cal_scale_error_EXP31_hadronic_mc (double scale[5], double pt, double tanl)
	Scale error for Exp.31 Hadron MC.
static void	cal_scale_error_EXP33_hadronic_mc (double scale[5], double pt, double tanl)
	Scale error for Exp.33,35 Hadron MC.
static void	cal_scale_error_EXP35_hadronic_mc (double scale[5], double pt, double tanl)
	Scale error for Exp.33,35 Hadron MC.
static void	cal_scale_error_EXP37_hadronic_mc (double scale[5], double pt, double tanl)
	Scale error for Exp.37 Hadron MC.
static void	cal_scale_error_EXP3941_hadronic_mc (double scale[5], double pt, double)
	Scale error for Exp.39,41 Hadron MC.
static void	cal_scale_error_EXP43_hadronic_mc (double scale[5], double pt, double)
	Scale error for Exp.43 Hadron MC.
static void	cal_scale_error_EXP45_hadronic_mc (double scale[5], double pt, double tanl)
	Scale error for Exp.45 Hadron MC.
static void	cal_scale_error_EXP47_hadronic_mc (double scale[5], double pt, double)
	Scale error for Exp.47 Hadron MC.
static void	cal_scale_error_EXP49_hadronic_mc (double scale[5], double pt, double)
	Scale error for Exp.49 Hadron MC.
static void	cal_scale_error_EXP51_hadronic_mc (double scale[5], double pt, double tanl)
	Scale error for Exp.51 Hadron MC.
static void	cal_scale_error_EXP53_hadronic_mc (double scale[5], double pt, double)
	Scale error for Exp.53 Hadron MC.
static void	cal_scale_error_EXP55_hadronic_mc (double scale[5], double pt, double)
	Scale error for Exp.55 Hadron MC.
static void	cal_scale_error_EXP6165_hadronic_mc (double scale[5], double pt, double)
	Scale error for Exp.61-65 Hadron MC.
static void	cal_scale_error_EXP67_hadronic_mc (double scale[5], double pt, double)
	Scale error for Exp.67 Hadron MC.
static void	cal_scale_error_EXP6971_hadronic_mc (double scale[5], double pt, double)
	Scale error for Exp.69,71 Hadron MC.
static cal_scale_error_func_set_t	get_scale_error_func_for_exprun (const int no_exp, const int)
	Get scale error fucntion for different Exp.
static void	cal_scale_error (double scale[5], const double pt, const double tanl, const int expmc, const int no_exp, const int no_run)
	Calculate scale error.
static int	scale_error_impl (const int message_level, const int reprocess_version)
	The implementation of B2BIIFixMdstModule::scale_error() More...
void	recsim_mdst_propgt_ (float *, float[], float[], float[], float[], float[], int *)
	recsim_mdst_propgt from legacy C code
static void	scale_err_ms (Belle::Mdst_trk_fit &fit, const double scale[])
	Scale error.
static void	smear_trk_ms (Belle::Mdst_trk_fit &fit, const double scale[])
	Smear MC tracks.
static void	smear_charged ()
	Smear tracks in Mdst_Charged.
	B2BIIConvertBeamParamsModule ()
	Create parameters.
void	initialize () override
	Initialize phanter banks.
void	beginRun () override
	Set run info in panther and load IPProfile/Benergy and convert the values to payloads.
	B2BIIFixMdstModule ()
	Constructor.
virtual void	initialize () override
	Initialize the module.
virtual void	terminate () override
	Terminates the module.
virtual void	event () override
	Called for each event.
virtual void	beginRun () override
	Called when a new run is started.

Variables
static int	SE_Message_Level
	Message level of scale_error(). More...
static int	SE_Reprocess_Version = 0
	Reprocess verison used in scale_error()
static const struct cal_scale_error_func_set_t	EXP0723_scale
	Scale error for Exp.7-23. More...
static const struct cal_scale_error_func_set_t	EXP2527_scale
	Scale error for Exp.25-27. More...
static const struct cal_scale_error_func_set_t	EXP31_scale
	Scale error for Exp.31. More...
static const struct cal_scale_error_func_set_t	EXP33_scale
	Scale error for Exp.33. More...
static const struct cal_scale_error_func_set_t	EXP35_scale
	Scale error for Exp.35. More...
static const struct cal_scale_error_func_set_t	EXP37_scale
	Scale error for Exp.37. More...
static const struct cal_scale_error_func_set_t	EXP3941_scale
	Scale error for Exp.39,41. More...
static const struct cal_scale_error_func_set_t	EXP43_scale
	Scale error for Exp.43. More...
static const struct cal_scale_error_func_set_t	EXP45_scale
	Scale error for Exp.45. More...
static const struct cal_scale_error_func_set_t	EXP47_scale
	Scale error for Exp.47. More...
static const struct cal_scale_error_func_set_t	EXP49_scale
	Scale error for Exp.49. More...
static const struct cal_scale_error_func_set_t	EXP51_scale
	Scale error for Exp.51. More...
static const struct cal_scale_error_func_set_t	EXP53_scale
	Scale error for Exp.53. More...
static const struct cal_scale_error_func_set_t	EXP55_scale
	Scale error for Exp.55. More...
static const struct cal_scale_error_func_set_t	EXP6165_scale
	Scale error for Exp.61-65. More...
static const struct cal_scale_error_func_set_t	EXP67_scale
	Scale error for Exp.67. More...
static const struct cal_scale_error_func_set_t	EXP6971_scale
	Scale error for Exp.69-71. More...
static const struct cal_scale_error_func_set_t	DUMMY_scale
	Dummy scale. More...
static double	s_benergy_value = NOMINAL_ENERGY
	Beam energy.

Muid user parameters (set via basf)
int	get_reprocess_version ()
	Get reprocess version of input mdst.
bool	good_event ()
	Check if event correspondes to the period when Belle detector not in good condition.
double	Benergy (int expnum=0, int runnum=0)
	Return Beam energy.
void	fix_pi0_brecon (void)
	Fix relation gamma<->pi0 in brecon table.
static double	pi0resol (double, double, const char *, bool, int, int)
	Treat pi0 mass width as a func. More...
void	correct_ecl (int, int)
	Correct photon's momenta and error matrix.
void	make_pi0 (int, double, double)
	Create Mdst_pi0 from Mdst_gamma and Mdst_ecl to let people get mass-constraint fitted momentum of pi0 after ad_hoc correction.
void	make_pi0_primary_vertex (int, double, double, const HepPoint3D &, const CLHEP::HepSymMatrix &)
	Fill Mdst_pi0 based on the fit result.
void	correct_ecl_primary_vertex (const HepPoint3D &, const CLHEP::HepSymMatrix &)
	Correct ecl using primary vertex.
void	shift_tof (const int mode)
	Shift tof times to account for residuals. More...
void	shift_tof_set (const int expno, const int runno, const int mode, const int im, const double pmom, const double sgn, double &shift)
	Return time shifts for different exp.
double	vee_mass_nofit (const Belle::Mdst_vee2 &vee2, float scale=1.0)
	Calculates V0 mass with non-constraint fit results.
void	scale_momenta (float scale_data=1.0, float scale_mc=1.0, int mode=0)
	Scale momenta of Mdst_trk.
void	scale_momenta_set (const int, const int, const int, double &)
	Return scale factors for 2001 summer confs. More...
void	scale_momenta_set_v1 (const int, const int, const int, double &)
	Return scale factors set_v1.
void	scale_momenta_set_v2 (const int, const int, const int, double &)
	Return scale factors set_v2.
void	scale_error (const int message_level=0)
	Apply scale error.
int	set_primary_vertex (HepPoint3D &v, CLHEP::HepSymMatrix &ve)
	Set primary vertex assuming all tracks are pions.
int	add_extra_trk_vee2 ()
	Add Mdst_trk_extra and Mdst_vee_extra to Mdst_trk and Mdst_vee2, respectively.
int	remove_extra_trk_vee2 ()
	Remove extra tracks from Mdst_trk and Mdst_vee2.
void	smear_trk ()
	Apply track smearing (on MC)

Detailed Description

Function Documentation

ecl_adhoc_corr()

static double Belle2::ecl_adhoc_corr ( int  Exp, int  Run, int  iflag05th, double  Energy, double    )

static

The function giving correction factor.

Corresponding Data/MC so that energy in data should be divided by this.

Definition at line 155 of file B2BIIFixMdstModule_ecl.cc.

  {
    double return_value;
 
    // Default.
    return_value = 1.0;
 
    double x = std::log10(Energy * 1000.0);
    switch (Exp) {
      case 7:
        //if( x < 2.5 )
        //{
        //   return_value = 0.85445 + 0.89162e-1*x
        //      -0.96202e-2*x*x -0.10944e-2*x*x*x;
        //}
        // New curve for b20020416 processed mdst.
        if (x < 2.6) {
          return_value = 0.85878 + 0.81947e-1 * x - 0.99708e-2 * x * x
                         - 0.28161e-3 * x * x * x;
        }
        break;
      case 9:
        //if( x < 3.0 )
        //{
        //   return_value = 0.83416 + 0.93534e-1*x
        //      -0.71632e-2*x*x -0.19072e-2*x*x*x;
        //}
        // New curve for b20020416 processed mdst.
        if (x < 3.2) {
          return_value = 0.86582 + 0.63194e-1 * x - 0.59391e-2 * x * x
                         + 0.17727e-2 * x * x * x - 0.62325e-3 * x * x * x * x;
        }
        break;
      case 11:
        // Proper curve for b20010507 processed mdst.
        //if( x < 3.0 )
        //{
        //   return_value = 0.84927 + 0.86308e-1*x
        //      -0.93898e-2*x*x -0.10557e-2*x*x*x;
        //}
        //
        // Proper curve for b20020416 processed mdst.
        if (x < 3.7) {
          if (x < 3.0) {
            return_value = 0.88989 + 0.41082e-1 * x
                           - 0.21919e-2 * x * x + 0.27116e-2 * x * x * x - 0.89113e-3 * x * x * x * x;
          } else { // Linear extrapolation to unity.
            return_value = 0.994441 + 4.90176e-3 * (x - 3.0);
          }
        }
        break;
      case 13:
        // Proper curve for b20010507 processed mdst.
        //if( x < 3.0 )
        //{
        //return_value = 0.78138 + 0.17899*x
        // -0.48821e-1*x*x +0.43692e-2*x*x*x;
        //}
        //
        // Proper curve for b20020416 processed mdst.
        if (x < 3.2) {
          return_value = 0.87434 + 0.62145e-1 * x
                         - 0.29691e-2 * x * x - 0.15843e-2 * x * x * x + 0.86858e-4 * x * x * x * x;
        }
        break;
      case 15:
        // Proper curve for b20020125 processed mdst.
        //if( x < 3.0 )
        //{
        //   return_value = 0.86858 + 0.75587e-1*x
        //      -0.11150e-1*x*x -0.21454e-4*x*x*x;
        //}
        //else if( x < 3.67 ) // Lenear extrapolation to be unity.
        //{
        //   return_value = 0.994412 + 8.11e-3*( x - 3.0 );
        //}
        // Proper curve for b20020416 processed mdst.
        if (x < 3.2) {
          return_value = 0.83073 + 0.10137 * x
                         - 0.59946e-2 * x * x - 0.56516e-2 * x * x * x + 0.87225e-3 * x * x * x * x;
        }
 
        break;
      case 17:
        // Proper curve for b20020416 processed mdst.
        if (x < 3.0) {
          return_value = 0.89260 + 0.54731e-1 * x
                         - 0.25736e-2 * x * x - 0.16493e-2 * x * x * x + 0.1032e-3 * x * x * x * x;
        } else if (x < 3.5) {
          return_value = 0.997459 + 0.0059039 * (x - 3.0);
        }
        break;
      case 19:
        // Correction curve obtained by Run001-299 HadronB
        // processed by b20020416.
        if (Run < 300) {
          if (x < 3.0) {
            return_value = 0.86432 + 0.87554e-1 * x
                           - 0.84182e-2 * x * x - 0.39880e-2 * x * x * x + 0.69435e-3 * x * x * x * x;
          }
        }
        if (299 < Run && Run < 900) { // For Run300-899 or later by b20020416.
          if (x < 2.9) {
            return_value = 0.92082 + 0.45896e-1 * x
                           - 0.68067e-2 * x * x + 0.94055e-3 * x * x * x - 0.27717e-3 * x * x * x * x;
          }
        }
        if (899 < Run) { // Till the end.
          if (x < 3.1) {
            return_value = 0.85154 + 0.97812e-1 * x
                           - 0.85774e-2 * x * x - 0.59092e-2 * x * x * x + 0.1121e-2 * x * x * x * x;
          }
        }
        break;
      case 21:
        // Proper curve for b20021205 processed mdst.
        if (x < 3.0) {
          return_value = 0.84940 + 0.86266e-1 * x
                         - 0.82204e-2 * x * x - 0.12912e-2 * x * x * x;
        }
        break;
      case 23:
        // Proper curve for b20021205 processed mdst.
        if (x < 3.1) {
          return_value = 0.85049 + 0.85418e-1 * x
                         - 0.45747e-2 * x * x - 0.40081e-2 * x * x * x + 0.52113e-3 * x * x * x * x;
        }
        break;
      case 25:
        if (x < 3.0) {
          return_value = 0.87001 + 0.73693e-1 * x
                         - 0.80094e-2 * x * x - 0.78527e-3 * x * x * x + 0.25888e-4 * x * x * x * x;
        }
        break;
      case 27: // It turned out no signif. change through all the runs.
        if (x < 2.85) {
          return_value = 0.90051 + 0.56094e-1 * x
                         - 0.14842e-1 * x * x + 0.55555e-2 * x * x * x - 0.10378e-2 * x * x * x * x;
        }
        break;
      case 31:
        if (iflag05th == 1) { // For 0.5MeV threshold.
          // For b20040429.
          if (x < 2.95) {
            return_value = 0.80295 + 0.13395 * x
                           - 0.10773e-1 * x * x - 0.85861e-2 * x * x * x + 0.15331e-2 * x * x * x * x;
          }
        } else { // For theta-dep. threshold.
          if (x < 3.25) {
            return_value = 0.57169 + 0.32548 * x
                           - 0.41157e-1 * x * x - 0.21971e-1 * x * x * x + 0.50114e-2 * x * x * x * x;
          }
        }
        break;
      case 33:
        if (iflag05th == 1) { // For 0.5MeV threshold.
          // For b20040429.
          if (x < 2.64) {
            return_value = 0.81153 + 0.10847 * x
                           - 0.24652e-2 * x * x - 0.54738e-2 * x * x * x + 0.41243e-3 * x * x * x * x;
          }
        } else { // For theta-dep. threshold, b20090127.
          if (x < 3.15) {
            return_value = 0.59869 + 0.29276 * x
                           - 0.15849e-1 * x * x - 0.31322e-1 * x * x * x + 0.62491e-2 * x * x * x * x;
          }
        }
        break;
      case 35:
        if (iflag05th == 1) { // For 0.5MeV threshold.
          // For b20040429.
          if (x < 2.54) {
            return_value = 0.83528 + 0.10402 * x
                           - 0.62047e-2 * x * x - 0.62411e-2 * x * x * x + 0.10312e-2 * x * x * x * x;
          }
        } else { // For theta-dep. threshold, b20090127.
          if (x < 3.2) {
            return_value =  0.58155 + 0.30642 * x
                            - 0.16981e-1 * x * x - 0.32609e-1 * x * x * x + 0.64874e-2 * x * x * x * x;
          }
        }
        break;
      case 37:
        if (iflag05th == 1) { // For 0.5MeV threshold.
          if (x < 2.5) {
            return_value = 0.83706 + 0.10726 * x
                           - 0.62423e-2 * x * x - 0.42425e-2 * x * x * x;
          }
        } else { // For theta-dep. threshold.
          if (x < 3.15) {
            return_value =  0.58801 + 0.30569 * x
                            - 0.18832e-1 * x * x - 0.32116e-1 * x * x * x + 0.64899e-2 * x * x * x * x;
          }
        }
        break;
      case 39:
        if (iflag05th == 1) { // For 0.5MeV threshold.
          if (x < 3.0) {
            if (x < 2.8) {
              return_value = 0.80827 + 0.095353 * x
                             - 0.14818e-2 * x * x - 0.23854e-2 * x * x * x - 0.22454e-3 * x * x * x * x;
            } else {
              return_value = 0.0112468 * (x - 2.8) + 0.997847;
            }
          }
        } else { // For theta-dep. threshold.
          if (x < 3.25) {
            return_value = 0.61133 + 0.27115 * x
                           - 0.12724e-1 * x * x - 0.28167e-1 * x * x * x + 0.54699e-2 * x * x * x * x;
          }
        }
        break;
      case 41:
        if (iflag05th == 1) { // For 0.5MeV threshold.
          if (Run < 800) {
            if (x < 2.5) {
              return_value = 0.90188 + 0.76563e-1 * x
                             - 0.16328e-1 * x * x + 0.56816e-3 * x * x * x;
            }
          }
          if (799 < Run) { // Run1200-1261 no signif. change from Run800-1199.
            if (x < 2.95) {
              return_value = 0.88077 + 0.71720e-1 * x
                             - 0.92197e-2 * x * x - 0.15932e-2 * x * x * x + 0.38133e-3 * x * x * x * x;
            }
          }
        } else { // For theta-dep. threshold.
          if (x < 3.2) {
            return_value =  0.57808 + 0.31083 * x
                            - 0.20247e-1 * x * x - 0.31658e-1 * x * x * x + 0.64087e-2 * x * x * x * x;
          }
        }
        break;
      case 43: // modified 20090829.
        if (iflag05th == 1) { // For 0.5MeV threshold.
          if (Run < 300) {
            if (x < 3.3) {
              return_value = 0.85592 + 0.93352e-1 * x
                             - 0.93144e-2 * x * x - 0.43681e-2 * x * x * x + 0.7971e-3 * x * x * x * x;
            }
          }
          if (299 < Run && Run < 700) { // added 20050630.
            if (x < 3.0) {
              return_value = 0.89169 + 0.73301e-1 * x
                             - 0.13856e-1 * x * x + 0.47303e-3 * x * x * x;
            }
          }
          if (699 < Run) { // added 20050727
            if (x < 3.1) {
              return_value = 0.90799 + 0.69815e-1 * x
                             - 0.17490e-1 * x * x + 0.14651e-2 * x * x * x;
            }
          }
        } else { // For theta-dep. threshold b20090127.
          if (x < 3.25) {
            return_value = 0.58176 + 0.30245 * x - 0.16390e-1 * x * x
                           - 0.32258e-1 * x * x * x + 0.64121e-2 * x * x * x * x;
          }
          //if( 999 < Run ){// For 5S runs, 2/fb, obsolete 20091209.
          //if( x < 3.3 ){
          //return_value = 0.58463 + 0.29780*x -0.14441e-1*x*x
          //-0.32323e-1*x*x*x + 0.63607e-2*x*x*x*x;
          //}
          //}
        }
        break;
      case 45:
        if (iflag05th == 1) { // For 0.5MeV threshold.
          if (x < 3.1) { // added 20060413
            return_value = 0.84823 + 0.10394 * x
                           - 0.92233e-2 * x * x - 0.72700e-2 * x * x * x + 0.14552e-2 * x * x * x * x;
          }
        } else { // For theta-dep. threshold b20090127.
          if (x < 3.25) {
            return_value = 0.62797 + 0.23897 * x
                           + 0.10261e-1 * x * x - 0.35700e-1 * x * x * x + 0.63846e-2 * x * x * x * x;
          }
        }
        break;
      case 47:
        if (iflag05th == 1) { // For 0.5MeV threshold.
          if (x < 3.1) { // added 20060420
            return_value = 0.86321 + 0.82987e-1 * x
                           - 0.12139e-1 * x * x - 0.12920e-3 * x * x * x;
          }
        } else { // For theta-dep. threshold b20090127.
          if (x < 3.25) {
            return_value = 0.58061 + 0.30399 * x
                           - 0.17520e-1 * x * x - 0.31559e-1 * x * x * x + 0.62681e-2 * x * x * x * x;
          }
        }
        break;
      case 49:
        if (iflag05th == 1) { // For 0.5MeV threshold.
          if (x < 3.13) { // added 20060511
            if (x < 2.80) {
              return_value = 0.85616 + 0.78041e-1 * x
                             - 0.38987e-2 * x * x - 0.31224e-2 * x * x * x + 0.33374e-3 * x * x * x * x;
            } else {
              return_value = ((1.0 - 0.99608) / (3.13 - 2.80)) * (x - 2.80) + 0.99608;
            }
          }
        } else { // For theta-dep. threshold b20090127.
          if (x < 3.3) {
            return_value = 0.57687 + 0.29948 * x
                           - 0.10594e-1 * x * x - 0.34561e-1 * x * x * x + 0.67064e-2 * x * x * x * x;
          }
        }
        break;
      case 51:
        if (iflag05th == 1) { // For 0.5MeV threshold, added 20070323.
          if (x < 3.07) {
            return_value = 0.89063 + 0.72152e-1 * x
                           - 0.14143e-1 * x * x + 0.73116e-3 * x * x * x;
          }
        } else { // For theta-dep. threshold b20090127.
          if (x < 3.3) {
            return_value = 0.59310 + 0.28618 * x
                           - 0.13858e-1 * x * x - 0.30230e-1 * x * x * x + 0.59173e-2 * x * x * x * x;
          }
        }
        break;
      case 53: // modified 20090807.
        if (iflag05th == 1) { // For 0.5MeV threshold.
          if (x < 2.63) {
            //std::cout<<"Exp.53 0.5MeV threshold."<<std::endl;
            return_value = 0.83425 + 0.10468 * x
                           - 0.53641e-2 * x * x - 0.60276e-2 * x * x * x + 0.77763e-3 * x * x * x * x;
          }
        } else { // For theta-dep. threshold b20090127.
          if (x < 3.1) {
            //std::cout<<"Exp.53 theta-dep. threshold."<<std::endl;
            return_value = 0.66100 + 0.25192 * x
                           - 0.16419e-1 * x * x - 0.25720e-1 * x * x * x + 0.52189e-2 * x * x * x * x;
          }
        }
        break;
      case 55:
        if (iflag05th == 1) { // For 0.5MeV threshold.
          if (x < 3.3) {
            return_value = 0.86202 + 0.79575e-1 * x
                           - 0.66721e-2 * x * x - 0.28609e-2 * x * x * x + 0.42784e-3 * x * x * x * x;
          }
        } else { // For theta-dep. threshold b20090127.
          if (x < 3.25) {
            return_value = 0.58789 + 0.29310 * x
                           - 0.15784e-1 * x * x - 0.30619e-1 * x * x * x + 0.60648e-2 * x * x * x * x;
          }
        }
        break;
      case 61:
        if (Run < 1210) { // On 4S by Th.-dep. threshold.
          if (x < 3.5) {
            return_value = 0.68839 + 0.18218 * x
                           - 0.23140e-2 * x * x - 0.17439e-1 * x * x * x + 0.29960e-2 * x * x * x * x;
          }
        }
        if (Run > 1209) { // 5S scan.
          if (iflag05th == 1) { // processed by 0.5MeV threshold.
            if (x < 3.3) {
              if (x < 2.8) {
                return_value = 0.94294 - 0.77497e-2 * x
                               - 0.43464e-3 * x * x + 0.99837e-2 * x * x * x - 0.23726e-2 * x * x * x * x;
              } else {
                return_value = 0.0162997 * (x - 2.80) + 0.991162;
              }
            }
          } else { // 5S scan. by Th.-dep. threshold.
            if (x < 3.5) {
              return_value = 0.64816 + 0.22492 * x
                             - 0.91745e-2 * x * x - 0.21736e-1 * x * x * x + 0.41333e-2 * x * x * x * x;
            }
          }
        }
        break;
      case 63:
        if (x < 3.4) {
          return_value = 0.69302 + 0.18393 * x
                         - 0.10983e-1 * x * x - 0.14148e-1 * x * x * x + 0.27298e-2 * x * x * x * x;
        }
        break;
      case 65:
        if (Run > 1000) { // 1S run turned out to be same as 5S scan.
          if (iflag05th == 1) { // For 0.5MeV threshold.
            if (x < 3.3) {
              if (x < 2.8) {
                return_value = 0.94294 - 0.77497e-2 * x
                               - 0.43464e-3 * x * x + 0.99837e-2 * x * x * x - 0.23726e-2 * x * x * x * x;
              } else {
                return_value = 0.0162997 * (x - 2.80) + 0.991162;
              }
            }
          } else { // For theta-dep. threshold, b20090127.
            if (x < 3.6) {
              return_value = 0.70987 + 0.16230 * x
                             - 0.64867e-2 * x * x - 0.12021e-1 * x * x * x + 0.20874e-2 * x * x * x * x;
            }
          }
        } else { // 4S runs processed with new lib.
          if (x < 3.4) {
            return_value = 0.66833 + 0.20602 * x
                           - 0.14322e-1 * x * x - 0.15712e-1 * x * x * x + 0.31114e-2 * x * x * x * x;
          }
        }
        break;
      case 67:
        if (x < 3.4) {
          return_value = 0.64196 + 0.23069 * x
                         - 0.20303e-1 * x * x - 0.15680e-1 * x * x * x + 0.31611e-2 * x * x * x * x;
        }
        break;
      case 69:
        if (x < 3.35) {
          //return_value = 0.64196 + 0.23752*x
          //-0.85197e-2*x*x -0.24366e-1*x*x*x +0.46397e-2*x*x*x*x;
          return_value = 0.99 * (0.64196 + 0.23752 * x
                                 - 0.85197e-2 * x * x - 0.24366e-1 * x * x * x
                                 + 0.46397e-2 * x * x * x * x);
        }
        break;
      case 71:
        if (x < 3.35) {
          return_value = 0.66541 + 0.12579 * x
                         + 0.89999e-1 * x * x - 0.58305e-1 * x * x * x + 0.86969e-2 * x * x * x * x;
        }
        break;
      case 73:
        if (x < 3.45) {
          return_value = 0.62368 + 0.24142 * x
                         - 0.11677e-1 * x * x - 0.22477e-1 * x * x * x + 0.42765e-2 * x * x * x * x;
        }
        break;
    }
    // This return value has to be used as division.
    // i.e. E_corrected = E_original/return_vale.
    return return_value;
  }

pi0resol()

double pi0resol ( double  p, double  theta, const char *  side, bool  mcdata, int  , int    )

staticprivate

Treat pi0 mass width as a func.

of pi0 momentum.

Definition at line 685 of file B2BIIFixMdstModule_ecl.cc.

scale_error_impl()

static int Belle2::scale_error_impl ( const int  message_level, const int  reprocess_version  )

static

The implementation of B2BIIFixMdstModule::scale_error()

Returns

0 if scale error is properly applied. 1 if scale error is already applied. 2 if scale error for this exp is not available.

See also

B2BIIFixMdstModule::scale_error()

Definition at line 2775 of file B2BIIFixMdstModule_trk.cc.

  {
    SE_Message_Level = message_level;
    SE_Reprocess_Version = reprocess_version;
 
    int no_exp{0}, no_run{0}, no_evt{0}, no_frm{0}, expmc{0};
    get_event_id(&no_exp, &no_run, &no_evt, &no_frm, &expmc);
 
    if (DUMMY_scale == get_scale_error_func_for_exprun(no_exp, no_run)) return 2;
    if (is_already_scaled()) return 1;
 
    static int first_exp73 = true;
    if (no_exp == 73 && first_exp73) {
      first_exp73 = false;
      B2ERROR(
        "scale_error: warning: scale parameters for exp#71 are tentatively used for exp#73.");
    }
 
    double scale[5] = { // ... temporary
      1.0, // dr
      1.0, // phi0
      1.0, // kappa
      1.0, // dz
      1.0, // tanl
    };
 
 
    // scale error matrices in mdst_trk_fit
    Belle::Mdst_trk_fit_Manager& fitmgr = Belle::Mdst_trk_fit_Manager::get_manager();
    for (std::vector<Belle::Mdst_trk_fit>::iterator it = fitmgr.begin(); it != fitmgr.end(); ++it) {
      Belle::Mdst_trk_fit& fit = *it;
      if (fit.helix(2) == 0.) continue;
 
      const double pt   = 1. / fabs(fit.helix(2));
      const double tanl = fit.helix(4);
 
      cal_scale_error(scale, pt, tanl, expmc, no_exp, no_run);
 
      fit.error(0, scale[0]*scale[0]*fit.error(0));
      fit.error(1, scale[1]*scale[0]*fit.error(1));
      fit.error(2, scale[1]*scale[1]*fit.error(2));
      fit.error(3, scale[2]*scale[0]*fit.error(3));
      fit.error(4, scale[2]*scale[1]*fit.error(4));
      fit.error(5, scale[2]*scale[2]*fit.error(5));
      fit.error(6, scale[3]*scale[0]*fit.error(6));
      fit.error(7, scale[3]*scale[1]*fit.error(7));
      fit.error(8, scale[3]*scale[2]*fit.error(8));
      fit.error(9, scale[3]*scale[3]*fit.error(9));
      fit.error(10, scale[4]*scale[0]*fit.error(10));
      fit.error(11, scale[4]*scale[1]*fit.error(11));
      fit.error(12, scale[4]*scale[2]*fit.error(12));
      fit.error(13, scale[4]*scale[3]*fit.error(13));
      fit.error(14, scale[4]*scale[4]*fit.error(14));
    }
 
 
    // scale error matrices in mdst_daughters
    Belle::Mdst_vee_daughters_Manager& daumgr = Belle::Mdst_vee_daughters_Manager::get_manager();
    for (std::vector<Belle::Mdst_vee_daughters>::iterator it = daumgr.begin(); it != daumgr.end(); ++it) {
      Belle::Mdst_vee_daughters& dau = *it;
      if (dau.helix_p(2) == 0. || dau.helix_m(2) == 0.) continue;
 
      // positive track
      const double pt_p   = 1. / fabs(dau.helix_p(2));
      const double tanl_p = dau.helix_p(4);
 
      cal_scale_error(scale, pt_p, tanl_p, expmc, no_exp, no_run);
 
      dau.error_p(0, scale[0]*scale[0]*dau.error_p(0));
      dau.error_p(1, scale[1]*scale[0]*dau.error_p(1));
      dau.error_p(2, scale[1]*scale[1]*dau.error_p(2));
      dau.error_p(3, scale[2]*scale[0]*dau.error_p(3));
      dau.error_p(4, scale[2]*scale[1]*dau.error_p(4));
      dau.error_p(5, scale[2]*scale[2]*dau.error_p(5));
      dau.error_p(6, scale[3]*scale[0]*dau.error_p(6));
      dau.error_p(7, scale[3]*scale[1]*dau.error_p(7));
      dau.error_p(8, scale[3]*scale[2]*dau.error_p(8));
      dau.error_p(9, scale[3]*scale[3]*dau.error_p(9));
      dau.error_p(10, scale[4]*scale[0]*dau.error_p(10));
      dau.error_p(11, scale[4]*scale[1]*dau.error_p(11));
      dau.error_p(12, scale[4]*scale[2]*dau.error_p(12));
      dau.error_p(13, scale[4]*scale[3]*dau.error_p(13));
      dau.error_p(14, scale[4]*scale[4]*dau.error_p(14));
 
      // negative track
      const double pt_m   = 1. / fabs(dau.helix_m(2));
      const double tanl_m = dau.helix_m(4);
 
      cal_scale_error(scale, pt_m, tanl_m, expmc, no_exp, no_run);
 
      dau.error_m(0, scale[0]*scale[0]*dau.error_m(0));
      dau.error_m(1, scale[1]*scale[0]*dau.error_m(1));
      dau.error_m(2, scale[1]*scale[1]*dau.error_m(2));
      dau.error_m(3, scale[2]*scale[0]*dau.error_m(3));
      dau.error_m(4, scale[2]*scale[1]*dau.error_m(4));
      dau.error_m(5, scale[2]*scale[2]*dau.error_m(5));
      dau.error_m(6, scale[3]*scale[0]*dau.error_m(6));
      dau.error_m(7, scale[3]*scale[1]*dau.error_m(7));
      dau.error_m(8, scale[3]*scale[2]*dau.error_m(8));
      dau.error_m(9, scale[3]*scale[3]*dau.error_m(9));
      dau.error_m(10, scale[4]*scale[0]*dau.error_m(10));
      dau.error_m(11, scale[4]*scale[1]*dau.error_m(11));
      dau.error_m(12, scale[4]*scale[2]*dau.error_m(12));
      dau.error_m(13, scale[4]*scale[3]*dau.error_m(13));
      dau.error_m(14, scale[4]*scale[4]*dau.error_m(14));
    }
 
    return 0;
  }

scale_momenta_set()

void scale_momenta_set ( const int  mode, const int  expno, const int  runno, double &  scale  )

private

Return scale factors for 2001 summer confs.

analyses, only for exp < 17

Definition at line 406 of file B2BIIFixMdstModule_trk.cc.

shift_tof()

void shift_tof ( const int  mode )

private

Shift tof times to account for residuals.

Based on scale_momenta code

Definition at line 59 of file B2BIIFixMdstModule_tof.cc.

  {
//===============================================================
// Shifts tof times to account for residuals.
// Based on scale_momenta code
//======================================================
//Check existence of belle_event
 
    Belle::Belle_event_Manager& evtmgr = Belle::Belle_event_Manager::get_manager();
    if (evtmgr.count() == 0) return;  //do nothing if no Belle::Belle_Event
 
    Belle::Mdst_event_add_Manager& addmgr = Belle::Mdst_event_add_Manager::get_manager();
    if (addmgr.count() == 0) return;  //do nothing if no Belle::Mdst_Event_Add
    if (addmgr[0].flag(4) != 0) return; //do nothing if already shifted
 
//check mode
    if (mode < 0 || mode > 2) {
      B2ERROR("shift_tof: invalid mode specified;");
      B2ERROR("mode must be 0, 1, 2");
      return;
    }
 
    int expmc = evtmgr[0].ExpMC();
    int expno = evtmgr[0].ExpNo();
    int runno = evtmgr[0].RunNo();
 
    if (expmc == 2) return;  //for now, apply no shift to MC
 
    addmgr[0].flag(4, 1);    //set flag
 
    // Loop over all charged tracks
 
    Belle::Mdst_charged_Manager& ChgMgr = Belle::Mdst_charged_Manager::get_manager();
    for (std::vector<Belle::Mdst_charged>::iterator it = ChgMgr.begin();
         it != ChgMgr.end(); ++it) {
      Belle::Mdst_charged& Chg = *it;
      if (Chg) {
        // Get mdst_tof table for the track
        Belle::Mdst_tof& Tof = Chg.tof();
        if (Tof) {
          // Get momentum and charge of track
          TVector3 Mom(Chg.px(), Chg.py(), Chg.pz());
          double pmom = Mom.Mag();
          double sgn = Chg.charge();
          // Loop over mass hypotheses
          for (int im = 0; im < 5; im++) {
            double shift = 0.;
            shift_tof_set(expno, runno, mode, im, pmom, sgn, shift);
            if (fabs(shift) < .005) continue;
            double ot = Tof.tof_exp(im);              //old t (expected)
            double odt = ot - Tof.tof();              //old dt
            double opid = Tof.pid(im);                //old pid
//         double ocl=Tof.cl(im);                    //old cl
            int ndf = (int)(Tof.ndf(im) + .001);      //ndf
            if (opid > 0. && opid < 1.) {
              double ct = ot + shift;                //corr t
              double cdt = ct - Tof.tof();           //corr dt
              double err = fabs(odt) / sqrt(-2.*log(opid)); //est error
              double cch = pow(cdt / err, 2);        //corr chisq
              double cpid = exp(-.5 * cch);          //corr pid
              double ccl = TMath::Prob(cch, ndf);    //corr cl
              Tof.tof_exp(im, ct);
              Tof.pid(im, cpid);
              Tof.cl(im, ccl);
            }
          }
        }
      }
    }
  }

Variable Documentation

DUMMY_scale

const struct cal_scale_error_func_set_t DUMMY_scale

static

Initial value:

= {
    null_scale,
    null_scale,
    null_scale
  }

Dummy scale.

Definition at line 2615 of file B2BIIFixMdstModule_trk.cc.

EXP0723_scale

const struct cal_scale_error_func_set_t EXP0723_scale

static

Initial value:

= {
    cal_scale_error_EXP0723_hadronic_mc,
    cal_scale_error_EXP0723_cosmic_mc,
    cal_scale_error_EXP0723_cosmic_data,
  }

Scale error for Exp.7-23.

Definition at line 2497 of file B2BIIFixMdstModule_trk.cc.

EXP2527_scale

const struct cal_scale_error_func_set_t EXP2527_scale

static

Initial value:

= {
    cal_scale_error_EXP2527_hadronic_mc,
    cal_scale_error_EXP2527_cosmic_mc,
    cal_scale_error_EXP2527_cosmic_data
  }

Scale error for Exp.25-27.

Definition at line 2519 of file B2BIIFixMdstModule_trk.cc.

EXP31_scale

const struct cal_scale_error_func_set_t EXP31_scale

static

Initial value:

= {
    cal_scale_error_EXP31_hadronic_mc,
    cal_scale_error_EXP31_cosmic_mc,
    cal_scale_error_EXP31_cosmic_data
  }

Scale error for Exp.31.

Definition at line 2525 of file B2BIIFixMdstModule_trk.cc.

EXP33_scale

const struct cal_scale_error_func_set_t EXP33_scale

static

Initial value:

= {
    cal_scale_error_EXP33_hadronic_mc,
    cal_scale_error_EXP33_cosmic_mc,
    cal_scale_error_EXP33_cosmic_data
  }

Scale error for Exp.33.

Definition at line 2531 of file B2BIIFixMdstModule_trk.cc.

EXP35_scale

const struct cal_scale_error_func_set_t EXP35_scale

static

Initial value:

= {
    cal_scale_error_EXP35_hadronic_mc,
    cal_scale_error_EXP35_cosmic_mc,
    cal_scale_error_EXP35_cosmic_data
  }

Scale error for Exp.35.

Definition at line 2537 of file B2BIIFixMdstModule_trk.cc.

EXP37_scale

const struct cal_scale_error_func_set_t EXP37_scale

static

Initial value:

= {
    cal_scale_error_EXP37_hadronic_mc,
    cal_scale_error_EXP37_cosmic_mc,
    cal_scale_error_EXP37_cosmic_data
  }

Scale error for Exp.37.

Definition at line 2543 of file B2BIIFixMdstModule_trk.cc.

EXP3941_scale

const struct cal_scale_error_func_set_t EXP3941_scale

static

Initial value:

= {
    cal_scale_error_EXP3941_hadronic_mc,
    cal_scale_error_EXP3941_cosmic_mc,
    cal_scale_error_EXP3941_cosmic_data
  }

Scale error for Exp.39,41.

Definition at line 2549 of file B2BIIFixMdstModule_trk.cc.

EXP43_scale

const struct cal_scale_error_func_set_t EXP43_scale

static

Initial value:

= {
    cal_scale_error_EXP43_hadronic_mc,
    cal_scale_error_EXP43_cosmic_mc,
    cal_scale_error_EXP43_cosmic_data
  }

Scale error for Exp.43.

Definition at line 2555 of file B2BIIFixMdstModule_trk.cc.

EXP45_scale

const struct cal_scale_error_func_set_t EXP45_scale

static

Initial value:

= {
    cal_scale_error_EXP45_hadronic_mc,
    cal_scale_error_EXP45_cosmic_mc,
    cal_scale_error_EXP4547_cosmic_data
  }

Scale error for Exp.45.

Definition at line 2561 of file B2BIIFixMdstModule_trk.cc.

EXP47_scale

const struct cal_scale_error_func_set_t EXP47_scale

static

Initial value:

= {
    cal_scale_error_EXP47_hadronic_mc,
    cal_scale_error_EXP47_cosmic_mc,
    cal_scale_error_EXP4547_cosmic_data
  }

Scale error for Exp.47.

Definition at line 2567 of file B2BIIFixMdstModule_trk.cc.

EXP49_scale

const struct cal_scale_error_func_set_t EXP49_scale

static

Initial value:

= {
    cal_scale_error_EXP49_hadronic_mc,
    cal_scale_error_EXP49_cosmic_mc,
    cal_scale_error_EXP49_cosmic_data
  }

Scale error for Exp.49.

Definition at line 2573 of file B2BIIFixMdstModule_trk.cc.

EXP51_scale

const struct cal_scale_error_func_set_t EXP51_scale

static

Initial value:

= {
    cal_scale_error_EXP51_hadronic_mc,
    cal_scale_error_EXP51_cosmic_mc,
    cal_scale_error_EXP51_cosmic_data
  }

Scale error for Exp.51.

Definition at line 2579 of file B2BIIFixMdstModule_trk.cc.

EXP53_scale

const struct cal_scale_error_func_set_t EXP53_scale

static

Initial value:

= {
    cal_scale_error_EXP53_hadronic_mc,
    cal_scale_error_EXP53_cosmic_mc,
    cal_scale_error_EXP53_cosmic_data
  }

Scale error for Exp.53.

Definition at line 2585 of file B2BIIFixMdstModule_trk.cc.

EXP55_scale

const struct cal_scale_error_func_set_t EXP55_scale

static

Initial value:

= {
    cal_scale_error_EXP55_hadronic_mc,
    cal_scale_error_EXP55_cosmic_mc,
    cal_scale_error_EXP55_cosmic_data
  }

Scale error for Exp.55.

Definition at line 2591 of file B2BIIFixMdstModule_trk.cc.

EXP6165_scale

const struct cal_scale_error_func_set_t EXP6165_scale

static

Initial value:

= {
    cal_scale_error_EXP6165_hadronic_mc,
    cal_scale_error_EXP6165_cosmic_mc,
    cal_scale_error_EXP6165_cosmic_data
  }

Scale error for Exp.61-65.

Definition at line 2597 of file B2BIIFixMdstModule_trk.cc.

EXP67_scale

const struct cal_scale_error_func_set_t EXP67_scale

static

Initial value:

= {
    cal_scale_error_EXP67_hadronic_mc,
    cal_scale_error_EXP67_cosmic_mc,
    cal_scale_error_EXP67_cosmic_data
  }

Scale error for Exp.67.

Definition at line 2603 of file B2BIIFixMdstModule_trk.cc.

EXP6971_scale

const struct cal_scale_error_func_set_t EXP6971_scale

static

Initial value:

= {
    cal_scale_error_EXP6971_hadronic_mc,
    cal_scale_error_EXP6971_cosmic_mc,
    cal_scale_error_EXP6971_cosmic_data
  }

Scale error for Exp.69-71.

Definition at line 2609 of file B2BIIFixMdstModule_trk.cc.

SE_Message_Level

int SE_Message_Level

static

Message level of scale_error().

0: none, 1: warning, 2: info

Definition at line 1434 of file B2BIIFixMdstModule_trk.cc.