#include "Record_Info.h"
#include <stdint.h>
#include "zdab_utils.h"
#include "pack2b.h"
#include <stdlib.h> /* for size_t */
#include <stdio.h> /* for fprintf() */

int isOrphan(aPmtEventRecord *pmtRecord)
{
    /* Returns non-zero if the specified event is an orphan. */
    int i;

    uint32_t *mtc_data = (uint32_t *) &pmtRecord->TriggerCardData;

    for (i = 0; i < 6; ++i) {
        if (*mtc_data != 0) return 0;
        ++mtc_data;
    }

    return 1;
}

// PH 04/23/98
// Swap 4-byte integer/floats between native and external format
void swap_int32(int32_t *val_pt, int count)
{
    int32_t *last = val_pt + count;
    while (val_pt < last) {
        *val_pt = ((*val_pt << 24) & 0xff000000) |
                  ((*val_pt <<  8) & 0x00ff0000) |
                  ((*val_pt >>  8) & 0x0000ff00) |
                  ((*val_pt >> 24) & 0x000000ff);
        ++val_pt;
    }
    return;
}

// Swap 2-byte integers between native and external format
void swap_int16(int16_t *val_pt, int count)
{
    char tmp;
    int i;
    for (i=0; i<count; ++i) {
        tmp = ((char *) val_pt)[0];
        ((char *) val_pt)[0] = ((char *) val_pt)[1];
        ((char *) val_pt)[1] = tmp;
        ++val_pt;
    }
    return;
}

void unpack_mcvx(uint32_t *data, MCVXBank *b)
{
    unpack((uint8_t *) data,      "l",&b->cls);
    unpack((uint8_t *) (data+1),  "l",&b->inc);
    unpack((uint8_t *) (data+2),  "f",&b->x);
    unpack((uint8_t *) (data+3),  "f",&b->y);
    unpack((uint8_t *) (data+4),  "f",&b->z);
    unpack((uint8_t *) (data+5),  "F",&b->tim);
    unpack((uint8_t *) (data+7),  "l",&b->rgn);
    unpack((uint8_t *) (data+8),  "l",&b->idm);
    unpack((uint8_t *) (data+9),  "l",&b->rg2);
    unpack((uint8_t *) (data+10), "l",&b->im2);
    unpack((uint8_t *) (data+12), "f",&b->bnx);
    unpack((uint8_t *) (data+13), "f",&b->bny);
    unpack((uint8_t *) (data+14), "f",&b->bnz);
    unpack((uint8_t *) (data+15), "l",&b->cer);
}

void unpack_mctk(uint32_t *data, MCTKBank *b)
{
    unpack((uint8_t *) data,      "l",&b->idp);
    unpack((uint8_t *) (data+1),  "f",&b->drx);
    unpack((uint8_t *) (data+2),  "f",&b->dry);
    unpack((uint8_t *) (data+3),  "f",&b->drz);
    unpack((uint8_t *) (data+4),  "f",&b->ene);
    unpack((uint8_t *) (data+5),  "l",&b->rgn);
    unpack((uint8_t *) (data+6),  "l",&b->idm);
    unpack((uint8_t *) (data+7),  "f",&b->plx);
    unpack((uint8_t *) (data+8),  "f",&b->ply);
    unpack((uint8_t *) (data+9),  "f",&b->plz);
    unpack((uint8_t *) (data+10), "f",&b->stp);
    unpack((uint8_t *) (data+11), "f",&b->near);
}

void unpack_ev(uint32_t *data, EVBank *b)
{
    unpack((uint8_t *) data,     "l",&b->run);
    unpack((uint8_t *) (data+1), "l",&b->evn);
    unpack((uint8_t *) (data+2), "l",&b->dtp);
    unpack((uint8_t *) (data+3), "l",&b->jdy);
    unpack((uint8_t *) (data+4), "l",&b->ut1);
    unpack((uint8_t *) (data+5), "l",&b->ut2);
    unpack((uint8_t *) (data+6), "l",&b->dte);
    unpack((uint8_t *) (data+7), "l",&b->hmsc);
    unpack((uint8_t *) (data+8), "l",&b->gtr1);
    unpack((uint8_t *) (data+9), "l",&b->gtr2);
    unpack((uint8_t *) (data+10),"l",&b->npm);
    unpack((uint8_t *) (data+11),"l",&b->nph);
    unpack((uint8_t *) (data+12),"l",&b->sub_run);
    unpack((uint8_t *) (data+13),"l",&b->mc_pck);
    unpack((uint8_t *) (data+14),"l",&b->rec);
    unpack((uint8_t *) (data+15),"l",&b->vpck);
    unpack((uint8_t *) (data+16),"l",&b->gtr_id);
    unpack((uint8_t *) (data+17),"l",&b->trg_type);
    unpack((uint8_t *) (data+18),"l",&b->peak);
    unpack((uint8_t *) (data+19),"l",&b->diff);
    unpack((uint8_t *) (data+20),"l",&b->integral);
    unpack((uint8_t *) (data+21),"l",&b->err);
    unpack((uint8_t *) (data+22),"l",&b->data_set);
    unpack((uint8_t *) (data+22),"lll",&b->spare1[0],
                                       &b->spare1[1],
                                       &b->spare1[2]);
    unpack((uint8_t *) (data+26),"l",&b->ncd_status);
    unpack((uint8_t *) (data+27),"l",&b->num_muxg);
    unpack((uint8_t *) (data+29),"l",&b->num_mux);
    unpack((uint8_t *) (data+29),"l",&b->num_scope);
    unpack((uint8_t *) (data+30),"lllll",&b->spare2[0],
                                         &b->spare2[1],
                                         &b->spare2[2],
                                         &b->spare2[3],
                                         &b->spare2[4]);
    unpack((uint8_t *) (data+35),"l",&b->ncd_clk_up);
    unpack((uint8_t *) (data+36),"l",&b->ncd_clk_lw);
    unpack((uint8_t *) (data+37),"l",&b->ncd_reg);
    unpack((uint8_t *) (data+38),"l",&b->ncd_gtid);
    unpack((uint8_t *) (data+39),"l",&b->ncd_sync);
    unpack((uint8_t *) (data+40),"l",&b->spare3[0],
                                     &b->spare3[1],
                                     &b->spare3[2],
                                     &b->spare3[3],
                                     &b->spare3[4],
                                     &b->spare3[5],
                                     &b->spare3[6],
                                     &b->spare3[7],
                                     &b->spare3[8],
                                     &b->spare3[9]);
}

void unpack_pmt(uint32_t *data, PMTBank *b)
{
    unpack((uint8_t *) data,"l",&b->pn);
    unpack((uint8_t *) (data+1),"l",&b->pf);
    unpack((uint8_t *) (data+2),"f",&b->pt);
    unpack((uint8_t *) (data+3),"f",&b->phl);
    unpack((uint8_t *) (data+4),"f",&b->phs);
    unpack((uint8_t *) (data+5),"f",&b->plx);
    unpack((uint8_t *) (data+6),"f",&b->pt0);
    unpack((uint8_t *) (data+7),"l",&b->pif);
    unpack((uint8_t *) (data+8),"f",&b->pit);
    unpack((uint8_t *) (data+9),"f",&b->pihl);
    unpack((uint8_t *) (data+10),"f",&b->pihs);
    unpack((uint8_t *) (data+11),"f",&b->pilx);
    unpack((uint8_t *) (data+12),"f",&b->pit0);
    unpack((uint8_t *) (data+13),"l",&b->cell);
    unpack((uint8_t *) (data+14),"l",&b->pin);
    unpack((uint8_t *) (data+15),"l",&b->tslh);
    unpack((uint8_t *) (data+16),"l",&b->hca);
    unpack((uint8_t *) (data+17),"l",&b->eca_val);
    unpack((uint8_t *) (data+18),"l",&b->pca_val);
    unpack((uint8_t *) (data+19),"l",&b->anxx);
    unpack((uint8_t *) (data+20),"l",&b->ept);
    unpack((uint8_t *) (data+21),"l",&b->ehl);
    unpack((uint8_t *) (data+22),"l",&b->ehs);
    unpack((uint8_t *) (data+23),"l",&b->elx);
    unpack((uint8_t *) (data+24),"l",&b->pt1);
    unpack((uint8_t *) (data+25),"l",&b->ptm);
    unpack((uint8_t *) (data+26),"l",&b->ptms);
    unpack((uint8_t *) (data+27),"l",&b->qm);
    unpack((uint8_t *) (data+28),"l",&b->qms);
    unpack((uint8_t *) (data+29),"l",&b->qrc);
}

int swap_PmtRecord(aPmtEventRecord *aPmtRecord, size_t size)
{
    /* Swap a Pmt Event Record. This function swaps both the Pmt event record
     * and the PMT hits and sub fields. Returns -1 if the PMT record has too
     * many hits. */
    SWAP_INT32(aPmtRecord, sizeof(aPmtEventRecord)/sizeof(uint32_t));

    int npmt = aPmtRecord->NPmtHit;

    if (npmt > MAX_NHIT) {
        fprintf(stderr, "Read error: Bad ZDAB -- %d pmt hit!", npmt);
        return -1;
    } else {
        if (size < sizeof(aPmtEventRecord) + 3*npmt*4) {
            fprintf(stderr, "swap_PmtRecord: size of record is %zu bytes, but there are %i PMT hits", size, npmt);
            return -1;
        }
        // swap the hit data
        SWAP_INT32(aPmtRecord + 1, 3*npmt);
        // swap the sub-fields
        uint32_t *sub_header = &aPmtRecord->CalPckType;
        while (*sub_header & SUB_NOT_LAST) {
            if (size < (sub_header - (uint32_t *) aPmtRecord)*4 + (*sub_header & SUB_LENGTH_MASK)*4 + 4) {
                fprintf(stderr, "swap_PmtRecord: size of record is %zu bytes, "
                    "but sub-field requires %lu bytes",
                    size,
                    (sub_header - (uint32_t *) aPmtRecord)*4 + (*sub_header & SUB_LENGTH_MASK)*4 + 4);
                return -1;
            }
            sub_header += (*sub_header & SUB_LENGTH_MASK);
            SWAP_INT32(sub_header, 1); // swap the sub-field header
            // get number of data words (-1 because we don't want to include header size)
            uint32_t data_words = (*sub_header & SUB_LENGTH_MASK) - 1;
            if (size < (sub_header - (uint32_t *) aPmtRecord)*4 + (*sub_header & SUB_LENGTH_MASK)*4) {
                fprintf(stderr, "swap_PmtRecord: size of record is %zu bytes, "
                    "but sub-field requires %lu bytes",
                    size,
                    (sub_header - (uint32_t *) aPmtRecord)*4 + (*sub_header & SUB_LENGTH_MASK)*4);
                return -1;
            }
            SWAP_INT32(sub_header+1, data_words);
        }
    }

    return 0;
}

void swap_TrigRecord(struct TriggerInfo *aTrigRecord)
{
    /* Byte swap a Trigger Record. */
    SWAP_INT32(aTrigRecord, sizeof(struct TriggerInfo)/sizeof(uint32_t));
}

void swap_RunRecord(struct RunRecord *aRunRecord)
{
    /* Byte swap a Run Record. */
    SWAP_INT32(aRunRecord, sizeof(struct RunRecord)/sizeof(uint32_t));
}


>197</a>
<a id='n198' href='#n198'>198</a>
<a id='n199' href='#n199'>199</a>
<a id='n200' href='#n200'>200</a>
<a id='n201' href='#n201'>201</a>
<a id='n202' href='#n202'>202</a>
<a id='n203' href='#n203'>203</a>
<a id='n204' href='#n204'>204</a>
<a id='n205' href='#n205'>205</a>
<a id='n206' href='#n206'>206</a>
<a id='n207' href='#n207'>207</a>
<a id='n208' href='#n208'>208</a>
<a id='n209' href='#n209'>209</a>
<a id='n210' href='#n210'>210</a>
<a id='n211' href='#n211'>211</a>
<a id='n212' href='#n212'>212</a>
<a id='n213' href='#n213'>213</a>
<a id='n214' href='#n214'>214</a>
<a id='n215' href='#n215'>215</a>
<a id='n216' href='#n216'>216</a>
<a id='n217' href='#n217'>217</a>
<a id='n218' href='#n218'>218</a>
<a id='n219' href='#n219'>219</a>
<a id='n220' href='#n220'>220</a>
<a id='n221' href='#n221'>221</a>
<a id='n222' href='#n222'>222</a>
<a id='n223' href='#n223'>223</a>
<a id='n224' href='#n224'>224</a>
<a id='n225' href='#n225'>225</a>
<a id='n226' href='#n226'>226</a>
<a id='n227' href='#n227'>227</a>
<a id='n228' href='#n228'>228</a>
<a id='n229' href='#n229'>229</a>
<a id='n230' href='#n230'>230</a>
<a id='n231' href='#n231'>231</a>
<a id='n232' href='#n232'>232</a>
<a id='n233' href='#n233'>233</a>
<a id='n234' href='#n234'>234</a>
<a id='n235' href='#n235'>235</a>
<a id='n236' href='#n236'>236</a>
<a id='n237' href='#n237'>237</a>
<a id='n238' href='#n238'>238</a>
<a id='n239' href='#n239'>239</a>
<a id='n240' href='#n240'>240</a>
<a id='n241' href='#n241'>241</a>
<a id='n242' href='#n242'>242</a>
<a id='n243' href='#n243'>243</a>
<a id='n244' href='#n244'>244</a>
<a id='n245' href='#n245'>245</a>
<a id='n246' href='#n246'>246</a>
<a id='n247' href='#n247'>247</a>
<a id='n248' href='#n248'>248</a>
<a id='n249' href='#n249'>249</a>
<a id='n250' href='#n250'>250</a>
<a id='n251' href='#n251'>251</a>
<a id='n252' href='#n252'>252</a>
<a id='n253' href='#n253'>253</a>
<a id='n254' href='#n254'>254</a>
<a id='n255' href='#n255'>255</a>
<a id='n256' href='#n256'>256</a>
<a id='n257' href='#n257'>257</a>
<a id='n258' href='#n258'>258</a>
<a id='n259' href='#n259'>259</a>
<a id='n260' href='#n260'>260</a>
<a id='n261' href='#n261'>261</a>
<a id='n262' href='#n262'>262</a>
<a id='n263' href='#n263'>263</a>
<a id='n264' href='#n264'>264</a>
<a id='n265' href='#n265'>265</a>
<a id='n266' href='#n266'>266</a>
<a id='n267' href='#n267'>267</a>
<a id='n268' href='#n268'>268</a>
<a id='n269' href='#n269'>269</a>
<a id='n270' href='#n270'>270</a>
<a id='n271' href='#n271'>271</a>
<a id='n272' href='#n272'>272</a>
<a id='n273' href='#n273'>273</a>
<a id='n274' href='#n274'>274</a>
<a id='n275' href='#n275'>275</a>
<a id='n276' href='#n276'>276</a>
<a id='n277' href='#n277'>277</a>
<a id='n278' href='#n278'>278</a>
<a id='n279' href='#n279'>279</a>
<a id='n280' href='#n280'>280</a>
<a id='n281' href='#n281'>281</a>
<a id='n282' href='#n282'>282</a>
<a id='n283' href='#n283'>283</a>
<a id='n284' href='#n284'>284</a>
<a id='n285' href='#n285'>285</a>
<a id='n286' href='#n286'>286</a>
<a id='n287' href='#n287'>287</a>
<a id='n288' href='#n288'>288</a>
<a id='n289' href='#n289'>289</a>
<a id='n290' href='#n290'>290</a>
<a id='n291' href='#n291'>291</a>
<a id='n292' href='#n292'>292</a>
<a id='n293' href='#n293'>293</a>
<a id='n294' href='#n294'>294</a>
<a id='n295' href='#n295'>295</a>
<a id='n296' href='#n296'>296</a>
<a id='n297' href='#n297'>297</a>
<a id='n298' href='#n298'>298</a>
<a id='n299' href='#n299'>299</a>
<a id='n300' href='#n300'>300</a>
<a id='n301' href='#n301'>301</a>
<a id='n302' href='#n302'>302</a>
<a id='n303' href='#n303'>303</a>
<a id='n304' href='#n304'>304</a>
<a id='n305' href='#n305'>305</a>
<a id='n306' href='#n306'>306</a>
<a id='n307' href='#n307'>307</a>
<a id='n308' href='#n308'>308</a>
<a id='n309' href='#n309'>309</a>
<a id='n310' href='#n310'>310</a>
<a id='n311' href='#n311'>311</a>
<a id='n312' href='#n312'>312</a>
<a id='n313' href='#n313'>313</a>
<a id='n314' href='#n314'>314</a>
<a id='n315' href='#n315'>315</a>
<a id='n316' href='#n316'>316</a>
<a id='n317' href='#n317'>317</a>
<a id='n318' href='#n318'>318</a>
<a id='n319' href='#n319'>319</a>
<a id='n320' href='#n320'>320</a>
<a id='n321' href='#n321'>321</a>
<a id='n322' href='#n322'>322</a>
<a id='n323' href='#n323'>323</a>
<a id='n324' href='#n324'>324</a>
<a id='n325' href='#n325'>325</a>
<a id='n326' href='#n326'>326</a>
<a id='n327' href='#n327'>327</a>
<a id='n328' href='#n328'>328</a>
<a id='n329' href='#n329'>329</a>
<a id='n330' href='#n330'>330</a>
<a id='n331' href='#n331'>331</a>
<a id='n332' href='#n332'>332</a>
<a id='n333' href='#n333'>333</a>
<a id='n334' href='#n334'>334</a>
<a id='n335' href='#n335'>335</a>
<a id='n336' href='#n336'>336</a>
<a id='n337' href='#n337'>337</a>
<a id='n338' href='#n338'>338</a>
<a id='n339' href='#n339'>339</a>
<a id='n340' href='#n340'>340</a>
<a id='n341' href='#n341'>341</a>
<a id='n342' href='#n342'>342</a>
<a id='n343' href='#n343'>343</a>
<a id='n344' href='#n344'>344</a>
<a id='n345' href='#n345'>345</a>
<a id='n346' href='#n346'>346</a>
<a id='n347' href='#n347'>347</a>
<a id='n348' href='#n348'>348</a>
<a id='n349' href='#n349'>349</a>
<a id='n350' href='#n350'>350</a>
<a id='n351' href='#n351'>351</a>
<a id='n352' href='#n352'>352</a>
<a id='n353' href='#n353'>353</a>
<a id='n354' href='#n354'>354</a>
<a id='n355' href='#n355'>355</a>
<a id='n356' href='#n356'>356</a>
<a id='n357' href='#n357'>357</a>
<a id='n358' href='#n358'>358</a>
<a id='n359' href='#n359'>359</a>
<a id='n360' href='#n360'>360</a>
<a id='n361' href='#n361'>361</a>
<a id='n362' href='#n362'>362</a>
<a id='n363' href='#n363'>363</a>
<a id='n364' href='#n364'>364</a>
<a id='n365' href='#n365'>365</a>
<a id='n366' href='#n366'>366</a>
<a id='n367' href='#n367'>367</a>
<a id='n368' href='#n368'>368</a>
<a id='n369' href='#n369'>369</a>
<a id='n370' href='#n370'>370</a>
<a id='n371' href='#n371'>371</a>
<a id='n372' href='#n372'>372</a>
<a id='n373' href='#n373'>373</a>
<a id='n374' href='#n374'>374</a>
<a id='n375' href='#n375'>375</a>
<a id='n376' href='#n376'>376</a>
<a id='n377' href='#n377'>377</a>
<a id='n378' href='#n378'>378</a>
<a id='n379' href='#n379'>379</a>
<a id='n380' href='#n380'>380</a>
<a id='n381' href='#n381'>381</a>
<a id='n382' href='#n382'>382</a>
<a id='n383' href='#n383'>383</a>
<a id='n384' href='#n384'>384</a>
<a id='n385' href='#n385'>385</a>
<a id='n386' href='#n386'>386</a>
<a id='n387' href='#n387'>387</a>
<a id='n388' href='#n388'>388</a>
<a id='n389' href='#n389'>389</a>
<a id='n390' href='#n390'>390</a>
<a id='n391' href='#n391'>391</a>
<a id='n392' href='#n392'>392</a>
<a id='n393' href='#n393'>393</a>
<a id='n394' href='#n394'>394</a>
<a id='n395' href='#n395'>395</a>
<a id='n396' href='#n396'>396</a>
<a id='n397' href='#n397'>397</a>
<a id='n398' href='#n398'>398</a>
<a id='n399' href='#n399'>399</a>
<a id='n400' href='#n400'>400</a>
<a id='n401' href='#n401'>401</a>
<a id='n402' href='#n402'>402</a>
<a id='n403' href='#n403'>403</a>
<a id='n404' href='#n404'>404</a>
<a id='n405' href='#n405'>405</a>
<a id='n406' href='#n406'>406</a>
<a id='n407' href='#n407'>407</a>
<a id='n408' href='#n408'>408</a>
<a id='n409' href='#n409'>409</a>
<a id='n410' href='#n410'>410</a>
<a id='n411' href='#n411'>411</a>
<a id='n412' href='#n412'>412</a>
<a id='n413' href='#n413'>413</a>
<a id='n414' href='#n414'>414</a>
<a id='n415' href='#n415'>415</a>
<a id='n416' href='#n416'>416</a>
<a id='n417' href='#n417'>417</a>
<a id='n418' href='#n418'>418</a>
<a id='n419' href='#n419'>419</a>
<a id='n420' href='#n420'>420</a>
<a id='n421' href='#n421'>421</a>
<a id='n422' href='#n422'>422</a>
<a id='n423' href='#n423'>423</a>
<a id='n424' href='#n424'>424</a>
<a id='n425' href='#n425'>425</a>
<a id='n426' href='#n426'>426</a>
<a id='n427' href='#n427'>427</a>
<a id='n428' href='#n428'>428</a>
<a id='n429' href='#n429'>429</a>
<a id='n430' href='#n430'>430</a>
<a id='n431' href='#n431'>431</a>
<a id='n432' href='#n432'>432</a>
<a id='n433' href='#n433'>433</a>
<a id='n434' href='#n434'>434</a>
</pre></td>
<td class='lines'><pre><code>