]>
Commit | Line | Data |
---|---|---|
50837721 | 1 | /************************************************************************** |
2 | * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. * | |
3 | * * | |
4 | * Author: The ALICE Off-line Project. * | |
5 | * Contributors are mentioned in the code where appropriate. * | |
6 | * * | |
7 | * Permission to use, copy, modify and distribute this software and its * | |
8 | * documentation strictly for non-commercial purposes is hereby granted * | |
9 | * without fee, provided that the above copyright notice appears in all * | |
10 | * copies and that both the copyright notice and this permission notice * | |
11 | * appear in the supporting documentation. The authors make no claims * | |
12 | * about the suitability of this software for any purpose. It is * | |
13 | * provided "as is" without express or implied warranty. * | |
14 | **************************************************************************/ | |
15 | ||
16 | /* $Id$ */ | |
17 | ||
88cb7938 | 18 | ========================================================== |
19 | Please add to this README file all information concerning | |
20 | config files, simulation, digitalization, clusterization, | |
21 | reconstruction and macro analysis | |
22 | ||
6b1e4b22 | 23 | ========================================================== |
24 | How to check that your aliroot is working well | |
25 | ========================================================== | |
26 | There is a script file AlirootRun_MUONtest.script which | |
27 | allows for simulating, reconstructing and making the | |
28 | invariant analysis of the generated Upsilon (1S). | |
29 | The used configuration file is Config.C in MUON | |
30 | directory. | |
31 | There you have to type : | |
32 | ||
33 | source $ALICE_ROOT/MUON/AlirootRun_MUONtest.script | |
34 | ||
35 | If you do not recover a few Upsilons in their mass region | |
36 | in the file MUONmassPlot.root | |
37 | ||
88cb7938 | 38 | ========================================================== |
7985603c | 39 | How to run a MUON generation |
88cb7938 | 40 | ========================================================== |
41 | aliroot | |
fc1fafa5 | 42 | root [0] gAlice->Run(10,"$ALICE_ROOT/MUON/Config.C"); |
88cb7938 | 43 | |
44 | 1 single muon of 7 GeV/c in the MUON spectrometer | |
45 | acceptance will be simulated using geant3. | |
46 | Hit information will be store in the root file in the | |
47 | execution directory. | |
88cb7938 | 48 | If you want to change the option or to define a new directory |
49 | for hits, you have to do the following before: | |
50 | root [0] gAlice->SetConfigFunction("Config( \"/home/martinez/aliroot/work_NewIO/test/\" , \"box\" );"); | |
51 | ||
a3448039 | 52 | |
88cb7938 | 53 | ============================================================ |
54 | How to run MUONCheck macro | |
55 | ============================================================ | |
56 | To check the content of a root data file, the MUONCheck | |
57 | provides a ascii output on screen. | |
58 | ||
59 | To compile MUONCheck.C | |
2058bf30 | 60 | .includepath $ALICE_ROOT/STEER |
61 | .includepath $ALICE_ROOT/MUON | |
62 | .L $ALICE_ROOT/MUON/MUONCheck.C++ | |
bf17dbfd | 63 | To Load |
64 | gSystem->Load("$ALICE_ROOT/MUON/MUONCheck_C.so") | |
2058bf30 | 65 | |
278a86d5 | 66 | To print Kine : (default file is galice.root ) |
2b32c661 | 67 | MUONkine() or MUONkine("galice.root",##) for the event number ## |
88cb7938 | 68 | |
69 | To print hits : (default file is galice.root if not MUONhits("toto.root""); ) | |
2b32c661 | 70 | MUONhits() or MUONhits("galice.root",##) for the event number ## |
88cb7938 | 71 | |
72 | To print digits : (default file is galice.root) | |
2b32c661 | 73 | MUONdigits() or MUONdigits("galice.root",##) for the event number ## |
88cb7938 | 74 | |
75 | To print rawcluster : (default file is galice.root) | |
2b32c661 | 76 | MUONrecpoints() or MUONrecpoints("galice.root",##) for the event number ## |
88cb7938 | 77 | |
78 | To print trigger : (default file is galice.root) | |
2b32c661 | 79 | MUONTestTrigger() or MUONTestTrigger("galice.root",##) for the event number ## |
80 | ||
81 | .... | |
88cb7938 | 82 | |
02d8f072 | 83 | |
84 | ============================================================ | |
85 | How to check the Geometry | |
86 | ============================================================ | |
fc1fafa5 | 87 | gAlice->Init("$ALICE_ROOT/MUON/Config.C"); |
02d8f072 | 88 | .L $ALICE/geant3/TGeant3/G3GUI.C |
89 | G3GUI() | |
90 | ||
91 | ||
7985603c | 92 | ============================================================ |
93 | How to run MUONdisplay | |
94 | ============================================================ | |
1eccde20 | 95 | First you need to perform a full simulation: |
7985603c | 96 | generation, digitalisation and clusterisation |
97 | .L $ALICE_ROOT/MUON/MUONdisplay.C | |
98 | MUONdisplay(0,"galice.root") | |
99 | ||
02d8f072 | 100 | ============================================================ |
101 | Tracking parameters, cuts, energy loss and physics processes | |
102 | ============================================================ | |
103 | Tracking parameters in MUON are automatically defined by GEANT | |
104 | MUON takes the default values of CUTs and physics processes | |
105 | defined by the Config files, except for the gas mixture medium | |
106 | of the tracking chambers. The CUT's and physics processes of | |
107 | the gas mixture medium is then defined in the galice.cuts file | |
108 | in the data directory. In particular ILOSS parameter MUST be | |
109 | equal unity (1) in order simulate a realistic energy loss | |
110 | distribution (mean value and fluctuations) in the active gas. | |
a88eb0d0 | 111 | |
112 | ============================================================ | |
113 | Tracking of particle in the magnetic field | |
114 | ============================================================ | |
115 | GEANT has two ways for tracking charged particles in the | |
116 | magnetic field: HELIX et RKUTA. | |
117 | HELIX is faster and works well if the gradient of magnetic | |
118 | field is small. | |
119 | For MUON, HELIX is a not a good approximation and we must | |
120 | use RKUTA to get the optimal mass resolution of the | |
121 | spectrometer. The choice of HELIX or RKUTA is done in the | |
122 | config file when the magnetic field is defined: | |
123 | AliMagFMaps* field = new AliMagFMaps("Maps","Maps", TRACKING, FACTOR, MAXB, AliMagFMaps::k4kG); | |
124 | gAlice->SetField(field); | |
125 | TRACKING must be 1 for RKUTA and 2 for HELIX (the default value for aliroot is 2 (HELIX)) | |
126 | FACTOR allows you to set the magnetic field to 0, just putting FACTOR=0. Default value is 1. | |
127 | MAXB is the maximum magnetic field which is 10.T | |
2b32c661 | 128 | |
f4f795ed | 129 | =========================================================== |
130 | MUON cocktail for physics .............. | |
131 | =========================================================== | |
132 | There is a MUON cocktail generator of the muon sources in the | |
133 | EVGEN directory. This class derives from AliGenCocktail. | |
134 | In the init of this class I have filled the cocktail with | |
135 | the muon sources: J/Psi, Upsilon, Open Charm, Open Beauty, | |
136 | Pion, Kaons. The code needs only the production cross section | |
137 | at 4pi (for the moment this values are in the code since I | |
138 | prefere them do not be modified), and the code calculates the | |
139 | rate of particles in the acceptance, making the scaling based | |
140 | on the number of collisions for the hard probes and on the | |
141 | number of participants for soft sources: Pions and Kaons. | |
142 | ||
143 | In the Genereate of this class all entries in the cocktail | |
144 | are called and we define a "primordial trigger" with requires | |
145 | a minimum number of muons above a Pt cut in the required acceptance. | |
146 | In order to normalized to the real number of simulated events, | |
147 | there are 2 data members in the class fNsuceeded adn fNGenerate | |
148 | which tell us what is the biais source. | |
149 | ||
150 | Enclose an example to use this generator: | |
151 | AliGenMUONCocktail * gener = new AliGenMUONCocktail(); | |
152 | gener->SetPtRange(1.,100.); // Transverse momentum range | |
153 | gener->SetPhiRange(0.,360.); // Azimuthal angle range | |
154 | gener->SetYRange(-4.0,-2.4); | |
155 | gener->SetMuonPtCut(1.); | |
156 | gener->SetMuonThetaCut(171.,178.); | |
157 | gener->SetMuonMultiplicity(2); | |
35e21dec | 158 | gener->SetImpactParameterRange(0.,5.); // 10% most centra PbPb collisions |
f4f795ed | 159 | gener->SetVertexSmear(kPerTrack); |
160 | gener->SetOrigin(0,0,0); // Vertex position | |
161 | gener->SetSigma(0,0,0.0); // Sigma in (X,Y,Z) (cm) on IP position | |
162 | gener->Init(); | |
163 | ||
2b32c661 | 164 | =========================================================== |
69be760c | 165 | csh Script for the full reconstruction with raw data generator |
2b32c661 | 166 | =========================================================== |
2b32c661 | 167 | |
9f532ceb | 168 | |
5473f16a | 169 | aliroot -b << EOF |
170 | AliSimulation MuonSim("YourConfig.C") | |
69be760c | 171 | MuonSim.SetWriteRawData("MUON"); |
5473f16a | 172 | MuonSim.Run(XXX) |
2b32c661 | 173 | .q |
174 | EOF | |
175 | ||
cac2eb58 | 176 | aliroot -b << EOF |
177 | TPluginManager* pluginManager = gROOT->GetPluginManager(); | |
178 | pluginManager->AddHandler("AliReconstructor", "MUON","AliMUONReconstructor", "MUON","AliMUONReconstructor()") | |
179 | AliReconstruction MuonRec("galice.root") | |
d59a4389 | 180 | MuonRec.SetRunTracking("") |
5473f16a | 181 | MuonRec.SetRunVertexFinder(kFALSE) |
182 | MuonRec.SetRunLocalReconstruction("MUON") | |
183 | MuonRec.SetFillESD("MUON") | |
cac2eb58 | 184 | MuonRec.Run() |
2b32c661 | 185 | .q |
186 | EOF | |
187 | ||
2b32c661 | 188 | aliroot -b << EOF |
189 | .includepath $ALICE_ROOT/STEER | |
190 | .includepath $ALICE_ROOT/MUON | |
5473f16a | 191 | .L $ALICE_ROOT/MUON/MUONmassPlot_ESD.C++ |
2b32c661 | 192 | MUONmassPlot("galice.root",0,99999); |
193 | .q | |
194 | EOF | |
a88eb0d0 | 195 | |
b8dc484b | 196 | |
197 | ============================================================ | |
198 | How to run MUONRecoCheck macro | |
199 | ============================================================ | |
200 | To check the muon reconstruction by comparing the reconstructed tracks | |
201 | with the reference tracks made of "AliTrackReference" for the hits and | |
202 | kinematic informations (TParticle) for the vertex. | |
203 | This macro can be used to check the track reconstruction e.g. efficiency, | |
204 | momentum resolution ... but also to make physics analysis whenever | |
205 | track identification is needed. | |
206 | ||
207 | To compile MUONRecoCheck.C | |
208 | .includepath $ALICE_ROOT/STEER | |
209 | .includepath $ALICE_ROOT/MUON | |
210 | .L $ALICE_ROOT/MUON/MUONRecoCheck.C+ | |
211 | ||
212 | // To run MUONRecoCheck | |
213 | MUONRecoCheck(nEvent,"galice.root"); // nEvent = nb of events | |
214 | ||
29fc2c86 | 215 | |
216 | ============================================================ | |
217 | How to run MUONTracker macro | |
218 | ============================================================ | |
219 | To make the track reconstruction directly from AliTrackReference hits | |
220 | which are recorded in TrackRefs.root during the simulation. | |
221 | It can be used to check the reconstruction without clusterization. | |
222 | ||
223 | To compile MUONTracker.C | |
224 | .includepath $ALICE_ROOT/STEER | |
225 | .includepath $ALICE_ROOT/MUON | |
226 | .L $ALICE_ROOT/MUON/MUONTracker.C+ | |
227 | ||
228 | // To run MUONTracker | |
229 | MUONTracker(iEventMin,iEventMax,"galice.root"); // iEventMin: first event | |
fc1fafa5 | 230 | |
231 | =========================================================== | |
232 | Macro MUONGenerateGeometryData.C | |
233 | =========================================================== | |
234 | ||
235 | Macro for generating the geometry data files: | |
236 | (transform_*.dat, svmap_*.dat). | |
237 | - MUON/data/transform_*dat files contain all the information | |
238 | about the transformation (traslation + rotation) needed for | |
239 | the positionning of each detection element of the MUON | |
240 | spectrometer. | |
241 | - MUON/data/svmap_*.dat files contain all the information to link | |
242 | each geant volume (it can be extended to other virtual MC) with | |
243 | a detection element. The point here is that a given detection | |
244 | element, i.e. a slat chamber consists of many geant volumes. | |
245 | the correspondence is then defined in an input file. | |
246 | Eanch time there is a change in the virtual MC geometry, these | |
247 | input files must be re-generated via the macro | |
248 | MUONGenerateGeometryData.C | |
249 | ||
250 | To be run from aliroot: | |
251 | .x MUONGenerateGeometryData.C | |
252 | ||
253 | The generated files do not replace the existing ones | |
254 | but have different names (with extension ".out"). | |
255 | To compare/replace the existing files with generated ones | |
256 | run the scripts in MUON/data: | |
257 | compare_data [transform] [svmaps ] | |
258 | reset_data [transform] [svmaps] | |
259 | Author: I. Hrivnacova, IPN Orsay | |
29fc2c86 | 260 | |
261 | ||
88cb7938 | 262 | =========================================================== |
263 | Still working .............. | |
264 | =========================================================== |