/*
 * TofDedxHeCut.h
 *
 *  Created on: 03/ago/2009
 *      Author: Nicola Mori
 */

/*! @file TofDedxHeCut.h  The TofDedxHeCut.h class definition file. */

#ifndef TOFDEDXHECUT_H_
#define TOFDEDXHECUT_H_

#include "../../PamCutBase/PamCutBase.h"

/*! @brief Cutting a band dE/dx vs. Beta for Helium.
 * This cut discards events where TOF dE/dx vs. Beta are outside the helium "corridor". This corridor is currently
 * defined by:
 *
 * dE/dx > 1 / (0.2 * (beta - 0.05)) - 2.
 * dE/dx < 1 / (0.25 * (beta - 0.1))
 * dE/dx > 11. * beta - 7.7
 *
 * dE/dx is the mean of the energy releases on S11, S12, S21, S22.
 */
class TofDedxHeCut: public PamCut {

public:
  /*! @brief Constructor.
   *
   * @param cutName The cut's name.
   * @param standAlone If true, the cut will use the standalone measure of beta; otherwise it will use the
   *                   physical track.
   */
  TofDedxHeCut(const char *cutName, bool standAlone = true) :
    PamCut(cutName), _standAlone(standAlone) {
  }
  /*! @brief Destructor. */
  ~TofDedxHeCut() {

  }

  /*! @brief The dE/dx vs. beta Helium check.
   *
   * @param event The event to analyze.
   * @return #CUTOK if event is inside helium "corridor".
   * @return 0 otherwise.
   */
  int Check(PamLevel2 *event);

private:

  bool _standAlone;
};

#endif /* TOFDEDXHECUT_H_ */