Subversion Repositories f9daq

Compare Revisions

Ignore whitespace Rev 263 → Rev 264

/cvi/apps/petdemo/PETProjDataMgr.c
0,0 → 1,645
//
// ../bin/FBP2D FBP2D.par
//
// ********************************************************************
// * License and Disclaimer *
// * *
// * The GAMOS software is copyright of the Copyright Holders of *
// * the GAMOS Collaboration. It is provided under the terms and *
// * conditions of the GAMOS Software License, included in the file *
// * LICENSE and available at http://fismed.ciemat.es/GAMOS/license .*
// * These include a list of copyright holders. *
// * *
// * Neither the authors of this software system, nor their employing *
// * institutes,nor the agencies providing financial support for this *
// * work make any representation or warranty, express or implied, *
// * regarding this software system or assume any liability for its *
// * use. Please see the license in the file LICENSE and URL above *
// * for the full disclaimer and the limitation of liability. *
// * *
// * This code implementation is the result of the scientific and *
// * technical work of the GAMOS collaboration. *
// * By using, copying, modifying or distributing the software (or *
// * any work based on the software) you agree to acknowledge its *
// * use in resulting scientific publications, and indicate your *
// * acceptance of all terms of the GAMOS Software license. *
// ********************************************************************
//
#include "PETProjDataMgr.h"
 
 
static struct PETProjDataMgr mgr;
 
void SetNPlanes(int n) {
mgr.m_NOfPlanes=n;
};
void SetNBin(int n) {
mgr.m_NOfBins=n;
};
void SetNAng(int n) {
mgr.m_NOfAngles=n;
};
void SetDebug(int n) {
mgr.m_Debug=n;
};
void SetNumberOfChannels(int n) {
mgr.m_nch=n;
};
 
void SetRingDiameter( double x) {
mgr.m_RingDiameter = x;
};
void SetAxialDistance( double x) {
mgr.m_AxialDistance = x;
};
 
//----------------------------------------------------------------------
 
//----------------------------------------------------------------------
struct PETProjDataMgr *GetInstance() {
 
 
return &mgr;
 
}
 
//-----------------------------------------------------------------------
int PETProjDataMgrInit() {
 
/*
// ARGUMENTS:
" cout << " -------------------------- \n"
" Arguments convention: \n"
" -a Axial FOV (mm), <theDist_axial=100.0> \n"
" -d Diameter Transaxial FOV (mm), <m_RingDiameter=300.0> \n"
" -i Type of the input file (by default: 0 = Arce_binary), <typeINfile=0> \n"
" \n"
" -m Maximum ring difference (by default: -1 = m_NOfPlanes), <m_MaxRingDifferenceiff> \n"
" -n Name of output file, <m_Filename> \n"
" -p Axial number of planes, <m_NOfPlanes> \n"
" -r Number of bins, \"distancias\", <m_NOfBins> \n"
// -s Span TO DO:span !!!!!
" -t Number of angular views, \"direcciones\", <m_NOfAngles> \n"
" -v Verbosity (by default: 0=silent, 3=debug), <verbos> \n"
" -x Maximum number of coincidences to be stored (by default: -1 = no limit), <Max_Coinci> \n"
" -o Output type (by default: 0 = mcc Interfile, 1 = STIR Interfile), <OutType> \n"
 
" \n"
" PET Reconstruction. CIEMAT 2009-11 \n"
" mario.canadas@ciemat.es \n"
" -------------------------- \n";
*/
 
mgr.m_AxialDistance = (9-1)*2.25; // Axial pixel dimension*NOfPlanes
mgr.m_RingDiameter = 100.0; // notranji premer peta
mgr.m_NOfPlanes = 9; // stevilo ravnin
mgr.m_NOfBins = 128; // stevilo binov v razdalji
mgr.m_nch = 128; // stevilo padov okoli in okoli
mgr.m_NOfAngles = fabs(mgr.m_nch)/2; // stevilo kotov = stevilo padov okoli in okoli /2
 
mgr.m_MaxRingDifference = -1;
//mgr.m_MaxRingDifference = 3; // najvecja razdalja med padi
// toDo: theSpan = int(GmParameterMgr::GetInstance()->GetNumericValue("PET:ProjData:Span",1));
 
 
mgr.m_OutFormat = 1; // 1.. projections 0 .. image
 
mgr.m_Debug=1;
 
if (mgr.m_MaxRingDifference==-1) mgr.m_MaxRingDifference=mgr.m_NOfPlanes-1;
 
 
mgr.m_TotalAxialPlanes=mgr.m_NOfPlanes*mgr.m_NOfPlanes;
if (mgr.m_OutFormat==1) mgr.m_TotalAxialPlanes= (2*mgr.m_NOfPlanes-1 - mgr.m_MaxRingDifference)*mgr.m_MaxRingDifference + mgr.m_NOfPlanes; // total number of Axial planes (segments*planes) in STIR format
 
/*--- Initialize sino3D ---*/
mgr.m_projections = (SINO_TYPE ** *) malloc(mgr.m_NOfBins*sizeof(SINO_TYPE **));
for(int i=0; i<mgr.m_NOfBins; i++) {
mgr.m_projections[i] = (SINO_TYPE **) malloc(mgr.m_NOfAngles*sizeof(SINO_TYPE *));
for(int j=0; j<mgr.m_NOfAngles; j++) {
mgr.m_projections[i][j] = (SINO_TYPE *) malloc(mgr.m_TotalAxialPlanes *sizeof(SINO_TYPE)); /// ! If mgr.m_OutFormat==1 (STIR output):Matrix size depends on the MAX_Ring_Difference
for(int k=0; k<mgr.m_TotalAxialPlanes; k++) {
mgr.m_projections[i][j][k]=0;
}
}
}
 
mgr.m_TotalProjectionCoincidences=0;
mgr.m_TotalCoincidences=0;
//OutputType = "pet";
return 0;
}
 
//-----------------------------------------------------------------------
void SetProjection( int axialplane, int id) {
for(int i=0; i<mgr.m_NOfBins; i++) {
for(int j=0; j<mgr.m_NOfAngles; j++) {
mgr.m_projections[i][j][axialplane]=H2D_GetBinContent(id,i+1,j+1);
}
}
 
}
//-----------------------------------------------------------------------
 
//-----------------------------------------------------------
// from Gate
//------------------------------------------------------------
double ComputeSinogramS(double X1, double Y1, double X2, double Y2) {
double s;
 
double denom = (Y1-Y2) * (Y1-Y2) + (X2-X1) * (X2-X1);
 
if (denom!=0.) {
denom = sqrt(denom);
s = ( X1 * (Y2-Y1) + Y1 * (X1-X2) ) / denom;
} else {
s = 0.;
}
 
double theta;
if ((X1-X2)!=0.) {
theta=atan((X1-X2) /(Y1-Y2));
} else {
theta=3.1416/2.;
}
if ((theta > 0.) && ((X1-X2) > 0.)) s = -s;
if ((theta < 0.) && ((X1-X2) < 0.)) s = -s;
if ( theta < 0.) {
theta = theta+3.1416;
s = -s;
}
return s;
}
 
 
void AddEvent( const HVector3 pos1 , const HVector3 pos2) {
int z1_i, z2_i;
//for discretization on the crystal: int x1_i, x2_i, y1_i, y2_i;
 
double z1_abs=pos1.x[2]+mgr.m_AxialDistance/2;
double z2_abs=pos2.x[2]+mgr.m_AxialDistance/2;
double a, b, phi, dis;
int phi_i, dis_i;
int ring_diff;
 
double _PI=2*asin(1);
 
mgr.m_TotalCoincidences++;
 
z1_i=(int)(mgr.m_NOfPlanes* z1_abs/mgr.m_AxialDistance); //round --> mgr.m_NOfPlanes+1 ...
z2_i=(int)(mgr.m_NOfPlanes* z2_abs/mgr.m_AxialDistance);
 
// control; if z_i out of range: return
 
if ( (pos1.x[0]==pos2.x[0]) && (pos1.x[1]==pos2.x[1]) ) {
#ifndef GAMOS_NO_VERBOSE
if( mgr.m_Debug ) {
printf( "AddEvent:WARNING! Event_1 == Event_2 ; x= %f y= %f z= %f\n", pos2.x[0], pos2.x[1], pos2.x[2] );
}
#endif
return;
}
 
if ( (z1_i<0) || (z2_i<0) || (z1_i>= mgr.m_NOfPlanes) || (z2_i>= mgr.m_NOfPlanes) ) {
#ifndef GAMOS_NO_VERBOSE
if( mgr.m_Debug ) {
printf("PETProjDataMgr::AddEvent:WARNING! Event out of bounds (Axial):");
printf( "x1= %f y1= %f z1= %f ", pos1.x[0], pos1.x[1], pos1.x[2] );
printf( "x2= %f y2= %f z2= %f \n", pos2.x[0], pos2.x[1], pos2.x[2] );
}
#endif
return;
}
 
ring_diff = (int)fabs(z1_i-z2_i);
 
// max ring difference; control:
if (ring_diff > mgr.m_MaxRingDifference) {
#ifndef GAMOS_NO_VERBOSE
if( mgr.m_Debug ) {
printf("PETProjDataMgr::AddEvent:WARNING! Event out of bounds (Max. Ring Diff.): %f>%f",ring_diff , mgr.m_MaxRingDifference );
printf( "x1= %f y1= %f z1= %f ", pos1.x[0], pos1.x[1], pos1.x[2] );
printf( "x2= %f y2= %f z2= %f \n", pos2.x[0], pos2.x[1], pos2.x[2] );
}
#endif
return;
}
 
a=(double)(pos2.x[1]- pos1.x[1]);
b=(double)(pos2.x[0]- pos1.x[0]);
 
if (a==0.0) {
phi=_PI*0.5;
} else {
phi=atan(b/a);
}
 
if (phi<0) phi = phi +_PI;
 
dis=pos1.x[0]*cos(phi) - pos1.x[1]*sin(phi);
//dis=ComputeSinogramS(pos1.x(), pos1.y(), pos2.x(), pos2.x());
// control; transaxial FOV
if ( fabs(dis) > mgr.m_RingDiameter*0.5 ) {
#ifndef GAMOS_NO_VERBOSE
if( mgr.m_Debug ) {
printf("PETProjDataMgr::AddEvent:WARNING! Event out of bounds (Transaxial):" );
printf( "x1= %f y1= %f z1= %f ", pos1.x[0], pos1.x[1], pos1.x[2] );
printf( "x2= %f y2= %f z2= %f \n", pos2.x[0], pos2.x[1], pos2.x[2] );
}
#endif
return;
}
 
dis = dis + mgr.m_RingDiameter*0.5;
 
// discret values:
phi_i=RoundRealToNearestInteger( (double)(mgr.m_NOfAngles-1)*phi/_PI );
dis_i=RoundRealToNearestInteger( (double)(mgr.m_NOfBins-1)*dis/(double)mgr.m_RingDiameter );
 
if ((phi_i>=mgr.m_NOfAngles) || (dis_i>=mgr.m_NOfBins)) return; // only possible "=" because 'round' check it..
 
// OLD: (SRRB included) sino3D[dis_i][phi_i][ (z1_i+z2_i)+ring_diff*(mgr.m_NOfPlanes-1) ]++;
 
int Zpos;
 
if (mgr.m_OutFormat==0) {
Zpos = (z1_i*mgr.m_NOfPlanes + z2_i);
} else {
 
if (z1_i>=z2_i) { // SIN Max Ring_Diff: Zpos= ( ((mgr.m_NOfPlanes-ring_diff)*(mgr.m_NOfPlanes-1-ring_diff))/2 + z2_i );
 
Zpos= ( ((2*mgr.m_NOfPlanes-1 - mgr.m_MaxRingDifference - ring_diff)*(mgr.m_MaxRingDifference - ring_diff))/2 + z2_i);
 
} else {
Zpos= ( (mgr.m_TotalAxialPlanes) - ((2*mgr.m_NOfPlanes-1 - mgr.m_MaxRingDifference - ring_diff +1)*(mgr.m_MaxRingDifference - ring_diff +1))/2 + z1_i );
 
}
}
 
mgr.m_projections[dis_i][phi_i][ Zpos ]++;
mgr.m_TotalProjectionCoincidences++;
 
#ifndef GAMOS_NO_VERBOSE
if( mgr.m_Debug >1) {
printf("PETProjDataMgr::AddEvent:");
printf( "x1= %f y1= %f z1= %f ", pos1.x[0], pos1.x[1], pos1.x[2] );
printf( "x2= %f y2= %f z2= %f \n", pos2.x[0], pos2.x[1], pos2.x[2] );
printf("PETProjDataMgr::AddEvent: Sinogram pos.: distance(s)= %f angular view(phi)= %f Zpos= %f Segment (Ring diff.) =%f\n" , dis_i , phi_i, Zpos , ring_diff );
}
#endif
 
 
}
 
//-----------------------------------------------------------------------
void PETProjDataMgrFree() {
int i,j;
if (mgr.m_projections){
for(i=0; i<mgr.m_NOfBins; i++) {
for(j=0; j<mgr.m_NOfAngles; j++) {
if (mgr.m_projections[i][j]) free(mgr.m_projections[i][j]);
mgr.m_projections[i][j]= NULL;
}
if (mgr.m_projections[i]) free(mgr.m_projections[i]);
mgr.m_projections[i]= NULL;
 
}
free(mgr.m_projections);
mgr.m_projections= NULL;
}
 
}
 
 
//-----------------------------------------------------------------------
void WriteInterfile(const char * Filename ) {
 
char name_hv[512];
char name_v[512];
FILE *fp;
if (mgr.m_OutFormat==0) {
 
strcpy(name_hv, Filename);
strcpy(name_v, Filename);
strcat(name_hv, ".hv");
strcat(name_v, ".v");
fp = fopen(name_hv, "w");
 
fprintf (fp, "!INTERFILE := \n");
fprintf (fp, "name of data file := %s\n", name_v);
fprintf (fp, "!GENERAL DATA := \n");
fprintf (fp, "!GENERAL IMAGE DATA :=\n");
fprintf (fp, "!type of data := tomographic\n");
fprintf (fp, "!version of keys := 3.3\n");
fprintf (fp, "!data offset in bytes := 0\n");
fprintf (fp, "imagedata byte order := littleendian\n");
fprintf (fp, "!PET STUDY (General) :=\n");
fprintf (fp, "!PET data type := 3D-Sinogram\n");
fprintf (fp, "process status := Reconstructed\n");
fprintf (fp, "!number format := unsigned short\n");
fprintf (fp, "!number of bytes per pixel := 2\n");
fprintf (fp, "number of dimensions := 3\n");
fprintf (fp, "matrix axis label [1] := x\n");
fprintf (fp, "!matrix size [1] := %i\n",mgr.m_NOfBins);
fprintf (fp, "scaling factor (mm/pixel) [1] := %f\n",(float)(mgr.m_RingDiameter/(mgr.m_NOfBins-1.0)));
 
fprintf (fp, "matrix axis label [2] := y\n");
fprintf (fp, "!matrix size [2] := %i\n",mgr.m_NOfAngles);
 
fprintf (fp, "scaling factor (degree/pixel) [2] := %f\n",(float)(360./mgr.m_NOfAngles));
 
fprintf (fp, "matrix axis label [3] := z\n");
fprintf (fp, "!matrix size [3] := %i\n",mgr.m_NOfPlanes*mgr.m_NOfPlanes);
fprintf (fp, "scaling factor (mm/pixel) [3] := %f\n",(float)(mgr.m_AxialDistance/(mgr.m_NOfPlanes-1.0)));
 
fprintf (fp, "number of slices := %i\n",mgr.m_NOfPlanes*mgr.m_NOfPlanes);
fprintf (fp, "number of time frames := 1\n");
fprintf (fp, "image scaling factor[1] := 1\n");
fprintf (fp, "data offset in bytes[1] := 0\n");
fprintf (fp, "quantification units := 1\n");
fprintf (fp, "!END OF INTERFILE := \n");
 
fclose(fp);
//(size_t)(mgr.m_NOfBins*mgr.m_NOfAngles*mgr.m_NOfPlanes*mgr.m_NOfPlanes);
 
} else {
 
strcpy(name_hv, Filename);
strcpy(name_v, Filename);
 
strcat(name_hv, ".hs"); // STIR extension: .hs .s
strcat(name_v, ".s");
fp =fopen(name_hv, "w");
 
fprintf (fp, "!INTERFILE := \n");
fprintf (fp, "name of data file := %s\n",name_v);
fprintf (fp, "!GENERAL DATA := \n");
fprintf (fp, "!GENERAL IMAGE DATA :=\n");
fprintf (fp, "!type of data := PET\n");
// fprintf (fp, "!version of keys := 3.3\n"); STIR format is not 3.3 (almost but not completely), ERROR in STIR if it is not removed
// fprintf (fp, "!data offset in bytes := 0\n");
fprintf (fp, "imagedata byte order := littleendian\n");
fprintf (fp, "!PET STUDY (General) :=\n");
fprintf (fp, "!PET data type := Emission\n");
fprintf (fp, "applied corrections := {arc correction}\n"); // {none}\n");
// fprintf (fp, "process status := Reconstructed\n");
fprintf (fp, "!number format := unsigned integer\n");
fprintf (fp, "!number of bytes per pixel := 2\n");
 
fprintf (fp, "number of dimensions := 4\n");
fprintf (fp, "matrix axis label [4] := segment\n");
fprintf (fp, "!matrix size [4] := %i\n",mgr.m_MaxRingDifference*2 + 1);
// fprintf (fp, "scaling factor (mm/pixel) [1] := %f\n",(float)(d_FOV/(mgr.m_NOfBins-1)));
fprintf (fp, "matrix axis label [3] := axial coordinate\n");
fprintf (fp, "!matrix size [3] := { ");
if (mgr.m_MaxRingDifference==0) {
fprintf (fp, "%i}\n", mgr.m_NOfPlanes);
} else {
int m;
for(m=mgr.m_NOfPlanes-mgr.m_MaxRingDifference; m<=mgr.m_NOfPlanes; m++) fprintf (fp, "%i,", m);
for(m=mgr.m_NOfPlanes-1; m>mgr.m_NOfPlanes-mgr.m_MaxRingDifference; m--) fprintf (fp, "%i,", m);
fprintf (fp, "%i}\n", mgr.m_NOfPlanes-mgr.m_MaxRingDifference);
}
fprintf (fp, "matrix axis label [2] := view\n");
fprintf (fp, "!matrix size [2] := %i\n",mgr.m_NOfAngles);
fprintf (fp, "matrix axis label [1] := tangential coordinate\n");
fprintf (fp, "!matrix size [1] := %i\n",mgr.m_NOfBins);
 
fprintf (fp, "minimum ring difference per segment := {"); // TO DO : add SPAN (mgr.m_MaxRingDifferenceiff per seg. variable)
fprintf (fp, "%i", -mgr.m_MaxRingDifference);
int m;
for(m=-mgr.m_MaxRingDifference+1; m<=mgr.m_MaxRingDifference; m++) fprintf (fp, ",%i", m);
fprintf (fp, "}\n");
fprintf (fp, "maximum ring difference per segment := {"); // TO DO : add SPAN (mgr.m_MaxRingDifferenceiff per seg. variable)
fprintf (fp, "%i", -mgr.m_MaxRingDifference);
for(m=-mgr.m_MaxRingDifference+1; m<=mgr.m_MaxRingDifference; m++) fprintf (fp, ",%i", m);
fprintf (fp, "}\n");
 
fprintf (fp, "inner ring diameter (cm) := %f\n", mgr.m_RingDiameter/10); // STIR Required parameter, now assigned to FOV (not detectors)
fprintf (fp, "average depth of interaction (cm) := 0.0001\n");
fprintf (fp, "default bin size (cm) := %f\n",0.1*((float)mgr.m_RingDiameter/((float)mgr.m_NOfBins-1.0)) );
fprintf (fp, "number of rings := %i\n",mgr.m_NOfPlanes );
fprintf (fp, "distance between rings (cm) := %f\n", 0.1*((float)mgr.m_AxialDistance/(float)(mgr.m_NOfPlanes-1)) ); // Axial pixel dimension
 
fprintf (fp, "number of detectors per ring := %i\n",mgr.m_NOfAngles*2 );
// fprintf (fp, "number of slices := %i\n",mgr.m_NOfPlanes*mgr.m_NOfPlanes);
fprintf (fp, "number of time frames := 1\n");
fprintf (fp, "image scaling factor[1] := 1\n");
fprintf (fp, "data offset in bytes[1] := 0\n");
fprintf (fp, "quantification units := 1\n");
fprintf (fp, "!END OF INTERFILE := \n");
 
fclose(fp);
 
}
int mbsize = mgr.m_NOfBins*mgr.m_NOfAngles*mgr.m_TotalAxialPlanes*sizeof(SINO_TYPE);
mgr.m_Buffer = (SINO_TYPE *) malloc( mbsize );
 
long unsigned int cont=0;
int i,j,k;
 
for(k=0; k<mgr.m_TotalAxialPlanes; k++) {
for(j=0; j<mgr.m_NOfAngles; j++) {
for(i=0; i<mgr.m_NOfBins; i++) {
mgr.m_Buffer[cont]=mgr.m_projections[i][j][k];
cont++;
}
}
}
 
fp=fopen(name_v, "wb");
 
//printf(4096*sizeof(SINO_TYPE) );
int nb=fwrite(mgr.m_Buffer,1,mbsize, fp);
fclose(fp);
 
#ifndef GAMOS_NO_VERBOSE
printf("PETProjDataMgr::WriteInterfile: File name: %s\n", Filename );
printf("PETProjDataMgr::WriteInterfile: Numer of bytes written: %d\n" , nb );
printf("PETProjDataMgr::WriteInterfile: Planes = %d bins = %d ang_views = %d\n", mgr.m_NOfPlanes, mgr.m_NOfBins, mgr.m_NOfAngles );
printf("PETProjDataMgr::WriteInterfile: Dimensions (mm): Transaxial FOV = %f ; Axial FOV = %f ; Transaxial_pix =%f ; Plane width = %f\n", mgr.m_RingDiameter , mgr.m_AxialDistance , mgr.m_RingDiameter/(mgr.m_NOfBins-1) , mgr.m_AxialDistance/(mgr.m_NOfPlanes-1) ); // Image Axial Pixel(ssrb) == 0.5*(Plane_Width);
printf("... " );
 
printf("PETProjDataMgr::WriteInterfile: Total Coinci: %d\n" , mgr.m_TotalCoincidences );
printf("PETProjDataMgr::WriteInterfile: Sino3D Coinci: %d\n" ,mgr.m_TotalProjectionCoincidences );
#endif
 
}
 
double Mag2(HVector3 x) {
return x.x[0]*x.x[0]+x.x[1]*x.x[1]+x.x[2]*x.x[2];
}
 
 
double Mag(HVector3 x) {
return sqrt(Mag2(x));
}
 
void SetPhi(HVector3 *r,double phi) {
 
}
void SetTheta(HVector3 *r,double phi) {
 
}
 
double Product(HVector3 a, HVector3 b) {
double c=0;
for (int i=0; i<3; i++) c += a.x[i] * b.x[i];
return c;
}
 
HVector3 Multiply(double a, HVector3 b) {
HVector3 c;
for (int i=0; i<3; i++) c.x[i]= a* b.x[i];
return c;
}
 
HVector3 Add(HVector3 a, HVector3 b) {
HVector3 c;
for (int i=0; i<3; i++) c.x[i]= a.x[i] + b.x[i];
return c;
}
 
HVector3 Subtract(HVector3 a, HVector3 b) {
HVector3 c;
for (int i=0; i<3; i++) c.x[i]= a.x[i] - b.x[i];
return c;
}
 
 
HVector3 Hits2Digits(const HVector3 r) {
if (!mgr.m_nch) return r;
float smear=0.5;
 
if (mgr.m_nch<0) smear=Random(0,1);
 
double angle = atan2(r.x[0],r.x[1]); // vrne kot med -pi in pi
double twopi=2*Pi();
if (angle<0) angle+=twopi;
 
angle= ((int)(angle/twopi*fabs(mgr.m_nch))+smear)*twopi/fabs(mgr.m_nch);
//(mgr.m_rnd->Rndm()-0.5)*mgr.m_AxialDistance;
HVector3 x;
x.x[0]=sin(angle);
x.x[1]=cos(angle);
x.x[2]=0;
return x; // z coordinata ni cisto v redu
 
}
 
int FwdProject(double x,double y, double z, int nmax, int h) {
HVector3 r;
r.x[0]=x;
r.x[1]=y;
r.x[2]=z;
int h2d=1;
double tfac=mgr.m_RingDiameter*mgr.m_RingDiameter/4-Mag2(r);
double rfac= mgr.m_AxialDistance/mgr.m_RingDiameter;
for (int i=0; i<nmax; i++) {
 
double phi= Random(0,Pi());
HVector3 s;
s.x[0]=1;
s.x[1]=0;
s.x[2]=0;
SetPhi(&s,phi);
double sign = (Random(0,1)>0.5)? 1 : 0;
double theta = acos(Random(0,rfac));
theta+=sign*Pi();
 
SetTheta(&s,theta);
double t=Product(r,s);
HVector3 rx=Add(r,Multiply(-t,s));
 
double d=sqrt(t*t+tfac);
 
HVector3 r1= Add(rx,Multiply(d,s));
HVector3 r2= Add(rx,Multiply(-d,s));
 
//r1=Hits2Digits(r1);
//r2=Hits2Digits(r2);
 
HVector3 s1=Subtract(r2,r1);
double s1len= Mag(s1);
int niter=(int) (100*s1len/mgr.m_RingDiameter);
for (int j=0; j<niter; j++) {
r2=Add(r1,Multiply(Random(0,1),s1));
if (h2d) H2D_Fill(h,r2.x[0],r2.x[1],1);
else H3D_Fill(h,r2.x[0],r2.x[1],r2.x[2],1);
}
}
return 0;
}
 
int FwdProject2d(int img, int h) {
 
for (int i=0; i<H2D_GetNbinsX(img); i++) {
double x_=H2D_GetXBinCenter( img, i+1 );
for (int j=0; j<H2D_GetNbinsY(img); j++) {
double y_=H2D_GetYBinCenter(img, j+1 );
double density= H2D_GetBinContent(img,i+1,j+1);
if (density>0) FwdProject(x_,y_,mgr.m_AxialDistance*(Random(-0.5,0.5)), density,h);
}
}
return 0;
}
 
 
int FwdProject3d(int img, int h) {
 
for (int i=0; i<H3D_GetNbinsX(img); i++) {
double x_=H3D_GetXBinCenter( img, i+1 );
for (int j=0; j<H3D_GetNbinsY(img); j++) {
double y_=H3D_GetYBinCenter(img, j+1 );
for (int k=0; k<H3D_GetNbinsZ(img); k++) {
double z_=H3D_GetZBinCenter( img, k+1 );
double density= H3D_GetBinContent(img, i+1,j+1,k+1);
if (density>0) FwdProject(x_,y_,z_, density,h);
}
}
}
return 0;
}
 
 
int Phantom(int kaj) {
int img = H2D_Init(1, "img","Original Image",100,-50,50,100,-50,50);
int i=0;
// izberi sliko 0: kroglice, 1: point source 2: central ball
switch (kaj) {
 
case 0:
for (i=0; i<H2D_GetNbinsX(img); i++) {
for (int j=0; j<H2D_GetNbinsY(img); j++) {
double x_=H2D_GetXBinCenter( img, i+1 );
double y_=H2D_GetYBinCenter( img, j+1 );
double density=1000;
if ((x_*x_+y_*y_)<6) H2D_SetBinContent(img, i+1,j+1,density);
 
density=500;
if ((x_-25)*(x_-25)+y_*y_<12) H2D_SetBinContent(img,i+1,j+1,density);
density=2000;
if ((y_-25)*(y_-25)+x_*x_<2) H2D_SetBinContent(img,i+1,j+1,density);
}
}
break;
 
case 2:
for (i=0; i<H2D_GetNbinsX(img); i++) {
for (int j=0; j<H2D_GetNbinsY(img); j++) {
double x_=H2D_GetXBinCenter( img, i+1 );
double y_=H2D_GetYBinCenter( img, j+1 );
double density=1000;
if ((x_*x_+y_*y_)<12.5) H2D_SetBinContent(img, i+1,j+1,density);
}
}
break;
 
case 1:
H2D_Fill(img, 25,25,10000);
break;
 
}
 
return img;
}
/cvi/apps/petdemo/PETProjDataMgr.h
0,0 → 1,110
//
// ********************************************************************
// * License and Disclaimer *
// * *
// * The GAMOS software is copyright of the Copyright Holders of *
// * the GAMOS Collaboration. It is provided under the terms and *
// * conditions of the GAMOS Software License, included in the file *
// * LICENSE and available at http://fismed.ciemat.es/GAMOS/license .*
// * These include a list of copyright holders. *
// * *
// * Neither the authors of this software system, nor their employing *
// * institutes,nor the agencies providing financial support for this *
// * work make any representation or warranty, express or implied, *
// * regarding this software system or assume any liability for its *
// * use. Please see the license in the file LICENSE and URL above *
// * for the full disclaimer and the limitation of liability. *
// * *
// * This code implementation is the result of the scientific and *
// * technical work of the GAMOS collaboration. *
// * By using, copying, modifying or distributing the software (or *
// * any work based on the software) you agree to acknowledge its *
// * use in resulting scientific publications, and indicate your *
// * acceptance of all terms of the GAMOS Software license. *
// ********************************************************************
//
 
#ifndef PETProjDataMgr_h
#define PETProjDataMgr_h
 
/*---------------------------------------------------------------------------
ClassName: PETProjDataMgr
Author: M. Canadas, P. Arce
Changes: 01/11: creation
 
-------------------------------------------------------------------------
// Description
PET output (List-mode X Y Z) into projection data (sinograms).
 
Output data: Sinograms for PET image reconstruction. Interfile format, STIR compatible (.hs/.s),
STIR Software for Tomographic Image Reconstruction: http://stir.sourceforge.net/main.htm
 
-----------------------------------------------------------------------*/
#include <stdlib.h>
#include <stdio.h>
#include <stdint.h>
#include "H1D.h"
#include "H2D.h"
#include "H3D.h"
 
#include "daq.h"
 
 
typedef unsigned short SINO_TYPE; //!!NOTE: Check "Write_sino3D" (.hv file) if SINO_TYPE changes !!
 
struct PETProjDataMgr {
int m_Debug;
double m_AxialDistance;
double m_RingDiameter;
int m_NOfPlanes;
int m_NOfBins;
int m_NOfAngles;
int m_MaxRingDifference;
int m_OutFormat;
int m_nch;
int m_TotalAxialPlanes;
char m_Filename;
 
//G4bool bDumpCout;
SINO_TYPE ***m_projections;
SINO_TYPE *m_Buffer;
unsigned long int m_TotalCoincidences;
unsigned long int m_TotalProjectionCoincidences;
 
 
};
//-----------------------------------------------------------------------
int PETProjDataMgrInit();
void PETProjDataMgrFree();
void AddEvent(const HVector3 pos1, const HVector3 pos2);
//void ReadFile();
//PETOutput ReadEvent( G4bool& bEof );
void SetProjection(int axialplane, int h2d);
void WriteInterfile(const char *fname);
void SetNPlanes(int n);
void SetNBin(int n);
void SetNAng(int n);
void SetDebug(int n);
void SetNumberOfChannels(int n);
 
void SetRingDiameter( double x);
void SetAxialDistance( double x);
int Phantom(int kaj);
int FwdProject2d(int img2d, int h2d);
int FwdProject3d(int img3d, int h3d);
 
int FwdProject(double x,double y, double z, int nmax, int h1d);
HVector3 Hits2Digits(const HVector3 r);
 
 
struct PETProjDataMgr* GetInstance();
#endif
 
/cvi/apps/petdemo/analyse.c
0,0 → 1,351
 
#include <stdlib.h>
#include <stdio.h>
#include <stdint.h>
#include "daq.h"
#include "H1D.h"
#include "H2D.h"
#include "PETProjDataMgr.h"
 
double gAdc[MAXCH]; // raw ADC ji
float gNtdata[MAXCH*3];
int m_write=0;
 
double gData[MAXCH]; // korigirani ADC ji
 
float gSum[MAXPMT];
float gRawSum[MAXPMT];
float gMax[MAXPMT];
float gSumCluster[MAXPMT];
int gNabove[MAXPMT];
 
 
 
 
 
double GetEnergy(int ch, int adc) {
return (adc-conf.apedestals[ch])/(conf.apeak[ch]-conf.apedestals[ch])*conf.peakscaling;
}
 
 
 
int GetGlobalPosition(int ipmt, double angle, float cx, float cy, float *gposition) {
 
double phi = angle + conf.module[ipmt].phi;
double r = conf.module[ipmt].r;
double sinphi = sin(phi);
double cosphi = cos(phi);
 
gposition[0]= cosphi* r + sinphi * cx ;
gposition[1]= sinphi* r - cosphi * cx ;
gposition[2]= cy;
 
//if (m_debug) printf( "[%d] phi=%f %f global ( %f , %f , %f )\n" ,ipmt, module[ipmt].phi,angle, rv.x(), rv.y(), rv.z());
return 0;
}
 
 
int GetRealPosition(int ipmt, float px,float py,float *cx,float *cy) {
// iz CoG izracuna pozicijo na kristalu
 
int binx = 0; //####m_crystalid[ipmt]->GetXaxis()->FindBin(px);
int biny = 0; //####m_crystalid[ipmt]->GetYaxis()->FindBin(py);
int crystalid = 0;//####m_crystalid[ipmt]->GetBinContent(binx,biny);
*cx= ( crystalid % conf.nofcrystalsx - conf.nofcrystalsy/2. + 0.5 ) * conf.crystalpitchx;
*cy= ( crystalid / conf.nofcrystalsx - conf.nofcrystalsx/2. + 0.5 ) * conf.crystalpitchy;
 
 
// razmazi pozicijo po povrsini kristala
double rndx = Random(-0.5,0.5)* conf.crystalpitchx;
double rndy = Random(-0.5,0.5)* conf.crystalpitchy;
*cx+= rndx;
*cy+= rndy;
return crystalid;
}
 
 
int HistogramsInit() {
 
if (conf.write) {
char varlist[1024], tmplist[1024];
sprintf( varlist,"sum0:sum1:sum2:sum3:nh0:nh1:nh2:nh3:c0:c1:c2:c3:px0:px1:px2:px3:py0:py1:py2:py3:rx0:rx1:rx2:rx3:ry0:ry1:ry2:ry3:max0:max1:max2:max3");
for (int i=0; i<MAXCH; i++) {
sprintf(tmplist,"%s",varlist);
sprintf(varlist,"%s:a%d",tmplist,i);
}
/*
energija .. (adc -pedestal) * scaling
 
sum0 .. sum3 vsota energij po vseh kanalih fotopomnozevalke
nh0 .. nh4 stevilo zadetkov na fotopomnozevalki (adc nad thresholdom)
c0 .. c4 vsota energij za kanale nad thresholdom
px0 .. px4 x koordinata COG na fotopomnozevalki
py0 .. py4 y koordinata COG na fotopomnozevalki
max0 .. max3 maximalna energija na kanalu na fotopomnozevalki
a0 .. a63 energija na i-tem kanalu
*/
//gNtuple = new TNtuple("nt","Pet RAW data",varlist);
printf ("#%s#\n",varlist);
}
m_Nhits = H2D_Init(301,"hnabove","Stevilo zadetkov nad nivojem diskriminacije",4,-0.5,3.5, 16,+0.5,16.5 ); // stevilo hitov nad thresholdom
m_GlobalPosition =H2D_Init(302,"globalxy","Rekonstruirana koordinata zadetka",200,-80,80,200,-80,80); // reconstruirana koordinata v globalnem sistemu
 
 
char hname[0xFF];
char hn[0xFF];
int i;
for (i=0; i<MAXCH; i++) {
sprintf(hname,"Raw ADC Ch. %d ;ADC;N",i);
sprintf(hn,"ach%d",i);
//m_Adc[i] = H1D_Init(i,hn,hname,4000,-0.5,3999.5); // osnovni adcji
m_Adc[i] = H1D_Init(i,hn,hname,400,500-0.5,899.5); // osnovni adcji
sprintf(hname,"cADC za dogodke z manj kot x hitov na pmtju Ch. %d ;ADC;N",i);
sprintf(hn,"cutch%d",i);
m_AdcCut[i] = H1D_Init(i+MAXCH, hn,hname,4000,-0.5,3999.5); // adcji za zadetke z manj kot 7 hiti na PMTju
sprintf(hname,"cADC proti stevilu zadetkov na pmtju Ch. %d ;ADC;N",i);
sprintf(hn,"singlech%d",i);
m_Adc_vs_Nhits[i] = H2D_Init(i,hn,hname,200,0,4000, 17,-0.5,16.5); // adc proti stevilu zadetkov na PMTju
sprintf(hname,"cADC proti vsoti kanalov na pmtju, Ch. %d ;ADC;ADCsum",i);
sprintf(hn,"corrch%d",i);
m_Adc_vs_Sum[i] = H2D_Init(i+MAXCH,hn,hname,200,0,4000, 200,0,12000 ); // raw adc proti vsoti kanalov na PMTju
sprintf(hname,"Raw ADC proti vsoti kanalov na pmtju, Ch. %d ;ADC;ADCsum",i);
sprintf(hn,"adcvssumch%d",i);
m_Adc_vs_Sum_Uncorected[i] =H2D_Init(i+MAXCH*2,hn,hname,100,0,3500, 100,5000,12000 ); // adc proti vsoti kanalov na PMTju
}
 
for (i=0; i<MAXPMT; i++) {
sprintf(hname,"Vsota cADC na PMTju %d ;ADC;N",i);
sprintf(hn,"pmt%d",i);
m_AdcSum[i] =H1D_Init(200+i,hn,hname,200,0,20000); // vsota adcjev na PMTju
sprintf(hname,"Vsota cADC na PMTju za dogodke z manj kot x zadetki %d ;ADC;N",i);
sprintf(hn,"clusterpmt%d",i);
m_AdcSumCluster[i] = H1D_Init(216+i,hn,hname,200,0,20000); // vsota adcjev na pmtju, ki so nad threshholdom, za dogodke z manj kot 5 zadetki na PMTju
 
sprintf(hname,"Center naboja CoG PMT %d ;x;y",i);
sprintf(hn,"pmt1%d",i);
m_CenterOfGravity[i] =H2D_Init(200+i,hn,hname,200,-1.25,4.25,200,-1.25,4.25); // center naboja na PMTju
 
sprintf(hname,"Center naboja CoG za zadetke nad thresholdom PMT %d ;x;y",i);
sprintf(hn,"pmt2%d",i);
m_CenterOfGravityforChAboveThreshold[i] = H2D_Init(216+i,hn,hname,200,0,3,200,0,3); // center naboja za zadetke nad thresholdom
 
sprintf(hname,"Rekonstruirana koordinata PMT %d ;x(mm);y(mm)",i);
sprintf(hn,"pmt3%d",i);
m_ReconstructedPosition[i] = H2D_Init(232+i,hn,hname,200,-12,12,200,-12,12); // korigirana koordinata na pmtju
 
sprintf(hname,"SumAdc PMT %d ;SumAdc(a.u.);crystal ID",i);
sprintf(hn,"sumadc%d",i);
m_SumAdc[i] =H2D_Init(248+i,hn,hname,200,0,12000,81,-0.5,80.5); // SumADC spekter za razlicne kristale
}
 
// store mapping
//TH1I *hx= new TH1I("mapping", "ch mapping", 16, -0.5,15.5);
//for (int j=0;j<16;j++) hx->Fill(j,m_geo->ch[j].idx);
 
return 0;
};
 
int HistogramsWrite(const char *fname){
FILE *fp = fopen(fname, "wb");
if (fp){
int i;
for (i=0;i<H1DMAX;i++){ if (H1D_Exist(i)) H1D_Write2File(i, fp); }
for (i=0;i<H2DMAX;i++){ if (H2D_Exist(i)) H2D_Write2File(i, fp); }
for (i=0;i<H3DMAX;i++){ if (H3D_Exist(i)) H3D_Write2File(i, fp); }
fclose(fp);
}
 
return 0;
}
 
 
int HistogramsFill() {
int i=0, ipmt;
for (i=0; i<MAXCH; i++) { // zanka cez vse elektronske kanale;
float x=gData[i];
int ipmt= i/16;
float s=gSum[ipmt];
// int idx= 2*(ipmt)+64;
// const float fSlope1=2;
// const float fSlope2=3;
// if ( ((s-gSumPedestals[ipmt] ) > fSlope1*(x-gPedestals[i])) && ((s-gSumPedestals[ipmt] ) < fSlope2*(x-gPedestals[ipmt])) ) {
if (gNabove[ipmt]<7) {
H2D_Fill(m_Adc_vs_Sum[i],x,s,1);
H2D_Fill(m_Adc_vs_Sum_Uncorected[i], gAdc[i],gRawSum[ipmt],1);
H1D_Fill(m_AdcCut[i],x,1);
}
}
 
HVector3 r_coincidence[2];
for (int k=0; k<2; k++) for (int i1=0; i1<3; i1++) r_coincidence[k].x[i1]=0;
int f_coincidence[2]= {0,0};
 
for (ipmt=0; ipmt<4; ipmt++) { // zanka preko pmtjev
int j2= (ipmt/2)*2+1-ipmt%2; // sosednja fotopomnozevalka
 
float posx[2]= {0,0};
float posy[2]= {0,0};
float sum[2]= {0,0};
for (int ich=0; ich<16; ich++) { // zanka preko elektronskih kanalov na fotopomnozevalki
int ch= ich+ipmt*16;
if (gMax[ipmt]>conf.adcthreshold) {
posx[0]+= gData[ch]*conf.channel[ich].ix;
posy[0]+= gData[ch]*conf.channel[ich].iy;
sum[0] += gData[ch];
 
if (gData[ch]> 0.1*gMax[ipmt]) { // pri racunanju pozicije upostevaj le kanale, ki imajo vrednost vecjo od ratio*maksimalna na tisti fotopomnozevalki ********
posx[1]+= gData[ch]*conf.channel[ich].ix;
posy[1]+= gData[ch]*conf.channel[ich].iy;
sum[1] += gData[ch];
}
}
}
 
if (m_write) {
gNtdata[12+ipmt]=posx[0]/sum[0];
gNtdata[16+ipmt]=posy[0]/sum[0];
gNtdata[20+ipmt]=(sum[1]>0)?posx[1]/sum[1]:-10;
gNtdata[24+ipmt]=(sum[1]>0)?posy[1]/sum[1]:-10;
gNtdata[28+ipmt]=sum[0];
}
if (gMax[ipmt]>gMax[j2]) { // izberi fotopomnozevalko, ki ima vec zadetkov kot soseda TRG= (P1|P2) & (P3|P4)
 
if ( sum[0] > 0 ) {
float px=posx[0]/sum[0];
float py=posy[0]/sum[0];
 
H2D_Fill(m_CenterOfGravity[ipmt],px,py,1);
 
float cx=0; //koordinata na pmtju
float cy=0;
int crystalID = GetRealPosition(ipmt,px,py,&cx,&cy);
H2D_Fill(m_ReconstructedPosition[ipmt],cx,cy,1);
H2D_Fill(m_SumAdc[ipmt],gSum[ipmt], crystalID,1 );
if (debug > 2) printf("ipmt=%d cx=%f cy=%f\n",ipmt,cx,cy);
GetGlobalPosition(ipmt, conf.rotation, cx, cy,r_coincidence[ipmt/2].x) ;
f_coincidence[ipmt/2] =1;
H2D_Fill(m_GlobalPosition,r_coincidence[ipmt/2].x[0],r_coincidence[ipmt/2].x[1],1);
 
}
if ( sum[1] > 0 ) {
 
float px=posx[1]/sum[1];
float py=posy[1]/sum[1];
// float px0=posx[0]/sum[0];
// float py0=posy[0]/sum[0];
// if (px0<0.25 || px0>2.75 || py0< 0.25 || py0> 2.75)
// if (((py0>0.75 && py0<1.25) || (py0>1.75 && py0<2.25)) && (px0< 0.25 || px0>2.75 ) )
H2D_Fill(m_CenterOfGravityforChAboveThreshold[ipmt],px,py,1);
}
 
}
}
 
for (i=0; i<MAXCH; i++) if (gData[i]>conf.adcthreshold) {
H2D_Fill(m_Adc_vs_Nhits[i],gData[i],gNabove[i/16],1);
}
for (i=0; i<MAXCH; i++) if (gData[i]>conf.adcthreshold && gNabove[i/16]<5 ) {
gSumCluster[i/16]+=gData[i];
}
for (ipmt=0; ipmt<4; ipmt++) {
H1D_Fill(m_AdcSum[ipmt],gSum[ipmt],1);
H1D_Fill(m_AdcSumCluster[ipmt],gSumCluster[ipmt],1);
}
 
for (i=0; i<4; i++) H1D_Fill(m_Nhits,i,gNabove[i]);
 
if (m_write) {
for (i=0; i<4; i++) {
gNtdata[i]=gSum[i];
gNtdata[4+i]=gNabove[i];
gNtdata[8+i]=gSumCluster[i];
}
for (i=0; i<MAXCH; i++) gNtdata[32+i]=gData[i];
//gNtuple->Fill(gNtdata);
}
//
// coincidences
//
if (f_coincidence[0]&&f_coincidence[1]) {
conf.coincidences++;
if (debug > 1) {
printf( "Coincidence ( %f , %f , %f ) ( %f , %f , %f )\n" ,
r_coincidence[0].x[0], r_coincidence[0].x[1], r_coincidence[0].x[2],
r_coincidence[1].x[0], r_coincidence[1].x[1], r_coincidence[1].x[2] );
}
AddEvent( r_coincidence[0] , r_coincidence[1]);
}
return 0;
};
 
 
 
int analyse(int nb, uint32_t *buf) {
 
int idx=1;
int neve=buf[0]/2;
//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
for (int nev=0; nev<neve; nev++) {
int len=buf[idx];
int sb =buf[idx+1];
unsigned int *pbuf=&buf[idx+2];
if (sb!=0xffab) {
printf("0x%04x!0xffab len=%d\n",sb,len);
break;
}
// postavi na nic
#define InitArrayWithAValue(arr,n,x) {for (int i=0;i<n;i++) arr[i]=x;}
 
//------------------------------------------------------------
// postavi na nic
 
InitArrayWithAValue( gSum , MAXPMT , 0);
InitArrayWithAValue( gRawSum , MAXPMT , 0);
InitArrayWithAValue( gMax , MAXPMT , 0);
InitArrayWithAValue( gSumCluster, MAXPMT , 0);
InitArrayWithAValue( gNabove , MAXPMT , 0);
InitArrayWithAValue( gData , MAXCH, 0);
InitArrayWithAValue( gAdc , MAXCH, 0);
for (int i0=0; i0<len-2; i0+=2) {
 
int data0 = pbuf[i0];
int data1 = pbuf[i0+1];
int geo = (data1 >> 11) & 0x1f;
int ch = (data1&0x1f) | (geo<<5);
int dtype = (data1>>9)&0x3;
int adc = data0&0xfff;
if (dtype == 0 & ch<MAXCH) {
if (debug) printf("%d ADC %d %d\n", geo, ch, adc);
if (ch<MAXCH) H1D_Fill(ch, adc,1);
int ipmt = ch/16;
 
gAdc[ch]=adc;
gRawSum[ipmt]+=adc;
 
gData[ch]=GetEnergy(ch,adc);
 
gSum[ipmt]+=gData[ch];
if (gData[ch] >gMax[ipmt] ) gMax[ipmt]= gData[ch];
if (gData[ch] >conf.adcthreshold ) gNabove[ipmt]++;
//#### m_Adc[ch]->Fill(adc);
}
};// for (int i0=0;i0<len-2;i0+=2)
//------------------------------------------------------------
 
idx+=len+1;
HistogramsFill();
if (m_write) {
gNtdata[0]=0;
for (int i=0; i<MAXCH; i++) gNtdata[i+1]=gAdc[i];
//####gNtuple->Fill(gNtdata);
}
 
} // for (int nev=0;nev<neve;nev++)
 
 
 
 
return 0;
}
/cvi/apps/petdemo/config.ini
0,0 → 1,20
[scintillator]
 
nofcrystalsx=9
nofcrystalsy=9
crystalpitchx=2.25
crystalpitchy=2.25
 
 
[sensor]
modules=ini/modules.map
channels=ini/m16.map
 
 
[calibration]
 
sumpedestals=ini/sumpedestals.dat
pedestals=ini/pedestals_new.dat
photopeak=ini/fotovrh.dat
channelcalibration=ini/calibration_all.root
adcthreshold=300
/cvi/apps/petdemo/daq.h
0,0 → 1,117
#ifndef _daq_h_
#define _daq_h_
 
#define MAXCH 64
#define MAXPMT 16
 
#define NSTEPS 4790 // en krog je 4790 korakov
 
typedef struct {
float x[3];
} HVector3;
 
struct Module {
double r;
double phi;
};
 
 
 
struct Channel {
int ix;
int iy;
};
 
 
struct Config {
int nofcrystalsx;
int nofcrystalsy;
double crystalpitchx;
double crystalpitchy;
char *modules;
char *channels;
char *sumpedestals;
char *pedestals;
char *photopeak;
char *channelcalibration;
int adcthreshold;
double peakscaling;
struct Module module[16];
struct Channel channel[16];
float apedestals[MAXCH];
float apeak[MAXCH];
float rotation;
int write;
int coincidences; // stevilo koincidenc
} conf;
 
 
 
#define NUM_CHANNELS 1
#define PEDESTAL 255
#define SCAN_TYPE 1
 
#define RUNREC_ID 1
#define ENDREC_ID 2
#define POSREC_ID 3
#define EVTREC_ID 4
 
typedef struct {
unsigned int id,length;
unsigned int fver,time;
unsigned int num_events,num_channels,pedestal,xy;
int nx,x0,dx,ny,y0,dy;
} RUNREC;
typedef struct {
unsigned int id,length;
unsigned int time,dum;
} ENDREC;
 
typedef struct {
unsigned int id,length;
unsigned int time;
int num_iter_x,mikro_pos_x,set_pos_x;
int num_iter_y,mikro_pos_y,set_pos_y;
} POSREC;
 
int debug;
int fPedestal;
int fThresholdEnable;
int fStop;
int fMode;
int clear();
int end();
int event(unsigned int *, int,int *, int );
int init();
int vmconnect();
int vmdisconnect();
int vmacquire (int neve,FILE *fp, char *fpedname);
 
 
int HistogramsInit();
int HistogramsFill();
int HistogramsWrite(const char *fname);
int GetGlobalPosition(int ipmt, double angle, float cx, float cy, float *gposition);
int analyse(int nb, uint32_t *data);
 
 
//-----------------
 
int m_Adc[MAXCH];
int m_AdcCut[MAXCH];
int m_Adc_vs_Nhits[MAXCH];
int m_Adc_vs_Sum[MAXCH];
int m_Adc_vs_Sum_Uncorected[MAXCH];
int m_Nhits;
int m_AdcSum[MAXPMT];
int m_CenterOfGravity[MAXPMT];
int m_ReconstructedPosition[MAXPMT];
int m_CenterOfGravityforChAboveThreshold[MAXPMT];
int m_GlobalPosition;
int m_AdcSumCluster[MAXPMT];
int m_MaxAdc[MAXPMT];
int m_SumAdc[MAXPMT];
 
//-----------------
 
#endif
/cvi/apps/petdemo/ini/calibration_all.root
Cannot display: file marked as a binary type.
svn:mime-type = application/octet-stream
Property changes:
Added: svn:mime-type
## -0,0 +1 ##
+application/octet-stream
\ No newline at end of property
Index: apps/petdemo/ini/config.xml
===================================================================
--- apps/petdemo/ini/config.xml (nonexistent)
+++ apps/petdemo/ini/config.xml (revision 264)
@@ -0,0 +1,24 @@
+<?xml version="1.0"?>
+<pet>
+
+<scintillator>
+<nofcrystalsx>9</nofcrystalsx>
+<nofcrystalsy>9</nofcrystalsy>
+<crystalpitchx>2.25</crystalpitchx>
+<crystalpitchy>2.25</crystalpitchy>
+</scintillator>
+
+<sensor>
+<modules>ini/modules.map</modules>
+<channels>ini/m16.map</channels>
+</sensor>
+
+<calibration>
+<sumpedestals>ini/sumpedestals.dat</sumpedestals>
+<pedestals>ini/pedestals_new.dat</pedestals>
+<photopeak>ini/fotovrh.dat</photopeak>
+<channelcalibration>ini/calibration_all.root</channelcalibration>
+<adcthreshold>300</adcthreshold>
+</calibration>
+
+</pet>
Index: apps/petdemo/ini/fotovrh.dat
===================================================================
--- apps/petdemo/ini/fotovrh.dat (nonexistent)
+++ apps/petdemo/ini/fotovrh.dat (revision 264)
@@ -0,0 +1,64 @@
+0 3196.160841
+1 1711.854801
+2 1635.564402
+3 1597.419203
+4 1711.854801
+5 1584.704136
+6 1355.832941
+7 1432.123339
+8 1533.843871
+9 1648.279469
+10 1597.419203
+11 1521.128804
+12 1470.268539
+13 1686.424668
+14 1622.849336
+15 1660.994535
+16 2522.262320
+17 2220.457458
+18 2423.898521
+19 2334.893056
+20 2551.049185
+21 2449.328654
+22 2703.629982
+23 2614.624517
+24 2487.473853
+25 2347.608122
+26 2156.882126
+27 2233.172524
+28 2589.194384
+29 2233.172524
+30 1864.435598
+31 1915.295864
+32 2852.854047
+33 1915.295864
+34 1915.295864
+35 2156.882126
+36 1991.586262
+37 1966.156129
+38 1915.295864
+39 1902.580797
+40 1940.725996
+41 1877.150664
+42 1940.725996
+43 1877.150664
+44 1928.010930
+45 1940.725996
+46 1991.586262
+47 1928.010930
+48 2852.854047
+49 2080.591727
+50 1813.575332
+51 1978.871196
+52 2245.887591
+53 2080.591727
+54 1813.575332
+55 1940.725996
+56 1902.580797
+57 1889.865731
+58 1673.709601
+59 1699.139734
+60 1813.575332
+61 1877.150664
+62 1508.413738
+63 1686.424668
Index: apps/petdemo/ini/m16.map
===================================================================
--- apps/petdemo/ini/m16.map (nonexistent)
+++ apps/petdemo/ini/m16.map (revision 264)
@@ -0,0 +1,17 @@
+#CH IX IY
+15 0 0
+14 0 1
+13 1 0
+12 1 1
+11 0 3
+10 0 2
+9 1 3
+8 1 2
+7 2 3
+6 2 2
+5 3 3
+4 3 2
+3 2 0
+2 2 1
+1 3 0
+0 3 1
Index: apps/petdemo/ini/m16_reverse.map
===================================================================
--- apps/petdemo/ini/m16_reverse.map (nonexistent)
+++ apps/petdemo/ini/m16_reverse.map (revision 264)
@@ -0,0 +1,17 @@
+#CH IX IY
+0 0 0
+1 0 1
+2 1 0
+3 1 1
+4 0 3
+5 0 2
+6 1 3
+7 1 2
+8 2 3
+9 2 2
+10 3 3
+11 3 2
+12 2 0
+13 2 1
+14 3 0
+15 3 1
Index: apps/petdemo/ini/modules.map
===================================================================
--- apps/petdemo/ini/modules.map (nonexistent)
+++ apps/petdemo/ini/modules.map (revision 264)
@@ -0,0 +1,5 @@
+#module ID (connector) r(mm) angle(deg)
+0 61 0
+1 61 22.5
+2 61 180
+3 61 202.5
Index: apps/petdemo/ini/pedestals.dat
===================================================================
--- apps/petdemo/ini/pedestals.dat (nonexistent)
+++ apps/petdemo/ini/pedestals.dat (revision 264)
@@ -0,0 +1,64 @@
+0 339.764
+1 356.526
+2 328.481
+3 283.352
+4 307.915
+5 365.996
+6 316.9
+7 316.175
+8 348.054
+9 338.349
+10 385.901
+11 342.318
+12 353.552
+13 348.903
+14 340.133
+15 335.819
+16 337.356
+17 334.646
+18 322.328
+19 320.371
+20 304.205
+21 347.213
+22 336.6
+23 325.069
+24 351.203
+25 347.073
+26 307.865
+27 285.498
+28 332.387
+29 325.668
+30 352.567
+31 343.887
+32 360.176
+33 328.103
+34 337.231
+35 324.995
+36 294.699
+37 354.216
+38 304.069
+39 369.572
+40 375.065
+41 347.03
+42 341.669
+43 311.298
+44 338.887
+45 322.734
+46 388.779
+47 320.116
+48 323.118
+49 350.691
+50 344.091
+51 334.575
+52 349.252
+53 351.313
+54 338.096
+55 330.029
+56 348.157
+57 304.036
+58 300.07
+59 304.476
+60 307.644
+61 296.693
+62 327.384
+63 345.163
Index: apps/petdemo/ini/pedestals_new.dat
===================================================================
--- apps/petdemo/ini/pedestals_new.dat (nonexistent)
+++ apps/petdemo/ini/pedestals_new.dat (revision 264)
@@ -0,0 +1,64 @@
+0 339.13
+1 354.26
+2 327.546
+3 282.69
+4 312.863
+5 365.795
+6 317.055
+7 316.198
+8 348.804
+9 337.947
+10 385.535
+11 341.414
+12 351.877
+13 346.92
+14 339.204
+15 334.106
+16 331.704
+17 333.163
+18 302.789
+19 319.067
+20 325.209
+21 348.417
+22 337.329
+23 326.671
+24 351.699
+25 348.799
+26 308.549
+27 286.519
+28 332.446
+29 324.031
+30 351.92
+31 342.048
+32 359.25
+33 325.419
+34 334.892
+35 322.25
+36 293.95
+37 354.971
+38 304.267
+39 370.495
+40 374.192
+41 347.656
+42 340.723
+43 311.266
+44 301.842
+45 320.301
+46 387.6
+47 318.137
+48 321.878
+49 348.859
+50 342.913
+51 332.14
+52 346.995
+53 350.787
+54 337.321
+55 329.6
+56 346.56
+57 303.365
+58 298.587
+59 302.891
+60 306.159
+61 293.716
+62 324.878
+63 341.973
Index: apps/petdemo/ini/ph511.dat
===================================================================
--- apps/petdemo/ini/ph511.dat (nonexistent)
+++ apps/petdemo/ini/ph511.dat (revision 264)
@@ -0,0 +1,84 @@
+0 1700
+1 1900
+2 1680
+3 1700
+4 1870
+5 1680
+6 1430
+7 1460
+8 1680
+9 1700
+10 1713
+11 1650
+12 1630
+13 1830
+14 1750
+15 1740
+16 2440
+17 2320
+18 2640
+19 2500
+20 2790
+21 2600
+22 2900
+23 2820
+24 2660
+25 2470
+26 2310
+27 2330
+28 2817
+29 2860
+30 2140
+31 2140
+32 2120
+33 2410
+34 1990
+35 2480
+36 2200
+37 2100
+38 2040
+39 2060
+40 2090
+41 1970
+42 2032
+43 1970
+44 2054
+45 2120
+46 2076
+47 1920
+48 2206
+49 2420
+50 1930
+51 2290
+52 2400
+53 2425
+54 1950
+55 2400
+56 2076
+57 2010
+58 1792
+59 1830
+60 2000
+61 2032
+62 1640
+63 1614
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
Index: apps/petdemo/ini/sumpedestals.dat
===================================================================
--- apps/petdemo/ini/sumpedestals.dat (nonexistent)
+++ apps/petdemo/ini/sumpedestals.dat (revision 264)
@@ -0,0 +1,4 @@
+0 5413.22
+1 5264.29
+2 5394.42
+3 5261.11
Index: apps/petdemo/libxxusb.cpp
===================================================================
--- apps/petdemo/libxxusb.cpp (nonexistent)
+++ apps/petdemo/libxxusb.cpp (revision 264)
@@ -0,0 +1,1651 @@
+
+// libxxusb.cpp : Defines the entry point for the DLL application.
+//
+
+
+
+#include <string.h>
+#include <malloc.h>
+#include "usb.h"
+#include "libxxusb.h"
+#include <time.h>
+
+
+
+
+// 03/09/06 Release 3.00 changes
+// 07/28/06 correction CAMAC write for F to be in range 16...23
+// 10/09/06 correction CAMAC read for F to be in range <16 OR >23
+// 10/16/06 CAMAC DGG corrected
+// 12/28/07 Open corrected for bug when calling register after opening
+/*
+******** xxusb_longstack_execute ************************
+
+ Executes stack array passed to the function and returns the data read from the VME bus
+
+ Paramters:
+ hdev: USB device handle returned from an open function
+ DataBuffer: pointer to the dual use buffer
+ when calling , DataBuffer contains (unsigned short) stack data, with first word serving
+ as a placeholder
+ upon successful return, DataBuffer contains (unsigned short) VME data
+ lDataLen: The number of bytes to be fetched from VME bus - not less than the actual number
+ expected, or the function will return -5 code. For stack consisting only of write operations,
+ lDataLen may be set to 1.
+ timeout: The time in ms that should be spent tryimg to write data.
+
+ Returns:
+ When Successful, the number of bytes read from xxusb.
+ Upon failure, a negative number
+
+ Note:
+ The function must pass a pointer to an array of unsigned integer stack data, in which the first word
+ is left empty to serve as a placeholder.
+ The function is intended for executing long stacks, up to 4 MBytes long, both "write" and "read"
+ oriented, such as using multi-block transfer operations.
+ Structure upon call:
+ DataBuffer(0) = 0(don't care place holder)
+ DataBuffer(1) = (unsigned short)StackLength bits 0-15
+ DataBuffer(2) = (unsigned short)StackLength bits 16-20
+ DataBuffer(3 - StackLength +2) (unsigned short) stack data
+ StackLength represents the number of words following DataBuffer(1) word, thus the total number
+ of words is StackLength+2
+ Structure upon return:
+ DataBuffer(0 - (ReturnValue/2-1)) - (unsigned short)array of returned data when ReturnValue>0
+*/
+
+int xxusb_longstack_execute(usb_dev_handle *hDev, void *DataBuffer, int lDataLen, int timeout)
+{
+ int ret;
+ char *cbuf;
+ unsigned short *usbuf;
+ int bufsize;
+
+ cbuf = (char *)DataBuffer;
+ usbuf = (unsigned short *)DataBuffer;
+ cbuf[0]=12;
+ cbuf[1]=0;
+ bufsize = 2*(usbuf[1]+0x10000*usbuf[2])+4;
+ ret=usb_bulk_write(hDev, XXUSB_ENDPOINT_OUT, cbuf, bufsize, timeout);
+ if (ret>0)
+ ret=usb_bulk_read(hDev, XXUSB_ENDPOINT_IN, cbuf, lDataLen, timeout);
+
+ return ret;
+}
+
+/*
+******** xxusb_bulk_read ************************
+
+ Reads the content of the usbfifo whenever "FIFO full" flag is set,
+ otherwise times out.
+
+ Paramters:
+ hdev: USB device handle returned from an open function
+ DataBuffer: pointer to an array to store data that is read from the VME bus;
+ the array may be declared as byte, unsigned short, or unsigned long
+ lDatalen: The number of bytes to read from xxusb
+ timeout: The time in ms that should be spent waiting for data.
+
+ Returns:
+ When Successful, the number of bytes read from xxusb.
+ Upon failure, a negative number
+
+ Note:
+ Depending upon the actual need, the function may be used to return the data in a form
+ of an array of bytes, unsigned short integers (16 bits), or unsigned long integers (32 bits).
+ The latter option of passing a pointer to an array of unsigned long integers is meaningful when
+ xxusb data buffering option is used (bit 7=128 of the global register) that requires data
+ 32-bit data alignment.
+
+*/
+int xxusb_bulk_read(usb_dev_handle *hDev, void *DataBuffer, int lDataLen, int timeout)
+{
+int ret;
+char *cbuf;
+cbuf = (char *)DataBuffer;
+ ret = usb_bulk_read(hDev, XXUSB_ENDPOINT_IN, cbuf, lDataLen, timeout);
+ return ret;
+}
+
+/*
+******** xxusb_bulk_write ************************
+
+ Writes the content of an array of bytes, unsigned short integers, or unsigned long integers
+ to the USB port fifo; times out when the USB fifo is full (e.g., when xxusb is busy).
+
+ Paramters:
+ hdev: USB device handle returned from an open function
+ DataBuffer: pointer to an array storing the data to be sent;
+ the array may be declared as byte, unsigned short, or unsigned long
+ lDatalen: The number of bytes to to send to xxusb
+ timeout: The time in ms that should be spent waiting for data.
+
+ Returns:
+ When Successful, the number of bytes passed to xxusb.
+ Upon failure, a negative number
+
+ Note:
+ Depending upon the actual need, the function may be used to pass to xxusb the data in a form
+ of an array of bytes, unsigned short integers (16 bits), or unsigned long integers (32 bits).
+*/
+int xxusb_bulk_write(usb_dev_handle *hDev, void *DataBuffer, int lDataLen, int timeout)
+{
+int ret;
+char *cbuf;
+cbuf = (char *)DataBuffer;
+ ret = usb_bulk_write(hDev, XXUSB_ENDPOINT_OUT, cbuf, lDataLen, timeout);
+ return ret;
+}
+
+/*
+******** xxusb_usbfifo_read ************************
+
+ Reads data stored in the xxusb fifo and packs them in an array of long integers.
+
+ Paramters:
+ hdev: USB device handle returned from an open function
+ DataBuffer: pointer to an array of long to store data that is read
+ the data occupy only the least significant 16 bits of the 32-bit data words
+ lDatalen: The number of bytes to read from the xxusb
+ timeout: The time in ms that should be spent waiting for data.
+
+ Returns:
+ When Successful, the number of bytes read from xxusb.
+ Upon failure, a negative number
+
+ Note:
+ The function is not economical as it wastes half of the space required for storing
+ the data received. Also, it is relatively slow, as it performs extensive data repacking.
+ It is recommended to use xxusb_bulk_read with a pointer to an array of unsigned short
+ integers.
+*/
+int xxusb_usbfifo_read(usb_dev_handle *hDev, int *DataBuffer, int lDataLen, int timeout)
+{
+int ret;
+char *cbuf;
+unsigned short *usbuf;
+int i;
+
+cbuf = (char *)DataBuffer;
+usbuf = (unsigned short *)DataBuffer;
+
+ ret = usb_bulk_read(hDev, XXUSB_ENDPOINT_IN, cbuf, lDataLen, timeout);
+ if (ret > 0)
+ for (i=ret/2-1; i >= 0; i=i-1)
+ {
+ usbuf[i*2]=usbuf[i];
+ usbuf[i*2+1]=0;
+ }
+ return ret;
+}
+
+
+//******************************************************//
+//******************* GENERAL XX_USB *******************//
+//******************************************************//
+// The following are functions used for both VM_USB & CC_USB
+
+
+/*
+******** xxusb_register_write ************************
+
+ Writes Data to the xxusb register selected by RedAddr. For
+ acceptable values for RegData and RegAddr see the manual
+ the module you are using.
+
+ Parameters:
+ hdev: usb device handle returned from open device
+ RegAddr: The internal address if the xxusb
+ RegData: The Data to be written to the register
+
+ Returns:
+ Number of bytes sent to xxusb if successful
+ 0 if the register is write only
+ Negative numbers if the call fails
+*/
+short xxusb_register_write(usb_dev_handle *hDev, short RegAddr, long RegData)
+{
+ long RegD;
+ char buf[8]={5,0,0,0,0,0,0,0};
+ int ret;
+ int lDataLen;
+ int timeout;
+ if ((RegAddr==0) || (RegAddr==12) || (RegAddr==15))
+ return 0;
+ buf[2]=(char)(RegAddr & 15);
+ buf[4]=(char)(RegData & 255);
+
+ RegD = RegData >> 8;
+ buf[5]=(char)(RegD & 255);
+ RegD = RegD >>8;
+ if (RegAddr==8)
+ {
+ buf[6]=(char)(RegD & 255);
+ lDataLen=8;
+ }
+ else
+ lDataLen=6;
+ timeout=10;
+ ret=xxusb_bulk_write(hDev, buf, lDataLen, timeout);
+ return ret;
+}
+
+/*
+******** xxusb_stack_write ************************
+
+ Writes a stack of VME/CAMAC calls to the VM_USB/CC_USB
+ to be executed upon trigger.
+
+ Parameters:
+ hdev: usb device handle returned from an open function
+ StackAddr: internal register to which the stack should be written
+ lpStackData: Pointer to an array holding the stack
+
+ Returns:
+ The number of Bytes written to the xxusb when successful
+ A negative number upon failure
+*/
+short xxusb_stack_write(usb_dev_handle *hDev, short StackAddr, long *intbuf)
+{
+ int timeout;
+ short ret;
+ short lDataLen;
+ char buf[2000];
+ short i;
+ int bufsize;
+
+ buf[0]=(char)((StackAddr & 51) + 4);
+ buf[1]=0;
+ lDataLen=(short)(intbuf[0] & 0xFFF);
+ buf[2]=(char)(lDataLen & 255);
+ lDataLen = lDataLen >> 8;
+ buf[3] = (char)(lDataLen & 255);
+ bufsize=intbuf[0]*2+4;
+ if (intbuf[0]==0)
+ return 0;
+ for (i=1; i <= intbuf[0]; i++)
+ {
+ buf[2+2*i] = (char)(intbuf[i] & 255);
+ buf[3+2*i] = (char)((intbuf[i] >>8) & 255);
+ }
+ timeout=50;
+ ret=usb_bulk_write(hDev, XXUSB_ENDPOINT_OUT, buf, bufsize, timeout);
+ return ret;
+}
+
+/*
+******** xxusb_stack_execute **********************
+
+ Writes, executes and returns the value of a DAQ stack.
+
+ Parameters:
+ hdev: USB device handle returned from an open function
+ intbuf: Pointer to an array holding the values stack. Upon return
+ Pointer value is the Data returned from the stack.
+
+ Returns:
+ When successful, the number of Bytes read from xxusb
+ Upon Failure, a negative number.
+*/
+short xxusb_stack_execute(usb_dev_handle *hDev, long *intbuf)
+{
+ int timeout;
+ short ret;
+ short lDataLen;
+ char buf[26700];
+ short i;
+ int bufsize;
+ int ii = 0;
+
+ buf[0]=12;
+ buf[1]=0;
+ lDataLen=(short)(intbuf[0] & 0xFFF);
+ buf[2]=(char)(lDataLen & 255);
+ lDataLen = lDataLen >> 8;
+ buf[3] = (char)(lDataLen & 15);
+ bufsize=intbuf[0]*2+4;
+ if (intbuf[0]==0)
+ return 0;
+ for (i=1; i <= intbuf[0]; i++)
+ {
+ buf[2+2*i] = (char)(intbuf[i] & 255);
+ buf[3+2*i] = (char)((intbuf[i] >>8) & 255);
+ }
+ timeout=2000;
+ ret=usb_bulk_write(hDev, XXUSB_ENDPOINT_OUT, buf, bufsize, timeout);
+ if (ret>0)
+ {
+ lDataLen=26700;
+ timeout=6000;
+ ret=usb_bulk_read(hDev, XXUSB_ENDPOINT_IN, buf, lDataLen, timeout);
+ if (ret>0)
+ for (i=0; i < ret; i=i+2)
+ intbuf[ii++]=(UCHAR)(buf[i]) +(UCHAR)( buf[i+1])*256;
+ }
+ return ret;
+}
+
+/*
+******** xxusb_stack_read ************************
+
+ Reads the current DAQ stack stored by xxusb
+
+ Parameters:
+ hdev: USB device handle returned by an open function
+ StackAddr: Indicates which stack to read, primary or secondary
+ intbuf: Pointer to a array where the stack can be stored
+
+ Returns:
+ Number of bytes read from xxusb when successful
+ Upon failure, a negative number
+*/
+short xxusb_stack_read(usb_dev_handle *hDev, short StackAddr, long *intbuf)
+{
+ int timeout;
+ short ret;
+ short lDataLen;
+ short bufsize;
+ char buf[1600];
+ int i;
+
+ buf[0]=(char)(StackAddr & 51);
+ buf[1]=0;
+ lDataLen = 2;
+ timeout=100;
+ ret=usb_bulk_write(hDev, XXUSB_ENDPOINT_OUT, buf, lDataLen, timeout);
+ if (ret < 0)
+ return ret;
+ else
+ bufsize=1600;
+ int ii=0;
+ {
+ ret=usb_bulk_read(hDev, XXUSB_ENDPOINT_IN, buf, bufsize, timeout);
+ if (ret>0)
+ for (i=0; i < ret; i=i+2)
+ intbuf[ii++]=(UCHAR)(buf[i]) + (UCHAR)(buf[i+1])*256;
+ return ret;
+
+ }
+}
+
+/*
+******** xxusb_register_read ************************
+
+ Reads the current contents of an internal xxusb register
+
+ Parameters:
+ hdev: USB device handle returned from an open function
+ RegAddr: The internal address of the register from which to read
+ RegData: Pointer to a long to hold the data.
+
+ Returns:
+ When Successful, the number of bytes read from xxusb.
+ Upon failure, a negative number
+*/
+short xxusb_register_read(usb_dev_handle *hDev, short RegAddr, long *RegData)
+{
+//long RegD;
+int timeout;
+ char buf[4]={1,0,0,0};
+ int ret;
+ int lDataLen;
+
+ buf[2]=(char)(RegAddr & 15);
+ timeout=10;
+ lDataLen=4;
+ ret=xxusb_bulk_write(hDev, buf, lDataLen, timeout);
+ if (ret < 0)
+ return (short)ret;
+ else
+ {
+ lDataLen=8;
+ timeout=100;
+ ret=xxusb_bulk_read(hDev, buf, lDataLen, timeout);
+ if (ret<0)
+ return (short)ret;
+ else
+ {
+ *RegData=(UCHAR)(buf[0])+256*(UCHAR)(buf[1]);
+ if (ret==4)
+ *RegData=*RegData+0x10000*(UCHAR)(buf[2]);
+ return (short)ret;
+ }
+ }
+}
+
+/*
+******** xxusb_reset_toggle ************************
+
+ Toggles the reset state of the FPGA while the xxusb in programming mode
+
+ Parameters
+ hdev: US B device handle returned from an open function
+
+ Returns:
+ Upon failure, a negative number
+*/
+short xxusb_reset_toggle(usb_dev_handle *hDev)
+{
+ short ret;
+ char buf[2] = {(char)255,(char)255};
+ int lDataLen=2;
+ int timeout=1000;
+ ret = usb_bulk_write(hDev, XXUSB_ENDPOINT_OUT, buf,lDataLen, timeout);
+ return (short)ret;
+}
+
+/*
+******** xxusb_devices_find ************************
+
+ Determines the number and parameters of all xxusb devices attched to
+ the computer.
+
+ Parameters:
+ xxdev: pointer to an array on which the device parameters are stored
+
+ Returns:
+ Upon success, returns the number of devices found
+ Upon Failure returns a negative number
+*/
+short xxusb_devices_find(xxusb_device_type *xxdev)
+{
+ short DevFound = 0;
+ usb_dev_handle *udev;
+ struct usb_bus *bus;
+ struct usb_device *dev;
+struct usb_bus *usb_busses;
+ char string[256];
+ short ret;
+ usb_init();
+ usb_find_busses();
+ usb_busses=usb_get_busses();
+ usb_find_devices();
+ for (bus=usb_busses; bus; bus = bus->next)
+ {
+ for (dev = bus->devices; dev; dev= dev->next)
+ {
+ if (dev->descriptor.idVendor==XXUSB_WIENER_VENDOR_ID)
+ {
+ udev = usb_open(dev);
+ if (udev)
+ {
+ ret = usb_get_string_simple(udev, dev->descriptor.iSerialNumber, string, sizeof(string));
+ if (ret >0 )
+ {
+ xxdev[DevFound].usbdev=dev;
+ strcpy(xxdev[DevFound].SerialString, string);
+ DevFound++;
+ }
+ usb_close(udev);
+ }
+ else return -1;
+ }
+ }
+ }
+ return DevFound;
+}
+
+/*
+******** xxusb_device_close ************************
+
+ Closes an xxusb device
+
+ Parameters:
+ hdev: USB device handle returned from an open function
+
+ Returns: 1
+*/
+short xxusb_device_close(usb_dev_handle *hDev)
+{
+ short ret;
+ ret=usb_release_interface(hDev,0);
+ usb_close(hDev);
+ return 1;
+}
+
+/*
+******** xxusb_device_open ************************
+
+ Opens an xxusb device found by xxusb_device_find
+
+ Parameters:
+ dev: a usb device
+
+ Returns:
+ A USB device handle
+*/
+usb_dev_handle* xxusb_device_open(struct usb_device *dev)
+{
+ short ret;
+ long val;
+ int count =0;
+ usb_dev_handle *udev;
+ udev = usb_open(dev);
+ ret = usb_set_configuration(udev,1);
+ ret = usb_claim_interface(udev,0);
+// RESET USB (added 10/16/06 Andreas Ruben)
+ ret=xxusb_register_write(udev, 10, 0x04);
+// Loop to find known state (added 12/28/07 TH / AR)
+ ret =-1;
+ while ((ret <0) && (count <10))
+ {
+ xxusb_register_read(udev, 0, &val);
+ count++;
+ }
+
+ return udev;
+}
+
+/*
+******** xxusb_flash_program ************************
+
+ --Untested and therefore uncommented--
+*/
+short xxusb_flash_program(usb_dev_handle *hDev, char *config, short nsect)
+{
+ int i=0;
+ int k=0;
+ short ret=0;
+ time_t t1,t2;
+
+ char *pconfig;
+ char *pbuf;
+ pconfig=config;
+ char buf[518] ={(char)0xAA,(char)0xAA,(char)0x55,(char)0x55,(char)0xA0,(char)0xA0};
+ while (*pconfig++ != -1);
+ for (i=0; i<nsect; i++)
+ {
+ pbuf=buf+6;
+ for (k=0; k<256; k++)
+ {
+ *(pbuf++)=*(pconfig);
+ *(pbuf++)=*(pconfig++);
+ }
+ ret = usb_bulk_write(hDev, XXUSB_ENDPOINT_OUT, buf, 518, 2000);
+ if (ret<0)
+ return ret;
+ t1=clock()+(time_t)(0.03*CLOCKS_PER_SEC);
+ while (t1>clock());
+ t2=clock();
+ }
+ return ret;
+}
+
+/*
+******** xxusb_flashblock_program ************************
+
+ --Untested and therefore uncommented--
+*/
+short xxusb_flashblock_program(usb_dev_handle *hDev, UCHAR *config)
+{
+ int k=0;
+ short ret=0;
+
+ UCHAR *pconfig;
+ char *pbuf;
+ pconfig=config;
+ char buf[518] ={(char)0xAA,(char)0xAA,(char)0x55,(char)0x55,(char)0xA0,(char)0xA0};
+ pbuf=buf+6;
+ for (k=0; k<256; k++)
+ {
+ *(pbuf++)=(UCHAR)(*(pconfig));
+ *(pbuf++)=(UCHAR)(*(pconfig++));
+ }
+ ret = usb_bulk_write(hDev, XXUSB_ENDPOINT_OUT, buf, 518, 2000);
+ return ret;
+}
+
+/*
+******** xxusb_serial_open ************************
+
+ Opens a xxusb device whose serial number is given
+
+ Parameters:
+ SerialString: a char string that gives the serial number of
+ the device you wish to open. It takes the form:
+ VM0009 - for a vm_usb with serial number 9 or
+ CC0009 - for a cc_usb with serial number 9
+
+ Returns:
+ A USB device handle
+*/
+usb_dev_handle* xxusb_serial_open(char *SerialString)
+{
+ short DevFound = 0;
+ usb_dev_handle *udev = NULL;
+ struct usb_bus *bus;
+ struct usb_device *dev;
+struct usb_bus *usb_busses;
+char string[7];
+ short ret;
+// usb_set_debug(4);
+ usb_init();
+ usb_find_busses();
+ usb_busses=usb_get_busses();
+ usb_find_devices();
+ for (bus=usb_busses; bus; bus = bus->next)
+ {
+ for (dev = bus->devices; dev; dev= dev->next)
+ {
+ if (dev->descriptor.idVendor==XXUSB_WIENER_VENDOR_ID)
+ {
+ udev = xxusb_device_open(dev);
+ if (udev)
+ {
+ ret = usb_get_string_simple(udev, dev->descriptor.iSerialNumber, string, sizeof(string));
+ if (ret >0 )
+ {
+ if (strcmp(string,SerialString)==0)
+ return udev;
+ }
+ usb_close(udev);
+ }
+ }
+ }
+ }
+ udev = NULL;
+ return udev;
+}
+
+
+//******************************************************//
+//****************** EZ_VME Functions ******************//
+//******************************************************//
+// The following are functions used to perform simple
+// VME Functions with the VM_USB
+
+/*
+******** VME_write_32 ************************
+
+ Writes a 32 bit data word to the VME bus
+
+ Parameters:
+ hdev: USB devcie handle returned from an open function
+ Address_Modifier: VME address modifier for the VME call
+ VME_Address: Address to write the data to
+ Data: 32 bit data word to be written to VME_Address
+
+ Returns:
+ Number of bytes read from xxusb when successful
+ Upon failure, a negative number
+*/
+short VME_write_32(usb_dev_handle *hdev, short Address_Modifier, long VME_Address, long Data)
+{
+ long intbuf[1000];
+ short ret;
+ intbuf[0]=7;
+ intbuf[1]=0;
+ intbuf[2]=Address_Modifier;
+ intbuf[3]=0;
+ intbuf[4]=(VME_Address & 0xffff);
+ intbuf[5]=((VME_Address >>16) & 0xffff);
+ intbuf[6]=(Data & 0xffff);
+ intbuf[7]=((Data >> 16) & 0xffff);
+ ret = xxusb_stack_execute(hdev, intbuf);
+ return ret;
+}
+
+/*
+******** VME_read_32 ************************
+
+
+ Reads a 32 bit data word from a VME address
+
+ Parameters:
+ hdev: USB devcie handle returned from an open function
+ Address_Modifier: VME address modifier for the VME call
+ VME_Address: Address to read the data from
+ Data: 32 bit data word read from VME_Address
+
+ Returns:
+ Number of bytes read from xxusb when successful
+ Upon failure, a negative number
+*/
+short VME_read_32(usb_dev_handle *hdev, short Address_Modifier, long VME_Address, long *Data)
+{
+ long intbuf[1000];
+ short ret;
+ intbuf[0]=5;
+ intbuf[1]=0;
+ intbuf[2]=Address_Modifier +0x100;
+ intbuf[3]=0;
+ intbuf[4]=(VME_Address & 0xffff);
+ intbuf[5]=((VME_Address >>16) & 0xffff);
+ ret = xxusb_stack_execute(hdev, intbuf);
+ *Data=intbuf[0] + (intbuf[1] * 0x10000);
+ return ret;
+}
+
+/*
+******** VME_write_16 ************************
+
+ Writes a 16 bit data word to the VME bus
+
+ Parameters:
+ hdev: USB devcie handle returned from an open function
+ Address_Modifier: VME address modifier for the VME call
+ VME_Address: Address to write the data to
+ Data: word to be written to VME_Address
+
+ Returns:
+ Number of bytes read from xxusb when successful
+ Upon failure, a negative number
+*/
+short VME_write_16(usb_dev_handle *hdev, short Address_Modifier, long VME_Address, long Data)
+{
+ long intbuf[1000];
+ short ret;
+ intbuf[0]=7;
+ intbuf[1]=0;
+ intbuf[2]=Address_Modifier;
+ intbuf[3]=0;
+ intbuf[4]=(VME_Address & 0xffff)+ 0x01;
+ intbuf[5]=((VME_Address >>16) & 0xffff);
+ intbuf[6]=(Data & 0xffff);
+ intbuf[7]=0;
+ ret = xxusb_stack_execute(hdev, intbuf);
+ return ret;
+}
+
+/*
+******** VME_read_16 ************************
+
+ Reads a 16 bit data word from a VME address
+
+ Parameters:
+ hdev: USB devcie handle returned from an open function
+ Address_Modifier: VME address modifier for the VME call
+ VME_Address: Address to read the data from
+ Data: word read from VME_Address
+
+ Returns:
+ Number of bytes read from xxusb when successful
+ Upon failure, a negative number
+*/
+short VME_read_16(usb_dev_handle *hdev,short Address_Modifier, long VME_Address, long *Data)
+{
+ long intbuf[1000];
+ short ret;
+ intbuf[0]=5;
+ intbuf[1]=0;
+ intbuf[2]=Address_Modifier +0x100;
+ intbuf[3]=0;
+ intbuf[4]=(VME_Address & 0xffff)+ 0x01;
+ intbuf[5]=((VME_Address >>16) & 0xffff);
+ ret = xxusb_stack_execute(hdev, intbuf);
+ *Data=intbuf[0];
+ return ret;
+}
+
+/*
+******** VME_BLT_read_32 ************************
+
+ Performs block transfer of 32 bit words from a VME address
+
+ Parameters:
+ hdev: USB devcie handle returned from an open function
+ Address_Modifier: VME address modifier for the VME call
+ count: number of data words to read
+ VME_Address: Address to read the data from
+ Data: pointer to an array to hold the data words
+
+ Returns:
+ Number of bytes read from xxusb when successful
+ Upon failure, a negative number
+*/
+short VME_BLT_read_32(usb_dev_handle *hdev, short Adress_Modifier, int count, long VME_Address, long Data[])
+{
+ long intbuf[1000];
+ short ret;
+ int i=0;
+ if (count > 255) return -1;
+ intbuf[0]=5;
+ intbuf[1]=0;
+ intbuf[2]=Adress_Modifier +0x100;
+ intbuf[3]=(count << 8);
+ intbuf[4]=(VME_Address & 0xffff);
+ intbuf[5]=((VME_Address >>16) & 0xffff);
+ ret = xxusb_stack_execute(hdev, intbuf);
+ int j=0;
+ for (i=0;i<(2*count);i=i+2)
+ {
+ Data[j]=intbuf[i] + (intbuf[i+1] * 0x10000);
+ j++;
+ }
+ return ret;
+}
+
+//******************************************************//
+//****************** VM_USB Registers ******************//
+//******************************************************//
+// The following are functions used to set the registers
+// in the VM_USB
+
+/*
+******** VME_register_write ************************
+
+ Writes to the vmusb registers that are accessible through
+ VME style calls
+
+ Parameters:
+ hdev: USB devcie handle returned from an open function
+ VME_Address: The VME Address of the internal register
+ Data: Data to be written to VME_Address
+
+ Returns:
+ Number of bytes read from xxusb when successful
+ Upon failure, a negative number
+*/
+short VME_register_write(usb_dev_handle *hdev, long VME_Address, long Data)
+{
+ long intbuf[1000];
+ short ret;
+
+ intbuf[0]=7;
+ intbuf[1]=0;
+ intbuf[2]=0x1000;
+ intbuf[3]=0;
+ intbuf[4]=(VME_Address & 0xffff);
+ intbuf[5]=((VME_Address >>16) & 0xffff);
+ intbuf[6]=(Data & 0xffff);
+ intbuf[7]=((Data >> 16) & 0xffff);
+ ret = xxusb_stack_execute(hdev, intbuf);
+ return ret;
+}
+
+/*
+******** VME_register_read ************************
+
+ Reads from the vmusb registers that are accessible trough VME style calls
+
+ Parameters:
+ hdev: USB devcie handle returned from an open function
+ VME_Address: The VME Address of the internal register
+ Data: Data read from VME_Address
+
+ Returns:
+ Number of bytes read from xxusb when successful
+ Upon failure, a negative number
+*/
+short VME_register_read(usb_dev_handle *hdev, long VME_Address, long *Data)
+{
+ long intbuf[1000];
+ short ret;
+
+ intbuf[0]=5;
+ intbuf[1]=0;
+ intbuf[2]=0x1100;
+ intbuf[3]=0;
+ intbuf[4]=(VME_Address & 0xffff);
+ intbuf[5]=((VME_Address >>16) & 0xffff);
+ ret = xxusb_stack_execute(hdev, intbuf);
+ *Data=intbuf[0] + (intbuf[1] * 0x10000);
+ return ret;
+}
+
+/*
+******** VME_LED_settings ************************
+
+ Sets the vmusb LED's
+
+ Parameters:
+ hdev: USB devcie handle returned from an open function
+ LED: The number which corresponds to an LED values are:
+ 0 - for Top YELLOW LED
+ 1 - for RED LED
+ 2 - for GREEN LED
+ 3 - for Bottom YELLOW LED
+ code: The LED aource selector code, valid values for each LED
+ are listed in the manual
+ invert: to invert the LED lighting
+ latch: sets LED latch bit
+
+ Returns:
+ Number of bytes read from xxusb when successful
+ Upon failure, a negative number
+*/
+short VME_LED_settings(usb_dev_handle *hdev, int LED, int code, int invert, int latch)
+{
+ short ret;
+// long internal;
+ long Data;
+
+ if( (LED <0) ||(LED > 3) || (code < 0) || (code > 7)) return -1;
+
+ VME_register_read(hdev,0xc,&Data);
+ if(LED == 0)
+ {
+ Data = Data & 0xFFFFFF00;
+ Data = Data | code;
+ if (invert == 1 && latch == 1) Data = Data | 0x18;
+ if (invert == 1 && latch == 0) Data = Data | 0x08;
+ if (invert == 0 && latch == 1) Data = Data | 0x10;
+ }
+ if(LED == 1)
+ {
+ Data = Data & 0xFFFF00FF;
+ Data = Data | (code * 0x0100);
+ if (invert == 1 && latch == 1) Data = Data | 0x1800;
+ if (invert == 1 && latch == 0) Data = Data | 0x0800;
+ if (invert == 0 && latch == 1) Data = Data | 0x1000;
+ }
+ if(LED == 2)
+ {
+ Data = Data & 0xFF00FFFF;
+ Data = Data | (code * 0x10000);
+ if (invert == 1 && latch == 1) Data = Data | 0x180000;
+ if (invert == 1 && latch == 0) Data = Data | 0x080000;
+ if (invert == 0 && latch == 1) Data = Data | 0x100000;
+ }
+ if(LED == 3)
+ {
+ Data = Data & 0x00FFFFFF;
+ Data = Data | (code * 0x10000);
+ if (invert == 1 && latch == 1) Data = Data | 0x18000000;
+ if (invert == 1 && latch == 0) Data = Data | 0x08000000;
+ if (invert == 0 && latch == 1) Data = Data | 0x10000000;
+ }
+ ret = VME_register_write(hdev, 0xc, Data);
+ return ret;
+}
+
+/*
+******** VME_DGG ************************
+
+ Sets the parameters for Gate & Delay channel A of vmusb
+
+ Parameters:
+ hdev: USB devcie handle returned from an open function
+ channel: Which DGG channel to use Valid Values are:
+ 0 - For DGG A
+ 1 - For DGG B
+ trigger: Determines what triggers the start of the DGG Valid values are:
+ 0 - Channel disabled
+ 1 - NIM input 1
+ 2 - NIM input 2
+ 3 - Event Trigger
+ 4 - End of Event
+ 5 - USB Trigger
+ 6 - Pulser
+ output: Determines which NIM output to use for the channel, Vaild values are:
+ 0 - for NIM O1
+ 1 - for NIM O2
+ delay: 32 bit word consisting of
+ lower 16 bits: Delay_fine in steps of 12.5ns between trigger and start of gate
+ upper 16 bits: Delay_coarse in steps of 81.7us between trigger and start of gate
+ gate: the time the gate should stay open in steps of 12.5ns
+ invert: is 1 if you wish to invert the DGG channel output
+ latch: is 1 if you wish to run the DGG channel latched
+
+ Returns:
+ Returns 1 when successful
+ Upon failure, a negative number
+*/
+short VME_DGG(usb_dev_handle *hdev, unsigned short channel, unsigned short trigger, unsigned short output,
+ long delay, unsigned short gate, unsigned short invert, unsigned short latch)
+{
+ long Data, DGData, Delay_ext;
+ long internal;
+ short ret;
+
+
+ ret = VME_register_read(hdev, 0x10, &Data);
+ // check and correct values
+ if(ret<=0) return -1;
+
+ if(channel >1) channel =1;
+ if(invert >1) invert =1;
+ if(latch >1) latch =1;
+ if(output >1) output =1;
+ if(trigger >6) trigger =0;
+
+ // define Delay and Gate data
+ DGData = gate * 0x10000;
+ DGData += (unsigned short) delay;
+
+ // Set channel, output, invert, latch
+ if (output == 0)
+ {
+ Data = Data & 0xFFFFFF00;
+ Data += 0x04 + channel +0x08*invert + 0x10*latch;
+ }
+ if (output == 1)
+ {
+ Data = Data & 0xFFFF00FF;
+ Data += (0x04 + channel +0x08*invert + 0x10*latch)*0x100;
+ }
+
+ // Set trigger, delay, gate
+
+ if(channel ==0) // CHANNEL DGG_A
+ {
+ internal = (trigger * 0x1000000) ;
+ Data= Data & 0xF0FFFFFF;
+ Data += internal;
+ ret = VME_register_write(hdev,0x10,Data);
+ if(ret<=0) return -1;
+ ret=VME_register_write(hdev,0x14,DGData);
+ if(ret<=0) return -1;
+ // Set coarse delay in DGG_Extended register
+ ret = VME_register_read(hdev,0x38,&Data);
+ Delay_ext= (Data & 0xffff0000);
+ Delay_ext+= ((delay/0x10000) & 0xffff);
+ ret = VME_register_write(hdev,0x38,Delay_ext);
+ }
+ if( channel ==1) // CHANNEL DGG_B
+ {
+ internal = (trigger * 0x10000000) ;
+ Data= Data & 0x0FFFFFFF;
+ Data += internal;
+ ret = VME_register_write(hdev,0x10,Data);
+ if(ret<=0) return -1;
+ ret=VME_register_write(hdev,0x18,DGData);
+ if(ret<=0) return -1;
+ // Set coarse delay in DGG_Extended register
+ ret = VME_register_read(hdev,0x38,&Data);
+ Delay_ext= (Data & 0x0000ffff);
+ Delay_ext+= (delay & 0xffff0000);
+ ret = VME_register_write(hdev,0x38,Delay_ext);
+ }
+ return 1;
+
+}
+
+/*
+******** VME_Output_settings ************************
+
+ Sets the vmusb NIM output register
+
+ Parameters:
+ hdev: USB devcie handle returned from an open function
+ Channel: The number which corresponds to an output:
+ 1 - for Output 1
+ 2 - for Output 2
+ code: The Output selector code, valid values
+ are listed in the manual
+ invert: to invert the output
+ latch: sets latch bit
+
+ Returns:
+ Number of bytes read from xxusb when successful
+ Upon failure, a negative number
+*/
+short VME_Output_settings(usb_dev_handle *hdev, int Channel, int code, int invert, int latch)
+{
+
+ short ret;
+// long internal;
+ long Data;
+
+ if( (Channel <1) ||(Channel > 2) || (code < 0) || (code > 7)) return -1;
+ VME_register_read(hdev,0x10,&Data);
+ if(Channel == 1)
+ {
+ Data = Data & 0xFFFF00;
+ Data = Data | code;
+ if (invert == 1 && latch == 1) Data = Data | 0x18;
+ if (invert == 1 && latch == 0) Data = Data | 0x08;
+ if (invert == 0 && latch == 1) Data = Data | 0x10;
+ }
+ if(Channel == 2)
+ {
+ Data = Data & 0xFF00FF;
+ Data = Data | (code * 0x0100);
+ if (invert == 1 && latch == 1) Data = Data | 0x1800;
+ if (invert == 1 && latch == 0) Data = Data | 0x0800;
+ if (invert == 0 && latch == 1) Data = Data | 0x1000;
+ }
+ ret = VME_register_write(hdev, 0x10, Data);
+ return ret;
+}
+
+
+//******************************************************//
+//****************** CC_USB Registers ******************//
+//******************************************************//
+// The following are functions used to set the registers
+// in the CAMAC_USB
+
+/*
+******** CAMAC_register_write *****************
+
+ Performs a CAMAC write to CC_USB register
+
+ Parameters:
+ hdev: USB device handle returned from an open function
+ A: CAMAC Subaddress
+ F: CAMAC Function
+ Data: data to be written
+
+ Returns:
+ Number of bytes written to xxusb when successful
+ Upon failure, a negative number
+*/
+short CAMAC_register_write(usb_dev_handle *hdev, int A, long Data)
+{
+ int F = 16;
+ int N = 25;
+ long intbuf[4];
+ int ret;
+
+ intbuf[0]=1;
+ intbuf[1]=(long)(F+A*32+N*512 + 0x4000);
+ intbuf[0]=3;
+ intbuf[2]=(Data & 0xffff);
+ intbuf[3]=((Data >>16) & 0xffff);
+ ret = xxusb_stack_execute(hdev, intbuf);
+
+ return ret;
+}
+
+/*
+******** CAMAC_register_read ************************
+
+ Performs a CAMAC read from CC_USB register
+
+ Parameters:
+ hdev: USB device handle returned from an open function
+ N: CAMAC Station Number
+ A: CAMAC Subaddress
+ F: CAMAC Function
+ Q: The Q response from the CAMAC dataway
+ X: The comment accepted response from CAMAC dataway
+
+ Returns:
+ Number of bytes read from xxusb when successful
+ Upon failure, a negative number
+*/
+short CAMAC_register_read(usb_dev_handle *hdev, int A, long *Data)
+{
+ int F = 0;
+ int N = 25;
+ long intbuf[4];
+ int ret;
+
+ intbuf[0]=1;
+ intbuf[1]=(long)(F+A*32+N*512 + 0x4000);
+ ret = xxusb_stack_execute(hdev, intbuf);
+ *Data=intbuf[0] + (intbuf[1] * 0x10000);
+
+ return ret;
+}
+
+/*
+******** CAMAC_DGG ************************
+
+ Sets the parameters for Gate & Delay channel A of CC_USB
+
+ Parameters:
+ hdev: USB devcie handle returned from an open function
+ channel: Which DGG channel to use Valid Values are:
+ 0 - For DGG A
+ 1 - For DGG B
+ trigger: Determines what triggers the start of the DGG Valid values are:
+ 0 - Channel disabled
+ 1 - NIM input 1
+ 2 - NIM input 2
+ 3 - NIM input 2
+ 4 - Event Trigger
+ 5 - End of Event
+ 6 - USB Trigger
+ 7 - Pulser
+ output: Determines which NIM output to use for the channel, Vaild values are:
+ 1 - for NIM O1
+ 2 - for NIM O2
+ 3 - for NIM O3
+ delay: Delay in steps of 12.5ns between trigger and start of gate
+ gate: the time the gate should stay open in steps of 12.5ns
+ invert: is 1 if you wish to invert the DGG channel output
+ latch: is 1 if you wish to run the DGG channel latched
+
+ Returns:
+ Returns 1 when successful
+ Upon failure, a negative number
+*/
+short CAMAC_DGG(usb_dev_handle *hdev, short channel, short trigger, short output,
+ int delay, int gate, short invert, short latch)
+
+
+
+{
+// short channel_ID;
+ long Data;
+ long internal;
+ short ret;
+ long Delay_ext;
+
+ ret = CAMAC_register_read(hdev,5,&Data);
+ //Set trigger
+ if((output < 1 ) || (output >3) || (channel < 0 ) || (channel > 1))
+ return -1;
+ if(output ==1)
+ {
+ if(channel ==0)
+ {
+ internal = 0x03;
+ } else {
+ internal = 0x04;
+ }
+ }
+ if(output ==2)
+ {
+ if(channel ==0)
+ {
+ internal = 0x04;
+ } else {
+ internal = 0x05;
+ }
+ }
+ if(output ==3)
+ {
+ if(channel ==0)
+ {
+ internal = 0x05;
+ } else {
+ internal = 0x06;
+ }
+ }
+
+
+// Set invert bit
+ if(invert ==1)
+ internal = internal | 0x10;
+ else
+ internal = internal & 0x0F;
+ // Set Latch Bit
+ if(latch==1)
+ internal = internal | 0x20;
+ else
+ internal = internal & 0x1F;
+// Add new data to old
+ if(output == 1)
+ {
+ Data = Data & 0xFFFF00;
+ Data = Data | internal;
+ }
+ if(output == 2)
+ {
+ Data = Data & 0xFF00FF;
+ Data = Data |(internal * 0x100);
+ }
+ if(output == 3)
+ {
+ Data = Data & 0x00FFFF;
+ Data = Data | (internal * 0x10000) ;
+ }
+ CAMAC_register_write(hdev, 5, Data);
+ ret = CAMAC_register_read(hdev,6,&Data);
+//Set Trigger
+ if(trigger <0 || trigger > 7)
+ return -1;
+ if(channel ==0)
+ {
+ Data = Data & 0xFF00FFFF;
+ internal = trigger * 0x10000;
+ Data = Data | internal;
+ } else {
+ Data = Data & 0x00FFFFFF;
+ internal = trigger * 0x1000000;
+ Data = Data | internal;
+ }
+ ret = CAMAC_register_write(hdev, 6, Data);
+ if(channel == 0)
+ {
+// Write Delay and Gate info
+ ret = CAMAC_register_read(hdev,13,&Data);
+ Delay_ext= (Data & 0xffff0000);
+ Delay_ext+= ((delay/0x10000) & 0xffff);
+ internal = gate * 0x10000;
+ Data = internal + (delay & 0xffff);
+ ret=CAMAC_register_write(hdev,7,Data);
+// Set coarse delay in DGG_Extended register
+ ret=CAMAC_register_write(hdev,13,Delay_ext);
+ }
+ else
+ {
+ ret=CAMAC_register_write(hdev,8,Data);
+ ret = CAMAC_register_read(hdev,13,&Data);
+ Delay_ext= (Data & 0x0000ffff);
+ Delay_ext+= (delay & 0xffff0000);
+ internal = gate * 0x10000;
+ Data = internal + (delay & 0xffff);
+// Set coarse delay in DGG_Extended register
+ ret=CAMAC_register_write(hdev,13,Delay_ext);
+ }
+ return 1;
+}
+
+/*
+******** CAMAC_LED_settings ************************
+
+ Writes a data word to the vmusb LED register
+
+ Parameters:
+ hdev: USB devcie handle returned from an open function
+ LED: The number which corresponds to an LED values are:
+ 1 - for RED LED
+ 2 - for GREEN LED
+ 3 - for Yellow LED
+ code: The LED aource selector code, valid values for each LED
+ are listed in the manual
+ invert: to invert the LED lighting
+ latch: sets LED latch bit
+
+ Returns:
+ Number of bytes read from xxusb when successful
+ Upon failure, a negative number
+*/
+short CAMAC_LED_settings(usb_dev_handle *hdev, int LED, int code, int invert, int latch)
+{
+
+ short ret;
+// long internal;
+ long Data;
+
+ if( (LED <1) ||(LED > 3) || (code < 0) || (code > 7))
+ return -1;
+
+ CAMAC_register_read(hdev,4,&Data);
+
+ if(LED == 1)
+ {
+ Data = Data & 0xFFFF00;
+ Data = Data | code;
+ if (invert == 1 && latch == 1)
+ Data = Data | 0x30;
+ if (invert == 1 && latch == 0)
+ Data = Data | 0x10;
+ if (invert == 0 && latch == 1)
+ Data = Data | 0x20;
+ }
+ if(LED == 2)
+ {
+ Data = Data & 0xFF00FF;
+ Data = Data | (code * 0x0100);
+ if (invert == 1 && latch == 1)
+ Data = Data | 0x3000;
+ if (invert == 1 && latch == 0)
+ Data = Data | 0x1000;
+ if (invert == 0 && latch == 1)
+ Data = Data | 0x2000;
+ }
+ if(LED == 3)
+ {
+ Data = Data & 0x00FFFF;
+ Data = Data | (code * 0x10000);
+ if (invert == 1 && latch == 1)
+ Data = Data | 0x300000;
+ if (invert == 1 && latch == 0)
+ Data = Data | 0x100000;
+ if (invert == 0 && latch == 1)
+ Data = Data | 0x200000;
+ }
+ ret = CAMAC_register_write(hdev, 4, Data);
+ return ret;
+}
+
+/*
+******** CAMAC_Output_settings ************************
+
+ Writes a data word to the vmusb LED register
+
+ Parameters:
+ hdev: USB devcie handle returned from an open function
+ Channel: The number which corresponds to an output:
+ 1 - for Output 1
+ 2 - for Output 2
+ 3 - for Output 3
+ code: The Output selector code, valid values
+ are listed in the manual
+ invert: to invert the output
+ latch: sets latch bit
+
+ Returns:
+ Number of bytes read from xxusb when successful
+ Upon failure, a negative number
+*/
+short CAMAC_Output_settings(usb_dev_handle *hdev, int Channel, int code, int invert, int latch)
+{
+ short ret;
+// long internal;
+ long Data;
+
+ if( (Channel <1) ||(Channel > 3) || (code < 0) || (code > 7))
+ return -1;
+
+ CAMAC_register_read(hdev,5,&Data);
+
+ if(Channel == 1)
+ {
+ Data = Data & 0xFFFF00;
+ Data = Data | code;
+ if (invert == 1 && latch == 1)
+ Data = Data | 0x30;
+ if (invert == 1 && latch == 0)
+ Data = Data | 0x10;
+ if (invert == 0 && latch == 1)
+ Data = Data | 0x20;
+ }
+ if(Channel == 2)
+ {
+ Data = Data & 0xFF00FF;
+ Data = Data | (code * 0x0100);
+ if (invert == 1 && latch == 1)
+ Data = Data | 0x3000;
+ if (invert == 1 && latch == 0)
+ Data = Data | 0x1000;
+ if (invert == 0 && latch == 1)
+ Data = Data | 0x2000;
+ }
+ if(Channel == 3)
+ {
+ Data = Data & 0x00FFFF;
+ Data = Data | (code * 0x10000);
+ if (invert == 1 && latch == 1)
+ Data = Data | 0x300000;
+ if (invert == 1 && latch == 0)
+ Data = Data | 0x100000;
+ if (invert == 0 && latch == 1)
+ Data = Data | 0x200000;
+ }
+ ret = CAMAC_register_write(hdev, 5, Data);
+ return ret;
+}
+
+/*
+******** CAMAC_write_LAM_mask ************************
+
+ Writes the data word to the LAM mask register
+
+ Parameters:
+ hdev: USB devcie handle returned from an open function
+ Data: LAM mask to write
+
+ Returns:
+ Number of bytes read from xxusb when successful
+ Upon failure, a negative number
+*/
+short CAMAC_write_LAM_mask(usb_dev_handle *hdev, long Data)
+{
+ short ret;
+ ret = CAMAC_register_write(hdev, 9, Data);
+
+ return ret;
+}
+
+/*
+******** CAMAC_read_LAM_mask ************************
+
+ Writes the data word to the LAM mask register
+
+ Parameters:
+ hdev: USB devcie handle returned from an open function
+ Data: LAM mask to write
+
+ Returns:
+ Number of bytes read from xxusb when successful
+ Upon failure, a negative number
+*/
+short CAMAC_read_LAM_mask(usb_dev_handle *hdev, long *Data)
+{
+ long intbuf[4];
+ int ret;
+ int N = 25;
+ int F = 0;
+ int A = 9;
+
+ // CAMAC direct read function
+ intbuf[0]=1;
+ intbuf[1]=(long)(F+A*32+N*512 + 0x4000);
+ ret = xxusb_stack_execute(hdev, intbuf);
+ *Data=intbuf[0] + (intbuf[1] & 255) * 0x10000;
+ return ret;
+}
+
+
+//******************************************************//
+//**************** EZ_CAMAC Functions ******************//
+//******************************************************//
+// The following are functions used to perform simple
+// CAMAC Functions with the CC_USB
+
+
+/*
+******** CAMAC_write ************************
+
+ Performs a CAMAC write using NAF comments
+
+ Parameters:
+ hdev: USB device handle returned from an open function
+ N: CAMAC Station Number
+ A: CAMAC Subaddress
+ F: CAMAC Function (16...23)
+ Q: The Q response from the CAMAC dataway
+ X: The comment accepted response from CAMAC dataway
+
+ Returns:
+ Number of bytes written to xxusb when successful
+ Upon failure, a negative number
+*/
+short CAMAC_write(usb_dev_handle *hdev, int N, int A, int F, long Data, int *Q, int *X)
+{
+ long intbuf[4];
+ int ret;
+// CAMAC direct write function
+ intbuf[0]=1;
+ intbuf[1]=(long)(F+A*32+N*512 + 0x4000);
+ if ((F > 15) && (F < 24))
+ {
+ intbuf[0]=3;
+ intbuf[2]=(Data & 0xffff);
+ intbuf[3]=((Data >>16) & 255);
+ ret = xxusb_stack_execute(hdev, intbuf);
+ *Q = (intbuf[0] & 1);
+ *X = ((intbuf[0] >> 1) & 1);
+ }
+ return ret;
+}
+
+/*
+******** CAMAC_read ************************
+
+ Performs a CAMAC read using NAF comments
+
+ Parameters:
+ hdev: USB device handle returned from an open function
+ N: CAMAC Station Number
+ A: CAMAC Subaddress
+ F: CAMAC Function (F<16 or F>23)
+ Q: The Q response from the CAMAC dataway
+ X: The comment accepted response from CAMAC dataway
+
+ Returns:
+ Number of bytes read from xxusb when successful
+ Upon failure, a negative number
+*/
+short CAMAC_read(usb_dev_handle *hdev, int N, int A, int F, long *Data, int *Q, int *X)
+{
+ long intbuf[4];
+ int ret;
+ // CAMAC direct read function
+ intbuf[0]=1;
+ intbuf[1]=(long)(F+A*32+N*512 + 0x4000);
+ ret = xxusb_stack_execute(hdev, intbuf);
+ if ((F < 16) || (F >23))
+ {
+ *Data=intbuf[0] + (intbuf[1] & 255) * 0x10000; //24-bit word
+ *Q = ((intbuf[1] >> 8) & 1);
+ *X = ((intbuf[1] >> 9) & 1);
+ }
+ return ret;
+}
+
+/*
+******** CAMAC_Z ************************
+
+ Performs a CAMAC init
+
+ Parameters:
+ hdev: USB device handle returned from an open function
+
+ Returns:
+ Number of bytes written to xxusb when successful
+ Upon failure, a negative number
+*/
+short CAMAC_Z(usb_dev_handle *hdev)
+{
+ long intbuf[4];
+ int ret;
+// CAMAC Z = N(28) A(8) F(29)
+ intbuf[0]=1;
+ intbuf[1]=(long)(29+8*32+28*512 + 0x4000);
+ ret = xxusb_stack_execute(hdev, intbuf);
+ return ret;
+}
+
+/*
+******** CAMAC_C ************************
+
+ Performs a CAMAC clear
+
+ Parameters:
+ hdev: USB device handle returned from an open function
+
+ Returns:
+ Number of bytes written to xxusb when successful
+ Upon failure, a negative number
+*/
+short CAMAC_C(usb_dev_handle *hdev)
+{
+ long intbuf[4];
+ int ret;
+ intbuf[0]=1;
+ intbuf[1]=(long)(29+9*32+28*512 + 0x4000);
+ ret = xxusb_stack_execute(hdev, intbuf);
+ return ret;
+}
+
+/*
+******** CAMAC_I ************************
+
+ Set CAMAC inhibit
+
+ Parameters:
+ hdev: USB device handle returned from an open function
+
+ Returns:
+ Number of bytes written to xxusb when successful
+ Upon failure, a negative number
+*/
+short CAMAC_I(usb_dev_handle *hdev, int inhibit)
+{
+ long intbuf[4];
+ int ret;
+ intbuf[0]=1;
+ if (inhibit) intbuf[1]=(long)(24+9*32+29*512 + 0x4000);
+ else intbuf[1]=(long)(26+9*32+29*512 + 0x4000);
+ ret = xxusb_stack_execute(hdev, intbuf);
+ return ret;
+}
+
+
Index: apps/petdemo/libxxusb.dll
===================================================================
Cannot display: file marked as a binary type.
svn:mime-type = application/octet-stream
/apps/petdemo/libxxusb.dll
Property changes:
Added: svn:mime-type
## -0,0 +1 ##
+application/octet-stream
\ No newline at end of property
Index: apps/petdemo/libxxusb.h
===================================================================
--- apps/petdemo/libxxusb.h (nonexistent)
+++ apps/petdemo/libxxusb.h (revision 264)
@@ -0,0 +1,111 @@
+#include <usb.h>
+
+
+#define XXUSB_WIENER_VENDOR_ID 0x16DC /* Wiener, Plein & Baus */
+#define XXUSB_VMUSB_PRODUCT_ID 0x000B /* VM-USB */
+#define XXUSB_CCUSB_PRODUCT_ID 0x0001 /* CC-USB */
+#define XXUSB_ENDPOINT_OUT 2 /* Endpoint 2 Out*/
+#define XXUSB_ENDPOINT_IN 0x86 /* Endpoint 6 In */
+#define XXUSB_FIRMWARE_REGISTER 0
+#define XXUSB_GLOBAL_REGISTER 1
+#define XXUSB_ACTION_REGISTER 10
+#define XXUSB_DELAYS_REGISTER 2
+#define XXUSB_WATCHDOG_REGISTER 3
+#define XXUSB_SELLEDA_REGISTER 6
+#define XXUSB_SELNIM_REGISTER 7
+#define XXUSB_SELLEDB_REGISTER 4
+#define XXUSB_SERIAL_REGISTER 15
+#define XXUSB_LAMMASK_REGISTER 8
+#define XXUSB_LAM_REGISTER 12
+#define XXUSB_READOUT_STACK 2
+#define XXUSB_SCALER_STACK 3
+#define XXUSB_NAF_DIRECT 12
+
+struct XXUSB_STACK
+{
+long Data;
+short Hit;
+short APatt;
+short Num;
+short HitMask;
+};
+
+struct XXUSB_CC_COMMAND_TYPE
+{
+short Crate;
+short F;
+short A;
+short N;
+long Data;
+short NoS2;
+short LongD;
+short HitPatt;
+short QStop;
+short LAMMode;
+short UseHit;
+short Repeat;
+short AddrScan;
+short FastCam;
+short NumMod;
+short AddrPatt;
+long HitMask[4];
+long Num;
+};
+
+struct xxusb_device_typ
+{
+ struct usb_device *usbdev;
+ char SerialString[7];
+};
+
+typedef struct xxusb_device_typ xxusb_device_type;
+typedef unsigned char UCHAR;
+typedef struct usb_bus usb_busx;
+
+
+int xxusb_longstack_execute(usb_dev_handle *hDev, void *DataBuffer, int lDataLen, int timeout);
+int xxusb_bulk_read(usb_dev_handle *hDev, void *DataBuffer, int lDataLen, int timeout);
+int xxusb_bulk_write(usb_dev_handle *hDev, void *DataBuffer, int lDataLen, int timeout);
+int xxusb_usbfifo_read(usb_dev_handle *hDev, int *DataBuffer, int lDataLen, int timeout);
+
+short xxusb_register_read(usb_dev_handle *hDev, short RegAddr, long *RegData);
+short xxusb_stack_read(usb_dev_handle *hDev, short StackAddr, long *StackData);
+short xxusb_stack_write(usb_dev_handle *hDev, short StackAddr, long *StackData);
+short xxusb_stack_execute(usb_dev_handle *hDev, long *StackData);
+short xxusb_register_write(usb_dev_handle *hDev, short RegAddr, long RegData);
+short xxusb_reset_toggle(usb_dev_handle *hDev);
+
+short xxusb_devices_find(xxusb_device_type *xxusbDev);
+short xxusb_device_close(usb_dev_handle *hDev);
+usb_dev_handle* xxusb_device_open(struct usb_device *dev);
+short xxusb_flash_program(usb_dev_handle *hDev, char *config, short nsect);
+short xxusb_flashblock_program(usb_dev_handle *hDev, UCHAR *config);
+usb_dev_handle* xxusb_serial_open(char *SerialString);
+
+short VME_register_write(usb_dev_handle *hdev, long VME_Address, long Data);
+short VME_register_read(usb_dev_handle *hdev, long VME_Address, long *Data);
+short VME_LED_settings(usb_dev_handle *hdev, int LED, int code, int invert, int latch);
+short VME_DGG(usb_dev_handle *hdev, unsigned short channel, unsigned short trigger,unsigned short output, long delay, unsigned short gate, unsigned short invert, unsigned short latch);
+
+short VME_Output_settings(usb_dev_handle *hdev, int Channel, int code, int invert, int latch);
+
+short VME_read_16(usb_dev_handle *hdev,short Address_Modifier, long VME_Address, long *Data);
+short VME_read_32(usb_dev_handle *hdev, short Address_Modifier, long VME_Address, long *Data);
+short VME_BLT_read_32(usb_dev_handle *hdev, short Address_Modifier, int count, long VME_Address, long Data[]);
+short VME_write_16(usb_dev_handle *hdev, short Address_Modifier, long VME_Address, long Data);
+short VME_write_32(usb_dev_handle *hdev, short Address_Modifier, long VME_Address, long Data);
+
+short CAMAC_DGG(usb_dev_handle *hdev, short channel, short trigger, short output, int delay, int gate, short invert, short latch);
+short CAMAC_register_read(usb_dev_handle *hdev, int A, long *Data);
+short CAMAC_register_write(usb_dev_handle *hdev, int A, long Data);
+short CAMAC_LED_settings(usb_dev_handle *hdev, int LED, int code, int invert, int latch);
+short CAMAC_Output_settings(usb_dev_handle *hdev, int Channel, int code, int invert, int latch);
+short CAMAC_read_LAM_mask(usb_dev_handle *hdev, long *Data);
+short CAMAC_write_LAM_mask(usb_dev_handle *hdev, long Data);
+
+short CAMAC_write(usb_dev_handle *hdev, int N, int A, int F, long Data, int *Q, int *X);
+short CAMAC_read(usb_dev_handle *hdev, int N, int A, int F, long *Data, int *Q, int *X);
+short CAMAC_Z(usb_dev_handle *hdev);
+short CAMAC_C(usb_dev_handle *hdev);
+short CAMAC_I(usb_dev_handle *hdev, int inhibit);
+
Index: apps/petdemo/libxxusb.lib
===================================================================
Cannot display: file marked as a binary type.
svn:mime-type = application/octet-stream
/apps/petdemo/libxxusb.lib
Property changes:
Added: svn:mime-type
## -0,0 +1 ##
+application/octet-stream
\ No newline at end of property
Index: apps/petdemo/petdemo.c
===================================================================
--- apps/petdemo/petdemo.c (nonexistent)
+++ apps/petdemo/petdemo.c (revision 264)
@@ -0,0 +1,511 @@
+#include "H2D.h"
+#include "H1D.h"
+#include "inifile.h"
+#include <utility.h>
+#include <ansi_c.h>
+#include <cvirte.h>
+#include <userint.h>
+#include "petdemo_uir.h"
+
+#include "daq.h"
+#include "PETProjDataMgr.h"
+#include <math.h>
+
+#define uSMC_USB
+#ifdef uSMC_USB
+# include "uSMC.h"
+# define uSMC_SERIAL "0000000000005660"
+const char serials[3][16]= {uSMC_SERIAL};
+#endif /* uSMC_USB */
+
+static int node;
+static int CurPos = 0;
+
+int p1;
+#define MAX_THREADS 10
+static CmtThreadPoolHandle poolHandle = 0;
+int ctrl_c=0;
+int daq_on=0;
+int gadcchannel;
+static int tfID;
+static int controlID;
+static int plothandle[0xFF];
+char strbuf[0xFF];
+
+int gLog=0;
+
+
+
+int printf(const char *format, ...) {
+ va_list aptr;
+ int ret;
+ FILE *flog;
+
+ va_start(aptr, format);
+ ret = vsprintf(strbuf, format, aptr);
+ va_end(aptr);
+ SetCtrlVal(p1,P1_STDIO,strbuf);
+
+ if (gLog) {
+ flog = fopen ("stdio.log", "a");
+ fprintf (flog, "%s", strbuf);
+ fclose (flog);
+ }
+ return(ret);
+}
+
+void CVICALLBACK EndOfThread ( CmtThreadPoolHandle poolhandle,
+ CmtThreadFunctionID functionID, unsigned int event,
+ int value, void *callbackData ) {
+
+ daq_on=0;
+ //SetDimming(0);
+ printf("End of Thread \n");
+ return ;
+
+}
+
+
+
+
+int __stdcall WinMain (HINSTANCE hInstance, HINSTANCE hPrevInstance,
+ LPSTR lpszCmdLine, int nCmdShow) {
+ if (InitCVIRTE (hInstance, 0, 0) == 0)
+ return -1; /* out of memory */
+ if ((p1 = LoadPanel (0, "petdemo_uir.uir", P1)) < 0)
+ return -1;
+ SetStdioPort (CVI_STDIO_WINDOW);
+ SetSleepPolicy(VAL_SLEEP_MORE);
+ CmtNewThreadPool (MAX_THREADS, &poolHandle);
+#ifdef uSMC_USB
+ uSMC_Open();
+ node=uSMC_FindSerial(serials)+1;
+ uSMC_Init(node,1);
+ SetCtrlVal(p1, P1_STAGELED,1);
+
+ uSMC_GetPosition(node,&CurPos);
+ SetCtrlVal(p1,P1_CPOSITION,CurPos);
+
+ SetCtrlVal(p1, P1_STAGELED,0);
+
+#endif /* uSMC_USB */
+
+ DisplayPanel (p1);
+ RunUserInterface ();
+ DiscardPanel (p1);
+ CmtDiscardThreadPool (poolHandle);
+
+#ifdef uSMC_USB
+ uSMC_PowerOff(node);
+ uSMC_Close();
+#endif /* uSMC_USB */
+
+ return 0;
+}
+
+
+
+
+int addheader(FILE *fp, int type, int *data) {
+
+
+ time_t t;
+
+
+
+ switch(type) {
+ case RUNREC_ID: {
+
+
+
+ RUNREC runrec; // start header: appears once at the start of the file
+
+ runrec.id = RUNREC_ID;
+ runrec.length = sizeof(runrec);
+ runrec.fver = 0x10000;
+ time(&t);
+ runrec.time=t;
+ runrec.num_events = data[6];
+ runrec.num_channels = NUM_CHANNELS;
+ runrec.pedestal = PEDESTAL;
+ runrec.xy = SCAN_TYPE; // 0 -> single data scan :: 1 -> XY position scan
+ runrec.nx = data[0];
+ runrec.x0 = data[1];
+ runrec.dx = data[2];
+ runrec.ny = data[3];
+ runrec.y0 = data[4];
+ runrec.dy = data[5];
+
+ printf("Writing header to file\n");
+ printf("RECID = %u\n",runrec.id);
+ printf("Length = %u\n",runrec.length);
+ printf("File version = %u\n",runrec.fver);
+ printf("Time = %u\n",runrec.time);
+ printf("Number of events per step = %u\n",runrec.num_events);
+ printf("Number of channels measured = %u\n",runrec.num_channels);
+ printf("Pedestal = %u\n",runrec.pedestal);
+ printf("Scan type = %u :: 0 -> single data scan :: 1 -> XY position scan\n",runrec.xy);
+ printf("Number of steps in X = %d\n",runrec.nx);
+ printf("Start position X = %d\n",runrec.x0);
+ printf("Step size direction X = %d\n",runrec.dx);
+ printf("Number of steps in Y = %d\n",runrec.ny);
+ printf("Start position Y = %d\n",runrec.y0);
+ printf("Step size direction Y = %d\n",runrec.dy);
+
+ if (fp) fwrite(&runrec, runrec.length,1,fp);
+
+ }
+ break;
+
+
+ case ENDREC_ID: {
+
+
+ ENDREC endrec; // end header: appears once at the end of the file
+
+ endrec.id = ENDREC_ID;
+ endrec.length = sizeof(endrec);
+ time(&t);
+ endrec.time=t;
+
+ printf("Writing header to file\n");
+ printf("RECID = %u\n",endrec.id);
+ printf("Length = %u\n",endrec.length);
+ printf("Time = %u\n",endrec.time);
+
+ if (fp) fwrite(&endrec, endrec.length,1,fp);
+
+
+ break;
+ }
+ case POSREC_ID: {
+
+
+ POSREC posrec; // position header: appears at every change of position
+
+ posrec.id = POSREC_ID;
+ posrec.length = sizeof(posrec);
+ time(&t);
+ posrec.time = t;
+ posrec.num_iter_x = data[0];
+ posrec.mikro_pos_x = data[1];
+ posrec.set_pos_x = data[2];
+ posrec.num_iter_y = data[3];
+ posrec.mikro_pos_y = data[4];
+ posrec.set_pos_y = data[5];
+
+ printf("Writing header to file\n");
+ printf("RECID = %u\n",posrec.id);
+ printf("Length = %u\n",posrec.length);
+ printf("Time = %u\n",posrec.time);
+ printf("Iteration X = %d\n",posrec.num_iter_x);
+ printf("MIKRO Position X = %d\n",posrec.mikro_pos_x);
+ printf("Set position X = %d\n",posrec.set_pos_x);
+ printf("Iteration Y = %d\n",posrec.num_iter_y);
+ printf("MIKRO Position Y = %d\n",posrec.mikro_pos_y);
+ printf("Set position Y = %d\n",posrec.set_pos_y);
+
+ if (fp) fwrite(&posrec, posrec.length,1,fp);
+
+
+
+ break;
+ }
+ }
+
+ return 0;
+}
+
+
+int CVICALLBACK ExitCB (int panel, int control, int event,
+ void *callbackData, int eventData1, int eventData2) {
+ switch (event) {
+ case EVENT_COMMIT:
+ QuitUserInterface (0);
+ break;
+ }
+ return 0;
+}
+
+
+int Fi2Pos(float x) {
+ // x in degrees
+ return (int) (x/360*NSTEPS);
+}
+
+
+int mdaq_init(const char *fname){
+
+ IniText iniText;
+char pathName[MAX_PATHNAME_LEN];
+char dirName[MAX_PATHNAME_LEN];
+
+
+/* set up the pathName for the .ini file */
+GetProjectDir (dirName);
+MakePathname (dirName, fname, pathName);
+
+/* create object for holding the value/tag pairs */
+iniText = Ini_New (TRUE); /* TRUE for automatic sorting */
+
+/* read in the tag/value pairs */
+Ini_ReadFromFile (iniText, pathName);
+
+/* create the in–memory tag/value pairs */
+
+Ini_GetInt (iniText, "scintillator", "nofcrystalsx", &conf.nofcrystalsx);
+Ini_GetInt (iniText, "scintillator", "nofcrystalsy", &conf.nofcrystalsy);
+
+Ini_GetDouble (iniText, "scintillator", "crystalpitchx", &conf.crystalpitchx);
+Ini_GetDouble (iniText, "scintillator", "crystalpitchy", &conf.crystalpitchy);
+
+
+Ini_GetStringCopy (iniText, "sensor", "modules", &conf.modules);
+Ini_GetStringCopy (iniText, "sensor", "channels", &conf.channels);
+
+
+Ini_GetStringCopy (iniText, "calibration", "sumpedestals", &conf.sumpedestals);
+Ini_GetStringCopy (iniText, "calibration", "pedestals", &conf.pedestals);
+Ini_GetStringCopy (iniText, "calibration", "photopeak", &conf.photopeak);
+Ini_GetStringCopy (iniText, "calibration", "channelcalibration", &conf.channelcalibration);
+Ini_GetInt (iniText, "calibration", "adcthreshold", &conf.adcthreshold);
+
+/* dispose of the in–memory tag/value pairs */
+Ini_Dispose (iniText);
+return 0;
+}
+
+
+void ReadModuleMap(const char *fname){
+ int id;
+ float r,phi;
+ char line[400];
+ const int ndim=400;
+ FILE *fp=fopen(fname,"r");
+ if (!fp) return;
+ while (fgets(line,ndim,fp)!=NULL) {
+ sscanf(line,"%d%f%f",&id,&r,&phi);
+ if (id<16) {
+ conf.module[id].r=r;
+ conf.module[id].phi=phi*Pi()/180.;
+ if (debug) printf("%s %d %f %f\n",fname,id, r, phi);
+ }
+ }
+ fclose(fp);
+
+}
+
+void ReadChannelMap(const char *fname){
+ int id;
+ int ix, iy;
+ char line[400];
+ const int ndim=400;
+ FILE *fp=fopen(fname,"r");
+ if (!fp) return;
+ while (fgets(line,ndim,fp)!=NULL) {
+ sscanf(line,"%d%d%d",&id,&ix,&iy);
+ if (id<16) {
+ conf.channel[id].ix=ix;
+ conf.channel[id].iy=iy;
+ if (debug) printf("%s %d %d %d\n",fname,id, ix, iy);
+ }
+ }
+ fclose(fp);
+
+}
+
+ int readfile(const char *fname, float *x, int nmax, int defaultvalue){
+ int id;
+ float ix;
+ char line[400];
+ const int ndim=400;
+ FILE *fp=fopen(fname,"r");
+ for (int i=0;i<nmax;i++){
+ x[i]=defaultvalue;
+ }
+ if (!fp) return -1;
+ while (fgets(line,ndim,fp)!=NULL) {
+ sscanf(line,"%d%f",&id,&ix);
+
+ if (id<nmax) x[id]=ix;
+ if (debug) printf("%s %d %f\n",fname,id, ix);
+
+ }
+ fclose(fp);
+ return 0;
+ }
+
+
+
+
+
+
+
+
+
+
+
+
+int Geometry(const char *fnameconfig){
+
+ mdaq_init(fnameconfig);
+
+ ReadModuleMap(conf.modules);
+ ReadChannelMap(conf.channels);
+
+ printf( "Reading ...%s\n " ,conf.channelcalibration );
+ /*
+ m_calibration = new TFile(m_calibrationrootName);
+ for (int i=0; i<4;i++) {
+ char hn[256];
+ sprintf(hn,"pmt1%d_calib",i);
+ m_crystalid[i]= (TH1I *) m_calibration->Get(hn);
+ m_crystalid[i]->ls();
+ }
+ */
+
+ conf.peakscaling=3000;
+
+ readfile(conf.pedestals,conf.apedestals,4*16, 0);
+ readfile(conf.photopeak,conf.apeak,4*16, conf.peakscaling);
+ return 0;
+ };
+
+
+int CVICALLBACK mdaq(void *functionData) {
+
+ double dfi, fi0;
+ unsigned int nsteps;
+ char filename[0xFF];
+ char fname[0xFF];
+ int daqtime;
+ daq_on = 1;
+ ctrl_c=0;
+ GetCtrlVal(p1,P1_DFI, &dfi);
+ GetCtrlVal(p1,P1_FI0, &fi0);
+ GetCtrlVal(p1,P1_NFI, &nsteps);
+ GetCtrlVal(p1,P1_DAQTIME, &daqtime);
+ GetCtrlVal(p1,P1_FILENAME, filename);
+ GetCtrlVal(p1,P1_DEBUG, &debug);
+ sprintf(fname,"%s.dat", filename);
+ Geometry("config.ini");
+ HistogramsInit();
+ PETProjDataMgrFree();
+ PETProjDataMgrInit();
+
+ SetDebug(debug);
+ SetRingDiameter(2*conf.module[0].r);
+
+
+
+
+ FILE *fp = fopen(fname, "wb");
+ int runhdr[7]= {nsteps, Fi2Pos(fi0), Fi2Pos(dfi), 1,0,1, -daqtime};
+ addheader(fp,RUNREC_ID, runhdr );
+
+
+ vmconnect();
+ for (int i=0; i<nsteps; i++) {
+
+ double fi=fi0+dfi*i;
+ conf.rotation = fi*Pi()/180;
+#ifdef uSMC_USB
+ SetCtrlVal(p1,P1_STAGELED,1);
+ int SetPos = Fi2Pos(fi);
+ uSMC_MoveTo(node,SetPos);
+ uSMC_GetPosition (node,&CurPos);
+ SetCtrlVal(p1,P1_CPOSITION,CurPos);
+
+ SetCtrlVal(p1,P1_STAGELED,0);
+ if (ctrl_c) break;
+ int poshdr[6]= {i, CurPos, SetPos, 0,0,0};
+ addheader(fp, POSREC_ID, poshdr );
+ vmacquire(-daqtime,fp, NULL);
+#endif /* uSMC_USB */
+ SetCtrlVal(p1,P1_CEVE,i);
+ printf("Step =%d angle =%f\n", i, fi);
+
+ }
+ addheader(fp, ENDREC_ID, NULL);
+ fclose(fp);
+ daq_on = 1;
+
+ WriteInterfile(filename);
+ sprintf(fname,"%s.sroot", filename);
+ HistogramsWrite(fname);
+
+
+ vmdisconnect();
+ return 0;
+}
+
+
+
+int CVICALLBACK StartCB (int panel, int control, int event,
+ void *callbackData, int eventData1, int eventData2) {
+ ThreadFunctionPtr mythread = NULL;
+ switch (event) {
+
+ case EVENT_COMMIT:
+
+ controlID=0;
+ if (panel == p1 && control == P1_START) {
+ mythread = mdaq;
+ controlID= control;
+ }
+ if (mythread!=NULL) {
+ printf("New Thread panel=%d button=%d\n", panel, control);
+
+ // SetDimming(1);
+
+
+ if (!daq_on) CmtScheduleThreadPoolFunctionAdv (poolHandle, mythread, &controlID,
+ DEFAULT_THREAD_PRIORITY,
+ EndOfThread,
+ EVENT_TP_THREAD_FUNCTION_END,
+ NULL, RUN_IN_SCHEDULED_THREAD,
+ &tfID);
+ }
+ break;
+ }
+ return 0;
+}
+
+int CVICALLBACK StopCB (int panel, int control, int event,
+ void *callbackData, int eventData1, int eventData2) {
+ switch (event) {
+ case EVENT_COMMIT:
+ ctrl_c=1;
+ break;
+ }
+ return 0;
+}
+
+int CVICALLBACK DebugCB (int panel, int control, int event,
+ void *callbackData, int eventData1, int eventData2) {
+ switch (event) {
+ case EVENT_COMMIT:
+ GetCtrlVal(p1,P1_DEBUG, &debug);
+ break;
+ }
+ return 0;
+}
+
+
+
+int CVICALLBACK RedrawCB (int panel, int control, int event,
+ void *callbackData, int eventData1, int eventData2) {
+ int pmtid;
+ switch (event) {
+ case EVENT_COMMIT:
+ case EVENT_TIMER_TICK:
+ GetCtrlVal(p1,P1_CH, &gadcchannel);
+ pmtid = gadcchannel/16;
+ if (H1D_Exist(m_Adc[gadcchannel])) H1D_Draw(m_Adc[gadcchannel], p1,P1_GADC, &plothandle[0]);
+ if (H2D_Exist( m_CenterOfGravity[pmtid])) H2D_Draw(m_CenterOfGravity[pmtid], p1,P1_GXY, &plothandle[1]);
+ if (H1D_Exist(m_AdcSum[pmtid])) H1D_Draw(m_AdcSum[pmtid], p1,P1_GSUMADC, &plothandle[2]);
+ break;
+ }
+ return 0;
+}
+
Index: apps/petdemo/petdemo.prj
===================================================================
--- apps/petdemo/petdemo.prj (nonexistent)
+++ apps/petdemo/petdemo.prj (revision 264)
@@ -0,0 +1,617 @@
+[Project Header]
+Version = 1302
+Pathname = "/c/home/cvi/apps/petdemo/petdemo.prj"
+CVI Dir = "/c/program files/national instruments/cvi2013"
+CVI Shared Dir = "/C/Program Files/National Instruments/Shared/CVI"
+CVI Pub Local Dir = "/C/ProgramData/National Instruments/CVI2013"
+CVI Pub Global Dir = "/C/ProgramData/National Instruments/CVI"
+IVI Standard Root Dir = "/C/Program Files/IVI Foundation/IVI"
+VXIplug&play Framework Dir = "/C/Program Files/IVI Foundation/VISA/winnt"
+IVI Standard Root 64-bit Dir = "/C/Program Files/IVI Foundation/IVI"
+VXIplug&play Framework 64-bit Dir = "/C/Program Files/IVI Foundation/VISA/win64"
+Number of Files = 17
+Target Type = "Executable"
+Flags = 2064
+Copied From Locked InstrDrv Directory = False
+Copied from VXIPNP Directory = False
+Locked InstrDrv Name = ""
+Don't Display Deploy InstrDrv Dialog = False
+
+[Folders]
+Folder 0 = "User Interface Files"
+FolderEx 0 = "User Interface Files"
+Folder 1 = "Source Files"
+FolderEx 1 = "Source Files"
+Folder 2 = "Instrument Files"
+FolderEx 2 = "Instrument Files"
+Folder 3 = "Library Files"
+FolderEx 3 = "Library Files"
+Folder 4 = "Include Files"
+FolderEx 4 = "Include Files"
+
+[File 0001]
+File Type = "User Interface Resource"
+Res Id = 1
+Path Is Rel = True
+Path Rel To = "Project"
+Path Rel Path = "petdemo_uir.uir"
+Path = "/c/home/cvi/apps/petdemo/petdemo_uir.uir"
+Exclude = False
+Project Flags = 0
+Folder = "User Interface Files"
+Folder Id = 0
+
+[File 0002]
+File Type = "CSource"
+Res Id = 2
+Path Is Rel = True
+Path Rel To = "Project"
+Path Rel Path = "analyse.c"
+Path = "/c/home/cvi/apps/petdemo/analyse.c"
+Exclude = False
+Compile Into Object File = False
+Project Flags = 0
+Folder = "Source Files"
+Folder Id = 1
+
+[File 0003]
+File Type = "CSource"
+Res Id = 3
+Path Is Rel = True
+Path Rel To = "Project"
+Path Rel Path = "petdemo.c"
+Path = "/c/home/cvi/apps/petdemo/petdemo.c"
+Exclude = False
+Compile Into Object File = False
+Project Flags = 0
+Folder = "Source Files"
+Folder Id = 1
+
+[File 0004]
+File Type = "CSource"
+Res Id = 4
+Path Is Rel = True
+Path Rel To = "Project"
+Path Rel Path = "PETProjDataMgr.c"
+Path = "/c/home/cvi/apps/petdemo/PETProjDataMgr.c"
+Exclude = False
+Compile Into Object File = False
+Project Flags = 0
+Folder = "Source Files"
+Folder Id = 1
+
+[File 0005]
+File Type = "CSource"
+Res Id = 5
+Path Is Rel = True
+Path Rel To = "Project"
+Path Rel Path = "vmusb.c"
+Path = "/c/home/cvi/apps/petdemo/vmusb.c"
+Exclude = False
+Compile Into Object File = False
+Project Flags = 0
+Folder = "Source Files"
+Folder Id = 1
+
+[File 0006]
+File Type = "CSource"
+Res Id = 6
+Path Is Rel = True
+Path Rel To = "Project"
+Path Rel Path = "../../instr/WUSBVME_DLL/wusbvme_dll.c"
+Path = "/c/home/cvi/instr/WUSBVME_DLL/wusbvme_dll.c"
+Exclude = False
+Compile Into Object File = False
+Project Flags = 0
+Folder = "Source Files"
+Folder Id = 1
+
+[File 0007]
+File Type = "Function Panel"
+Res Id = 7
+Path Is Rel = True
+Path Rel To = "Project"
+Path Rel Path = "../../instr/HISTO/H1D.fp"
+Path = "/c/home/cvi/instr/HISTO/H1D.fp"
+Exclude = False
+Project Flags = 0
+Folder = "Instrument Files"
+Folder Id = 2
+
+[File 0008]
+File Type = "Function Panel"
+Res Id = 8
+Path Is Rel = True
+Path Rel To = "Project"
+Path Rel Path = "../../instr/HISTO/H2D.fp"
+Path = "/c/home/cvi/instr/HISTO/H2D.fp"
+Exclude = False
+Project Flags = 0
+Folder = "Instrument Files"
+Folder Id = 2
+
+[File 0009]
+File Type = "Function Panel"
+Res Id = 9
+Path Is Rel = True
+Path Rel To = "Project"
+Path Rel Path = "../../instr/HISTO/H3D.fp"
+Path = "/c/home/cvi/instr/HISTO/H3D.fp"
+Exclude = False
+Project Flags = 0
+Folder = "Instrument Files"
+Folder Id = 2
+
+[File 0010]
+File Type = "Function Panel"
+Res Id = 10
+Path Is Rel = True
+Path Rel To = "CVI"
+Path Rel To Override = "CVI"
+Path Rel Path = "toolslib/toolbox/inifile.fp"
+Path = "/c/Program Files/National Instruments/CVI2013/toolslib/toolbox/inifile.fp"
+Exclude = False
+Project Flags = 0
+Folder = "Instrument Files"
+Folder Id = 2
+
+[File 0011]
+File Type = "Function Panel"
+Res Id = 11
+Path Is Rel = True
+Path Rel To = "Project"
+Path Rel Path = "../../instr/uSMC/uSMC.fp"
+Path = "/c/home/cvi/instr/uSMC/uSMC.fp"
+Exclude = False
+Project Flags = 0
+Folder = "Instrument Files"
+Folder Id = 2
+
+[File 0012]
+File Type = "Library"
+Res Id = 12
+Path Is Rel = True
+Path Rel To = "Project"
+Path Rel Path = "../../instr/vme/msvc/libusb.lib"
+Path = "/c/home/cvi/instr/vme/msvc/libusb.lib"
+Exclude = False
+Project Flags = 0
+Folder = "Library Files"
+Folder Id = 3
+
+[File 0013]
+File Type = "Library"
+Res Id = 13
+Path Is Rel = True
+Path Rel To = "Project"
+Path Rel Path = "libxxusb.lib"
+Path = "/c/home/cvi/apps/petdemo/libxxusb.lib"
+Exclude = False
+Project Flags = 0
+Folder = "Library Files"
+Folder Id = 3
+
+[File 0014]
+File Type = "Library"
+Res Id = 14
+Path Is Rel = True
+Path Rel To = "Project"
+Path Rel Path = "../../instr/uSMC/msvc/USMCDLL.lib"
+Path = "/c/home/cvi/instr/uSMC/msvc/USMCDLL.lib"
+Exclude = False
+Project Flags = 0
+Folder = "Library Files"
+Folder Id = 3
+
+[File 0015]
+File Type = "Include"
+Res Id = 15
+Path Is Rel = True
+Path Rel To = "Project"
+Path Rel Path = "daq.h"
+Path = "/c/home/cvi/apps/petdemo/daq.h"
+Exclude = False
+Project Flags = 0
+Folder = "Include Files"
+Folder Id = 4
+
+[File 0016]
+File Type = "Include"
+Res Id = 16
+Path Is Rel = True
+Path Rel To = "Project"
+Path Rel Path = "PETProjDataMgr.h"
+Path = "/c/home/cvi/apps/petdemo/PETProjDataMgr.h"
+Exclude = False
+Project Flags = 0
+Folder = "Include Files"
+Folder Id = 4
+
+[File 0017]
+File Type = "Include"
+Res Id = 17
+Path Is Rel = True
+Path Rel To = "Project"
+Path Rel Path = "../../instr/vme/vme.h"
+Path = "/c/home/cvi/instr/vme/vme.h"
+Exclude = False
+Project Flags = 0
+Folder = "Include Files"
+Folder Id = 4
+
+[Custom Build Configs]
+Num Custom Build Configs = 0
+
+[Default Build Config Debug]
+Config Name = "Debug"
+Is 64-Bit = False
+Is Release = False
+Default Calling Convention = "cdecl"
+Optimization Level = "Optimize for speed (level 2)"
+Require Prototypes = True
+Show Warning IDs in Build Output = False
+Selected Warning Level = "None"
+Warning List None = "4,9,84,105,106,107,108,109,110,111"
+Warning List Common = ""
+Warning List Extended = ""
+Warning List All = ""
+Warning Mode = 0
+Enable Unreferenced Identifiers Warning = False
+Enable Pointer Mismatch Warning = False
+Enable Unreachable Code Warning = False
+Enable Assignment In Conditional Warning = False
+Uninitialized Locals Compile Warning = "Aggressive"
+Require Return Values = True
+Enable C99 Extensions = False
+Enable OpenMP Extensions = False
+Stack Size = 250000
+Stack Reserve = 1048576
+Stack Commit = 4096
+Image Base Address = 4194304
+Image Base Address x64 = 4194304
+Compiler Defines = "/DWIN32_LEAN_AND_MEAN /DXXDLL"
+Sign = False
+Sign Store = ""
+Sign Certificate = ""
+Sign Timestamp URL = ""
+Sign URL = ""
+Manifest Embed = False
+Icon File Is Rel = False
+Icon File = ""
+Application Title = ""
+Use IVI Subdirectories for Import Libraries = False
+Use VXIPNP Subdirectories for Import Libraries = False
+Use Dflt Import Lib Base Name = True
+Where to Copy DLL = "Do not copy"
+Custom Directory to Copy DLL Is Rel = False
+Custom Directory to Copy DLL = ""
+Generate Source Documentation = "None"
+Runtime Support = "Full Runtime Support"
+Runtime Binding = "Shared"
+Embed Project .UIRs = False
+Generate Map File = False
+Embed Timestamp = True
+Create Console Application = False
+Using LoadExternalModule = False
+DLL Exports = "Include File Symbols"
+Register ActiveX Server = False
+Numeric File Version = "1,0,0,0"
+Numeric Prod Version = "1,0,0,0"
+Comments = ""
+Comments Ex = ""
+Company Name = ""
+Company Name Ex = "%company"
+File Description = "petdemo (Debug x86)"
+File Description Ex = "%application (%rel_dbg %arch)"
+File Version = "1.0"
+File Version Ex = "%f1.%f2"
+Internal Name = "petdemo"
+Internal Name Ex = "%basename"
+Legal Copyright = "Copyright © 2017"
+Legal Copyright Ex = "Copyright © %company %Y"
+Legal Trademarks = ""
+Legal Trademarks Ex = ""
+Original Filename = "petdemo.exe"
+Original Filename Ex = "%filename"
+Private Build = ""
+Private Build Ex = ""
+Product Name = " petdemo"
+Product Name Ex = "%company %application"
+Product Version = "1.0"
+Product Version Ex = "%p1.%p2"
+Special Build = ""
+Special Build Ex = ""
+Add Type Lib To DLL = False
+Include Type Lib Help Links = False
+TLB Help Style = "HLP"
+Type Lib FP File Is Rel = False
+Type Lib FP File = ""
+
+[Default Build Config Release]
+Config Name = "Release"
+Is 64-Bit = False
+Is Release = True
+Default Calling Convention = "cdecl"
+Optimization Level = "Optimize for speed (level 2)"
+Require Prototypes = True
+Show Warning IDs in Build Output = False
+Selected Warning Level = "None"
+Warning List None = "4,9,84,105,106,107,108,109,110,111"
+Warning List Common = ""
+Warning List Extended = ""
+Warning List All = ""
+Warning Mode = 0
+Enable Unreferenced Identifiers Warning = False
+Enable Pointer Mismatch Warning = False
+Enable Unreachable Code Warning = False
+Enable Assignment In Conditional Warning = False
+Uninitialized Locals Compile Warning = "Aggressive"
+Require Return Values = True
+Enable C99 Extensions = False
+Enable OpenMP Extensions = False
+Stack Size = 250000
+Stack Reserve = 1048576
+Stack Commit = 4096
+Image Base Address = 4194304
+Image Base Address x64 = 4194304
+Compiler Defines = "/DWIN32_LEAN_AND_MEAN"
+Sign = False
+Sign Store = ""
+Sign Certificate = ""
+Sign Timestamp URL = ""
+Sign URL = ""
+Manifest Embed = False
+Icon File Is Rel = False
+Icon File = ""
+Application Title = ""
+Use IVI Subdirectories for Import Libraries = False
+Use VXIPNP Subdirectories for Import Libraries = False
+Use Dflt Import Lib Base Name = True
+Where to Copy DLL = "Do not copy"
+Custom Directory to Copy DLL Is Rel = False
+Custom Directory to Copy DLL = ""
+Generate Source Documentation = "None"
+Runtime Support = "Full Runtime Support"
+Runtime Binding = "Shared"
+Embed Project .UIRs = False
+Generate Map File = False
+Embed Timestamp = True
+Create Console Application = False
+Using LoadExternalModule = False
+DLL Exports = "Include File Symbols"
+Register ActiveX Server = False
+Add Type Lib To DLL = False
+Include Type Lib Help Links = False
+TLB Help Style = "HLP"
+Type Lib FP File Is Rel = False
+Type Lib FP File = ""
+
+[Default Build Config Debug64]
+Config Name = "Debug64"
+Is 64-Bit = True
+Is Release = False
+Default Calling Convention = "cdecl"
+Optimization Level = "Optimize for speed (level 2)"
+Require Prototypes = True
+Show Warning IDs in Build Output = False
+Selected Warning Level = "None"
+Warning List None = "4,9,84,105,106,107,108,109,110,111"
+Warning List Common = ""
+Warning List Extended = ""
+Warning List All = ""
+Warning Mode = 0
+Enable Unreferenced Identifiers Warning = False
+Enable Pointer Mismatch Warning = False
+Enable Unreachable Code Warning = False
+Enable Assignment In Conditional Warning = False
+Uninitialized Locals Compile Warning = "Aggressive"
+Require Return Values = True
+Enable C99 Extensions = False
+Enable OpenMP Extensions = False
+Stack Size = 250000
+Stack Reserve = 1048576
+Stack Commit = 4096
+Image Base Address = 4194304
+Image Base Address x64 = 4194304
+Compiler Defines = "/DWIN32_LEAN_AND_MEAN"
+Sign = False
+Sign Store = ""
+Sign Certificate = ""
+Sign Timestamp URL = ""
+Sign URL = ""
+Manifest Embed = False
+Icon File Is Rel = False
+Icon File = ""
+Application Title = ""
+Use IVI Subdirectories for Import Libraries = False
+Use VXIPNP Subdirectories for Import Libraries = False
+Use Dflt Import Lib Base Name = True
+Where to Copy DLL = "Do not copy"
+Custom Directory to Copy DLL Is Rel = False
+Custom Directory to Copy DLL = ""
+Generate Source Documentation = "None"
+Runtime Support = "Full Runtime Support"
+Runtime Binding = "Shared"
+Embed Project .UIRs = False
+Generate Map File = False
+Embed Timestamp = True
+Create Console Application = False
+Using LoadExternalModule = False
+DLL Exports = "Include File Symbols"
+Register ActiveX Server = False
+Add Type Lib To DLL = False
+Include Type Lib Help Links = False
+TLB Help Style = "HLP"
+Type Lib FP File Is Rel = False
+Type Lib FP File = ""
+
+[Default Build Config Release64]
+Config Name = "Release64"
+Is 64-Bit = True
+Is Release = True
+Default Calling Convention = "cdecl"
+Optimization Level = "Optimize for speed (level 2)"
+Require Prototypes = True
+Show Warning IDs in Build Output = False
+Selected Warning Level = "None"
+Warning List None = "4,9,84,105,106,107,108,109,110,111"
+Warning List Common = ""
+Warning List Extended = ""
+Warning List All = ""
+Warning Mode = 0
+Enable Unreferenced Identifiers Warning = False
+Enable Pointer Mismatch Warning = False
+Enable Unreachable Code Warning = False
+Enable Assignment In Conditional Warning = False
+Uninitialized Locals Compile Warning = "Aggressive"
+Require Return Values = True
+Enable C99 Extensions = False
+Enable OpenMP Extensions = False
+Stack Size = 250000
+Stack Reserve = 1048576
+Stack Commit = 4096
+Image Base Address = 4194304
+Image Base Address x64 = 4194304
+Compiler Defines = "/DWIN32_LEAN_AND_MEAN"
+Sign = False
+Sign Store = ""
+Sign Certificate = ""
+Sign Timestamp URL = ""
+Sign URL = ""
+Manifest Embed = False
+Icon File Is Rel = False
+Icon File = ""
+Application Title = ""
+Use IVI Subdirectories for Import Libraries = False
+Use VXIPNP Subdirectories for Import Libraries = False
+Use Dflt Import Lib Base Name = True
+Where to Copy DLL = "Do not copy"
+Custom Directory to Copy DLL Is Rel = False
+Custom Directory to Copy DLL = ""
+Generate Source Documentation = "None"
+Runtime Support = "Full Runtime Support"
+Runtime Binding = "Shared"
+Embed Project .UIRs = False
+Generate Map File = False
+Embed Timestamp = True
+Create Console Application = False
+Using LoadExternalModule = False
+DLL Exports = "Include File Symbols"
+Register ActiveX Server = False
+Add Type Lib To DLL = False
+Include Type Lib Help Links = False
+TLB Help Style = "HLP"
+Type Lib FP File Is Rel = False
+Type Lib FP File = ""
+
+[Compiler Options]
+Default Calling Convention = "cdecl"
+Require Prototypes = True
+Require Return Values = True
+Enable Pointer Mismatch Warning = False
+Enable Unreachable Code Warning = False
+Enable Unreferenced Identifiers Warning = False
+Enable Assignment In Conditional Warning = False
+O Option Compatible With 5.0 = False
+Enable C99 Extensions = False
+Uninitialized Locals Compile Warning = "Aggressive"
+Precompile Prefix Header = False
+Prefix Header File = ""
+
+[Run Options]
+Stack Size = 250000
+Stack Commit = 4096
+Image Base Address = 4194304
+Image Base Address x64 = 4194304
+
+[Compiler Defines]
+Compiler Defines = "/DWIN32_LEAN_AND_MEAN"
+
+[Include Paths]
+Include Path 1 Is Rel = True
+Include Path 1 Rel To = "Project"
+Include Path 1 Rel Path = "../../instr/CAENV965"
+Include Path 1 = "/c/home/cvi/instr/CAENV965"
+Include Path 2 Is Rel = True
+Include Path 2 Rel To = "Project"
+Include Path 2 Rel Path = "../../instr/WUSBVME_DLL"
+Include Path 2 = "/c/home/cvi/instr/WUSBVME_DLL"
+
+[Create Executable]
+Executable File_Debug Is Rel = True
+Executable File_Debug Rel To = "Project"
+Executable File_Debug Rel Path = "petdemo.exe"
+Executable File_Debug = "/c/home/cvi/apps/petdemo/petdemo.exe"
+Executable File_Release Is Rel = True
+Executable File_Release Rel To = "Project"
+Executable File_Release Rel Path = "petdemo.exe"
+Executable File_Release = "/c/home/cvi/apps/petdemo/petdemo.exe"
+Executable File_Debug64 Is Rel = True
+Executable File_Debug64 Rel To = "Project"
+Executable File_Debug64 Rel Path = "petdemo.exe"
+Executable File_Debug64 = "/c/home/cvi/apps/petdemo/petdemo.exe"
+Executable File_Release64 Is Rel = True
+Executable File_Release64 Rel To = "Project"
+Executable File_Release64 Rel Path = "petdemo.exe"
+Executable File_Release64 = "/c/home/cvi/apps/petdemo/petdemo.exe"
+Icon File Is Rel = False
+Icon File = ""
+Application Title = ""
+DLL Exports = "Include File Symbols"
+Use IVI Subdirectories for Import Libraries = False
+Use VXIPNP Subdirectories for Import Libraries = False
+Use Dflt Import Lib Base Name = True
+Where to Copy DLL = "Do not copy"
+Custom Directory to Copy DLL Is Rel = False
+Custom Directory to Copy DLL = ""
+Generate Source Documentation = "None"
+Add Type Lib To DLL = False
+Include Type Lib Help Links = False
+TLB Help Style = "HLP"
+Type Lib FP File Is Rel = False
+Type Lib FP File = ""
+Type Lib Guid = ""
+Runtime Support = "Full Runtime Support"
+Instrument Driver Support Only = False
+Embed Project .UIRs = False
+Generate Map File = False
+
+[External Compiler Support]
+UIR Callbacks File Option = 0
+Using LoadExternalModule = False
+Create Project Symbols File = True
+UIR Callbacks Obj File Is Rel = False
+UIR Callbacks Obj File = ""
+Project Symbols H File Is Rel = False
+Project Symbols H File = ""
+Project Symbols Obj File Is Rel = False
+Project Symbols Obj File = ""
+
+[ActiveX Server Options]
+Specification File Is Rel = False
+Specification File = ""
+Source File Is Rel = False
+Source File = ""
+Include File Is Rel = False
+Include File = ""
+IDL File Is Rel = False
+IDL File = ""
+Register ActiveX Server = False
+
+[Signing Info]
+Sign = False
+Sign Debug Build = False
+Store = ""
+Certificate = ""
+Timestamp URL = ""
+URL = ""
+
+[Manifest Info]
+Embed = False
+
+[tpcSection]
+tpcEnabled = 0
+tpcOverrideEnvironment = 0
+tpcEnabled x64 = 0
+tpcOverrideEnvironment x64 = 0
+
Index: apps/petdemo/petdemo_uir.h
===================================================================
--- apps/petdemo/petdemo_uir.h (nonexistent)
+++ apps/petdemo/petdemo_uir.h (revision 264)
@@ -0,0 +1,61 @@
+/**************************************************************************/
+/* LabWindows/CVI User Interface Resource (UIR) Include File */
+/* */
+/* WARNING: Do not add to, delete from, or otherwise modify the contents */
+/* of this include file. */
+/**************************************************************************/
+
+#include <userint.h>
+
+#ifdef __cplusplus
+ extern "C" {
+#endif
+
+ /* Panels and Controls: */
+
+#define P1 1
+#define P1_EXIT 2 /* control type: command, callback function: ExitCB */
+#define P1_STOP 3 /* control type: command, callback function: StopCB */
+#define P1_START 4 /* control type: command, callback function: StartCB */
+#define P1_CPOSITION 5 /* control type: numeric, callback function: (none) */
+#define P1_NFI 6 /* control type: numeric, callback function: (none) */
+#define P1_DFI 7 /* control type: numeric, callback function: (none) */
+#define P1_FI0 8 /* control type: numeric, callback function: (none) */
+#define P1_CEVE 9 /* control type: numeric, callback function: (none) */
+#define P1_DAQTIME 10 /* control type: numeric, callback function: (none) */
+#define P1_GXY 11 /* control type: graph, callback function: (none) */
+#define P1_GSUMADC 12 /* control type: graph, callback function: (none) */
+#define P1_GADC 13 /* control type: graph, callback function: (none) */
+#define P1_GSINOGRAM 14 /* control type: graph, callback function: (none) */
+#define P1_PROGRESS 15 /* control type: slide, callback function: (none) */
+#define P1_STAGELED 16 /* control type: LED, callback function: (none) */
+#define P1_STDIO 17 /* control type: textBox, callback function: (none) */
+#define P1_FILENAME 18 /* control type: string, callback function: (none) */
+#define P1_DEBUG 19 /* control type: radioButton, callback function: DebugCB */
+#define P1_CH 20 /* control type: numeric, callback function: RedrawCB */
+#define P1_TIMER 21 /* control type: timer, callback function: RedrawCB */
+#define P1_DAQRATE 22 /* control type: strip, callback function: (none) */
+
+
+ /* Control Arrays: */
+
+ /* (no control arrays in the resource file) */
+
+
+ /* Menu Bars, Menus, and Menu Items: */
+
+ /* (no menu bars in the resource file) */
+
+
+ /* Callback Prototypes: */
+
+int CVICALLBACK DebugCB(int panel, int control, int event, void *callbackData, int eventData1, int eventData2);
+int CVICALLBACK ExitCB(int panel, int control, int event, void *callbackData, int eventData1, int eventData2);
+int CVICALLBACK RedrawCB(int panel, int control, int event, void *callbackData, int eventData1, int eventData2);
+int CVICALLBACK StartCB(int panel, int control, int event, void *callbackData, int eventData1, int eventData2);
+int CVICALLBACK StopCB(int panel, int control, int event, void *callbackData, int eventData1, int eventData2);
+
+
+#ifdef __cplusplus
+ }
+#endif
Index: apps/petdemo/petdemo_uir.uir
===================================================================
Cannot display: file marked as a binary type.
svn:mime-type = application/octet-stream
/apps/petdemo/petdemo_uir.uir
Property changes:
Added: svn:mime-type
## -0,0 +1 ##
+application/octet-stream
\ No newline at end of property
Index: apps/petdemo/vmusb.c
===================================================================
--- apps/petdemo/vmusb.c (nonexistent)
+++ apps/petdemo/vmusb.c (revision 264)
@@ -0,0 +1,433 @@
+//#define XXDLL
+//#include "vme.h"
+
+
+
+typedef unsigned short ADDRESS_MODIFIER;
+#define Std_NoPriv_Data (ADDRESS_MODIFIER)0x39
+
+#include "wusbvme_dll.h"
+#include "toolbox.h"
+#include <stdio.h>
+#include <stdlib.h>
+#include <string.h>
+
+
+#include <errno.h>
+#include <ctype.h>
+#include <time.h>
+
+
+#include <signal.h>
+
+
+#include "CAENV965_DEF.h"
+
+#include "petdemo_uir.h"
+#include "daq.h"
+
+extern int p1;
+extern int ctrl_c;
+
+
+#define VERSION 1.0
+#define TIMEOUT 3
+
+int *fThreshold = NULL;
+
+/* VME modules */
+
+#define CAEN_V792 0x340000 // IJS V792
+#define CAEN_V792_1 0x530000 // FMF1 V792
+#define CAEN_V792_2 0x630000 // FMF2 V792
+#define CAEN_V965 0x350000 // IJS V965
+int addr[3]= {CAEN_V792,CAEN_V792_1,CAEN_V965 };
+const int nadc=2;
+#define BUFF_L 2048
+static int stackwrite[10000];
+static int stackdata[10000],stackdump[27000];
+
+
+/************************************************/
+int weight_while(int num) {
+ int i, tmp;
+
+ for ( i =0; i <num; i++ ) tmp = 0;
+ return 0;
+}
+#define WWHILE weight_while(0)
+
+#define TRUE 1
+#define FALSE 0
+
+int timer_out;
+
+void timerast (int signumber) {
+ timer_out = TRUE;
+ fprintf(stderr,"TIMEOUT !!!\n");
+}
+
+void tmlnk (int tout) {
+
+}
+
+void tmulk () {
+
+}
+
+int fexist( char *path) {
+
+
+ if(!path || !*path) return 0;
+ ssize_t fileSize;
+ int res = FileExists (path, &fileSize);
+
+
+ if (res>0) {
+ return 1;
+ } else {
+ return 0;
+ }
+}
+
+
+
+
+int init() {
+
+ xxusb_register_write(WUSB_udev,1,0x0); // Stop DAQ mode
+ while (xxusb_usbfifo_read(WUSB_udev,(int *) stackdump,BUFF_L,100)>0);
+
+
+ int rate=1000;
+ int freq=80000000/rate-40; // 80 MHz
+ if (freq<72) freq=72;
+ // Set DGG channel A as a pulser, output on O1,
+ // with delay =500 x 12.5ns,
+ // and width = 500 x 12.5ns,
+ // not latching or inverting
+ // VME_DGG(udev,0,6,0,24000,6000,0,0);
+ VME_DGG(WUSB_udev,0,6,0,freq,40,0,0);
+ // Set DGG channel B to trigger on NIM1, output on O2,
+ // with delay =200 x 12.5ns,
+ // and width = 200 x 12.5ns,
+ // not latching or inverting
+ VME_DGG(WUSB_udev,1,1,1,0,10,0,1);
+ printf("CAEN V965 Pedestal set to %d\n", fPedestal);
+
+// INIT stackdata
+ int fPedestal=255;
+ WIENER_VMUSB_StackInit();
+
+ for (int i=0; i<nadc; i++) {
+ WIENER_VMUSB_StackWriteA24D16( addr[i] + CAENV965_CRN , 0x0);
+ WIENER_VMUSB_StackWriteA24D16( addr[i] + CAENV965_GEO , i);
+ WIENER_VMUSB_StackReadA24D16(addr[i] + CAENV965_GEO);
+ for (int j=0; j<32; j++) {
+ WIENER_VMUSB_StackWriteA24D16(addr[i] + CAENV965_THM + 0x02*j, fThreshold[j+i*32]); // threshold/kill for 32 channels
+ }
+ WIENER_VMUSB_StackWriteA24D16( addr[i] + CAENV965_BS1, 0x80 ); // soft reset
+ WIENER_VMUSB_StackWriteA24D16( addr[i] + CAENV965_BC1, 0x80 ); // soft reset
+ WIENER_VMUSB_StackWriteA24D16( addr[i] + CAENV965_PED, fPedestal ); // pedestal
+ WIENER_VMUSB_StackWriteA24D16( addr[i] + CAENV965_BS2,0x5000);
+ WIENER_VMUSB_StackWriteA24D16( addr[i] + CAENV965_BS2,0x4); // clear module
+ WIENER_VMUSB_StackWriteA24D16( addr[i] + CAENV965_BC2,0x4);
+ }
+ WIENER_VMUSB_StackMarker(0xFAFC);
+ //WIENER_VMUSB_StackPrint();
+
+
+ VME_LED_settings(WUSB_udev, 0,0,0,0); // Set Yellow LED to light with with USB out buffer not empty
+ VME_LED_settings(WUSB_udev, 1,1,0,0); // Set Red LED to light with NIM1
+ VME_LED_settings(WUSB_udev,2,0,0,0); // Set Green LED to light when stack is not empty
+
+ unsigned int vmereg;
+ VME_register_read(WUSB_udev,0x00,&vmereg);
+ printf("VMUSB Firmware ID -> 0x%08X\n",vmereg);
+
+ VME_register_read(WUSB_udev,0x04,&vmereg);
+ printf("VMUSB Global Mode -> 0x%08X\n",vmereg);
+
+ vmereg=(vmereg&0xF000)|0x0004;
+ VME_register_write(WUSB_udev,0x04,vmereg);
+ VME_register_write(WUSB_udev,0x08,0x00000080);
+ VME_register_write(WUSB_udev,0x28,0x0);
+ VME_register_write(WUSB_udev,0x2C,0x0);
+ VME_register_write(WUSB_udev,0x30,0x0);
+ VME_register_write(WUSB_udev,0x34,0x0);
+ VME_register_write(WUSB_udev,0x3C,0x000);
+
+ int nb = WIENER_VMUSB_StackGetUint32(10000,stackdata);
+ int ret= xxusb_stack_execute(WUSB_udev,(uint32_t *)stackdata);
+ printf("Init::%d ret=%d\n",nb,ret);
+ if (ret>0) for (int k=0; k<ret/2; k++) printf ("stackdata=0x%08X\n",stackdata[k]);
+
+
+ WIENER_VMUSB_StackClear();
+ WIENER_VMUSB_StackMarker(0xFFAB);
+ for (int j=0; j<36; j++) WIENER_VMUSB_StackReadA24D32(addr[0] + CAENV965_OB);
+ WIENER_VMUSB_StackMarker(0xFAFB);
+ for (int k0=0; k0<nadc; k0++) {
+ WIENER_VMUSB_StackWriteA24D16(addr[k0] + CAENV965_BS2,0x4); // clear module
+ WIENER_VMUSB_StackWriteA24D16(addr[k0] + CAENV965_BC2,0x4);
+ }
+
+ int nb0= WIENER_VMUSB_StackGetUint32(10000,&stackwrite[0]);
+ if (nb0>768) {
+ fprintf(stderr,"nb0=%d > 768 error xxusb_stack_write\n", nb0);
+ exit(-1);
+ }
+ nb =xxusb_stack_write(WUSB_udev,0x2,(uint32_t *) stackwrite);
+ nb0=xxusb_stack_read(WUSB_udev,0x2,(uint32_t *) stackdata);
+ for (int k1=0; k1<stackwrite[0]+1; k1++) {
+ if (stackdata[k1]!=stackwrite[k1]) printf("%d %d init err %x %x\n",nb,nb0,
+ stackwrite[k1], stackdata[k1]);
+ }
+
+ if (fMode==2) xxusb_register_write(WUSB_udev,1,0x1); // Start DAQ mode
+
+ printf("daq::init() \n");
+ return 0;
+}
+
+int vmconnect() {
+ WIENER_VMUSB_VME_START("VM0120");
+ //VME_START(WIENER_VMUSB);
+
+ printf("daq::connect()\n");
+ return 0;
+}
+
+int vmdisconnect() {
+ /* zakljuci */
+ WIENER_VMUSB_VME_STOP();
+ printf("daq::disconnect()\n");
+ return 0;
+}
+
+int clear() {
+ return 0;
+}
+
+int module_header(int recid,uint32_t *data,int len) {
+ data[0] = recid;
+ data[1] = (len >0)? len : 0 ;
+ return data[1]+2;
+}
+
+
+int event(unsigned int *data, int maxn, int *ctr, int print) {
+ int tout=200; /* 1/100 of a second */
+ const int lsize=sizeof(uint32_t);
+
+
+ ctr[0]++;
+ ctr[1]++;
+
+ int count=0;
+ switch (fMode) {
+ case 0: { // normal calls
+ unsigned short clr= 0x4;
+ unsigned int status=0;
+ uint32_t mdata;
+
+ for (int i=0; i<nadc; i++) {
+ // wait for trg
+ tmlnk (tout);
+ do WIENER_VMUSB_VME_A24D16_R( addr[i] + CAENV965_SR1, &status);
+ while ( (status&0x1)==0 && timer_out==0 && ctrl_c==0 );
+ tmulk();
+ // readout data
+ if (timer_out) return 0;
+ if (ctrl_c) return 0;
+ int len=0;
+
+ do {
+ WIENER_VMUSB_VME_A24D32_R(addr[i] + CAENV965_OB, &mdata);
+ mdata=data[count++];
+ len++;
+ } while ( (mdata & 0x4000000)==0 && timer_out==0) ; // bit 26 EOB or not valid datum
+ // clear
+ WIENER_VMUSB_VME_A24D16_W( addr[i] + CAENV965_BS2, &clr);
+ WIENER_VMUSB_VME_A24D16_W( addr[i] + CAENV965_BC2, &clr);
+
+ if (count+2<maxn) {
+ if (print) printf("V965 %3d\n",len);
+ count+=module_header(0x130+i,&data[count],len);
+ ctr[2]++;
+ ctr[3]+=len;
+ }
+
+ timer_out=0;
+ }
+ }
+ break;
+ case 1: { // stack execute
+ WIENER_VMUSB_StackGetUint32(10000,(int *)data);
+ int ret=xxusb_stack_execute(WUSB_udev,(uint32_t *) data); //The first element of the array is the number of bytes.
+ if (ret< 0 ) {
+ printf ("xxusb_stack_execute error err=%d\n",ret); \
+ count = 0;
+ } else count= ret/lsize;
+
+ }
+ break;
+ case 2: { // stack load
+ int ret=xxusb_usbfifo_read(WUSB_udev,(int *) data,BUFF_L,100);
+ if (ret< 0 ) {
+ if (ret!=-110) {
+ printf ("xxusb_usbfifo_read error err=%d\n",ret);
+ end();
+ init();
+ }
+ count = 0;
+ } else {
+ if (debug) {
+ for (int i=0; i<100; i++) {
+ printf ("%4d fifodata=0x%08X\n",i, data[i]);
+ if (data[i]==0xFAFB) break;
+ }
+ /*
+
+ 0 fifodata=0x0000000D
+ 1 fifodata=0x00000049
+ 2 fifodata=0x0000FFAB
+ 3 fifodata=0x00002000
+ 4 fifodata=0x00000200
+ 5 fifodata=0x00004141
+ 6 fifodata=0x00000000
+ 7 fifodata=0x00004057
+
+ 70 fifodata=0x00000400
+ 71 fifodata=0x00000035
+ 72 fifodata=0x00000600
+ 73 fifodata=0x00000035
+ 74 fifodata=0x0000FAFB
+
+
+ */
+ }
+ static double t0=0;
+ static int nc =0;
+ nc += data[0];
+ if (print) {
+ //printf("##------------------ret=%d data[0]=%d\n",ret,(int)data[0]);
+ double t1 = Timer();
+ double dt = t1-t0;
+ t0 = t1;
+ float rate = (dt>0)?nc/dt:0;
+
+ //printf("--------ret=%d data[0]=%f nc=%d dt=%f\n",ret,rate,nc,dt);
+ PlotStripChart (p1, P1_DAQRATE, &rate, 1, 0, 0,VAL_FLOAT);
+ nc =0;
+ }
+ count= ret/lsize;
+ ctr[2]+=data[0];
+ ctr[3]+=count;
+ }
+ }
+ break;
+ }
+ return count*lsize;
+}
+
+int end() {
+
+ xxusb_register_write(WUSB_udev,1,0x0); // Stop DAQ mode
+ while (xxusb_usbfifo_read(WUSB_udev,(int *) stackdata,BUFF_L,30)>0);
+ printf("daq::end()\n");
+ return 0;
+}
+
+
+
+
+int vmacquire (int neve,FILE *fp, char *fpedname) {
+ // neve .... negative argument time ( in s )limited event loop
+ if (fThreshold== NULL) fThreshold = (int *) malloc (sizeof(int)*128);
+ // print welcome message
+ time_t t, told;
+ double tstart, tstop;
+ time(&t);
+ fThresholdEnable=0;
+
+#define BSIZE 10000
+ uint32_t data[10000];
+
+ fMode = 2;
+ fPedestal=255;
+ for (int ki=0; ki<128; ki++) {
+ if (ki<72) fThreshold[ki]=0;
+ else fThreshold[ki]=0x1<<8; // samo 4 kanali na zadnjem modulu so enablani
+ }
+ fThresholdEnable=0;
+ fStop=0;
+
+
+
+ if (fpedname !=NULL) {
+
+ FILE *fped=fopen(fpedname,"r");
+
+ int j=0;
+ int ndim=400;
+ char line[ndim];
+ int val=0;
+ while (fgets(line,ndim,fped)!=NULL) {
+ sscanf(line,"%d",&val);
+ fThreshold[j++]=val;
+ }
+ fThresholdEnable=1;
+ fclose(fped);
+
+ }
+
+
+
+ init();
+
+
+ int hdr[4]= {2}; // recid
+ int i=0;
+ int ntotal=0;
+ int counters[30]= {0,0,0,0,0, 0,0,0,0,0,0,0};
+ char names[10][20]= {"TRG","CAEN V965"};
+ time(&t);
+ tstart=Timer();
+ told =t;
+ tstop=tstart+360000;
+ if (neve<-1) {
+ tstop=tstart-neve;
+ neve=-1;
+ }
+ for (i=0; i!=neve && !ctrl_c && Timer()<tstop; i++) {
+ time(&t);
+ if (t!=told ) printf("%d in %2.2f min daq::event() %s\n",i, (double)(Timer()-tstart)/60., ctime(&t));
+ int nb=event(data,BSIZE, counters,t!=told);
+ analyse(nb, data);
+ if (nb>0) {
+ // zapis v datoteko
+ hdr[1]=nb+4*sizeof(int);
+ hdr[2]=time(NULL);
+ hdr[3]=i;
+
+ fwrite(hdr, 1, sizeof(int)*4, fp); //gzip
+ ntotal += fwrite(data,1, nb, fp);
+ told=t;
+ } else i--;
+ }
+
+ end();
+
+ printf("Number of Events: %d\n",i);
+ if (ctrl_c) printf("User Program termination CTRL-C\n");
+ if (Timer()>tstop ) printf("Timeout termination tstart# t>tstop: %d# %d >%d\n",(int)t, (int)tstart, (int) tstop);
+
+
+ printf("%d bytes written to file\nCounts:\n", (int) (ntotal*sizeof(int)));
+ for (i=0; i<2; i++) printf("%s\t%d\t%d\n",names[i],counters[2*i],counters[2*i+1]) ;
+
+
+ return 0;
+}
+
Index: instr/CAENV965/CAENV965.c
===================================================================
--- instr/CAENV965/CAENV965.c (revision 263)
+++ instr/CAENV965/CAENV965.c (revision 264)
@@ -1,13 +1,10 @@
#include <ansi_c.h>
+
+
+#include "vme.h"
#include "CAENV965.h"
#include "CAENV965_DEF.h"
-#ifdef SISVME
-#include "sisvme_dll.h"
-#endif
-#ifdef WIENVME
-#include "wienvme_dll.h"
-#endif
static unsigned long ModuleAddress[10];
/cvi/instr/HISTO/H1D.c
11,7 → 11,7
 
#include "H1D.h"
 
#define H1DMAX 500
 
H1D *h1[H1DMAX];
//int Printf(char *format, ...);
 
146,7 → 146,7
H1D_Clear(h1d);
H1D_Print(h1d);
//Printf("InitH1D 0x%x\n", h1d );
return 0;
return h1d;
 
}
 
/cvi/instr/HISTO/H1D.h
13,6 → 13,7
 
 
/* 2d histogramming */
#define H1DMAX 500
#define H1D_ID 0x11
#pragma pack(4)
typedef struct {
/cvi/instr/HISTO/H2D.c
8,7 → 8,7
 
#include "H2D.h"
 
#define H2DMAX 500
 
H2D *h2[H2DMAX];
//int Printf(char *format, ...);
 
128,7 → 128,7
H2D_Clear(h2d);
H2D_Print(h2d);
//Printf("InitH2D 0x%x\n", h2d );
return 0;
return h2d;
}
 
double _VI_FUNC H2D_GetXBinCenter(int h2d,int xbin) {
/cvi/instr/HISTO/H2D.h
36,6 → 36,7
 
/* 2d histogramming */
#define H2DMAX 500
#define H2D_ID 0x10
#pragma pack(4)
typedef struct {
59,7 → 60,7
#pragma pack()
 
double _VI_FUNC H2D_GetYBinCenter(int h2d,int ybin);
double _VI_FUNC H2D_GetYBinCenter(int h2d,int xbin);
double _VI_FUNC H2D_GetXBinCenter(int h2d,int xbin);
 
int _VI_FUNC H2D_Clear(int h2d);
int _VI_FUNC H2D_Print(int h2d);
/cvi/instr/HISTO/H3D.c
10,7 → 10,7
#include "H2D.h"
#include "H1D.h"
 
#define H3DMAX 500
 
H3D *h3[H3DMAX];
//int Printf(char *format, ...);
 
169,6 → 169,8
}
 
int hid=499;
int i=0;
int j=0;
switch (direction){
case 0:// xy
H2D_Init(hid,"projectionXY","projectionXY",
175,8 → 177,8
H3D_GetNbinsX(histogram), H3D_GetMinX(histogram),H3D_GetMaxX(histogram),
H3D_GetNbinsY(histogram), H3D_GetMinY(histogram),H3D_GetMaxY(histogram)
);
for (int i=0; i < H3D_GetNbinsX(histogram); i++ )
for (int j=0; j < H3D_GetNbinsY(histogram); j++ )
for (i=0; i < H3D_GetNbinsX(histogram); i++ )
for (j=0; j < H3D_GetNbinsY(histogram); j++ )
H2D_SetBinContent(hid,i,j,H3D_GetBinContent(histogram,i,j,slice));
break;
case 1:// xz
184,8 → 186,8
H3D_GetNbinsX(histogram), H3D_GetMinX(histogram),H3D_GetMaxX(histogram),
H3D_GetNbinsZ(histogram), H3D_GetMinZ(histogram),H3D_GetMaxZ(histogram)
);
for (int i=0; i < H3D_GetNbinsX(histogram); i++ )
for (int j=0; j < H3D_GetNbinsZ(histogram); j++ )
for (i=0; i < H3D_GetNbinsX(histogram); i++ )
for (j=0; j < H3D_GetNbinsZ(histogram); j++ )
H2D_SetBinContent(hid,i,j,H3D_GetBinContent(histogram,i,slice,j));
break;
case 2:// yz
194,8 → 196,8
H3D_GetNbinsY(histogram), H3D_GetMinY(histogram),H3D_GetMaxY(histogram),
H3D_GetNbinsZ(histogram), H3D_GetMinZ(histogram),H3D_GetMaxZ(histogram)
);
for (int i=0; i < H3D_GetNbinsY(histogram); i++ )
for (int j=0; j < H3D_GetNbinsZ(histogram); j++ )
for (i=0; i < H3D_GetNbinsY(histogram); i++ )
for (j=0; j < H3D_GetNbinsZ(histogram); j++ )
H2D_SetBinContent(hid,i,j,H3D_GetBinContent(histogram,slice,i,j));
break;
/cvi/instr/HISTO/H3D.h
23,6 → 23,7
#include <stdio.h>
 
/* 2d histogramming */
#define H3DMAX 500
#define H3D_ID 0x12
#pragma pack(4)
typedef struct
54,8 → 55,9
 
int _VI_FUNC H3D_SliXY(int histogram,int slice, int direction);
 
double _VI_FUNC H3D_GetZBinCenter(int h3d,int zbin);
double _VI_FUNC H3D_GetYBinCenter(int h3d,int ybin);
double _VI_FUNC H3D_GetYBinCenter(int h3d,int xbin);
double _VI_FUNC H3D_GetXBinCenter(int h3d,int xbin);
 
int _VI_FUNC H3D_Clear(int h3d);
int _VI_FUNC H3D_Print(int h3d);
/cvi/instr/WUSBVME_DLL/wusbvme_dll.c
13,6 → 13,8
#endif
 
#include "wusbvme_dll.h"
#include <userint.h>
#include <utility.h>
 
usb_dev_handle *WUSB_udev;
 
69,16 → 71,21
 
 
void _VI_FUNC WIENER_VMUSB_load (const char* module_path) {
if (module_path == NULL) {
DLLHandle = LoadLibrary("c:\\Windows\\system32\\libxxusb.dll");
} else {
DLLHandle = LoadLibrary(module_path);
printf("Module Path %s\n", module_path);
}
if (!DLLHandle) {
printf ("\n\nFailed to Open libxxusb.dll \n");
printf ("\n\nFailed to Open libxxusb.dll \n");
printf("LastError 0x%x\n",GetLastError() );
MessagePopup ("ERROR", "Failed to load libxxusb.dll ");
return;
}
} else {
printf ("\n\nSuccesfully Opened libxxusb.dll \n");
}
if (!(xxusb_register_read_Ptr = (void *) GetProcAddress(DLLHandle,"xxusb_register_read"))) exit(1);
if (!(xxusb_stack_read_Ptr = (void*) GetProcAddress(DLLHandle,"xxusb_stack_read"))) exit(1);
if (!(xxusb_stack_write_Ptr = (void*) GetProcAddress(DLLHandle,"xxusb_stack_write"))) exit(1);
120,9 → 127,9
 
xxusb_device_type devices[100];
struct usb_device *dev;
if (serial != NULL)
if (serial != NULL) {
WUSB_udev = xxusb_serial_open((char *) serial);
else {
} else {
//Find XX_USB devices and open the first one found
xxusb_devices_find(devices);
dev = devices[0].usbdev;
157,7 → 164,7
WIENER_VMUSB_StackInit();
}
 
int WIENER_VMUSB_StackGetUInt16(int maxn, uint16_t *stack) {
int WIENER_VMUSB_StackGetUint16(int maxn, uint16_t *stack) {
int i;
uint16_t * u16Stack = (uint16_t *) fStack;
for (i=0; i< fStackSize*2; i++) {
167,7 → 174,7
return fStackSize*2;
}
 
int WIENER_VMUSB_StackGetUInt32(int maxn, uint32_t *stack) {
int WIENER_VMUSB_StackGetUint32(int maxn, uint32_t *stack) {
int i;
uint16_t * u16Stack = (uint16_t *) fStack;
for (i=0; i< fStackSize*2; i++) {
235,7 → 242,7
}
uint16_t * exstack = (uint16_t *) wusbcc_stack;
exstack[0]=0;
nb = WIENER_VMUSB_StackGetUInt16(MAXSTACKSIZE, &exstack[1] );
nb = WIENER_VMUSB_StackGetUint16(MAXSTACKSIZE, &exstack[1] );
//printf("WUSBXX_stack_execute stack size %d\n", nb);
//WIENER_VMUSB_StackPrint();
nb = xxusb_longstack_execute (dev, exstack, nb, 1000);
248,6 → 255,8
 
 
 
 
 
#ifndef VME_D32
 
#define VME_D8 0x1
/cvi/instr/WUSBVME_DLL/wusbvme_dll.h
38,10 → 38,10
void WIENER_VMUSB_StackClear();
 
 
#define WIENER_VMUSB_StackReadA24D16( ADDR, DATA ) WIENER_VMUSB_StackAppend(CMD_READ | CMD_A24 , CMD_D16 | (ADDR) )
#define WIENER_VMUSB_StackReadA24D32( ADDR, DATA ) WIENER_VMUSB_StackAppend(CMD_READ | CMD_A24 , CMD_D32 | (ADDR) )
#define WIENER_VMUSB_StackReadA32D16( ADDR, DATA ) WIENER_VMUSB_StackAppend(CMD_READ | CMD_A32 , CMD_D16 | (ADDR) )
#define WIENER_VMUSB_StackReadA32D32( ADDR, DATA ) WIENER_VMUSB_StackAppend(CMD_READ | CMD_A32 , CMD_D32 | (ADDR) )
#define WIENER_VMUSB_StackReadA24D16( ADDR ) WIENER_VMUSB_StackAppend(CMD_READ | CMD_A24 , CMD_D16 | (ADDR) )
#define WIENER_VMUSB_StackReadA24D32( ADDR ) WIENER_VMUSB_StackAppend(CMD_READ | CMD_A24 , CMD_D32 | (ADDR) )
#define WIENER_VMUSB_StackReadA32D16( ADDR ) WIENER_VMUSB_StackAppend(CMD_READ | CMD_A32 , CMD_D16 | (ADDR) )
#define WIENER_VMUSB_StackReadA32D32( ADDR ) WIENER_VMUSB_StackAppend(CMD_READ | CMD_A32 , CMD_D32 | (ADDR) )
 
#define WIENER_VMUSB_StackWriteA24D16( ADDR, DATA) WIENER_VMUSB_StackAppendData(CMD_WRITE | CMD_A24,CMD_D16 | (ADDR), (DATA))
#define WIENER_VMUSB_StackWriteA24D32( ADDR, DATA) WIENER_VMUSB_StackAppendData(CMD_WRITE | CMD_A24,CMD_D32 | (ADDR), (DATA))
65,6 → 65,7
void _VI_FUNC WIENER_VMUSB_load (const char *module_path);
void _VI_FUNC WIENER_VMUSB_open (const char *serial);
void _VI_FUNC WIENER_VMUSB_close (void);
extern usb_dev_handle *WUSB_udev;
 
#define WIENER_VMUSB_VME_START(NODE) WIENER_VMUSB_load(NULL);WIENER_VMUSB_open((NODE));
#define WIENER_VMUSB_VME_STOP() WIENER_VMUSB_close()