trgReader20.cxx

#include <stdio.h>
#include <string.h>
#include <arpa/inet.h>

#include <rtsLog.h>


#include <rtsSystems.h>

// one needs to diffuse the crummy trigger stuff first
#define TRG_VERSION     0x20
#include <daqFormats.h>


#include <rts.h>  // fror the sake of swap32 entry

#include "daq_trg.h"

// this maps from the raw data to the following representation
// sector1,slat1,ch1 ... ch5, sector1,slat2,ch1.... ch5, sector2,slat1,ch1...

static unsigned char ctbMap[240] = {
7, 6, 5, 4, 3, 23, 22, 21, 20, 19, 
2, 1, 0, 15, 14, 18, 17, 16, 31, 30, 
13, 12, 11, 10, 9, 29, 28, 27, 26, 25, 
39, 38, 37, 36, 35, 55, 54, 53, 52, 51, 
34, 33, 32, 47, 46, 50, 49, 48, 63, 62, 
45, 44, 43, 42, 41, 61, 60, 59, 58, 57, 
71, 70, 69, 68, 67, 87, 86, 85, 84, 83, 
66, 65, 64, 79, 78, 82, 81, 80, 95, 94, 
77, 76, 75, 74, 73, 93, 92, 91, 90, 89, 
103, 102, 101, 100, 99, 119, 118, 117, 116, 115, 
98, 97, 96, 111, 110, 114, 113, 112, 127, 126, 
109, 108, 107, 106, 105, 125, 124, 123, 122, 121, 
135, 134, 133, 132, 131, 151, 150, 149, 148, 147, 
130, 129, 128, 143, 142, 146, 145, 144, 159, 158, 
141, 140, 139, 138, 137, 157, 156, 155, 154, 153, 
167, 166, 165, 164, 163, 183, 182, 181, 180, 179, 
162, 161, 160, 175, 174, 178, 177, 176, 191, 190, 
173, 172, 171, 170, 169, 189, 188, 187, 186, 185, 
199, 198, 197, 196, 195, 215, 214, 213, 212, 211, 
194, 193, 192, 207, 206, 210, 209, 208, 223, 222, 
205, 204, 203, 202, 201, 221, 220, 219, 218, 217, 
231, 230, 229, 228, 227, 247, 246, 245, 244, 243, 
226, 225, 224, 239, 238, 242, 241, 240, 255, 254, 
237, 236, 235, 234, 233, 253, 252, 251, 250, 249, 
} ;


static unsigned char mwcMap[96] = {
71, 70, 69, 68, 67, 66, 65, 64, 79, 78, 77, 76, 
95, 94, 93, 92, 87, 86, 85, 84, 83, 82, 81, 80, 
99, 98, 97, 96, 111, 110, 109, 108, 103, 102, 101, 100, 
119, 118, 117, 116, 115, 114, 113, 112, 127, 126, 125, 124, 
7, 6, 5, 4, 3, 2, 1, 0, 15, 14, 13, 12, 
31, 30, 29, 28, 23, 22, 21, 20, 19, 18, 17, 16, 
35, 34, 33, 32, 47, 46, 45, 44, 39, 38, 37, 36, 
55, 54, 53, 52, 51, 50, 49, 48, 63, 62, 61, 60, 
} ;




// read the Trigger RAW data
int trgReader20(char *arg, struct trg_t *trg)
{
        int i ;
        int sdes, ssum, sraw ;
        int len ;

        struct TRGD *trgd = (struct TRGD *)arg ;

        int swap = 0;
        if(trgd->bh.byte_order != DAQ_RAW_FORMAT_ORDER) swap=1;

        len = swap ? swap32(trgd->bh.length)*4 : trgd->bh.length;

        sdes = sizeof(trgd->bh) + sizeof(trgd->desc) ;
        ssum = sdes + sizeof(trgd->sum) ;
        sraw = ssum + sizeof(trgd->raw[0]) ;

        if(len < ssum) {
                LOG(NOTE,"No TRG Summaries, or Raw...",0,0,0,0,0) ;
                return 0 ;
        }
        else if(len < sraw) {
                LOG(NOTE,"TRG: Summaries ONLY!",0,0,0,0,0) ;
                return 0 ;
        }

        LOG(DBG,"len %d, desc len %d, sum len %d, raw trg 0 len %d, all raw len %d",
            len,sizeof(trgd->desc),sizeof(trgd->sum),sizeof(trgd->raw[0]),sizeof(trgd->raw)) ;

        LOG(DBG,"evt desc bytes %d, sum bytes %d, raw bytes %d",trgd->desc.TCUdataBytes,
            trgd->sum.TrgSumBytes,trgd->raw[0].RawDetBytes,0,0) ;

        int pre, post ;

        pre = swap ? swap16(trgd->desc.npre) : trgd->desc.npre ;
        post = swap ? swap16(trgd->desc.npost) : trgd->desc.npost ;

        int trg_length = len - sizeof(struct bankHeader) - 40 ; // 40 is supurious
        int exp_length = sizeof(EvtDescData) + sizeof(trgd->sum) + sizeof(trgd->raw[0]) + 
                        sizeof(trgd->raw[0])*pre + sizeof(trgd->raw[0])*post ;

        if(trg_length != exp_length) {
                LOG(NOTE,"Trigger data: is %d, expect %d bytes, pre %d, post %d",
                    trg_length, exp_length,pre,post,0) ;
        }
        else {
                LOG(DBG,"Trigger data: is %d, expect %d bytes, pre %d, post %d",
                    trg_length, exp_length,pre,post,0) ;
        }

        // get the last DSM aka TCU bits
        trg->tcubits = swap ? swap16(trgd->desc.DSMInput) : trgd->desc.DSMInput;
        trg->detlive = swap ? swap16(trgd->desc.externalBusy) : trgd->desc.externalBusy ;


        LOG(DBG,"TrgDataFmtVer 0x%X, 0x%04X %c %c", trgd->desc.TrgDataFmtVer,
            trgd->sum.TrgSumHeader,
            trgd->raw[0].CTBdataHeader[0], 
            trgd->raw[0].CTBdataHeader[1],0) ;
LOG(DBG,"TrgSumBytes %d, 0x%04X %c %c %c", trgd->sum.TrgSumBytes,
            trgd->sum.L0SumHeader,
            trgd->raw[0].RawDetHeader[0], 
            trgd->raw[0].RawDetHeader[1],0) ;

        for(i=0;i<240;i++) {
                trg->CTB[i] = trgd->raw[0].CTB[ctbMap[i]] ;
        }

        for(i=0;i<96;i++) {
                trg->MWC[i] = trgd->raw[0].MWC[mwcMap[i]] ;
        }

        // east
        for(i=0;i<240;i++) {
                trg->BEMC[0][i] = trgd->raw[0].BEMCEast[i] ;
        }
        // west
        for(i=0;i<240;i++) {
                trg->BEMC[1][i] = trgd->raw[0].BEMCWest[i] ;
        }
        // layer1
        for(i=0;i<48;i++) {
                trg->BEMC_l1[i] = swap ? swap16(trgd->raw[0].BEMClayer1[i]) : trgd->raw[0].BEMClayer1[i] ;
        }



        for(i=0;i<144;i++) {
                trg->EEMC[i] = trgd->raw[0].EEMC[i] ;
        }
        for(i=0;i<16;i++) {
                trg->EEMC_l1[i] = swap ? swap16(trgd->raw[0].EEMClayer1[i]) :  trgd->raw[0].EEMClayer1[i];
        }

        // FPD:
        // first index is 0 East, 1 West
        // second index is 0 NortSouth, 1 TopBottom
        // East, NS
        for(i=0;i<112;i++) {
                trg->FPD[0][0][i] = trgd->raw[0].FPDEastNSLayer0[i] ;
        }
        for(i=0;i<8;i++) {
          trg->FPD_l1[0][0][i] = swap ? swap16(trgd->raw[0].FPDEastNSLayer1[i]) : trgd->raw[0].FPDEastNSLayer1[i] ;
        }
        // East, TB
        for(i=0;i<64;i++) {
                trg->FPD[0][1][i] = trgd->raw[0].FPDEastTBLayer0[i] ;
        }
        for(i=0;i<8;i++) {
                trg->FPD_l1[0][1][i] = swap ? swap16(trgd->raw[0].FPDEastTBLayer1[i]) : trgd->raw[0].FPDEastTBLayer1[i] ;
        }
        // West, NS
        for(i=0;i<112;i++) { 
                trg->FPD[1][0][i] = trgd->raw[0].FPDWestNSLayer0[i] ;
        }
        for(i=0;i<8;i++) {
                trg->FPD_l1[1][0][i] = swap ? swap16(trgd->raw[0].FPDWestNSLayer1[i]) : trgd->raw[0].FPDWestNSLayer1[i];
        }
        // West, TB
        for(i=0;i<64;i++) {
                trg->FPD[1][1][i] = trgd->raw[0].FPDWestTBLayer0[i] ;
        }
        for(i=0;i<8;i++) {
                trg->FPD_l1[1][1][i] = swap ? swap16(trgd->raw[0].FPDWestTBLayer1[i]) : trgd->raw[0].FPDWestTBLayer1[i];
        }




        for(i=0;i<80;i++) {
                trg->BBC[i] = trgd->raw[0].BBC[i] ;
        }
        for(i=0;i<16;i++) {
                trg->BBC_l1[i] = swap ? swap16(trgd->raw[0].BBClayer1[i]) : trgd->raw[0].BBClayer1[i] ;
        }


        for(i=0;i<16;i++) {
                trg->ZDC[i] = trgd->raw[0].ZDC[i] ;
        }
        for(i=0;i<8;i++) {
                trg->ZDC_l1[i] = swap ? swap16(trgd->raw[0].ZDClayer1[i]) : trgd->raw[0].ZDClayer1[i];
        }

        trg->trg_sum = NULL ;


/*
        int j ;
        for(j=0;j<11;j++) {
        for(i=0;i<16;i++) {
                printf("pre %d, ZDC %d = %u\n",j,i,trgd->raw[j].ZDC[i]) ;
        }
        for(i=0;i<8;i++) {
                printf("pre %d, ZDC_l1 %d = %u\n",j,i,b2h16(trgd->raw[j].ZDClayer1[i])) ;
        }
        }
*/

        return len ;
}