#include <iostream>
#include <fstream>
#include <sstream>
#include <cinttypes>
#include <iomanip>
#include <cstdlib>
#include <cstdio>
#include <cstring>

Functions
int	main (int argc, char *argv[])

Variables
uint16_t	mem [nChs][nEvts]

const int	nChs = 68

const int	nEvts = 2048

Detailed Description

The TCC memory for FE data emulation takes a fixed number, 2048, of events. This standalone application completes a FE emulation data file with an arbitrary number of events (<=2048) in order to have the required 2048 events. The N initial events are repeated till having 2048 events. In general the number of initial events is choosen as a divider of 2048.

Definition in file recycleTccEmu.cc.

Function Documentation

◆ main()

int main	(	int	argc,
		char *	argv[]
	)

===============================================================================================================================================================================================

variant2: for each run define phi-averaged A for normalization channel (Dref,16) and then, divide Rijk on it, i.e. get RRijk

eta=27

eta=25

eta=23

eta=22

eta=21

eta=26

eta=24

eta=19

eta=17

eta=25

eta=23

eta=22

eta=21

eta=26

eta=24

eta=20

eta=19

eta=18

eta=27 L1=1

eta=25 L1=1

eta=23 L1=1

eta=22 L1=1

eta=21 L1=1

eta=29 L1=1

eta=26 L1=1

eta=24 L1=1

eta=20 L1=1

eta=19 L1=1

eta=18 L1=1

eta=17 L1=1

eta=28 L7=1

eta=27 L7=1

eta=25 L7=1

eta=23 L7=1

eta=22 L7=1

eta=21 L7=1

eta=26 L7=1

eta=24 L7=1

eta=20 L7=1

eta=19 L7=1

eta=18 L7=1

eta=17 L7=1

eta=27

eta=25

eta=23

eta=22

eta=21

eta=26

eta=24

eta=19

eta=17

eta=25

eta=23

eta=22

eta=21

eta=26

eta=24

eta=20

eta=19

eta=18

eta=27 L1=1

eta=25 L1=1

eta=23 L1=1

eta=22 L1=1

eta=21 L1=1

eta=26 L1=1

eta=24 L1=1

eta=20 L1=1

eta=19 L1=1

eta=18 L1=1

eta=17 L1=1

eta=28 L7=1

eta=27 L7=1

eta=25 L7=1

eta=23 L7=1

eta=22 L7=1

eta=21 L7=1

eta=26 L7=1

eta=24 L7=1

eta=20 L7=1

eta=19 L7=1

eta=18 L7=1

eta=17 L7=1

eta=27

eta=28

errA with average Amplitudes

Summed Amplitude Plots:

Summed Amplitude Plots:

Summed Amplitude Plots:

Sticking ADC counts in neughbouring TS for HB:

Sticking ADC counts in neughbouring TS for HE:

RBX:

errA with average Amplitudes

Summed Amplitude Plots:

Summed Amplitude Plots:

Summed Amplitude Plots:

RBX:

errA with average Amplitudes

Summed Amplitude Plots:

Summed Amplitude Plots:

Summed Amplitude Plots:

RBX:

Prepare maps of good/bad channels:

CALO JETS

PF JETS

Definition at line 23 of file recycleTccEmu.cc.

References dir2webdir::argc, GCPpyPlots::argv, gather_cfg::cout, TauDecayModes::dec, beamvalidation::exit(), recoMuon::in, mem, nChs, nEvts, MillePedeFileConverter_cfg::out, alignCSCRings::s, and heppy_batch::val.

                                  {
   if ((argc >= 2 && ((strcmp(argv[1], "-h") == 0) || (strcmp(argv[1], "--help") == 0))) || argc != 3) {
     cout << "Usage: recycleTccEmu infile outfile\n";
     return 1;
   }
 
   string ifilename = argv[1];
   string ofilename = argv[2];
 
   for (int iCh = 0; iCh < nChs; ++iCh) {
     for (int iEvts = 0; iEvts < nEvts; ++iEvts) {
       mem[iCh][iEvts] = 0xFFFF;
     }
   }
 
   ifstream in(ifilename.c_str());
   int chnb;
   int bcnb;
   int val;
   int dummy;
   int oldLineCnt = 0;
 
   //reads input file:
   if (in) {
     while (!in.eof()) {
       in >> dec >> chnb >> bcnb >> hex >> val >> dummy;
       mem[chnb - 1][bcnb] = val & 0x7FF;
       if (mem[chnb - 1][bcnb] != val) {
         cout << "Invalid Et value at line " << oldLineCnt + 1 << ".\n";
         exit(1);
       }
       // cout<<"Channel: "<< dec <<chnb <<", BX: "
       // << dec << bcnb << " filled with val:"<< hex<< mem[chnb-1][bcnb]
       // << dec << endl;
       ++oldLineCnt;
     }
   } else {
     cout << "Failed to open file " << ifilename << "\n";
   }
 
   in.close();
   ofstream out(ofilename.c_str());
 
   if (!out) {
     cout << "Failed to open file '" << ofilename << "' in write mode.\n";
     return 1;
   }
 
   bool singleOldEventCnt = true;
   int oldEventCnt[nChs];
   //fills output file:
   for (int iCh = 0; iCh < nChs; ++iCh) {
     int evtcnt = 0;
     //find first not initialized events:
     while (evtcnt < nEvts && mem[iCh][evtcnt] != 0xFFFF) {
       ++evtcnt;
     }
     //cout << "ch " << iCh << " event count: " << evtcnt << "\n";
     oldEventCnt[iCh] = evtcnt;
     if (oldEventCnt[0] != oldEventCnt[iCh])
       singleOldEventCnt = false;
     if (evtcnt == 0) {
       cout << "Error: no data found for channel " << iCh << "\n";
     }
     //clones data of channel iCh
     for (int ievt = evtcnt; ievt < nEvts; ++ievt) {
       if (mem[iCh][ievt] != 0xFFFF) {
         cout << "Error: memory offset of channel " << iCh << " events are not contiguous.\n";
         exit(1);
       }
       mem[iCh][ievt] = mem[iCh][ievt % evtcnt];
     }
 
     for (int ievt = 0; ievt < nEvts; ++ievt) {
       out << iCh + 1 << "\t" << ievt << "\t" << hex << "0x" << setfill('0') << setw(4) << mem[iCh][ievt] << setfill(' ')
           << dec << "\t0"
           << "\n";
     }
   }
 
   //warning for aperiodic case:
   if (singleOldEventCnt && (nEvts % oldEventCnt[0] != 0)) {
     cout << "Warning: ouput event count (2048) is not a mulitple of input "
             "event counts\n";
   }
   if (!singleOldEventCnt) {
     stringstream s;
     for (int iCh = 0; iCh < nChs; ++iCh) {
       if (nEvts % oldEventCnt[iCh]) {
         s << (s.str().empty() ? "" : ", ") << iCh;
       }
     }
     if (!s.str().empty())
       cout << "Warning: ouput event count (2048) for channel" << (s.str().size() > 1 ? "s" : "") << " " << s.str()
            << " is not a mulitple of input event counts\n";
   }
 
   if (!singleOldEventCnt) {
     cout << "Info: in the input file the event count depends on the channel";
   }
 }