mRICH.h

#ifndef __MRICH_H__
#define __MRICH_H__
	
//
//  Hello mRICH Fans:
//
//  This is an example class that inherits from the PID base class.
//  We shall implement all the necessary functions as required by the
//  base class.
//
//  The UNIQUE identifiers for this class are radius, eta extent, and 
//  time resolution. 
//
//  Note that in keeping with any well written base class, we shall
//  give default values to all the constructor arguments.
//
//  This routine assumes units of cm for distances and picoSeconds for time.
//

#include "PID.h"
#include "TVector3.h"
#include "TRandom.h"
#include <iostream>
#include <TH1.h>
#include <TH2.h>

	
class mRICH: public PID
{
public:
  mRICH(double trackResolution=0.5, double timePrecision=1.0, double pixS=0.5, double p=5.0);
  virtual ~mRICH() {}
	
  //bool   valid   (double eta, double p) {return (((eta>-8&&eta<-1) || (eta>1&&eta<8)) && (p>pLow && p<pHigh));}
  bool   valid   (double eta, double p) {return (((eta>-8&&eta<8)) && (p>pLow && p<pHigh));}
  double maxP    (double eta, double p, double numSigma, PID::type PID);
  double minP    (double eta, double numSigma, PID::type PID) {return 0;}
  string name    () {return myName;}
  void   description ();
  void ReadMap(TString name);
  double getAng(double mass);
  double getdAng(double mass);
  double getNgamma(double t, double mass);
  double T_Aer(double t, double lam);
  double T_QE(double lam);

  double SetLensFocalLength(double f);
  double SetAerogelThickness(double t);
  double GetmRICHParameters(); 

 protected:
  std::string myName;

  double Nsigma;

  // Physical constants (should come from elsewhere!)
  double mPion;    // GeV/c^2
  double mKaon;    // GeV/c^2
  double mProton;  // GeV/c^2
  double mom;      // GeV/c
  double c;        // cm/picosecond;
  double n;
  double f;        //mm
  double a;        //mm
  double N_gam;
  double pi;
  double alpha;
  double L;
  double th0;

  double fTrackResolution;
  double fTimePrecision;

  double pLow;
  double pHigh;

  TH2F *fpixMap;

};
	
#endif /* __PID_H__ */