forked from sPHENIX-Collaboration/coresoftware
/
PHTpcTrackerUtil.cc
66 lines (58 loc) · 2.11 KB
/
PHTpcTrackerUtil.cc
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
/*!
* \file PHTpcTrackerUtil.cc
* \brief
* \author Dmitry Arkhipkin <arkhipkin@gmail.com>
*/
#include "PHTpcTrackerUtil.h"
#include <trackbase/TrkrCluster.h>
#include <trackbase/TrkrClusterContainer.h> // for TrkrClusterContainer
#include <trackbase/TrkrHitSetContainer.h>
#include <trackbase/TrkrHitSet.h>
#include <trackbase/TrkrDefs.h>
#include <phool/PHLog.h>
#include <cstdint> // for uint64_t
#include <cstring>
#include <map> // for _Rb_tree_const_iterator
#include <utility> // for pair
namespace PHTpcTrackerUtil
{
std::vector<std::vector<double> > convert_clusters_to_hits(TrkrClusterContainer* cluster_map, TrkrHitSetContainer* hitsets )
{
//***** convert clusters to kdhits
std::vector<std::vector<double> > kdhits;
if (!cluster_map)
{
LOG_WARN("tracking.PHTpcTrackerUtil.convert_clusters_to_hits") << "cluster map is not provided";
return kdhits;
}
auto hitsetrange = hitsets->getHitSets(TrkrDefs::TrkrId::tpcId);
for (auto hitsetitr = hitsetrange.first;
hitsetitr != hitsetrange.second;
++hitsetitr){
std::string separator = "";
auto range = cluster_map->getClusters(hitsetitr->first);
for( auto it = range.first; it != range.second; ++it )
{
// TrkrDefs::cluskey cluskey = it->first;
TrkrCluster* cluster = it->second;
if ((std::pow((double) cluster->getPosition(0), 2) +
std::pow((double) cluster->getPosition(1), 2) +
std::pow((double) cluster->getPosition(2), 2)) > (25.0 * 25.0))
{
std::vector<double> kdhit(4);
kdhit[0] = cluster->getPosition(0);
kdhit[1] = cluster->getPosition(1);
kdhit[2] = cluster->getPosition(2);
uint64_t key = cluster->getClusKey();
std::memcpy(&kdhit[3], &key, sizeof(key));
// HINT: way to get original uint64_t value from double:
//
// LOG_DEBUG("tracking.PHTpcTrackerUtil.convert_clusters_to_hits")
// << "orig: " << cluster->getClusKey() << ", readback: " << (*((int64_t*)&kdhit[3]));
kdhits.push_back(kdhit);
}
}
}
return kdhits;
}
} // namespace PHTpcTrackerUtil