1 / 11

MC simulation of test beam events

MC simulation of test beam events. Fabrizio Salvatore (RHUL). Simulation of Desy tb events - I. First attempt to simulate MC events in test beam mode (Desy06 test beam) at three different institutions: RHUL, B’ham, Manc.

cholliday
Download Presentation

MC simulation of test beam events

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. MC simulation of test beam events Fabrizio Salvatore (RHUL)

  2. Simulation of Desy tb events - I • First attempt to simulate MC events in test beam mode (Desy06 test beam) at three different institutions: RHUL, B’ham, Manc. • Use same version of Mokka (v06-01) and same version of packages against which Mokka is linked: • Geant4 8.0.p01 • CLHEP 2.0.2.3 • Note: use same versions of g++ compiler (gcc 3.2.3), but different updates • RHUL: 3.2.3-56; Bham: 3.2.3-49; Man: 3.2.3-49/56 • LCIO v01-07 • Gear 00-03 • GSL 1.7

  3. Simulation of Desy tb events - II • Small modification needed in Geant4 and in Mokka/Kernel in order to define initial seeds used by the random numbers generator (David Bailey): • G4: in source/intercoms/include/G4UIcommand.hh • Mokka: Added class ‘G4UIcmdWithVariableIntList’ defined in Kernel/src/G4UIcmdWithVariableIntList.cc • In src/ControlMessenger.cc: G4int DoIt(G4String parameterList); → virtual G4int DoIt(G4String parameterList); //Random number commands _randomSeedsCommand = new G4UIcmdWithVariableIntList("/Mokka/init/RandomNumberSeeds",this); _randomSeedsCommand->SetGuidance("List of integer seeds to be passed to the random number generator."); _randomSeedsCommand->AvailableForStates(G4State_PreInit);

  4. Simulation of Desy tb events - III • Generate 10K 1GeV electron events using TBDesy0506, at normal incidence • Use same initial random numbers seeds • Use caliceMarlin to extract distribution of variables related to ECAL and tracking to compare the events (lcio output) generated in the three different sites (thanks to Michele and David W. for the help on setting up caliceMarlin) • Plot variable distribution and ratio of variable_at_site_i – variable_at_site_ j variable_at_site_i

  5. N(hits) vs Plane RHUL B’ham RHUL Man Man B’ham

  6. Energy vs Plane RHUL B’ham RHUL Man Man B’ham

  7. N(hits) vs i ECAL RHUL Man RHUL B’ham Man B’ham

  8. N(hits) vs j ECAL RHUL B’ham RHUL Man Man B’ham

  9. r(hit) and E(hit/mips) Man B’ham RHUL B’ham RHUL Man Man B’ham RHUL B’ham Man B’ham RHUL Man

  10. Summary - I • Full comparison plots available at: http://www.pp.rhul.ac.uk/~calice/fab/Mokka/test/ • Although same versions of all programs and same initial seeds for the random numbers generator are used, resulting MC events are not exactly the same (as one would expect) • Differences in the distributions of several kinematic variables are almost all within statistical uncertainties • No big systematic bias is clearly visible

  11. Summary - II • Modified code to handle random number seeds not yet in official Mokka release • Modification in Geant4 needs to be included in standard G4 release • Need to setup a strategy for future massive MC production • Compile executable at one single site and distribute the executable to all institutes involved in MC production • Use different executables and quote differences as systematic effect on results

More Related