[u/mrichter/AliRoot.git] / STEER / AliGenInfo.h



////////////////////////////////////////////////////////////////////////
//
// Start of implementation of the class digitRow
//
////////////////////////////////////////////////////////////////////////
const Int_t kgRowBytes = 32;

class digitRow: public TObject {

public:
  digitRow();
  virtual ~digitRow(){;}
  void SetRow(Int_t row);
  Bool_t TestRow(Int_t row);
  digitRow & operator=(const digitRow &digOld);
  Int_t RowsOn(Int_t upto=8*kgRowBytes);
  Int_t Last();
  Int_t First();
  void Reset();

//private:
  UChar_t fDig[kgRowBytes];
  ClassDef(digitRow,1)  // container for digit pattern
};
ClassImp(digitRow)


////////////////////////////////////////////////////////////////////////
//
// Start of implementation of the class AliMCInfo
//
////////////////////////////////////////////////////////////////////////

class AliMCInfo: public TObject {

public:
  AliMCInfo();
  ~AliMCInfo();
  void Update();

  AliTrackReference  fTrackRef;      // track reference saved in the output tree
  AliTrackReference  fTrackRefOut;   // decay track reference saved in the output tree
  AliTrackReference  fTRdecay;       // track reference at decay point
  //
  Int_t fPrimPart;                   // index of primary particle in TreeH
  TParticle fParticle;           // generated particle 
  Float_t   fMass;               // mass of the particle
  Float_t fCharge;               //
  Int_t fLabel;                   // track label
  Int_t fEventNr;                 // event number
  Int_t fMCtracks;                // indication of how many times the track is retuturned back
  Int_t fPdg;                     //pdg code
  Float_t fDecayCoord[3];         // position of particle decay
  Double_t fVDist[4];             //distance of the particle vertex from primary vertex
  Bool_t fTPCdecay;               //indicates decay in TPC
  Int_t fRowsWithDigitsInn;    // number of rows with digits in the inner sectors
  Int_t fRowsWithDigits;       // number of rows with digits in the outer sectors
  Int_t fRowsTrackLength;      // last - first row with digit
  Float_t fPrim;               // theoretical dedx in tpc according particle momenta and mass
  digitRow fTPCRow;               // information about digits row pattern
  Int_t fNTPCRef;                 // tpc references counter
  Int_t fNITSRef;                 // ITS references counter
  Int_t fNTRDRef;                  // TRD references counter
  Int_t fNTOFRef;                  // TOF references counter
  TClonesArray * fTPCReferences;  //containner with all track references -in the TPC
  TClonesArray * fITSReferences;  //container with ITS references
  TClonesArray * fTRDReferences;  //container with TRD references  
  TClonesArray * fTOFReferences;  //container with TRD references  
  //
  ClassDef(AliMCInfo,3)  // container for 
};
ClassImp(AliMCInfo)


class AliGenV0Info: public TObject {
public:
  AliMCInfo fMCd;      //info about daughter particle - second particle for V0
  AliMCInfo fMCm;      //info about mother particle   - first particle for V0
  TParticle fMotherP;   //particle info about mother particle
  void Update(Float_t vertex[3]);        // put some derived info to special field 
  Double_t    fMCDist1;    //info about closest distance according closest MC - linear DCA
  Double_t    fMCDist2;    //info about closest distance parabolic DCA
  //
  Double_t     fMCPdr[3];    //momentum at vertex daughter  - according approx at DCA
  Double_t     fMCPd[4];     //exact momentum from MC info
  Double_t     fMCX[3];      //exact position of the vertex
  Double_t     fMCXr[3];     //rec. position according helix
  //
  Double_t     fMCPm[3];    //momentum at the vertex mother
  Double_t     fMCAngle[3]; //three angels
  Double_t     fMCRr;       // rec position of the vertex 
  Double_t     fMCR;        //exact r position of the vertex
  Int_t        fPdg[2];   //pdg code of mother and daugter particles
  Int_t        fLab[2];   //MC label of the partecle  
  //
  Double_t       fInvMass;  //reconstructed invariant mass -
  Float_t        fPointAngleFi; //point angle fi
  Float_t        fPointAngleTh; //point angle theta
  Float_t        fPointAngle;   //point angle full
  //
  ClassDef(AliGenV0Info,1)  // container for  
};
ClassImp(AliGenV0Info)


class AliGenKinkInfo: public TObject {
public:
  AliMCInfo fMCd;      //info about daughter particle - second particle for V0
  AliMCInfo fMCm;      //info about mother particle   - first particle for V0
  void Update();        // put some derived info to special field 
  Float_t GetQt();      //
  Double_t    fMCDist1;    //info about closest distance according closest MC - linear DCA
  Double_t    fMCDist2;    //info about closest distance parabolic DCA
  //
  Double_t     fMCPdr[3];    //momentum at vertex daughter  - according approx at DCA
  Double_t     fMCPd[4];     //exact momentum from MC info
  Double_t     fMCX[3];      //exact position of the vertex
  Double_t     fMCXr[3];     //rec. position according helix
  //
  Double_t     fMCPm[3];    //momentum at the vertex mother
  Double_t     fMCAngle[3]; //three angels
  Double_t     fMCRr;       // rec position of the vertex 
  Double_t     fMCR;        //exact r position of the vertex
  Int_t        fPdg[2];   //pdg code of mother and daugter particles
  Int_t        fLab[2];   //MC label of the partecle
  ClassDef(AliGenKinkInfo,1)  // container for  
};
ClassImp(AliGenKinkInfo)


////////////////////////////////////////////////////////////////////////
// 
// Start of implementation of the class AliGenInfoMaker
//
////////////////////////////////////////////////////////////////////////

class AliGenInfoMaker {

public:
  AliGenInfoMaker();
  AliGenInfoMaker(const char * fnGalice, const char* fnRes    ="genTracks.root",
		   Int_t nEvents=1, Int_t firstEvent=0);
  virtual ~AliGenInfoMaker();
  void Reset();
  Int_t Exec();
  Int_t Exec(Int_t nEvents, Int_t firstEventNr);
  void CreateTreeGenTracks();
  void CloseOutputFile();
  Int_t TreeKLoop();
  Int_t TreeTRLoop();
  Int_t TreeDLoop();
  Int_t BuildKinkInfo();  // build information about MC kinks
  Int_t BuildV0Info();  // build information about MC kinks
  Int_t BuildHitLines();  // build information about MC kinks
  //void  FillInfo(Int_t iParticle);
  void SetFirstEventNr(Int_t i) {fFirstEventNr = i;}
  void SetNEvents(Int_t i) {fNEvents = i;}
  void SetDebug(Int_t level) {fDebug = level;}
  Int_t SetIO(Int_t eventNr);
  Int_t CloseIOEvent();
  Int_t CloseIO();
  Int_t SetIO();
  Float_t TR2LocalX(AliTrackReference *trackRef,
		    AliTPCParam *paramTPC);
  AliMCInfo * GetInfo(UInt_t i){return (i<fNParticles)? fGenInfo[i]:0;}
  AliMCInfo * MakeInfo(UInt_t i);

public:
  Int_t fDebug;                   //! debug flag  
  Int_t fEventNr;                 //! current event number
  Int_t fLabel;                   //! track label
  Int_t fNEvents;                 //! number of events to process
  Int_t fFirstEventNr;            //! first event to process
  UInt_t fNParticles;              //! number of particles in TreeK
  TTree * fTreeGenTracks;          //! output tree with generated tracks
  TTree * fTreeKinks;             //!  output tree with Kinks
  TTree * fTreeV0;                //!  output tree with V0
  TTree * fTreeHitLines;          //! tree with hit lines
  char  fFnRes[1000];             //! output file name with stored tracks
  TFile *fFileGenTracks;          //! output file with stored fTreeGenTracks
  //
  AliRunLoader * fLoader;         //! pointer to the run loader
  TTree * fTreeD;                 //! current tree with digits
  TTree * fTreeTR;                //! current tree with TR
  AliStack *fStack;               //! current stack
  // 
  AliMCInfo **   fGenInfo;    //! array with pointers to gen info
  Int_t fNInfos;                  //! number of tracks with infos
  //
  AliTPCParam* fParamTPC;         //! AliTPCParam
  Float_t fVPrim[3];             //! primary vertex position
                                  // the fVDist[3] contains size of the 3-vector

private:

// some constants for the original non-pareller tracking (by Y.Belikov)
  static const Int_t seedRow11 = 158;  // nRowUp - 1
  static const Int_t seedRow12 = 139;  // nRowUp - 1 - (Int_t) 0.125*nRowUp
  static const Int_t seedRow21 = 149;  // seedRow11 - shift
  static const Int_t seedRow22 = 130;  // seedRow12 - shift
  static const Double_t kRaddeg = 180./3.14159265358979312;
  // 
  static const Double_t fgTPCPtCut = 0.03; // do not store particles with generated pT less than this
  static const Double_t fgITSPtCut = 0.2; // do not store particles with generated pT less than this
  static const Double_t fgTRDPtCut = 0.2; // do not store particles with generated pT less than this
  static const Double_t fgTOFPtCut = 0.15; // do not store particles with generated pT less than this
 
  ClassDef(AliGenInfoMaker,1)    // class which creates and fills tree with TPCGenTrack objects
};
ClassImp(AliGenInfoMaker)


class AliComparisonDraw: public TObject{
public:
  AliComparisonDraw(){fPoints=0; fView=0;}
  void InitView();
  TH1F * DrawXY(const char * chx, const char *chy, const char* selection, 
		const char * quality,Int_t nbins, Float_t minx, Float_t maxx, 
		Float_t miny, Float_t maxy, Int_t nBinsRes=30);
  TH1F * DrawLogXY(const char * chx, const char *chy, const char* selection, 
		   const char * quality, Int_t nbins,Float_t minx, Float_t maxx, 
		   Float_t miny, Float_t maxy, Int_t nBinsRes=30); 
  TH1F * Eff(const char *variable, const char* selection, const char * quality, 
	     Int_t nbins,Float_t min, Float_t max); 
  TH1F * EffLog(const char *variable, const char* selection, const char * quality, 
	     Int_t nbins,Float_t min, Float_t max);
  //
  static void   AliLabelAxes(TH1* histo, const char* xAxisTitle, const char* yAxisTitle);
  static Double_t* CreateLogBins(Int_t nBins, Double_t xMin, Double_t xMax);
  //
  static TH1F*  CreateEffHisto(TH1F* hGen, TH1F* hRec);
  static TH1F*  CreateResHisto(TH2F* hRes2, TH1F **phMean, 
				Bool_t drawBinFits = kTRUE,Bool_t overflowBinFits = kFALSE);
  void   DrawFriend2D(const char * chx, const char *chy, const char* selection, TTree * tfriend);
  void   GetPoints3D(const char * label, const char * chpoints, const char* selection, TTree * tpoints, Int_t color=6, Float_t rmin=4.);
  void   Draw3D(Int_t min=0, Int_t max = 10000);
  void SavePoints(const char* name);
 public: 
  TTree * fTree;
  TH1F  * fRes;  //temporary file
  TH1F  * fMean;  //temporary file
  TView * fView;  //3D view
  TCanvas *fCanvas; //canvas
  TObjArray *fPoints;
  ClassDef(AliComparisonDraw,1)
};
ClassImp(AliComparisonDraw)


class AliPointsMI: public TObject{
 public:
  AliPointsMI();
  AliPointsMI(Int_t n, Float_t *x,Float_t *y, Float_t *z);
  void Reset();
  void Reset(AliDetector * det, Int_t particle);  //load points for given particle
  ~AliPointsMI();
  Int_t   fN;  //number of points;
  Float_t *fX; //[fN] pointer to x
  Float_t *fY; //[fN] pointer to y
  Float_t *fZ; //[fN] pointer to Z  
  Int_t   fCapacity; //!allocated size of the x,y,x
  Int_t   fLabel0; //label
  Int_t   fLabel1; //label
  ClassDef(AliPointsMI,1)
};
ClassImp(AliPointsMI)


AliTPCParam * GetTPCParam();
Float_t TPCBetheBloch(Float_t bg);
Commit	Line	Data
ae7d73d2	1
	2
	3
	4	////////////////////////////////////////////////////////////////////////
	5	//
	6	// Start of implementation of the class digitRow
	7	//
	8	////////////////////////////////////////////////////////////////////////
	9	const Int_t kgRowBytes = 32;
	10
	11	class digitRow: public TObject {
	12
	13	public:
	14	digitRow();
	15	virtual ~digitRow(){;}
	16	void SetRow(Int_t row);
	17	Bool_t TestRow(Int_t row);
	18	digitRow & operator=(const digitRow &digOld);
	19	Int_t RowsOn(Int_t upto=8*kgRowBytes);
	20	Int_t Last();
	21	Int_t First();
	22	void Reset();
	23
	24	//private:
	25	UChar_t fDig[kgRowBytes];
	26	ClassDef(digitRow,1) // container for digit pattern
	27	};
	28	ClassImp(digitRow)
	29
	30
	31	////////////////////////////////////////////////////////////////////////
	32	//
	33	// Start of implementation of the class AliMCInfo
	34	//
	35	////////////////////////////////////////////////////////////////////////
	36
	37	class AliMCInfo: public TObject {
	38
	39	public:
	40	AliMCInfo();
	41	~AliMCInfo();
	42	void Update();
	43
	44	AliTrackReference fTrackRef; // track reference saved in the output tree
	45	AliTrackReference fTrackRefOut; // decay track reference saved in the output tree
	46	AliTrackReference fTRdecay; // track reference at decay point
	47	//
51ad6848	48	Int_t fPrimPart; // index of primary particle in TreeH
ae7d73d2	49	TParticle fParticle; // generated particle
	50	Float_t fMass; // mass of the particle
	51	Float_t fCharge; //
	52	Int_t fLabel; // track label
	53	Int_t fEventNr; // event number
	54	Int_t fMCtracks; // indication of how many times the track is retuturned back
	55	Int_t fPdg; //pdg code
	56	Float_t fDecayCoord[3]; // position of particle decay
	57	Double_t fVDist[4]; //distance of the particle vertex from primary vertex
	58	Bool_t fTPCdecay; //indicates decay in TPC
	59	Int_t fRowsWithDigitsInn; // number of rows with digits in the inner sectors
	60	Int_t fRowsWithDigits; // number of rows with digits in the outer sectors
	61	Int_t fRowsTrackLength; // last - first row with digit
	62	Float_t fPrim; // theoretical dedx in tpc according particle momenta and mass
	63	digitRow fTPCRow; // information about digits row pattern
	64	Int_t fNTPCRef; // tpc references counter
51ad6848	65	Int_t fNITSRef; // ITS references counter
	66	Int_t fNTRDRef; // TRD references counter
	67	Int_t fNTOFRef; // TOF references counter
ae7d73d2	68	TClonesArray * fTPCReferences; //containner with all track references -in the TPC
	69	TClonesArray * fITSReferences; //container with ITS references
	70	TClonesArray * fTRDReferences; //container with TRD references
51ad6848	71	TClonesArray * fTOFReferences; //container with TRD references
ae7d73d2	72	//
	73	ClassDef(AliMCInfo,3) // container for
	74	};
	75	ClassImp(AliMCInfo)
	76
	77
	78
	79	class AliGenV0Info: public TObject {
	80	public:
51ad6848	81	AliMCInfo fMCd; //info about daughter particle - second particle for V0
51ad6848	82	AliMCInfo fMCm; //info about mother particle - first particle for V0
81e97e0d	83	TParticle fMotherP; //particle info about mother particle
51ad6848	84	void Update(Float_t vertex[3]); // put some derived info to special field
ae7d73d2	85	Double_t fMCDist1; //info about closest distance according closest MC - linear DCA
	86	Double_t fMCDist2; //info about closest distance parabolic DCA
	87	//
	88	Double_t fMCPdr[3]; //momentum at vertex daughter - according approx at DCA
	89	Double_t fMCPd[4]; //exact momentum from MC info
	90	Double_t fMCX[3]; //exact position of the vertex
	91	Double_t fMCXr[3]; //rec. position according helix
	92	//
	93	Double_t fMCPm[3]; //momentum at the vertex mother
	94	Double_t fMCAngle[3]; //three angels
	95	Double_t fMCRr; // rec position of the vertex
	96	Double_t fMCR; //exact r position of the vertex
	97	Int_t fPdg[2]; //pdg code of mother and daugter particles
d6355c8c	98	Int_t fLab[2]; //MC label of the partecle
51ad6848	99	//
	100	Double_t fInvMass; //reconstructed invariant mass -
	101	Float_t fPointAngleFi; //point angle fi
	102	Float_t fPointAngleTh; //point angle theta
	103	Float_t fPointAngle; //point angle full
	104	//
ae7d73d2	105	ClassDef(AliGenV0Info,1) // container for
	106	};
	107	ClassImp(AliGenV0Info)
	108
	109
	110
	111
51ad6848	112	class AliGenKinkInfo: public TObject {
	113	public:
	114	AliMCInfo fMCd; //info about daughter particle - second particle for V0
	115	AliMCInfo fMCm; //info about mother particle - first particle for V0
	116	void Update(); // put some derived info to special field
	117	Float_t GetQt(); //
	118	Double_t fMCDist1; //info about closest distance according closest MC - linear DCA
	119	Double_t fMCDist2; //info about closest distance parabolic DCA
	120	//
	121	Double_t fMCPdr[3]; //momentum at vertex daughter - according approx at DCA
	122	Double_t fMCPd[4]; //exact momentum from MC info
	123	Double_t fMCX[3]; //exact position of the vertex
	124	Double_t fMCXr[3]; //rec. position according helix
	125	//
	126	Double_t fMCPm[3]; //momentum at the vertex mother
	127	Double_t fMCAngle[3]; //three angels
	128	Double_t fMCRr; // rec position of the vertex
	129	Double_t fMCR; //exact r position of the vertex
	130	Int_t fPdg[2]; //pdg code of mother and daugter particles
	131	Int_t fLab[2]; //MC label of the partecle
	132	ClassDef(AliGenKinkInfo,1) // container for
	133	};
	134	ClassImp(AliGenKinkInfo)
	135
	136
	137
	138
ae7d73d2	139	////////////////////////////////////////////////////////////////////////
	140	//
	141	// Start of implementation of the class AliGenInfoMaker
	142	//
	143	////////////////////////////////////////////////////////////////////////
	144
	145	class AliGenInfoMaker {
	146
	147	public:
	148	AliGenInfoMaker();
	149	AliGenInfoMaker(const char * fnGalice, const char* fnRes ="genTracks.root",
	150	Int_t nEvents=1, Int_t firstEvent=0);
	151	virtual ~AliGenInfoMaker();
	152	void Reset();
	153	Int_t Exec();
	154	Int_t Exec(Int_t nEvents, Int_t firstEventNr);
	155	void CreateTreeGenTracks();
	156	void CloseOutputFile();
	157	Int_t TreeKLoop();
	158	Int_t TreeTRLoop();
	159	Int_t TreeDLoop();
51ad6848	160	Int_t BuildKinkInfo(); // build information about MC kinks
	161	Int_t BuildV0Info(); // build information about MC kinks
	162	Int_t BuildHitLines(); // build information about MC kinks
ae7d73d2	163	//void FillInfo(Int_t iParticle);
	164	void SetFirstEventNr(Int_t i) {fFirstEventNr = i;}
	165	void SetNEvents(Int_t i) {fNEvents = i;}
	166	void SetDebug(Int_t level) {fDebug = level;}
	167	Int_t SetIO(Int_t eventNr);
	168	Int_t CloseIOEvent();
	169	Int_t CloseIO();
	170	Int_t SetIO();
	171	Float_t TR2LocalX(AliTrackReference *trackRef,
	172	AliTPCParam *paramTPC);
	173	AliMCInfo * GetInfo(UInt_t i){return (i<fNParticles)? fGenInfo[i]:0;}
	174	AliMCInfo * MakeInfo(UInt_t i);
	175
	176	public:
	177	Int_t fDebug; //! debug flag
	178	Int_t fEventNr; //! current event number
	179	Int_t fLabel; //! track label
	180	Int_t fNEvents; //! number of events to process
	181	Int_t fFirstEventNr; //! first event to process
	182	UInt_t fNParticles; //! number of particles in TreeK
51ad6848	183	TTree * fTreeGenTracks; //! output tree with generated tracks
	184	TTree * fTreeKinks; //! output tree with Kinks
	185	TTree * fTreeV0; //! output tree with V0
	186	TTree * fTreeHitLines; //! tree with hit lines
	187	char fFnRes[1000]; //! output file name with stored tracks
ae7d73d2	188	TFile *fFileGenTracks; //! output file with stored fTreeGenTracks
	189	//
	190	AliRunLoader * fLoader; //! pointer to the run loader
	191	TTree * fTreeD; //! current tree with digits
	192	TTree * fTreeTR; //! current tree with TR
	193	AliStack *fStack; //! current stack
	194	//
	195	AliMCInfo ** fGenInfo; //! array with pointers to gen info
	196	Int_t fNInfos; //! number of tracks with infos
	197	//
	198	AliTPCParam* fParamTPC; //! AliTPCParam
51ad6848	199	Float_t fVPrim[3]; //! primary vertex position
ae7d73d2	200	// the fVDist[3] contains size of the 3-vector
	201
	202	private:
	203
	204	// some constants for the original non-pareller tracking (by Y.Belikov)
	205	static const Int_t seedRow11 = 158; // nRowUp - 1
	206	static const Int_t seedRow12 = 139; // nRowUp - 1 - (Int_t) 0.125*nRowUp
	207	static const Int_t seedRow21 = 149; // seedRow11 - shift
	208	static const Int_t seedRow22 = 130; // seedRow12 - shift
	209	static const Double_t kRaddeg = 180./3.14159265358979312;
	210	//
51ad6848	211	static const Double_t fgTPCPtCut = 0.03; // do not store particles with generated pT less than this
ae7d73d2	212	static const Double_t fgITSPtCut = 0.2; // do not store particles with generated pT less than this
ae7d73d2	213	static const Double_t fgTRDPtCut = 0.2; // do not store particles with generated pT less than this
51ad6848	214	static const Double_t fgTOFPtCut = 0.15; // do not store particles with generated pT less than this
ae7d73d2	215
	216	ClassDef(AliGenInfoMaker,1) // class which creates and fills tree with TPCGenTrack objects
	217	};
	218	ClassImp(AliGenInfoMaker)
	219
	220
	221
	222	class AliComparisonDraw: public TObject{
	223	public:
51ad6848	224	AliComparisonDraw(){fPoints=0; fView=0;}
51ad6848	225	void InitView();
ae7d73d2	226	TH1F * DrawXY(const char * chx, const char chy, const char selection,
ae7d73d2	227	const char * quality,Int_t nbins, Float_t minx, Float_t maxx,
f007cb5f	228	Float_t miny, Float_t maxy, Int_t nBinsRes=30);
ae7d73d2	229	TH1F * DrawLogXY(const char * chx, const char chy, const char selection,
ae7d73d2	230	const char * quality, Int_t nbins,Float_t minx, Float_t maxx,
d6355c8c	231	Float_t miny, Float_t maxy, Int_t nBinsRes=30);
ae7d73d2	232	TH1F * Eff(const char variable, const char selection, const char * quality,
	233	Int_t nbins,Float_t min, Float_t max);
	234	TH1F * EffLog(const char variable, const char selection, const char * quality,
	235	Int_t nbins,Float_t min, Float_t max);
	236	//
	237	static void AliLabelAxes(TH1* histo, const char* xAxisTitle, const char* yAxisTitle);
	238	static Double_t* CreateLogBins(Int_t nBins, Double_t xMin, Double_t xMax);
	239	//
	240	static TH1F* CreateEffHisto(TH1F* hGen, TH1F* hRec);
	241	static TH1F* CreateResHisto(TH2F* hRes2, TH1F **phMean,
	242	Bool_t drawBinFits = kTRUE,Bool_t overflowBinFits = kFALSE);
51ad6848	243	void DrawFriend2D(const char * chx, const char chy, const char selection, TTree * tfriend);
	244	void GetPoints3D(const char * label, const char * chpoints, const char* selection, TTree * tpoints, Int_t color=6, Float_t rmin=4.);
	245	void Draw3D(Int_t min=0, Int_t max = 10000);
	246	void SavePoints(const char* name);
ae7d73d2	247	public:
	248	TTree * fTree;
	249	TH1F * fRes; //temporary file
	250	TH1F * fMean; //temporary file
51ad6848	251	TView * fView; //3D view
	252	TCanvas *fCanvas; //canvas
	253	TObjArray *fPoints;
ae7d73d2	254	ClassDef(AliComparisonDraw,1)
	255	};
	256	ClassImp(AliComparisonDraw)
	257
	258
51ad6848	259	class AliPointsMI: public TObject{
	260	public:
	261	AliPointsMI();
	262	AliPointsMI(Int_t n, Float_t x,Float_t y, Float_t *z);
	263	void Reset();
	264	void Reset(AliDetector * det, Int_t particle); //load points for given particle
	265	~AliPointsMI();
	266	Int_t fN; //number of points;
	267	Float_t *fX; //[fN] pointer to x
	268	Float_t *fY; //[fN] pointer to y
	269	Float_t *fZ; //[fN] pointer to Z
	270	Int_t fCapacity; //!allocated size of the x,y,x
	271	Int_t fLabel0; //label
	272	Int_t fLabel1; //label
	273	ClassDef(AliPointsMI,1)
	274	};
	275	ClassImp(AliPointsMI)
	276
ae7d73d2	277
	278	AliTPCParam * GetTPCParam();
	279	Float_t TPCBetheBloch(Float_t bg);