void ntuple_200j::MCInit( FInteger ID ) {
hbname(ID,"mcgnrl",_ievtyp,"$set");
hbname(ID,"mcevsh",_tp,"$set");
hbname(ID,"mcjets",_nxjdp,"$set");
hbname(ID,"mcpart",_ntrkp,"$set");
hbname(ID,"mcctrk",_ntrkh,"$set");
}
int main(int argc, char **argv)
{
int input;
int lrec=65536;
int status;
int cycle=1; // initialize to 1 just to avoid compiler warnings
hlimit(PAWC_SIZE);
hbset("BSIZE",&lrec,&status);
quest_.iq[9] = 256000; // extend RZ quota to 2^32 bits
hropen(50,"RZfile","bcal2nt.hbook","NQE",&lrec,&status);
hbnt(1,"BCal diagnostic ntuple"," ");
hbname(1,"bcalnt",&bcalnt,BCALNT_FORM);
hddm_s::HDDM record;
for (input=1; input<argc; input++) {
std::ifstream ifs(argv[input]);
if (!ifs.is_open()) {
std::cerr << "Error - could not open input file "
<< argv[input] << std::endl;
exit(1);
}
hddm_s::istream istr(ifs);
while (ifs.good()) {
istr >> record;
process_event(record);
record.clear();
}
}
hrout(1,&cycle," ");
hrend("RZfile");
}
void ntuple_200j::DAInit( FInteger ID ) {
hbname(ID,"dagnrl",_irun,"$set");
hbname(ID,"jdcuts",_ino,"$set");
hbname(ID,"daevsh",_tdtc,"$set");
hbname(ID,"dajets",_nxjdtc,"$set");
hbname(ID,"dactrk",_ntrk,"$set");
hbname(ID,"daclus",_nclus,"$set");
}
