PndFsmMvdPid.h

Go to the documentation of this file.

//--------------------------------------------------------------------------
// File and Version Information:
//      $Id: PndFsmMvdPid.hh,v 1.4 2006/08/10 09:58:41 klausg Exp $
//
// Description:
//      Class PndFsmMvdPid
//
//  Implementation of the MVD for the Fast Sim Detectors
//
//  This software was developed for the PANDA collaboration.  If you
//  use all or part of it, please give an appropriate acknowledgement.
//
// Author List:
//      Klaus Goetzen                    Original Author
//
// Copyright Information:
//      Copyright (C) 2006              GSI
//
//------------------------------------------------------------------------

#ifndef PNDFSMMVDPID_H
#define PNDFSMMVDPID_H

//----------------------
// Base Class Headers --
//----------------------

//-------------------------------
// Collaborating class Headers --
//-------------------------------
#include <iosfwd>
#include "ArgList.h"
#include "PndFsmAbsDet.h"

//--------------------------------------------
// Collaborating class forward declarations --
// -------------------------------------------
class PndFsmTrack;
class PndFsmResponse;

class PndFsmMvdPid : public PndFsmAbsDet {
 public:
  //--------------------
  // Public interface --
  //--------------------

  //
  // Constructors
  //

  PndFsmMvdPid();
  PndFsmMvdPid(ArgList &par);
  virtual ~PndFsmMvdPid();

  virtual PndFsmResponse *respond(PndFsmTrack *t);

 private:
  bool detected(PndFsmTrack *t) const;

  void initParameters();
  void print(std::ostream &o);
  bool setParameter(std::string &name, double value);

  double _thtMin;
  double _thtMax;
  double _ptmin;
  double _dEdxResMulti;

  enum { nPidType = 5 };
  enum PidType { electron, muon, pion, kaon, proton };

  double MeanEnergyLoss(PidType particle);
  double LandauGaus(double s_mpv, double width1, double width2);
  double Likelihood(PidType particle);
  double mpv(PidType particle);
  double width1(PidType particle);
  double width2(PidType particle);

  double _momentum;
  double _energyloss;
};

#endif