--- /dev/null
+void lumi ()
+{
+// Dynamically link some shared libs
+ gSystem->Load("$(ALICE_ROOT)/lib/tgt_$(ALICE_TARGET)/libLHC");
+
+//
+// LHC
+//
+ AliLHC* lhc = new AliLHC();
+ lhc->SetRadius(2665887./2./TMath::Pi());
+ lhc->SetAverageBeta(5050.);
+ lhc->SetAverageDisp(129.);
+ lhc->SetSetUpTime(1./3. * 3600.);
+ lhc->SetFillingTime(3. * 3600.);
+//
+// The Interaction regions
+//
+//
+
+ Float_t bstar = 50.;
+//
+ AliLhcIRegion* alice = new AliLhcIRegion(lhc, "ALICE", "ALICE");
+ alice->SetBetaStar(bstar);
+ AliLhcIRegion* cms = new AliLhcIRegion(lhc, "CMS" , "CMS" );
+ cms->SetBetaStar(bstar);
+ AliLhcIRegion* atlas = new AliLhcIRegion(lhc, "ATLAS", "ATLAS");
+ atlas->SetBetaStar(bstar);
+
+ lhc->AddIRegion(alice);
+// lhc->AddIRegion(cms);
+// lhc->AddIRegion(atlas);
+//
+// The beams
+//
+ Float_t n = 6.8e7;
+ Float_t epsH = 1.5e-4;
+ Float_t epsL = 2.5e-9;
+ Float_t de = 1.14e-4;
+ Float_t energy = 7000;
+ Int_t a1 = 208;
+ Int_t z1 = 82;
+ Int_t a2 = 208;
+ Int_t z2 = 82;
+
+
+//
+ AliLhcBeam* beam1 = new AliLhcBeam(lhc);
+ beam1->SetN(n);
+ beam1->SetNEmittance(epsH); // (cm)
+ beam1->SetLongEmittance(epsL); // (GeV s)
+ beam1->SetEnergy(energy); // (GeV)
+ beam1->SetParticle(a1,z1);
+ beam1->SetEnergySpread(de);
+
+ AliLhcBeam* beam2 = new AliLhcBeam(lhc);
+ beam2->SetN(n);
+ beam2->SetNEmittance(epsH); // (cm)
+ beam2->SetLongEmittance(epsL); // (GeV s)
+ beam2->SetEnergy(energy); // (GeV)
+ beam2->SetParticle(a2,z2);
+ beam2->SetEnergySpread(de);
+
+
+ lhc->SetBeams(beam1, beam2);
+//
+// The Processes
+//
+ AliLhcProcessBB* bb = new AliLhcProcessBB(lhc, "BB", "Beam-Beam Losses");
+ bb->SetCrossSection(505.);
+ AliLhcProcessIBS* ibs = new AliLhcProcessIBS(lhc, "IBS", "Intra Beam Scattering");
+ AliLhcProcessBT* bt = new AliLhcProcessBT(lhc, "BT", "Beta* Tuning");
+ bt->SetBetaMin(50.);
+
+ lhc->AddProcess(ibs);
+ lhc->AddProcess(bb);
+// lhc->AddProcess(bt);
+//
+// Run the collider
+//
+ lhc->SetTime(100., 20.*3600.);
+ lhc->Init();
+ lhc->EvolveTime();
+ lhc->Evaluate();
+}
+
+
+
+
+
+
+
+
+
+