+// Base class AliAODPairBaseCut:
+// This class defines the range of some property - pure virtual
+// Property is coded by AliAODCutTypes type
+// Derived classes:
+// AliAODQInvCut
+// AliAODKtCut
+// AliAODKStarCut
+// AliAODKStarOutCut
+// AliAODKStarSideCut
+// AliAODKStarLongCut
+// AliAODQSideLCMSCut
+// AliAODQOutLCMSCut
+// AliAODQLongLCMSCut
+// AliAODDeltaECut
+// AliAODDeltaPCut
+// AliAODDeltaPvectorCut
+// AliAODDeltaPhiCut
+// AliAODDeltaThetaCut
+// AliAODCluterOverlapCut
+// AliAODAvSeparationCut
+// AliAODSeparationCut
+// AliAODITSSeparationCut
+// AliAODOutSideSameSignCut
+// AliAODOutSideDiffSignCut
+// AliAODLogicalOperPairCut
+// AliAODOrPairCut
+// AliAODAndPairCut
+// Author: Piotr.Skowronski@cern.ch
#include "AliAODPairBaseCut.h"
#include "AliTrackPoints.h"
#ifndef ALIAODPAIRBASECUT_H
#define ALIAODPAIRBASECUT_H
-
-#include <TNamed.h>
+// Base class AliAODPairBaseCut:
+// This class defines the range of some property - pure virtual
+// Property is coded by AliAODCutTypes type
+// Derived classes:
+// AliAODQInvCut
+// AliAODKtCut
+// AliAODKStarCut
+// AliAODKStarOutCut
+// AliAODKStarSideCut
+// AliAODKStarLongCut
+// AliAODQSideLCMSCut
+// AliAODQOutLCMSCut
+// AliAODQLongLCMSCut
+// AliAODDeltaECut
+// AliAODDeltaPCut
+// AliAODDeltaPvectorCut
+// AliAODDeltaPhiCut
+// AliAODDeltaThetaCut
+// AliAODCluterOverlapCut
+// AliAODAvSeparationCut
+// AliAODSeparationCut
+// AliAODITSSeparationCut
+// AliAODOutSideSameSignCut
+// AliAODOutSideDiffSignCut
+// AliAODLogicalOperPairCut
+// AliAODOrPairCut
+// AliAODAndPairCut
+
+#include <TObject.h>
#include "AliAODPair.h"
class AliAODPairBaseCut: public TObject
{
- //This class defines the range of some property - pure virtual
- //Property is coded by AliAODCutTypes type
-
public:
enum EAODPairCutProperty
#include "AliAODPairCut.h"
/* $Id$ */
-//____________________________________
+
/////////////////////////////////////////////////////////////////////////
//
// Class AliAODPairCut:
-//
// implements cut on the pair of particles
// more info: http://aliweb.cern.ch/people/skowron/analyzer/index.html
// Author: Piotr.Skowronski@cern.ch
//-------------------------------------------------------------------
+
#include "AliAODPair.h"
#include "AliAODParticleCut.h"
-#include "AliTrackPoints.h"
-#include "AliClusterMap.h"
+//#include "AliTrackPoints.h"
+//#include "AliClusterMap.h"
ClassImp(AliAODPairCut)
const Int_t AliAODPairCut::fgkMaxCuts = 50;
//Class implements cut on the pair of particles
//
//more info: http://aliweb.cern.ch/people/skowron/analyzer/index.html
+
#include <TNamed.h>
#include "AliAODPairBaseCut.h"
/////////////////////////////////////////////////////////////
#include "AliVAODParticle.h"
-//#include <TLorentzVector.h>
-//#include <TMath.h>
+#include <TLorentzVector.h>
+#include <TMath.h>
#include <TDatabasePDG.h>
//-------------------------------------------------------------------------
// Implementation of the Analysis Oriented Data (AOD) V0 vertex class
-//
// Origin: B.Hippolyte, IReS, hippolyt@in2p3.fr
// G.Van Buren, BNL, gene@bnl.gov (original STAR MuDsts)
-//
// Purpose: Having observables for physics available for V0s
//-------------------------------------------------------------------------
-#include <Riostream.h>
+
#include <TMath.h>
#include "AliESD.h"
this->Fill(rV0Vertex,rEvent);
}
-// AliAODv0::~AliAODv0(){
-// }
-
-void AliAODv0::Fill(AliESDv0* rV0Vertex ,AliESD* rEvent){// Filling method
+void AliAODv0::Fill(AliESDv0* rV0Vertex ,AliESD* rEvent){
+ // Fills the data memebers of the AOD
fEvent=rEvent;
Double_t tDecayVertexV0[3]; rV0Vertex->GetXYZ(tDecayVertexV0[0],tDecayVertexV0[1],tDecayVertexV0[2]);
fDecayVertexV0X = tDecayVertexV0[0];
fDcaNegToPrimVertex = TMath::Sqrt(tDcaNegToPrimVertex[0]*tDcaNegToPrimVertex[0]+tDcaNegToPrimVertex[1]*tDcaPosToPrimVertex[1]);
}
-void AliAODv0::ResetV0(){// Reset method
+void AliAODv0::ResetV0(){
+ // Sets the default values of the AOD data members
fDecayVertexV0X = 999;
fDecayVertexV0Y = 999;
fDecayVertexV0Z = 999;
//-------------------------------------------------------------------------
// Implementation of the Analysis Oriented Data (AOD) V0 vertex class
-//
// Origin: B.Hippolyte, IReS, hippolyt@in2p3.fr
// G.Van Buren, BNL, gene@bnl.gov (original STAR MuDsts)
-//
// Purpose: Having observables for physics available for V0s
//-------------------------------------------------------------------------
/* $Id$ */
-//_________________________________________________
-///////////////////////////////////////////////////
+///////////////////////////////////////////////////////////////////////////
//
// class AliClusterMap
-//
// class that describes cluster occupation at TPC
// Each padraw has a corresponding bit in fPadRawMap
-//
-//
// more info: http://aliweb.cern.ch/people/skowron/analyzer/index.html
// Piotr.Skowronski@cern.ch
//
#include "AliESDtrack.h"
#include "AliLog.h"
#include "AliTPCtrack.h"
-#include "AliVAODParticle.h"
-const Int_t AliClusterMap::fNPadRows = 159;
+const Int_t AliClusterMap::fgkNPadRows = 159;
ClassImp(AliClusterMap)
AliClusterMap::AliClusterMap():
- fPadRawMap(fNPadRows)
+ fPadRawMap(fgkNPadRows)
{
//ctor
}
/***********************************************************************/
AliClusterMap::AliClusterMap(AliESDtrack* track):
- fPadRawMap( (track)?track->GetTPCClusterMap():fNPadRows )
+ fPadRawMap( (track)?track->GetTPCClusterMap():fgkNPadRows )
{
//ctor
/***********************************************************************/
AliClusterMap::AliClusterMap(AliTPCtrack* track):
- fPadRawMap(fNPadRows)
+ fPadRawMap(fgkNPadRows)
{
//ctor
{
//Prints the bit map
TString msg;
- for ( Int_t i = 0; i < fNPadRows; i++)
+ for ( Int_t i = 0; i < fgkNPadRows; i++)
{
if ( fPadRawMap.TestBitNumber(i) )
{
Int_t nh = 0;
Int_t an = 0;
- for ( Int_t i = 0; i < fNPadRows; i++)
+ for ( Int_t i = 0; i < fgkNPadRows; i++)
{
Bool_t x = HasClAtPadRow(i);
Bool_t y = clmap.HasClAtPadRow(i);
#ifndef ALICLUSTERMAP_H
#define ALICLUSTERMAP_H
-//_________________________________________________
-///////////////////////////////////////////////////
+
+///////////////////////////////////////////////////////////////////////////
//
// class AliClusterMap
-//
// class that describes cluster occupation at TPC
// Each padraw has a corresponding bit in fPadRawMap
-//
-//
// more info: http://aliweb.cern.ch/people/skowron/analyzer/index.html
// Piotr.Skowronski@cern.ch
//
protected:
private:
TBits fPadRawMap;//bit vector of length 150 correspondind to total number of padraws in TPC
- static const Int_t fNPadRows;
+ static const Int_t fgkNPadRows; // Number of pad rows
ClassDef(AliClusterMap,1)
};
#include "AliEventCut.h"
-//________________________________
+
///////////////////////////////////////////////////////////
//
-// class AliRunAnalysis
-//
-//
-//
+// class AliEventCut
//
+// Event cut. It has list of base event cuts.
+// Each of base event cut checks only one property.
+// Logical base cuts also exists that point to other base cuts.
+// Using them one can build complicated cut with binary tree structure
+// Author: Piotr.Skowronski@cern.ch
///////////////////////////////////////////////////////////
-#include <TObjArray.h>
-
#include "AliEventBaseCut.h"
ClassImp(AliEventCut)
#ifndef ALIEVENTCUT_H
#define ALIEVENTCUT_H
-//________________________________
+
///////////////////////////////////////////////////////////
//
// class AliEventCut
// Each of base event cut checks only one property.
// Logical base cuts also exists that point to other base cuts.
// Using them one can build complicated cut with binary tree structure
-//
+// Author: Piotr.Skowronski@cern.ch
///////////////////////////////////////////////////////////
#include <TObject.h>
protected:
AliEventBaseCut* FindCut(AliEventBaseCut::EEventCutProperty prop);
- TObjArray fBaseCuts;
+ TObjArray fBaseCuts; // Array of cuts
private:
ClassDef(AliEventCut,1)
};
#include "AliFlowAnalysis.h"
-//________________________________
-///////////////////////////////////////////////////////////
-//
+
+//*********************************************************
// class AliFlowAnalysis
-//
// Flow Analysis
-//
-//
// S.Radomski@gsi.de
// Piotr.Skowronski@cern.ch
-//
-///////////////////////////////////////////////////////////
-/*********************************************************/
-
-#include <TString.h>
-#include <TParticle.h>
+//*********************************************************
-#include <AliStack.h>
#include <AliAOD.h>
#include <AliVAODParticle.h>
#include <AliAODParticleCut.h>
Int_t AliFlowAnalysis::ProcessEvent(AliAOD* aodrec, AliAOD* aodsim)
{
-
+ // Process AOD events containing the reconstructed and simulated information
Info("ProcessEvent","Sim AOD address %#x",aodsim);
Double_t psi = 0, v2 = 0;
if (aodrec)
#ifndef ALIFLOWANALYSIS_H
#define ALIFLOWANALYSIS_H
-//________________________________
-///////////////////////////////////////////////////////////
-//
+
+//*********************************************************
// class AliFlowAnalysis
-//
// Flow Analysis
-//
-//
// S.Radomski@gsi.de
// Piotr.Skowronski@cern.ch
-//
-///////////////////////////////////////////////////////////
+//*********************************************************
#include "AliAnalysis.h"
#include <TFile.h>
#include <TGliteXmlEventlist.h>
#include <TH1.h>
-#include <TKey.h>
-#include <TObjString.h>
#include <TPDGCode.h>
#include <TParticle.h>
#include <TString.h>
-#include <TTree.h>
-
-#include "AliAnalysis.h"
-#include "AliAODRun.h"
#include "AliAOD.h"
#include "AliAODParticle.h"
-#include "AliAODParticleCut.h"
-#include "AliAODRun.h"
-#include "AliAnalysis.h"
#include "AliClusterMap.h"
#include "AliESD.h"
#include "AliESDtrack.h"
-#include "AliKalmanTrack.h"
#include "AliLog.h"
#include "AliReaderESD.h"
-#include "AliRun.h"
#include "AliRunLoader.h"
#include "AliStack.h"
Int_t AliReaderESD::ReadESDCentral(AliESD* esd)
{
//****** Tentative particle type "concentrations"
- static const Double_t concentr[5]={0.05, 0., 0.85, 0.10, 0.05};
+ static const Double_t kConcentr[5]={0.05, 0., 0.85, 0.10, 0.05};
Double_t pidtable[kNSpecies];//array used for reading pid probabilities from ESD track
Double_t w[kNSpecies];
//Here we apply Bayes' formula
Double_t rc=0.;
- for (Int_t s=0; s<AliPID::kSPECIES; s++) rc+=concentr[s]*pidtable[s];
+ for (Int_t s=0; s<AliPID::kSPECIES; s++) rc+=kConcentr[s]*pidtable[s];
if (rc==0.0)
{
AliDebug(3,"Particle rejected since total bayessian PID probab. is zero.");
continue;
}
- for (Int_t s=0; s<AliPID::kSPECIES; s++) w[s]=concentr[s]*pidtable[s]/rc;
+ for (Int_t s=0; s<AliPID::kSPECIES; s++) w[s]=kConcentr[s]*pidtable[s]/rc;
if (AliDebugLevel() > 4)
{
/**********************************************************/
Int_t AliReaderESD::ReadESDMuon(AliESD* esd)
{
-
+ // Reads the muon tracks from the ESD
Double_t vertexpos[3];//vertex position, assuming no secondary decay
const AliESDVertex* vertex = esd->GetVertex();
AliDebug(1,Form("Reading Event %d \nFound %d tracks.",fCurrentEvent,nTracks));
// settings
- Float_t Chi2Cut = 100.;
- Float_t PtCutMin = 1.;
- Float_t PtCutMax = 10000.;
+ Float_t chi2Cut = 100.;
+ Float_t ptCutMin = 1.;
+ Float_t ptCutMax = 10000.;
Float_t muonMass = 0.105658389;
Int_t pdgcode = -13;
Double_t thetaX, thetaY, pYZ;
- Double_t pxRec1, pyRec1, pzRec1, E1;
+ Double_t pxRec1, pyRec1, pzRec1, e1;
Int_t charge;
Int_t ntrackhits;
pxRec1 = pzRec1 * TMath::Tan(thetaX);
pyRec1 = pzRec1 * TMath::Tan(thetaY);
charge = Int_t(TMath::Sign(1.,muonTrack->GetInverseBendingMomentum()));
- E1 = TMath::Sqrt(muonMass * muonMass + pxRec1 * pxRec1 + pyRec1 * pyRec1 + pzRec1 * pzRec1);
- fV1.SetPxPyPzE(pxRec1, pyRec1, pzRec1, E1);
+ e1 = TMath::Sqrt(muonMass * muonMass + pxRec1 * pxRec1 + pyRec1 * pyRec1 + pzRec1 * pzRec1);
+ fV1.SetPxPyPzE(pxRec1, pyRec1, pzRec1, e1);
ntrackhits = muonTrack->GetNHit();
fitfmin = muonTrack->GetChi2();
// chi2 per d.o.f.
Float_t ch1 = fitfmin / (2.0 * ntrackhits - 5);
- if ((ch1 < Chi2Cut) && (pt1 > PtCutMin) && (pt1 < PtCutMax)) {
+ if ((ch1 < chi2Cut) && (pt1 > ptCutMin) && (pt1 < ptCutMax)) {
AliAODParticle* track = new AliAODParticle(pdgcode*charge,1,iTrack,
- pxRec1, pyRec1,pzRec1, E1,
+ pxRec1, pyRec1,pzRec1, e1,
vertexpos[0], vertexpos[1], vertexpos[2], 0.);
fEventRec->AddParticle(track);
}
/* $Id$ */
-//_______________________________________________________________________
-/////////////////////////////////////////////////////////////////////////
-//
+//***********************************************************************
// class AliReaderESDTree
-//
// Reader for MUON ESD Tree (only for rec)
-//
// finck@subatech.in2p3.fr
-//
-/////////////////////////////////////////////////////////////////////////
+//***********************************************************************
#include <TString.h>
#include <TTree.h>
#include "AliESD.h"
#include "AliLog.h"
#include "AliReaderESDTree.h"
-#include "AliRun.h"
#include "AliRunLoader.h"
+
ClassImp(AliReaderESDTree)
AliReaderESDTree::AliReaderESDTree(const Char_t* esdfilename, const Char_t* galfilename):
#ifndef AliReaderESDTree_H
#define AliReaderESDTree_H
-//_______________________________________________________________________
-/////////////////////////////////////////////////////////////////////////
-//
+
+//***********************************************************************
// class AliReaderESDTree
-//
// Reader for ESD Tree
-//
// Ch. Finck
-//
-/////////////////////////////////////////////////////////////////////////
+//***********************************************************************
+
#include "AliReaderESD.h"
-#include <TString.h>
class TFile;
class TTree;
/* $Id$ */
+//********************************************************
+// class AliRunAnalysis *
+// Analysis manager *
+// Author: Piotr.Skowronski@cern.ch *
+//********************************************************
+
#include "AliRunAnalysis.h"
#include "AliLog.h"
-
-//________________________________
-///////////////////////////////////////////////////////////
-//
-// class AliRunAnalysis
-//
-//
-//
-// Piotr.Skowronski@cern.ch
-//
-///////////////////////////////////////////////////////////
-
-#include <stdlib.h>
-
-#include <TString.h>
-#include <TObjString.h>
-#include <TClass.h>
-#include <TFile.h>
-#include <TKey.h>
-#include <TObjArray.h>
-
-#include <AliRun.h>
-#include <AliRunLoader.h>
-#include <AliStack.h>
-#include <AliESDtrack.h>
-#include <AliESD.h>
-
-
+#include "AliAnalysis.h"
#include "AliEventCut.h"
#include "AliReader.h"
-#include "AliVAODParticle.h"
ClassImp(AliRunAnalysis)
#ifndef ALIRUNANALYSIS_H
#define ALIRUNANALYSIS_H
-//________________________________
+
///////////////////////////////////////////////////////////
//
// class AliRunAnalysis
-//
// Analysis manager
-//
-//
-// Piotr.Skowronski@cern.ch
+// Author: Piotr.Skowronski@cern.ch
//
///////////////////////////////////////////////////////////
#include <TTask.h>
#include <TObjArray.h>
-#include "AliAnalysis.h"
class AliEventCut;
-class TFile;
class AliReader;
+class AliAnalysis;
+class AliAOD;
class AliRunAnalysis: public TTask
{
/* $Id$ */
-#include "AliTrackPoints.h"
-#include "AliLog.h"
-
-//_________________________________
////////////////////////////////////////////////////////////
// //
// class AliTrackPoints //
// //
////////////////////////////////////////////////////////////
+#include "AliTrackPoints.h"
+#include "AliLog.h"
+
#include <TClonesArray.h>
#include <TFile.h>
#include <TMath.h>
#include "AliTrackReference.h"
#include "AliITStrackV2.h"
+#include "AliRun.h"
+#include "AliESD.h"
+#include "AliRunLoader.h"
+#include "AliTPCtrack.h"
+#include "TTree.h"
+#include "TBranch.h"
+#include "TH2D.h"
+
ClassImp(AliTrackPoints)
Int_t AliTrackPoints::fgDebug = 0;
{
//constructor
}
+
/***************************************************************/
AliTrackPoints::AliTrackPoints(AliTrackPoints::ETypes type, AliESDtrack* track, Float_t mf):
}
}
+
/***************************************************************/
AliTrackPoints::AliTrackPoints(Int_t n, AliESDtrack* track, Float_t mf, Float_t dr, Float_t r0):
MakePoints(dr,r0,x,par,c,alpha);
}
+
/***************************************************************/
AliTrackPoints::AliTrackPoints(Int_t n, AliTPCtrack* track, Float_t dr, Float_t r0):
Double_t alpha = track->GetAlpha();
Double_t c=track->GetC();
MakePoints(dr,r0,x,par,c,alpha);
-}
+}
+
/***************************************************************/
AliTrackPoints::~AliTrackPoints()
delete [] fY;
delete [] fZ;
}
+
/***************************************************************/
void AliTrackPoints::MakePoints( Float_t dr, Float_t r0, Double_t x, Double_t* par, Double_t c, Double_t alpha)
fX[i],fY[i],fZ[i],rc,TMath::Hypot(fX[i],fY[i])));
}
}
+
/***************************************************************/
void AliTrackPoints::MakeITSPoints(AliESDtrack* track)
// z=R*Pz/Pt
AliITStrackV2 itstrack(*track,kTRUE);
Double_t x,y,z;
- static const Double_t r[6] = {4.0, 7.0, 14.9, 23.8, 39.1, 43.6};
+ static const Double_t kR[6] = {4.0, 7.0, 14.9, 23.8, 39.1, 43.6};
for (Int_t i = 0; i < 6; i++)
{
- itstrack.GetGlobalXYZat(r[i],x,y,z);
+ itstrack.GetGlobalXYZat(kR[i],x,y,z);
fX[i] = x;
fY[i] = y;
fZ[i] = z;
}
/***************************************************************/
+
void AliTrackPoints::MakeITSPointsInnerFromVertexOuterFromTPC(AliESDtrack* track, Float_t mf)
{
//makes trackpoints for ITS
//for 3 inner layers calculates out of the vector at vertex
//for 3 outer ---------------//------------------ at inner TPC
- static const Double_t r[6] = {4.0, 7.0, 14.9, 23.8, 39.1, 43.6};
+ static const Double_t kR[6] = {4.0, 7.0, 14.9, 23.8, 39.1, 43.6};
AliITStrackV2 itstrack(*track,kTRUE);
Double_t x,y,z;
for (Int_t i = 0; i < 3; i++)
{
- itstrack.GetGlobalXYZat(r[i],x,y,z);
+ itstrack.GetGlobalXYZat(kR[i],x,y,z);
fX[i] = x;
fY[i] = y;
fZ[i] = z;
AliDebug(3,Form("X %f Y %f Z %f R asked %f R obtained %f",
- fX[i],fY[i],fZ[i],r[i],TMath::Hypot(fX[i],fY[i])));
+ fX[i],fY[i],fZ[i],kR[i],TMath::Hypot(fX[i],fY[i])));
}
for (Int_t i = 3; i < 6; i++)
{
Float_t ax,ay,az;
- AliTrackPoints tmptp(1,track,mf,0,r[i]);
+ AliTrackPoints tmptp(1,track,mf,0,kR[i]);
tmptp.PositionAt(0,ax,ay,az);
fX[i] = ax;
fY[i] = ay;
fZ[i] = az;
AliDebug(3,Form("X %f Y %f Z %f R asked %f R obtained %f",
- fX[i],fY[i],fZ[i],r[i],TMath::Hypot(fX[i],fY[i])));
+ fX[i],fY[i],fZ[i],kR[i],TMath::Hypot(fX[i],fY[i])));
}
}
-
/***************************************************************/
void AliTrackPoints::PositionAt(Int_t n, Float_t &x,Float_t &y,Float_t &z)
AliDebug(2,Form("n %d; X %f; Y %f; Z %f",n,x,y,z));
}
+
/***************************************************************/
void AliTrackPoints::Move(Float_t x, Float_t y, Float_t z)
fZ[i]+=z;
}
}
+
/***************************************************************/
Double_t AliTrackPoints::AvarageDistance(const AliTrackPoints& tr)
return retval;
}
+
/***************************************************************/
void AliTrackPoints::Print(Option_t* /*option*/) const
{
-
+ // Prints the coordinates of the track points
Info("Print","There is %d points",fN);
for(Int_t i = 0; i < fN; i++)
{
}
/***************************************************************/
-/***************************************************************/
-/***************************************************************/
-/***************************************************************/
-/***************************************************************/
-/***************************************************************/
-
-#include "AliRun.h"
-#include "AliESD.h"
-#include "AliRunLoader.h"
-#include "AliTPCtrack.h"
-#include "TTree.h"
-#include "TBranch.h"
-#include "TH2D.h"
-#include "TCanvas.h"
-#include "AliMagF.h"
-
-
-
void AliTrackPoints::Testesd(Int_t entr,const char* fname )
{
+ // This is a former macro converted to member function.
+ // It is used to test the functionality of the class using ESD tracks
delete gAlice;
gAlice = 0x0;
AliRunLoader* rl = AliRunLoader::Open();
}
- Int_t N = 170;
- AliTrackPoints* tp = new AliTrackPoints(N,t,mf,1.);
+ Int_t nTrackPoints = 170;
+ AliTrackPoints* tp = new AliTrackPoints(nTrackPoints,t,mf,1.);
Float_t xmin = -250;
Float_t xmax = 250;
Float_t x,y,z;
- for (Int_t i = 0;i<N;i++)
+ for (Int_t i = 0;i<nTrackPoints;i++)
{
Double_t r = 84.1+i;
tp->PositionAt(i,x,y,z);
delete rl;
}
-/***************************************************************/
-/***************************************************************/
/***************************************************************/
void AliTrackPoints::Testtpc(Int_t entr)
{
+ // This is a former macro converted to member function.
+ // It is used to test the functionality of the class using TPC tracks
delete gAlice;
gAlice = 0x0;
AliRunLoader* rl = AliRunLoader::Open();
TBranch* b=l->TreeT()->GetBranch("tracks");
b->SetAddress(&t);
l->TreeT()->GetEntry(entr);
- Int_t N = 160;
- AliTrackPoints* tp = new AliTrackPoints(N,t,1.);
+ Int_t nTrackPoints = 160;
+ AliTrackPoints* tp = new AliTrackPoints(nTrackPoints,t,1.);
Float_t xmin = -250;
Float_t xmax = 250;
Float_t x,y,z;
- for (Int_t i = 0;i<N;i++)
+ for (Int_t i = 0;i<nTrackPoints;i++)
{
Double_t r = 84.1+i;
tp->PositionAt(i,x,y,z);
/* $Id$ */
-/////////////////////////////////////////////////////////////
-//
+//***********************************************************
+// Class AliVAODParticle
// base class for AOD particles
-//
-/////////////////////////////////////////////////////////////
+// Author: Piotr.Skowronski@cern.ch
+//***********************************************************
#include <TObject.h>
-#include <TLorentzVector.h>
+//#include <TLorentzVector.h>
#include <TVector3.h>
class AliTrackPoints;
//twise, f.g. as a pion and as a kaon than both have the same UID)
virtual void SetUID(Int_t /*id*/){/* *this */}
// type information
- virtual Bool_t IsSimulated() {return kFALSE;};
- virtual Bool_t IsTrack() {return kFALSE;};
- virtual Bool_t IsCluster() {return kFALSE;};
+ virtual Bool_t IsSimulated() const {return kFALSE;};
+ virtual Bool_t IsTrack() const {return kFALSE;};
+ virtual Bool_t IsCluster() const {return kFALSE;};
//HBT specific
virtual AliTrackPoints* GetTPCTrackPoints() const {return 0x0;}