[u/mrichter/AliRoot.git] / PHOS / AliPHOSPIDv0.h

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

/* $Id$ */


//_________________________________________________________________________
// Implementation version v0 of the PHOS particle identifier 
// Identification is based on information from PPSD and EMC
// Oh yeah                 
//*-- Author: Yves Schutz (SUBATECH)

// --- ROOT system ---
class TFormula ;
class TVector3 ;

// --- Standard library ---

// --- AliRoot header files ---
class AliPHOSEmcRecPoint ;
class AliPHOSRecPoint ;

#include "AliPHOSPID.h"

class  AliPHOSPIDv0 : public AliPHOSPID {

public:

  AliPHOSPIDv0() ;          // ctor            
  AliPHOSPIDv0(const char* headerFile, const char * tsBranch = "Default", const Bool_t toSplit=kFALSE) ;
  AliPHOSPIDv0(AliPHOSPIDv0 & pid) ;          // cpy ctor            
  virtual ~AliPHOSPIDv0() ; // dtor

  virtual void Exec(Option_t * option);
  //  virtual char * GetRecParticlesBranch()const {return (char*) fRecParticlesTitle.Data() ;}      
  //  virtual char * GetTrackSegmentsBranch()const{return (char*) fTrackSegmentsTitle.Data(); }
  virtual const Int_t GetRecParticlesInRun() const  {return fRecParticlesInRun ;}  

  virtual void PlotDispersionCuts()const ;
  virtual void Print(Option_t * option)const ; 
  virtual void SetIdentificationMethod(char * option = "CPV DISP" ){fIDOptions = option ;} 
  virtual void SetShowerProfileCut(char * formula = "0.35*0.35 - (x-1.386)*(x-1.386) - 1.707*1.707*(y-1.008)*(y-1.008)") ;
  virtual void SetDispersionCut(Float_t cut){fDispersion = cut ; } 
  virtual void SetCpvtoEmcDistanceCut(Float_t cut )      {fCpvEmcDistance = cut ;}
  virtual void SetTimeGate(Float_t gate)                 {fTimeGate = gate ;}
  //  virtual void SetTrackSegmentsBranch(const char* title) { fTrackSegmentsTitle = title;}
  //  virtual void SetRecParticlesBranch (const char* title) { fRecParticlesTitle = title;} 
  virtual const char * Version() const { return "pid-v0" ; }  
  AliPHOSPIDv0 & operator = (const  AliPHOSPIDv0 & pid) { return *this ; }                 
 private:
  
  virtual void Init() ;
  void     MakeRecParticles(void ) ;
  Float_t  GetDistance(AliPHOSEmcRecPoint * emc, AliPHOSRecPoint * cpv, Option_t * Axis)const ; // Relative Distance CPV-EMC
  TVector3 GetMomentumDirection(AliPHOSEmcRecPoint * emc, AliPHOSRecPoint * cpv)const ;
  void     PrintRecParticles(Option_t * option) ;
  virtual void WriteRecParticles(Int_t event) ; 

 private:

  TString                fHeaderFileName ;    // file name with event header
  TString                fTrackSegmentsTitle; // branch name with track segments
  TString                fRecPointsTitle ;    // branch name with rec points
  TString                fRecParticlesTitle ; // branch name with rec particles
  TString                fIDOptions ;         // PID option
  Int_t                  fNEvent ;            // current event number

  AliPHOSClusterizer   * fClusterizer ;       // !
  AliPHOSTrackSegmentMaker * fTSMaker ;       // !

  TFormula             * fFormula ;           // formula to define cut on the shower elips axis
  Float_t                fDispersion ;        // dispersion cut
  Float_t                fCpvEmcDistance ;    // Max EMC-CPV distance
  Float_t                fTimeGate ;          // Time of the latest EmcRecPoint accepted as EM
  Int_t                  fRecParticlesInRun ; //! Total number of recparticles in one run

  ClassDef( AliPHOSPIDv0,1)  // Particle identifier implementation version 1

};

#endif // AliPHOSPIDV0_H
Commit	Line	Data
b91d88dc	1	#ifndef ALIPHOSPIDV0_H
	2	#define ALIPHOSPIDV0_H
	3	/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
	4	* See cxx source for full Copyright notice */
	5
	6	/* $Id$ */
	7
	8
	9	//_________________________________________________________________________
	10	// Implementation version v0 of the PHOS particle identifier
	11	// Identification is based on information from PPSD and EMC
	12	// Oh yeah
	13	//*-- Author: Yves Schutz (SUBATECH)
	14
	15	// --- ROOT system ---
	16	class TFormula ;
	17	class TVector3 ;
	18
	19	// --- Standard library ---
	20
	21	// --- AliRoot header files ---
	22	class AliPHOSEmcRecPoint ;
	23	class AliPHOSRecPoint ;
	24
	25	#include "AliPHOSPID.h"
	26
	27	class AliPHOSPIDv0 : public AliPHOSPID {
	28
	29	public:
	30
	31	AliPHOSPIDv0() ; // ctor
fbf811ec	32	AliPHOSPIDv0(const char* headerFile, const char * tsBranch = "Default", const Bool_t toSplit=kFALSE) ;
531f44b7	33	AliPHOSPIDv0(AliPHOSPIDv0 & pid) ; // cpy ctor
b91d88dc	34	virtual ~AliPHOSPIDv0() ; // dtor
	35
	36	virtual void Exec(Option_t * option);
fbf811ec	37	// virtual char * GetRecParticlesBranch()const {return (char*) fRecParticlesTitle.Data() ;}
fbf811ec	38	// virtual char * GetTrackSegmentsBranch()const{return (char*) fTrackSegmentsTitle.Data(); }
b91d88dc	39	virtual const Int_t GetRecParticlesInRun() const {return fRecParticlesInRun ;}
	40
	41	virtual void PlotDispersionCuts()const ;
	42	virtual void Print(Option_t * option)const ;
	43	virtual void SetIdentificationMethod(char * option = "CPV DISP" ){fIDOptions = option ;}
	44	virtual void SetShowerProfileCut(char * formula = "0.350.35 - (x-1.386)(x-1.386) - 1.7071.707(y-1.008)*(y-1.008)") ;
	45	virtual void SetDispersionCut(Float_t cut){fDispersion = cut ; }
	46	virtual void SetCpvtoEmcDistanceCut(Float_t cut ) {fCpvEmcDistance = cut ;}
	47	virtual void SetTimeGate(Float_t gate) {fTimeGate = gate ;}
fbf811ec	48	// virtual void SetTrackSegmentsBranch(const char* title) { fTrackSegmentsTitle = title;}
fbf811ec	49	// virtual void SetRecParticlesBranch (const char* title) { fRecParticlesTitle = title;}
b91d88dc	50	virtual const char * Version() const { return "pid-v0" ; }
531f44b7	51	AliPHOSPIDv0 & operator = (const AliPHOSPIDv0 & pid) { return *this ; }
b91d88dc	52	private:
	53
	54	virtual void Init() ;
	55	void MakeRecParticles(void ) ;
	56	Float_t GetDistance(AliPHOSEmcRecPoint * emc, AliPHOSRecPoint * cpv, Option_t * Axis)const ; // Relative Distance CPV-EMC
	57	TVector3 GetMomentumDirection(AliPHOSEmcRecPoint * emc, AliPHOSRecPoint * cpv)const ;
	58	void PrintRecParticles(Option_t * option) ;
	59	virtual void WriteRecParticles(Int_t event) ;
	60
	61	private:
	62
	63	TString fHeaderFileName ; // file name with event header
	64	TString fTrackSegmentsTitle; // branch name with track segments
	65	TString fRecPointsTitle ; // branch name with rec points
	66	TString fRecParticlesTitle ; // branch name with rec particles
	67	TString fIDOptions ; // PID option
	68	Int_t fNEvent ; // current event number
	69
	70	AliPHOSClusterizer * fClusterizer ; // !
	71	AliPHOSTrackSegmentMaker * fTSMaker ; // !
	72
fbf811ec	73	TFormula * fFormula ; // formula to define cut on the shower elips axis
b91d88dc	74	Float_t fDispersion ; // dispersion cut
	75	Float_t fCpvEmcDistance ; // Max EMC-CPV distance
	76	Float_t fTimeGate ; // Time of the latest EmcRecPoint accepted as EM
	77	Int_t fRecParticlesInRun ; //! Total number of recparticles in one run
	78
	79	ClassDef( AliPHOSPIDv0,1) // Particle identifier implementation version 1
	80
	81	};
	82
	83	#endif // AliPHOSPIDV0_H