{
//////////////////////////////////////////////////////////
//   This file has been automatically generated 
//     (Fri Sep 17 15:27:03 2004 by ROOT version3.05/07)
//   from TTree AncSvx/SVX Ancestors
//   found on file: ancsvx_01.root
//////////////////////////////////////////////////////////


//Reset ROOT and connect tree file
   gROOT->Reset();
   TFile *f = (TFile*)gROOT->GetListOfFiles()->FindObject("ancsvx_list.dat");
   if (!f) {
      f = new TFile("ancsvx_softlink.root");
   }
   TTree *AncSvx = (TTree*)gDirectory->Get("AncSvx");

//Declaration of leaves types
   Float_t         TRACK;
   Float_t         NFILE;
   Float_t         PTOT;
   Float_t         PTHETA;
   Float_t         PPHI;
   Float_t         R_VERTEX;
   Float_t         Z_VERTEX;
   Float_t         THET_VER;
   Float_t         PHI_VER;
   Float_t         ITPARENT;
   Float_t         IDPARENT;
   Float_t         IDPART;
   Float_t         ITORIGIN;
   Float_t         IDORIGIN;
   Float_t         IHIT;
   Float_t         LAYER;
   Float_t         THETA;
   Float_t         PHI;
   Float_t         XGLOBAL;
   Float_t         YGLOBAL;
   Float_t         ZGLOBAL;
   Float_t         PMOMX;
   Float_t         PMOMY;
   Float_t         PMOMZ;
   Float_t         IPC123;
   Float_t         PC1THETA;
   Float_t         PC1PHI;
   Float_t         PC1RDIST;
   Float_t         PC1ZDIST;
   Float_t         PC2THETA;
   Float_t         PC2PHI;
   Float_t         PC2RDIST;
   Float_t         PC2ZDIST;
   Float_t         PC3THETA;
   Float_t         PC3PHI;
   Float_t         PC3RDIST;
   Float_t         PC3ZDIST;
   Float_t         DELE;
   Float_t         XLOCIN;
   Float_t         YLOCIN;
   Float_t         ZLOCIN;
   Float_t         XLOCOUT;
   Float_t         YLOCOUT;
   Float_t         ZLOCOUT;
   Float_t         HITVOL0;
   Float_t         HITVOL1;
   Float_t         HITVOL2;
   Float_t         HITVOL3;
   Float_t         HITVOL4;
   Float_t         HITVOL5;
   Float_t         NHIT;
   Float_t         X0_EVENT;
   Float_t         Y0_EVENT;
   Float_t         Z0_EVENT;
   Float_t         B_IMPACT;
   Float_t         EVENT;

//Set branch addresses
   AncSvx->SetBranchAddress("TRACK",&TRACK);
   AncSvx->SetBranchAddress("NFILE",&NFILE);
   AncSvx->SetBranchAddress("PTOT",&PTOT);
   AncSvx->SetBranchAddress("PTHETA",&PTHETA);
   AncSvx->SetBranchAddress("PPHI",&PPHI);
   AncSvx->SetBranchAddress("R_VERTEX",&R_VERTEX);
   AncSvx->SetBranchAddress("Z_VERTEX",&Z_VERTEX);
   AncSvx->SetBranchAddress("THET_VER",&THET_VER);
   AncSvx->SetBranchAddress("PHI_VER",&PHI_VER);
   AncSvx->SetBranchAddress("ITPARENT",&ITPARENT);
   AncSvx->SetBranchAddress("IDPARENT",&IDPARENT);
   AncSvx->SetBranchAddress("IDPART",&IDPART);
   AncSvx->SetBranchAddress("ITORIGIN",&ITORIGIN);
   AncSvx->SetBranchAddress("IDORIGIN",&IDORIGIN);
   AncSvx->SetBranchAddress("IHIT",&IHIT);
   AncSvx->SetBranchAddress("LAYER",&LAYER);
   AncSvx->SetBranchAddress("THETA",&THETA);
   AncSvx->SetBranchAddress("PHI",&PHI);
   AncSvx->SetBranchAddress("XGLOBAL",&XGLOBAL);
   AncSvx->SetBranchAddress("YGLOBAL",&YGLOBAL);
   AncSvx->SetBranchAddress("ZGLOBAL",&ZGLOBAL);
   AncSvx->SetBranchAddress("PMOMX",&PMOMX);
   AncSvx->SetBranchAddress("PMOMY",&PMOMY);
   AncSvx->SetBranchAddress("PMOMZ",&PMOMZ);
   AncSvx->SetBranchAddress("IPC123",&IPC123);
   AncSvx->SetBranchAddress("PC1THETA",&PC1THETA);
   AncSvx->SetBranchAddress("PC1PHI",&PC1PHI);
   AncSvx->SetBranchAddress("PC1RDIST",&PC1RDIST);
   AncSvx->SetBranchAddress("PC1ZDIST",&PC1ZDIST);
   AncSvx->SetBranchAddress("PC2THETA",&PC2THETA);
   AncSvx->SetBranchAddress("PC2PHI",&PC2PHI);
   AncSvx->SetBranchAddress("PC2RDIST",&PC2RDIST);
   AncSvx->SetBranchAddress("PC2ZDIST",&PC2ZDIST);
   AncSvx->SetBranchAddress("PC3THETA",&PC3THETA);
   AncSvx->SetBranchAddress("PC3PHI",&PC3PHI);
   AncSvx->SetBranchAddress("PC3RDIST",&PC3RDIST);
   AncSvx->SetBranchAddress("PC3ZDIST",&PC3ZDIST);
   AncSvx->SetBranchAddress("DELE",&DELE);
   AncSvx->SetBranchAddress("XLOCIN",&XLOCIN);
   AncSvx->SetBranchAddress("YLOCIN",&YLOCIN);
   AncSvx->SetBranchAddress("ZLOCIN",&ZLOCIN);
   AncSvx->SetBranchAddress("XLOCOUT",&XLOCOUT);
   AncSvx->SetBranchAddress("YLOCOUT",&YLOCOUT);
   AncSvx->SetBranchAddress("ZLOCOUT",&ZLOCOUT);
   AncSvx->SetBranchAddress("HITVOL0",&HITVOL0);
   AncSvx->SetBranchAddress("HITVOL1",&HITVOL1);
   AncSvx->SetBranchAddress("HITVOL2",&HITVOL2);
   AncSvx->SetBranchAddress("HITVOL3",&HITVOL3);
   AncSvx->SetBranchAddress("HITVOL4",&HITVOL4);
   AncSvx->SetBranchAddress("HITVOL5",&HITVOL5);
   AncSvx->SetBranchAddress("NHIT",&NHIT);
   AncSvx->SetBranchAddress("X0_EVENT",&X0_EVENT);
   AncSvx->SetBranchAddress("Y0_EVENT",&Y0_EVENT);
   AncSvx->SetBranchAddress("Z0_EVENT",&Z0_EVENT);
   AncSvx->SetBranchAddress("B_IMPACT",&B_IMPACT);
   AncSvx->SetBranchAddress("EVENT",&EVENT);

  Int_t nentries = AncSvx->GetEntries();
  cout <<" AncSvx nentries = "<<nentries<<endl;

//xxx arrays of floats
Float_t event[nentries];
Float_t track[nentries];
Float_t ptot[nentries];
Float_t ptheta[nentries];
Float_t pphi[nentries];
Float_t theta[nentries];
Float_t phi[nentries];
Float_t x[nentries];
Float_t y[nentries];
Float_t z[nentries];
Float_t px[nentries];
Float_t py[nentries];
Float_t pz[nentries];
Float_t layer[nentries];
Float_t id[nentries];
Int_t ievent[10001];

//     This is the loop skeleton
//       To read only selected branches, Insert statements like:
// AncSvx->SetBranchStatus("*",0);  // disable all branches
// TTreePlayer->SetBranchStatus("branchname",1);  // activate branchname

  Int_t nbytes = 0;
  Float_t oldevent = 1.;
  Int_t iev = 0;
  ievent[1]=0;

  for (Int_t i=0; i<nentries;i++) {
    nbytes += AncSvx->GetEntry(i);
    event[i]=EVENT;           // watch for new event number:
    if (EVENT>oldevent) {
      oldevent=EVENT;
      iev=EVENT;
      ievent[iev]=i;
      cout<<"iev:"<<iev<<", ievent[iev]="<<ievent[iev]<<endl; 
      if (iev%10==0) {cout <<"record "<<i<<", event "<<iev<<" "<<endl;}
    }
    track[i]=TRACK;           // fill hit arrays
    layer[i] =  LAYER;
    ptot[i]=PTOT;
    ptheta[i]=PTHETA;
    pphi[i]=PPHI;
    theta[i]=THETA;
    phi[i]=PHI;
    x[i]=XGLOBAL;
    y[i]=YGLOBAL;
    z[i]=ZGLOBAL;
    px[i]=PMOMX;
    py[i]=PMOMY;
    pz[i]=PMOMZ;
    id[i]=IDPART;  

  }                           // end loop to stuff arrays
  cout <<"Total records "<<i<<", events "<<iev<<" "<<endl;
//---------------------------------------------------------------------------------

const int max_tracks=200;
Float_t tracks[max_tracks][13][7];    // track#(1-2) layer(1-12) (x,y,z,px,py,pz)
Int_t idpart  [max_tracks];           // pid   (1-2)
Float_t vertex[max_tracks][7];        // track#(1-2) (xvertex,yvertex,zvertex,ptot,ptheta,pphi)

Int_t itrack = 0;
Int_t ilayer = 0;
Float_t pzvertex = 0.;
Float_t ptvertex = 0.;
Float_t pxvertex = 0.;
Float_t pyvertex = 0.;
Float_t degtorad = 3.14159/180.;

//define histograms
TH1F *g01 = new TH1F("g01","layer", 13, 0.0, 13.0);
TH2F *g02 = new TH2F("g02","layer vs r",13,0.0,13.0, 80, 0.0,20.0);
TH2F *g03 = new TH2F("g03","z vs r",80, -40.0,40.0, 80, 0.0,20.0);
TH2F *g04 = new TH2F("g04","z vs r, layer 5",80, -40.0,40.0, 80, 0.0,20.0);
TH2F *g05 = new TH2F("g05","z vs r, layer 5",80, -40.0,40.0, 80, 0.0,20.0);
TH2F *g06 = new TH2F("g06","z vs r, layer 5",80, -40.0,40.0, 80, 0.0,20.0);
TH2F *g07 = new TH2F("g07","z vs r, layer 5",80, -40.0,40.0, 80, 0.0,20.0);
TH2F *g08 = new TH2F("g08","z vs r, layer 5",80, -40.0,40.0, 80, 0.0,20.0);
g04->SetMarkerColor(2);
g05->SetMarkerColor(3);
g06->SetMarkerColor(4);
g07->SetMarkerColor(5);
g08->SetMarkerColor(6);

g04->SetLineColor(1);
g05->SetLineColor(2);
g06->SetLineColor(4);
g07->SetLineColor(6);
g08->SetLineColor(3);

g04->SetLineWidth(5);
g05->SetLineWidth(5);
g06->SetLineWidth(5);
g07->SetLineWidth(5);
g08->SetLineWidth(5);

int nmess = 1;                          // messages every nmess events (typically 1 or 100)
iev = 100;                              //xxx <10K uncomment for testing
cout <<"Now processing "<<iev<<" events:"<<endl;
for (Int_t ii=1; ii<iev; ii++) {         // loop over events

  for (Int_t i=1; i<=6; i++) {            // first clear track and vertex arrays
    vertex[1][i]=0;
    vertex[2][i]=0;
    for (Int_t j=1; j<=12;j++) {
      tracks[1][j][i]=0.;
      tracks[2][j][i]=0.;
    }
  }

  //work on each event's hits, filling tracks and vertex arrays
  for (Int_t i=ievent[ii]; i<ievent[ii+1]; i++) {
    itrack=track[i];
    ilayer=layer[i];
    // cout<<" itrack, ilayer= "<<itrack<<"  "<<ilayer<<endl;
    if (itrack>=200) {
      cout <<" Track " <<itrack <<" ";
      cout <<" layer " <<ilayer <<" (too many tracks) "<<endl;
      //break;   should be 'continue' ?
    }
    else {                                      // store track info 
      idpart[itrack]=id[i];

      tracks[itrack][ilayer][1] = x[i];
      tracks[itrack][ilayer][2] = y[i];
      tracks[itrack][ilayer][3] = z[i];
      tracks[itrack][ilayer][4] = px[i];
      tracks[itrack][ilayer][5] = py[i];
      tracks[itrack][ilayer][6] = pz[i];

      // store vertex info
      pzvertex = ptot[i]*cos(ptheta[i]*degtorad);
      ptvertex = ptot[i]*sin(ptheta[i]*degtorad);
      pxvertex = ptvertex*cos(pphi[i]*degtorad);
      pyvertex = ptvertex*sin(pphi[i]*degtorad);
      //vertex[itrack][1] = vx[ii][6];            note vx,y,z are note defined
      //vertex[itrack][2] = vy[ii][6];            or set
      //vertex[itrack][3] = vz[ii][6];
      vertex[itrack][4] = pxvertex;
      vertex[itrack][5] = pyvertex;
      vertex[itrack][6] = pzvertex;

      //printf("%d %d %d %d %6.1f = i,itrack,ilayer,idpart[i],z[i]\n",i,itrack,ilayer,idpart[i],z[i]);
      g01->Fill(layer[i]);
      g02->Fill(layer[i],sqrt(x[i]**2+y[i]**2));
      g03->Fill(z[i],sqrt(x[i]**2+y[i]**2));
      if (ilayer == 5) g04->Fill(z[i],sqrt(x[i]**2+y[i]**2));
      if (ilayer == 6) g05->Fill(z[i],sqrt(x[i]**2+y[i]**2));
      if (ilayer == 7) g06->Fill(z[i],sqrt(x[i]**2+y[i]**2));
      if (ilayer == 8) g07->Fill(z[i],sqrt(x[i]**2+y[i]**2));
      if (ilayer  < 5) g08->Fill(z[i],sqrt(x[i]**2+y[i]**2));
    }      // <200 track in the event
  }
}
  //-------------------------------------------------------------------------//

}