GetDetectorBins.cxx

Go to the documentation of this file.

#include "GetDetectorBins.h"

#include <UgliGeometry/UgliScintPlnHandle.h>
#include <UgliGeometry/UgliSteelPlnHandle.h>
#include <UgliGeometry/UgliStripHandle.h>

#include <Conventions/PlaneCoverage.h>
#include <Conventions/Detector.h>

#include <vector>
using namespace std;


void DataUtil::GetDetectorBinsUV(UgliGeomHandle ugh, PlaneView::EPlaneView view,
                                 int &nbins, double &min, double &max)
{
    vector<UgliScintPlnHandle> scints = ugh.GetScintPlnHandleVector();

    // u/v bins
    double binper=0;
    switch (ugh.GetVldRange().GetDetectorMask()) {
    case Detector::kNear:
        for (size_t ind = 0; ind < scints.size(); ++ind) {
            if (scints[ind].GetPlaneCoverage() == PlaneCoverage::kNearFull) {
                if (scints[ind].GetPlaneView() == view) {
                    vector<UgliStripHandle> v = scints[ind].GetStripHandleVector();
                    min = v.front().GetTPos()-v.front().GetHalfWidth();
                    max = v.back().GetTPos()+v.back().GetHalfWidth();
                    binper = v.size()/(max-min);
                }
                if (binper) break;
            }
        }
        break;
    case Detector::kFar: case Detector::kCalib: {
        vector<UgliStripHandle> v = scints[0].GetStripHandleVector();
        min = v.front().GetTPos()-v.front().GetHalfWidth();
        max =  v.back().GetTPos() +v.back().GetHalfWidth();
        binper = v.size()/(max-min);
        break;}
    default:
        break;
    }
    nbins = (int)(((max-min)*binper)+0.5);
}


void DataUtil::GetDetectorBinsZ(UgliGeomHandle ugh,
                                int &nbins, double &zmin, double &zmax)
{
    vector<UgliSteelPlnHandle> steel = ugh.GetSteelPlnHandleVector();

    // Z bins
    double zbinper=0;
    switch (ugh.GetVldRange().GetDetectorMask()) {
    case Detector::kNear: case Detector::kCalib:
        ugh.GetZExtent(zmin,zmax);
        zbinper = steel.size()/(zmax-zmin);
        break;
    case Detector::kFar:{
        double max,min;
        ugh.GetZExtent(zmin,min,0);
        ugh.GetZExtent(max,zmax,1);
        zbinper = steel.size()/((zmax-zmin)-(max-min));
        break;}
    default:
        break;
    }
    nbins = (int)(((zmax-zmin)*zbinper)+0.5);
}

---