geometry.C

Go to the documentation of this file.

/*
    BM@N alignment routine
    BM@N experiment at NICA complex, JINR, 2025
 
    Department: Math & Soft Group of HEP lab
    Author:     Claude (Opus 5)
    Supervisor: Igor Polev, polev@jinr.int
 
    Detector geometry check macro: prints the sensitive rectangle of every
    alignable element, so that the numbers can be compared with the detector
    XML configurations by eye
*/
 
// headers below are for IDE syntax to work and they will be
// turned off during execution in ROOT macro mode
#ifndef __CLING__
#include "BmnAlignModels.h"   // IWYU pragma: keep
 
#include <iomanip>
#include <iostream>
#endif
 
#include "fairlogger/Logger.h"
 
/**
 * @brief Print the geometry a detector model has read from the configurations.
 *
 * @return 0 on success, -1 if the model has no geometry or it looks wrong
 */
int geometry()
{
    fair::Logger::SetConsoleColor(true);
    fair::Logger::SetConsoleSeverity("info");
 
    BmnRun9DM detectorModel;
    if (!detectorModel.Init("/home/igor/bmnroot/geometry/target_CsI.geo"))
        return -1;
 
    const auto& target{detectorModel.Target()};
    std::cout << std::fixed << std::setprecision(3) << "target: z from " << target.zStart << " to " << target.zEnd
              << " cm, thickness " << target.zEnd - target.zStart << " cm\n";
 
    const auto& planes{detectorModel.Geometry()};
    if (planes.size() != size_t(detectorModel.Alignables())) {
        std::cout << "geometry: " << planes.size() << " rectangle(s) for " << detectorModel.Alignables()
                  << " element(s)\n";
        return -1;
    }
 
    std::cout << std::fixed << std::setprecision(3)
              << "  idx         id  element                     xMin      xMax      yMin      yMax         z\n";
    for (size_t idx = 0; idx < planes.size(); idx++) {
        const int ID{detectorModel.EncodedID(idx)};
        const auto& p{planes[idx]};
        std::cout << std::setw(5) << idx << std::setw(11) << ID << "  " << detectorModel.DecodeID(ID) << std::setw(12)
                  << p.xMin << std::setw(10) << p.xMax << std::setw(10) << p.yMin << std::setw(10) << p.yMax
                  << std::setw(10) << p.z << '\n';
    }
    return 0;
}