void
iAIDA::AIDA_HBookStore::HBook::releaseZebraDirectory(const std::string& name)
{
std::string tname = name;
int len = tname.length();
hrend_(const_cast<char*>(tname.c_str()), len);
std::string zebraDir;
if ( name[0] == '/' ) zebraDir = name;
else zebraDir = "//" + name;
char opt[1] = {' '};
rzclos_( const_cast< char*>( zebraDir.c_str() ), opt, static_cast<int>( zebraDir.size() ), 1 );
}
main(int argc, char **argv)
{
int nr,sec,layer,strip,wire,nl,ncol,nrow,i,j,l,l1,l2,ichan,nn,iev,ii;
int ind,ind1,ind2,status,status1,handle,handle1,k,m,rocid;
int scal,nw,scaler,scaler_old;
unsigned int hel,str,strob,helicity,strob_old,helicity_old, tgbi1;
int tmp1 = 1, tmp2 = 2, iret, bit1, bit2, leading[2000], trailing[2000], nleading, ntrailing;
float *bcsfl, rndm_();
char strr[1000], bankname[5], histname[128];
static int syn[32], id;
MTDCHead *mtdchead;
MTDC *mtdc;
char *HBOOKfile = "schist.his";
int nwpawc,lun,lrec,istat,icycle,idn,nbins,nbins1,igood,offset;
float x1,x2,y1,y2,ww,tmpx,tmpy,ttt,ref;
unsigned short *tlv1, *tlv2, *buf16;
printf(" dchist reached !\n");
if(argc != 2)
{
printf("Usage: dchist <fpack_filename>\n");
exit(1);
}
bcsfl = (float*)bcs_.iw;
bosInit(bcs_.iw,NBCS);
nwpawc = NWPAWC;
hlimit_(&nwpawc);
lun = 11;
lrec = LREC;
hropen_(&lun,"NTUPEL",HBOOKfile,"N",&lrec,&istat,
strlen("NTUPEL"),strlen(HBOOKfile),1);
if(istat)
{
printf("\aError: cannot open RZ file %s for writing.\n", HBOOKfile);
exit(0);
}
else
{
printf("RZ file >%s< opened for writing, istat = %d\n\n", HBOOKfile, istat);
}
nbins=200;
x1 = 0.;
x2 = 2000.;
ww = 0.;
for(k=1; k<=6; k++)
{
for(j=1; j<=2; j++)
{
for(i=1; i<=57; i++)
{
if(j==1) sprintf(histname,"ADC: sec %1d slab %2d left",k,i);
else sprintf(histname,"ADC: sec %1d slab %2d right",k,i);
idn = k*1000+j*100+i;
hbook1_(&idn,histname,&nbins,&x1,&x2,&ww,strlen(histname));
}
}
}
/*
nbins=100;
x1 = 0.;
x2 = 10000.;
nbins1=100;
y1 = 0.;
y2 = 10000.;
ww = 0.;
idn=3;
sprintf(histname,"Reference: channel1 vs channel2",idn);
hbook2_(&idn,histname,&nbins,&x1,&x2,&nbins1,&y1,&y2,&ww,strlen(histname));
idn=4;
sprintf(histname,"Reference: channel1 vs channel2 (sync)",idn);
hbook2_(&idn,histname,&nbins,&x1,&x2,&nbins1,&y1,&y2,&ww,strlen(histname));
*/
sprintf(strr,"OPEN INPUT UNIT=1 FILE='%s' ",argv[1]);
printf("fparm string: >%s<\n",strr);
status = fparm_(strr,strlen(strr));
for(iev=0; iev<1000000; iev++)
{
frbos_(bcs_.iw,&tmp1,"E",&iret,1);
if(!(iev%10000)) printf("===== Event no. %d\n",iev);
for(k=1; k<=6; k++)
{
if((ind1=bosNlink(bcs_.iw,"SC ",k)) > 0)
{
unsigned int *fb, *fbend;
unsigned short *fb16;
int crate,slot,channel,edge,data,ncol1,nrow1,adc1,adc2;
ncol1 = bcs_.iw[ind1-6];
nrow1 = bcs_.iw[ind1-5];
nw = nrow1;
offset = 0;
/*
printf("ncol=%d nrow=%d nw=%d\n",ncol1,nrow1,nw);
*/
ww = 1.0;
/*printf("\n");*/
fb16 = (unsigned short *)&bcs_.iw[ind1];
for(ii=0; ii<nrow1; ii++)
{
wire = (fb16[0])&0xFF;
adc1 = fb16[2];
adc2 = fb16[4];
/*if(k==1) printf("sec=%d wire=%d adc1=%d adc2=%d (0x%08x)\n",
k,wire,adc1,adc2,fb16);*/
if(wire>0&&wire<58)
{
tmpx = adc1;
layer = 1;
idn=k*1000+layer*100+wire;
if(tmpx>0.)
{
/*if(k==1) printf("wire=%d idn=%d,tmpx=%f\n",wire,idn,tmpx);*/
hf1_(&idn,&tmpx,&ww);
}
tmpx = adc2;
layer = 2;
idn=k*1000+layer*100+wire;
if(tmpx>0.) hf1_(&idn,&tmpx,&ww);
}
else
{
printf("ERROR: sec=%d layer=%d wire=%d (0x%08x)\n",k,layer,wire,*fb);
}
fb16 += 5;
}
}
}
/* */
if(iret == -1 || iret > 0)
{
printf(" End-of-File flag, iret =%d\n",iret);
goto a111111;
}
else if(iret < 0)
{
printf(" Error1, iret =%d\n",iret);
goto a111111;
}
if(iret != 0)
{
printf(" Error2, iret =%d\n",iret);
}
bdrop_(bcs_.iw,"E",1);
bgarb_(bcs_.iw);
}
a111111:
fparm_("CLOSE",5);
/* closing HBOOK file */
idn = 0;
hrout_(&idn,&icycle," ",1);
/*hprint_(&idn);*/
hrend_("NTUPEL", 6);
exit(0);
/* end of v1190 test */
}
void hrend(const char*filename)
{
void hrend_(const char *,int);
hrend_(filename, strlen(filename));
return;
}
void hrend(char*filename)
{
hrend_(filename, strlen(filename));
return;
}