]>
Commit | Line | Data |
---|---|---|
7bd443c7 | 1 | ======== Continuous readout toy simulation ======= |
526ddf0e | 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 | ||
8a218638 | 5 | AliToyMCTrack, AliToyMCEvent //classes for storing event and track parameters |
526ddf0e | 6 | |
8a218638 | 7 | AliToyMCEventGenerator //base class for event generators, contains DistortTrack method to distort track using method borrowed from AliTPCCorrection |
526ddf0e | 8 | |
8a218638 | 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 |
526ddf0e | 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+ | |
8597ab27 | 16 | makeTree(Double_t collFreq/*kHz*/, Double_t bunchFreq/*MHz*/, Int_t nEvents) |
17 | ||
18 | ||
7bd443c7 | 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 |