testMuonArrivalTime.cc

Go to the documentation of this file.

#include <iostream>
#include "MuonArrivalTime.h"

#include <tst/Verify.h>
#include <cppunit/extensions/HelperMacros.h>

using namespace std;
using namespace utl;
using namespace tst;

class MuonArrivalTimeTest : public CppUnit::TestFixture {

  CPPUNIT_TEST_SUITE(MuonArrivalTimeTest);
  CPPUNIT_TEST(testTimeCDF);
  CPPUNIT_TEST(testApproximateTimeCDF);
  CPPUNIT_TEST(testFirstTimePDF);
  CPPUNIT_TEST(testApproximateFirstTimePDF);
  CPPUNIT_TEST_SUITE_END();

  MuonArrivalTime fTM;

public:

  void
  setUp()
  {
    fTM.SetTheta(60*degree);
    fTM.SetCoordinates(1000*meter, 0);
  }

  void tearDown() { }

  void
  testTimeCDF()
  {
    const double res = fTM.TimeCDF(1e9*ns);
    CPPUNIT_ASSERT(Verify<CloseTo>(res, 1.0));
  }

  void
  testApproximateTimeCDF()
  {
    const double res = fTM.ApproximateTimeCDF(1e9*ns);
    CPPUNIT_ASSERT(Verify<CloseTo>(res, 1.0));
  }

  void
  testFirstTimePDF()
  {
    const unsigned int nmuon = 10;
    const double lntmin = std::log(1e1*ns);
    const double lntmax = std::log(1e3*ns);
    const double n = 100;

    double tsum = 0;
    double ttsum = 0;
    for (int i = 0; i < n; ++i)
    {
      const double z = double(i)/(n-1);
      const double t = std::exp((1.0-z)*lntmin + z*lntmax);
      const double dp = fTM.FirstTimePDF(t,nmuon)*t;
      tsum += dp*t;
      ttsum += dp*t*t;
    }

    const double dlnt = (lntmax-lntmin)/n;
    const double tmean = tsum*dlnt;
    const double tsigma = std::sqrt(ttsum*dlnt-tmean*tmean);

    CPPUNIT_ASSERT(Verify<CloseTo>(tmean, 76.*ns, 1e-1));
    CPPUNIT_ASSERT(Verify<CloseTo>(tsigma, 15.*ns, 1e-1));
  }

  void
  testApproximateFirstTimePDF()
  {
    const unsigned int nmuon = 10;
    const double lntmin = std::log(1e1*ns);
    const double lntmax = std::log(1e3*ns);
    const double n = 100;

    double tsum = 0;
    double ttsum = 0;
    for (int i = 0; i < n; ++i)
    {
      const double z = double(i)/(n-1);
      const double t = std::exp((1.0-z)*lntmin + z*lntmax);
      const double dp = fTM.ApproximateFirstTimePDF(t,nmuon)*t;
      tsum += dp*t;
      ttsum += dp*t*t;
    }

    const double dlnt = (lntmax-lntmin)/n;
    const double tmean = tsum*dlnt;
    const double tsigma = std::sqrt(ttsum*dlnt-tmean*tmean);

    CPPUNIT_ASSERT(Verify<CloseTo>(tmean, 76.*ns, 1e-1));
    CPPUNIT_ASSERT(Verify<CloseTo>(tsigma, 19.*ns, 1e-1));
  }

};


CPPUNIT_TEST_SUITE_REGISTRATION(MuonArrivalTimeTest);


// Configure (x)emacs for this file ...
// Local Variables:
// mode: c++
// End: