00001 void bg(int i, double p){
00002 double mass = 0;
00003 switch(i){
00004 case 0:
00005 mass = 0.511e-3;
00006 cout << "mass " << mass << " betagamma " << p/mass << endl;
00007 break;
00008 case 1:
00009 mass = 105.658e-3;
00010 cout << "mass " << mass << " betagamma " << p/mass << endl;
00011 break;
00012 case 2:
00013 mass = 139.570e-3;
00014 cout << "mass " << mass << " betagamma " << p/mass << endl;
00015 break;
00016 case 3:
00017 mass = 493.677e-3;
00018 cout << "mass " << mass << " betagamma " << p/mass << endl;
00019 break;
00020 case 4:
00021 mass = 938.272e-3;
00022 cout << "mass " << mass << " betagamma " << p/mass << endl;
00023 break;
00024 default:
00025 cout << "bg(index, momentum) index:0/1/2/3/4; e/mu/pi/k/p" << endl;
00026 }
00027 }
00028
00029
00030 void rebg(int i, double bg){
00031 double mass = 0;
00032 switch(i){
00033 case 0:
00034 mass = 0.511e-3;
00035 cout << "mass " << mass << " momentum " << bg*mass << endl;
00036 break;
00037 case 1:
00038 mass = 105.658e-3;
00039 cout << "mass " << mass << " momentum " << bg*mass << endl;
00040 break;
00041 case 2:
00042 mass = 139.570e-3;
00043 cout << "mass " << mass << " momentum " << bg*mass << endl;
00044 break;
00045 case 3:
00046 mass = 493.677e-3;
00047 cout << "mass " << mass << " momentum " << bg*mass << endl;
00048 break;
00049 case 4:
00050 mass = 938.272e-3;
00051 cout << "mass " << mass << " momentum " << bg*mass << endl;
00052 break;
00053 default:
00054 cout << "bg(index, betagamma) index:0/1/2/3/4; e/mu/pi/k/p" << endl;
00055 }
00056 }
