LHC/lumi.C

   1 void lumi ()
   2 {
   3 // Dynamically link some shared libs
   4     gSystem->Load("$(ALICE_ROOT)/lib/tgt_$(ALICE_TARGET)/libLHC");
   5
   6 //
   7 //  LHC
   8 //
   9     AliLHC* lhc = new AliLHC();
  10     lhc->SetRadius(2665887./2./TMath::Pi());
  11     lhc->SetAverageBeta(5050.);
  12     lhc->SetAverageDisp(129.);
  13     lhc->SetSetUpTime(1./3. * 3600.);
  14     lhc->SetFillingTime(3.  * 3600.);
  15 //
  16 //  The Interaction regions
  17 //
  18 //
  19
  20     Float_t bstar = 50.;
  21 //
  22     AliLhcIRegion* alice =  new AliLhcIRegion(lhc, "ALICE", "ALICE");
  23     alice->SetBetaStar(bstar);
  24     AliLhcIRegion* cms   =  new AliLhcIRegion(lhc, "CMS"  , "CMS"  );
  25     cms->SetBetaStar(bstar);
  26     AliLhcIRegion* atlas =  new AliLhcIRegion(lhc, "ATLAS", "ATLAS");
  27     atlas->SetBetaStar(bstar);
  28
  29     lhc->AddIRegion(alice);
  30 //    lhc->AddIRegion(cms);
  31 //    lhc->AddIRegion(atlas);
  32 //
  33 //  The beams
  34 //
  35     Float_t n         = 6.8e7;
  36     Float_t epsH      = 1.5e-4;
  37     Float_t epsL      = 2.5e-9;
  38     Float_t de        = 1.14e-4;
  39     Float_t energy    = 7000;
  40     Int_t a1          = 208;
  41     Int_t z1          =  82;
  42     Int_t a2          = 208;
  43     Int_t z2          =  82;
  44
  45
  46 //
  47     AliLhcBeam* beam1 =  new AliLhcBeam(lhc);
  48     beam1->SetN(n);
  49     beam1->SetNEmittance(epsH);            // (cm)
  50     beam1->SetLongEmittance(epsL);         // (GeV s)
  51     beam1->SetEnergy(energy);                 // (GeV)
  52     beam1->SetParticle(a1,z1);
  53     beam1->SetEnergySpread(de);
  54
  55     AliLhcBeam* beam2 =  new AliLhcBeam(lhc);
  56     beam2->SetN(n);
  57     beam2->SetNEmittance(epsH);            // (cm)
  58     beam2->SetLongEmittance(epsL);         // (GeV s)
  59     beam2->SetEnergy(energy);                 // (GeV)
  60     beam2->SetParticle(a2,z2);
  61     beam2->SetEnergySpread(de);
  62
  63
  64     lhc->SetBeams(beam1, beam2);
  65 //
  66 //  The Processes
  67 //
  68     AliLhcProcessBB*  bb  = new AliLhcProcessBB(lhc, "BB", "Beam-Beam Losses");
  69     bb->SetCrossSection(505.);
  70     AliLhcProcessIBS* ibs = new AliLhcProcessIBS(lhc, "IBS", "Intra Beam Scattering");
  71     AliLhcProcessBT*  bt  = new AliLhcProcessBT(lhc, "BT", "Beta* Tuning");
  72     bt->SetBetaMin(50.);
  73
  74     lhc->AddProcess(ibs);
  75     lhc->AddProcess(bb);
  76 //    lhc->AddProcess(bt);
  77 //
  78 //  Run the collider
  79 //
  80     lhc->SetTime(100., 20.*3600.);
  81     lhc->Init();
  82     lhc->EvolveTime();
  83     lhc->Evaluate();
  84 }
  85
  86
  87
  88
  89
  90
  91
  92
  93
  94