Public Types | Public Member Functions | Private Member Functions | Private Attributes | List of all members
HalfApaChannelRanges Class Reference

#include <HalfApaChannelRanges.h>

Inheritance diagram for HalfApaChannelRanges:
IndexRangeTool

Public Types

using Name = std::string
 
using NameVector = std::vector< Name >
 
using Index = IndexRange::Index
 
using IndexVector = std::vector< Index >
 
using IndexRangeMap = std::map< Name, IndexRange >
 
- Public Types inherited from IndexRangeTool
using Name = IndexRange::Name
 

Public Member Functions

 HalfApaChannelRanges (fhicl::ParameterSet const &ps)
 
 ~HalfApaChannelRanges () override=default
 
IndexRange get (Name nam) const override
 
- Public Member Functions inherited from IndexRangeTool
virtual ~IndexRangeTool ()=default
 

Private Member Functions

void insertLen (Name nam, Index begin, Index len, Name lab, Name lab1="", Name lab2="")
 

Private Attributes

Index m_LogLevel
 
IndexVector m_ApaNumbers
 
Name m_ExtraRanges
 
IndexRangeMap m_Ranges
 
const IndexRangeTool * m_pExtraRanges =nullptr
 

Detailed Description

Definition at line 43 of file HalfApaChannelRanges.h.

Member Typedef Documentation

using HalfApaChannelRanges::Index = IndexRange::Index

Definition at line 49 of file HalfApaChannelRanges.h.

using HalfApaChannelRanges::IndexRangeMap = std::map<Name, IndexRange>

Definition at line 51 of file HalfApaChannelRanges.h.

using HalfApaChannelRanges::IndexVector = std::vector<Index>

Definition at line 50 of file HalfApaChannelRanges.h.

using HalfApaChannelRanges::Name = std::string

Definition at line 47 of file HalfApaChannelRanges.h.

using HalfApaChannelRanges::NameVector = std::vector<Name>

Definition at line 48 of file HalfApaChannelRanges.h.

Constructor & Destructor Documentation

HalfApaChannelRanges::HalfApaChannelRanges ( fhicl::ParameterSet const &  ps)

Definition at line 22 of file HalfApaChannelRanges_tool.cc.

23 : m_LogLevel(ps.get<Index>("LogLevel")),
24  m_ApaNumbers(ps.get<IndexVector>("ApaNumbers")),
25  m_ExtraRanges(ps.get<Name>("ExtraRanges")) {
26  const Name myname = "HalfApaChannelRanges::ctor: ";
27  const Index ntps = m_ApaNumbers.size();
28  Index nchaApa = 1280;
29  Index nfchau = 40;
30  Index nfchav = 40;
31  Index nfchax = 48;
32  Index nchau = 400;
33  Index nchav = 400;
34  Index nchaz = 240;
35  Index nchax = 2*nchaz;
36  Index nchai = nchau + nchav;
37  const IndexVector& apaIdx = m_ApaNumbers;
38  insertLen("all", 0, ntps*nchaApa, "All", "");
39  for ( Index itps=0; itps<ntps; ++itps ) {
40  string sitps;
41  string siapa;
42  if ( ntps > 1 ) {
43  sitps = std::to_string(itps);
44  siapa = std::to_string(apaIdx[itps]);
45  }
46  string labTps = "TPC set " + sitps;
47  string labApa = "APA " + siapa;
48  string sloc = "";
49  Index ch0 = itps*nchaApa;
50  string stps = "tps" + sitps;
51  string sapa = "apa" + siapa;
52  insertLen( stps, ch0, nchaApa, labTps, labApa);
53  insertLen( sapa, ch0, nchaApa, labApa);
54  string stpp = "tpp" + sitps;
55  Index chu0 = ch0;
56  Index chv0 = chu0 + nchau;
57  Index chx10 = chv0 + nchav;
58  Index chx20 = chx10 + nchaz;
59  Index chx0 = chx10;
60  bool beamRight = 2*(itps/2) == itps;
61  Index chz0 = beamRight ? chx20 : chx10;
62  Index chc0 = beamRight ? chx10 : chx20;
63  insertLen(stpp + "u", chu0, nchau, "TPC plane " + sitps + "u", labApa);
64  insertLen(sapa + "u", chu0, nchau, "APA plane " + siapa + "u");
65  insertLen(stpp + "v", chv0, nchav, "TPC plane " + sitps + "v", labApa);
66  insertLen(sapa + "v", chv0, nchav, "APA plane " + siapa + "v");
67  insertLen(stpp + "z1", chc0, nchaz, "TPC plane " + sitps + "z1", labApa);
68  insertLen(sapa + "z1", chc0, nchaz, "APA plane " + siapa + "z1");
69  insertLen(stpp + "z2", chz0, nchaz, "TPC plane " + sitps + "z2", labApa);
70  insertLen(sapa + "z2", chz0, nchaz, "APA plane " + siapa + "z2");
71  insertLen(stpp + "i", chu0, nchai, "TPC planes " + sitps + "i", labApa);
72  insertLen(sapa + "i", chu0, nchai, "APA planes " + siapa + "i");
73  insertLen(stpp + "x", chx0, nchax, "TPC planes " + sitps + "x", labApa);
74  insertLen(sapa + "x", chx0, nchax, "APA planes " + siapa + "x");
75  insertLen(stpp + "z", chx0, nchax, "collection planes " + sitps, labApa);
76  insertLen(sapa + "z", chx0, nchax, "collection planes " + siapa);
77  Index fchu0 = chu0;
78  Index fchv0 = chv0;
79  Index fchx0 = chx10;
80  Index ifmbu = 6;
81  Index ifmbv = 10;
82  Index ifmbx = 10;
83  // Loop over FEMBS in offline order.
84  for ( Index ifmbOff=0; ifmbOff<20; ++ifmbOff ) {
85  ostringstream ssnamu;
86  ssnamu << siapa << setfill('0') << setw(2) << ifmbu << "u";
87  string namu = ssnamu.str();
88  insertLen("femb" + namu, fchu0, nfchau, "FEMB block " + namu);
89  ostringstream ssnamv;
90  ssnamv << siapa << setfill('0') << setw(2) << ifmbv << "v";
91  string namv = ssnamv.str();
92  insertLen("femb" + namv, fchv0, nfchav, "FEMB block " + namv);
93  ostringstream ssnamx;
94  ssnamx << siapa << setfill('0') << setw(2) << ifmbx << "x";
95  string namx = ssnamx.str();
96  insertLen("femb" + namx, fchx0, nfchax, "FEMB block " + namx);
97  fchu0 += nfchau;
98  fchv0 += nfchav;
99  fchx0 += nfchax;
100  ifmbu += 1;
101  if ( ifmbu > 10 ) ifmbu = 1;
102  ifmbv -= 1;
103  if ( ifmbv == 0 ) ifmbv = 10;
104  if ( ifmbOff < 4 ) ifmbx -= 1;
105  else if ( ifmbOff == 4 ) ifmbx = 1;
106  else ifmbx += 1;
107 
108  }
109  }
110  if ( m_ExtraRanges.size() ) {
111  DuneToolManager* ptm = DuneToolManager::instance();
112  m_pExtraRanges = ptm->getShared<IndexRangeTool>(m_ExtraRanges);
113  if ( m_pExtraRanges == nullptr ) {
114  cout << myname << "WARNING: Extra range tool not found: " << m_ExtraRanges << endl;
115  }
116  }
117  if ( m_LogLevel >= 1 ) {
118  cout << myname << " LogLevel: " << m_LogLevel << endl;
119  cout << myname << " Extra range tool: " << m_ExtraRanges << endl;
120  }
121 }
IndexVector
std::vector< Index > IndexVector
Definition: test_FclIndexVectorMapTool.cxx:26
HalfApaChannelRanges::insertLen
void insertLen(Name nam, Index begin, Index len, Name lab, Name lab1="", Name lab2="")
Definition: HalfApaChannelRanges_tool.cc:183
Name
ChannelGroupService::Name Name
Definition: FixedChannelGroupService_service.cc:19
Index
unsigned int Index
Definition: sspmapmaker_v1.c:14
genie::units::ps
static constexpr double ps
Definition: Units.h:99
setw
Q_EXPORT QTSManip setw(int w)
Definition: qtextstream.h:331
HalfApaChannelRanges::m_pExtraRanges
const IndexRangeTool * m_pExtraRanges
Definition: HalfApaChannelRanges.h:70
setfill
Q_EXPORT QTSManip setfill(int f)
Definition: qtextstream.h:337
art::to_string
std::string to_string(ModuleType const mt)
Definition: ModuleType.h:34
DuneToolManager::instance
static DuneToolManager * instance(std::string fclname="", int dbg=1)
Definition: DuneToolManager.cxx:95
DuneToolManager::getShared
T * getShared(std::string name)
Definition: DuneToolManager.h:101
endl
QTextStream & endl(QTextStream &s)
Definition: qtextstream.cpp:2030
HalfApaChannelRanges::~HalfApaChannelRanges ( )
overridedefault

Member Function Documentation

IndexRange HalfApaChannelRanges::get ( Name  nam) const
overridevirtual

Implements IndexRangeTool.

Definition at line 125 of file HalfApaChannelRanges_tool.cc.

125  {
126  const Name myname = "HalfApaChannelRanges::runData: ";
127  if ( m_pExtraRanges != nullptr ) {
128  IndexRange rout = m_pExtraRanges->get(nam);
129  if ( rout.isValid() ) return rout;
130  }
131  // Special handling for online wire specifier fembAFFVCC
132  // We look up FEMB block fembAFFV and pick the channel corresponding to CC,
133  if ( nam.size() == 10 && nam.substr(0,4) == "femb" ) {
134  IndexRange fbran = get(nam.substr(0,8));
135  Index ifchmax = 0;
136  if ( fbran.isValid() ) {
137  char cpla = nam[7];
138  bool dirSame = true; // Is channel numbering dir same in FEMB and offline?
139  if ( cpla == 'u' ) {
140  ifchmax = 40;
141  dirSame = false;
142  } else if ( cpla == 'v' ) {
143  ifchmax = 40;
144  } else if ( cpla == 'x' || cpla == 'w') {
145  ifchmax = 48;
146  istringstream ssapa(nam.substr(4,1));
147  int iapa = 0;
148  ssapa >> iapa;
149  istringstream ssfmb(nam.substr(5,2));
150  int ifmb = 0;
151  ssfmb >> ifmb;
152  bool beamL = ifmb >= 1 && ifmb <= 10;
153  bool beamR = ifmb >= 11 && ifmb <= 20;
154  if ( beamL ) dirSame = false;
155  else if ( beamR ) dirSame = true;
156  else ifchmax = 0;
157  }
158  if ( ifchmax != fbran.size() ) {
159  cout << myname << "WARNING: FEMB block has unexpected size: "
160  << ifchmax << " != " << fbran.size() << endl;
161  ifchmax = 0;
162  }
163  string sch = nam.substr(8,2);
164  istringstream ssfch(sch);
165  Index ifch;
166  ssfch >> ifch;
167  if ( ifch > 0 && ifch <= ifchmax ) {
168  if ( dirSame ) --ifch;
169  else ifch = ifchmax - ifch;
170  Index ich = fbran.begin + ifch;
171  return IndexRange(nam, ich, ich+1, fbran.label(0) + sch, fbran.label(1));
172  }
173  }
174  }
175  IndexRangeMap::const_iterator iran = m_Ranges.find(nam);
176  if ( iran == m_Ranges.end() ) return IndexRange();
177  return iran->second;
178 }
IndexRange::begin
Index begin
Definition: IndexRange.h:34
Name
ChannelGroupService::Name Name
Definition: FixedChannelGroupService_service.cc:19
IndexRange::size
Index size() const
Definition: IndexRange.h:88
const_iterator
intermediate_table::const_iterator const_iterator
Definition: intermediate_table.cc:28
Index
unsigned int Index
Definition: sspmapmaker_v1.c:14
IndexRange::isValid
bool isValid() const
Definition: IndexRange.h:94
IndexRange::label
Name label(Index ilab=0) const
Definition: IndexRange.h:106
HalfApaChannelRanges::m_pExtraRanges
const IndexRangeTool * m_pExtraRanges
Definition: HalfApaChannelRanges.h:70
endl
QTextStream & endl(QTextStream &s)
Definition: qtextstream.cpp:2030
IndexRangeTool::get
virtual IndexRange get(Name nam) const =0
void HalfApaChannelRanges::insertLen ( Name  nam,
Index  begin,
Index  len,
Name  lab,
Name  lab1 = "",
Name  lab2 = "" 
)
private

Definition at line 183 of file HalfApaChannelRanges_tool.cc.

183  {
184  m_Ranges[nam] = IndexRange(nam, begin, begin+len, lab, lab1, lab2);
185 }
util::begin
decltype(auto) constexpr begin(T &&obj)
ADL-aware version of std::begin.
Definition: StdUtils.h:72

Member Data Documentation

IndexVector HalfApaChannelRanges::m_ApaNumbers
private

Definition at line 66 of file HalfApaChannelRanges.h.

Name HalfApaChannelRanges::m_ExtraRanges
private

Definition at line 67 of file HalfApaChannelRanges.h.

Index HalfApaChannelRanges::m_LogLevel
private

Definition at line 65 of file HalfApaChannelRanges.h.

const IndexRangeTool* HalfApaChannelRanges::m_pExtraRanges =nullptr
private

Definition at line 70 of file HalfApaChannelRanges.h.

IndexRangeMap HalfApaChannelRanges::m_Ranges
private

Definition at line 69 of file HalfApaChannelRanges.h.

The documentation for this class was generated from the following files:
Generated by   doxygen 1.8.11