]>
Commit | Line | Data |
---|---|---|
1 | ======== Continuous readout toy simulation ======= | |
2 | ||
3 | loadlibs.C //modified version of the one used for GEMtest code. I haven't looked into this so it probably contains a lot of stuff that is not needed | |
4 | ||
5 | AliToyMCTrack, AliToyMCEvent //classes for storing event and track parameters | |
6 | ||
7 | AliToyMCEventGenerator //base class for event generators, contains DistortTrack method to distort track using method borrowed from AliTPCCorrection | |
8 | ||
9 | AliToyMCEventGeneratorSimple //simple first MC event generator, borrowing from Marians macro. Inherits from AliToyMCEventGenerator to get DistortTrack method. Generate(Double_t time) generates and returns an event at time time | |
10 | ||
11 | makeTree.C //short macro to generate events from a poisson distr and storing in a ROOT tree. | |
12 | ||
13 | ||
14 | to run : root -l loadlibs.C | |
15 | .L makeTree.C+ | |
16 | makeTree(Double_t collFreq/*kHz*/, Double_t bunchFreq/*MHz*/, Int_t nEvents) | |
17 | ||
18 | ||
19 | ======== Continuous readout toy reconstruction ======= | |
20 | AliToyMCReconstruction -- Class to reconstruct | |
21 | ||
22 | ||
23 | ========== Simple example ========= | |
24 | 1. simultate | |
25 | ||
26 | AliToyMCEventGeneratorSimple s | |
27 | s.SetSpaceCharge(AliToyMCEventGeneratorSimple::kEps5) | |
28 | s.SetOutputFileName("toyMC_eps05_10ev_50tr.NoMat.root"); | |
29 | s.RunSimulation(10,50); | |
30 | ||
31 | 2. run reconstruction | |
32 | const char *file="toyMC_eps05_50ev_500tr.root" | |
33 | AliToyMCReconstruction rec | |
34 | rec.SetUseMaterialBudget(kTRUE) | |
35 | ||
36 | rec.SetRecoSettings(1,0,AliToyMCReconstruction::kNoCorrection) | |
37 | rec.RunReco(file,-1); > /tmp/log_1_0_0 | |
38 | ||
39 | rec.SetRecoSettings(1,1,AliToyMCReconstruction::kIdeal) | |
40 | rec.RunReco(file,-1); > /tmp/log_1_1_3 | |
41 | ||
42 | rec.SetRecoSettings(0,1,AliToyMCReconstruction::kIdeal) | |
43 | rec.RunReco(file,-1); > /tmp/log_0_1_3 | |
44 | ||
45 | rec.SetRecoSettings(0,1,AliToyMCReconstruction::kAverageEta) | |
46 | rec.RunReco(file,-1); > /tmp/log_0_1_2 | |
47 | ||
48 | rec.SetRecoSettings(0,1,AliToyMCReconstruction::kNoCorrection) | |
49 | rec.RunReco(file,-1); > /tmp/log_0_1_0 | |
50 | ||
51 | rec.SetRecoSettings(0,0,AliToyMCReconstruction::kNoCorrection) | |
52 | rec.RunReco(file,-1); > /tmp/log_0_0_0 | |
53 | ||
54 | 3. visualisation | |
55 | .L $ALICE_ROOT/TPC/Upgrade/macros/toyMCRecPlots.C | |
56 | toyMCRecPlots("toyMC_eps20_50ev_500tr.0.1_0_0_140_10.debug.root") | |
57 |