| | 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 | |
git://git.uio.no / u / mrichter / AliRoot.git / blame_incremental