MUSIC unknown
Analysis for the MUSIC active-target ionization chamber
Loading...
Searching...
No Matches
RegionCuts.cpp File Reference
#include "RegionCuts.hpp"
#include "Paths.hpp"
#include "PlottingUtils.hpp"
#include <TAxis.h>
#include <TCanvas.h>
#include <TFile.h>
#include <TH1D.h>
#include <TMath.h>
#include <TParameter.h>
#include <TSystem.h>
#include <cmath>
#include <iostream>
#include <vector>
Include dependency graph for RegionCuts.cpp:

Go to the source code of this file.

Namespaces

namespace  RegionCutStore
 On-disk home of the per-strip (a,n) and (a,a') region cuts.
namespace  RegionCutFinder
 Finding regions by fitting a two-component mixture to the scatter.

Functions

TString RegionCutStore::Dir ()
 Directory holding the region cut files.
TString RegionCutStore::Path (const char *name, Int_t reac)
 Path to one cut's file.
void RegionCutStore::Save (Int_t reac, TCutG *cut_an, TCutG *cut_aa, Double_t n_an_assigned=-1.0)
 Write a strip's two region cuts.
void RegionCutStore::SaveFit (Int_t reac, const RegionFit &fit)
 Store the fitted components beside the (a,n) cut.
Bool_t RegionCutStore::LoadFit (Int_t reac, RegionFit &fit)
 Read back a stored fit.
Double_t RegionCutStore::LoadAssigned (const char *name, Int_t reac)
 The attributed count stored with a cut.
TCutG * RegionCutStore::Load (const char *name, Int_t reac)
 Load a cut, from either storage generation.
TCutG * RegionCutStore::LoadDrawn (const char *name, Int_t reac)
 Load only a hand-drawn cut, never a fitted one.
RegionFit RegionCutFinder::FitMixture (TH2F *scatter, Int_t reac, Double_t x_lo, Double_t x_hi, Double_t y_lo, Double_t y_hi)
 Fit the two-component mixture to one strip's scatter.
Double_t RegionCutFinder::AboveRidge (const Gauss2D &beam, Double_t x, Double_t y)
 How far a point sits above the beam ridge, in conditional sigma.
RegionFit RegionCutFinder::FitBeam (TH2F *scatter, Int_t reac, Double_t x_lo, Double_t x_hi, Double_t y_lo, Double_t y_hi)
 Fit the beam component alone, with no reaction component.
TCutG * RegionCutFinder::RidgeBandCut (const char *name, const Gauss2D &beam, Double_t nsig_lo, Double_t nsig_hi, Double_t x_lo, Double_t x_hi, Double_t y_lo, Double_t y_hi)
 Closed polygon of a band above the beam ridge.
TCutG * RegionCutFinder::EllipseCut (const char *name, const Gauss2D &g, Double_t nsigma, Int_t npts=64)
 Closed polygon of a component's Mahalanobis contour.
Double_t RegionCutFinder::CountInside (TH2F *scatter, TCutG *cut, const RegionFit &fit)
 Count scatter events inside a cut.
void RegionCutFinder::SaveFigures (TH2F *scatter, Int_t reac, const RegionFit &fit, TCutG *an, TCutG *aa, const TString &subdir)
 Save the diagnostic figures for a strip's regions.