ARICHAeroTilesInfo.cc

/**************************************************************************
 * basf2 (Belle II Analysis Software Framework)                           *
 * Author: The Belle II Collaboration                                     *
 *                                                                        *
 * See git log for contributors and copyright holders.                    *
 * This file is licensed under LGPL-3.0, see LICENSE.md.                  *
 **************************************************************************/
 
#include <framework/logging/Logger.h>
 
#include <arich/dbobjects/ARICHAeroTilesInfo.h>
#include <iostream>
 
using namespace std;
using namespace Belle2;
 
float ARICHAeroTilesInfo::getRefractiveIndex(int slot, int layer) const
{
 
  if (slot > 124 || slot < 1) { B2WARNING("ARICHAeroTilesInfo::getTransmissionLength: Slot number " << slot << " not valid!"); return 0.;}
  if (layer > 1 || layer < 0) { B2WARNING("ARICHAeroTilesInfo::getTransmissionLength: Layer number " << layer << " not valid!"); return 0.;}
 
  std::vector<int> slotID{slot, layer};
  std::vector<float> opticalProperties = m_aeroPos2optical.find(slotID)->second;
 
  return opticalProperties[0];
}
 
float ARICHAeroTilesInfo::getTransmissionLength(int slot, int layer) const
{
 
  if (slot > 124 || slot < 1) { B2WARNING("ARICHAeroTilesInfo::getTransmissionLength: Slot number " << slot << " not valid!"); return 0.;}
  if (layer > 1 || layer < 0) { B2WARNING("ARICHAeroTilesInfo::getTransmissionLength: Layer number " << layer << " not valid!"); return 0.;}
 
  std::vector<int> slotID{slot, layer};
  std::vector<float> opticalProperties = m_aeroPos2optical.find(slotID)->second;
 
  return opticalProperties[1];
}
 
float ARICHAeroTilesInfo::getThickness(int slot, int layer) const
{
 
  if (slot > 124 || slot < 1) { B2WARNING("ARICHAeroTilesInfo::getThickness: Slot number " << slot << " not valid!"); return 0.;}
  if (layer > 1 || layer < 0) { B2WARNING("ARICHAeroTilesInfo::getThickness: Layer number " << layer << " not valid!"); return 0.;}
 
  std::vector<int> slotID{slot, layer};
  std::vector<float> opticalProperties = m_aeroPos2optical.find(slotID)->second;
 
  return opticalProperties[2];
}
 
void ARICHAeroTilesInfo::addMapping(int slot, int layer, float refractiveIndex, float transmissionLength, float thickness)
{
 
  if (slot > 124 || slot < 1) { B2WARNING("ARICHAeroTilesInfo::getTransmissionLength: Slot number " << slot << " not valid!"); }
  if (layer > 1 || layer < 0) { B2WARNING("ARICHAeroTilesInfo::getTransmissionLength: Layer number " << layer << " not valid!"); }
 
  std::vector<int> slotID{slot, layer};
  std::vector<float> opticalProperties{refractiveIndex, transmissionLength, thickness};
  m_aeroPos2optical.insert(std::pair<std::vector<int>, std::vector<float>>(slotID, opticalProperties));
}
 
void ARICHAeroTilesInfo::print() const
{
  B2INFO("AEROGEL SLOT MAPPING with optical properties of tiles");
  for (int tile = 1; tile < N_TILES + 1; tile++) {
    B2INFO("Slot number: " << tile << "\n"
           << "       downstream tile : refIndex    = " << getRefractiveIndex(tile, 0) << "\n"
           << "                         transLength = " << getTransmissionLength(tile, 0) << "\n"
           << "                         thickness   = " << getThickness(tile, 0) << "\n"
           << "       upstream tile   : refIndex    = " << getRefractiveIndex(tile, 1) << "\n"
           << "                         transLength = " << getTransmissionLength(tile, 1) << "\n"
           << "                         thickness   = " << getThickness(tile, 1) << "\n");
    //B2INFO("Slot number: " << tile << "\n"
    //     << "       downstream tile : refIndex    = " << getRefractiveIndex(tile, 0) << "\n"
    //     << "                         transLength = " << getTransmissionLength(tile, 0) << "\n"
    //     << "       upstream tile   : refIndex    = " << getRefractiveIndex(tile, 1) << "\n"
    //     << "                         transLength = " << getTransmissionLength(tile, 1) << "\n");
  }
}