This File Converts MYOWN Data Files into a Canned Format for Later Use
     It Creates Files lengthinfo.tst,objsinfo.tst,eveinfo.tst

<<Graphics`Graphics`
nameoffile="~/pheno/hidval/hv1006/hv306001/HV.had.dat"
blacklist =Import[nameoffile,"Table"];
If[Part[Part[blacklist,1],1]==1,shift=0;Nrec=Length[blacklist],shift=1;Nrec=Length[blacklist]-1]
Length[blacklist]
shift
Nrec
DeleteCases[Table[{i,Length[Part[blacklist,i]]},{i,Length[blacklist]/1}],{_,15}]

theta[eta_]:=2 ArcTan[Exp[-eta]]+.00001
px[eta_,phi_,pt_]:=Round[pt Cos[phi]*100]*.01
py[eta_,phi_,pt_]:=Round[pt Sin[phi]*100]*.01
pz[eta_,phi_,pt_]:=Round[pt/Tan[theta[eta]]*100]*.01
ee[eta_,phi_,pt_,jmass_]:=Round[Sqrt[pt^2(1+1/Tan[theta[eta]]^2)+jmass^2]*100]*.01
newname[a_,b_,c_] :=If[a-1==0,IntegerPart[1*b]*a-a,0]+If[a-2==0,IntegerPart[1*b]*a-a,0]+If[a-3==0,IntegerPart[1*b]*a-a,0]+If[a-4==0,IntegerPart[1*(1-2c)]*a-a,0]+a

fvecmass[r1_]:=Re[Sqrt[Part[r1,9]^2-Part[r1,6]^2-Part[r1,7]^2-Part[r1,8]^2]]
fvecmassboost[r1_]:=Re[Sqrt[Part[r1,13]^2-Part[r1,6]^2-Part[r1,7]^2-Part[r1,12]^2]]
pzzboost[pz_,ee_,beta_]:=1/Sqrt[1-beta^2](pz+ee beta)
eezboost[pz_,ee_,beta_]:=1/Sqrt[1-beta^2](ee+pz beta)
zboost4vec[px_,py_,pz_,ee_,beta_]:={px,py,1/Sqrt[1-beta^2](pz+ee beta),1/Sqrt[1-beta^2](ee+pz beta)}
boostrec[r1_,zb_]:=r1-{0,0,0,0,0,Part[r1,6],Part[r1,7],Part[r1,8],Part[r1,9],0,0}+{0,0,0,0,0,Part[zb,1],Part[zb,2],Part[zb,3],Part[zb,4],0,0}
getbetaz[r1_]:=If[Part[r1,9]>0,Part[r1,8]/Part[r1,9],0]
rapidity[r1_]:=If[Part[r1,9]>0,ArcTanh[Part[r1,8]/Part[r1,9]],0]
psrapidity[r1_]:=- Log[Tan[ArcTan[Sqrt[Part[r1,6]^2+Part[r1,7]^2]/(.0001+pz[Part[r1,3],Part[r1,4],Part[r1,5]])]/2]]

                                                                                                  
convert[record_] := {Part[record,1],
newname[Part[record,2],Part[record,6],Part[record,8]],
Part[record,3],Part[record,4],Part[record,5],
px[Part[record,3],Part[record,4],Part[record,5]],
py[Part[record,3],Part[record,4],Part[record,5]],
pz[Part[record,3],Part[record,4],Part[record,5]],
ee[Part[record,3],Part[record,4],Part[record,5],If[Part[record,2]==4,Part[record,6],0]],
If[Part[record,2]==4,Part[record,6],0],Part[record,7]}


Do[Do[whitelist[i]=Part[blacklist,i+shift],{i,1,Nrec}]]
n=1;Do[If[Part[whitelist[i],1]==1,pointer[n]=i;n=n+1],{i,1,Nrec}]
pointer[n]=Nrec+1
Nev=n-1

Nev=25
Nrec=pointer[26]-1

Clear[Tmvis,Tbetavis,Tptot0,Tmvis0,Timbalance,Tetavis,Tptot,Tht,Tmvis2,TMET]
Tbetavis=Table[getbetaz[Sum[whitelist[i]UnitStep[Part[whitelist[i],5]-25],{i,pointer[n],pointer[n+1]-2}]],{n,1,Nev}];

Do[If[Round[n/10000]==n/10000,Print[n]];Do[bluelist[i]=If[Part[whitelist[i],2]!=6,Flatten[Append[whitelist[i],{pzzboost[Part[whitelist[i],8],Part[whitelist[i],9],-Part[Tbetavis,n]],eezboost[Part[whitelist[i],8],Part[whitelist[i],9],-Part[Tbetavis,n]]}]],Flatten[Append[whitelist[i],{Part[whitelist[i],8],Part[whitelist[i],9]}]]],{i,pointer[n],pointer[n+1]-1}],{n,1,Nev}]

Tptot0=Table[Sum[whitelist[i],{i,pointer[n],pointer[n+1]-1}],{n,1,Nev}];
Tmvis0=Table[Re[fvecmass[Part[Tptot0,n]]],{n,1,Nev}];

Tptot=Table[Sum[bluelist[i]UnitStep[Part[bluelist[i],5]-25],{i,pointer[n],pointer[n+1]-1}],{n,1,Nev}];
Tmvis=Table[Re[fvecmassboost[Part[Tptot,n]]],{n,1,Nev}];

Tetavis=Table[rapidity[Part[Tptot,n]],{n,1,Nev}];
TMET=Table[If[Part[whitelist[pointer[n+1]-1],2]==6,Part[whitelist[pointer[n+1]-1],5],0],{n,1,Nev}];
Timbalance=Table[Sqrt[(Part[Part[Tptot,n],6])^2+(Part[Part[Tptot,n],7])^2],{n,1,Nev}];

Do[If[Round[n/10000]==n/10000,Print[n]];eventinfo[n]=Flatten[{n,pointer[n],Part[Timbalance,n],Part[Tmvis0,n],Part[Tbetavis,n],Part[Tmvis,n],If[Part[whitelist[pointer[n+1]-1],2]==6,Part[whitelist[pointer[n+1]-1],5],0],Table[Count[Table[If[Part[whitelist[j],5]>3,whitelist[j],{9,9,9,9,9,9,9,9,9,9,9,9,9,9,9}],{j,pointer[n],pointer[n+1]-1}],{_,k,_,_,_,_,_,_,_,_,_,_,_,_,_}],{k,-4,6}],Table[Count[Table[If[Part[whitelist[j],5]>-1,whitelist[j],{9,9,9,9,9,9,9,9,9,9,9,9,9,9,9}],{j,pointer[n],pointer[n+1]-1}],{_,k,_,_,_,_,_,_,_,_,_,_,_,_,_}],{k,-4,6}]}],{n,1,Nev}]
OpenWrite["lengthinfo.tst"]
Write["lengthinfo.tst",Nrec]
Write["lengthinfo.tst",Nev]
Write["lengthinfo.tst",nameoffile]
Close["lengthinfo.tst"]

OpenWrite["objinfo.tst"]
Do[Write["objinfo.tst",bluelist[i]],{i,1,Nrec}]
Close["objinfo.tst"]
OpenWrite["eveinfo.tst"]
Do[Write["eveinfo.tst",eventinfo[n]],{n,1,Nev}]
Close["eveinfo.tst"]

~/pheno/hidval/hv1006/hv306001/HV.had.dat

385669

385669

0

385669

{}

385670

1000

25

10725

General :: spell1 : Possible spelling error: new symbol name \"Tetavis\" is similar to existing symbol \"Tbetavis\".  More…

OutputStream[lengthinfo.tst, 9]

lengthinfo.tst

OutputStream[objinfo.tst, 10]

objinfo.tst

OutputStream[eveinfo.tst, 11]

eveinfo.tst

"

<<"startup.m"

<<Graphics`MultipleListPlot`
<<Graphics`Graphics3D`

Nev
Nev=25

red[ii_]:=(KroneckerDelta[id[ii],4]-KroneckerDelta[id[ii],1])Min[pt[ii]/6+.5,1]+KroneckerDelta[id[ii],1]
green[ii_]:=(KroneckerDelta[id[ii],2]-KroneckerDelta[id[ii],1])Min[pt[ii]/6+.5,1]+KroneckerDelta[id[ii],1]
blue[ii_]:=(KroneckerDelta[id[ii],0]-KroneckerDelta[id[ii],1])Min[pt[ii]/6+.5,1]+KroneckerDelta[id[ii],1]
cal[ii_]:=If[id[ii]==6,0,If[vertx[ii]^2+verty[ii]^2<9,Max[cc /. Solve[(vertx[ii]/1000.+cc*(px[ii]+.001))^2+(verty[ii]/1000.+cc*(py[ii]+.001))^2==9,cc],0],0]]
startx[ii_]:=If[id[ii]==0||id[ii]==4,vertx[ii]/1000.+cal[ii] px[ii],vertx[ii]/1000.]
starty[ii_]:=If[id[ii]==0||id[ii]==4,verty[ii]/1000.+cal[ii] py[ii],verty[ii]/1000.]
endx[ii_]:=If[id[ii]==0||id[ii]==4,startx[i]+px[i]/2,startx[ii]+UnitStep[pt[ii]-3]px[i]/2]
endy[ii_]:=If[id[ii]==0||id[ii]==4,starty[i]+py[i]/2,starty[ii]+UnitStep[pt[ii]-3]py[i]/2]

vertx[ii_]:=Part[bluelist[ii],13]
verty[ii_]:=Part[bluelist[ii],14]
vertz[ii_]:=Part[bluelist[ii],15]

A0[0]=0;A0[1]=0;A0[2]=0;A0[4]=0;A0[6]=0;B0[0]=1;B0[1]=1;B0[2]=1;B0[4]=1;B0[6]=1;CC0[0]={RGBColor[.7,0,1],Thickness[0.02]};C0[1]={RGBColor[.7,.2,1],Thickness[0.02]};C0[2]={RGBColor[.7,.4,1],Thickness[0.02]};C0[4]={RGBColor[1,0,0],Thickness[0.02]};C0[6]={RGBColor[0,0,0],Thickness[0.02]};


Pxy[i_]:=Sqrt[(px[i])^2+(py[i])^2+.0000001]
Pxz[i_]:=Sqrt[(px[i])^2+(pz[i])^2+.0000001]
Pyz[i_]:=Sqrt[(pz[i])^2+(py[i])^2+.0000001]


TrkXY[i_]:={{A0[id[i]]*px[i]/Pxy[i],A0[id[i]]*py[i]/Pxy[i]},{px[i](B0[id[i]]+A0[id[i]]/Pxy[i]),py[i](B0[id[i]]+A0[id[i]]/Pxy[i])}}
TrkXZ[i_]:={{A0[id[i]]*px[i]/Pxy[i],A0[id[i]]*pz[i]/Pxy[i]},{px[i](B0[id[i]]+A0[id[i]]/Pxy[i]),pz[i](B0[id[i]]+A0[id[i]]/Pxy[i])}}
TrkYZ[i_]:={{A0[id[i]]*py[i]/Pxy[i],A0[id[i]]*pz[i]/Pxy[i]},{py[i](B0[id[i]]+A0[id[i]]/Pxy[i]),pz[i](B0[id[i]]+A0[id[i]]/Pxy[i])}}
Clear[lego,Tlego]
Do[
lego=SparseArray[{j_,k_}->0,{31,31}];
Tlego=SparseArray[{j_,k_}->0,{31,31}];
Do[If[Abs[eta[i]]<3,lego[[IntegerPart[(eta[i]+3)*5+1],IntegerPart[1+phi[i]/2/Pi*30]]]+=pt[i]],{i,pointer[nn],pointer[nn+1]-2}];
Do[If[lego[[j,k]]<1,lego[[j,k]]=0],{j,1,31},{k,1,31}];
Do[Tlego[[j,k]]=lego[[k,j]],{j,1,31},{k,1,31}];
LEGO1[nn]=Tlego;
LEGO2[nn]=lego
    ,{nn,1,Nev}]
    
Do[
HadPlots[nn]=Table[ {AbsoluteThickness[Min[Pxy[i]/4,2]],Hue[.3+.3*pz[i]/(ee[i]+.001)], Dashing[{0.01-0.008*(KroneckerDelta[id[i],2]),0.02*(.5 KroneckerDelta[id[i],2]+ KroneckerDelta[id[i],0])}], Line[{{startx[i], starty[i]}, {endx[i],endy[i]}}]},{i,pointer[nn], pointer[nn+1]-2}],{nn,1,Nev}]

25

25

General :: spell1 : Possible spelling error: new symbol name \"red\" is similar to existing symbol \"Red\".  More…

General :: spell1 : Possible spelling error: new symbol name \"green\" is similar to existing symbol \"Green\".  More…

General :: spell1 : Possible spelling error: new symbol name \"blue\" is similar to existing symbol \"Blue\".  More…

General :: spell1 : Possible spelling error: new symbol name \"verty\" is similar to existing symbol \"vertx\".  More…

General :: spell1 : Possible spelling error: new symbol name \"starty\" is similar to existing symbol \"startx\".  More…

General :: spell1 : Possible spelling error: new symbol name \"endy\" is similar to existing symbol \"endx\".  More…

General :: spell : Possible spelling error: new symbol name \"vertz\" is similar to existing symbols  {vertx, verty} .  More…

General :: spell1 : Possible spelling error: new symbol name \"TrkXZ\" is similar to existing symbol \"TrkXY\".  More…

General :: spell1 : Possible spelling error: new symbol name \"TrkYZ\" is similar to existing symbol \"TrkXZ\".  More…

General :: spell1 : Possible spelling error: new symbol name \"Tlego\" is similar to existing symbol \"lego\".  More…


nameoffile="~/pheno/hidval/atlas919/HV.bb.dat"
blacklist =Import[nameoffile,"Table"];
If[Part[Part[blacklist,1],1]==1,shift=0;Nrec=Length[blacklist],shift=1;Nrec=Length[blacklist]-1]
Length[blacklist]
shift
Nrec
DeleteCases[Table[{i,Length[Part[blacklist,i]]},{i,Length[blacklist]/1}],{_,15}]

theta[eta_]:=2 ArcTan[Exp[-eta]]+.00001
px[eta_,phi_,pt_]:=Round[pt Cos[phi]*100]*.01
py[eta_,phi_,pt_]:=Round[pt Sin[phi]*100]*.01
pz[eta_,phi_,pt_]:=Round[pt/Tan[theta[eta]]*100]*.01
ee[eta_,phi_,pt_,jmass_]:=Round[Sqrt[pt^2(1+1/Tan[theta[eta]]^2)+jmass^2]*100]*.01
newname[a_,b_,c_] :=If[a-1==0,IntegerPart[1*b]*a-a,0]+If[a-2==0,IntegerPart[1*b]*a-a,0]+If[a-3==0,IntegerPart[1*b]*a-a,0]+If[a-4==0,IntegerPart[1*(1-2c)]*a-a,0]+a

fvecmass[r1_]:=Re[Sqrt[Part[r1,9]^2-Part[r1,6]^2-Part[r1,7]^2-Part[r1,8]^2]]
fvecmassboost[r1_]:=Re[Sqrt[Part[r1,13]^2-Part[r1,6]^2-Part[r1,7]^2-Part[r1,12]^2]]
pzzboost[pz_,ee_,beta_]:=1/Sqrt[1-beta^2](pz+ee beta)
eezboost[pz_,ee_,beta_]:=1/Sqrt[1-beta^2](ee+pz beta)
zboost4vec[px_,py_,pz_,ee_,beta_]:={px,py,1/Sqrt[1-beta^2](pz+ee beta),1/Sqrt[1-beta^2](ee+pz beta)}
boostrec[r1_,zb_]:=r1-{0,0,0,0,0,Part[r1,6],Part[r1,7],Part[r1,8],Part[r1,9],0,0}+{0,0,0,0,0,Part[zb,1],Part[zb,2],Part[zb,3],Part[zb,4],0,0}
getbetaz[r1_]:=If[Part[r1,9]>0,Part[r1,8]/Part[r1,9],0]
rapidity[r1_]:=If[Part[r1,9]>0,ArcTanh[Part[r1,8]/Part[r1,9]],0]
psrapidity[r1_]:=- Log[Tan[ArcTan[Sqrt[Part[r1,6]^2+Part[r1,7]^2]/(.0001+pz[Part[r1,3],Part[r1,4],Part[r1,5]])]/2]]

                                                                                                  
convert[record_] := {Part[record,1],
newname[Part[record,2],Part[record,6],Part[record,8]],
Part[record,3],Part[record,4],Part[record,5],
px[Part[record,3],Part[record,4],Part[record,5]],
py[Part[record,3],Part[record,4],Part[record,5]],
pz[Part[record,3],Part[record,4],Part[record,5]],
ee[Part[record,3],Part[record,4],Part[record,5],If[Part[record,2]==4,Part[record,6],0]],
If[Part[record,2]==4,Part[record,6],0],Part[record,7]}


Do[Do[whitelist[i]=Part[blacklist,i+shift],{i,1,Nrec}]]
n=1;Do[If[Part[whitelist[i],1]==1,pointer[n]=i;n=n+1],{i,1,Nrec}]
pointer[n]=Nrec+1
Nev=n-1


Clear[Tmvis,Tbetavis,Tptot0,Tmvis0,Timbalance,Tetavis,Tptot,Tht,Tmvis2,TMET]
Tbetavis=Table[getbetaz[Sum[whitelist[i]UnitStep[Part[whitelist[i],5]-25],{i,pointer[n],pointer[n+1]-2}]],{n,1,Nev}];

Do[If[Round[n/10000]==n/10000,Print[n]];Do[bluelist[i]=If[Part[whitelist[i],2]!=6,Flatten[Append[whitelist[i],{pzzboost[Part[whitelist[i],8],Part[whitelist[i],9],-Part[Tbetavis,n]],eezboost[Part[whitelist[i],8],Part[whitelist[i],9],-Part[Tbetavis,n]]}]],Flatten[Append[whitelist[i],{Part[whitelist[i],8],Part[whitelist[i],9]}]]],{i,pointer[n],pointer[n+1]-1}],{n,1,Nev}]

Tptot0=Table[Sum[whitelist[i],{i,pointer[n],pointer[n+1]-1}],{n,1,Nev}];
Tmvis0=Table[Re[fvecmass[Part[Tptot0,n]]],{n,1,Nev}];

Tptot=Table[Sum[bluelist[i]UnitStep[Part[bluelist[i],5]-25],{i,pointer[n],pointer[n+1]-1}],{n,1,Nev}];
Tmvis=Table[Re[fvecmassboost[Part[Tptot,n]]],{n,1,Nev}];

Tetavis=Table[rapidity[Part[Tptot,n]],{n,1,Nev}];
TMET=Table[If[Part[whitelist[pointer[n+1]-1],2]==6,Part[whitelist[pointer[n+1]-1],5],0],{n,1,Nev}];
Timbalance=Table[Sqrt[(Part[Part[Tptot,n],6])^2+(Part[Part[Tptot,n],7])^2],{n,1,Nev}];

Do[If[Round[n/10000]==n/10000,Print[n]];eventinfo[n]=Flatten[{n,pointer[n],Part[Timbalance,n],Part[Tmvis0,n],Part[Tbetavis,n],Part[Tmvis,n],If[Part[whitelist[pointer[n+1]-1],2]==6,Part[whitelist[pointer[n+1]-1],5],0],Table[Count[Table[If[Part[whitelist[j],5]>3,whitelist[j],{9,9,9,9,9,9,9,9,9,9,9,9,9,9,9}],{j,pointer[n],pointer[n+1]-1}],{_,k,_,_,_,_,_,_,_,_,_,_,_,_,_}],{k,-4,6}],Table[Count[Table[If[Part[whitelist[j],5]>-1,whitelist[j],{9,9,9,9,9,9,9,9,9,9,9,9,9,9,9}],{j,pointer[n],pointer[n+1]-1}],{_,k,_,_,_,_,_,_,_,_,_,_,_,_,_}],{k,-4,6}]}],{n,1,Nev}]
OpenWrite["lengthinfo.tst"]
Write["lengthinfo.tst",Nrec]
Write["lengthinfo.tst",Nev]
Write["lengthinfo.tst",nameoffile]
Close["lengthinfo.tst"]

OpenWrite["objinfo.tst"]
Do[Write["objinfo.tst",bluelist[i]],{i,1,Nrec}]
Close["objinfo.tst"]
OpenWrite["eveinfo.tst"]
Do[Write["eveinfo.tst",eventinfo[n]],{n,1,Nev}]
Close["eveinfo.tst"]


Do[
BPlots[nn]=Table[ {AbsoluteThickness[Min[Pxy[i]/30,4]],Hue[Part[bluelist[i],10]/15.],
Dashing[{0.01,If[Abs[id[i]-100]==15,.05,0]}],Line[{{startx[i], starty[i]}, {endx[i],endy[i]}}]},{i,pointer[nn], pointer[nn+1]-2}],{nn,1,Nev}]


Do[
Show[
GraphicsArray[
{{Graphics[HadPlots[nn],PlotRange->{{-5,5},{-5,5}},AspectRatio->1,Frame->True,PlotLabel->nn],Graphics[BPlots[nn],PlotRange->{{-5,5},{-5,5}},AspectRatio->1,Frame->True,PlotLabel->nn]},{
Histogram3D[LEGO1[nn],ViewPoint -> {0,-2.6,3},FrequencyData->True,DisplayFunction→Identity,DefaultColor->RGBColor[0,0,0],TextStyle→{FontFamily→"Times",FontWeight->"Bold",FontSlant->"Italic",FontSize→12},Frame->True],DensityGraphics[LEGO2[nn]/Max[LEGO2[nn]],TextStyle→{FontFamily→"Times",FontWeight->"Bold",
    FontSlant->"Italic",FontSize→12},PlotLabel->"LEGO plot,   eta  vs. phi  in .2 bins",ColorFunctionScaling->False,ColorFunction->(RGBColor[Min[4#,1]^.2,Min[4#,1]^.5,Min[4#,1]^.1]&),DefaultColor->RGBColor[0,0,0],PlotRegion->Automatic]}},ImageSize->{600,600}]],{nn,1,Nev}]


-----------------------------------------------------------------------------------------------------------------------------------------------------------------

Do[
Show[
GraphicsArray[
{Graphics[HadPlots[nn],PlotRange->{{-5,5},{-5,5}},AspectRatio->1,Frame->True,PlotLabel->nn],Graphics[HadPlots[nn],PlotRange->{{-.1,.1},{-.1,.1}},AspectRatio->1,Frame->True,PlotLabel->nn],
Histogram3D[LEGO1[nn],ViewPoint -> {0,-2.6,3},FrequencyData->True,DisplayFunction→Identity,DefaultColor->RGBColor[0,0,0],TextStyle→{FontFamily→"Times",FontWeight->"Bold",FontSlant->"Italic",FontSize→12},Frame->True],DensityGraphics[LEGO2[nn]/Max[LEGO2[nn]],TextStyle→{FontFamily→"Times",FontWeight->"Bold",
    FontSlant->"Italic",FontSize→12},PlotLabel->"LEGO plot,   eta  vs. phi  in .2 bins",ColorFunctionScaling->False,ColorFunction->(RGBColor[Min[4#,1]^.2,Min[4#,1]^.5,Min[4#,1]^.1]&),DefaultColor->RGBColor[0,0,0],PlotRegion->Automatic]},ImageSize->{1200,300}]],{nn,1,25}]

[Graphics:HTMLFiles/hv306001_30.gif]

[Graphics:HTMLFiles/hv306001_31.gif]

                                     1 Power :: infy : Infinite expression  - encountered. More…                                      0

∞ :: indet : Indeterminate expression 0 ComplexInfinity encountered. More…

[Graphics:HTMLFiles/hv306001_34.gif]

[Graphics:HTMLFiles/hv306001_35.gif]

[Graphics:HTMLFiles/hv306001_36.gif]

[Graphics:HTMLFiles/hv306001_37.gif]

[Graphics:HTMLFiles/hv306001_38.gif]

[Graphics:HTMLFiles/hv306001_39.gif]

[Graphics:HTMLFiles/hv306001_40.gif]

[Graphics:HTMLFiles/hv306001_41.gif]

[Graphics:HTMLFiles/hv306001_42.gif]

[Graphics:HTMLFiles/hv306001_43.gif]

[Graphics:HTMLFiles/hv306001_44.gif]

[Graphics:HTMLFiles/hv306001_45.gif]

[Graphics:HTMLFiles/hv306001_46.gif]

[Graphics:HTMLFiles/hv306001_47.gif]

[Graphics:HTMLFiles/hv306001_48.gif]

[Graphics:HTMLFiles/hv306001_49.gif]

[Graphics:HTMLFiles/hv306001_50.gif]

[Graphics:HTMLFiles/hv306001_51.gif]

[Graphics:HTMLFiles/hv306001_52.gif]

[Graphics:HTMLFiles/hv306001_53.gif]

[Graphics:HTMLFiles/hv306001_54.gif]

[Graphics:HTMLFiles/hv306001_55.gif]

[Graphics:HTMLFiles/hv306001_56.gif]





HadPlots[nn]=Table[ {AbsoluteThickness[Min[Pxy[i]/4,2]],Hue[.3+.3*pz[i]/(ee[i]+.001)], Dashing[{0.01-0.008*(KroneckerDelta[id[i],2]),0.02*(.5 KroneckerDelta[id[i],2]+ KroneckerDelta[id[i],0])}], Line[{{startx[i], starty[i]}, {endx[i],endy[i]}}]},{i,pointer[nn], pointer[nn+1]-2}],{nn,1,25}]


Created by Mathematica  (October 8, 2006) Valid XHTML 1.1!