03192624e89cba3387dd9ac235ea12a716b8a6a0
[u/mrichter/AliRoot.git] / PHOS / PHOS-HTML / simulation.html
1 <!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN"
2                "http://www.w3.org/TR/REC-html40/loose.dtd">
3 <html>
4 <head>
5    <meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
6    <meta name="GENERATOR" content="Mozilla/4.51 [en] (X11; I; Linux 2.2.5-15 i686) [Netscape]">
7    <title>PHOS Reconstruction Proposal</title>
8 <link REL="stylesheet" href="http://www-subatech.in2p3.fr/~photons/gps_alice.css" type="text/css">
9 </head>
10 <body>
11
12 <h1 class="page-header">
13 Simulation</h1>
14
15 <h2>Event generation</h2>
16
17 <P>
18 One of the main ingredients to simulate the physic performance of PHOS is the
19 generation of a realistic background in central Pb+Pb collisions at LHC energies. <BR> <BR>
20 However, most of ALICE subdetectors are sensitive to charged particles and, in
21 consequence, the charged hadronic background is the most important one (see for
22 instance <A HREF="http://www.cern.ch/ALICE/Projects/offline/evgen/Welcome.html">AliGenHIJINGPara class</A>). <BR>
23 However, PHOS detector is sensitive to most of particles: electrons, muon, photons, charged and neutral hadrons. In
24 particular the main background is generated by the electromagnetic decay of the neutral pion, which was omitted 
25 in the background generators classes defined in aliroot package 3.01.  <BR>
26 <BR>
27 Within this respect, we have defined a new class: 
28 <A HREF="ttp://www.cern.ch/ALICE/Projects/offline/aliroot/roothtml/src/AliGenPHOSlib.cxx.html">AliGenPHOSlib class</A>,
29 and adapted the 
30 <A HREF="http://www.cern.ch/ALICE/Projects/offline/aliroot/roothtml/src/AliGenParam.cxx.html">AliGenParam class</A> class in
31 order to be able to simulate the full neutral and charged background which is present in heavy-ion collisions.
32 All these stuff has been included in the aliroot version 3.02.<BR>
33 <BR>
34
35 Enclosed, an exemple of the aliroot config file is reported. 
36 This neutral and charged background contains pi+, pi-,pi0,K+,K-,kshort,Klong,eta,omega,p. p bar, n and n bar with 
37 relative yields given by hijing in Pb+Pb central collisions (results from venus or shaker are close to the hijing ones).
38 Pion transverse momentum (Pt) distribution is taken from AliGenMUONlib class, version 3.01 of aliroot: 
39 Pt Parameterization CDF (PRL 61(88) 1819):  Power law  for Pt > 500 MeV and Mt scaling below (slope T=160 MeV).
40 Meson Pt distributions are calculated by mt-scaling with pions (see 
41 <A HREF="http://www.cern.ch/ALICE/Projects/offline/evgen/Welcome.html">AliGenMUONlib or AliGenPHOSlib class</A>). 
42
43 <PRE class="code">
44 //=======================================================================
45 // ************* STEERING parameters FOR ALICE SIMULATION **************
46 // --- Specify event type to be tracked through the ALICE setup
47 // --- All positions are in cm, angles in degrees, and P and E in GeV
48 //
49 // The following Cocktail generator is defined to simulate the neutral and
50 // charged background in the ALICE detector. This background is important 
51 // in the case of photon detector as PHOS. We simulated a cocktail of 
52 // pions (pi+, pi- and pi0) , kaons (K+, K-, Kshort and Klong), eta mesons, 
53 // omega mesons and main baryons (protons, antiprotons, neutrons and
54 // antineutrons) 
55 //
56 // 1-Nov-1999 Gines MARTINEZ, GPS @ SUBATECH, Nantes, France  
57 //
58 Int_t ParticleDensity =  8000    // Number of particles created in the selected rapidity range
59 Int_t NumberOfPions   =  0.77 * ParticleDensity ;  // ~77% of pions: pi^+, pi^- and pi^0 with equal probability 
60 Int_t NumberOfKaons   =  0.11 * ParticleDensity ;  // ~11% of kaons: K^+, K^-, K_short and K_long with equal probability
61 Int_t NumberOfEtas    =  0.04 * ParticleDensity ;  // ~4.0% of eta mesons
62 Int_t NumberOfOmegas  =  0.01 * ParticleDensity ;  // ~1.0% of omega mesons
63 Int_t NumberOfBaryons =  0.07 * ParticleDensity ;  // ~7.0% of baryons: p, pbar, n and nbar with equal probability
64
65 AliGenCocktail *gener = new AliGenCocktail();                                \\ Cocktail class
66  gener->SetPtRange(.5,5.);       // Transverse momentum range   
67  gener->SetPhiRange(0.,360.);    // Azimuthal angle range  
68  gener->SetYRange(-0.5,0.5);     // Pseudorapidity range
69  gener->SetOrigin(0,0,0);        // Vertex position
70  gener->SetSigma(0,0,5.6);       // Sigma in (X,Y,Z) (cm) on IP position
71 //
72 // P  I  O  N  S
73     AliGenParam *generpion = new AliGenParam(NumberOfPions,Pion,
74                  AliGenPHOSlib::GetPt(Pion), AliGenPHOSlib::GetY(Pion), AliGenPHOSlib::GetIp(Pion)  );
75     generpion->SetWeighting(non_analog);   // Selecting the Pt distribution provided by AliGenPHOSlib::GetPt
76     generpion->SetForceDecay(nodecay);     // Requiring that the generated particle is directly tracked by GEANT
77 //
78 //  K  A  O  N  S
79    AliGenParam *generkaon = new AliGenParam(NumberOfKaons,Kaon, 
80                 AliGenPHOSlib::GetPt(Kaon), AliGenPHOSlib::GetY(Kaon), AliGenPHOSlib::GetIp(Kaon)   );
81     generkaon->SetWeighting(non_analog);   // Selecting the Pt distribution provided by AliGenPHOSlib::GetPt
82     generkaon->SetForceDecay(nodecay);     // Requiring that the generated particle is directly tracked by GEANT
83 //
84 // E  T  A  S
85     AliGenParam *genereta = new AliGenParam(NumberOfEtas,Eta,            
86                  AliGenPHOSlib::GetPt(Eta), AliGenPHOSlib::GetY(Eta), AliGenPHOSlib::GetIp(Eta) );
87     genereta->SetWeighting(non_analog);   // Selecting the Pt distribution provided by AliGenPHOSlib::GetPt
88     genereta->SetForceDecay(nodecay);     // Requiring that the generated particle is directly tracked by GEANT
89 //
90 // O M E G A S
91     AliGenParam *generomega = new AliGenParam(NumberOfOmegas,Omega,            
92                  AliGenPHOSlib::GetPt(Omega), AliGenPHOSlib::GetY(Omega), AliGenPHOSlib::GetIp(Omega) );
93     generomega->SetWeighting(non_analog);   // Selecting the Pt distribution provided by AliGenPHOSlib::GetPt
94     generomega->SetForceDecay(nodecay);    // Requiring that the generated particle is directly tracked by GEANT
95 //
96 // B A R Y O N S
97
98     AliGenParam *generbaryon = new AliGenParam(NumberOfBaryons,Baryon,            
99                  AliGenPHOSlib::GetPt(Baryon), AliGenPHOSlib::GetY(Baryon), AliGenPHOSlib::GetIp(Baryon) );
100     generbaryon->SetWeighting(non_analog);   // Selecting the Pt distribution provided by AliGenPHOSlib::GetPt
101     generbaryon->SetForceDecay(nodecay);    // Requiring that the generated particle is directly tracked by GEANT
102     
103     
104   gener->AddGenerator(generpion,"pion",1.);
105   gener->AddGenerator(generkaon,"kaon",1.);
106   gener->AddGenerator(genereta,"eta",1.);
107   gener->AddGenerator(generomega,"omega",1.);
108   gener->AddGenerator(generbaryon,"baryon",1.);
109   gener->Init();
110 </pre>
111
112
113 The resulting hit and digit TTree's (one per event) will be stored
114 in the root file on disk.
115
116 <hr>
117 <address class="left">
118 &copy; <a href="mailto:aphecetc@in2p3.fr">Groupe Photons Subatech</a> <a href="http://www-subatech.in2p3.fr/~photons/subatech/en_index.shtml">[Go
119 to the GPS Home Page]</a></address>
120
121 <address class="right">
122 <!-- Created: Tue Oct 26 19:52:56 CEST 1999 -->
123 <!-- hhmts start -->
124 Last modified: Mon Nov 15 18:55:06 CET 1999
125 <!-- hhmts end --></address>
126
127 <div align=right><a href="http://validator.w3.org/check/referer"><img SRC="vh40.gif" ALT="Valid HTML 4.0!" BORDER=0 height=31 width=88></a><a href="http://jigsaw.w3.org/css-validator"><img SRC="vcss.gif" ALT="Valid CSS!" BORDER=0 height=31 width=88></a></div>
128
129 </body>
130 </html>