STEER/AliMagFC.h

   1 #ifndef ALIMAGFC_H
   2 #define ALIMAGFC_H
   3 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
   4 *See cxx source for full Copyright notice                               */
   5
   6 /* $Id$ */
   7
   8 //-------------------------------------------------------------------------
   9 //     Constant magnetic field class
  10 //     Used by AliRun class
  11 //     Author:
  12 //-------------------------------------------------------------------------
  13
  14 #include "AliMagF.h"
  15
  16 enum BeamType_t {kBeamTypeAA, kBeamTypepp};
  17 class AliMagFC  : public AliMagF
  18 {
  19   //Alice Constant Magnetic Field
  20
  21 public:
  22   AliMagFC();
  23   AliMagFC(const char *name, const char *title, Int_t integ,
  24            Float_t factor, Float_t fmax);
  25   virtual ~AliMagFC(){}
  26   virtual void Field(Float_t *x, Float_t *b) const;
  27   virtual void ReadField() {}
  28   virtual void ZDCField(Float_t *x, Float_t *b) const;
  29   virtual void SetBeamType(BeamType_t type)      {fBeamType    = type;}
  30   virtual void SetBeamEnergy(Float_t energy)     {fBeamEnergy  = energy;}
  31   virtual void SetCompensatorMagnet(Bool_t flag) {fCompensator = flag;}
  32
  33 private:
  34   Bool_t     fCompensator; // Flag for compensator magnetic field (kTrue -> ON)
  35   BeamType_t fBeamType;    // Beam type: A-A (fBeamType=0) or p-p (fBeamType=1)
  36   Float_t    fBeamEnergy;  // Beam energy in GeV
  37   Float_t    fQuadGradient;// Gradient field for inner triplet quadrupoles
  38   Float_t    fDipoleField; // Field value for D1 and D2 dipoles
  39   Float_t    fCCorrField;  // Side C 2nd compensator field
  40   Float_t    fACorr1Field; // Side A 1st compensator field
  41   Float_t    fACorr2Field; // Side A 2nd compensator field
  42
  43   ClassDef(AliMagFC,3)  //Class for all Alice Constant MagField
  44 };
  45
  46
  47 // %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%  ZDC part  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  48
  49 // ************************ LHC optics v6.5 *****************************
  50 // ---- Position of the magnetic elements of LHC beam optics ----
  51 // -> SIDE C
  52 static const Float_t kCCorrBegin = -1972.5, kCCorrEnd = kCCorrBegin - 153., kCCorrSqRadius = 4.5*4.5;
  53 //
  54 static const Float_t kCTripletBegin  = -2296.5;
  55 static const Float_t kCQ1Begin = kCTripletBegin,        kCQ1End = kCQ1Begin-637., kCQ1SqRadius = 3.5*3.5;
  56 static const Float_t kCQ2Begin = kCTripletBegin-908.5,  kCQ2End = kCQ2Begin-550., kCQ2SqRadius = 3.5*3.5;
  57 static const Float_t kCQ3Begin = kCTripletBegin-1558.5, kCQ3End = kCQ3Begin-550., kCQ3SqRadius = 3.5*3.5;
  58 static const Float_t kCQ4Begin = kCTripletBegin-2400.,  kCQ4End = kCQ4Begin-637., kCQ4SqRadius = 3.5*3.5;
  59 //
  60 static const Float_t kCD1Begin = -5838.3,  kCD1End = kCD1Begin-945., kCD1SqRadius = 4.5*4.5;
  61 static const Float_t kCD2Begin = -12167.8, kCD2End = kCD2Begin-945., kCD2SqRadius = 4.5*4.5;
  62 static const Float_t kCD2XCentre1 = -9.7;
  63 static const Float_t kCD2XCentre2 =  9.7;
  64 //
  65 // -> SIDE A
  66 // NB -> kACorr1Begin = 919. to be checked
  67 static const Float_t kACorr1Begin = 919., kACorr1End = kACorr1Begin+260., kCCorr1SqRadius = 4.*4.;
  68 static const Float_t kACorr2Begin = 1972.5, kACorr2End = kACorr2Begin+153., kCCorr2SqRadius = 4.5*4.5;
  69 static const Float_t kATripletBegin  = 2296.5;
  70 static const Float_t kAQ1Begin = kATripletBegin,        kAQ1End = kAQ1Begin+637., kAQ1SqRadius = 3.5*3.5;
  71 static const Float_t kAQ2Begin = kATripletBegin+908.5,  kAQ2End = kAQ2Begin+550., kAQ2SqRadius = 3.5*3.5;
  72 static const Float_t kAQ3Begin = kATripletBegin+1558.5, kAQ3End = kAQ3Begin+550., kAQ3SqRadius = 3.5*3.5;
  73 static const Float_t kAQ4Begin = kATripletBegin+2400.,  kAQ4End = kAQ4Begin+637., kAQ4SqRadius = 3.5*3.5;
  74 //
  75 static const Float_t kAD1Begin = 5838.3,  kAD1End = kAD1Begin+945., kAD1SqRadius = 3.375*3.375;
  76 static const Float_t kAD2Begin = 12167.8, kAD2End = kAD2Begin+945., kAD2SqRadius = 3.75*3.75;
  77 static const Float_t kAD2XCentre1 = -9.4;
  78 static const Float_t kAD2XCentre2 =  9.4;
  79
  80 // %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%  ZDC part  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  81
  82 #endif
  83