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

  • These instructions are based on release 20.20.X-VAL.
    • asetup 20.20.X-VAL,rel1,here

Particle Generation

  • See an example code: /afs/cern.ch/atlas/software/releases/20.20.5/AtlasSimulation/20.20.5/Generators/ParticleGun/examples/jobOption.ParticleGun_flatpt_2particle.py
  • Particle Generation files(Single muon with pT=100 GeV, Flat distribution in eta and phi plane): /afs/cern.ch/user/t/tsaito/public/SingleMuonSample/

Detector Simulation with NSW

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

  • Run
    • AtlasG4_tf.py --conditionsTag='OFLCOND-MC12b-SIM-00' --inputEVNTFile '/afs/cern.ch/user/t/tsaito/public/SingleMuonSample/MC12.NSW.107206.atlasG4.01.0000.Pt100.generator.pool.root' --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

MM Digitization (SVN)

  • 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/m/mkhader/public/mm_asetup2020/NSW.config.digi.py,/afs/cern.ch/user/m/mkhader/public/mm_asetup2020/NSWPRDValAlg.digi.py --preExec 'from AthenaCommon.DetFlags import DetFlags;DetFlags.ID_setOff();DetFlags.Calo_setOff();DetFlags.Muon_setOn();DetFlags.MDT_setOff();DetFlags.CSC_setOff();DetFlags.TGC_setOff();DetFlags.RPC_setOff();DetFlags.sTGC_setOff();DetFlags.Micromegas_setOn();DetFlags.Truth_setOn()' --inputHITSFile /afs/cern.ch/user/m/mkhader/public/mm_asetup2020/mm_asetup2020_samples/hits/MC12.NSW.107206.atlasG4.01.0001.Pt100.asetup2020.hits.pool.root --outputRDOFile test.rdos.pool.root --maxEvents 10

Pileup Digitization

  • Run
    • Digi_tf.py --inputHitsFile='%s' --outputRDOFile='%s' --preExec='HITtoRDO:from DigitizationFlags import digitizationFlags;digitizationFlags.overrideMetadata=["PhysicsList","SimLayout","MagneticField","SimulatedDetectors"];from DetFlags import DetFlags;DetFlags.Truth_setOn();DetFlags.BCM_setOn();DetFlags.ID_setOff();DetFlags.Calo_setOff();DetFlags.Muon_setOn(); DetFlags.MDT_setOn();DetFlags.CSC_setOn();DetFlags.TGC_setOn();DetFlags.RPC_setOn(); DetFlags.sTGC_setOn();DetFlags.Micromegas_setOn();DetFlags.digitize.MDT_setOff(); DetFlags.digitize.TGC_setOff();DetFlags.digitize.RPC_setOff();DetFlags.digitize.CSC_setOff(); DetFlags.digitize.Micromegas_setOn();DetFlags.digitize.sTGC_setOn();' --maxEvents=XX --skipEvents=XX --digiSeedOffset1=seed1 --digiSeedOffset2=seed2 --geometryVersion='ATLAS-R2-2015-03-01-00_VALIDATION' --conditionsTag='OFLCOND-MC12b-SIM-00' --preInclude='Digitization/ForceUseOfPileUpTools.py,/afs/cern.ch/user/t/tsaito/preInclude.PileUpBunchTrainsUpgrade_mu160cavern7_Myconfig.py' --LowPtMinbiasHitsFile='filename.root' --HighPtMinbiasHitsFile='filename.root' --inputCavernHitsFile='filename.root' --postExec="mm=job.StandardPileUpToolsAlg.PileUpTools[\\"MmDigitizationTool\\"];mm.RndmSvc=jobproperties.Digitization.rndmSvc.get_Value();mm.qThreshold=0.00000001;mm.SaveInternalHistos=True;mm.qThresholdForTrigger=2;mm.ValidationSetup=True;streamRDO.ItemList+=[\\"MuonSimDataCollection#MM_SDO\\",\\"MuonSimDataCollection#sTGC_SDO\\",\\"GenericMuonSimHitCollection#InputMicroMegasHits\\"];from GaudiSvcConf import THistSvc;ServiceMgr += THistSvc();ServiceMgr.THistSvc.Output=[\\"NSWPRDValAlg DATAFILE='NSWPRDValAlg.root' OPT='RECREATE'\\"]" --postInclude /afs/cern.ch/user/m/mkhader/public/mm_asetup2020/NSW.config.digi.py,/afs/cern.ch/user/m/mkhader/public/mm_asetup2020/NSWPRDValAlg.digi.py

MM Trigger Simulation (SVN)

  • 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

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.

Background Samples (produced by release 20.3)

  • Cavern background: /afs/cern.ch/work/i/iravinov/public/geant4/geant_hits/cavern_background/
  • Minimum bias: /afs/cern.ch/work/i/iravinov/public/geant4/geant_hits/minbias/

-- TomoyukiSaito - 2016-05-10

Edit | Attach | Watch | Print version | History: r11 < r10 < r9 < r8 < r7 | Backlinks | Raw View | WYSIWYG | More topic actions
Topic revision: r11 - 2017-06-27 - 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