Updates for mixing.
[u/mrichter/AliRoot.git] / EVGEN / AliGenHalo.h
CommitLineData
f87cfe57 1#ifndef ALIGENHALO_H
2#define ALIGENHALO_H
3da30618 3/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * See cxx source for full Copyright notice */
5
6/* $Id$ */
7
f87cfe57 8
fe4da5cc 9#include "AliGenerator.h"
bd1cece1 10#include <TString.h>
675e9664 11
d12ebe74 12// Read beam halo background particles from a boundary source
13// Boundary source is in the LHCb format
14// http://www.hep.manchester.ac.uk/u/robert/LHC_backgrounds/Note-MIBStudies.pdf
15// and has been provided by Robert Appleby
675e9664 16// Author: andreas.morsch@cern.ch
fe4da5cc 17
18class AliGenHalo : public AliGenerator
19{
d12ebe74 20public:
21 enum constants{kY1Day0, kY1Day70, kY2D0, kY2D10, kY3D90, kLHCPR674Startup, kLHCPR674Conditioned};
fe4da5cc 22 AliGenHalo();
23 AliGenHalo(Int_t npart);
24 virtual ~AliGenHalo();
25 virtual void Init();
bd1cece1 26 virtual void SetFileName(TString filename) {fFileName=TString(filename);}
fe4da5cc 27 virtual void Generate();
d12ebe74 28 virtual Float_t GasPressureWeight(Float_t zPrimary);
29 virtual void SetSide(Int_t flag = 1) {fSide = flag;}
30 virtual void SetNskip(Int_t nskip) {fNskip = nskip;}
31 virtual void SetRunPeriod(Int_t t = kY3D90) {fRunPeriod = t;}
32 virtual void SetTimePerEvent(Float_t t = 1.e-4) {fTimePerEvent = t;}
33 virtual void Draw(Option_t * opt="");
73d273a4 34 virtual void CountEvents();
d12ebe74 35 private:
c321951c 36 virtual void SkipEvents();
d12ebe74 37 virtual Int_t ReadNextParticle();
198bb1c7 38 protected:
d12ebe74 39 FILE* fFile; // ! Pointer to file
198bb1c7 40 TString fFileName; // Choose the file
d12ebe74 41 Int_t fSide; // Muon arm side (1) / Castor side (-1)
42 Int_t fRunPeriod; // LHC Running Period
43 Float_t fTimePerEvent; // Time corresponding to one event [s]
44 Int_t fNskip; // Number of entries to skip
45 Float_t* fZ1; // ! z-positions for gas pressure tables
46 Float_t* fZ2; // ! z-positions for gas pressure tables
47 Float_t* fG1; // ! gas pressures
48 Float_t* fG2; // ! gas pressures
49 Int_t fGPASize; // ! Size of arrays
50 Int_t fLossID; // ! unique loss ID
51 Int_t fLossA; // ! atomic number of scatterer
52 Int_t fPdg; // ! pdg code
53 Float_t fLossT0; // ! relative time
54 Float_t fLossZ; // ! scaterring position
55 Float_t fLossW; // ! weight of proton loss
56 Float_t fXS; // ! x-position on scoring plane
57 Float_t fYS; // ! y-position on scoring plane
58 Float_t fZS; // ! z-position on scoring plane
59 Float_t fDX; // ! direction cosine x
60 Float_t fDY; // ! direction cosine y
61 Float_t fEkin; // ! kinetic energy
62 Float_t fTS; // ! relative arrival time
63 Float_t fWS; // ! weight
ac3faee4 64 private:
93a2041b 65 AliGenHalo(const AliGenHalo &Halo);
66 AliGenHalo & operator=(const AliGenHalo & rhs);
d12ebe74 67
68 ClassDef(AliGenHalo,1) // LHC beam halo boundary source
fe4da5cc 69};
70#endif
71
72
73
74
75
76