tests/src/MapsSensorTest.cpp

Go to the documentation of this file.

#include "MapsSensor.hh"
#include "TRandom2.h"
#include "TFile.h"
#include <iostream>
int main() {
        MapsSensor* mySensor = new MapsSensor(0, 0);
        TRandom2 t2;
        unsigned thresh(100);
        for (unsigned h(0); h < 100000; h++) {
                double x, y, z;
                t2.Sphere(x, y, z, 1);
                double t(t2.Uniform());
                MapsSensor::coord c(x*168, y*168);
                if (z > 0) {
                        mySensor->registerTrackConfirm(thresh, c, static_cast<unsigned>(t
                                        *8000));

                } else {
                        std::cout << __PRETTY_FUNCTION__ << __PRETTY_LINE__ << ": disabled!\n
                        //mySensor->registerTrackMiss(thresh, static_cast<unsigned>(t*8000));

                }
        }
        
        thresh = 130;
        for (unsigned h(0); h < 100000; h++) {
                double x, y, z;
                t2.Sphere(x, y, z, 1);
                double t(t2.Uniform());
                MapsSensor::coord c(x*168, y*168);
                if (z > 0) {
                        mySensor->registerTrackConfirm(thresh, c, static_cast<unsigned>(t
                                        *8000));

                } else {
                        std::cout << __PRETTY_FUNCTION__ << __PRETTY_LINE__ << ": disabled!\n
                        //mySensor->registerTrackMiss(thresh, static_cast<unsigned>(t*8000));

                }
        }
        
        
        thresh = 160;
        for (unsigned h(0); h < 100000; h++) {
                double x, y, z;
                t2.Sphere(x, y, z, 1);
                double t(t2.Uniform());
                MapsSensor::coord c(x*168, y*168);
                if (z > 0) {
                        mySensor->registerTrackConfirm(thresh, c, static_cast<unsigned>(t
                                        *8000));

                } else {
                        std::cout << __PRETTY_FUNCTION__ << __PRETTY_LINE__ << ": disabled!\n";
                        //mySensor->registerTrackMiss(thresh, static_cast<unsigned>(t*8000));

                }
                std::pair<double, double> resid(t2.Gaus(0, 5), t2.Gaus(2, 5));
                mySensor->setResidual(resid);
                
        }
        
        TFile* file = new TFile("TestMapsSensor.root", "recreate");
        std::cout << "\n" << *mySensor << " efficiency at th = 100: \t"
                        << mySensor->getEfficiency(thresh) << "\n";
        TH1F effTimes;
        std::cout << "Calling get efficiency timestamps...\n";
        mySensor->getEfficiencyTimestamps(effTimes);
        std::cout << "Done!\n";
        TH1F ineffTimes;
        mySensor->getInefficiencyTimestamps(ineffTimes);
        effTimes.Write();
        ineffTimes.Write();
        
        std::cout << "Getting efficiency curve...\n";
        TH1F efficiency;
        mySensor->getEfficiencyCurve(efficiency, 10, 100, 200);
        efficiency.Write();
        
        std::cout << "Getting residual curve...\n";
        TH2F residuals;
        mySensor->getResidualXYPlot(residuals);
        residuals.Write();
        
        std::cout << "Getting efficiency xy plot...\n";
        TH2F effxy;
        mySensor->getEfficiencyXYPlot(effxy);
        effxy.Write();
        

        std::cout << "Getting inefficiency xy plot...\n";
        TH2F ineffxy;
        mySensor->getInefficiencyXYPlot(ineffxy);
        ineffxy.Write();
        
        file->Write();
        file->Close();
}

---