/*
$Log$
+Revision 1.14 2001/01/26 19:58:48 hristov
+Major upgrade of AliRoot code
+
+Revision 1.13 2000/12/21 15:30:18 fca
+Correcting coding convention violations
+
+Revision 1.12 2000/12/12 13:18:59 hristov
+Protection against FPE
+
+Revision 1.11 2000/11/30 07:12:48 alibrary
+Introducing new Rndm and QA classes
+
+Revision 1.10 2000/10/02 21:28:14 fca
+Removal of useless dependecies via forward declarations
+
Revision 1.9 2000/07/13 16:19:09 fca
Mainly coding conventions + some small bug fixes
// //
//////////////////////////////////////////////////////////////////////////
-#include <TROOT.h>
#include <TTree.h>
#include <TButton.h>
#include <TCanvas.h>
#include <TView.h>
-#include <TText.h>
-#include <TPolyMarker3D.h>
#include <TPaveLabel.h>
#include <TPaveText.h>
-#include <TList.h>
#include <TDiamond.h>
-#include <TNode.h>
#include <TArc.h>
-#include <TTUBE.h>
#include <TSlider.h>
#include <TSliderBox.h>
#include <TGaxis.h>
#include <TVirtualX.h>
#include <TMath.h>
-#include <X3DBuffer.h>
#include "AliRun.h"
#include "AliDetector.h"
#include "AliDisplay.h"
#include "AliPoints.h"
#include "TParticle.h"
-#include "AliGeometry.h"
#include "TGeometry.h"
+#include "AliHeader.h"
static const Float_t kptcutmax = 2;
static const Float_t ketacutmax = 1.5;
}
//----------------------------------------------------------------------------
-void AliDisplay::ShowTrack(Int_t idx) {
+void AliDisplay::ShowTrack(Int_t idx)
+{
+ //
+ // Display track idx
+ //
AliDetector *mTPC=(AliDetector*)gAlice->GetModule("TPC");
TObjArray *points=mTPC->Points();
int ntracks=points->GetEntriesFast();
pm->Draw("same");
// fPad->Update();
// fPad->Modified();
- TClonesArray *particles=gAlice->Particles();
- TParticle *p = (TParticle*)particles->UncheckedAt(idx);
+ // TClonesArray *particles=gAlice->Particles();
+ // TParticle *p = (TParticle*)particles->UncheckedAt(idx);
+ TParticle *p = gAlice->Particle(idx);
printf("\nTrack index %d\n",idx);
printf("Particle ID %d\n",p->GetPdgCode());
printf("Parent %d\n",p->GetFirstMother());
pxyz = pm->GetP();
r = TMath::Sqrt(pxyz[0]*pxyz[0] + pxyz[1]*pxyz[1]);
theta = TMath::ATan2(r,TMath::Abs(pxyz[2]));
- if(theta) eta = -TMath::Log(TMath::Tan(0.5*theta)); else eta = 1e10;
+ if(theta) eta = -TMath::Log(TMath::Abs(TMath::Tan(0.5*theta))); else eta = 1e10;
if (pxyz[2] < 0) eta = -eta;
if (eta < etamin || eta > etamax) continue;
pm->Draw();