Instructions of NSW Simulation in lxplus or lxatut.
  1. Particle Generation
  2. Detector Simulation with NSW
  3. MM Digitization
  4. MM Trigger Simulation

Particle Generation

  • See /home/tsaito/work/MM_Digitization/ParticleGenerator/script_20151028.py

Detector Simulation with NSW

  • Setup
    • asetup 20.3.X-VAL, rel1, here (old: asetup 20.3.0.3,here)
  • Run
    • AtlasG4_tf.py --inputEVNTFile 'your input file' --preExec simFlags.SimulateNewSmallWheel=True --postInclude /afs/cern.ch/user/j/jomeyer/public/NSW/NSW.config.simu.py --postInclude /afs/cern.ch/user/j/jomeyer/public/NSW/NSW.config.simu.py --geometryVersion ATLAS-R2-2015-03-01-00_VALIDATION --outputHITSFile test.hits.pool.root --maxEvents 250
  • Input files
    • Single muon with pT=100 GeV, Flat distribution in eta and phi plane within NSW acceptance: /data/maxi162/atljphys/tsaito/MM_Digitization/signal/Generator/Flat_2.9_20151028/MC12.107206.ParticleGenerator_singlemu_Pt100/MC12.107206.ParticleGenerator_singlemu_Pt100.generator.pool.root (lxatut)
    • Other inputs: /afs/cern.ch/user/j/jomeyer/public/NSW/EVNT/
  • Output files
    • /data//maxi16s/atljphys/tsaito/MM_Digitization/NEWGEOMETRY_2015/Hits/Flat_2.9_20151029/output/*pool.root

MM Digitization (SVN)

  • Setup
    • cmt co svn+ssh://svn.cern.ch/reps/atlasoff/MuonSpectrometer/MuonDigitization/MM_Digitization/trunk MuonSpectrometer/MuonDigitization/MM_Digitization
    • cmt co svn+ssh://svn.cern.ch/reps/atlasoff/MuonSpectrometer/MuonValidation/MuonPRDTest/trunk MuonSpectrometer/MuonValidation/MuonPRDTest
    • To fix geometry errors ( needed in the latest one?)
      • cmt co svn+ssh://svn.cern.ch/reps/atlasoff/DetectorDescription/IdDictParser/trunk DetectorDescription/IdDictParser
      • In DetectorDescription/IdDictParser/data/IdDictMuonSpectrometer_S._02.xml, increased “mmChannel” maxvalue from 2500 to 8192

  • Compile
    • setupWorkArea.py
    • cd $TestArea/WorkArea/cmt
    • cmt br cmt config
    • cmt br gmake
    • cd $TestArea

  • Run
    • Reco_tf.py --postInclude /afs/cern.ch/user/j/jomeyer/public/NSW/NSW.config.digi.py,/afs/cern.ch/user/j/jomeyer/public/NSW/NSWPRDValAlg.digi.py --inputHITSFile 'your input file' --outputRDOFile test.rdos.pool.root
  • Input files
    • /data//maxi16s/atljphys/tsaito/MM_Digitization/NEWGEOMETRY_2015/Hits/Flat_2.9_20151029/output/*pool.root
  • Output files
    • /data//maxi16s/atljphys/tsaito/MM_Digitization/NEWGEOMETRY_2015/Digit_stereo_20150918/Flat_2.9/output_val/merged_NSW_atlasG4_Pt100_HITS_Digi_newGeometry_stereo_pitchokay_20150918.root

MM Trigger Simulation (SVN)

(Latest version, April 2016)

Instructions to get the MM trigger code running in lxplus or lxatut. There are a lot of tweaks needed to get it running in athena 20.1.0.2.

  • Setup
    • asetup 20.1.0.2,here,opt,slc6,64,gcc48
    • cmt co -r TrigT1NSWSimTools-00-00-13-branch Trigger/TrigT1/TrigT1NSWSimTools/Trigger/TrigT1/TrigT1NSWSimTools/
    • pkgco.py MuonGeoModelTest-04-00-05
    • cp Trigger/TrigT1/TrigT1NSWSimTools/jobOptions/NSWGeoSetup.py MuonSpectrometer/MuonGeoModelTest/share/
    • pkgco.py MuonReadoutGeometry-03-02-01
    • Change at the end (~lines 944-948) of MuonSpectrometer/MuonDetDescr/MuonReadoutGeometry/MuonReadoutGeometry/NSWgeometry.h
      • static const float MM_STEREO_ANGLE_ML1[4] = {-0.02618, 0.02618,0.,0.};
      • static const float MM_STEREO_ANGLE_ML2[4] = {0.,0.,-0.02618,0.02618};
      • //static const float MM_STEREO_ANGLE_ML1[4] = {0.,0.,0.,0.};
      • //static const float MM_STEREO_ANGLE_ML2[4] = {0.,0.,0.,0.};
    • cmt co DetectorDescription/IdDictParser/trunk DetectorDescription/IdDictParser
    • cmt co Trigger/TrigT1/TrigT1NSW/trunk Trigger/TrigT1/TrigT1NSW
    • In Trigger/TrigT1/TrigT1NSW/src/NSWL1Simulation.cxx, Change m_doMM = false to m_doMM = true.

  • Compile
    • setupWorkArea.py
    • cd $TestArea/WorkArea/cmt
    • cmt br cmt config
    • cmt br gmake
    • cd $TestArea

  • Run
    • cp Trigger/TrigT1/TrigT1NSWSimTools/jobOptions/NSWL1.py .
    • athena NSWL1.py

  • Preliminary output files (to be validated)
    • /data//maxi16s/atljphys/tsaito/TrigSimulation/TriggerOutput/20160407/output/merged_MC12.NSW.107206.atlasG4.01.0000.Pt.100.generator.HITS.root

Quick note about the directly ported code from Harvard: MMT_Finder.cxx and MMT_Fitter.cxx do all the machinery of the algorithm, while MMT_struct.cxx has all the hard-coded numbers about the strip positioning and wedge geometry.

-- TomoyukiSaito - 2016-04-08

Edit | Attach | Watch | Print version | History: r3 < r2 < r1 | Backlinks | Raw View | WYSIWYG | More topic actions
Topic revision: r3 - 2016-04-26 - TomoyukiSaito
 
    • Cern Search Icon Cern Search
    • TWiki Search Icon TWiki Search
    • Google Search Icon Google Search

    Main All webs login

This site is powered by the TWiki collaboration platform Powered by PerlCopyright &© 2008-2019 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding TWiki? Send feedback