3 // Dynamically link some shared libs
4 gSystem->Load("$(ALICE_ROOT)/lib/tgt_$(ALICE_TARGET)/libLHC");
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.);
16 // The Interaction regions
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);
29 lhc->AddIRegion(alice);
30 // lhc->AddIRegion(cms);
31 // lhc->AddIRegion(atlas);
36 Float_t epsH = 1.5e-4;
37 Float_t epsL = 2.5e-9;
39 Float_t energy = 7000;
47 AliLhcBeam* beam1 = new AliLhcBeam(lhc);
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);
55 AliLhcBeam* beam2 = new AliLhcBeam(lhc);
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);
64 lhc->SetBeams(beam1, beam2);
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");
76 // lhc->AddProcess(bt);
80 lhc->SetTime(100., 20.*3600.);