Belle II Software  release-08-01-10
TrackSegmentFinder.cc
1 /**************************************************************************
2  * basf2 (Belle II Analysis Software Framework) *
3  * Author: The Belle II Collaboration *
4  * *
5  * See git log for contributors and copyright holders. *
6  * This file is licensed under LGPL-3.0, see LICENSE.md. *
7  **************************************************************************/
8 
9 //-----------------------------------------------------------------------------
10 // Description : A class to find Track Segments
11 //-----------------------------------------------------------------------------
12 
13 #define TRG_SHORT_NAMES
14 #define TRGCDC_SHORT_NAMES
15 
16 #include <bitset>
17 #include "TVectorD.h"
18 #include "trg/trg/Debug.h"
19 #include "trg/trg/State.h"
20 #include "trg/trg/Channel.h"
21 #include "trg/trg/Utilities.h"
22 #include "trg/cdc/TRGCDC.h"
23 #include "trg/cdc/Cell.h"
24 #include "trg/cdc/Wire.h"
25 #include "trg/cdc/WireHit.h"
26 #include "trg/cdc/Segment.h"
27 #include "trg/cdc/SegmentHit.h"
28 #include "trg/cdc/TrackSegmentFinder.h"
29 #include "framework/datastore/StoreArray.h"
30 #include "framework/datastore/RelationArray.h"
31 #include <framework/geometry/B2Vector3.h>
32 #include "cdc/dataobjects/CDCSimHit.h"
33 #include "mdst/dataobjects/MCParticle.h"
34 #include "trg/cdc/FrontEnd.h"
35 #include "trg/cdc/Merger.h"
36 
37 using namespace std;
38 
39 namespace Belle2 {
45 //...For debug : yi...
46 // vector<TRGSignalVector * > dbgIn;
47 // vector<TRGSignalVector * > dbgOut;
48 
49  TRGCDCTrackSegmentFinder::TRGCDCTrackSegmentFinder(const TRGCDC& TRGCDC,
50  bool makeRootFile,
51  bool logicLUTFlag)
52  :
53  TRGBoard("", TRGClock("", 0, 0), TRGClock("", 0, 0), TRGClock("", 0, 0),
54  TRGClock("", 0, 0)),
55  _cdc(TRGCDC),
56  m_logicLUTFlag(logicLUTFlag),
57  m_makeRootFile(makeRootFile),
58  _type(), _tosbE(), _tosbT() // 2019/07/31 by ytlai
59  {
60 
61  m_hitPatternInformation = new TClonesArray("TVectorD");
62  m_particleEfficiency = new TClonesArray("TVectorD");
63  m_tsInformation = new TClonesArray("TVectorD");
64  m_nnPatternInformation = new TClonesArray("TVectorD");
65 
66  // For ROOT file
67  TDirectory* currentDir = gDirectory;
68  if (m_makeRootFile) {
69  m_fileTSF = new TFile("TSF.root", "RECREATE");
70 
71  m_treeInputTSF = new TTree("m_treeInputTSF", "InputTSF");
72  m_treeInputTSF->Branch("hitPatternInformation", &m_hitPatternInformation,
73  32000, 0);
74 
75  m_treeOutputTSF = new TTree("m_treeOutputTSF", "OutputTSF");
76  m_treeOutputTSF->Branch("particleEfficiency", &m_particleEfficiency, 32000, 0);
77  m_treeOutputTSF->Branch("tsInformation", &m_tsInformation, 32000, 0);
78 
79  // For neural network TSF. Filled only when TSF and priority is hit.
80  m_treeNNTSF = new TTree("m_treeNNTSF", "NNTSF");
81  m_treeNNTSF->Branch("nnPatternInformation", &m_nnPatternInformation, 32000, 0);
82  }
83  currentDir->cd();
84  }
85 
86 
87  TRGCDCTrackSegmentFinder::TRGCDCTrackSegmentFinder(const TRGCDC& TRGCDC,
88  const std::string& name,
89  boardType type,
90  const TRGClock& systemClock,
91  const TRGClock& dataClock,
92  const TRGClock& userClockInput,
93  const TRGClock& userClockOutput,
94  const std::vector<TRGCDCSegment*>& tsSL)
95  : TRGBoard(name, systemClock, dataClock, userClockInput, userClockOutput),
96  _cdc(TRGCDC),
97  m_logicLUTFlag(),
98  m_rootTSFFilename(),
99  m_fileTSF(),
100  _type(type),
101  // _tisb(0),
102  _tosbE(0),
103  _tosbT(0),
104  _tsSL(tsSL)
105  {
106 
107  }
108  TRGCDCTrackSegmentFinder::~TRGCDCTrackSegmentFinder()
109  {
110 
111  delete m_nnPatternInformation;
112  delete m_tsInformation;
113  delete m_particleEfficiency;
114  delete m_hitPatternInformation;
115  delete m_treeInputTSF;
116  if (m_makeRootFile) delete m_fileTSF;
117 
118  }
119 
120  TRGCDCTrackSegmentFinder::boardType
121  TRGCDCTrackSegmentFinder:: type(void) const
122  {
123  return _type;
124  }
125  void TRGCDCTrackSegmentFinder::terminate(void)
126  {
127  if (m_makeRootFile) {
128  m_fileTSF->Write();
129  m_fileTSF->Close();
130  }
131  }
132 
133  void TRGCDCTrackSegmentFinder::doit(std::vector<TRGCDCSegment*>& tss,
134  const bool trackSegmentClockSimulation,
135  std::vector<TRGCDCSegmentHit*>& segmentHits,
136  std::vector<TRGCDCSegmentHit*>* segmentHitsSL)
137  {
138  TRGDebug::enterStage("Track Segment Finder");
139 
140  // Saves TS information
141  //saveTSInformation(tss);
142 
143  //...Store TS hits...
144  const unsigned n = tss.size();
145  for (unsigned i = 0; i < n; i++) {
146  TCSegment& s = * tss[i];
147  s.simulate(trackSegmentClockSimulation, m_logicLUTFlag,
148  TRGCDC::getTRGCDC()->getCDCHitCollectionName());
149  if (s.signal().active()) {
150  TCSHit* th = new TCSHit(s);
151  s.hit(th);
152  segmentHits.push_back(th);
153  segmentHitsSL[s.layerId()].push_back(th);
154  }
155  }
156 
157  // Save TSF results
158  saveTSFResults(segmentHitsSL);
159 
160  // Saves NNTS information. Only when ts is hit.
161  //saveNNTSInformation(tss);
162 
163  if (TRGDebug::level() > 1) {
164  cout << TRGDebug::tab() << "TS hit list" << endl;
165  string dumpOption = "trigger";
166  if (TRGDebug::level() > 2)
167  dumpOption = "detail";
168  for (unsigned i = 0; i < _cdc.nSegments(); i++) {
169  const TCSegment& s = _cdc.segment(i);
170  if (s.signal().active())
171  s.dump(dumpOption, TRGDebug::tab(4));
172  }
173 
174  cout << TRGDebug::tab() << "TS hit list (2)" << endl;
175  if (TRGDebug::level() > 2)
176  dumpOption = "detail";
177  for (unsigned i = 0; i < segmentHits.size(); i++) {
178  const TCSHit& s = * segmentHits[i];
179  s.dump(dumpOption, TRGDebug::tab(4));
180  }
181 
182  cout << TRGDebug::tab() << "TS hit list (3)" << endl;
183  if (TRGDebug::level() > 2)
184  dumpOption = "detail";
185  for (unsigned j = 0; j < _cdc.nSuperLayers(); j++) {
186  for (unsigned i = 0; i < segmentHitsSL[j].size(); i++) {
187  const vector<TCSHit*>& s = segmentHitsSL[j];
188  for (unsigned k = 0; k < s.size(); k++)
189  s[k]->dump(dumpOption, TRGDebug::tab(4));
190  }
191  }
192  }
193 
194  TRGDebug::leaveStage("Track Segment Finder");
195  }
196 
197 
198  vector <TRGSignalVector*>
199  TRGCDCTrackSegmentFinder::findTSHit(TRGSignalVector* eachInput, int tsid)
200  {
201 
202  //variables for common
203  const string na = "TSCandidate" + TRGUtilities::itostring(
204  tsid) + " in " + name();
205  TCSegment* tsi = _tsSL[tsid];
206  vector <TRGSignalVector*> result;
207 
208  //variables for EvtTime & Low pT
209  vector<bool> fTimeVect;
210  // int tmpFTime = 0 ;
211 
212  //variables for Tracker & N.N
213  vector <bool> tmpOutBool;
214 
215  TRGSignalVector* resultT = new TRGSignalVector(na, eachInput->clock(), 22);
216  TRGSignalVector* resultE = new TRGSignalVector(na, eachInput->clock(), 10);
217  TRGSignalVector* Hitmap = new TRGSignalVector(na + "HitMap",
218  eachInput->clock(), 0);
219  TRGSignalVector pTime(na + "PriorityTime", eachInput->clock(), 0);
220  TRGSignalVector fTime(na + "FastestTime", eachInput->clock(), 0);
221  for (unsigned i = 0; i < 12; i++) {
222  Hitmap->push_back((*eachInput)[i]);
223  (*Hitmap)[i].widen(16);
224  }
225  for (unsigned i = 0; i < 4; i++) {
226  pTime.push_back((*eachInput)[i + 12]);
227  fTime.push_back((*eachInput)[i + 16]);
228  }
229  for (unsigned i = 0; i < 5; i++) {
230  pTime.push_back((*eachInput)[i + 20]);
231  fTime.push_back((*eachInput)[i + 20]);
232  }
233 
234  vector <int> changeTime = Hitmap->stateChanges();
235 
236  int* LUTValue = new int[changeTime.size()];
237  if (changeTime.size()) {
238  int hitPosition = 0;
239  bool fTimeBool[10];
240  int tmpPTime = 0 ;
241  int tmpCTime = 0 ;
242  int tmpOutInt;
243  fTime.state(changeTime[0]).copy2bool(fTimeBool);
244  fTimeBool[9] = true;
245  fTimeVect.insert(fTimeVect.begin(), fTimeBool, fTimeBool + 10);
246  //tmpFTime = mkint(fTime.state(changeTime[0]));
247  bool eOUT = true;
248  for (unsigned i = 0; i < changeTime.size(); i++) {
249  LUTValue[i] = tsi->LUT()->getValue(mkint(Hitmap->state(changeTime[i])));
250 
252  if ((LUTValue[i]) && (eOUT)) {
253  resultE->set(fTimeVect, changeTime[i]);
254  eOUT = false;
255  }
256 
257  bool priority1rise = (*Hitmap)[6].riseEdge(changeTime[i]);
258  bool priority2rise = ((*Hitmap)[7].riseEdge(changeTime[i]) or
259  (*Hitmap)[8].riseEdge(changeTime[i]));
260 
262  //ready for output
263  if (priority1rise) {
264  hitPosition = 3;
265  tmpPTime = mkint(pTime.state(changeTime[i]));
266  tmpCTime = changeTime[i];
267  } else if (priority2rise) {
268  if (!hitPosition) {
269  tmpPTime = mkint(pTime.state(changeTime[i]));
270  tmpCTime = changeTime[i];
271  if ((*Hitmap)[0].state(changeTime[i])) hitPosition = 2;
272  else hitPosition = 1;
273  }
274  }
275 
276  // output selection
277  if ((hitPosition) && (LUTValue[i]) && ((changeTime[i] - tmpCTime) < 16)) {
278  tmpOutInt = tsid * pow(2, 13) + tmpPTime * pow(2, 4) + LUTValue[i] * pow(2,
279  2) + hitPosition;
280  tmpOutBool = mkbool(tmpOutInt, 22);
281  if (hitPosition == 3) {
282  if (priority1rise) resultT->set(tmpOutBool, changeTime[i]);
283  else {
284  if ((LUTValue[i] == 1) | (LUTValue[i] == 2)) {
285  if (!((LUTValue[i - 1] == 1) |
286  (LUTValue[i - 1] == 2)))
287  resultT->set(tmpOutBool, changeTime[i]);
288  } else {
289  if (!(LUTValue[i - 1])) resultT->set(tmpOutBool, changeTime[i]);
290  }
291  }
292  } else {
293  if (priority2rise) resultT->set(tmpOutBool, changeTime[i]);
294  else {
295  if ((LUTValue[i] == 1) | (LUTValue[i] == 2)) {
296  if (!((LUTValue[i - 1] == 1) |
297  (LUTValue[i - 1] == 2)))
298  resultT->set(tmpOutBool, changeTime[i]);
299  } else {
300  if (!(LUTValue[i])) resultT->set(tmpOutBool, changeTime[i]);
301  }
302  }
303  }
304  }
305 
306  }
307 
308  }
309  result.push_back(resultT);
310  result.push_back(resultE);
311 
312  // cppcheck-suppress uninitdata
313  delete[] LUTValue;
314  delete Hitmap;
315 
316  return result;
317  }
318 
319  TRGSignalVector*
320  TRGCDCTrackSegmentFinder::packerOuterTracker(vector<TRGSignalVector*>& hitList,
321  vector<int>& cList,
322  const unsigned maxHit)
323  {
324  TRGSignalVector* result =
325  new TRGSignalVector("", (* hitList[0]).clock(), 21 * maxHit);
326 
327  for (unsigned ci = 0; ci < cList.size(); ci++) {
328  unsigned cntHit = 0;
329  for (unsigned hi = 0; hi < hitList.size(); hi++) {
330  TRGState s = (* hitList[hi]).state(cList[ci]);
331  if (s.active()) {
332  if (cntHit >= maxHit) continue;
333  for (unsigned j = 0; j < 21; j++) {
334  if ((* hitList[hi])[j].state(cList[ci])) {
335  (* result)[21 * (maxHit - 1) - (cntHit * 21) + j]
336  .set(cList[ci], cList[ci] + 1);
337  }
338  }
339  if (TRGDebug::level()) {
340  TRGState t = hitList[hi]->state(cList[ci]).subset(13, 9);
341  cout << TRGDebug::tab() << " hit found : TSF out local ID="
342  << unsigned(t) << "(" << t << ")" << endl;
343  }
344 
345  ++cntHit;
346  // result->dump("", "??? ");
347  }
348  }
349  }
350 
351  return result;
352  }
353 
354  TRGSignalVector*
355  TRGCDCTrackSegmentFinder::packerOuterEvt(vector<TRGSignalVector*> hitList, vector<int> cList,
356  int maxHit)
357  {
358 
359  //TRGSignalVector * result = new TRGSignalVector("",(*hitList[0]).clock() ,N+9*maxHit);
360  TRGSignalVector* result = new TRGSignalVector("", (*hitList[0]).clock(),
361  hitList.size() + 9 * maxHit);
362 
363  for (unsigned ci = 0; ci < cList.size(); ci++) {
364  int cntHit = 0;
365  for (unsigned hi = 0; hi < hitList.size(); hi++) {
366  if ((*hitList[hi]).state(cList[ci]).active()) {
367  (*result)[9 * maxHit + hi].set(cList[ci], cList[ci] + 1);
368  if (cntHit >= maxHit) continue;
369  for (unsigned j = 0; j < (((*hitList[hi])).size() - 1); j++) {
370  if ((*hitList[hi])[j].state(cList[ci]))
371  (*result)[9 * (maxHit - 1) - (cntHit * 9) + j].set(cList[ci], cList[ci] + 1);
372  }
373  cntHit++;
374  }
375  }
376  }
377  return result;
378  }
379 
380  double
381  TRGCDCTrackSegmentFinder::mkint(TRGState bitInput)
382  {
383  double r = 0;
384  bool* binput = new bool[bitInput.size()];
385  bitInput.copy2bool(binput);
386  for (unsigned i = 0; i < bitInput.size(); i++) {
387  if (binput[i])
388  r += pow(2, i);
389  }
390  delete[] binput;
391  return r;
392  }
393 
394  vector<bool>
395  TRGCDCTrackSegmentFinder::mkbool(int N, int bitSize)
396  {
397  vector<bool> boutput;
398  boutput.resize(bitSize);
399  int tmpint = N;
400  for (unsigned i = 0; tmpint; i++) {
401  if (tmpint % 2) boutput[i] = true;
402  else boutput[i] = false;
403  tmpint = tmpint / 2;
404  }
405 
406  return boutput;
407  }
408 
409 //
410 
411  void TRGCDCTrackSegmentFinder::saveTSInformation(std::vector<TRGCDCSegment*>&
412  tss)
413  {
414 
415  TClonesArray& hitPatternInformation = *m_hitPatternInformation;
416  hitPatternInformation.Clear();
417 
418  StoreArray<CDCSimHit> SimHits;
419  if (! SimHits.getEntries()) {
420  //cout << "CDCTRGTSF !!! can not access to CDCSimHits" << endl;
421  return;
422  }
423  StoreArray<MCParticle> mcParticles;
424  if (! mcParticles.getEntries()) {
425  //cout << "CDCTRGTSF !!! can not access to MCParticles" << endl;
426  return;
427  }
428 
429  //cout<<"Save TS information"<<endl;
430  // Loop over all TSs. wireHit has only one wireSimHit. Meaning there is a particle overlap hit bug.
431  const unsigned nTSs = tss.size();
432  unsigned nHitTSs = 0;
433  for (unsigned iTS = 0; iTS < nTSs; iTS++) {
434  const TCSegment& ts = * tss[iTS];
435  const TRGCDCWire* priority;
436  const TRGCDCWire* secondPriorityL;
437  const TRGCDCWire* secondPriorityR;
438  // Find priority wires
439  if (ts.wires().size() == 15) {
440  priority = ts.wires()[0];
441  secondPriorityR = ts.wires()[1];
442  secondPriorityL = ts.wires()[2];
443  } else {
444  priority = ts.wires()[5];
445  secondPriorityR = ts.wires()[6];
446  secondPriorityL = ts.wires()[7];
447  }
448  // Find L/R, phi of priority wires
449  vector<int> priorityLRs(3);
450  vector<float> priorityPhis(3);
451  B2Vector3D posOnTrack;
452  B2Vector3D posOnWire;
453  if (priority->hit() != 0) {
454  int iSimHit = priority->hit()->iCDCSimHit();
455  priorityLRs[0] = SimHits[iSimHit]->getPosFlag();
456  posOnTrack = SimHits[iSimHit]->getPosTrack();
457  posOnWire = SimHits[iSimHit]->getPosWire();
458  priorityPhis[0] = (posOnTrack - posOnWire).Phi() + M_PI_2 -
459  posOnWire.Phi();
460  //cout<<ts.name()<<endl;
461  //cout<<"Track: "<<posOnTrack.x()<<" "<<posOnTrack.y()<<" "<<posOnTrack.z()<<endl;
462  //cout<<"Wire: "<<posOnWire.x()<<" "<<posOnWire.y()<<" "<<posOnWire.z()<<endl;
463  //cout<<"Before Phi: "<<(posOnTrack - posOnWire).Phi()<<" PosOnWirePhi: "<<posOnWire.Phi()<<" After Phi: "<<priorityPhis[0]<<endl;
464  //cout<<"LR: "<<priorityLRs[0]<<endl;
465  } else {
466  priorityLRs[0] = -1;
467  priorityPhis[0] = 9999;
468  }
469  if (secondPriorityR->hit() != 0) {
470  int iSimHit = secondPriorityR->hit()->iCDCSimHit();
471  priorityLRs[1] = SimHits[iSimHit]->getPosFlag();
472  posOnTrack = SimHits[iSimHit]->getPosTrack();
473  posOnWire = SimHits[iSimHit]->getPosWire();
474  priorityPhis[1] = (posOnTrack - posOnWire).Phi() + M_PI_2 -
475  posOnWire.Phi();
476  } else {
477  priorityLRs[1] = -1;
478  priorityPhis[1] = 9999;
479  }
480  if (secondPriorityL->hit() != 0) {
481  int iSimHit = secondPriorityL->hit()->iCDCSimHit();
482  priorityLRs[2] = SimHits[iSimHit]->getPosFlag();
483  posOnTrack = SimHits[iSimHit]->getPosTrack();
484  posOnWire = SimHits[iSimHit]->getPosWire();
485  priorityPhis[2] = (posOnTrack - posOnWire).Phi() + M_PI_2 -
486  posOnWire.Phi();
487  } else {
488  priorityLRs[2] = -1;
489  priorityPhis[2] = 9999;
490  }
491 
492  const unsigned nWires = ts.wires().size();
493  unsigned nHitWires = 0;
494  // Find TSPatternInformation for each particle
495  map<int, unsigned> particleTSPattern;
496  // Loop over wires in TS
497  for (unsigned iWire = 0; iWire < nWires; iWire++) {
498  //...Copy signal from a wire...
499  const TRGSignal& wireSignal = ts.wires()[iWire]->signal();
500  if (wireSignal.active()) ++nHitWires;
501  // Find MC particle of hit wire
502  const TRGCDCWireHit* wireHit = ts.wires()[iWire]->hit();
503  if (wireHit != 0) {
504  MCParticle* particle = SimHits[wireHit->iCDCSimHit()]->getRelatedFrom<MCParticle>();
505  int iMCParticle = (particle) ? particle->getArrayIndex() : mcParticles.getEntries();
506  // If new particle
507  if (particleTSPattern[iMCParticle] == 0) {
508  particleTSPattern[iMCParticle] = 1 << iWire;
509  } else {
510  particleTSPattern[iMCParticle] |= 1 << iWire;
511  }
512  }
513  } // Loop over wires in TS
514 
515  // Print pattern for each particle
516  //for(map<int, unsigned >::const_iterator it = particleTSPattern.begin(); it != particleTSPattern.end(); it++ ) {
517  // bitset<15> printPattern((*it).second);
518  // cout<<ts.name()<<" MC Particle: "<< (*it).first <<" pattern: "<<printPattern<<endl;
519  //}
520 
521  if (nHitWires != 0) {
522  // Ignore TSPatterns that have 2 particles passing TS.
523  if (particleTSPattern.size() == 1) {
524  map<int, unsigned>::const_iterator it = particleTSPattern.begin();
525  bitset<15> printPattern((*it).second);
526  //cout<<ts.name()<<" MC Particle: "<< (*it).first <<" pattern: "<<printPattern<<endl;
527  //cout<<(*it).first<<" "<<ts.superLayerId()<<" "<<printPattern<<endl;
528  // Save TS Pattern information
529  TVectorD tsPatternInformation(9);
530  tsPatternInformation[0] = (*it).first;
531  tsPatternInformation[1] = ts.superLayerId();
532  tsPatternInformation[2] = double((*it).second);
533  tsPatternInformation[3] = priorityLRs[0];
534  tsPatternInformation[4] = priorityPhis[0];
535  tsPatternInformation[5] = priorityLRs[1];
536  tsPatternInformation[6] = priorityPhis[1];
537  tsPatternInformation[7] = priorityLRs[2];
538  tsPatternInformation[8] = priorityPhis[2];
539  new (hitPatternInformation[nHitTSs++]) TVectorD(tsPatternInformation);
540  }
541  //cout<<ts.name()<<" has "<<nHitWires<<" hit wires."<<endl;
542  }
543  } // End of loop over all TSs
544  if (m_makeRootFile) m_treeInputTSF->Fill();
545  //cout<<"End saving TS information"<<endl;
546  } // End of save function
547 
548 
549 
550  void TRGCDCTrackSegmentFinder::saveTSFResults(std::vector<TRGCDCSegmentHit*>*
551  segmentHitsSL)
552  {
553 
554  TClonesArray& particleEfficiency = *m_particleEfficiency;
555  TClonesArray& tsInformation = *m_tsInformation;
556  particleEfficiency.Clear();
557  tsInformation.Clear();
558 
559  int mcInformation = 1;
561  StoreArray<CDCSimHit> SimHits;
562  if (! SimHits) {
563  //cout << "CDCTRGTSF !!! can not access to CDCSimHits" << endl;
564  mcInformation = 0;
565  }
566  StoreArray<MCParticle> mcParticles;
567  if (! mcParticles) {
568  //cout << "CDCTRGTSF !!! can not access to MCParticles" << endl;
569  mcInformation = 0;
570  }
571  if (mcInformation) {
572  RelationArray cdcSimHitRel(mcParticles, SimHits);
573  // Make map for hit to mcParticle
574  map<int, int> simHitsMCParticlesMap;
575  // Loop over all particles
576  for (int iPart = 0; iPart < cdcSimHitRel.getEntries(); iPart++) {
577  // Loop over all hits for particle
578  for (unsigned iHit = 0; iHit < cdcSimHitRel[iPart].getToIndices().size();
579  iHit++) {
580  //cout<<"From: "<<cdcSimHitRel[iPart].getFromIndex()<<" To: "<<cdcSimHitRel[iPart].getToIndex(iHit)<<endl;
581  simHitsMCParticlesMap[cdcSimHitRel[iPart].getToIndex(iHit)] =
582  cdcSimHitRel[iPart].getFromIndex();
583  }
584  }
586  //for(map<int, int >::iterator it = simHitsMCParticlesMap.begin(); it != simHitsMCParticlesMap.end(); it++) {
587  // cout<<"SimHit Index: "<<(*it).first<<" MCParticle Index: "<<(*it).second<<endl;
588  //}
589  // Find efficiency for each particle
590  // Find hit TS SuperLayer for particles
591  // particleNHiTS[iMCParticle] = hitTSSL
592  map<int, unsigned> particleNHitTS;
593  // Loop over all hit TSs
594  for (int iSuperLayer = 0; iSuperLayer < 9; iSuperLayer++) {
595  // map<int, bool> particleHitTS;
596  for (unsigned iTS = 0; iTS < segmentHitsSL[iSuperLayer].size(); iTS++) {
597  const TCSegment& ts = segmentHitsSL[iSuperLayer][iTS]->segment();
598  unsigned nWires = ts.wires().size();
599  for (unsigned iWire = 0; iWire < nWires; iWire++) {
600  const TRGCDCWireHit* wireHit = ts.wires()[iWire]->hit();
601  if (wireHit != 0) {
602  int iMCParticle = simHitsMCParticlesMap[wireHit->iCDCSimHit()];
603  if (particleNHitTS[iMCParticle] == 0) {
604  unsigned hitTSSL;
605  hitTSSL = 1 << iSuperLayer;
606  particleNHitTS[iMCParticle] = hitTSSL;
607  } else
608  particleNHitTS[iMCParticle] |= 1 << iSuperLayer;
609  //cout<<ts.name()<<" "<<ts.wires()[iWire]->name()<<" was hit.";
610  //cout<<" Particle was "<<simHitsMCParticlesMap[wireHit->iCDCSimHit()]<<endl;
611  } // If wire is hit
612  } // End of loop over all wires in TS
613  } // End of loop over all TS in super layer
614  } // End of loop over all hit TSs
616  //for( map<int, unsigned>::const_iterator it = particleNHitTS.begin(); it != particleNHitTS.end(); it++) {
617  // bitset<9> printSuperLayers((*it).second);
618  // cout<<"MC particle: "<<(*it).first<<" HitSuperLayers: "<<printSuperLayers<<" nCount: "<<printSuperLayers.count()<<endl;
619  //}
620 
621  // Find last CDC hit for each MC particle
622  // tsEfficiency[i][tsEfficiency,particle pT,#MCTS]
623  vector<vector<float>> tsEfficiency;
624  // Loop over all particles
625  for (int iPart = 0; iPart < cdcSimHitRel.getEntries(); iPart++) {
626  int lastWireHit = -1;
627  // Loop over all hits for particle
628  for (unsigned iHit = 0; iHit < cdcSimHitRel[iPart].getToIndices().size();
629  iHit++) {
630  int iSimHit = cdcSimHitRel[iPart].getToIndex(iHit);
631  if (SimHits[iSimHit]->getWireID().getICLayer() > lastWireHit) lastWireHit =
632  SimHits[iSimHit]->getWireID().getICLayer();
633  //cout<<"Particle: "<<cdcSimHitRel[iPart].getFromIndex()<<" CDCSimHit: "<<iSimHit<<endl;
634  //cout<<"SuperLayer: "<<SimHits[iSimHit]->getWireID().getISuperLayer()<<" wireLayer: "<<SimHits[iSimHit]->getWireID().getICLayer()<<endl;
635  }
636  //cout<<"iMCParticle: "<<cdcSimHitRel[iPart].getFromIndex()<<" Last wire Hit: "<<lastWireHit<<endl;
637  // Calculate last superlayer
638  int lastSLHit = 0;
639  if (lastWireHit >= 53) lastSLHit = 9;
640  else if (lastWireHit >= 47) lastSLHit = 8;
641  else if (lastWireHit >= 41) lastSLHit = 7;
642  else if (lastWireHit >= 35) lastSLHit = 6;
643  else if (lastWireHit >= 29) lastSLHit = 5;
644  else if (lastWireHit >= 23) lastSLHit = 4;
645  else if (lastWireHit >= 17) lastSLHit = 3;
646  else if (lastWireHit >= 11) lastSLHit = 2;
647  else if (lastWireHit >= 5) lastSLHit = 1;
648  // Get number of hit TS for particle
649  int iMCParticle = cdcSimHitRel[iPart].getFromIndex();
650  bitset<9> hitSuperLayers(particleNHitTS[iMCParticle]);
651  int numberHitSuperLayers = hitSuperLayers.count();
652  //cout<<"iMCParticle: "<< iMCParticle << " # hit TS: "<<numberHitSuperLayers<<" MC # TS: "<<lastSLHit<<endl;
653  float mcPt = mcParticles[iMCParticle]->getMomentum().Rho();
654  float efficiency;
655  if (lastSLHit == 0) efficiency = -1;
656  else efficiency = float(numberHitSuperLayers) / lastSLHit;
657  //cout<<"Efficiency: "<<float(numberHitSuperLayers)/lastSLHit<<" MC pT: "<<mcPt<<endl;
658  vector<float> tempEfficiency;
659  tempEfficiency.resize(3);
660  tempEfficiency[0] = efficiency;
661  tempEfficiency[1] = mcPt;
662  tempEfficiency[2] = lastSLHit;
663  tsEfficiency.push_back(tempEfficiency);
664  } // End of looping over all particles
666  //for(unsigned iEfficiency=0; iEfficiency<tsEfficiency.size(); iEfficiency++) {
667  // cout<<"Efficiency: "<<tsEfficiency[iEfficiency][0]<<" Pt: "<<tsEfficiency[iEfficiency][1]<<" #MCTS: "<<tsEfficiency[iEfficiency][2]<<endl;
668  //}
669 
670  // Save TS efficiency for each particle
671  for (unsigned iEfficiency = 0; iEfficiency < tsEfficiency.size();
672  iEfficiency++) {
673  TVectorD t_particleEfficiency(3);
674  t_particleEfficiency[0] = tsEfficiency[iEfficiency][0];
675  t_particleEfficiency[1] = tsEfficiency[iEfficiency][1];
676  t_particleEfficiency[2] = tsEfficiency[iEfficiency][2];
677  new (particleEfficiency[iEfficiency]) TVectorD(t_particleEfficiency);
678  }
679 
680  } // End of no MC information
681 
682  // Save TS information
683  // [FIXME] Doesn't work when second only priority is hit.
684  // Loop over all hit TSs
685  int iHitTS = 0;
686  for (int iSuperLayer = 0; iSuperLayer < 9; iSuperLayer++) {
687  for (unsigned iTS = 0; iTS < segmentHitsSL[iSuperLayer].size(); iTS++) {
688  const TCSegment& ts = segmentHitsSL[iSuperLayer][iTS]->segment();
689  //unsigned nWires = ts.wires().size();
690  unsigned iWire = 5;
691  if (iSuperLayer == 0) iWire = 0;
692  const TRGCDCWireHit* wireHit = ts.wires()[iWire]->hit();
693  if (wireHit != 0) {
694  //cout<<"[TSF]: "<<ts.name()<<" was hit at ";
695  unsigned nHits = ts.wires()[iWire]->signal().nSignals();
696  for (unsigned iHit = 0; iHit < nHits; iHit++) {
697  if (iHit % 2 == 1) continue;
698  TVectorD tempTSInformation(3);
699  tempTSInformation[0] = iSuperLayer;
700  tempTSInformation[1] = ts.localId();
701  tempTSInformation[2] = ts.wires()[iWire]->signal().stateChanges()[iHit];
702  new (tsInformation[iHitTS++]) TVectorD(tempTSInformation);
703  //cout<<ts.wires()[iWire]->signal().stateChanges()[iHit]<<", ";
704  //iHit++;
705  }
706  //cout<<endl;
707  } // If wire is hit
708  } // End of loop over all TS in super layer
709  } // End of loop over all hit TSs
710 
711  if (m_makeRootFile) m_treeOutputTSF->Fill();
712 
713  } // End of saving TSF results
714 
715 
716 
717  void TRGCDCTrackSegmentFinder::saveNNTSInformation(
718  std::vector<TRGCDCSegment*>& tss)
719  {
720 
721  StoreArray<CDCSimHit> SimHits;
722  if (! SimHits) {
723  //cout << "CDCTRGTSF !!! can not access to CDCSimHits" << endl;
724  return;
725  }
726 
727  TClonesArray& nnPatternInformation = *m_nnPatternInformation;
728  nnPatternInformation.Clear();
729 
730  // Save Timing information in ROOT file. Fill for each TS.
731  // Loop over all TSs. wireHit has only one wireSimHit. Meaning there is a particle overlap hit bug.
732  const unsigned nTSs = tss.size();
733  unsigned indexSaving = 0;
734  for (unsigned iTS = 0; iTS < nTSs; iTS++) {
735  // If TS is hit
736  const TCSegment& ts = * tss[iTS];
737  const TCSHit* tsHit = ts.hit();
738  if (tsHit) {
739  const TRGCDCWire* priority;
740  if (ts.wires().size() == 15) priority = ts.wires()[0];
741  else priority = ts.wires()[5];
742  // If priority wire is hit
743  if (priority->hit()) {
744 
745  // Calculate timeWires
746  // Fill wire timing. Not hit = 9999. Unit is ns.
747  vector<float> wireTime;
748  if (ts.superLayerId() == 0)
749  wireTime.resize(15);
750  else
751  wireTime.resize(11);
752  // Loop over all wires
753  //cout<<ts.name();
754  const unsigned nWires = ts.wires().size();
755  for (unsigned iWire = 0; iWire < nWires; iWire++) {
756  const TRGCDCWire* wire = ts.wires()[iWire];
757  const TRGCDCWireHit* wireHit = ts.wires()[iWire]->hit();
758  // If wire is hit
759  if (wireHit) {
760  // Only check first change. This could become a bug.
761  wireTime[iWire] = wire->signal().stateChanges()[0];
762  //cout<<ts.wires()[iWire]->name()<<" Clock: "<< ts.wires()[iWire]->signal().clock().frequency()<<" Drift lenght: "<<wireHit->drift()<<" Drift time: "<<ts.wires()[iWire]->signal().stateChanges()[0]<<endl;
763  } else
764  wireTime[iWire] = 9999;
765  //cout<<" "<<wire->name();
766  } // End loop over all wires
767  //cout<<endl;
769  //for(unsigned iWire = 0; iWire < wireTime.size(); iWire++) {
770  // cout<<"Wire: "<<iWire<<" Time: "<<wireTime[iWire]<<endl;
771  //}
772  // Get additional wire information for 6th layer
773  //cout<<" JB: "<<ts.wires().back()->layerId()<<" "<<ts.wires().back()->localId()<<endl;
774  int lastLayer = ts.wires().back()->layerId();
775  int lastWire = ts.wires().back()->localId();
776  int nWiresLayer = _cdc.layer(lastLayer + 1)->nCells();
777  if (nWires == 15) {
778  for (unsigned iWire = 0; iWire < 6; iWire++) {
779  int wireIndex = lastWire - 4 + iWire;
780  if (wireIndex < 0) wireIndex += nWiresLayer;
781  if (wireIndex >= nWiresLayer) wireIndex -= nWiresLayer;
782  //cout<<"Call: "<<(*_cdc.layer(lastLayer+1))[wireIndex]->localId()<<endl;
783  const TRGCDCCell* wire = (*_cdc.layer(lastLayer + 1))[wireIndex];
784  const TRGCDCCellHit* wireHit = wire->hit();
785  // If wire is hit
786  if (wireHit) {
787  // Only check first change. This could become a bug.
788  wireTime.push_back(wire->signal().stateChanges()[0]);
789  //cout<<ts.wires()[iWire]->name()<<" Clock: "<< ts.wires()[iWire]->signal().clock().frequency()<<" Drift lenght: "<<wireHit->drift()<<" Drift time: "<<ts.wires()[iWire]->signal().stateChanges()[0]<<endl;
790  } else
791  wireTime.push_back(9999);
792  } // Loop over all extra wires
793  } else {
794  for (unsigned iWire = 0; iWire < 5; iWire++) {
795  int wireIndex = lastWire - 3 + iWire;
796  if (wireIndex < 0) wireIndex += nWiresLayer;
797  if (wireIndex >= nWiresLayer) wireIndex -= nWiresLayer;
798  //cout<<"Call: "<<(*_cdc.layer(lastLayer+1))[wireIndex]->localId()<<endl;
799  const TRGCDCCell* wire = (*_cdc.layer(lastLayer + 1))[wireIndex];
800  const TRGCDCCellHit* wireHit = wire->hit();
801  // If wire is hit
802  if (wireHit) {
803  // Only check first change. This could become a bug.
804  wireTime.push_back(wire->signal().stateChanges()[0]);
805  //cout<<ts.wires()[iWire]->name()<<" Clock: "<< ts.wires()[iWire]->signal().clock().frequency()<<" Drift lenght: "<<wireHit->drift()<<" Drift time: "<<ts.wires()[iWire]->signal().stateChanges()[0]<<endl;
806  } else
807  wireTime.push_back(9999);
808  } // Loop over all extra wires
809  }
810  //ts.wires()[14]->layerId(), localId()
811  //_cdc._layers[0][localId]->signal().stateChanges()[0]
812 
813 
815  //float fastestTime = 9999;
816  //for(unsigned iWire=0; iWire<wireTime.size();iWire++) {
817  // if (fastestTime > wireTime[iWire]) fastestTime = wireTime[iWire];
818  //}
820  //for(unsigned iWire=0; iWire<wireTime.size();iWire++) {
821  // if(wireTime[iWire] != 9999) wireTime[iWire] -= fastestTime;
822  //}
823 
825  //for(unsigned iWire = 0; iWire < wireTime.size(); iWire++) {
826  // cout<<"Wire: "<<iWire<<" Time: "<<wireTime[iWire]<<endl;
827  //}
828  //cout<<ts.name()<<" is found and priority wire is hit"<<endl;
829  // Calculate mc result
830  float mcLRDriftTime = priority->signal().stateChanges()[0];
831  if (priority->hit()->mcLR()) mcLRDriftTime *= -1;
832  //cout<<ts.name()<<" LRDriftTime: "<<mcLRDriftTime<<endl;
833 
834  // Save timing information in Root
835  TVectorD t_nnPatternInformation;
836  if (ts.superLayerId() == 0) {
837  //t_nnPatternInformation.ResizeTo(17);
838  t_nnPatternInformation.ResizeTo(23);
839  t_nnPatternInformation[0] = ts.superLayerId();
840  t_nnPatternInformation[1] = mcLRDriftTime;
841  //for(unsigned iWire=0; iWire<15; iWire++){
842  for (unsigned iWire = 0; iWire < 21; iWire++)
843  t_nnPatternInformation[iWire + 2] = wireTime[iWire];
844  new (nnPatternInformation[indexSaving++]) TVectorD(t_nnPatternInformation);
845  } else {
846  //t_nnPatternInformation.ResizeTo(13);
847  t_nnPatternInformation.ResizeTo(17);
848  t_nnPatternInformation[0] = ts.superLayerId();
849  t_nnPatternInformation[1] = mcLRDriftTime;
850  //for(unsigned iWire=0; iWire<11; iWire++){
851  for (unsigned iWire = 0; iWire < 15; iWire++)
852  t_nnPatternInformation[iWire + 2] = wireTime[iWire];
853  new (nnPatternInformation[indexSaving++]) TVectorD(t_nnPatternInformation);
854  }
855 
856  } // End of if priority cell is hit
857  } // End of if TS is hit
858  } // End loop of all TSs
859 
860  if (m_makeRootFile) m_treeNNTSF->Fill();
861 
862  } // End of save function
863 
864  void
865  TRGCDCTrackSegmentFinder::push_back(const TRGCDCMerger* a)
866  {
867  std::vector<const TRGCDCMerger*>::push_back(a);
868  }
869 
870 
871  void
872  TRGCDCTrackSegmentFinder::simulateOuter(void)
873  {
874 
875  const string sn = "TSF::simulateOuter : " + name();
876  TRGDebug::enterStage(sn);
877 
878  //...Loop over mergers to create a super layer hit map...
879  for (unsigned m = 0; m < nInput(); m++) {
880  TRGSignalBundle* b = input(m)->signal();
881 
882  for (unsigned i = 0; i < 16; i++) {
883  _secMap.push_back(& ((* b)[0][0][208 + i]));
884  for (unsigned j = 0; j < 5; j++)
885  _hitMap[j].push_back(& ((* b)[0][0][j * 16 + i]));
886  for (unsigned j = 0; j < 4; j++)
887  _priMap.push_back(& ((* b)[0][0][80 + i * 4 + j]));
888  for (unsigned j = 0; j < 4; j++)
889  _fasMap.push_back(& ((* b)[0][0][144 + i * 4 + j]));
890  }
891 
892  for (unsigned i = 0; i < 4; i++)
893  _edg0Map.push_back(& ((* b)[0][0][224 + i]));
894  for (unsigned i = 0; i < 4; i++)
895  _edg1Map.push_back(& ((* b)[0][0][224 + 4 + i]));
896  for (unsigned i = 0; i < 4; i++)
897  _edg2Map.push_back(& ((* b)[0][0][224 + 8 + i]));
898  }
899 
900  //...Storage preparation...
901  const unsigned nTSF = nInput() * 16;
902  vector<TRGSignalVector*> trker[4];
903  vector<int> tsfStateChanges;
904 
905  //...Form a TSF...
906  for (unsigned t = 0; t < nTSF; t++) {
907 
908  const string n = name() + "-" + TRGUtilities::itostring(t);
909  TRGSignalVector* s = new TRGSignalVector(n, clockData());
910  _tsfIn.push_back(s);
911 
912  s->push_back(* _secMap[t]);
913 
914  if (t == 0) {
915  s->push_back(* (_hitMap[0][nTSF - 1]));
916  s->push_back(* (_hitMap[0][0]));
917  s->push_back(* (_hitMap[0][1]));
918  s->push_back(* (_hitMap[1][nTSF - 1]));
919  s->push_back(* (_hitMap[1][0]));
920  s->push_back(* (_hitMap[2][0]));
921  s->push_back(* (_hitMap[3][nTSF - 1]));
922  s->push_back(* (_hitMap[3][0]));
923  s->push_back(* (_hitMap[4][nTSF - 1]));
924  s->push_back(* (_hitMap[4][0]));
925  s->push_back(* (_hitMap[4][1]));
926  } else if (t == (nTSF - 1)) {
927  s->push_back(* (_hitMap[0][nTSF - 2]));
928  s->push_back(* (_hitMap[0][nTSF - 1]));
929  s->push_back(* (_hitMap[0][0]));
930  s->push_back(* (_hitMap[1][nTSF - 2]));
931  s->push_back(* (_hitMap[1][nTSF - 1]));
932  s->push_back(* (_hitMap[2][nTSF - 1]));
933  s->push_back(* (_hitMap[3][nTSF - 2]));
934  s->push_back(* (_hitMap[3][nTSF - 1]));
935 
936  s->push_back(* (_hitMap[4][nTSF - 2]));
937  s->push_back(* (_hitMap[4][nTSF - 1]));
938  s->push_back(* (_hitMap[4][0]));
939  } else {
940  s->push_back(* (_hitMap[0][t - 1]));
941  s->push_back(* (_hitMap[0][t]));
942  s->push_back(* (_hitMap[0][t + 1]));
943  s->push_back(* (_hitMap[1][t - 1]));
944  s->push_back(* (_hitMap[1][t]));
945  s->push_back(* (_hitMap[2][t]));
946  s->push_back(* (_hitMap[3][t - 1]));
947  s->push_back(* (_hitMap[3][t]));
948  s->push_back(* (_hitMap[4][t - 1]));
949  s->push_back(* (_hitMap[4][t]));
950  s->push_back(* (_hitMap[4][t + 1]));
951  }
952 
953  //...Priority timing...
954  s->push_back(* _priMap[t * 4 + 0]);
955  s->push_back(* _priMap[t * 4 + 1]);
956  s->push_back(* _priMap[t * 4 + 2]);
957  s->push_back(* _priMap[t * 4 + 3]);
958 
959  //...Priority timing...
960  priorityTiming(t, nTSF, * s, (* s)[6], (* s)[7], (* s)[8]);
961 
962  //...Fastest timing...
963  fastestTimingOuter(t, nTSF, * s);
964 
965  //...Clock counter is omitted...
966 
967  //...Simulate TSF...
968  vector<TRGSignalVector*> result = simulateOuter(s, t);
969  _toBeDeleted.push_back(result[1]); // Event timing omitted now
970 
971  TRGSignalVector* forTracker0 = result[0];
972  TRGSignalVector* forTracker1 = new TRGSignalVector(* forTracker0);
973 
974  _tsfOut.push_back(forTracker0);
975  _tsfOut.push_back(forTracker1);
976 
977 
978  //...State change list...
979  vector<int> sc = forTracker0->stateChanges();
980  for (unsigned i = 0; i < sc.size(); i++) {
981  bool skip = false;
982  for (unsigned j = 0; j < tsfStateChanges.size(); j++) {
983  if (tsfStateChanges[j] == sc[i]) {
984  skip = true;
985  break;
986  }
987  }
988  if (! skip) tsfStateChanges.push_back(sc[i]);
989  }
990 
991  //...Store it for each tracker division...
992  const unsigned pos = t / (nTSF / 4);
993  if (pos == 0) {
994  addID(* forTracker0, t + (nTSF / 4));
995  addID(* forTracker1, t);
996  trker[3].push_back(forTracker0);
997  trker[0].push_back(forTracker1);
998  } else if (pos == 1) {
999  addID(* forTracker0, t);
1000  addID(* forTracker1, t - 1 * (nTSF / 4));
1001  trker[0].push_back(forTracker0);
1002  trker[1].push_back(forTracker1);
1003  } else if (pos == 2) {
1004  addID(* forTracker0, t - 1 * (nTSF / 4));
1005  addID(* forTracker1, t - 2 * (nTSF / 4));
1006  trker[1].push_back(forTracker0);
1007  trker[2].push_back(forTracker1);
1008  } else {
1009  addID(* forTracker0, t - 2 * (nTSF / 4));
1010  addID(* forTracker1, t - 3 * (nTSF / 4));
1011  trker[2].push_back(forTracker0);
1012  trker[3].push_back(forTracker1);
1013  }
1014 
1015  if (TRGDebug::level())
1016  if (forTracker0->active())
1017  cout << TRGDebug::tab() << "hit found : TSF out="
1018  << t << endl;
1019  }
1020 
1021  //...Sort state changes...
1022  std::sort(tsfStateChanges.begin(), tsfStateChanges.end());
1023 
1024  //...Output for 2D...
1025  for (unsigned i = 0; i < 4; i++) {
1026  string n = name() + "-trker" + TRGUtilities::itostring(i);
1027  TRGSignalVector* tOut = packerOuterTracker(trker[i],
1028  tsfStateChanges,
1029  20);
1030  tOut->name(n);
1031  TRGSignalBundle* b = new TRGSignalBundle(n, clockData());
1032  b->push_back(tOut);
1033  output(i)->signal(b);
1034  _toBeDeleted.push_back(tOut);
1035 
1036  if (TRGCDC::getTRGCDC()->firmwareSimulationMode() & 0x4)
1037  b->dumpCOE("",
1038  TRGCDC::getTRGCDC()->firmwareSimulationStartDataClock(),
1039  TRGCDC::getTRGCDC()->firmwareSimulationStopDataClock());
1040 
1041  // b->dump();
1042  }
1043 
1044  TRGDebug::leaveStage(sn);
1045  }
1046 
1047  vector <TRGSignalVector*>
1048  TRGCDCTrackSegmentFinder::simulateOuter(TRGSignalVector* in, unsigned tsid)
1049  {
1050 
1051  //variables for common
1052  const string na = "TSCandidate" + TRGUtilities::itostring(tsid) + " in " +
1053  name();
1054  TCSegment* tsi = _tsSL[tsid];
1055  vector <TRGSignalVector*> result;
1056 
1057  //variables for EvtTime & Low pT
1058  vector<bool> fTimeVect;
1059  // int tmpFTime = 0 ;
1060 
1061  //variables for Tracker & N.N
1062  vector <bool> tmpOutBool;
1063 
1064  //iwTRGSignalVector* resultT = new TRGSignalVector(na, in->clock(),22);
1065  TRGSignalVector* resultT = new TRGSignalVector(na, in->clock(), 13);
1066  TRGSignalVector* resultE = new TRGSignalVector(na, in->clock(), 10);
1067  TRGSignalVector* Hitmap = new TRGSignalVector(na + "HitMap", in->clock(), 0);
1068  TRGSignalVector pTime(na + "PriorityTime", in->clock(), 0);
1069  TRGSignalVector fTime(na + "FastestTime", in->clock(), 0);
1070  for (unsigned i = 0; i < 12; i++) {
1071  Hitmap->push_back((*in)[i]);
1072  (*Hitmap)[i].widen(16);
1073  }
1074  for (unsigned i = 0; i < 4; i++) {
1075  pTime.push_back((*in)[i + 12]);
1076  fTime.push_back((*in)[i + 16]);
1077  }
1078 
1079  //...Clock counter...
1080  const TRGSignalVector& cc = in->clock().clockCounter();
1081  for (unsigned i = 0; i < 5; i++) {
1082  pTime.push_back(cc[i]);
1083  fTime.push_back(cc[i]);
1084  }
1085 
1086  vector <int> changeTime = Hitmap->stateChanges();
1087 
1088  int* LUTValue = new int[changeTime.size()];
1089  if (changeTime.size()) {
1090  int hitPosition = 0;
1091  bool fTimeBool[10];
1092  int tmpPTime = 0 ;
1093  int tmpCTime = 0 ;
1094  int tmpOutInt;
1095  fTime.state(changeTime[0]).copy2bool(fTimeBool);
1096  fTimeBool[9] = true;
1097  fTimeVect.insert(fTimeVect.begin(), fTimeBool, fTimeBool + 10);
1098  //tmpFTime = mkint(fTime.state(changeTime[0]));
1099  bool eOUT = true;
1100  for (unsigned i = 0; i < changeTime.size(); i++) {
1101  LUTValue[i] = tsi->LUT()->getValue(mkint(Hitmap->state(changeTime[i])));
1102 
1104  if ((LUTValue[i]) && (eOUT)) {
1105  resultE->set(fTimeVect, changeTime[i]);
1106  eOUT = false;
1107  }
1108 
1109  bool priority1rise = (*Hitmap)[6].riseEdge(changeTime[i]);
1110  bool priority2rise = ((*Hitmap)[7].riseEdge(changeTime[i]) or
1111  (*Hitmap)[8].riseEdge(changeTime[i]));
1112 
1114  //ready for output
1115  if (priority1rise) {
1116  hitPosition = 3;
1117  tmpPTime = mkint(pTime.state(changeTime[i]));
1118  tmpCTime = changeTime[i];
1119  } else if (priority2rise) {
1120  if (!hitPosition) {
1121  tmpPTime = mkint(pTime.state(changeTime[i]));
1122  tmpCTime = changeTime[i];
1123  if ((*Hitmap)[0].state(changeTime[i])) hitPosition = 2;
1124  else hitPosition = 1;
1125  }
1126  }
1127 
1128  // output selection
1129  if ((hitPosition) && (LUTValue[i]) && ((changeTime[i] - tmpCTime) < 16)) {
1130  //iw tmpOutInt = tsid * pow(2, 13) + tmpPTime * pow(2, 4) +
1131  //iw LUTValue[i] * pow(2,2) + hitPosition;
1132  tmpOutInt = tmpPTime * pow(2, 4) +
1133  LUTValue[i] * pow(2, 2) + hitPosition;
1134  tmpOutBool = mkbool(tmpOutInt, 13); // ID removed : iw
1135  if (hitPosition == 3) {
1136  if (priority1rise) resultT->set(tmpOutBool, changeTime[i]);
1137  else {
1138  if ((LUTValue[i] == 1) | (LUTValue[i] == 2)) {
1139  if (!((LUTValue[i - 1] == 1) |
1140  (LUTValue[i - 1] == 2)))
1141  resultT->set(tmpOutBool, changeTime[i]);
1142  } else {
1143  if (!(LUTValue[i - 1])) resultT->set(tmpOutBool, changeTime[i]);
1144  }
1145  }
1146  } else {
1147  if (priority2rise) resultT->set(tmpOutBool, changeTime[i]);
1148  else {
1149  if ((LUTValue[i] == 1) | (LUTValue[i] == 2)) {
1150  if (!((LUTValue[i - 1] == 1) |
1151  (LUTValue[i - 1] == 2)))
1152  resultT->set(tmpOutBool, changeTime[i]);
1153  } else {
1154  if (!(LUTValue[i])) resultT->set(tmpOutBool, changeTime[i]);
1155  }
1156  }
1157  }
1158  }
1159 
1160  }
1161 
1162  }
1163  result.push_back(resultT);
1164  result.push_back(resultE);
1165 
1166  // cppcheck-suppress uninitdata
1167  delete[] LUTValue;
1168  delete Hitmap;
1169 
1170  return result;
1171  }
1172 
1173  void
1174  TRGCDCTrackSegmentFinder::priorityTiming(unsigned t,
1175  const unsigned nTSF,
1176  TRGSignalVector& s,
1177  const TRGSignal& center,
1178  const TRGSignal& right,
1179  const TRGSignal& left) const
1180  {
1181 
1182  unsigned rem = t % 16;
1183 
1184  if (rem != 0) {
1185  s.push_back(* _priMap[t * 4 + 0]);
1186  s.push_back(* _priMap[t * 4 + 1]);
1187  s.push_back(* _priMap[t * 4 + 2]);
1188  s.push_back(* _priMap[t * 4 + 3]);
1189  return;
1190  }
1191 
1192  //...Prepare signal vector for timing cells...
1193  TRGSignalVector tc("timing cell hit", clockData());
1194  tc.push_back(center);
1195  tc.push_back(right);
1196  tc.push_back(left);
1197 
1198  //...No timing hit case...
1199  if (! tc.active()) {
1200  s.push_back(* _priMap[t * 4 + 0]);
1201  s.push_back(* _priMap[t * 4 + 1]);
1202  s.push_back(* _priMap[t * 4 + 2]);
1203  s.push_back(* _priMap[t * 4 + 3]);
1204  return;
1205  }
1206 
1207  const string sn = "TSF priority timing outer";
1208  TRGDebug::enterStage(sn);
1209 
1210  //...Prepare signal vector for each timing cells...
1211  TRGSignalVector t0("t0", clockData()); // center
1212  t0.push_back(* _priMap[t * 4 + 0]);
1213  t0.push_back(* _priMap[t * 4 + 1]);
1214  t0.push_back(* _priMap[t * 4 + 2]);
1215  t0.push_back(* _priMap[t * 4 + 3]);
1216 
1217  TRGSignalVector t1("t1", clockData()); // left
1218  unsigned p = t / 16;
1219  if (t == 0)
1220  p = nTSF / 16 - 1;
1221  else p = p - 1;
1222  t1.push_back(* _edg0Map[p * 4 + 0]);
1223  t1.push_back(* _edg0Map[p * 4 + 1]);
1224  t1.push_back(* _edg0Map[p * 4 + 2]);
1225  t1.push_back(* _edg0Map[p * 4 + 3]);
1226 
1227  //...Loop over state changes...
1228  vector<TRGState*> outputTimingStates;
1229  vector<int> clkStates;
1230  vector<int> sc = tc.stateChanges();
1231  vector<bool> sb;
1232  for (unsigned i = 0; i < sc.size(); i++) {
1233  int clk = sc[i];
1234  TRGState stt = tc.state(clk);
1235 
1236  //...Priority cell hit case...
1237  bool sel = true; // true=t0, false=t1
1238  if (stt[0])
1239  sel = true;
1240  else if (stt[1] && (! stt[2]))
1241  sel = false;
1242  else if ((! stt[1]) && stt[2])
1243  sel = true;
1244  else { // (stt[1] && stt[2])
1245  if (unsigned(t0.state(clk)) < unsigned(t1.state(clk)))
1246  sel = true;
1247  else
1248  sel = false;
1249  }
1250 
1251  if (sel) {
1252  outputTimingStates.push_back(new TRGState(t0.state(clk)));
1253  sb.push_back(false);
1254  } else {
1255  outputTimingStates.push_back(new TRGState(t1.state(clk)));
1256  sb.push_back(true);
1257  }
1258  clkStates.push_back(clk);
1259  }
1260 
1261  //...Create signals...
1262  const unsigned n = outputTimingStates.size();
1263  TRGSignalVector st("pri.timing", clockData(), 4);
1264  for (unsigned i = 0; i < n; i++) {
1265  if (sb[i])
1266  s[0].unset(clkStates[i], clkStates[i] + 1);
1267  st.set(* outputTimingStates[i], clkStates[i]);
1268  delete outputTimingStates[i];
1269  }
1270 
1271  //...Store signals...
1272  s.push_back(st[0]);
1273  s.push_back(st[1]);
1274  s.push_back(st[2]);
1275  s.push_back(st[3]);
1276 
1277  if (TRGDebug::level()) {
1278  tc.dump("", TRGDebug::tab() + " ");
1279  t0.dump("", TRGDebug::tab() + " ");
1280  t1.dump("", TRGDebug::tab() + " ");
1281  st.dump("", TRGDebug::tab() + " ");
1282  s.dump("", TRGDebug::tab() + " ");
1283  }
1284 
1285  TRGDebug::leaveStage(sn);
1286  }
1287 
1288  void
1289  TRGCDCTrackSegmentFinder::fastestTimingOuter(unsigned t,
1290  const unsigned nTSF,
1291  TRGSignalVector& s) const
1292  {
1293 
1294  const unsigned rem = t % 16;
1295 
1296  if ((rem != 0) && (rem != 15)) {
1297  s.push_back(* _fasMap[t * 4 + 0]);
1298  s.push_back(* _fasMap[t * 4 + 1]);
1299  s.push_back(* _fasMap[t * 4 + 2]);
1300  s.push_back(* _fasMap[t * 4 + 3]);
1301  return;
1302  }
1303 
1304  //...Check hit map if there is a hit...
1305  bool hit = false;
1306  for (unsigned i = 0; i < 11; i++) {
1307  if (s[i + 1].active()) {
1308  hit = true;
1309  break;
1310  }
1311  }
1312 
1313  //...No hit case...
1314  if (! hit) {
1315  s.push_back(* _fasMap[t * 4 + 0]);
1316  s.push_back(* _fasMap[t * 4 + 1]);
1317  s.push_back(* _fasMap[t * 4 + 2]);
1318  s.push_back(* _fasMap[t * 4 + 3]);
1319  return;
1320  }
1321 
1322  const string sn = "TSF fastest timing outer";
1323  TRGDebug::enterStage(sn);
1324 
1325  //...Prepare timing signal vectors
1326  TRGSignalVector t0("t0", clockData()); // main
1327  t0.push_back(* _fasMap[t * 4 + 0]);
1328  t0.push_back(* _fasMap[t * 4 + 1]);
1329  t0.push_back(* _fasMap[t * 4 + 2]);
1330  t0.push_back(* _fasMap[t * 4 + 3]);
1331 
1332  TRGSignalVector t1("t1", clockData()); // edge
1333  TRGSignal ht0("t0 hit", clockData());
1334  TRGSignal ht1("t1 hit", clockData());
1335  if (rem == 0) {
1336  unsigned n = t / 16;
1337  if (n == 0)
1338  n = nTSF / 16 - 1;
1339  else
1340  --n;
1341 
1342  t1.push_back(* _edg2Map[n * 4 + 0]);
1343  t1.push_back(* _edg2Map[n * 4 + 1]);
1344  t1.push_back(* _edg2Map[n * 4 + 2]);
1345  t1.push_back(* _edg2Map[n * 4 + 3]);
1346 
1347  ht0 = s[2];
1348  ht0 |= s[3];
1349  ht0 |= s[5];
1350  ht0 |= s[6];
1351  ht0 |= s[8];
1352  ht0 |= s[10];
1353  ht0 |= s[11];
1354 
1355  ht1 = s[1];
1356  ht1 |= s[4];
1357  ht1 |= s[7];
1358  ht1 |= s[9];
1359  } else {
1360  unsigned n = t / 16 + 1;
1361  if (n >= nTSF / 16)
1362  n = 0;
1363 
1364  t1.push_back(* _edg1Map[n * 4 + 0]);
1365  t1.push_back(* _edg1Map[n * 4 + 1]);
1366  t1.push_back(* _edg1Map[n * 4 + 2]);
1367  t1.push_back(* _edg1Map[n * 4 + 3]);
1368 
1369  ht0 = s[1];
1370  ht0 |= s[2];
1371  ht0 |= s[4];
1372  ht0 |= s[5];
1373  ht0 |= s[6];
1374  ht0 |= s[7];
1375  ht0 |= s[8];
1376  ht0 |= s[9];
1377  ht0 |= s[10];
1378 
1379  ht1 = s[3];
1380  ht1 |= s[11];
1381  }
1382 
1383  //...State changes...
1384  vector<int> sc = t0.stateChanges();
1385  vector<int> tmp = t1.stateChanges();
1386  sc.insert(sc.end(), tmp.begin(), tmp.end());
1387  std::sort(sc.begin(), sc.end());
1388 
1389  //...Loop over state changes...
1390  TRGSignalVector tm("fastest", clockData(), 4);
1391  int last = clockData().max();
1392  for (unsigned i = 0; i < sc.size(); i++) {
1393  if (sc[i] == last)
1394  continue;
1395 
1396  int clk = sc[i];
1397 
1398  TRGState ts0 = t0.state(clk);
1399  TRGState ts1 = t1.state(clk);
1400  unsigned tm0 = unsigned(ts0);
1401  unsigned tm1 = unsigned(ts1);
1402  bool th0 = ht0.state(clk);
1403  bool th1 = ht1.state(clk);
1404 
1405  if ((! th0) && (! th1))
1406  continue;
1407  else if (th0 && th1) {
1408  if (tm1 < tm0)
1409  tm.set(ts1, clk);
1410  else
1411  tm.set(ts0, clk);
1412  } else if (th0)
1413  tm.set(ts0, clk);
1414  else if (th1)
1415  tm.set(ts1, clk);
1416 
1417  last = clk;
1418  }
1419 
1420  //...Store signals...
1421  s.push_back(tm[0]);
1422  s.push_back(tm[1]);
1423  s.push_back(tm[2]);
1424  s.push_back(tm[3]);
1425 
1426  if (TRGDebug::level()) {
1427  ht0.name("t0 hit:" + ht0.name());
1428  ht0.dump("", TRGDebug::tab() + " ");
1429  t0.dump("", TRGDebug::tab() + " ");
1430  ht1.name("t1 hit:" + ht1.name());
1431  ht1.dump("", TRGDebug::tab() + " ");
1432  t1.dump("", TRGDebug::tab() + " ");
1433  tm.dump("", TRGDebug::tab() + " ");
1434  s.dump("", TRGDebug::tab() + " ");
1435  }
1436 
1437  TRGDebug::leaveStage(sn);
1438  }
1439 
1440  void
1441  TRGCDCTrackSegmentFinder::simulateInner(void)
1442  {
1443 
1444  const string sn = "TSF::simulateInner : " + name();
1445  TRGDebug::enterStage(sn);
1446 
1447  //...Output for 2D : empty bundle temporary...
1448  for (unsigned i = 0; i < 4; i++) {
1449  string n = name() + "-trker" + TRGUtilities::itostring(i);
1450  TRGSignalVector* dummy = new TRGSignalVector(n, clockData(), 420);
1451  TRGSignalBundle* b = new TRGSignalBundle(n, clockData());
1452  b->push_back(dummy);
1453  output(i)->signal(b);
1454  _toBeDeleted.push_back(dummy);
1455 
1456  if (TRGCDC::getTRGCDC()->firmwareSimulationMode() & 0x4)
1457  b->dumpCOE("",
1458  TRGCDC::getTRGCDC()->firmwareSimulationStartDataClock(),
1459  TRGCDC::getTRGCDC()->firmwareSimulationStopDataClock());
1460  }
1461 
1462  TRGDebug::leaveStage(sn);
1463  return;
1464 
1465 
1466  //...Loop over mergers to create a super layer hit map...
1467  for (unsigned m = 0; m < nInput(); m++) {
1468  TRGSignalBundle* b = input(m)->signal();
1469 
1470  for (unsigned i = 0; i < 16; i++) {
1471  _secMap.push_back(& ((* b)[0][0][208 + i]));
1472  for (unsigned j = 0; j < 5; j++)
1473  _hitMap[j].push_back(& ((* b)[0][0][j * 16 + i]));
1474  for (unsigned j = 0; j < 4; j++)
1475  _priMap.push_back(& ((* b)[0][0][80 + i * 4 + j]));
1476  for (unsigned j = 0; j < 4; j++)
1477  _fasMap.push_back(& ((* b)[0][0][144 + i * 4 + j]));
1478  }
1479 
1480  for (unsigned i = 0; i < 4; i++)
1481  _edg0Map.push_back(& ((* b)[0][0][224 + i]));
1482  for (unsigned i = 0; i < 4; i++)
1483  _edg1Map.push_back(& ((* b)[0][0][224 + 4 + i]));
1484  for (unsigned i = 0; i < 4; i++)
1485  _edg2Map.push_back(& ((* b)[0][0][224 + 8 + i]));
1486  for (unsigned i = 0; i < 4; i++)
1487  _edg3Map.push_back(& ((* b)[0][0][224 + 12 + i]));
1488  for (unsigned i = 0; i < 4; i++)
1489  _edg4Map.push_back(& ((* b)[0][0][224 + 16 + i]));
1490  }
1491 
1492  //...Storage preparation...
1493  const unsigned nTSF = nInput() * 16;
1494  vector<TRGSignalVector*> trker[4];
1495  vector<int> tsfStateChanges;
1496 
1497  //...Form a TSF...
1498  for (unsigned t = 0; t < nTSF; t++) {
1499 
1500  const string n = name() + "-" + TRGUtilities::itostring(t);
1501  TRGSignalVector* s = new TRGSignalVector(n, clockData());
1502  _tsfIn.push_back(s);
1503 
1504  s->push_back(* _secMap[t]);
1505 
1506  if (t == 0) {
1507  s->push_back(* (_hitMap[0][0]));
1508  s->push_back(* (_hitMap[1][nTSF - 1]));
1509  s->push_back(* (_hitMap[1][0]));
1510  s->push_back(* (_hitMap[2][nTSF - 1]));
1511  s->push_back(* (_hitMap[2][0]));
1512  s->push_back(* (_hitMap[2][1]));
1513  s->push_back(* (_hitMap[3][nTSF - 2]));
1514  s->push_back(* (_hitMap[3][nTSF - 1]));
1515  s->push_back(* (_hitMap[3][0]));
1516  s->push_back(* (_hitMap[3][1]));
1517  s->push_back(* (_hitMap[4][nTSF - 2]));
1518  s->push_back(* (_hitMap[4][nTSF - 1]));
1519  s->push_back(* (_hitMap[4][0]));
1520  s->push_back(* (_hitMap[4][1]));
1521  s->push_back(* (_hitMap[4][2]));
1522  } else if (t == 1) {
1523  s->push_back(* (_hitMap[0][1]));
1524  s->push_back(* (_hitMap[1][0]));
1525  s->push_back(* (_hitMap[1][1]));
1526  s->push_back(* (_hitMap[2][0]));
1527  s->push_back(* (_hitMap[2][1]));
1528  s->push_back(* (_hitMap[2][2]));
1529  s->push_back(* (_hitMap[3][nTSF - 1]));
1530  s->push_back(* (_hitMap[3][0]));
1531  s->push_back(* (_hitMap[3][1]));
1532  s->push_back(* (_hitMap[3][2]));
1533  s->push_back(* (_hitMap[4][nTSF - 1]));
1534  s->push_back(* (_hitMap[4][0]));
1535  s->push_back(* (_hitMap[4][1]));
1536  s->push_back(* (_hitMap[4][2]));
1537  s->push_back(* (_hitMap[4][3]));
1538  } else if (t == (nTSF - 2)) {
1539  s->push_back(* (_hitMap[0][t]));
1540  s->push_back(* (_hitMap[1][t - 11]));
1541  s->push_back(* (_hitMap[1][t]));
1542  s->push_back(* (_hitMap[2][t - 1]));
1543  s->push_back(* (_hitMap[2][t]));
1544  s->push_back(* (_hitMap[2][0]));
1545  s->push_back(* (_hitMap[3][t - 2]));
1546  s->push_back(* (_hitMap[3][t - 1]));
1547  s->push_back(* (_hitMap[3][t]));
1548  s->push_back(* (_hitMap[3][0]));
1549  s->push_back(* (_hitMap[4][t - 2]));
1550  s->push_back(* (_hitMap[4][t - 1]));
1551  s->push_back(* (_hitMap[4][t]));
1552  s->push_back(* (_hitMap[4][t + 1]));
1553  s->push_back(* (_hitMap[4][0]));
1554  } else if (t == (nTSF - 1)) {
1555  s->push_back(* (_hitMap[0][t]));
1556  s->push_back(* (_hitMap[1][t - 11]));
1557  s->push_back(* (_hitMap[1][t]));
1558  s->push_back(* (_hitMap[2][t - 1]));
1559  s->push_back(* (_hitMap[2][t]));
1560  s->push_back(* (_hitMap[2][0]));
1561  s->push_back(* (_hitMap[3][t - 2]));
1562  s->push_back(* (_hitMap[3][t - 1]));
1563  s->push_back(* (_hitMap[3][t]));
1564  s->push_back(* (_hitMap[3][0]));
1565  s->push_back(* (_hitMap[4][t - 2]));
1566  s->push_back(* (_hitMap[4][t - 1]));
1567  s->push_back(* (_hitMap[4][t]));
1568  s->push_back(* (_hitMap[4][0]));
1569  s->push_back(* (_hitMap[4][1]));
1570  } else {
1571  s->push_back(* (_hitMap[0][t]));
1572  s->push_back(* (_hitMap[1][t - 11]));
1573  s->push_back(* (_hitMap[1][t]));
1574  s->push_back(* (_hitMap[2][t - 1]));
1575  s->push_back(* (_hitMap[2][t]));
1576  s->push_back(* (_hitMap[2][0]));
1577  s->push_back(* (_hitMap[3][t - 2]));
1578  s->push_back(* (_hitMap[3][t - 1]));
1579  s->push_back(* (_hitMap[3][t]));
1580  s->push_back(* (_hitMap[3][0]));
1581  s->push_back(* (_hitMap[4][t - 2]));
1582  s->push_back(* (_hitMap[4][t - 1]));
1583  s->push_back(* (_hitMap[4][t]));
1584  s->push_back(* (_hitMap[4][t + 1]));
1585  s->push_back(* (_hitMap[4][t + 2]));
1586  }
1587 
1588  //...Priority timing...
1589  s->push_back(* _priMap[t * 4 + 0]);
1590  s->push_back(* _priMap[t * 4 + 1]);
1591  s->push_back(* _priMap[t * 4 + 2]);
1592  s->push_back(* _priMap[t * 4 + 3]);
1593 
1594  //...Priority timing...
1595  priorityTiming(t, nTSF, * s, (* s)[1], (* s)[2], (* s)[3]);
1596 
1597  //...Fastest timing...
1598  fastestTimingInner(t, nTSF, * s);
1599 
1600  //...Clock counter is omitted...
1601 
1602  //...Simulate TSF...
1603  vector<TRGSignalVector*> result = simulateInner(* s, t);
1604 
1605  TRGSignalVector* forTracker = result[0];
1606  _tsfOut.push_back(forTracker);
1607  _toBeDeleted.push_back(result[1]);
1608 
1609  //...State change list...
1610  vector<int> sc = forTracker->stateChanges();
1611  for (unsigned i = 0; i < sc.size(); i++) {
1612  bool skip = false;
1613  for (unsigned j = 0; j < tsfStateChanges.size(); j++) {
1614  if (tsfStateChanges[j] == sc[i]) {
1615  skip = true;
1616  break;
1617  }
1618  }
1619  if (! skip) tsfStateChanges.push_back(sc[i]);
1620  }
1621 
1622  //...Store it for each tracker division...
1623  const unsigned pos = t / (nTSF / 4);
1624  if (pos == 0) {
1625  trker[0].push_back(forTracker);
1626  trker[1].push_back(forTracker);
1627  } else if (pos == 1) {
1628  trker[1].push_back(forTracker);
1629  trker[2].push_back(forTracker);
1630  } else if (pos == 2) {
1631  trker[2].push_back(forTracker);
1632  trker[3].push_back(forTracker);
1633  } else {
1634  trker[3].push_back(forTracker);
1635  trker[0].push_back(forTracker);
1636  }
1637 
1638  //...Test...
1639  // if (TRGDebug::level() && t < 16) {
1640  // bool ok = ((* result[0]) == (* dbgOut[t])); // track part only
1641 
1642  // if (ok) {
1643  // cout << TRGDebug::tab() << name() << "...Comparison OK"
1644  // << endl;
1645  // }
1646  // else {
1647  // cout << TRGDebug::tab() << name()
1648  // << "...Comparison is not OK" << endl;
1649  // dbgOut[t]->dump("", "kt :");
1650  // result[0]->dump("", "2d :");
1651  // dbgIn[t]->dump("", "kt i:");
1652  // s->dump("", "2d i:");
1653  // }
1654  // }
1655  }
1656 
1657  //...Sort state changes...
1658  std::sort(tsfStateChanges.begin(), tsfStateChanges.end());
1659 
1660  //...Output for 2D...
1661  for (unsigned i = 0; i < 4; i++) {
1662  TRGSignalVector* tOut = packerOuterTracker(trker[i],
1663  tsfStateChanges,
1664  20);
1665  string n = name() + "trker" + TRGUtilities::itostring(i);
1666  TRGSignalBundle* b = new TRGSignalBundle(n, clockData());
1667  b->push_back(tOut);
1668  output(i)->signal(b);
1669  }
1670 
1671  TRGDebug::leaveStage(sn);
1672  }
1673 
1674  void
1675  TRGCDCTrackSegmentFinder::fastestTimingInner(unsigned t,
1676  const unsigned nTSF,
1677  TRGSignalVector& s) const
1678  {
1679 
1680  const unsigned rem = t % 16;
1681 
1682  if ((rem > 1) && (rem < 14)) {
1683  s.push_back(* _fasMap[t * 4 + 0]);
1684  s.push_back(* _fasMap[t * 4 + 1]);
1685  s.push_back(* _fasMap[t * 4 + 2]);
1686  s.push_back(* _fasMap[t * 4 + 3]);
1687  return;
1688  }
1689 
1690  //...Check hit map if there is a hit...
1691  bool hit = false;
1692  for (unsigned i = 0; i < 11; i++) {
1693  if (s[i + 1].active()) {
1694  hit = true;
1695  break;
1696  }
1697  }
1698 
1699  //...No hit case...
1700  if (! hit) {
1701  s.push_back(* _fasMap[t * 4 + 0]);
1702  s.push_back(* _fasMap[t * 4 + 1]);
1703  s.push_back(* _fasMap[t * 4 + 2]);
1704  s.push_back(* _fasMap[t * 4 + 3]);
1705  return;
1706  }
1707 
1708  const string sn = "TSF fastest timing inner";
1709  TRGDebug::enterStage(sn);
1710 
1711  //...Prepare timing signal vectors
1712  TRGSignalVector t0("t0", clockData()); // main
1713  t0.push_back(* _fasMap[t * 4 + 0]);
1714  t0.push_back(* _fasMap[t * 4 + 1]);
1715  t0.push_back(* _fasMap[t * 4 + 2]);
1716  t0.push_back(* _fasMap[t * 4 + 3]);
1717 
1718  TRGSignalVector t1("t1", clockData()); // edge
1719  TRGSignal ht0("t0 hit", clockData());
1720  TRGSignal ht1("t1 hit", clockData());
1721  if (rem == 0) {
1722  unsigned n = t / 16;
1723  if (n == 0)
1724  n = nTSF / 16 - 1;
1725  else
1726  --n;
1727 
1728  t1.push_back(* _edg3Map[n * 4 + 0]);
1729  t1.push_back(* _edg3Map[n * 4 + 1]);
1730  t1.push_back(* _edg3Map[n * 4 + 2]);
1731  t1.push_back(* _edg3Map[n * 4 + 3]);
1732 
1733  ht0 = s[0 + 1];
1734  ht0 |= s[2 + 1];
1735  ht0 |= s[4 + 1];
1736  ht0 |= s[5 + 1];
1737  ht0 |= s[8 + 1];
1738  ht0 |= s[9 + 1];
1739  ht0 |= s[12 + 1];
1740  ht0 |= s[13 + 1];
1741 
1742  ht1 = s[1 + 1];
1743  ht1 |= s[3 + 1];
1744  ht1 |= s[6 + 1];
1745  ht1 |= s[7 + 1];
1746  ht1 |= s[10 + 1];
1747  ht1 |= s[11 + 1];
1748  } else if (rem == 1) {
1749  unsigned n = t / 16;
1750  if (n == 0)
1751  n = nTSF / 16 - 1;
1752  else
1753  --n;
1754 
1755  t1.push_back(* _edg4Map[n * 4 + 0]);
1756  t1.push_back(* _edg4Map[n * 4 + 1]);
1757  t1.push_back(* _edg4Map[n * 4 + 2]);
1758  t1.push_back(* _edg4Map[n * 4 + 3]);
1759 
1760  // ht0 = s[0 + 1];
1761  ht0 = s[1 + 1];
1762  ht0 |= s[2 + 1];
1763  ht0 |= s[3 + 1];
1764  ht0 |= s[4 + 1];
1765  ht0 |= s[5 + 1];
1766  ht0 |= s[7 + 1];
1767  ht0 |= s[8 + 1];
1768  ht0 |= s[9 + 1];
1769  ht0 |= s[10 + 1];
1770  ht0 |= s[12 + 1];
1771  ht0 |= s[13 + 1];
1772 
1773  ht1 |= s[6 + 1];
1774  ht1 |= s[11 + 1];
1775  } else if (rem == 14) {
1776  unsigned n = t / 16 + 1;
1777  if (n >= nTSF / 16)
1778  n = 0;
1779 
1780  t1.push_back(* _edg1Map[n * 4 + 0]);
1781  t1.push_back(* _edg1Map[n * 4 + 1]);
1782  t1.push_back(* _edg1Map[n * 4 + 2]);
1783  t1.push_back(* _edg1Map[n * 4 + 3]);
1784 
1785  ht0 = s[0 + 1];
1786  ht0 |= s[1 + 1];
1787  ht0 |= s[2 + 1];
1788  ht0 |= s[3 + 1];
1789  ht0 |= s[4 + 1];
1790  ht0 |= s[5 + 1];
1791  ht0 |= s[6 + 1];
1792  ht0 |= s[7 + 1];
1793  ht0 |= s[8 + 1];
1794  ht0 |= s[9 + 1];
1795  ht0 |= s[10 + 1];
1796  ht0 |= s[11 + 1];
1797  ht0 |= s[12 + 1];
1798 
1799  ht1 = s[13 + 1];
1800  } else {
1801  unsigned n = t / 16 + 1;
1802  if (n >= nTSF / 16)
1803  n = 0;
1804 
1805  t1.push_back(* _edg2Map[n * 4 + 0]);
1806  t1.push_back(* _edg2Map[n * 4 + 1]);
1807  t1.push_back(* _edg2Map[n * 4 + 2]);
1808  t1.push_back(* _edg2Map[n * 4 + 3]);
1809 
1810  ht0 = s[0 + 1];
1811  ht0 |= s[1 + 1];
1812  ht0 |= s[2 + 1];
1813  ht0 |= s[3 + 1];
1814  ht0 |= s[4 + 1];
1815  ht0 |= s[6 + 1];
1816  ht0 |= s[7 + 1];
1817  ht0 |= s[8 + 1];
1818  ht0 |= s[10 + 1];
1819  ht0 |= s[11 + 1];
1820 
1821  ht1 = s[5 + 1];
1822  ht1 |= s[9 + 1];
1823  ht1 |= s[12 + 1];
1824  ht1 |= s[13 + 1];
1825  }
1826 
1827  //...State changes...
1828  vector<int> sc = t0.stateChanges();
1829  vector<int> tmp = t1.stateChanges();
1830  sc.insert(sc.end(), tmp.begin(), tmp.end());
1831  std::sort(sc.begin(), sc.end());
1832 
1833  //...Loop over state changes...
1834  TRGSignalVector tm("fastest", clockData(), 4);
1835  int last = clockData().max();
1836  for (unsigned i = 0; i < sc.size(); i++) {
1837  if (sc[i] == last)
1838  continue;
1839 
1840  int clk = sc[i];
1841 
1842  TRGState ts0 = t0.state(clk);
1843  TRGState ts1 = t1.state(clk);
1844  unsigned tm0 = unsigned(ts0);
1845  unsigned tm1 = unsigned(ts1);
1846  bool th0 = ht0.state(clk);
1847  bool th1 = ht1.state(clk);
1848 
1849  if ((! th0) && (! th1))
1850  continue;
1851  else if (th0 && th1) {
1852  if (tm1 < tm0)
1853  tm.set(ts1, clk);
1854  else
1855  tm.set(ts0, clk);
1856  } else if (th0)
1857  tm.set(ts0, clk);
1858  else if (th1)
1859  tm.set(ts1, clk);
1860 
1861  last = clk;
1862  }
1863 
1864  //...Store signals...
1865  s.push_back(tm[0]);
1866  s.push_back(tm[1]);
1867  s.push_back(tm[2]);
1868  s.push_back(tm[3]);
1869 
1870  if (TRGDebug::level()) {
1871  ht0.name("t0 hit:" + ht0.name());
1872  ht0.dump("", TRGDebug::tab() + " ");
1873  t0.dump("", TRGDebug::tab() + " ");
1874  ht1.name("t1 hit:" + ht1.name());
1875  ht1.dump("", TRGDebug::tab() + " ");
1876  t1.dump("", TRGDebug::tab() + " ");
1877  tm.dump("", TRGDebug::tab() + " ");
1878  s.dump("", TRGDebug::tab() + " ");
1879  }
1880 
1881  TRGDebug::leaveStage(sn);
1882  }
1883 
1884  vector<TRGSignalVector*>
1885  TRGCDCTrackSegmentFinder::simulateInner(const TRGSignalVector& s, unsigned)
1886  {
1887 
1888  // This is just a simple coincidence logic. Should be replaced by
1889  // real logic.
1890 
1891  //const int width = 10;
1892 
1893  //...Layer hit...
1894  //TRGSignal l0 = s[0 + 1].widen(width);
1895  //TRGSignal l1 = s[1 + 1].widen(width) | s[2 + 1].widen(width);
1896  //TRGSignal l2 = s[3 + 1].widen(width) | s[4 + 1].widen(width) |
1897  // s[5 + 1].widen(width);
1898  //TRGSignal l3 = s[6 + 1].widen(width) | s[7 + 1].widen(width) |
1899  // s[8 + 1].widen(width) | s[9 + 1].widen(width);
1900  //TRGSignal l4 = s[10 + 1].widen(width) | s[11 + 1].widen(width) |
1901  // s[12 + 1].widen(width) | s[13 + 1].widen(width);
1902 
1903  //...Layer coincidence...
1904  //TRGSignal a0 = l1 & l2 & l3 & l4;
1905  //TRGSignal a1 = l0 & l2 & l3 & l4;
1906  //TRGSignal a2 = l0 & l1 & l3 & l4;
1907  //TRGSignal a3 = l0 & l1 & l2 & l4;
1908  //TRGSignal a4 = l0 & l1 & l2 & l3;
1909 
1910  //...Final hit... (detail hit pattern changes are ignored here)
1911  //TRGSignal a = a0 | a1 | a2 | a3 | a4; TODO should it be used somewhere?
1912 
1913  //...Check timing cells...
1914  vector<int> sc = s.stateChanges();
1915  for (unsigned i = 0; i < sc.size(); i++) {
1916  int clk = sc[i];
1917  TRGState st = s.state(clk).subset(1, 3);
1918 
1919  if (st[0]) {
1920  } else {
1921  if (st[1] && st[2]) {
1922  } else if (st[1]) {
1923  } else if (st[2]) {
1924  }
1925  }
1926  }
1927 
1928  TRGSignalVector* r0 = new TRGSignalVector("tmp0", clockData(), 21);
1929  TRGSignalVector* r1 = new TRGSignalVector("tmp0", clockData(), 21);
1930  vector<TRGSignalVector*> v;
1931  v.push_back(r0);
1932  v.push_back(r1);
1933  return v;
1934  }
1935 
1936  void
1937  TRGCDCTrackSegmentFinder::addID(TRGSignalVector& s, unsigned id)
1938  {
1939  const TRGState sid(8, id);
1940  TRGSignalVector idv("TSFid", s.clock(), 8);
1941 
1942  //...Encode ID only when s is active...
1943  vector<int> sc = s.stateChanges();
1944  for (unsigned i = 0; i < sc.size(); i++) {
1945  if (s.active(sc[i]))
1946  idv.set(sid, sc[i]);
1947  }
1948 
1949  //...Merge ID vectors...
1950  s += idv;
1951  }
1952 
1953  void
1954  TRGCDCTrackSegmentFinder::simulate(void)
1955  {
1956 
1957  //...Delete old objects...
1958  for (unsigned i = 0; i < nOutput(); i++) {
1959  if (output(i))
1960  if (output(i)->signal())
1961  delete output(i)->signal();
1962  }
1963  for (unsigned i = 0; i < _tsfIn.size(); i++)
1964  delete _tsfIn[i];
1965  _tsfIn.clear();
1966  for (unsigned i = 0; i < _tsfOut.size(); i++)
1967  delete _tsfOut[i];
1968  _tsfOut.clear();
1969  for (unsigned i = 0; i < _toBeDeleted.size(); i++)
1970  delete _toBeDeleted[i];
1971  _toBeDeleted.clear();
1972 
1973  //...Clear old pointers...
1974  for (unsigned i = 0; i < 5; i++)
1975  _hitMap[i].clear();
1976  _priMap.clear();
1977  _fasMap.clear();
1978  _secMap.clear();
1979  _edg0Map.clear();
1980  _edg1Map.clear();
1981  _edg2Map.clear();
1982  _edg3Map.clear();
1983  _edg4Map.clear();
1984 
1985  if (_type == innerType)
1986  simulateInner();
1987  else
1988  simulateOuter();
1989  }
1990 
1991  void
1992  TRGCDCTrackSegmentFinder::simulate2(void)
1993  {
1994 
1995  const string sn = "TSF::simulate2 : " + name();
1996  TRGDebug::enterStage(sn);
1997 
1998  // if (_type == innerType) {
1999  // cout << "??? tmp skip" << endl;
2000  // TRGDebug::leaveStage(sn);
2001  // return;
2002  // }
2003 
2004  //...Delete old objects...
2005  for (unsigned i = 0; i < nOutput(); i++) {
2006  if (output(i))
2007  if (output(i)->signal())
2008  delete output(i)->signal();
2009  }
2010  for (unsigned i = 0; i < _tsfIn.size(); i++)
2011  delete _tsfIn[i];
2012  _tsfIn.clear();
2013  for (unsigned i = 0; i < _tsfOut.size(); i++)
2014  delete _tsfOut[i];
2015  _tsfOut.clear();
2016  for (unsigned i = 0; i < _toBeDeleted.size(); i++)
2017  delete _toBeDeleted[i];
2018  _toBeDeleted.clear();
2019 
2020  //...Clear old pointers...
2021  for (unsigned i = 0; i < 5; i++)
2022  _hitMap[i].clear();
2023  _priMap.clear();
2024  _fasMap.clear();
2025  _secMap.clear();
2026  _edg0Map.clear();
2027  _edg1Map.clear();
2028  _edg2Map.clear();
2029  _edg3Map.clear();
2030  _edg4Map.clear();
2031 
2032  //...Storage preparation...
2033  const unsigned nTSF = nInput() * 16;
2034  vector<TRGSignalVector*> trker[4];
2035  vector<int> tsfStateChanges;
2036 
2037  //...Creates hit maps...
2038  if (_type == innerType)
2039  hitMapInner();
2040  else
2041  hitMapOuter();
2042 
2043  //...Form a TSF...
2044  for (unsigned t = 0; t < nTSF; t++) {
2045 
2046  const string n = name() + "-" + TRGUtilities::itostring(t);
2047  TRGSignalVector* s = new TRGSignalVector(n, clockData());
2048  _tsfIn.push_back(s);
2049 
2050  //...Make input signals...
2051  if (_type == innerType)
2052  inputInner(t, nTSF, s);
2053  else
2054  inputOuter(t, nTSF, s);
2055 
2056  //...Simulate TSF...
2057  vector<TRGSignalVector*> result = simulateTSF(s, t);
2058  vector<TRGSignalVector*> result2 = simulateTSF2(s, t);
2059  _toBeDeleted.push_back(result[1]); // Event timing omitted now
2060 
2061  if (TRGDebug::level()) {
2062  if (result.size() != result2.size()) {
2063  cout << "TSF::simulateTSF2 has different response(size)"
2064  << endl;
2065  } else {
2066  for (unsigned i = 0; i < result.size(); i++) {
2067  if ((* result[i]) != (* result2[i]))
2068  cout << "TSF::simulateTSF2 has different response"
2069  << "(contents)" << endl;
2070  }
2071  }
2072  }
2073 
2074  TRGSignalVector* forTracker0 = result[0];
2075  TRGSignalVector* forTracker1 = new TRGSignalVector(* forTracker0);
2076 
2077  _tsfOut.push_back(forTracker0);
2078  _tsfOut.push_back(forTracker1);
2079 
2080  //...State change list...
2081  vector<int> sc = forTracker0->stateChanges();
2082  for (unsigned i = 0; i < sc.size(); i++) {
2083  bool skip = false;
2084  for (unsigned j = 0; j < tsfStateChanges.size(); j++) {
2085  if (tsfStateChanges[j] == sc[i]) {
2086  skip = true;
2087  break;
2088  }
2089  }
2090  if (! skip) tsfStateChanges.push_back(sc[i]);
2091  }
2092 
2093  //...Store it for each tracker division...
2094  const unsigned pos = t / (nTSF / 4);
2095  if (pos == 0) {
2096  addID(* forTracker0, t - 3 * (nTSF / 4));
2097  addID(* forTracker1, t);
2098  trker[3].push_back(forTracker0);
2099  trker[0].push_back(forTracker1);
2100  } else if (pos == 1) {
2101  addID(* forTracker0, t);
2102  addID(* forTracker1, t - 1 * (nTSF / 4));
2103  trker[0].push_back(forTracker0);
2104  trker[1].push_back(forTracker1);
2105  } else if (pos == 2) {
2106  addID(* forTracker0, t - 1 * (nTSF / 4));
2107  addID(* forTracker1, t - 2 * (nTSF / 4));
2108  trker[1].push_back(forTracker0);
2109  trker[2].push_back(forTracker1);
2110  } else {
2111  addID(* forTracker0, t - 2 * (nTSF / 4));
2112  addID(* forTracker1, t - 3 * (nTSF / 4));
2113  trker[2].push_back(forTracker0);
2114  trker[3].push_back(forTracker1);
2115  }
2116 
2117  if (TRGDebug::level())
2118  if (forTracker0->active())
2119  cout << TRGDebug::tab() << name() << " : TSF out="
2120  << t << endl;
2121  }
2122 
2123  //...Sort state changes...
2124  std::sort(tsfStateChanges.begin(), tsfStateChanges.end());
2125 
2126  //...Output for 2D...
2127  for (unsigned i = 0; i < 4; i++) {
2128  string n = name() + "-trker" + TRGUtilities::itostring(i);
2129  TRGSignalVector* tOut = packerForTracker(trker[i],
2130  tsfStateChanges,
2131  20);
2132  tOut->name(n);
2133  TRGSignalBundle* b = new TRGSignalBundle(n, clockData());
2134  b->push_back(tOut);
2135  output(i)->signal(b);
2136  _toBeDeleted.push_back(tOut);
2137 
2138  if (TRGCDC::getTRGCDC()->firmwareSimulationMode() & 0x4)
2139  b->dumpCOE("",
2140  TRGCDC::getTRGCDC()->firmwareSimulationStartDataClock(),
2141  TRGCDC::getTRGCDC()->firmwareSimulationStopDataClock());
2142  // b->dump();
2143  }
2144 
2145  TRGDebug::leaveStage(sn);
2146  }
2147 
2148  void
2149  TRGCDCTrackSegmentFinder::hitMapInner(void)
2150  {
2151 
2152  // dump("detail","??? ");
2153 
2154  //...Loop over mergers to create a super layer hit map...
2155  for (unsigned m = 0; m < nInput(); m++) {
2156  TRGSignalBundle* b = input(m)->signal();
2157 
2158  // b->dump("", "??? ");
2159 
2160  for (unsigned i = 0; i < 16; i++) {
2161  _secMap.push_back(& ((* b)[0][0][208 + i]));
2162  for (unsigned j = 0; j < 5; j++) {
2163  _hitMap[j].push_back(& ((* b)[0][0][j * 16 + i]));
2164 
2165  // _hitMap[j][i]->dump("", "??? " + TRGUtilities::itostring(i) + "-" + TRGUtilities::itostring(j));
2166 
2167  }
2168  for (unsigned j = 0; j < 4; j++)
2169  _priMap.push_back(& ((* b)[0][0][80 + i * 4 + j]));
2170  for (unsigned j = 0; j < 4; j++)
2171  _fasMap.push_back(& ((* b)[0][0][144 + i * 4 + j]));
2172  }
2173 
2174  for (unsigned i = 0; i < 4; i++)
2175  _edg0Map.push_back(& ((* b)[0][0][224 + i]));
2176  for (unsigned i = 0; i < 4; i++)
2177  _edg1Map.push_back(& ((* b)[0][0][224 + 4 + i]));
2178  for (unsigned i = 0; i < 4; i++)
2179  _edg2Map.push_back(& ((* b)[0][0][224 + 8 + i]));
2180  for (unsigned i = 0; i < 4; i++)
2181  _edg3Map.push_back(& ((* b)[0][0][224 + 12 + i]));
2182  for (unsigned i = 0; i < 4; i++)
2183  _edg4Map.push_back(& ((* b)[0][0][224 + 16 + i]));
2184  }
2185  }
2186 
2187  void
2188  TRGCDCTrackSegmentFinder::hitMapOuter(void)
2189  {
2190 
2191  //...Loop over mergers to create a super layer hit map...
2192  for (unsigned m = 0; m < nInput(); m++) {
2193  TRGSignalBundle* b = input(m)->signal();
2194 
2195  for (unsigned i = 0; i < 16; i++) {
2196  _secMap.push_back(& ((* b)[0][0][208 + i]));
2197  for (unsigned j = 0; j < 5; j++)
2198  _hitMap[j].push_back(& ((* b)[0][0][j * 16 + i]));
2199  for (unsigned j = 0; j < 4; j++)
2200  _priMap.push_back(& ((* b)[0][0][80 + i * 4 + j]));
2201  for (unsigned j = 0; j < 4; j++)
2202  _fasMap.push_back(& ((* b)[0][0][144 + i * 4 + j]));
2203  }
2204 
2205  for (unsigned i = 0; i < 4; i++)
2206  _edg0Map.push_back(& ((* b)[0][0][224 + i]));
2207  for (unsigned i = 0; i < 4; i++)
2208  _edg1Map.push_back(& ((* b)[0][0][224 + 4 + i]));
2209  for (unsigned i = 0; i < 4; i++)
2210  _edg2Map.push_back(& ((* b)[0][0][224 + 8 + i]));
2211  }
2212  }
2213 
2214  void
2215  TRGCDCTrackSegmentFinder::inputInner(const unsigned t,
2216  const unsigned nTSF,
2217  TRGSignalVector* s)
2218  {
2219  s->push_back(* _secMap[t]);
2220 
2221  if (t == 0) {
2222  s->push_back(* (_hitMap[0][0]));
2223  s->push_back(* (_hitMap[1][nTSF - 1]));
2224  s->push_back(* (_hitMap[1][0]));
2225  s->push_back(* (_hitMap[2][nTSF - 1]));
2226  s->push_back(* (_hitMap[2][0]));
2227  s->push_back(* (_hitMap[2][1]));
2228  s->push_back(* (_hitMap[3][nTSF - 2]));
2229  s->push_back(* (_hitMap[3][nTSF - 1]));
2230  s->push_back(* (_hitMap[3][0]));
2231  s->push_back(* (_hitMap[3][1]));
2232  s->push_back(* (_hitMap[4][nTSF - 2]));
2233  s->push_back(* (_hitMap[4][nTSF - 1]));
2234  s->push_back(* (_hitMap[4][0]));
2235  s->push_back(* (_hitMap[4][1]));
2236  s->push_back(* (_hitMap[4][2]));
2237  } else if (t == 1) {
2238  s->push_back(* (_hitMap[0][1]));
2239  s->push_back(* (_hitMap[1][0]));
2240  s->push_back(* (_hitMap[1][1]));
2241  s->push_back(* (_hitMap[2][0]));
2242  s->push_back(* (_hitMap[2][1]));
2243  s->push_back(* (_hitMap[2][2]));
2244  s->push_back(* (_hitMap[3][nTSF - 1]));
2245  s->push_back(* (_hitMap[3][0]));
2246  s->push_back(* (_hitMap[3][1]));
2247  s->push_back(* (_hitMap[3][2]));
2248  s->push_back(* (_hitMap[4][nTSF - 1]));
2249  s->push_back(* (_hitMap[4][0]));
2250  s->push_back(* (_hitMap[4][1]));
2251  s->push_back(* (_hitMap[4][2]));
2252  s->push_back(* (_hitMap[4][3]));
2253  } else if (t == (nTSF - 2)) {
2254  s->push_back(* (_hitMap[0][t]));
2255  s->push_back(* (_hitMap[1][t - 11]));
2256  s->push_back(* (_hitMap[1][t]));
2257  s->push_back(* (_hitMap[2][t - 1]));
2258  s->push_back(* (_hitMap[2][t]));
2259  s->push_back(* (_hitMap[2][0]));
2260  s->push_back(* (_hitMap[3][t - 2]));
2261  s->push_back(* (_hitMap[3][t - 1]));
2262  s->push_back(* (_hitMap[3][t]));
2263  s->push_back(* (_hitMap[3][0]));
2264  s->push_back(* (_hitMap[4][t - 2]));
2265  s->push_back(* (_hitMap[4][t - 1]));
2266  s->push_back(* (_hitMap[4][t]));
2267  s->push_back(* (_hitMap[4][t + 1]));
2268  s->push_back(* (_hitMap[4][0]));
2269  } else if (t == (nTSF - 1)) {
2270  s->push_back(* (_hitMap[0][t]));
2271  s->push_back(* (_hitMap[1][t - 11]));
2272  s->push_back(* (_hitMap[1][t]));
2273  s->push_back(* (_hitMap[2][t - 1]));
2274  s->push_back(* (_hitMap[2][t]));
2275  s->push_back(* (_hitMap[2][0]));
2276  s->push_back(* (_hitMap[3][t - 2]));
2277  s->push_back(* (_hitMap[3][t - 1]));
2278  s->push_back(* (_hitMap[3][t]));
2279  s->push_back(* (_hitMap[3][0]));
2280  s->push_back(* (_hitMap[4][t - 2]));
2281  s->push_back(* (_hitMap[4][t - 1]));
2282  s->push_back(* (_hitMap[4][t]));
2283  s->push_back(* (_hitMap[4][0]));
2284  s->push_back(* (_hitMap[4][1]));
2285  } else {
2286  // _hitMap[1][t - 1]->dump("", "??? " + TRGUtilities::itostring(1) + "-" + TRGUtilities::itostring(t - 1));
2287 
2288  TRGSignal& ts = * _hitMap[1][t - 1];
2289 
2290  s->push_back(* (_hitMap[0][t]));
2291  // s->push_back(* (_hitMap[1][t - 11])); // Why this makes SegV?
2292  s->push_back(ts);
2293  s->push_back(* (_hitMap[1][t]));
2294  s->push_back(* (_hitMap[2][t - 1]));
2295  s->push_back(* (_hitMap[2][t]));
2296  s->push_back(* (_hitMap[2][0]));
2297  s->push_back(* (_hitMap[3][t - 2]));
2298  s->push_back(* (_hitMap[3][t - 1]));
2299  s->push_back(* (_hitMap[3][t]));
2300  s->push_back(* (_hitMap[3][0]));
2301  s->push_back(* (_hitMap[4][t - 2]));
2302  s->push_back(* (_hitMap[4][t - 1]));
2303  s->push_back(* (_hitMap[4][t]));
2304  s->push_back(* (_hitMap[4][t + 1]));
2305  s->push_back(* (_hitMap[4][t + 2]));
2306  }
2307 
2308  //...Priority timing...
2309 // s->push_back(* _priMap[t * 4 + 0]);
2310 // s->push_back(* _priMap[t * 4 + 1]);
2311 // s->push_back(* _priMap[t * 4 + 2]);
2312 // s->push_back(* _priMap[t * 4 + 3]);
2313 
2314  //...Priority timing...
2315  priorityTiming(t, nTSF, * s, (* s)[1], (* s)[2], (* s)[3]);
2316 
2317  //...Fastest timing...
2318  fastestTimingInner(t, nTSF, * s);
2319 
2320  //...Clock counter is omitted...
2321  }
2322 
2323  void
2324  TRGCDCTrackSegmentFinder::inputOuter(const unsigned t,
2325  const unsigned nTSF,
2326  TRGSignalVector* s)
2327  {
2328  s->push_back(* _secMap[t]);
2329 
2330  if (t == 0) {
2331  s->push_back(* (_hitMap[0][nTSF - 1]));
2332  s->push_back(* (_hitMap[0][0]));
2333  s->push_back(* (_hitMap[0][1]));
2334  s->push_back(* (_hitMap[1][nTSF - 1]));
2335  s->push_back(* (_hitMap[1][0]));
2336  s->push_back(* (_hitMap[2][0]));
2337  s->push_back(* (_hitMap[3][nTSF - 1]));
2338  s->push_back(* (_hitMap[3][0]));
2339  s->push_back(* (_hitMap[4][nTSF - 1]));
2340  s->push_back(* (_hitMap[4][0]));
2341  s->push_back(* (_hitMap[4][1]));
2342  } else if (t == (nTSF - 1)) {
2343  s->push_back(* (_hitMap[0][nTSF - 2]));
2344  s->push_back(* (_hitMap[0][nTSF - 1]));
2345  s->push_back(* (_hitMap[0][0]));
2346  s->push_back(* (_hitMap[1][nTSF - 2]));
2347  s->push_back(* (_hitMap[1][nTSF - 1]));
2348  s->push_back(* (_hitMap[2][nTSF - 1]));
2349  s->push_back(* (_hitMap[3][nTSF - 2]));
2350  s->push_back(* (_hitMap[3][nTSF - 1]));
2351 
2352  s->push_back(* (_hitMap[4][nTSF - 2]));
2353  s->push_back(* (_hitMap[4][nTSF - 1]));
2354  s->push_back(* (_hitMap[4][0]));
2355  } else {
2356  s->push_back(* (_hitMap[0][t - 1]));
2357  s->push_back(* (_hitMap[0][t]));
2358  s->push_back(* (_hitMap[0][t + 1]));
2359  s->push_back(* (_hitMap[1][t - 1]));
2360  s->push_back(* (_hitMap[1][t]));
2361  s->push_back(* (_hitMap[2][t]));
2362  s->push_back(* (_hitMap[3][t - 1]));
2363  s->push_back(* (_hitMap[3][t]));
2364  s->push_back(* (_hitMap[4][t - 1]));
2365  s->push_back(* (_hitMap[4][t]));
2366  s->push_back(* (_hitMap[4][t + 1]));
2367  }
2368 
2369  //...Priority timing...
2370 // s->push_back(* _priMap[t * 4 + 0]);
2371 // s->push_back(* _priMap[t * 4 + 1]);
2372 // s->push_back(* _priMap[t * 4 + 2]);
2373 // s->push_back(* _priMap[t * 4 + 3]);
2374 
2375  //...Priority timing...
2376  priorityTiming(t, nTSF, * s, (* s)[6], (* s)[7], (* s)[8]);
2377 
2378  //...Fastest timing...
2379  fastestTimingOuter(t, nTSF, * s);
2380 
2381  //...Clock counter is omitted...
2382  }
2383 
2384  TRGSignalVector*
2385  TRGCDCTrackSegmentFinder::packerForTracker(vector<TRGSignalVector*>& hitList,
2386  vector<int>& cList,
2387  const unsigned maxHit)
2388  {
2389 
2390  TRGSignalVector* result =
2391  new TRGSignalVector("", (* hitList[0]).clock(), 21 * maxHit);
2392 
2393  for (unsigned ci = 0; ci < cList.size(); ci++) {
2394  unsigned cntHit = 0;
2395  for (unsigned hi = 0; hi < hitList.size(); hi++) {
2396  TRGState s = (* hitList[hi]).state(cList[ci]);
2397  if (s.active()) {
2398  if (cntHit >= maxHit) continue;
2399  for (unsigned j = 0; j < 21; j++) {
2400  if ((* hitList[hi])[j].state(cList[ci])) {
2401  (* result)[21 * (maxHit - 1) - (cntHit * 21) + j]
2402  .set(cList[ci], cList[ci] + 1);
2403  }
2404  }
2405  if (TRGDebug::level()) {
2406  TRGState t = hitList[hi]->state(cList[ci]).subset(13, 8);
2407  cout << TRGDebug::tab() << " hit found : TSF out local ID="
2408  << unsigned(t) << "(" << t << ")" << endl;
2409  }
2410 
2411  ++cntHit;
2412  // result->dump("", "??? ");
2413  }
2414  }
2415  }
2416 
2417  return result;
2418  }
2419 
2420  vector <TRGSignalVector*>
2421  TRGCDCTrackSegmentFinder::simulateTSFOld(TRGSignalVector* in, unsigned tsid)
2422  {
2423 
2424  //variables for common
2425  const string na = "TSF" + TRGUtilities::itostring(tsid) + " in " +
2426  name();
2427  TCSegment* tsi = _tsSL[tsid];
2428  vector <TRGSignalVector*> result;
2429 
2430  //variables for EvtTime & Low pT
2431  vector<bool> fTimeVect;
2432  // int tmpFTime = 0 ;
2433 
2434  //variables for Tracker & N.N
2435  vector <bool> tmpOutBool;
2436 
2437  //iwTRGSignalVector* resultT = new TRGSignalVector(na, in->clock(),22);
2438  TRGSignalVector* resultT = new TRGSignalVector(na, in->clock(), 13);
2439  TRGSignalVector* resultE = new TRGSignalVector(na, in->clock(), 10);
2440  TRGSignalVector* Hitmap = new TRGSignalVector(na + "HitMap", in->clock(), 0);
2441  TRGSignalVector pTime(na + "PriorityTime", in->clock(), 0);
2442  TRGSignalVector fTime(na + "FastestTime", in->clock(), 0);
2443  TRGSignal* pri0 = 0;
2444  TRGSignal* pri1 = 0;
2445  TRGSignal* pri2 = 0;
2446  if (_type == innerType) {
2447  for (unsigned i = 0; i < 16; i++) {
2448  Hitmap->push_back((* in)[i]);
2449  (* Hitmap)[i].widen(16);
2450  }
2451  for (unsigned i = 0; i < 4; i++) {
2452  pTime.push_back((* in)[i + 16]);
2453  fTime.push_back((* in)[i + 20]);
2454  }
2455  pri0 = & (*Hitmap)[1];
2456  pri1 = & (*Hitmap)[2];
2457  pri2 = & (*Hitmap)[3];
2458  } else {
2459  for (unsigned i = 0; i < 12; i++) {
2460  Hitmap->push_back((* in)[i]);
2461  (* Hitmap)[i].widen(16);
2462  }
2463  for (unsigned i = 0; i < 4; i++) {
2464  pTime.push_back((* in)[i + 12]);
2465  fTime.push_back((* in)[i + 16]);
2466  }
2467  pri0 = & (*Hitmap)[6];
2468  pri1 = & (*Hitmap)[7];
2469  pri2 = & (*Hitmap)[8];
2470  }
2471 
2472  //...Clock counter...
2473  const TRGSignalVector& cc = in->clock().clockCounter();
2474  for (unsigned i = 0; i < 5; i++) {
2475  pTime.push_back(cc[i]);
2476  fTime.push_back(cc[i]);
2477  }
2478 
2479  vector <int> changeTime = Hitmap->stateChanges();
2480 
2481  int* LUTValue = new int[changeTime.size()];
2482  if (changeTime.size()) {
2483 
2484  const string fn = "TSF::simulateTSF:tsid=" + to_string(tsid);
2485  TRGDebug::enterStage(fn);
2486 
2487  int hitPosition = 0;
2488  bool fTimeBool[10];
2489  int tmpPTime = 0 ;
2490  int tmpCTime = 0 ;
2491  int tmpOutInt;
2492  fTime.state(changeTime[0]).copy2bool(fTimeBool);
2493  fTimeBool[9] = true;
2494  fTimeVect.insert(fTimeVect.begin(), fTimeBool, fTimeBool + 10);
2495  //tmpFTime = mkint(fTime.state(changeTime[0]));
2496  bool eOUT = true;
2497  for (unsigned i = 0; i < changeTime.size(); i++) {
2498  int ct = changeTime[i];
2499 
2500 // LUTValue[i] = tsi->LUT()->getValue(mkint(Hitmap->state(ct)));
2501  LUTValue[i] = tsi->LUT()->getValue(unsigned(Hitmap->state(ct)));
2502 
2504  if ((LUTValue[i]) && (eOUT)) {
2505  resultE->set(fTimeVect, ct);
2506  eOUT = false;
2507  }
2508 
2509  bool priority1rise = pri0->riseEdge(ct);
2510  bool priority2rise = pri1->riseEdge(ct) or pri2->riseEdge(ct);
2511 
2513  //ready for output
2514  if (priority1rise) {
2515  hitPosition = 3;
2516 // tmpPTime = mkint(pTime.state(ct));
2517  tmpPTime = unsigned(pTime.state(ct));
2518  tmpCTime = ct;
2519  } else if (priority2rise) {
2520  if (!hitPosition) {
2521 // tmpPTime = mkint(pTime.state(ct));
2522  tmpPTime = unsigned(pTime.state(ct));
2523  tmpCTime = ct;
2524  if ((*Hitmap)[0].state(ct)) hitPosition = 2;
2525  else hitPosition = 1;
2526  }
2527  }
2528 
2529  // output selection
2530  if ((hitPosition) && (LUTValue[i]) && ((ct - tmpCTime) < 16)) {
2531  tmpOutInt = tmpPTime * pow(2, 4) +
2532  LUTValue[i] * pow(2, 2) + hitPosition;
2533  tmpOutBool = mkbool(tmpOutInt, 13); // ID removed : iw
2534 
2535  if (hitPosition == 3) {
2536  if (priority1rise) {
2537  resultT->set(tmpOutBool, ct);
2538  } else {
2539  if ((LUTValue[i] == 1) | (LUTValue[i] == 2)) {
2540  if (!((LUTValue[i - 1] == 1) |
2541  (LUTValue[i - 1] == 2)))
2542  resultT->set(tmpOutBool, ct);
2543  } else {
2544  if (!(LUTValue[i - 1])) resultT->set(tmpOutBool, ct);
2545  }
2546  }
2547  } else {
2548  if (priority2rise) resultT->set(tmpOutBool, ct);
2549  else {
2550  if ((LUTValue[i] == 1) | (LUTValue[i] == 2)) {
2551  if (!((LUTValue[i - 1] == 1) |
2552  (LUTValue[i - 1] == 2)))
2553  resultT->set(tmpOutBool, ct);
2554  } else {
2555  if (!(LUTValue[i])) resultT->set(tmpOutBool, ct);
2556  }
2557  }
2558  }
2559  }
2560 
2561  if (TRGDebug::level() > 1) {
2562  cout << TRGDebug::tab() << "clk=" << ct
2563  << ", pattern=" << Hitmap->state(ct)
2564  << ", LUT=" << LUTValue[i]
2565  << ", pri=" << priority1rise
2566  << ", sec=" << priority2rise
2567  << ", hitPos=" << hitPosition
2568  << ", (clk-tmpCTime)=" << (ct - tmpCTime)
2569  << endl;
2570  }
2571  }
2572 
2573  TRGDebug::leaveStage(fn);
2574  }
2575 
2576  result.push_back(resultT);
2577  result.push_back(resultE);
2578 
2579  // cppcheck-suppress uninitdata
2580  delete[] LUTValue;
2581  delete Hitmap;
2582 
2583  return result;
2584  }
2585 
2586  vector <TRGSignalVector*>
2587  TRGCDCTrackSegmentFinder::simulateTSF(TRGSignalVector* in, unsigned tsid)
2588  {
2589 
2590  //variables for common
2591  const string na = "TSF" + TRGUtilities::itostring(tsid) + " in " +
2592  name();
2593  TCSegment* tsi = _tsSL[tsid];
2594  vector <TRGSignalVector*> result;
2595 
2596  //variables for EvtTime & Low pT
2597  vector<bool> fTimeVect;
2598  // int tmpFTime = 0 ;
2599 
2600  //variables for Tracker & N.N
2601  vector <bool> tmpOutBool;
2602 
2603  TRGSignalVector* resultT = new TRGSignalVector(na, in->clock(), 13);
2604  TRGSignalVector* resultE = new TRGSignalVector(na, in->clock(), 10);
2605  TRGSignalVector* Hitmap = new TRGSignalVector(na + "HitMap",
2606  in->clock(),
2607  0);
2608  TRGSignalVector pTime(na + "PriorityTime", in->clock(), 0);
2609  TRGSignalVector fTime(na + "FastestTime", in->clock(), 0);
2610  TRGSignal* pri0 = 0;
2611  TRGSignal* pri1 = 0;
2612  TRGSignal* pri2 = 0;
2613  if (_type == innerType) {
2614  for (unsigned i = 0; i < 16; i++) {
2615  Hitmap->push_back((* in)[i]);
2616  (* Hitmap)[i].widen(16);
2617  }
2618  for (unsigned i = 0; i < 4; i++) {
2619  pTime.push_back((* in)[i + 16]);
2620  fTime.push_back((* in)[i + 20]);
2621  }
2622  pri0 = & (*Hitmap)[1];
2623  pri1 = & (*Hitmap)[2];
2624  pri2 = & (*Hitmap)[3];
2625  } else {
2626  for (unsigned i = 0; i < 12; i++) {
2627  Hitmap->push_back((* in)[i]);
2628  (* Hitmap)[i].widen(16);
2629  }
2630  for (unsigned i = 0; i < 4; i++) {
2631  pTime.push_back((* in)[i + 12]);
2632  fTime.push_back((* in)[i + 16]);
2633  }
2634  pri0 = & (*Hitmap)[6];
2635  pri1 = & (*Hitmap)[7];
2636  pri2 = & (*Hitmap)[8];
2637  }
2638 
2639  //...Clock counter...
2640  const TRGSignalVector& cc = in->clock().clockCounter();
2641  for (unsigned i = 0; i < 5; i++) {
2642  pTime.push_back(cc[i]);
2643  fTime.push_back(cc[i]);
2644  }
2645 
2646  vector<int> changeTime = Hitmap->stateChanges();
2647 
2648  int* LUTValue = new int[changeTime.size()];
2649  //int oldLUT = 0;
2650  int lastFastHit = in->clock().min();
2651  if (changeTime.size()) {
2652 
2653  int oldLUT = 0; // 2019/07/31 by ytlai, moved to here to reduce scope
2654  const string fn = "TSF::simulateTSF:tsid=" + to_string(tsid);
2655  TRGDebug::enterStage(fn);
2656 
2657  int hitPosition = 0;
2658  int tmpPTime = 0 ;
2659  int tmpCTime = 0 ;
2660  int tmpOutInt;
2661  //TRGState ft(10); //TODO variable not used at all
2662  bool eOut = false;
2663  for (unsigned i = 0; i < changeTime.size(); i++) {
2664 
2665  int ct = changeTime[i];
2666  TRGState st = Hitmap->state(ct);
2667 
2668  //...LUT...
2669  // 0:no hit, 1:right, 2:left, 3:LR unknown hit
2670  LUTValue[i] = tsi->LUT()->getValue(unsigned(st));
2671 
2672  //...Any wire hit (for the fastest timing)...
2673  bool active = st.active();
2674  if (active) {
2675  const int timeCounter = ct - lastFastHit;
2676 
2677  //...Record as the fastest timing hit...
2678  if (timeCounter > 15) {
2679  lastFastHit = ct;
2680  TRGState ftnow = fTime.state(ct);
2681  ftnow += TRGState(1, 1);
2682  //ft = ftnow;
2683  eOut = true;
2684  }
2685  }
2686 
2688  if ((LUTValue[i])) {
2689  if (eOut) {
2690  resultE->set(fTimeVect, ct);
2691  eOut = false;
2692  }
2693  }
2694 
2695  bool priority1rise = pri0->riseEdge(ct);
2696  bool priority2rise = pri1->riseEdge(ct) or pri2->riseEdge(ct);
2697 
2699  //ready for output
2700  if (priority1rise) {
2701  hitPosition = 3;
2702  tmpPTime = unsigned(pTime.state(ct));
2703  tmpCTime = ct;
2704  } else if (priority2rise) {
2705  if (!hitPosition) {
2706  tmpPTime = unsigned(pTime.state(ct));
2707  tmpCTime = ct;
2708  if ((*Hitmap)[0].state(ct)) hitPosition = 2;
2709  else hitPosition = 1;
2710  }
2711  }
2712 
2713  // output selection
2714  if ((ct - tmpCTime) < 16) {
2715  if ((hitPosition) && (LUTValue[i])) {
2716  //iw tmpOutInt = tsid * pow(2, 13) + tmpPTime * pow(2, 4) +
2717  //iw LUTValue[i] * pow(2,2) + hitPosition;
2718  tmpOutInt = tmpPTime * pow(2, 4) +
2719  LUTValue[i] * pow(2, 2) + hitPosition;
2720  tmpOutBool = mkbool(tmpOutInt, 13); // ID removed : iw
2721 
2722  if (hitPosition == 3) {
2723  if (priority1rise) {
2724  resultT->set(tmpOutBool, ct);
2725  oldLUT = LUTValue[i];
2726  } else {
2727  if ((LUTValue[i] == 1) | (LUTValue[i] == 2)) {
2728  if (!((oldLUT == 1) |
2729  (oldLUT == 2))) {
2730  resultT->set(tmpOutBool, ct);
2731  oldLUT = LUTValue[i];
2732  }
2733  } else {
2734  if (!(LUTValue[i - 1])) {
2735  resultT->set(tmpOutBool, ct);
2736  oldLUT = LUTValue[i];
2737  }
2738  }
2739  }
2740  } else {
2741  if ((priority2rise) && (!oldLUT)) resultT->set(tmpOutBool, ct);
2742  else {
2743  if ((LUTValue[i] == 1) | (LUTValue[i] == 2)) {
2744  if (!((oldLUT == 1) |
2745  (oldLUT == 2))) {
2746  resultT->set(tmpOutBool, ct);
2747  oldLUT = LUTValue[i];
2748  }
2749  } else {
2750  if (!(oldLUT)) {
2751  resultT->set(tmpOutBool, ct);
2752  oldLUT = LUTValue[i];
2753  }
2754  }
2755  }
2756  }
2757  }
2758  } else {
2759  oldLUT = 0;
2760  hitPosition = 0;
2761  }
2762 
2763  if (TRGDebug::level() > 1) {
2764  cout << TRGDebug::tab() << "clk=" << ct
2765  << ", pattern=" << st
2766  << ", LUT=" << LUTValue[i]
2767  << ", pri=" << priority1rise
2768  << ", sec=" << priority2rise
2769  << ", hitPos=" << hitPosition
2770  << ", (clk-tmpCTime)=" << (ct - tmpCTime)
2771  << endl;
2772  }
2773  }
2774  TRGDebug::leaveStage(fn);
2775  }
2776 
2777  if (TRGDebug::level()) {
2778  vector<int> sc = resultT->stateChanges();
2779  vector<unsigned> lv;
2780  for (unsigned i = 0; i < sc.size(); i++) {
2781  //...LUT value...
2782  unsigned l = unsigned(resultT->state(sc[i]).subset(2, 2));
2783  lv.push_back(l);
2784  cout << "clk=" << sc[i] << " LUT output[" << i << "]=" << l << endl;
2785  }
2786  bool found1or2 = false;
2787  bool found3 = false;
2788  for (unsigned i = 0; i < lv.size(); i++) {
2789  if (found1or2 && (lv[i] > 0))
2790  cout << "!!! simulateTSF something wrong(found1or2)" << endl;
2791  if (found3 && (lv[i] == 3))
2792  cout << "!!! simulateTSF something wrong(found3)" << endl;
2793 
2794  if ((lv[i] == 1) or (lv[i] == 2))
2795  found1or2 = true;
2796  else if (lv[i] == 3)
2797  found3 = true;
2798  }
2799  }
2800 
2801  result.push_back(resultT);
2802  result.push_back(resultE);
2803 
2804  // cppcheck-suppress uninitdata
2805  delete[] LUTValue;
2806  delete Hitmap;
2807 
2808  return result;
2809  }
2810  vector <TRGSignalVector*>
2811  TRGCDCTrackSegmentFinder::simulateTSF2(TRGSignalVector* in, unsigned tsid)
2812  {
2813 
2814  //variables for common
2815  const string na = "TSF" + TRGUtilities::itostring(tsid) + " in " +
2816  name();
2817  TCSegment* tsi = _tsSL[tsid];
2818  vector <TRGSignalVector*> result;
2819 
2820  //variables for EvtTime & Low pT
2821  vector<bool> fTimeVect;
2822  // int tmpFTime = 0 ;
2823 
2824  //variables for Tracker & N.N
2825  TRGSignalVector* resultT = new TRGSignalVector(na, in->clock(), 13);
2826  TRGSignalVector* resultE = new TRGSignalVector(na, in->clock(), 10);
2827  TRGSignalVector* Hitmap = new TRGSignalVector(na + "HitMap",
2828  in->clock(),
2829  0);
2830  TRGSignalVector pTime(na + "PriorityTime", in->clock(), 0);
2831  TRGSignalVector fTime(na + "FastestTime", in->clock(), 0);
2832  TRGSignal* pri0 = 0;
2833  TRGSignal* pri1 = 0;
2834  TRGSignal* pri2 = 0;
2835 
2836  if (_type == innerType) {
2837  for (unsigned i = 0; i < 16; i++) {
2838  Hitmap->push_back((* in)[i]);
2839  (* Hitmap)[i].widen(16);
2840  }
2841  for (unsigned i = 0; i < 4; i++) {
2842  pTime.push_back((* in)[i + 16]);
2843  fTime.push_back((* in)[i + 20]);
2844  }
2845  pri0 = & (*Hitmap)[1];
2846  pri1 = & (*Hitmap)[2];
2847  pri2 = & (*Hitmap)[3];
2848  } else {
2849  for (unsigned i = 0; i < 12; i++) {
2850  Hitmap->push_back((* in)[i]);
2851  (* Hitmap)[i].widen(16);
2852  }
2853  for (unsigned i = 0; i < 4; i++) {
2854  pTime.push_back((* in)[i + 12]);
2855  fTime.push_back((* in)[i + 16]);
2856  }
2857  pri0 = & (*Hitmap)[6];
2858  pri1 = & (*Hitmap)[7];
2859  pri2 = & (*Hitmap)[8];
2860  }
2861 
2862  //...Clock counter...
2863  const TRGSignalVector& cc = in->clock().clockCounter();
2864  for (unsigned i = 0; i < 5; i++) {
2865  pTime.push_back(cc[i]);
2866  fTime.push_back(cc[i]);
2867  }
2868 
2869  vector<int> changeTime = Hitmap->stateChanges();
2870 
2871  vector<unsigned> luts;
2872  int lastFastHit = in->clock().min();
2873 
2874  if (changeTime.size()) {
2875 
2876  const string fn = "TSF::simulateTSF2:tsid=" + to_string(tsid);
2877  TRGDebug::enterStage(fn);
2878 
2879  //TRGState ft(10); //not used, should be?
2880  bool eOut = false;
2881 
2882  unsigned stateHitPos = 0; // 0:wait, 1:active1st, 2:active2nd
2883  unsigned hitPos = 0;
2884  unsigned timing = 0;
2885  int timePosFound = 0;
2886  int counterPos = 0;
2887 
2888  unsigned stateLR = 0; // 0:wait, 1:LUT12, 2:LUT3
2889 // int timeLRFound = 0;
2890 // int counterLR = 0;
2891 
2892  for (unsigned i = 0; i < changeTime.size(); i++) {
2893 
2894  const int ct = changeTime[i];
2895  const TRGState st = Hitmap->state(ct);
2896 
2897  //...LUT...
2898  // 0:no hit, 1:right, 2:left, 3:LR unknown hit
2899  const unsigned lut = tsi->LUT()->getValue(unsigned(st));
2900  luts.push_back(lut);
2901 
2902  //...Any wire hit (for the fastest timing)...
2903  const bool active = st.active();
2904  if (active) {
2905  const int timeCounter = ct - lastFastHit;
2906 
2907  //...Record as the fastest timing hit...
2908  if (timeCounter > 15) {
2909  lastFastHit = ct;
2910  TRGState ftnow = fTime.state(ct);
2911  ftnow += TRGState(1, 1);
2912  //ft = ftnow;
2913  eOut = true;
2914  }
2915  }
2916 
2917  // output for EvtTime & Low pT tracker module
2918  if (lut) {
2919  if (eOut) {
2920  resultE->set(fTimeVect, ct);
2921  eOut = false;
2922  }
2923  }
2924 
2925  //...Hit position...
2926  const bool hit1st = pri0->riseEdge(ct);
2927  const bool hit2nd0 = pri1->riseEdge(ct);
2928  const bool hit2nd1 = pri2->riseEdge(ct);
2929  const bool hit2nd = hit2nd0 || hit2nd1;
2930 
2931  //...Hit position state machine...
2932 
2933  //...State : wait...
2934  if (stateHitPos == 0) {
2935  hitPos = 0;
2936  timing = 0;
2937  timePosFound = 0;
2938  if (hit1st) {
2939  hitPos = 3;
2940  timing = unsigned(pTime.state(ct));
2941  timePosFound = ct;
2942  stateHitPos = 1;
2943  } else if (hit2nd) {
2944  if ((* Hitmap)[0].state(ct))
2945  hitPos = 2;
2946  else
2947  hitPos = 1;
2948  timing = unsigned(pTime.state(ct));
2949  timePosFound = ct;
2950  stateHitPos = 2;
2951  }
2952  }
2953 
2954  //...State : active1st...
2955  else if (stateHitPos == 1) {
2956  counterPos = ct - timePosFound;
2957  if (counterPos > 15)
2958  stateHitPos = 0;
2959  }
2960 
2961  //...State : active2nd...
2962  else if (stateHitPos == 2) {
2963  if (hit1st) {
2964  hitPos = 3;
2965  timing = unsigned(pTime.state(ct));
2966  timePosFound = ct;
2967  stateHitPos = 1;
2968  } else {
2969  counterPos = ct - timePosFound;
2970  if (counterPos > 15)
2971  stateHitPos = 0;
2972  }
2973  }
2974 
2975  //...State unknown...
2976  else {
2977  cout << "TSF::simulateTSF2 !!! strange state in hit position"
2978  << endl;
2979  }
2980 
2981  //...L/R decision state machine...
2982 
2983  //...State : wait...
2984  if (stateLR == 0) {
2985  if (lut != 0) {
2986  const unsigned val = (timing << 4) | (lut << 2) | hitPos;
2987  const TRGState output(13, val);
2988  resultT->set(output, ct);
2989 // timeLRFound = ct;
2990  if ((lut == 1) || (lut == 2))
2991  stateLR = 1;
2992  else
2993  stateLR = 2;
2994 
2995  if (TRGDebug::level()) {
2996  if (((lut == 1) || (lut == 2)) && (stateHitPos == 0))
2997  cout << TRGDebug::tab()
2998  << "!!! state machines incosistent" << endl;
2999  }
3000  }
3001  }
3002 
3003  //...State : LUT12...
3004  else if (stateLR == 1) {
3005 // counterLR = ct - timeLRFound;
3006  if (counterPos > 15)
3007  stateLR = 0;
3008  }
3009 
3010  //...State : LUT3...
3011  else if (stateLR == 2) {
3012  if (lut) {
3013  if ((lut == 1) || (lut == 2)) {
3014  const unsigned val = (timing << 4) | (lut << 2) |
3015  hitPos;
3016  const TRGState output(13, val);
3017  resultT->set(output, ct);
3018 // timeLRFound = ct;
3019  stateLR = 1;
3020  } else {
3021 // counterLR = ct - timeLRFound;
3022  if (counterPos > 15)
3023  stateLR = 0;
3024  }
3025  }
3026  }
3027 
3028  //...State unknown...
3029  else {
3030  cout << "TSF::simulateTSF2 !!! strange state in L/R decision"
3031  << endl;
3032  }
3033 
3034  if (TRGDebug::level() > 1) {
3035  cout << TRGDebug::tab() << "clk=" << ct
3036  << ", ptn=" << st
3037  << ", LUT=" << lut
3038  << ", pri=" << hit1st
3039  << ", sec=" << hit2nd
3040  << ", hPos=" << hitPos
3041  << ", ctr=" << counterPos // << "," << counterLR
3042  << ", states=" << stateHitPos << "," << stateLR
3043  << endl;
3044  }
3045  }
3046  TRGDebug::leaveStage(fn);
3047  }
3048 
3049  if (TRGDebug::level()) {
3050  vector<int> sc = resultT->stateChanges();
3051  vector<unsigned> lv;
3052  for (unsigned i = 0; i < sc.size(); i++) {
3053  //...LUT value...
3054  unsigned l = unsigned(resultT->state(sc[i]).subset(2, 2));
3055  lv.push_back(l);
3056  cout << "clk=" << sc[i] << " LUT output[" << i << "]=" << l << endl;
3057  }
3058  }
3059 
3060  result.push_back(resultT);
3061  result.push_back(resultE);
3062 
3063  delete Hitmap;
3064 
3065  return result;
3066  }
3067 
3069 } // namespace Belle2
Belle2::B2Vector3::Phi
DataType Phi() const
The azimuth angle.
Definition: B2Vector3.h:151
Belle2::MCParticle
A Class to store the Monte Carlo particle information.
Definition: MCParticle.h:32
Belle2::RelationArray
Low-level class to create/modify relations between StoreArrays.
Definition: RelationArray.h:62
Belle2::RelationArray::getEntries
int getEntries() const
Get the number of elements.
Definition: RelationArray.h:242
Belle2::StoreArray
Accessor to arrays stored in the data store.
Definition: StoreArray.h:113
Belle2::StoreArray::getEntries
int getEntries() const
Get the number of objects in the array.
Definition: StoreArray.h:216
Belle2::TRGBoard
A class to represent a trigger board.
Definition: Board.h:25
Belle2::TRGCDCCellHit
A class to represent a wire hit in CDC.
Definition: CellHit.h:74
Belle2::TRGCDCCell
A class to represent a wire in CDC.
Definition: Cell.h:40
Belle2::TRGCDCCell::signal
virtual const TRGSignal & signal(void) const =0
returns trigger output. Null will returned if no signal.
Belle2::TRGCDCMerger
A class to represent a CDC merger board.
Definition: Merger.h:36
Belle2::TRGCDCTrackSegmentFinder::_priMap
std::vector< TRGSignal * > _priMap
Internal data of the priority hit timing.
Definition: TrackSegmentFinder.h:274
Belle2::TRGCDCTrackSegmentFinder::m_hitPatternInformation
TClonesArray * m_hitPatternInformation
Stores hitpattern information.
Definition: TrackSegmentFinder.h:101
Belle2::TRGCDCTrackSegmentFinder::_edg4Map
std::vector< TRGSignal * > _edg4Map
Internal data of the edge timing information.
Definition: TrackSegmentFinder.h:295
Belle2::TRGCDCTrackSegmentFinder::_edg0Map
std::vector< TRGSignal * > _edg0Map
Internal data of the edge timing information.
Definition: TrackSegmentFinder.h:283
Belle2::TRGCDCTrackSegmentFinder::m_treeInputTSF
TTree * m_treeInputTSF
ROOT TTree for input.
Definition: TrackSegmentFinder.h:95
Belle2::TRGCDCTrackSegmentFinder::m_makeRootFile
bool m_makeRootFile
make ROOT file or not
Definition: TrackSegmentFinder.h:117
Belle2::TRGCDCTrackSegmentFinder::m_particleEfficiency
TClonesArray * m_particleEfficiency
[Efficiency, Pt, # MC TS] Efficiency = -1 means that # MC TS is 0.
Definition: TrackSegmentFinder.h:107
Belle2::TRGCDCTrackSegmentFinder::_toBeDeleted
std::vector< TRGSignalVector * > _toBeDeleted
One time info. to be deleted in next event;.
Definition: TrackSegmentFinder.h:268
Belle2::TRGCDCTrackSegmentFinder::_edg1Map
std::vector< TRGSignal * > _edg1Map
Internal data of the edge timing information.
Definition: TrackSegmentFinder.h:286
Belle2::TRGCDCTrackSegmentFinder::m_nnPatternInformation
TClonesArray * m_nnPatternInformation
[superlayer id, lrDriftTime, timeWire0, timeWire1, ..., ...]
Definition: TrackSegmentFinder.h:114
Belle2::TRGCDCTrackSegmentFinder::m_treeOutputTSF
TTree * m_treeOutputTSF
ROOT TTree for output.
Definition: TrackSegmentFinder.h:104
Belle2::TRGCDCTrackSegmentFinder::_edg2Map
std::vector< TRGSignal * > _edg2Map
Internal data of the edge timing information.
Definition: TrackSegmentFinder.h:289
Belle2::TRGCDCTrackSegmentFinder::m_tsInformation
TClonesArray * m_tsInformation
[SuperLayer Id, Wire Id, Priority Timing]
Definition: TrackSegmentFinder.h:109
Belle2::TRGCDCTrackSegmentFinder::_tsfOut
std::vector< TRGSignalVector * > _tsfOut
TSF response storeage.
Definition: TrackSegmentFinder.h:265
Belle2::TRGCDCTrackSegmentFinder::_hitMap
std::vector< TRGSignal * > _hitMap[5]
Internal data of wire hit map.
Definition: TrackSegmentFinder.h:271
Belle2::TRGCDCTrackSegmentFinder::m_treeNNTSF
TTree * m_treeNNTSF
ROOT Tree for NNTSF.
Definition: TrackSegmentFinder.h:112
Belle2::TRGCDCTrackSegmentFinder::_secMap
std::vector< TRGSignal * > _secMap
Internal data of the priority cell hit position flag.
Definition: TrackSegmentFinder.h:280
Belle2::TRGCDCTrackSegmentFinder::_tsSL
std::vector< TCSegment * > _tsSL
list of TSF
Definition: TrackSegmentFinder.h:259
Belle2::TRGCDCTrackSegmentFinder::_tsfIn
std::vector< TRGSignalVector * > _tsfIn
TSF input storage.
Definition: TrackSegmentFinder.h:262
Belle2::TRGCDCTrackSegmentFinder::_edg3Map
std::vector< TRGSignal * > _edg3Map
Internal data of the edge timing information.
Definition: TrackSegmentFinder.h:292
Belle2::TRGCDCTrackSegmentFinder::m_logicLUTFlag
bool m_logicLUTFlag
0 is Logic. 1 is LUT.
Definition: TrackSegmentFinder.h:87
Belle2::TRGCDCTrackSegmentFinder::_fasMap
std::vector< TRGSignal * > _fasMap
Internal data of the fastest hit timing.
Definition: TrackSegmentFinder.h:277
Belle2::TRGCDCTrackSegmentFinder::boardType
boardType
enum of boardType of TrackSegmentFinder
Definition: TrackSegmentFinder.h:46
Belle2::TRGCDCWireHit
A class to represent a wire hit in CDC.
Definition: WireHit.h:33
Belle2::TRGCDCWire
A class to represent a wire in CDC.
Definition: Wire.h:56
Belle2::TRGCDC
The instance of TRGCDC is a singleton.
Definition: TRGCDC.h:69
Belle2::TRGClock
A class to represent a digitized signal. Unit is nano second.
Definition: Clock.h:38
Belle2::TRGSignalBundle
A class to represent a bundle of SignalVectors.
Definition: SignalBundle.h:26
Belle2::TRGSignalVector
A class to represent a bundle of digitized signals.
Definition: SignalVector.h:26
Belle2::TRGSignal
A class to represent a digitized signal. Unit is nano second.
Definition: Signal.h:23
Belle2::TRGState
A class to represent a state of multi bits.
Definition: State.h:24
Belle2::TRGCDCTrackSegmentFinder::simulateTSFOld
vector< TRGSignalVector * > simulateTSFOld(TRGSignalVector *in, unsigned id)
Simulate TSF response (unified version)
Definition: TrackSegmentFinder.cc:2421
Belle2::TRGState::size
unsigned size(void) const
returns bit size.
Definition: State.h:149
Belle2::TRGDebug::tab
static std::string tab(void)
returns tab spaces.
Definition: Debug.cc:47
Belle2::TRGCDCTrackSegmentFinder::~TRGCDCTrackSegmentFinder
~TRGCDCTrackSegmentFinder()
Destructor.
Definition: TrackSegmentFinder.cc:108
Belle2::TRGSignalVector::set
const TRGSignalVector & set(const TRGState &, int clockPosition)
sets state at given clock.
Definition: SignalVector.cc:181
Belle2::TRGCDC::getTRGCDC
static TRGCDC * getTRGCDC(void)
returns TRGCDC object.
Definition: TRGCDC.cc:192
Belle2::TRGCDCTrackSegmentFinder::fastestTimingOuter
void fastestTimingOuter(unsigned tsfID, const unsigned nTSF, TRGSignalVector &s) const
Make TRGSignals for the fastest timing bits.
Definition: TrackSegmentFinder.cc:1289
Belle2::TRGCDCCellHit::mcLR
int mcLR(void) const
returns mc left/right information
Definition: CellHit.h:245
Belle2::TRGBoard::input
const TRGChannel * input(unsigned i) const
returns input channel i.
Definition: Board.h:158
Belle2::TRGCDCTrackSegmentFinder::simulateOuter
void simulateOuter(void)
Firmware simulation for the outers. yi.
Definition: TrackSegmentFinder.cc:872
Belle2::TRGState::subset
TRGState subset(unsigned i, unsigned n) const
returns subset from i with n bits.
Definition: State.cc:356
Belle2::TRGState::copy2bool
void copy2bool(bool *array) const
returns a filled array.
Definition: State.h:255
Belle2::TRGSignalVector::state
TRGState state(int clockPosition) const
returns state at given clock position.
Definition: SignalVector.cc:158
Belle2::TRGCDCTrackSegmentFinder::packerOuterTracker
TRGSignalVector * packerOuterTracker(vector< TRGSignalVector * > &, vector< int > &, const unsigned)
Packing output for tracker.
Definition: TrackSegmentFinder.cc:320
Belle2::TRGSignalVector::clock
const TRGClock & clock(void) const
returns clock.
Definition: SignalVector.h:125
Belle2::TRGCDCTrackSegmentFinder::priorityTiming
void priorityTiming(unsigned tsfID, const unsigned nTSF, TRGSignalVector &s, const TRGSignal &center, const TRGSignal &right, const TRGSignal &left) const
Make TRGSignals for the priority timing bits.
Definition: TrackSegmentFinder.cc:1174
Belle2::TRGCDCTrackSegmentFinder::terminate
void terminate(void)
terminate
Definition: TrackSegmentFinder.cc:125
Belle2::TRGBoard::clockData
const TRGClock & clockData(void) const
returns data clock.
Definition: Board.h:123
Belle2::TRGCDCTrackSegmentFinder::packerForTracker
TRGSignalVector * packerForTracker(vector< TRGSignalVector * > &, vector< int > &, const unsigned)
Output packer for tracker.
Definition: TrackSegmentFinder.cc:2385
Belle2::TRGSignal::riseEdge
bool riseEdge(int clockPosition) const
returns true if signal is active and rising edge in give clock position.
Definition: Signal.h:304
Belle2::TRGCDCCellHit::hit
const CDCHit * hit(void) const
Access to CDCHit.
Definition: CellHit.cc:133
Belle2::TRGCDCTrackSegmentFinder::push_back
void push_back(const TRGCDCMerger *)
push back the Mergers of this TSF
Definition: TrackSegmentFinder.cc:865
Belle2::TRGCDCTrackSegmentFinder::findTSHit
vector< TRGSignalVector * > findTSHit(TRGSignalVector *eachInput, int)
Use LUT for find TSHit.
Definition: TrackSegmentFinder.cc:199
Belle2::TRGCDCTrackSegmentFinder::simulateTSF
vector< TRGSignalVector * > simulateTSF(TRGSignalVector *in, unsigned id)
Simulate TSF response (unified version, 2016/07/12)
Definition: TrackSegmentFinder.cc:2587
Belle2::TRGCDCTrackSegmentFinder::inputOuter
void inputOuter(const unsigned id, const unsigned nTSF, TRGSignalVector *s)
Creates input signals to TSF for the outer.
Definition: TrackSegmentFinder.cc:2324
Belle2::TRGBoard::name
const std::string & name(void) const
returns name.
Definition: Board.h:109
Belle2::TRGSignalVector::stateChanges
std::vector< int > stateChanges(void) const
returns a list of clock position of state change.
Definition: SignalVector.cc:130
Belle2::TRGDebug::enterStage
static void enterStage(const std::string &stageName)
Declare that you enter new stage.
Definition: Debug.cc:24
Belle2::TRGBoard::nInput
unsigned nInput(void) const
returns input channels.
Definition: Board.h:172
Belle2::TRGClock::min
int min(void) const
returns min. clock point.
Definition: Clock.h:194
Belle2::TRGCDCTrackSegmentFinder::TRGCDCTrackSegmentFinder
TRGCDCTrackSegmentFinder(const TRGCDC &, bool makeRootFile, bool logicLUTFlag)
Constructor.
Definition: TrackSegmentFinder.cc:49
Belle2::TRGCDCCell::hit
const TRGCDCCellHit * hit(void) const
returns a pointer to a TRGCDCCellHit.
Definition: Cell.h:362
Belle2::TRGCDCWire::signal
const TRGSignal & signal(void) const override
returns an input to the trigger. This is sync'ed to 1GHz clock.
Definition: Wire.h:226
Belle2::TRGUtilities::itostring
static std::string itostring(int i)
converts int to string. (Use boost::lexical_cast)
Definition: Utilities.cc:50
Belle2::TRGCDCTrackSegmentFinder::simulateInner
void simulateInner(void)
Firmware simulation for the outers. yi.
Definition: TrackSegmentFinder.cc:1441
Belle2::TRGClock::max
int max(void) const
returns max. clock point.
Definition: Clock.h:201
Belle2::TRGCDC::nSegments
unsigned nSegments(void) const
returns # of track segments.
Definition: TRGCDC.h:954
Belle2::TRGCDCTrackSegmentFinder::saveTSFResults
void saveTSFResults(std::vector< TRGCDCSegmentHit * > *segmentHitsSL)
save result of TSF
Definition: TrackSegmentFinder.cc:550
Belle2::TRGCDCTrackSegmentFinder::addID
static void addID(TRGSignalVector &s, unsigned id)
Add TSF ID to timing signal vector for the output.
Definition: TrackSegmentFinder.cc:1937
Belle2::TRGCDCWire::hit
const TRGCDCWireHit * hit(void) const
returns a pointer to a TRGCDCWireHit.
Definition: Wire.h:144
Belle2::TRGCDC::nSuperLayers
unsigned nSuperLayers(void) const
returns # of super layers.
Definition: TRGCDC.h:870
Belle2::TRGCDCTrackSegmentFinder::saveTSInformation
void saveTSInformation(std::vector< TRGCDCSegment * > &tss)
save the TS info
Definition: TrackSegmentFinder.cc:411
Belle2::TRGCDC::segment
const TRGCDCSegment & segment(unsigned id) const
returns a track segment.
Definition: TRGCDC.h:933
Belle2::TRGCDCTrackSegmentFinder::hitMapOuter
void hitMapOuter(void)
Creates the hit maps for the outer.
Definition: TrackSegmentFinder.cc:2188
Belle2::TRGChannel::signal
TRGSignalBundle * signal(void) const
returns signal.
Definition: Channel.h:93
Belle2::TRGCDCTrackSegmentFinder::simulate2
void simulate2(void)
Firmware simulation. Unified version of inner and outer : yi.
Definition: TrackSegmentFinder.cc:1992
Belle2::TRGSignal::state
bool state(int clockPosition) const
returns true if signal is active in given clock position.
Definition: Signal.h:286
Belle2::TRGCDCTrackSegmentFinder::doit
void doit(std::vector< TRGCDCSegment * > &tss, const bool trackSegmentClockSimulation, std::vector< TRGCDCSegmentHit * > &segmentHits, std::vector< TRGCDCSegmentHit * > *segmentHitsSL)
Member functions of doing TSF.
Definition: TrackSegmentFinder.cc:133
Belle2::TRGDebug::level
static int level(void)
returns the debug level.
Definition: Debug.cc:67
Belle2::TRGSignal::active
bool active(void) const
returns true if there is a signal.
Definition: Signal.h:277
Belle2::TRGClock::clockCounter
const TRGSignalVector & clockCounter(void) const
returns the clock counter.
Definition: Clock.cc:166
Belle2::TRGCDCTrackSegmentFinder::inputInner
void inputInner(const unsigned id, const unsigned nTSF, TRGSignalVector *s)
Creates input signals to TSF for the inner.
Definition: TrackSegmentFinder.cc:2215
Belle2::TRGCDC::layer
const TRGCDCLayer * layer(unsigned id) const
returns a pointer to a layer. 0 will be returned if 'id' is invalid.
Definition: TRGCDC.h:828
Belle2::TRGDebug::leaveStage
static void leaveStage(const std::string &stageName)
Declare that you leave a stage.
Definition: Debug.cc:34
Belle2::TRGCDCTrackSegmentFinder::simulateTSF2
vector< TRGSignalVector * > simulateTSF2(TRGSignalVector *in, unsigned id)
Simulate TSF response (unified version, state machine)
Definition: TrackSegmentFinder.cc:2811
Belle2::TRGCDCLayer::nCells
unsigned nCells(void) const
returns # of cells.
Definition: Layer.h:194
Belle2::TRGCDCWireHit::iCDCSimHit
unsigned iCDCSimHit(void) const
return index of CDCSimHit
Definition: WireHit.h:111
Belle2::TRGCDCTrackSegmentFinder::saveNNTSInformation
void saveNNTSInformation(std::vector< TRGCDCSegment * > &tss)
Saves NNTS information. Only when ts is hit.
Definition: TrackSegmentFinder.cc:717
Belle2::TRGCDCTrackSegmentFinder::simulate
void simulate(void)
Firmware simulation. yi.
Definition: TrackSegmentFinder.cc:1954
Belle2::TRGCDCTrackSegmentFinder::hitMapInner
void hitMapInner(void)
Creates the hit maps for the inner.
Definition: TrackSegmentFinder.cc:2149
Belle2::TRGCDCTrackSegmentFinder::packerOuterEvt
TRGSignalVector * packerOuterEvt(vector< TRGSignalVector * >, vector< int >, int)
Packing output for evtTime & Low pT.
Definition: TrackSegmentFinder.cc:355
Belle2::TRGCDCTrackSegmentFinder::mkint
double mkint(TRGState)
tranformatoin into integer
Definition: TrackSegmentFinder.cc:381
Belle2::TRGBoard::nOutput
unsigned nOutput(void) const
returns output channels.
Definition: Board.h:180
Belle2::TRGCDCTrackSegmentFinder::type
boardType type(void) const
board type of TSF
Definition: TrackSegmentFinder.cc:121
Belle2::TRGBoard::output
TRGChannel * output(unsigned i) const
returns output channel i.
Definition: Board.h:165
Belle2::TRGCDCTrackSegmentFinder::fastestTimingInner
void fastestTimingInner(unsigned tsfID, const unsigned nTSF, TRGSignalVector &s) const
Make TRGSignals for the fastest timing bits.
Definition: TrackSegmentFinder.cc:1675
Belle2::TRGCDCTrackSegmentFinder::mkbool
vector< bool > mkbool(int, int)
tranformatoin into bool
Definition: TrackSegmentFinder.cc:395
Belle2::TRGBoard::dump
void dump(const std::string &message="", const std::string &pre="") const
dumps contents.
Definition: Board.cc:42
Belle2
Abstract base class for different kinds of events.
Definition: MillepedeAlgorithm.h:17