#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:

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";
   }
 }