bug fixed for alignment, removed alignment database access from AliPMDUtility class
[u/mrichter/AliRoot.git] / LHC / lumi.C
CommitLineData
b6cd8e3f 1void 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