sPhenix code displayed by LXR master/​coresoftware/​offline/​packages/​tpc/​TpcMap.cc	Source navigation Diff markup Identifier search General search
	Version: master Revert   Change

File indexing completed on 2025-12-16 09:20:24

 #include "TpcMap.h"
 
 #include <cmath>
 #include <cstdlib>
 #include <fstream>
 #include <iostream>
 #include <limits>
 #include <sstream>
 #include <utility>
 
 void TpcMap::setMapNames(const std::string &r1, const std::string &r2, const std::string &r3)
 {
   const char *calibrationroot = getenv("CALIBRATIONROOT");
   std::string full_path;
   if (calibrationroot)
   {
     full_path = std::string(calibrationroot) + "/TPC/Mapping/PadPlane/";
   }
   else
   {
     calibrationroot = getenv("TPCCALIB");
     if (calibrationroot)
     {
       full_path = std::string(calibrationroot) + "/";
     }
     else
     {
       full_path = "./";
     }
   }
 
   std::string full_path_r1 = full_path + r1;
   std::string full_path_r2 = full_path + r2;
   std::string full_path_r3 = full_path + r3;
   int status;
   status = digest_map(full_path_r1, 0);
   if (status)
   {
     std::cout << "reading " << full_path_r1 << " failed" << std::endl;
   }
   status = digest_map(full_path_r2, 1);
   if (status)
   {
     std::cout << "reading " << full_path_r2 << " failed" << std::endl;
   }
   status = digest_map(full_path_r3, 2);
   if (status)
   {
     std::cout << "reading " << full_path_r3 << " failed" << std::endl;
   }
 }
 
 int TpcMap::digest_map(const std::string &fileName, const unsigned int section_offset)
 {
   std::ifstream infile(fileName, std::ios::in);
 
   if (!infile.is_open())
   {
     std::cout << "Could not open file: " << fileName << std::endl;
     _broken = 1;
     return -1;
   }
 
   std::string line;
   getline(infile, line);  // throwaway - we skip the first line
   //  cout << __FILE__<< " " << __LINE__ << ": " << line << endl;
 
   int abs_pad = std::numeric_limits<int>::max();
   int Radius = std::numeric_limits<int>::max();
   // int Pad;
   // int  U;
   // int  G;
   // string  Pin;
   // int PinColID;
   // int PinRowID;
   // string PadName;
   int FEE = std::numeric_limits<int>::max();
   //  int FEE_Connector;
   int FEE_Chan = std::numeric_limits<int>::max();
   // double phi;
   // double x;
   // double y;
   // double PadX;
   // double PadY;
   double PadR = NAN;
   double PadPhi = NAN;
 
   while (getline(infile, line))
   {
     //      cout << line<< endl;
     std::stringstream ss(line);
     std::string next;
 
     //  0  26    26
     //  1   0    0
     //  2   0    0
     //  3   1    1
     //  4   1    1
     //  5   0    C5
     //  6   2    2
     //  7   5    5
     //  8   0    ZZ.00.000
     //  9   5    5.0
     // 10   0    J2
     // 11 147    147
     // 12   0    0.005570199740407434
     // 13  69    69.99891405342764
     // 14   0    0.38991196551332985
     // 15  77    77.86043476908294
     // 16 305    305.14820499531316
     // 17 314    314.9248709046211
     // 18   0    0.24982557215053805
 
     int index = 0;
     while (ss.good())
     {
       getline(ss, next, ',');
       if (index == 2)
       {
         abs_pad = std::stoul(next);
       }
       else if (index == 1)
       {
         Radius = stoul(next);
       }
       else if (index == 9)
       {
         FEE = stoul(next);
       }
       else if (index == 11)
       {
         FEE_Chan = stoul(next);
       }
       else if (index == 17)
       {
         PadR = stod(next);
       }
       else if (index == 18)
       {
         PadPhi = stod(next);
       }
       index++;
     }
 
     if (section_offset == 1)
     {
       FEE += 6;
     }
     if (section_offset == 2)
     {
       FEE += 14;
     }
 
     struct tpc_map x
     {
     };
     x.padnr = abs_pad;
     x.layer = Radius + section_offset * 16 + 7;
     x.FEE = FEE;
     x.FEEChannel = FEE_Chan;
     x.PadR = PadR;
     x.PadPhi = PadPhi;
 
     unsigned int key = 256 * (FEE) + FEE_Chan;
     tmap[key] = x;
     //if(Radius<0){
     //  std::cout << " " << key << " " << FEE << " " << FEE_Chan << " " << PadR << "  " << PadPhi << " " << Radius << std::endl;
     //}
   }
   return 0;
 }
 
 unsigned int TpcMap::getLayer(const unsigned int FEE, const unsigned int FEEChannel, const unsigned int /* packetid */) const
 {
   if (FEE >= 26 || FEEChannel > 255)
   {
     return 0.;
   }
   unsigned int key = 256 * FEE + FEEChannel;
 
   std::map<unsigned int, struct tpc_map>::const_iterator itr = tmap.find(key);
   if (itr == tmap.end())
   {
     return 0;
   }
   return itr->second.layer;
 }
 
 unsigned int TpcMap::getPad(const unsigned int FEE, const unsigned int FEEChannel, const unsigned int /* packetid */) const
 {
   if (FEE >= 26 || FEEChannel > 255)
   {
     return 0.;
   }
   unsigned int key = 256 * FEE + FEEChannel;
 
   std::map<unsigned int, struct tpc_map>::const_iterator itr = tmap.find(key);
   if (itr == tmap.end())
   {
     return -100;
   }
   return itr->second.padnr;
 }
 
 double TpcMap::getR(const unsigned int FEE, const unsigned int FEEChannel, const unsigned int /* packetid */) const
 {
   if (FEE >= 26 || FEEChannel > 255)
   {
     return 0.;
   }
   unsigned int key = 256 * FEE + FEEChannel;
 
   std::map<unsigned int, struct tpc_map>::const_iterator itr = tmap.find(key);
   if (itr == tmap.end())
   {
     return -100;
   }
   return itr->second.PadR;
 }
 
 double TpcMap::getPhi(const unsigned int FEE, const unsigned int FEEChannel, const unsigned int /* packetid */) const
 {
   if (FEE > 25 || FEEChannel > 255)
   {
     return 0.;
   }
   unsigned int key = 256 * FEE + FEEChannel;
 
   std::map<unsigned int, struct tpc_map>::const_iterator itr = tmap.find(key);
   if (itr == tmap.end())
   {
     return -100;
   }
   return itr->second.PadPhi;
 }

This page was automatically generated by the 2.3.7 LXR engine. The LXR team