InterCalibConfig.h

Go to the documentation of this file.

// Copyright 2019-2020 CERN and copyright holders of ALICE O2.
// See https://alice-o2.web.cern.ch/copyright for details of the copyright holders.
// All rights not expressly granted are reserved.
//
// This software is distributed under the terms of the GNU General Public
// License v3 (GPL Version 3), copied verbatim in the file "COPYING".
//
// In applying this license CERN does not waive the privileges and immunities
// granted to it by virtue of its status as an Intergovernmental Organization
// or submit itself to any jurisdiction.
 
#ifndef O2_ZDC_INTERCALIBCONFIG_H
#define O2_ZDC_INTERCALIBCONFIG_H
 
#include "ZDCBase/Constants.h"
#include <Rtypes.h>
#include <array>
#include <string>
#include <limits>
 
 
namespace o2
{
namespace zdc
{
struct InterCalibConfig {
  static constexpr int NH = 9; 
  bool enabled[NH] = {true, true, true, true, true, true, true, true, true};
  double cutLow[NH] = {-std::numeric_limits<float>::infinity(), -std::numeric_limits<float>::infinity(), -std::numeric_limits<float>::infinity(), -std::numeric_limits<float>::infinity(), -std::numeric_limits<float>::infinity(), -std::numeric_limits<float>::infinity(), -std::numeric_limits<float>::infinity(), -std::numeric_limits<float>::infinity(), -std::numeric_limits<float>::infinity()};
  double cutHigh[NH] = {std::numeric_limits<float>::infinity(), std::numeric_limits<float>::infinity(), std::numeric_limits<float>::infinity(), std::numeric_limits<float>::infinity(), std::numeric_limits<float>::infinity(), std::numeric_limits<float>::infinity(), std::numeric_limits<float>::infinity(), std::numeric_limits<float>::infinity(), std::numeric_limits<float>::infinity()};
 
  // Distance from calorimeter side close to the beam (always positive)
  // Meaningful values are in the range of tower x centers i.e. from
  // 2.8 to 19.6 If one puts less than 2.8 then the computation will be
  // the same as for ZPA/ZPC with no cuts
  double xcut_ZPA = 0;
  double xcut_ZPC = 0;
  double rms_cut_ZP = 0;                                                                                                                                                                         // RMS of ZP centroid can go from 0 to 8.4 cm
  double towerCutLow_ZPA[4] = {0, 0, 0, 0};                                                                                                                                                      // Applied to all ZP fits except ZPI
  double towerCutLow_ZPC[4] = {0, 0, 0, 0};                                                                                                                                                      // Applied to all ZP fits except ZPI
  double towerCutHigh_ZPA[4] = {std::numeric_limits<float>::infinity(), std::numeric_limits<float>::infinity(), std::numeric_limits<float>::infinity(), std::numeric_limits<float>::infinity()}; // Applied to all ZP fits except ZPI
  double towerCutHigh_ZPC[4] = {std::numeric_limits<float>::infinity(), std::numeric_limits<float>::infinity(), std::numeric_limits<float>::infinity(), std::numeric_limits<float>::infinity()}; // Applied to all ZP fits except ZPI
  bool cross_check = false;
 
  int nb1[NH] = {0};      
  double amin1[NH] = {0}; 
  double amax1[NH] = {0}; 
  int nb2[NH] = {0};      
  double amin2[NH] = {0}; 
  double amax2[NH] = {0}; 
  double start[NH] = {0.75, 0.75, 0.75, 0.75, 1., 1., 1., 0.75, 0.75};
  double l_bnd[NH] = {0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1};
  double u_bnd[NH] = {10., 10., 10., 10., 10., 10., 10., 10., 10.};
  double l_bnd_o[NH] = {-20., -20., -20., -20., -20., -20., -20., -20., -20.};
  double u_bnd_o[NH] = {20., 20., 20., 20., 20., 20., 20., 20., 20.};
  double step_o[NH] = {0., 0., 0., 0., 0., 0., 0., 0., 0.};
  double min_e[NH] = {0., 0., 0., 0., 0., 0., 0., 0., 0.};
  std::string desc = "";
 
  void print() const;
  void resetCuts();
  void resetCutLow();
  void resetCutHigh();
  void resetCutLow(int ih);
  void resetCutHigh(int ih);
  void setMinEntries(double val);
  void setMinEntries(int ih, double val);
  void setCutLow(double val);
  void setCutHigh(double val);
  void setCutLow(int ih, double val);
  void setCutHigh(int ih, double val);
  void setCuts(double low, double high);
  void setCuts(int ih, double low, double high);
  void setBinning1D(int nb, double amin, double amax);
  void setBinning2D(int nb, double amin, double amax);
  void setBinning1D(int ih, int nb, double amin, double amax);
  void setBinning2D(int ih, int nb, double amin, double amax);
  void setDescription(std::string d) { desc = d; }
  void enable(bool c0, bool c1, bool c2, bool c3, bool c4, bool c5, bool c6, bool c7, bool c8)
  {
    enabled[0] = c0;
    enabled[1] = c1;
    enabled[2] = c2;
    enabled[3] = c3;
    enabled[4] = c4;
    enabled[5] = c5;
    enabled[6] = c6;
    enabled[7] = c7;
    enabled[8] = c8;
  }
  ClassDefNV(InterCalibConfig, 5);
};
} // namespace zdc
} // namespace o2
 
#endif