---

TMuiHVMask.cxx

Go to the documentation of this file.

#include<TMuiHVMask.h>
#include<TMuiKeyGen.h>
#include<fstream>
#include<sstream>
#include<MUIOO.h>
#include<boost/array.hpp>
#include<string>
#include<MutGeom.h>
#include<MuiGeomClasses.hh>

ClassImp(TMuiHVMask)

// Static data members definitions
//
sgi::hash_map<ULong_t, double> TMuiHVMask::_map;
ULong_t TMuiHVMask::_run_number = 0;
bool TMuiHVMask::_use_override = false;
float TMuiHVMask::_one_tube_effic = 0.65;
float TMuiHVMask::_two_tube_effic = 0.93;
bool TMuiHVMask::_disabled = true;
int TMuiHVMask::_mask_plane = -1;
std::string TMuiHVMask::_filename_north("tube_eff_north_default.txt");
std::string TMuiHVMask::_filename_south("tube_eff_south_default.txt");                  


bool
TMuiHVMask::initialize()
{  
  // Default filenames
  //
  boost::array<const char*,2> files = {{_filename_north.c_str(),_filename_south.c_str()}};
  MUTOO::PRINT(std::cout,"TMuiHVMask::Initialize");
  std::cout << _filename_north.c_str() << std::endl;
  std::cout << _filename_south.c_str() << std::endl;
  MUTOO::PRINT(std::cout,"**");
  
  for(size_t i=0; i<files.size(); ++i){

    ifstream lookup_file;    

    // Input file
    //
    lookup_file.open(files[i]);
    if(!lookup_file.is_open()) {
      MUIOO::TRACE("TMuiHVMask: can't find input file");
      continue;
    }
    
    // Fill the map
    //
    while(lookup_file) {

      UShort_t iarm=0, iplane=0, ipanel=0, iorient=0, itwopack=0;
      float effic=0;
      lookup_file >> iarm;
      lookup_file >> iplane;
      lookup_file >> ipanel;
      lookup_file >> iorient;
      lookup_file >> itwopack;
      lookup_file >> effic;    

      // Figure out how to to this correctly
      //
      if(iarm==0&&iplane==0&&ipanel==0&&iorient==0&&itwopack==0&&effic==0.0) continue;

      ULong_t key = TMuiKeyGen::get_key(iarm,iplane,ipanel,iorient,itwopack);
      
      // The number read is is 0, 65, 90 -- We reserve the right to adjust these
      // numbers hence the condition below that selects either the one tube or
      // two tube efficiency numbers from class data members
      //
      if(_use_override) {
        if(effic == 90) {
          _map.insert(std::make_pair(key,_two_tube_effic));
        } else if(effic == 65){
          _map.insert(std::make_pair(key,_one_tube_effic));
        } else if(effic == 0){
          _map.insert(std::make_pair(key,0));
        }
      } else {
        _map.insert(std::make_pair(key,effic));
      }
    }
    lookup_file.close();
  }
  
  // Message for the masses
  //
  MUIOO::PRINT(std::cout, "TMuiHVMask Initialization");
  std::cout << "TMuiHVMask read " << _map.size() << " entries" << std::endl;
  MUIOO::PRINT(std::cout, "**");
  
  return true;
}


double
TMuiHVMask::get_effic_twopack(UShort_t arm, UShort_t plane, UShort_t panel, UShort_t orient, UShort_t twopack)
{
  if(_disabled) return 1.0;

  static bool init = initialize();
  if(!init) return 1.0;

  private_map::const_iterator iter = _map.find(TMuiKeyGen::get_key(arm,plane,panel,orient,twopack));
  if(iter == _map.end()) {    
    std::ostringstream msg;
    msg << "TMuiHVMask no entry " << arm << " " << panel << " " << orient << " " << twopack << " "  << " entries";
    MUIOO::TRACE(msg.str().c_str());
    return _two_tube_effic;    
  } else {
    return iter->second;
  }
}

double
TMuiHVMask::get_effic_hvgroup(UShort_t arm, UShort_t plane, UShort_t panel, UShort_t orient,  UShort_t hvgroup)
{
  TMuiPanelGeo* panel_geo = TMuiGeometry::Geom()->getPanel(arm,plane,panel);    
  EOrient_t muid_orient = (orient==0) ? kHORIZ : kVERT;
  UShort_t ntwopack = panel_geo->getTwoPackCount(muid_orient);

  // Loop over two packs and return the effic for the first hv group
  //
  for(int twopack = 0; twopack<ntwopack; ++twopack){
    TMuiChannelId channel_id(arm,plane,panel,muid_orient,twopack);
    if(channel_id.get_HV_chain_group() == hvgroup){
      return get_effic_twopack(arm,plane,panel,orient,twopack);
    }
  }
  return 0.0;
}

---