proc init() { #define STRAHLZEIT 16 enable("ClusterMatch"); enable("VertexFit"); enable("KinFit"); //enable("FuzzyFilter"); broker("EDDAruninfo:auto", 2); broker("EDDAruninfo:strahlzeit", STRAHLZEIT); select ALPHA; access BEAM_INFO; access FIRST_MATCH; access KINFIT; access VERTEX_FT; access CLPATTERN; access ISPATTERN; access USER; //////////////////////////////////////// // // do not edit the histo-definitions! // alph=histo2d("ALPHA0", 100, 0, 30, 18, 400, 2800); dim firs[24]; firs[0]=histo2d("FIRST_MATCH0", 100, 0, 90, 18, 400, 2800); firs[1]=histo2d("FIRST_MATCH1", 100, 0, 90, 18, 400, 2800); firs[2]=histo2d("FIRST_MATCH2", 100, 0, 90, 18, 400, 2800); firs[3]=histo2d("FIRST_MATCH3", 100, 0, 180, 18, 400, 2800); firs[4]=histo2d("FIRST_MATCH4", 100, 0, 180, 18, 400, 2800); firs[5]=histo2d("FIRST_MATCH5", 100, 0, 1000, 18, 400, 2800); firs[6]=histo2d("FIRST_MATCH6", 100, 0, 1000, 18, 400, 2800); firs[7]=histo2d("FIRST_MATCH7", 100, 0, 360, 18, 400, 2800); firs[8]=histo2d("FIRST_MATCH8", 100, 0, 360, 18, 400, 2800); firs[9]=histo2d("FIRST_MATCH9", 100, 0, 10, 18, 400, 2800); firs[10]=histo2d("FIRST_MATCH10", 100, 0, 500, 18, 400, 2800); firs[11]=histo2d("FIRST_MATCH11", 100, 0, 500, 18, 400, 2800); firs[12]=histo2d("FIRST_MATCH12", 100, 0, 500, 18, 400, 2800); firs[13]=histo2d("FIRST_MATCH13", 100, 0, 500, 18, 400, 2800); firs[14]=histo2d("FIRST_MATCH14", 100, 0, 500, 18, 400, 2800); firs[15]=histo2d("FIRST_MATCH15", 100, 0, 500, 18, 400, 2800); firs[16]=histo2d("FIRST_MATCH16", 100, 0, 360, 18, 400, 2800); firs[17]=histo2d("FIRST_MATCH17", 100, 0, 360, 18, 400, 2800); firs[18]=histo2d("FIRST_MATCH18", 100, 0, 500, 18, 400, 2800); firs[19]=histo2d("FIRST_MATCH19", 100, -50, 50, 18, 400, 2800); firs[20]=histo2d("FIRST_MATCH20", 100, 0, 1000, 18, 400, 2800); firs[21]=histo2d("FIRST_MATCH21", 100, 0, 1000, 18, 400, 2800); firs[22]=histo2d("FIRST_MATCH22", 100, 0, 500, 18, 400, 2800); dim kinf[6]; kinf[0]=histo2d("KINFIT0", 100, 0, 90, 18, 400, 2800); kinf[1]=histo2d("KINFIT1", 100, 0, 360, 18, 400, 2800); kinf[2]=histo2d("KINFIT2", 100, -30, 30, 18, 400, 2800); kinf[3]=histo2d("KINFIT3", 100, -30, 30, 18, 400, 2800); kinf[4]=histo2d("KINFIT4", 100, 0, 20, 18, 400, 2800); kinf[5]=histo2d("KINFIT5", 100, 0, 20, 18, 400, 2800); dim vert[13]; vert[0]=histo2d("VERTEX_FT0", 100, -30, 30, 18, 400, 2800); vert[1]=histo2d("VERTEX_FT1", 100, -30, 30, 18, 400, 2800); vert[2]=histo2d("VERTEX_FT2", 100, -50, 200, 18, 400, 2800); vert[3]=histo2d("VERTEX_FT3", 100, 0, 25, 18, 400, 2800); vert[5]=histo2d("VERTEX_FT5", 100, 0, 20, 18, 400, 2800); vert[6]=histo2d("VERTEX_FT6", 100, 0, 180, 18, 400, 2800); vert[7]=histo2d("VERTEX_FT7", 100, 0, 180, 18, 400, 2800); vert[8]=histo2d("VERTEX_FT8", 100, 0, 180, 18, 400, 2800); vert[9]=histo2d("VERTEX_FT9", 100, 0, 180, 18, 400, 2800); vert[10]=histo2d("VERTEX_FT10", 100, 0, 180, 18, 400, 2800); vert[11]=histo2d("VERTEX_FT11", 100, 0, 180, 18, 400, 2800); vert[12]=histo2d("VERTEX_FT12", 100, 0, 30, 18, 400, 2800); dim user[6]; user[0]=histo2d("USER0", 100, 0, 100, 18, 400, 2800); user[1]=histo2d("USER1", 100, 0, 100, 18, 400, 2800); user[2]=histo2d("USER2", 100, 0, 100, 18, 400, 2800); user[3]=histo2d("USER3", 100, 0, 100, 18, 400, 2800); user[4]=histo2d("USER4", 100, 0, 100, 18, 400, 2800); //////////////////////////////////////// } proc analyze() { hfill2d(alph,ALPHA[0], BEAM_INFO[2]); hfill2d(firs[0],FIRST_MATCH[0], BEAM_INFO[2]); hfill2d(firs[1],FIRST_MATCH[1], BEAM_INFO[2]); hfill2d(firs[2],FIRST_MATCH[2], BEAM_INFO[2]); hfill2d(firs[3],FIRST_MATCH[3], BEAM_INFO[2]); hfill2d(firs[4],FIRST_MATCH[4], BEAM_INFO[2]); hfill2d(firs[5],FIRST_MATCH[5], BEAM_INFO[2]); hfill2d(firs[6],FIRST_MATCH[6], BEAM_INFO[2]); hfill2d(firs[7],FIRST_MATCH[7], BEAM_INFO[2]); hfill2d(firs[8],FIRST_MATCH[8], BEAM_INFO[2]); hfill2d(firs[9],FIRST_MATCH[9], BEAM_INFO[2]); hfill2d(firs[10],FIRST_MATCH[10], BEAM_INFO[2]); hfill2d(firs[11],FIRST_MATCH[11], BEAM_INFO[2]); hfill2d(firs[12],FIRST_MATCH[12], BEAM_INFO[2]); hfill2d(firs[13],FIRST_MATCH[13], BEAM_INFO[2]); hfill2d(firs[14],FIRST_MATCH[14], BEAM_INFO[2]); hfill2d(firs[15],FIRST_MATCH[15], BEAM_INFO[2]); hfill2d(firs[16],FIRST_MATCH[16], BEAM_INFO[2]); hfill2d(firs[17],FIRST_MATCH[17], BEAM_INFO[2]); hfill2d(firs[18],FIRST_MATCH[18], BEAM_INFO[2]); hfill2d(firs[19],FIRST_MATCH[19], BEAM_INFO[2]); hfill2d(firs[20],FIRST_MATCH[20], BEAM_INFO[2]); hfill2d(firs[21],FIRST_MATCH[21], BEAM_INFO[2]); hfill2d(firs[22],FIRST_MATCH[22], BEAM_INFO[2]); hfill2d(kinf[0],KINFIT[0], BEAM_INFO[2]); hfill2d(kinf[1],KINFIT[1], BEAM_INFO[2]); hfill2d(kinf[2],KINFIT[2], BEAM_INFO[2]); hfill2d(kinf[3],KINFIT[3], BEAM_INFO[2]); hfill2d(kinf[4],KINFIT[4], BEAM_INFO[2]); hfill2d(kinf[5],KINFIT[5], BEAM_INFO[2]); hfill2d(vert[0],VERTEX_FT[0], BEAM_INFO[2]); hfill2d(vert[1],VERTEX_FT[1], BEAM_INFO[2]); hfill2d(vert[2],VERTEX_FT[2], BEAM_INFO[2]); hfill2d(vert[3],VERTEX_FT[3], BEAM_INFO[2]); hfill2d(vert[5],VERTEX_FT[5], BEAM_INFO[2]); hfill2d(vert[6],VERTEX_FT[6], BEAM_INFO[2]); hfill2d(vert[7],VERTEX_FT[7], BEAM_INFO[2]); hfill2d(vert[8],VERTEX_FT[8], BEAM_INFO[2]); hfill2d(vert[9],VERTEX_FT[9], BEAM_INFO[2]); hfill2d(vert[10],VERTEX_FT[10], BEAM_INFO[2]); hfill2d(vert[11],VERTEX_FT[11], BEAM_INFO[2]); hfill2d(vert[12],VERTEX_FT[12], BEAM_INFO[2]); for (i=0; i