aboutsummaryrefslogtreecommitdiff
path: root/src/zdab_utils.c
diff options
context:
space:
mode:
Diffstat (limited to 'src/zdab_utils.c')
-rw-r--r--src/zdab_utils.c199
1 files changed, 199 insertions, 0 deletions
diff --git a/src/zdab_utils.c b/src/zdab_utils.c
new file mode 100644
index 0000000..d946eef
--- /dev/null
+++ b/src/zdab_utils.c
@@ -0,0 +1,199 @@
+#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_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));
+}
+
+
generated by cgit (git 2.34.1) at 2025-05-06 21:22:44 +0000