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

/* $Id $ */


#include "AliGenerator.h"
#include <TString.h>

// Read background particles from a boundary source
// Very specialized generator to simulate background from beam halo.
// Author: andreas.morsch@cern.ch

class AliGenHaloProtvino : public AliGenerator
{
public:
    enum constants{kY1Day0, kY1Day70, kY2D0, kY2D10, kY3D90};
    AliGenHaloProtvino();
    AliGenHaloProtvino(Int_t npart);
    AliGenHaloProtvino(const AliGenHaloProtvino &HaloProtvino);
    virtual ~AliGenHaloProtvino();
    virtual void Init();
    virtual void SetFileName(TString filename) {fFileName=TString(filename);}
    virtual void Generate();
    virtual Float_t GassPressureWeight(Float_t zPrimary);
    virtual void SetSide(Int_t flag = 1) {fSide = flag;}
    virtual void SetNskip(Int_t nskip) {fNskip = nskip;}
    virtual void SetRunPeriod(Int_t t = kY3D90) {fRunPeriod = t;}
    virtual void SetTimePerEvent(Float_t t = 1.e-4) {fTimePerEvent = t;}
    

    AliGenHaloProtvino & operator=(const AliGenHaloProtvino & rhs);

protected:
  FILE*    fFile;                       // ! Pointer to file
  TString  fFileName;                   //   Choose the file
  Int_t    fSide;                       //   Muon arm side (1) / Castor side (-1)
  Int_t    fRunPeriod;                  //   LHC Running Period
  Float_t  fTimePerEvent;               //   Time corresponding to one event [s]
  Int_t    fNskip;                      //   Number of entries to skip
  Float_t  fZ1[21],    fZ2[21];         // ! z-positions for gas pressure tables
  Float_t  fG1[21][5], fG2[21][5];      // ! gas pressures
 private:
  void Copy(TObject&) const;
  ClassDef(AliGenHaloProtvino,1)        //   LHC background boundary source (Protvino Group results)
      

};
#endif