add evt2root.cpp, this directly convert evt file to root, no event build
This commit is contained in:
parent
e9343e0168
commit
c93b2a76a0
2
.gitignore
vendored
2
.gitignore
vendored
|
|
@ -9,9 +9,11 @@ xia2ev2*
|
|||
pixie2root
|
||||
scan
|
||||
pxi-time-order
|
||||
evt2root
|
||||
|
||||
*.so
|
||||
*.d
|
||||
*.pcm
|
||||
|
||||
ti74pt7a
|
||||
fsu_run_2021
|
||||
|
|
|
|||
184
evt2root.cpp
Normal file
184
evt2root.cpp
Normal file
|
|
@ -0,0 +1,184 @@
|
|||
#include <stdio.h>
|
||||
#include <iostream>
|
||||
#include <fstream>
|
||||
#include <stdlib.h>
|
||||
#include <string.h>
|
||||
#include "TFile.h"
|
||||
#include "TTree.h"
|
||||
#include "TString.h"
|
||||
#include "TMath.h"
|
||||
#include "TBenchmark.h"
|
||||
#include <vector>
|
||||
|
||||
#define MAX_CRATES 2
|
||||
#define MAX_BOARDS_PER_CRATE 13
|
||||
#define MAX_CHANNELS_PER_BOARD 16
|
||||
#define BOARD_START 2
|
||||
|
||||
|
||||
class measurment{
|
||||
|
||||
public:
|
||||
UShort_t ch;
|
||||
UShort_t slot;
|
||||
UShort_t crate;
|
||||
UShort_t headerLength; /// headerLength > 4, more data except tarce.
|
||||
UShort_t eventLength; /// eventLength = headerLength + trace
|
||||
Bool_t pileup;
|
||||
ULong64_t time;
|
||||
UShort_t cfd;
|
||||
UShort_t energy;
|
||||
UShort_t trace_length;
|
||||
Bool_t trace_out_of_range;
|
||||
|
||||
UShort_t id;
|
||||
|
||||
measurment(){};
|
||||
|
||||
void Clear(){
|
||||
ch = 0;
|
||||
slot = 0;
|
||||
crate = 0;
|
||||
eventLength = 0;
|
||||
pileup = false;
|
||||
time = 0;
|
||||
cfd = 0;
|
||||
energy = 0;
|
||||
trace_length = 0;
|
||||
trace_out_of_range = 0;
|
||||
id = 0;
|
||||
}
|
||||
|
||||
void Print(){
|
||||
printf("Crate: %d, Slot: %d, Ch: %d | id: %d \n", crate, slot, ch, id);
|
||||
printf("HeaderLength: %d, Event Length: %d, energy: %d, timeStamp: %llu\n", headerLength, eventLength, energy, time);
|
||||
printf("trace_length: %d, pile-up:%d\n", trace_length, pileup);
|
||||
}
|
||||
|
||||
|
||||
};
|
||||
|
||||
int main(int argn, char **argv) {
|
||||
|
||||
if (argn != 2 && argn != 3 ) {
|
||||
printf("Usage :\n");
|
||||
printf("%s [evt File] [timeWindow] \n", argv[0]);
|
||||
printf(" timeWindow : number of tick, 1 tick = 10 ns. default = 100 \n");
|
||||
return 1;
|
||||
}
|
||||
|
||||
TString inFileName = argv[1];
|
||||
TString outFileName = inFileName;
|
||||
outFileName.Remove(inFileName.First('.'));
|
||||
outFileName.Append("_raw.root");
|
||||
|
||||
long int fprpos;
|
||||
TBenchmark gClock;
|
||||
gClock.Reset();
|
||||
gClock.Start("timer");
|
||||
|
||||
ULong64_t measureID = 0;
|
||||
|
||||
measurment data;
|
||||
|
||||
printf("====================================\n");
|
||||
|
||||
TFile * outFile = new TFile(outFileName, "recreate");
|
||||
TTree * tree = new TTree("tree", "tree");
|
||||
|
||||
tree->Branch("evID", &measureID, "data_ID/l");
|
||||
tree->Branch("detID", &data.id, "det_ID/s");
|
||||
tree->Branch("e", &data.energy, "energy/s");
|
||||
tree->Branch("t", &data.time, "time_stamp/l");
|
||||
|
||||
|
||||
FILE * inFile = fopen(inFileName, "r");
|
||||
if( inFile == NULL ){
|
||||
printf("Cannot read file : %s \n", inFileName.Data());
|
||||
return -404;
|
||||
}
|
||||
|
||||
//get file size
|
||||
fseek(inFile, 0L, SEEK_END);
|
||||
long int fprsize = ftell(inFile);
|
||||
rewind(inFile); ///back to the File begining
|
||||
|
||||
unsigned int header[4]; //read 4 header, unsigned int = 4 byte = 32 bits.
|
||||
unsigned long long nWords = 0;
|
||||
|
||||
//=============== Read File
|
||||
while ( ! feof(inFile) ){
|
||||
|
||||
fread(header, sizeof(header), 1, inFile);
|
||||
measureID ++;
|
||||
|
||||
/// see the Pixie-16 user manual, Table4-2
|
||||
data.ch = header[0] & 0xF ;
|
||||
data.slot = (header[0] >> 4) & 0xF;
|
||||
data.crate = (header[0] >> 8) & 0xF;
|
||||
data.headerLength = (header[0] >> 12) & 0x1F;
|
||||
data.eventLength = (header[0] >> 17) & 0x3FFF;
|
||||
data.pileup = header[0] >> 31 ;
|
||||
data.time = ((ULong64_t)(header[2] & 0xFFFF) << 32) + header[1];
|
||||
data.cfd = header[2] >> 16 ;
|
||||
data.energy = header[3] & 0xFF;
|
||||
data.trace_length = (header[3] >> 16) & 0x7FFF;
|
||||
data.trace_out_of_range = header[3] >> 31;
|
||||
|
||||
data.id = data.crate*MAX_BOARDS_PER_CRATE*MAX_CHANNELS_PER_BOARD + (data.slot-BOARD_START)*MAX_CHANNELS_PER_BOARD + data.ch;
|
||||
|
||||
nWords += data.eventLength;
|
||||
|
||||
///printf("----------------------%llu\n", nWords);
|
||||
///for(int i = 0; i < 4; i++){
|
||||
/// printf(" %x\n", header[i]);
|
||||
///}
|
||||
///data.Print();
|
||||
|
||||
|
||||
//=== jump to next measurement
|
||||
if( data.eventLength > 4 ){
|
||||
fseek(inFile, sizeof(int) * (data.eventLength-4), SEEK_CUR);
|
||||
}
|
||||
|
||||
//if( nWords > 200 ) break;
|
||||
|
||||
|
||||
//event stats, print status every 10000 events
|
||||
if ( measureID % 10000 == 0 ) {
|
||||
fprpos = ftell(inFile);
|
||||
float tempf = (float)fprsize/(1024.*1024.*1024.);
|
||||
gClock.Stop("timer");
|
||||
double time = gClock.GetRealTime("timer");
|
||||
gClock.Start("timer");
|
||||
printf("Total measurements: \x1B[32m%llu \x1B[0m\nPercent Complete: \x1B[32m%ld%% of %.3f GB\x1B[0m\nTime used:%3.0f min %5.2f sec\033[A\033[A\r",
|
||||
measureID, (100*fprpos/fprsize), tempf, TMath::Floor(time/60.), time - TMath::Floor(time/60.)*60.);
|
||||
}
|
||||
|
||||
//cern fill tree
|
||||
outFile->cd();
|
||||
tree->Fill();
|
||||
|
||||
}
|
||||
|
||||
fprpos = ftell(inFile);
|
||||
gClock.Stop("timer");
|
||||
double time = gClock.GetRealTime("timer");
|
||||
gClock.Start("timer");
|
||||
float tempf = (float)fprsize/(1024.*1024.*1024.);
|
||||
printf("Total measurements: \x1B[32m%llu \x1B[0m\nPercent Complete: \x1B[32m%ld%% of %.3f GB\x1B[0m\nTime used:%3.0f min %5.2f sec\033[A\r",
|
||||
measureID, (100*fprpos/fprsize), tempf, TMath::Floor(time/60.), time - TMath::Floor(time/60.)*60.);
|
||||
|
||||
fclose(inFile);
|
||||
|
||||
//cern save root
|
||||
outFile->cd();
|
||||
tree->Write();
|
||||
outFile->Close();
|
||||
|
||||
gClock.Stop("timer");
|
||||
time = gClock.GetRealTime("timer");
|
||||
printf("\n==================== finished.\r\n");
|
||||
printf("Total time spend : %3.0f min %5.2f sec\n", TMath::Floor(time/60.), time - TMath::Floor(time/60.)*60.);
|
||||
|
||||
}
|
||||
10
makefile
10
makefile
|
|
@ -1,6 +1,7 @@
|
|||
CC=g++
|
||||
|
||||
all: xia2root xia2ev2_nopart pixie2root scan pxi-time-order
|
||||
#all: xia2root xia2ev2_nopart pixie2root scan pxi-time-order
|
||||
all: xia2root xia2ev2_nopart pixie2root scan evt2root
|
||||
|
||||
xia2root: xia2root.cpp
|
||||
$(CC) xia2root.cpp -o xia2root `root-config --cflags --glibs`
|
||||
|
|
@ -11,8 +12,11 @@ xia2ev2_nopart: xia2ev2_nopart.cpp
|
|||
pixie2root: pixie2root.cpp
|
||||
$(CC) pixie2root.cpp -o pixie2root `root-config --cflags --glibs`
|
||||
|
||||
evt2root: evt2root.cpp
|
||||
$(CC) evt2root.cpp -o evt2root `root-config --cflags --glibs`
|
||||
|
||||
scan: scan.c
|
||||
$(CC) scan.c -o scan
|
||||
|
||||
pxi-time-order: pxi-time-order.c
|
||||
$(CC) pxi-time-order.c -o pxi-time-order
|
||||
#pxi-time-order: pxi-time-order.c
|
||||
# $(CC) pxi-time-order.c -o pxi-time-order
|
||||
|
|
|
|||
|
|
@ -328,11 +328,11 @@ int main(int argc, char **argv) {
|
|||
}
|
||||
else {
|
||||
tdif = subevt[sevtmult].time - etime;
|
||||
if (tdif < 0) {
|
||||
printf("SEVERE ERROR: tdiff < 0, file must be time sorted\n");
|
||||
printf("etime = %lld, time = %lld, and tdif = %lld\n", etime, subevt[sevtmult].time, tdif);
|
||||
return 0;
|
||||
}
|
||||
//if (tdif < 0) {
|
||||
// printf("SEVERE ERROR: tdiff < 0, file must be time sorted\n");
|
||||
// printf("etime = %lld, time = %lld, and tdif = %lld\n", etime, subevt[sevtmult].time, tdif);
|
||||
// return 0;
|
||||
//}
|
||||
}
|
||||
|
||||
//Check for end of event, rewind, and break out of while loop
|
||||
|
|
|
|||
Loading…
Reference in New Issue
Block a user