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

• Setup
  • cmt co NSW_Sim MuonSpectrometer/MuonG4/NSW_Sim
  • cmt co MuonReadoutGeometry MuonSpectrometer/MuonDetDescr/MuonReadoutGeometry

• 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)

• Setup
  • cmt co -r MM_Digitization-00-05-27 MuonSpectrometer/MuonDigitization/MM_Digitization
  • cmt co -r MuonPRDTest-01-01-10 MuonSpectrometer/MuonValidation/MuonPRDTest
  • cmt co -r MuonDigitContainer-01-02-23 MuonSpectrometer/MuonDigitContainer
  • cmt co -r MuonReadoutGeometry-03-03-11 MuonSpectrometer/MuonDetDescr/MuonReadoutGeometry
  • cmt co -r MuonAGDDBase-00-00-17 MuonSpectrometer/MuonDetDescr/MuonAGDDBase
  • cmt co -r MuonAGDDDescription-00-00-11 MuonSpectrometer/MuonDetDescr/MuonAGDDDescription

• 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)

• Setup
• cmt co TrigT1NSWSimTools Trigger/TrigT1/TrigT1NSWSimTools/
• cmt co Trigger/TrigT1/TrigT1NSW/trunk Trigger/TrigT1/TrigT1NSW
• cmt co MuonReadoutGeometry MuonSpectrometer/MuonDetDescr/MuonReadoutGeometry
• cmt co MuonAGDDBase MuonSpectrometer/MuonDetDescr/MuonAGDDBase
• cmt co MuonAGDDDescription MuonSpectrometer/MuonDetDescr/MuonAGDDDescription
• cmt co NSW_Sim MuonSpectrometer/MuonG4/NSW_Sim (Added, 28/02/2017)
• in Trigger/TrigT1/TrigT1NSW/src/NSWL1Simulation.cxx update
  • 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

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