[u/mrichter/AliRoot.git] / TPC / AliTPCcalibCosmic.h

#ifndef ALITPCCALIBCOSMIC_H
#define ALITPCCALIBCOSMIC_H

/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
 * See cxx source for full Copyright notice                               */

#include "AliTPCcalibBase.h"
#include "AliTPCCalPad.h"
#include "TH2F.h"
class TH1F;
class TList;
class AliESDEvent;
class AliESDtrack;

#include "TTreeStream.h"


class AliTPCcalibCosmic:public AliTPCcalibBase {
public:
  AliTPCcalibCosmic(); 
  AliTPCcalibCosmic(const Text_t *name, const Text_t *title);
  virtual ~AliTPCcalibCosmic();
  
  virtual void      Process(AliESDEvent *event);
  virtual Long64_t  Merge(TCollection *li);
  virtual void      Analyze();
  //
  void              FindPairs(AliESDEvent *event);
  Bool_t            IsPair(AliExternalTrackParam *tr0, AliExternalTrackParam *tr1);
  void              SetGainMap(AliTPCCalPad *GainMap){fGainMap = GainMap;};
  static void       CalculateBetheParams(TH2F *hist, Double_t * initialParam);
  static Double_t   CalculateMIPvalue(TH1F * hist);
  AliExternalTrackParam *Invert(AliExternalTrackParam *input);
  //
  TH1F   *          GetHistNTracks(){return fHistNTracks;};
  TH1F   *          GetHistClusters(){return fClusters;};
  TH2F   *          GetHistAcorde(){return fModules;};
  TH1F   *          GetHistPt(){return fHistPt;};
  TH2F   *          GetHistDeDx(){return fDeDx;};
  TH1F   *          GetHistMIP(){return fDeDxMIP;};
  //
  Double_t          GetMIPvalue(){return fMIPvalue;};
  //
  static void       BinLogX(TH1 * h);   // method for correct histogram binning


private:

  void              FillAcordeHist(AliESDtrack *upperTrack);

  AliTPCCalPad *fGainMap;         //  gain map from Krypton calibration
  TH1F  *fHistNTracks;            //  histogram showing number of ESD tracks per event
  TH1F  *fClusters;               //  histogram showing the number of clusters per track
  TH2F  *fModules;                //  2d histogram of tracks which are propagated to the ACORDE scintillator array
  TH1F  *fHistPt;                 //  Pt histogram of reconstructed tracks
  TH2F  *fDeDx;                   //  dEdx spectrum showing the different particle types
  TH1F  *fDeDxMIP;                //  TPC signal close to the MIP region of muons 0.4 < p < 0.45 GeV

  Double_t fMIPvalue;             //  MIP value calculated via a fit to fDeDxMIP

  // cuts
  //
  Float_t fCutMaxD;     // maximal distance in rfi ditection
  Float_t fCutTheta;    // maximal distance in theta ditection
  Float_t fCutMinDir;   // direction vector products

  AliTPCcalibCosmic(const AliTPCcalibCosmic&); 
  AliTPCcalibCosmic& operator=(const AliTPCcalibCosmic&); 

  ClassDef(AliTPCcalibCosmic, 1); 
};

#endif
Commit	Line	Data
f7f33dec	1	#ifndef ALITPCCALIBCOSMIC_H
	2	#define ALITPCCALIBCOSMIC_H
	3
	4	/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
	5	* See cxx source for full Copyright notice */
	6
	7	#include "AliTPCcalibBase.h"
d6920632	8	#include "AliTPCCalPad.h"
54b76c13	9	#include "TH2F.h"
f7f33dec	10	class TH1F;
f7f33dec	11	class TList;
f7f33dec	12	class AliESDEvent;
54b76c13	13	class AliESDtrack;
f7f33dec	14
	15	#include "TTreeStream.h"
	16
	17
	18	class AliTPCcalibCosmic:public AliTPCcalibBase {
	19	public:
	20	AliTPCcalibCosmic();
	21	AliTPCcalibCosmic(const Text_t name, const Text_t title);
	22	virtual ~AliTPCcalibCosmic();
	23
54b76c13	24	virtual void Process(AliESDEvent *event);
	25	virtual Long64_t Merge(TCollection *li);
	26	virtual void Analyze();
	27	//
	28	void FindPairs(AliESDEvent *event);
	29	Bool_t IsPair(AliExternalTrackParam tr0, AliExternalTrackParam tr1);
	30	void SetGainMap(AliTPCCalPad *GainMap){fGainMap = GainMap;};
	31	static void CalculateBetheParams(TH2F hist, Double_t initialParam);
	32	static Double_t CalculateMIPvalue(TH1F * hist);
428b47be	33	AliExternalTrackParam Invert(AliExternalTrackParam input);
54b76c13	34	//
	35	TH1F * GetHistNTracks(){return fHistNTracks;};
	36	TH1F * GetHistClusters(){return fClusters;};
	37	TH2F * GetHistAcorde(){return fModules;};
	38	TH1F * GetHistPt(){return fHistPt;};
	39	TH2F * GetHistDeDx(){return fDeDx;};
	40	TH1F * GetHistMIP(){return fDeDxMIP;};
	41	//
	42	Double_t GetMIPvalue(){return fMIPvalue;};
	43	//
	44	static void BinLogX(TH1 * h); // method for correct histogram binning
	45
	46
	47	private:
	48
	49	void FillAcordeHist(AliESDtrack *upperTrack);
	50
	51	AliTPCCalPad *fGainMap; // gain map from Krypton calibration
	52	TH1F *fHistNTracks; // histogram showing number of ESD tracks per event
	53	TH1F *fClusters; // histogram showing the number of clusters per track
	54	TH2F *fModules; // 2d histogram of tracks which are propagated to the ACORDE scintillator array
	55	TH1F *fHistPt; // Pt histogram of reconstructed tracks
	56	TH2F *fDeDx; // dEdx spectrum showing the different particle types
	57	TH1F *fDeDxMIP; // TPC signal close to the MIP region of muons 0.4 < p < 0.45 GeV
	58
	59	Double_t fMIPvalue; // MIP value calculated via a fit to fDeDxMIP
	60
9b27d39b	61	// cuts
	62	//
	63	Float_t fCutMaxD; // maximal distance in rfi ditection
	64	Float_t fCutTheta; // maximal distance in theta ditection
	65	Float_t fCutMinDir; // direction vector products
54b76c13	66
f7f33dec	67	AliTPCcalibCosmic(const AliTPCcalibCosmic&);
	68	AliTPCcalibCosmic& operator=(const AliTPCcalibCosmic&);
	69
	70	ClassDef(AliTPCcalibCosmic, 1);
	71	};
	72
	73	#endif
	74