[u/mrichter/AliRoot.git] / PWGPP / TPC / AliPerformanceMC.h

#ifndef ALIPERFORMANCEMC_H
#define ALIPERFORMANCEMC_H

//------------------------------------------------------------------------------
// Class to keep information for MC  tracks 
// to check propagation algorithm, B-field and material budget.   
// 
// Author: J.Otwinowski 09/06/2009 
//------------------------------------------------------------------------------

class TString;
class TNamed;
class TCanvas;
class TH1F;
class TH2F;
class TParticle;

class AliESDVertex;
class AliESDtrack;
class AliMCEvent;
class AliStack;
class AliTrackReference;
class AliESDEvent; 
class AliESDfriend; 
class AliESDfriendTrack; 
class AliMCEvent;
class AliMCParticle;
class AliMCInfoCuts;
class AliRecInfoCuts;
class AliExternalTrackParam;

#include "THnSparse.h"
#include "AliPerformanceObject.h"

class AliPerformanceMC : public AliPerformanceObject {
public :
  AliPerformanceMC(const Char_t* name="AliPerformanceMC", const Char_t* title="AliPerformanceMC",Int_t analysisMode=0,Bool_t hptGenerator=kFALSE);
  virtual ~AliPerformanceMC();

  // Init data members
  virtual void  Init();

  // Execute analysis
  virtual void  Exec(AliMCEvent* const mcEvent, AliESDEvent *const esdEvent, AliESDfriend *const esdFriend, const Bool_t bUseMC, const Bool_t bUseESDfriend);

  // Merge output objects (needed by PROOF) 
  virtual Long64_t Merge(TCollection* const list);

  // Analyse output histograms
  virtual void Analyse();

  // Get analysis folder
  virtual TFolder* GetAnalysisFolder() const {return fAnalysisFolder;}

  // Process events
  void ProcessTPC(AliTrackReference* const refIn, TParticle* const particle);
  void ProcessInnerTPC(AliTrackReference* const refIn,  AliTrackReference* const refOut, TParticle* const particle);
  void ProcessOuterTPCExt(TParticle *const  part, TClonesArray * const trefs);
  AliExternalTrackParam * MakeTrack(const AliTrackReference* ref, TParticle *const part);

  // Create folder for analysed histograms
  TFolder *CreateFolder(TString folder = "folderMC", TString title = "Analysed MC histograms");

  // Export objects to folder
  TFolder *ExportToFolder(TObjArray * array=0);

  // Selection cuts
  void SetAliRecInfoCuts(AliRecInfoCuts* const cuts=0) {fCutsRC = cuts;}   
  void SetAliMCInfoCuts(AliMCInfoCuts* const cuts=0) {fCutsMC = cuts;}  
   
  AliRecInfoCuts*  GetAliRecInfoCuts() const {return fCutsRC;}  
  AliMCInfoCuts*   GetAliMCInfoCuts()  const {return fCutsMC;}  

  TH1F*  MakeResol(TH2F * his, Int_t integ=0, Bool_t type=kFALSE, Int_t cut=0); 

  // getters
  //
  THnSparse *GetResolHisto() const  { return fResolHisto; }
  THnSparse *GetPullHisto()  const  { return fPullHisto; }
 
private:
  //
  // Control histograms
  // 5 track parameters (details in STEER/AliExternalTrackParam.h)
  //

  // resolution histogram
  THnSparseF *fResolHisto; //-> res_y:res_z:res_phi:res_lambda:res_pt:y:z:phi:eta:pt

  // pull histogram
  THnSparseF *fPullHisto;  //-> pull_y:pull_z:pull_snp:pull_tgl:pull_1pt:y:z:snp:tgl:1pt

  // Global cuts objects
  AliRecInfoCuts*  fCutsRC;      // selection cuts for reconstructed tracks
  AliMCInfoCuts*  fCutsMC;       // selection cuts for MC tracks

  // analysis folder 
  TFolder *fAnalysisFolder; // folder for analysed histograms

  AliPerformanceMC(const AliPerformanceMC&); // not implemented
  AliPerformanceMC& operator=(const AliPerformanceMC&); // not implemented

  ClassDef(AliPerformanceMC,1);
};

#endif
Commit	Line	Data
7cc34f08	1	#ifndef ALIPERFORMANCEMC_H
	2	#define ALIPERFORMANCEMC_H
	3
	4	//------------------------------------------------------------------------------
	5	// Class to keep information for MC tracks
	6	// to check propagation algorithm, B-field and material budget.
	7	//
	8	// Author: J.Otwinowski 09/06/2009
	9	//------------------------------------------------------------------------------
	10
	11	class TString;
	12	class TNamed;
	13	class TCanvas;
	14	class TH1F;
	15	class TH2F;
	16	class TParticle;
	17
	18	class AliESDVertex;
	19	class AliESDtrack;
	20	class AliMCEvent;
	21	class AliStack;
	22	class AliTrackReference;
	23	class AliESDEvent;
	24	class AliESDfriend;
	25	class AliESDfriendTrack;
	26	class AliMCEvent;
	27	class AliMCParticle;
	28	class AliMCInfoCuts;
	29	class AliRecInfoCuts;
	30	class AliExternalTrackParam;
	31
	32	#include "THnSparse.h"
	33	#include "AliPerformanceObject.h"
	34
	35	class AliPerformanceMC : public AliPerformanceObject {
	36	public :
4823e4d4	37	AliPerformanceMC(const Char_t* name="AliPerformanceMC", const Char_t* title="AliPerformanceMC",Int_t analysisMode=0,Bool_t hptGenerator=kFALSE);
7cc34f08	38	virtual ~AliPerformanceMC();
	39
	40	// Init data members
	41	virtual void Init();
	42
	43	// Execute analysis
	44	virtual void Exec(AliMCEvent* const mcEvent, AliESDEvent const esdEvent, AliESDfriend const esdFriend, const Bool_t bUseMC, const Bool_t bUseESDfriend);
	45
	46	// Merge output objects (needed by PROOF)
	47	virtual Long64_t Merge(TCollection* const list);
	48
	49	// Analyse output histograms
	50	virtual void Analyse();
	51
	52	// Get analysis folder
	53	virtual TFolder* GetAnalysisFolder() const {return fAnalysisFolder;}
	54
	55	// Process events
	56	void ProcessTPC(AliTrackReference* const refIn, TParticle* const particle);
	57	void ProcessInnerTPC(AliTrackReference* const refIn, AliTrackReference* const refOut, TParticle* const particle);
	58	void ProcessOuterTPCExt(TParticle const part, TClonesArray const trefs);
	59	AliExternalTrackParam * MakeTrack(const AliTrackReference* ref, TParticle *const part);
	60
	61	// Create folder for analysed histograms
	62	TFolder *CreateFolder(TString folder = "folderMC", TString title = "Analysed MC histograms");
	63
	64	// Export objects to folder
	65	TFolder ExportToFolder(TObjArray array=0);
	66
	67	// Selection cuts
	68	void SetAliRecInfoCuts(AliRecInfoCuts* const cuts=0) {fCutsRC = cuts;}
	69	void SetAliMCInfoCuts(AliMCInfoCuts* const cuts=0) {fCutsMC = cuts;}
	70
	71	AliRecInfoCuts* GetAliRecInfoCuts() const {return fCutsRC;}
	72	AliMCInfoCuts* GetAliMCInfoCuts() const {return fCutsMC;}
	73
	74	TH1F* MakeResol(TH2F * his, Int_t integ=0, Bool_t type=kFALSE, Int_t cut=0);
	75
	76	// getters
	77	//
	78	THnSparse *GetResolHisto() const { return fResolHisto; }
	79	THnSparse *GetPullHisto() const { return fPullHisto; }
	80
	81	private:
	82	//
	83	// Control histograms
	84	// 5 track parameters (details in STEER/AliExternalTrackParam.h)
	85	//
	86
	87	// resolution histogram
	88	THnSparseF *fResolHisto; //-> res_y:res_z:res_phi:res_lambda:res_pt:y:z:phi:eta:pt
	89
	90	// pull histogram
	91	THnSparseF *fPullHisto; //-> pull_y:pull_z:pull_snp:pull_tgl:pull_1pt:y:z:snp:tgl:1pt
	92
	93	// Global cuts objects
	94	AliRecInfoCuts* fCutsRC; // selection cuts for reconstructed tracks
	95	AliMCInfoCuts* fCutsMC; // selection cuts for MC tracks
	96
	97	// analysis folder
	98	TFolder *fAnalysisFolder; // folder for analysed histograms
	99
	100	AliPerformanceMC(const AliPerformanceMC&); // not implemented
	101	AliPerformanceMC& operator=(const AliPerformanceMC&); // not implemented
102
103	ClassDef(AliPerformanceMC,1);
104	};
105
106	#endif