 |
ATLAS Offline Software
|
This class return all neighbours for given hash ID of the cell.
More...
#include <TileNeighbour.h>
|
int | initialize_prev_eta (unsigned int ind, const std::vector< short int > &all_cells) |
|
int | initialize_next_eta (unsigned int ind, const std::vector< short int > &all_cells) |
|
int | initialize_prev_samp (unsigned int ind, const std::vector< short int > &all_cells) |
|
int | initialize_next_samp (unsigned int ind, const std::vector< short int > &all_cells) |
|
int | initialize_prev_samp_wide (unsigned int ind, const std::vector< short int > &all_cells) |
|
int | initialize_next_samp_wide (unsigned int ind, const std::vector< short int > &all_cells) |
|
int | fill_phi_vec (std::set< std::pair< IdentifierHash, int > > &ids, unsigned int hash_max, unsigned int max_phi, MsgStream &log) |
|
|
static void | get_id (std::string &strName, Identifier &id, const Tile_Base_ID *tileID) |
|
static void | get_name (Identifier &id, std::string §ion, int &module, std::string &cell, const Tile_Base_ID *tileID, MsgStream &log, const char *end) |
|
static void | print_list (std::vector< IdentifierHash > &nb_list, const Tile_Base_ID *tileID, const IdContext *context, MsgStream &log, const char *pref, const char *suff) |
|
This class return all neighbours for given hash ID of the cell.
- Author
- Alexander Solodkov
Definition at line 41 of file TileNeighbour.h.
◆ TileNeighbour()
TileNeighbour::TileNeighbour |
( |
void |
| ) |
|
◆ ~TileNeighbour()
TileNeighbour::~TileNeighbour |
( |
void |
| ) |
|
|
virtualdefault |
◆ fill_phi_vec()
int TileNeighbour::fill_phi_vec |
( |
std::set< std::pair< IdentifierHash, int > > & |
ids, |
|
|
unsigned int |
hash_max, |
|
|
unsigned int |
max_phi, |
|
|
MsgStream & |
log |
|
) |
| |
|
private |
Definition at line 367 of file TileNeighbour.cxx.
371 if(
ids.size() != hash_max) {
372 log << MSG::ERROR <<
"fill_phi_vec "
373 <<
" Error: set size NOT EQUAL to hash max. size " <<
ids.size()
374 <<
" hash max " << hash_max
383 int dphi = hash_max / max_phi;
385 int ind_min=0, ind_max=-1,
phi, phi_prev=-1,
index=0;
387 std::set<std::pair<IdentifierHash,int> >::const_iterator
first =
ids.begin();
388 std::set<std::pair<IdentifierHash,int> >::const_iterator last =
ids.end();
391 phi = (*first).second;
392 if (
phi != phi_prev) {
405 IdentifierHash hash_id = (*first).first;
◆ get_id()
Definition at line 562 of file TileNeighbour.cxx.
564 const char *
name = strName.c_str();
570 std::string::size_type
pos = strName.find(
'-', 0 );
571 if ( std::string::npos !=
pos )
sd = -1;
else sd = 1;
573 sscanf(
name+2,
"%80d",&tw);
577 case 'a':
case 'A': sm = 0; tw -= 1;
578 if (tw<10)
se = 1;
else se = 2;
580 case 'b':
case 'B': sm = 1; tw -= 1;
581 if (tw<10)
se = 1;
else se = 2;
583 case 'c':
case 'C': sm = 1; tw -= 1;
586 case 'd':
case 'D': sm = 2; tw *= 2;
587 if (tw<10)
se = 1;
else se = 2;
590 case 'e':
case 'E': sm = 3;
592 case 1: tw = 10;
break;
593 case 2: tw = 11;
break;
594 case 3: tw = 13;
break;
595 case 4: tw = 15;
break;
599 case 's':
case 'S': sm = 0; tw -= 1;
600 if (tw<9)
se = 1;
else se = 2;
602 case 'v':
case 'V': sm = 2; tw *= 2;
603 if (tw<10)
se = 1;
else se = 2;
605 default:
se = 0; sm = -1; tw = -1;
break;
◆ get_name()
void TileNeighbour::get_name |
( |
Identifier & |
id, |
|
|
std::string & |
section, |
|
|
int & |
module, |
|
|
std::string & |
cell, |
|
|
const Tile_Base_ID * |
tileID, |
|
|
MsgStream & |
log, |
|
|
const char * |
end |
|
) |
| |
|
staticprivate |
Definition at line 611 of file TileNeighbour.cxx.
618 int tower = tileID->
tower (
id);
623 case 0: strSection =
"Online";
break;
624 case 1: strSection =
"Barrel";
break;
625 case 2: strSection =
"Ext.Barrel";
break;
626 case 3: strSection = (
sample < 3) ?
"ITC":
"Gap.Scin";
break;
627 default: strSection =
"Unknown";
631 char sd = (
side<0) ?
'-' :
'+';
635 case 0: sm=(supercell?
'S':
'A');
643 case 2: sm=(supercell?
'V':
'D');
657 else snprintf(
name,
sizeof(
name),
"%c%c%d", sm,
sd,tw);
662 << strSection <<
" " <<
module <<
" " << strCell <<
end;
◆ initialize()
int TileNeighbour::initialize |
( |
const Tile_Base_ID * |
tileID, |
|
|
const std::string & |
filename = "TileNeighbour.txt" |
|
) |
| |
Definition at line 57 of file TileNeighbour.cxx.
76 log << MSG::ERROR <<
"Could not open file " <<
file <<
endmsg;
83 unsigned int line=0, record=0;
88 std::vector<Cell> allCells;
90 while (
fin >> token ) {
92 if ( token[0] ==
'#' || token[0] ==
'!' || ( token[0] ==
'/' && token[1] ==
'/' ) ) {
97 unsigned int toklen = strlen(token) - 1 ;
102 newCell.name = token;
103 for (
unsigned int i=0;
i<4; ++
i) newCell.neighbours[
i].clear();
106 while (
fin >> token) {
107 toklen = strlen(token) - 1;
108 bool lastelem = ( token[toklen] ==
';' );
111 if (strcmp(token,
"none") != 0)
112 newCell.neighbours[
column].emplace_back(token);
122 allCells.emplace_back(std::move(newCell));
131 unsigned int curSize = allCells.size();
132 for (
unsigned int i=0;
i<curSize; ++
i) {
134 std::string::size_type
pos = allCells[
i].name.find(
'-', 0 );
135 if ( std::string::npos !=
pos ) {
139 tmpName = allCells[
i].name;
140 tmpName.replace(
pos,1,
"+");
141 newCell.name = tmpName;
143 for (
unsigned int j=0; j<4; ++j) {
144 unsigned int j1 = (j<2) ? 1-j : j;
145 unsigned int nb_size=allCells[
i].neighbours[j].size();
146 for (
unsigned int k=0;
k<nb_size; ++
k) {
147 tmpName = allCells[
i].neighbours[j][
k];
148 pos = tmpName.find(
'-', 0 );
149 if ( std::string::npos !=
pos ) {
150 tmpName.replace(
pos,1,
"+");
152 pos = tmpName.find(
'+', 0 );
153 if ( std::string::npos !=
pos ) {
154 tmpName.replace(
pos,1,
"-");
157 newCell.neighbours[
j1].push_back(tmpName);
160 allCells.push_back(std::move(newCell));
164 unsigned int nids = 0;
165 unsigned int max_phi = 0,
phi=0;
166 std::set<std::pair<IdentifierHash,int> >
ids;
171 auto first = rit.begin();
172 auto last = rit.end();
176 if (
phi > max_phi) max_phi =
phi;
183 if(!(
ids.insert(std::make_pair(hashid,
phi))).second){
184 log << MSG::ERROR <<
"init_hashes "
185 <<
" Error: duplicated id for cell id. nids= " << nids
200 curSize = allCells.size();
202 std::string nb_name[4] = {
"Prev_eta_",
"Next_eta_",
"Prev_smp_",
"Next_smp_"};
204 for (
unsigned int i=0;
i<curSize; ++
i) {
207 for (
unsigned int ii=0; ii<
m_length; ++ii) {
209 allCells[
i].index = ii;
214 for (
unsigned int j=0; j<4; ++j) {
215 unsigned int nb_size=allCells[
i].neighbours[j].size();
216 allCells[
i].neighbours_ind[j].resize(nb_size);
217 for (
unsigned int k=0;
k<nb_size; ++
k) {
218 for (
unsigned int ii=0; ii<curSize; ++ii) {
219 if ( allCells[ii].
name == allCells[
i].neighbours[j][
k]) {
220 allCells[
i].neighbours_ind[j][
k] = ii;
223 << tileID->
to_string(allCells[
i].
id,-2) <<
" "
224 << allCells[
i].hash_id <<
" "
225 << allCells[
i].index <<
" "
226 << nb_name[j] <<
k <<
" "
227 << allCells[
i].neighbours[j][
k] <<
" "
228 << allCells[
i].neighbours_ind[j][
k] <<
" "
245 std::vector<short int> neighbours[4];
246 for (
unsigned int i=0;
i<curSize; ++
i) {
247 unsigned int ind = allCells[
i].index;
249 for (
unsigned int j=0; j<4; ++j) {
250 unsigned int nb_size=allCells[
i].neighbours[j].size();
251 neighbours[j].resize(nb_size);
252 for (
unsigned int k=0;
k<nb_size; ++
k) {
253 neighbours[j][
k] = allCells[allCells[
i].neighbours_ind[j][
k]].index;
265 int sd = tileID->
side(allCells[
i].
id);
266 int tw = tileID->
tower(allCells[
i].
id);
267 int sm = tileID->
sample(allCells[
i].
id);
270 for (
unsigned int j=2; j<4; ++j) {
273 unsigned int nb_size=allCells[
i].neighbours[j].size();
274 for (
unsigned int k=0;
k<nb_size; ++
k) {
275 unsigned int new_i = allCells[
i].neighbours_ind[j][
k];
276 for (
unsigned int new_j=0; new_j<2; ++new_j) {
277 unsigned int new_nb_size=allCells[new_i].neighbours[new_j].size();
278 for (
unsigned int new_k=0; new_k<new_nb_size; ++new_k) {
279 unsigned int new_nb_i = allCells[new_i].neighbours_ind[new_j][new_k];
280 short new_nb_index = allCells[new_nb_i].index;
281 int n=neighbours[j].size()-1;
283 if (neighbours[j][
n] == new_nb_index)
287 int new_sd = tileID->
side(allCells[new_nb_i].
id);
288 int new_tw = tileID->
tower(allCells[new_nb_i].
id);
289 int new_sm = tileID->
sample(allCells[new_nb_i].
id);
290 int new_se = tileID->
section(allCells[new_nb_i].
id);
291 if ( ( new_sd !=
sd && new_tw == tw ) ||
293 ( new_se ==
se && abs(new_tw - tw) < 3 ) ||
295 ( new_sm == sm+1 && abs(new_tw - tw) < 3 ) ||
296 ( new_sm != sm+1 && abs(new_tw - tw) < 4 ) ))))) {
297 neighbours[j].push_back(new_nb_index);
312 std::vector<IdentifierHash> nb_list;
315 for (
unsigned int i=0;
i<hash_max; ++
i) {
318 tileID->
get_id (
i,
id, &context);
319 if ( tileID->
module(
id) != 0 )
continue;
321 IdentifierHash hash_id=
i;
322 nb_list.clear(); nb_list.push_back(hash_id);
326 print_list(nb_list, tileID, &context,
log,
"Prev phi ",
" ; ");
329 print_list(nb_list, tileID, &context,
log,
"Next phi ",
" ; ");
332 print_list(nb_list, tileID, &context,
log,
"Prev eta ",
" ; ");
335 print_list(nb_list, tileID, &context,
log,
"Next eta ",
" ; ");
338 print_list(nb_list, tileID, &context,
log,
"Prev smp ",
" ; ");
341 print_list(nb_list, tileID, &context,
log,
"Next smp ",
" ; ");
344 print_list(nb_list, tileID, &context,
log,
"Face 2D ",
" ; ");
347 print_list(nb_list, tileID, &context,
log,
"Corn 2D ",
" ; ");
356 print_list(nb_list, tileID, &context,
log,
"All 3DC ",
" ; ");
359 print_list(nb_list, tileID, &context,
log,
"Super 3D ",
" ; ");
◆ initialize_next_eta()
int TileNeighbour::initialize_next_eta |
( |
unsigned int |
ind, |
|
|
const std::vector< short int > & |
all_cells |
|
) |
| |
|
private |
◆ initialize_next_samp()
int TileNeighbour::initialize_next_samp |
( |
unsigned int |
ind, |
|
|
const std::vector< short int > & |
all_cells |
|
) |
| |
|
private |
◆ initialize_next_samp_wide()
int TileNeighbour::initialize_next_samp_wide |
( |
unsigned int |
ind, |
|
|
const std::vector< short int > & |
all_cells |
|
) |
| |
|
private |
◆ initialize_prev_eta()
int TileNeighbour::initialize_prev_eta |
( |
unsigned int |
ind, |
|
|
const std::vector< short int > & |
all_cells |
|
) |
| |
|
private |
◆ initialize_prev_samp()
int TileNeighbour::initialize_prev_samp |
( |
unsigned int |
ind, |
|
|
const std::vector< short int > & |
all_cells |
|
) |
| |
|
private |
◆ initialize_prev_samp_wide()
int TileNeighbour::initialize_prev_samp_wide |
( |
unsigned int |
ind, |
|
|
const std::vector< short int > & |
all_cells |
|
) |
| |
|
private |
◆ next_eta()
◆ next_phi()
◆ next_samp()
◆ next_samp_wide()
◆ prev_eta()
◆ prev_phi()
◆ prev_samp()
◆ prev_samp_wide()
◆ print_list()
Definition at line 665 of file TileNeighbour.cxx.
675 std::string strSection;
677 std::sort(nb_list.begin(), nb_list.end());
678 unsigned int size=nb_list.size();
684 for (
unsigned int j=0; j<
size; ++j) {
687 tileID->
get_id (nb_list[j],
id, context);
◆ m_cell_index
std::vector<short int> TileNeighbour::m_cell_index |
|
private |
◆ m_debug
unsigned int TileNeighbour::m_debug |
|
private |
◆ m_hashid_vec
◆ m_length
unsigned int TileNeighbour::m_length |
|
private |
◆ m_maxHash
unsigned int TileNeighbour::m_maxHash |
|
private |
◆ m_next_eta
std::vector<std::vector<short int> > TileNeighbour::m_next_eta |
|
private |
◆ m_next_samp
std::vector<std::vector<short int> > TileNeighbour::m_next_samp |
|
private |
◆ m_next_samp_wide
std::vector<std::vector<short int> > TileNeighbour::m_next_samp_wide |
|
private |
◆ m_phid_index
std::vector<short int> TileNeighbour::m_phid_index |
|
private |
◆ m_prev_eta
std::vector<std::vector<short int> > TileNeighbour::m_prev_eta |
|
private |
◆ m_prev_samp
std::vector<std::vector<short int> > TileNeighbour::m_prev_samp |
|
private |
◆ m_prev_samp_wide
std::vector<std::vector<short int> > TileNeighbour::m_prev_samp_wide |
|
private |
The documentation for this class was generated from the following files:
int initialize_next_samp(unsigned int ind, const std::vector< short int > &all_cells)
int initialize_prev_eta(unsigned int ind, const std::vector< short int > &all_cells)
Scalar phi() const
phi method
size_type cell_hash_max() const
int initialize_prev_samp_wide(unsigned int ind, const std::vector< short int > &all_cells)
std::vector< IdentifierHash > m_hashid_vec
static void get_id(std::string &strName, Identifier &id, const Tile_Base_ID *tileID)
logLevel
If HLT PSK is set on command line read it from DB instead of COOL (ATR-25974)
int side(const Identifier &id) const
int sample(const Identifier &id) const
int initialize_next_eta(unsigned int ind, const std::vector< short int > &all_cells)
int fill_phi_vec(std::set< std::pair< IdentifierHash, int > > &ids, unsigned int hash_max, unsigned int max_phi, MsgStream &log)
IMessageSvc * m_msgSvc
pointer to the message service
int tower(const Identifier &id) const
std::vector< short int > m_cell_index
size_type m_SECTION_INDEX
virtual int get_id(const IdentifierHash &hash_id, Identifier &id, const IdContext *context=0) const
create compact id from hash id (return == 0 for OK)
int initialize_next_samp_wide(unsigned int ind, const std::vector< short int > &all_cells)
int module(const Identifier &id) const
std::vector< std::vector< short int > > m_prev_eta
IdentifierHash hash(Identifier id) const
Look up the hash code corresponding to an Identifier.
std::vector< std::vector< short int > > m_next_samp
std::vector< std::vector< short int > > m_prev_samp
A Range describes the possible ranges for the field values of an ExpandedIdentifier.
int get_neighbours(const IdentifierHash &id, const LArNeighbours::neighbourOption &option, std::vector< IdentifierHash > &neighbourList) const
access to hashes for neighbours return == 0 for neighbours found option = prevInPhi,...
std::vector< std::vector< short int > > m_next_samp_wide
const HashGroup & channels() const
Return the HashGroup for channels (cells).
MultiRange m_full_cell_range
IdContext cell_context() const
id for cells
std::vector< std::vector< short int > > m_prev_samp_wide
std::string show_to_string(Identifier id, const IdContext *context=0, char sep='.') const
or provide the printout in string form
std::string to_string(const Identifier &id, int level=0) const
void sort(typename std::reverse_iterator< DataModel_detail::iterator< DVL > > beg, typename std::reverse_iterator< DataModel_detail::iterator< DVL > > end, const Compare &comp)
Specialization of sort for DataVector/List.
static std::string find_file(const std::string &logical_file_name, const std::string &search_path)
std::vector< std::vector< short int > > m_next_eta
bool is_supercell(const Identifier &id) const
static void get_name(Identifier &id, std::string §ion, int &module, std::string &cell, const Tile_Base_ID *tileID, MsgStream &log, const char *end)
std::vector< short int > m_phid_index
int section(const Identifier &id) const
size_type cardinality() const
Computes a possible cardinality from all ranges.
Identifier cell_id(const Identifier &any_id) const
This class saves the "context" of an expanded identifier (ExpandedIdentifier) for compact or hash ver...
static void print_list(std::vector< IdentifierHash > &nb_list, const Tile_Base_ID *tileID, const IdContext *context, MsgStream &log, const char *pref, const char *suff)
int initialize_prev_samp(unsigned int ind, const std::vector< short int > &all_cells)