Constants.hpp

Go to the documentation of this file.

#ifndef CONSTANTS_HPP
#define CONSTANTS_HPP
 
#include "DedupStrategy.hpp"
#include "PulseHistoryGroups.hpp"
#include "RunEpoch.hpp"
#include "SlotLayout.hpp"
#include <Rtypes.h>
#include <RtypesCore.h>
#include <TString.h>
#include <iostream>
#include <map>
#include <utility>
#include <vector>


#ifndef MUSIC_HOT_PATH_LOGGING
#define MUSIC_HOT_PATH_LOGGING 0
#endif

const Int_t GATE_AXIS_GRID = -1;
 
struct StripSumScatterConfig {
  enum PureBeamGate {
    PURE_BEAM_GATE_S0_S1,
    PURE_BEAM_GATE_S1_S2,
    PURE_BEAM_GATE_S1_GRID
  };
  PureBeamGate PURE_BEAM_GATE;
  enum TagResolve { TAG_RESOLVE_FIRST_STRIP, TAG_RESOLVE_LARGEST_JUMP };
  TagResolve TAG_RESOLVE;
  Double_t PURE_BEAM_NSIGMA;

  Int_t POST_TRIGGER_SUM_STRIPS;
  Int_t POST_WINDOW_LAST_STRIP;

  Int_t MAX_STRIP_SUM_WORKERS;
 
  Int_t REACTION_STRIP_MIN;
  Int_t REACTION_STRIP_MAX;

  Double_t SMOOTHNESS_NSIGMA;

  Int_t TAIL_FALL_FROM_STRIP;
  Double_t TAIL_RISE_NSIGMA;

  Double_t TAIL_RETURN_NSIGMA;
  Double_t TAIL_RERISE_NSIGMA;

  Double_t POST_ABOVE_NSIGMA;
  Int_t POST_ABOVE_STRIPS;

  Int_t POST_CROSS_MIN_STRIP;

  Double_t TAIL_CLIFF_MAX_FRACTION;

  Bool_t CUT_VARIATION;
  Double_t CUT_VARIATION_NSIGMA_STEP;
  Double_t CUT_VARIATION_CLIFF_STEP;

  Bool_t Y_RATIO_TO_UPSTREAM;

  Double_t REAC_JUMP_NSIGMA;

  Double_t END_STRIP_NSIGMA;

  Double_t PILEUP_NSIGMA;
  Int_t PILEUP_MIN_STRIPS;
  Double_t NOISE_NSIGMA;
  Int_t NOISE_MIN_STRIPS;

  Bool_t REGION_CUT_REDRAW;

  Double_t AN_REGION_NSIGMA;
  Double_t AA_REGION_NSIGMA;

  enum AnRegionMode { AN_REGION_MIXTURE, AN_REGION_ALL_TAGGED };
  AnRegionMode AN_REGION_MODE;

  Bool_t REQUIRE_BEAM_UPSTREAM_OF_REAC;
  Double_t BEAM_UPSTREAM_NSIGMA;

  Int_t BOTH_MULT_MAX;
  Int_t BOTH_MULT_COUNT_TO;
 
  Int_t SEED_HALF_BINS;
 
  Int_t TRACES_PER_CLASS;
 
  Int_t X_LO;
  Int_t X_HI;

  Int_t GATE_STRIP;
  Double_t GATE_NSIGMA;
  Double_t GATE_CENTER;
  Double_t GATE_MIN;
  Double_t GATE_MAX;
  Int_t GATE_BINS;

  Double_t X_DISPLAY_MIN;
  Double_t X_DISPLAY_MAX;
  Int_t XBINS;
  Double_t Y_DISPLAY_MIN;
  Double_t Y_DISPLAY_MAX;
  Int_t YBINS;

  std::map<Int_t, std::pair<Double_t, Double_t>> Y_DISPLAY_RANGE;
 
  Long64_t SAMPLE_MAX_POINTS;
 
  Bool_t RERUN_SIM;
  Int_t CANDIDATE_REAC_STRIP;

  Bool_t SKIP_SAVGOL_PLOTS;

  Bool_t PLOT_REGION_MEAN_TRACES;
  Bool_t PLOT_ADC_TRACES;
  Bool_t REQUIRE_STRIP_16_BELOW_BEAM;
 
  void SetDefaults();
};


struct CrossSectionChannel {
  TString name;
  TString label;
  std::vector<TString> talys_exits;
  std::vector<TString> subtract_exits;
  std::vector<std::vector<Double_t>> reference_xs;
  TString reference_label;
  Int_t fit_model = -1;
};

struct TalysModel {
  TString label;                 
  std::vector<TString> keywords; 
};

enum TargetGas { kHELIUM  };

Int_t TargetGasA(TargetGas gas);

Double_t TargetGasAtomsPerMolecule(TargetGas gas);

struct CrossSectionConfig {
  TargetGas TARGET_GAS;
  Double_t GAS_PRESSURE_TORR;
  Double_t GAS_PRESSURE_TORR_ERR;

  Int_t BEAM_A;
  Int_t BEAM_Z;
  TString BEAM_ELEMENT;

  std::vector<TalysModel> TALYS_MODELS;
  std::vector<TString> TALYS_COMMON_KEYWORDS;

  TString BEAM_SIM_FILE;

  Int_t XS_STRIP_MIN;
  Int_t XS_STRIP_MAX;

  Int_t FELDMAN_COUSINS_MAX_COUNT;

  Bool_t EFFECTIVE_ENERGY;

  Bool_t PRELIMINARY;

  std::vector<TString> EPOCHS;

  std::vector<CrossSectionChannel> CHANNELS;
 
  void SetDefaults();
};
 
class DatasetConfig {
public:
  Bool_t USE_SOLARIS_DATA;
  TString SOL_BASE_DIR;
  TString SOL_SPLIT_DIR;
  Double_t SOL_SPLIT_CHUNK_SECONDS;
  Int_t SOL_N_SPLIT_WORKERS;
 
  TString COMPASS_BASE_DIR;

  std::vector<Int_t> RUN_NUMBERS;
  std::vector<RunEpoch> EPOCHS;
  Int_t N_CHUNKS;
 
  TString SIM_BEAM_FILE;

  Int_t N_BOARDS;
  Int_t N_CHANNELS;
  UShort_t TIMING_REF_BOARD;
  std::vector<UShort_t> TIMING_REF_BOARD_CHANNELS;

  Double_t TIMING_MIN_ENERGY;
  Double_t TIMING_MAX_ENERGY;
  Double_t TIMING_OVERLAP_MARGIN_S;
  Double_t TIMING_THRESH_DT_US;
  Double_t TIMING_MAX_ABS_SHIFT_S;
  Double_t TIMING_SHIFT_COARSE_STEP_US;
  Double_t TIMING_SHIFT_FINE_STEP_US;
  Double_t TIMING_SHIFT_FINE_HALF_WIDTH_US;
  Int_t TIMING_SHIFT_MIN_NPTS;
  Int_t TIMING_SHIFT_MAX_SCAN_CANDIDATES;
  Bool_t TIMING_DO_BOARD_SYNC;
  Bool_t TIMING_DO_SORT;
 
  Bool_t REJECT_FLAGGED_EVENTS;

  Bool_t PULSE_HISTORY_CORRECTION;
  Long64_t PULSE_HISTORY_MIN_EVENTS;
  Double_t PULSE_HISTORY_BEAM_LO;
  Double_t PULSE_HISTORY_BEAM_HI;

  Double_t PULSE_HISTORY_OWN_LO;
  Double_t PULSE_HISTORY_OWN_HI;
 
  // How far back to check and apply a correction. The trapezoid baseline is an
  // average over a
  Double_t PULSE_HISTORY_APPLY_MAX_US;

  Int_t PULSE_HISTORY_AMP_BINS;

  PulseHistoryGroups PULSE_HISTORY_GROUPS;

  Double_t PULSE_HISTORY_TRAP_RISE_US;
  Double_t PULSE_HISTORY_TRAP_FLAT_US;
  Double_t PULSE_HISTORY_PEAKING_US;

  Bool_t IGNORE_SHORT_STRIPS;
  Bool_t IGNORE_STRIP_0;
  Bool_t IGNORE_STRIP_17;
  Bool_t REQUIRE_STRIP_0;
 
  Bool_t HAS_CATHODE;
  Bool_t HAS_GRID;
  Bool_t HAS_STRIP0;
  Bool_t HAS_STRIP17;
 
  Bool_t SKIP_EXISTING;

  Bool_t SAVE_FULL_PLOTS;
  Int_t PLOT_SAMPLE_FILES;

  Int_t SAVE_SAMPLE_TRACES;
 
  Int_t MAX_FUSED_WORKERS;
 
  TString REFERENCE_CHANNEL;
  Double_t EVENT_TIME_WINDOW_US;
  Double_t SEED_HOLDOFF_US;
  Double_t SEED_HOLDOFF_MAX_RATIO;
  DedupStrategy DEDUP_STRATEGY;
 
  Bool_t USE_GPU_ACCELERATION;
  Int_t MAX_GPU_CONCURRENT_SORTS;

  Double_t BEAM_GATE_NSIGMA;
 
  Double_t STRIP_DE_MIN_NORMED;
  Double_t STRIP_DE_MAX_NORMED;
 
  StripSumScatterConfig STRIP_SUM_SCATTER_CONFIG;
  CrossSectionConfig CROSS_SECTION_CONFIG;
 
  Double_t STRIP_E_MIN_ADC;
  Double_t STRIP_E_MAX_ADC;
 
  Double_t GRID_MAX_ADC;
 
  Double_t REFERENCE_CHANNEL_MIN_ADC;
  Double_t REFERENCE_CHANNEL_MAX_ADC;
 
  Double_t STRIP0_MAX_ADC;
  Double_t STRIP17_MAX_ADC;
  Double_t CATHODE_MAX_ADC;
 
  Double_t LEFT_EVEN_MAX_ADC;
  Double_t LEFT_ODD_MAX_ADC;
  Double_t RIGHT_EVEN_MAX_ADC;
  Double_t RIGHT_ODD_MAX_ADC;
 
  std::map<std::pair<Int_t, Int_t>, TString> channelMap;
  std::map<std::pair<Int_t, Int_t>, TString> channelMap64;
 
  DatasetConfig();
};
 
namespace Constants {
extern const DatasetConfig &cfg;

void SetActiveEpoch(const RunEpoch *epoch);

const RunEpoch *GetActiveEpoch();
const RunEpoch *EpochForRun(Int_t run);
RunEpoch MakeEpoch(const TString &name, const std::vector<Int_t> &runs);
std::vector<Int_t> RunRange(Int_t first, Int_t last);
const TString &ActiveFileTag();

Int_t ActiveNBoards();
Int_t ActiveNChannels();
UShort_t ActiveTimingRefBoard();
const std::vector<UShort_t> &ActiveTimingRefBoardChannels();
Bool_t ActiveDoBoardSync();
Bool_t ActiveDoSort();
Bool_t ActiveHasCathode();
Bool_t ActiveUseSolarisData();
Double_t ActiveEventTimeWindowUs();
Double_t ActiveSeedHoldoffUs();
Double_t ActiveSeedHoldoffMaxRatio();
const TString &ActiveReferenceChannel();
Double_t ActiveReferenceChannelMinAdc();
Double_t ActiveReferenceChannelMaxAdc();
DedupStrategy ActiveDedupStrategy();
Double_t ActiveStripEMinAdc();
Double_t ActiveStripEMaxAdc();
Double_t ActiveCathodeMaxAdc();
Double_t ActiveGridMaxAdc();
Double_t ActiveStrip0MaxAdc();
Double_t ActiveStrip17MaxAdc();
Double_t ActiveLeftEvenMaxAdc();
Double_t ActiveLeftOddMaxAdc();
Double_t ActiveRightEvenMaxAdc();
Double_t ActiveRightOddMaxAdc();

const std::vector<Int_t> &ActiveRunNumbers();

Int_t ActiveMaxFiles();

Double_t ActiveSplitChunkSeconds();
const PulseHistoryGroups &ActivePulseHistoryGroups();
Bool_t ActiveIgnoreShortStrips();
void ActivateAnalysisEpoch();

Bool_t FileInSample();
Bool_t SavePlots();
std::ostream &Detail();
std::ostream &DetailErr();
Bool_t InPlotSample(Int_t k);
void SetPlotsThisFile(Bool_t on);

const std::map<std::pair<Int_t, Int_t>, TString> &ActiveChannelMap();
} // namespace Constants
 
#endif