Question

一个前提，我不是程序员，我是物理学家，我使用c ++作为分析数据的工具（ROOT包）。我的知识可能有限！

我有这种情况，我从文件中读取数据并将它们存储在一个向量中（没问题）

vector<double> data;

有了这些数据，我想绘制一个相关图，所以我需要将它们分成两个不同的子集，其中一个是2D直方图的X条目，另一个是Y条目。

分裂必须如下，我有这个表（我只复制它的一小部分只是为了解释问题）

************* LBA - LBC **************
--------------------------------------
Cell Name  |  Channel   |    PMT     |
D0         |          0 |          1 |
A1-L       |          1 |          2 |
BC1-R      |          2 |          3 |
BC1-L      |          3 |          4 |
A1-R       |          4 |          5 |
A2-L       |          5 |          6 |
BC2-R      |          6 |          7 |
BC2-L      |          7 |          8 |
A2-R       |          8 |          9 |
A3-L       |          9 |         10 |
A3-R       |         10 |         11 |
BC3-L      |         11 |         12 |
BC3-R      |         12 |         13 |
D1-L       |         13 |         14 |
D1-R       |         14 |         15 |
A4-L       |         15 |         16 |
BC4-R      |         16 |         17 |
BC4-L      |         17 |         18 |
A4-R       |         18 |         19 |
A5-L       |         19 |         20 |
...
None       |         31 |         32 |

正如您所看到的，A1-L和A1-R等条目对应于一个单元格的左侧和右侧，此左侧和右侧相关联int对应于频道，在本例中为1和4.我希望这些左侧和右侧位于我的2D直方图的X轴和Y轴上。

然后问题是以某种方式将这个表与数据向量相关联，这样我就可以选择属于同一个单元格的通道，并将它们放在X轴上，另一个放在Y轴上。为了使事情复杂化，还有一些没有合作伙伴的单元格，如本示例D0和没有像通道31那样关联的单元格的通道。

我尝试的解决方案是创建索引向量

vector<int> indexing = (0, 1, 4, ....);

和有序数据向量

vector<double> data_ordered;

并用

之类的东西填充有序矢量

for( vector<int> iterator it = indexing.begin(); it != indexing.end(); ++it)
    data_ordered.push_back(data.at(*it));

然后将data_ordered的偶数索引放在X轴上，奇数值放在Y轴上，但我遇到D0单元格和空单元格的问题！

我的另一个想法是创建一个struct之类的

struct cell{
    string cell_name;
    int left_channel;
    int right_channel;
    double data;
    ....
    other informations
};

然后尝试使用它，但在那里我缺乏c ++知识！有人能给我一个如何解决这个问题的提示吗？我希望我的问题足够清楚，并且尊重本网站的规则！

EDIT ----------

为了澄清问题，我尝试用一个例子来解释它

vector<double> data = (data0, data1, data2, data3, data4, ...);

do data0有索引0，如果我去表格，我看到它对应于没有其他伙伴的cell D0，让我们说现在可以忽略不计。 data1有索引1，它对应于cell A1（A1-L）的左侧部分，所以我需要在表格中找到合适的合作伙伴，其中有索引4，理想情况下我会选择来自包含数据的向量的data4。我希望这至少可以澄清一下情况！

Answer 1

这是一个可以做到你想要的引擎，粗略地说：

#include <vector>
#include <map>
#include <string>
#include <iostream>

enum sub_entry { left, right, only };

struct DataType {
  std::string cell;
  sub_entry sub;
  DataType( DataType const& o ): cell(o.cell), sub(o.sub) {};
  DataType( const char* c, sub_entry s=only ):
    cell( c ),
    sub( s )
  {}
  DataType(): cell("UNUSED"), sub(only) {};
  // lexographic weak ordering:
  bool operator<( DataType const& o ) const {
    if (cell != o.cell)
      return cell < o.cell;
    return sub < o.sub;
  }
};

typedef std::vector< double > RawData;
typedef std::vector< DataType > LookupTable;
typedef std::map< DataType, double > OrganizedData;

OrganizedData organize( RawData const& raw, LookupTable const& table )
{
  OrganizedData retval;
  for( unsigned i = 0; i < raw.size() && i < table.size(); ++i ) {
    DataType d = table[i];
    retval[d] = raw[i];
  }
  return retval;
}

void PrintOrganizedData( OrganizedData const& data ) {
  for (OrganizedData::const_iterator it = data.begin(); it != data.end(); ++it ) {
    std::cout << (*it).first.cell;
    switch( (*it).first.sub ) {
      case left: {
        std::cout << "-L";
      } break;
      case right: {
        std::cout << "-R";
      } break;
      case only: {
      } break;
    }
    std::cout << " is " << (*it).second << "\n";
  }
}

int main() {
  RawData test;
  test.push_back(3.14);
  test.push_back(2.8);
  test.push_back(-1);
  LookupTable table;
  table.resize(3);
  table[0] = DataType("A1", left);
  table[1] = "D0";
  table[2] = DataType("A1", right);
  OrganizedData org = organize( test, table );
  PrintOrganizedData( org );
}

查找表存储哪些通道映射到哪个单元名称和侧面。

查找表中未使用的条目应设置为DataType()，这会将其值标记为存储在"UNUSED"位置。（它仍将被存储，但你可以在之后丢弃它。）

结果是从(CellName, Side)到double数据的地图。我包括一个简单的打印机，只是转储数据。如果你有图形软件，你可以找到一种方法来制作图表。跳过"UNUSED"是一项涉及在该打印循环中检查(*it).first.cell == "UNUSED"的练习。

我相信一切都符合C ++ 03标准。如果你有一个C ++ 11编译器，那么上面的一堆会变得更漂亮。

根据表分裂矢量

1 个答案: