TagEfficiency.hpp

Go to the documentation of this file.

#ifndef TAG_EFFICIENCY_HPP
#define TAG_EFFICIENCY_HPP
#include "RegionCuts.hpp"
#include <Rtypes.h>
#include <TCutG.h>
#include <TString.h>
#include <map>
#include <vector>
 
class TFile;
class TH2F;
class TTree;
struct CrossSectionChannel;

struct TagEfficiencyRecord {
  Int_t reac = -1; 
  Double_t n_counted = 0.0;
  Double_t eff = 0.0;
  Double_t eff_err = 0.0;
  Double_t migrate = 0.0;
  Double_t tag_eff = 0.0;
};

namespace TagEfficiencyStore {

TString Path();
void Write(const TString &channel,
           const std::vector<TagEfficiencyRecord> &records,
           const TString &method);
Bool_t Load(const TString &channel, Int_t reac, TagEfficiencyRecord &record);
TString Method();
} // namespace TagEfficiencyStore

class TagEfficiency {
public:
  Bool_t Run();

  static const char *MethodLabel() {
    return "bootstrap: mean trace inside the hand-drawn cut, resampled with "
           "the measured per-strip beam widths";
  }
 
private:
  struct Strip {
    TCutG *cut = nullptr;       // the channel's hand-drawn cut at this strip
    std::vector<Double_t> mean; // mean trace of the events inside it
    Double_t count = 0.0;       // how many
  };
  struct Outcome {
    Long64_t n = 0, tagged = 0, inside = 0, next_inside = 0;
  };
 
  Bool_t LoadSigmas(TFile &cache);
  // One channel end to end; kTRUE when it produced at least one record.
  Bool_t RunChannel(const CrossSectionChannel &ch, TTree *reservoir);
  void LoadDrawnCuts(const TString &region);
  void MeanTraces(TTree *reservoir);
  Outcome Bootstrap(Int_t reac, TH2F *plane);
  void DrawMeanTrace(const TString &subdir, Int_t reac, const Strip &st);
  void DrawPlane(const TString &subdir, Int_t reac, TH2F *plane, TCutG *cut);
 
  std::map<Int_t, Strip> strips_; // the channel being run
};
 
#endif