TankResponseSpeed.cc

Go to the documentation of this file.

#include <tls/TankResponseFactory.h>
#include <tls/VTankResponse.h>

#include <iostream>
#include <cstdlib>
#include <cstdio>
#include <cmath>

#include <utl/Reader.h>
#include <utl/AugerUnits.h>
#include <utl/AugerException.h>
#include <utl/ErrorLogger.h>

#include <TRandom.h>
#include <TBenchmark.h>

using namespace tls;
using namespace utl;
using namespace std;

int main(int argc, char* argv[]) {

  const int nCycles = 10000;

  // Get the configuration. Which implementation is used is determined there.
  if (argc != 2)
  {
    cerr << "wrong number of arguments, need exactly one which is path to XML configuration file" << endl;
    exit(1);
  }

  Reader reader(argv[1]);
  Branch branch = reader.GetTopBranch().GetFirstChild();
  vector<Branch> branchList;
  while (branch) { branchList.push_back(branch); branch = branch.GetNextSibling(); }

  TBenchmark bench;
  
  cout << "No. of calls " << nCycles << endl;

  vector<double> thetaList;
  vector<double> signalList;

  for (int i=0;i<nCycles;++i)
  {
    thetaList.push_back(gRandom->Uniform(60.1*degree,87.9*degree));
    signalList.push_back(gRandom->Exp(10.));
  }
    
  cout << "Calls to PDF(...)" << endl;
  for (unsigned short ib=0;ib<branchList.size();++ib)
  {
    Branch& branch = branchList[ib];

    try {
      VTankResponse& tankResponse = TankResponseFactory::GetInstance().GetTankResponse(branch);
      ostringstream msg1;
      msg1 << branch.GetName() << "1";

      bench.Start(msg1.str().c_str());
      for (int i=0;i<nCycles;++i) tankResponse.PDF(signalList[i], thetaList[i], 0, 1);
      bench.Stop(msg1.str().c_str());
  
      ostringstream msg2;
      msg2 << branch.GetName() << "20";

      bench.Start(msg2.str().c_str());
      for (int i=0;i<nCycles;++i) tankResponse.PDF(signalList[i], thetaList[i], 0, 20);
      bench.Stop(msg2.str().c_str());
  
      char buf[1024];
      sprintf(buf,"%20s -  1 muon : real time %5.2f",branch.GetName().c_str(),bench.GetRealTime(msg1.str().c_str()));
      cout << buf << endl;
      sprintf(buf,"%20s - 20 muons: real time %5.2f",branch.GetName().c_str(),bench.GetRealTime(msg2.str().c_str()));
      cout << buf << endl;
    } catch (NonExistentComponentException& e) {
      ERROR(e.GetMessage());
    }
  }

  bench.Reset();
  cout << "Calls to CDF(...)" << endl;
  for (unsigned short ib=0;ib<branchList.size();++ib)
  {
    Branch& branch = branchList[ib];

    try {
      VTankResponse& tankResponse = TankResponseFactory::GetInstance().GetTankResponse(branch);

      ostringstream msg1;
      msg1 << branch.GetName() << "1";

      bench.Start(msg1.str().c_str());
      for (int i=0;i<nCycles;++i) tankResponse.CDF(signalList[i], thetaList[i], 0, 1);
      bench.Stop(msg1.str().c_str());
  
      ostringstream msg2;
      msg2 << branch.GetName() << "20";
      
      bench.Start(msg2.str().c_str());
      for (int i=0;i<nCycles;++i) tankResponse.CDF(signalList[i], thetaList[i], 0, 20);
      bench.Stop(msg2.str().c_str());
  
      char buf[1024];
      sprintf(buf,"%20s -  1 muon : real time %5.2f",branch.GetName().c_str(),bench.GetRealTime(msg1.str().c_str()));
      cout << buf << endl;
      sprintf(buf,"%20s - 20 muons: real time %5.2f",branch.GetName().c_str(),bench.GetRealTime(msg2.str().c_str()));
      cout << buf << endl;
    } catch (NonExistentComponentException& e) {
      ERROR(e.GetMessage());
    }
  }
  
  return 0;
}