CHeTGlobalSettings.hh

Go to the documentation of this file.

#ifndef CHETGLOBALSETTINGS_HH
#define CHETGLOBALSETTINGS_HH
 
#include <algorithm> // for std::find if needed
#include <cmath>
#include <iostream>
#include <map>
#include <string>
#include <vector>
 
namespace CHeT
{
 
namespace Config
{
 
// --- Configuration Structures ---
 
struct LayerConfig
{
    int nBundles; 
    double radius; 
    double thickness; 
    double phiOffset; 
    int direction; 
    double nominalStereoAngle; 
    int color; 
};
 
struct CylinderConfig
{
    int id; 
    double nominalRadius; 
    LayerConfig inner; 
    LayerConfig outer; 
};
 
struct FiberProp
{
    int cylinderId;
    int layerId; 
    double r; 
    double phi0; 
    int dir; 
    int color; 
};
 
struct BundlesIntersection
{
    double z; 
    double x; 
    double y; 
    double z_loc; 
    double x_loc; 
    double y_loc; 
    int cylinderId; 
    double widthZ; 
};
 
// --- Global Physical Constants ---
 
// Kept as constexpr in header for compile-time optimization
constexpr double L_HALF = 150.0; 
constexpr double FIBER_WIDTH = 0.5; 
constexpr int FIBERS_PER_SIPM = 4; 
constexpr double BUNDLE_WIDTH = FIBER_WIDTH * FIBERS_PER_SIPM; 
 
inline double GetStereoAngle(double radius, double deltaPhi = M_PI, double length = 2.0 * L_HALF)
{
    return std::atan((radius * deltaPhi) / length);
}
 
// --- Experimental Offsets ---
 
double GetOffsetExp();
 
void SetOffsetExp(double val);
 
double GetDeltaI(int cylIdx);
 
void SetDeltaI(int cylIdx, double val);
 
void SetDeltas(const std::vector<double> &deltas);
 
std::vector<double> GetDeltas();
 
// --- Rotation Settings ---
 
void SetRotation(double rx, double ry, double rz);
 
void GetRotation(double &rx, double &ry, double &rz);
 
void SetTranslation(double tx, double ty, double tz);
 
void GetTranslation(double &tx, double &ty, double &tz);
 
void ApplyRotation(double &x, double &y, double &z);
 
void ApplyInverseRotation(double &x, double &y, double &z);
 
void ApplyTransformation(double &x, double &y, double &z);
 
void ApplyInverseTransformation(double &x, double &y, double &z);
 
// --- Function Declarations ---
 
int GetBoardGlobalOffset(int board_id);
 
void SetActiveCylinders(const std::vector<int> &active_ids);
 
std::vector<int> GetActiveCylinders();
 
const std::vector<CylinderConfig> &GetCylinders();
 
double wrap0_2pi(double angle);
 
FiberProp GetFiberProp(int b_id);
 
int GetGlobalBundleId(int board_id, int channel_id);
 
std::vector<BundlesIntersection> FindIntersections(const std::vector<int> &hit_ids);
 
// --- Debug / Helper Functions ---
 
int GetGlobalIdFromGeometry(int cyl_id, int layer_id, int layer_idx);
 
void PrintBundleMapping(int global_id);
 
void MapExplorer();
 
} // namespace Config
 
} // namespace CHeT
 
#endif // CHETGLOBALSETTINGS_HH