Constants.hpp

Go to the documentation of this file.

#ifndef CONSTANTS_HPP
#define CONSTANTS_HPP
 
#include "DedupStrategy.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 1
#endif

struct StripSumScatterConfig {
  enum PureBeamGate { PURE_BEAM_GATE_S0_S1, PURE_BEAM_GATE_S1_S2 };
  PureBeamGate PURE_BEAM_GATE;

  Int_t POST_TRIGGER_SUM_STRIPS;
  Int_t POST_WINDOW_LAST_STRIP;
  std::map<Int_t, Int_t> POST_WINDOW_STRIPS;
  Int_t MAX_STRIP_SUM_WORKERS;
 
  Int_t REACTION_STRIP_MIN;
  Int_t REACTION_STRIP_MAX;
 
  Int_t REQUIRE_SMOOTHNESS_END_STRIP;
  Double_t REQUIRE_SMOOTHNESS_MAX_STEP;

  Double_t REAC_JUMP_NSIGMA;
  Double_t REAC_JUMP_MAX;
  Double_t END_STRIP_MAX;
 
  Double_t PILEUP_THRESHOLD;
  Double_t NOISE_THRESHOLD;
 
  Bool_t REJECT_NOISE;
  Double_t NOISE_THRESH_PY;
  Int_t NOISE_MIN_STRIPS;
 
  Bool_t REJECT_PILEUP;
  Double_t PILEUP_THRESH_PY;
  Int_t PILEUP_MIN_STRIPS;

  Bool_t REGION_CUT_REDRAW;

  Double_t AN_REGION_NSIGMA;
  Double_t AA_REGION_NSIGMA;
  enum AnRegionMode { AN_REGION_MIXTURE, AN_REGION_RIDGE_BAND };
  AnRegionMode AN_REGION_MODE;
  Double_t AN_RIDGE_NSIGMA_LO;
  Double_t AN_RIDGE_NSIGMA_HI;

  Bool_t REQUIRE_BEAM_UPSTREAM_OF_REAC;
  Double_t BEAM_UPSTREAM_NSIGMA;

  Int_t BOTH_MULT_MAX;
  Int_t BOTH_MULT_COUNT_TO;
 
  Bool_t REJECT_OFFBEAM;
  Double_t OFFBEAM_DIST;
  Int_t OFFBEAM_MIN_STRIPS;
  Double_t PARITY_ASYM_MAX;

  Bool_t PLOT_PARITY_REJECTED_GRID;
 
  Double_t TRIGGER_NSIGMA;
  Double_t TRIGGER_CFD_FRAC;
  Int_t PLATEAU_POST;
  Int_t CLUSTER_SMOOTH_WINDOW;
  Int_t SEED_HALF_BINS;
 
  Int_t SAVGOL_HALF;
 
  Int_t TRACES_PER_CLASS;
 
  Int_t X_LO;
  Int_t X_HI;
 
  Int_t GATE_STRIP_X;
  Int_t GATE_STRIP_Y;
  Double_t GATE_NSIGMA_X;
  Double_t GATE_NSIGMA_Y;
  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 REQUIRE_SMOOTHNESS;
  Bool_t REQUIRE_GATE_S3_S4;
  Bool_t REQUIRE_GATE_S5_S6;
  Bool_t SKIP_SAVGOL_PLOTS;
  Bool_t REQUIRE_STRIP_16_BELOW_BEAM;

  Bool_t ALT_DECODE_REGION_TRACES;
 
  void SetDefaults();
};
 
// One TALYS calculation: what the plot calls it, and the input lines that
// select it (see CrossSectionConfig::TALYS_MODELS).
// One measured reaction channel; see CrossSectionConfig::CHANNELS.
struct CrossSectionChannel {
  TString name;
  TString label;
  std::vector<TString> talys_exits;
  std::vector<std::vector<Double_t>> reference_xs;
  TString reference_label;
};

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;

  Int_t BEAM_A;
  Int_t BEAM_Z;
  TString BEAM_ELEMENT;

  std::vector<TalysModel> TALYS_MODELS;

  TString BEAM_SIM_FILE;

  Int_t XS_STRIP_MIN;
  Int_t XS_STRIP_MAX;

  Bool_t EFFECTIVE_ENERGY;

  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;
  Double_t PULSE_HISTORY_APPLY_MAX_US;
  Int_t PULSE_HISTORY_AMP_BINS;
 
  Bool_t IGNORE_SHORT_STRIPS;
  Bool_t IGNORE_STRIP_0;
  Bool_t IGNORE_STRIP_17;
 
  Bool_t HAS_CATHODE;
  Bool_t HAS_GRID;
  Bool_t HAS_STRIP0;
  Bool_t HAS_STRIP17;
 
  Bool_t SKIP_EXISTING;
  Bool_t SAVE_PLOTS;
 
  Bool_t SKIP_CALIBRATION;

  Int_t SAVE_SAMPLE_TRACES;
 
  Int_t MAX_FUSED_WORKERS;
 
  TString REFERENCE_CHANNEL;
  Double_t EVENT_TIME_WINDOW_US;
  DedupStrategy DEDUP_STRATEGY;
 
  Bool_t USE_GPU_ACCELERATION;
  Int_t MAX_GPU_CONCURRENT_SORTS;
 
  Double_t STRIP_DE_OVERVIEW_MIN_NORMED;
  Double_t STRIP_DE_OVERVIEW_MAX_NORMED;
  Double_t BEAM_GATE_NSIGMA_X;
  Double_t BEAM_GATE_NSIGMA_Y;
 
  Double_t STRIP_DE_MIN_NORMED;
  Double_t STRIP_DE_MAX_NORMED;
  Double_t CATHODE_E_MAX_NORMED;
  Double_t TOTAL_E_MIN_NORMED;
  Double_t TOTAL_E_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 TOTAL_E_MIN_ADC;
  Double_t TOTAL_E_MAX_ADC;
 
  Double_t GRID_MIN_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;
  std::map<std::pair<Int_t, Int_t>, Long64_t> ttfOffsetPs;
 
  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();
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();

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