File indexing completed on 2024-04-06 12:22:46
0001 #include <iostream>
0002
0003 #include "MuonAnalysis/MomentumScaleCalibration/interface/MuScleFitUtils.h"
0004
0005 int main()
0006 {
0007 std::cout << "Testing MuScleFitUtils" << std::endl;
0008 std::cout << std::endl;
0009
0010 int iY = 0;
0011 double mass = 90.;
0012 double massResol = 1.;
0013 int iRes = 0;
0014
0015
0016 for (int iy=0; iy<=MuScleFitUtils::nbins; ++iy) {
0017 MuScleFitUtils::GLZNorm[iY][iy] = 0.;
0018 for (int ix=0; ix<=MuScleFitUtils::nbins; ++ix) {
0019 MuScleFitUtils::GLZValue[iY][ix][iy] = 1.;
0020 MuScleFitUtils::GLZNorm[iY][iy] += MuScleFitUtils::GLZValue[iY][ix][iy];
0021 }
0022 }
0023
0024 double prob = MuScleFitUtils::probability( mass, massResol,
0025 MuScleFitUtils::GLZValue, MuScleFitUtils::GLZNorm,
0026 iRes, iY );
0027
0028 std::cout << "Probability = " << prob << std::endl;
0029
0030 }