ATLAS Memo

EVNTの生成

  • mc15_13TeV.424016.ParticleGun_single_mu_Pt200to2000.evgen.EVNT.e5771
    • https://bigpanda.cern.ch/task/10704665/
    • https://ami.in2p3.fr/app/?subapp=amiTags_show&userdata=e5771
      asetup AtlasProduction,19.2.5.17
      Generate_tf.py --evgenJobOpts 'MC15JobOpts-00-05-77_v0.tar.gz' --ecmEnergy '13000'
      
      --AMITag=e5771
      --ecmEnergy=13000
      --evgenJobOpts=MC15JobOpts-00-05-77_v0.tar.gz
      --jobConfig=MC15JobOptions/MC15.424016.ParticleGun_single_mu_Pt200to2000.py
      --runNumber=424016
      
      Generate_tf.py --AMITag e5771 --ecmEnergy "13000" --evgenJobOpts MC15JobOpts-00-05-77_v0.tar.gz --jobConfig MC15JobOptions/MC15.424016.ParticleGun_single_mu_Pt200to2000.py --runNumber 424016 --outputEVNTFile EVNT.pool.root --firstEvent 1 --randomSeed 1 --maxEvents 10000
    • MC15JobOptions/MC15.424016.ParticleGun_single_mu_Pt200to2000.py
      evgenConfig.description = "Single mu with flat phi, eta in [-3.0, 3.0], and uniform pT in [200,2000] GeV"
      evgenConfig.keywords = ["singleParticle", "muon"]
      
      include("MC15JobOptions/ParticleGun_Common.py")
      
      import ParticleGun as PG
      genSeq.ParticleGun.sampler.pid = (-13, 13)
      genSeq.ParticleGun.sampler.mom = PG.PtEtaMPhiSampler(pt=[200000, 2000000], eta=[-3.0, 3.0])

以上を参考にして、

  • MC15.500000.ParticleGun_single_piminus_Pt1000.py
    evgenConfig.description = "Single piminus with flat phi, eta in [-1.25, 1.25], and pT of 1 GeV"
    evgenConfig.keywords = ["singleParticle", "pi-"]
    
    include("MC15JobOptions/ParticleGun_Common.py")
    
    import ParticleGun as PG
    genSeq.ParticleGun.sampler.n = 100
    genSeq.ParticleGun.sampler.pid = -211
    genSeq.ParticleGun.sampler.mom = PG.PtEtaMPhiSampler(pt=1000, eta=[-1.25, 1.25])

を作った。1イベントに負パイオンを100個。そして

asetup AtlasProduction,19.2.5.17
Generate_tf.py --ecmEnergy "13000" --jobConfig MC15.500000.ParticleGun_single_piminus_Pt1000.py --runNumber 500000 --outputEVNTFile EVNT.pool.root --firstEvent 1 --randomSeed 1 --maxEvents 10000 

を走らせることで、single pionを1万イベント作れた。

EVNT->DAOD_TRUTH0

DAOD_TRUTH0の作り方。

asetup AthDerivation 21.2.44.0
Reco_tf.py --AMI p3655 --inputEVNTFile EVNT.pool.root --outputDAODFile pool.root --reductionConf TRUTH0

EVNT->HITS

  • mc16_13TeV.424016.ParticleGun_single_mu_Pt200to2000.simul.HITS.e5771_s3112

を参考にしてMC16の設定でGeant4を走らせた。

asetup AtlasOffline 21.0.15
Sim_tf.py --AMI s3126 --inputEVNTFile EVNT.pool.root --outputHITSFile HITS.pool.root --randomSeed 0

HITS->RDO->HIST without pileup

  • mc16_13TeV.424016.ParticleGun_single_mu_Pt200to2000.recon.ESD.e5771_s3112_r10201

上記のサンプルはパイルアップありだが、ないものが欲しいのと、出力ファイルはHISTだけで十分なので、MC16dの設定を書き換えた。

asetup Athena 21.0.53

Reco_tf.py \
--digiSteeringConf 'StandardSignalOnlyTruth' \
--conditionsTag 'default:OFLCOND-MC16-SDR-20' \
--valid 'True' \
--autoConfiguration 'everything' \
--steering 'doRDO_TRIG' \
--postExec \
'all:CfgMgr.MessageSvc().setError+=["HepMcParticleLink"]' \
"ESDtoAOD:fixedAttrib=[s if 'CONTAINER_SPLITLEVEL = \"99\"' not in s else '' for s in svcMgr.AthenaPoolCnvSvc.PoolAttributes];svcMgr.AthenaPoolCnvSvc.PoolAttributes=fixedAttrib" \
--postInclude 'default:PyJobTransforms/UseFrontier.py' \
--preExec 'all:rec.Commissioning.set_Value_and_Lock(True);from AthenaCommon.BeamFlags import jobproperties;jobproperties.Beam.numberOfCollisions.set_Value_and_Lock(20.0);from LArROD.LArRODFlags import larRODFlags;larRODFlags.NumberOfCollisions.set_Value_and_Lock(20);larRODFlags.nSamples.set_Value_and_Lock(4);larRODFlags.doOFCPileupOptimization.set_Value_and_Lock(True);larRODFlags.firstSample.set_Value_and_Lock(0);larRODFlags.useHighestGainAutoCorr.set_Value_and_Lock(True); from LArDigitization.LArDigitizationFlags import jobproperties;jobproperties.LArDigitizationFlags.useEmecIwHighGain.set_Value_and_Lock(False)' \
'ESDtoAOD:from TriggerJobOpts.TriggerFlags import TriggerFlags;TriggerFlags.AODEDMSet.set_Value_and_Lock("AODSLIM");' \
'all:from ParticleBuilderOptions.AODFlags import AODFlags; AODFlags.ThinInDetForwardTrackParticles.set_Value_and_Lock(True)' \
--triggerConfig 'RDOtoRDOTrigger=MCRECO:DBF:TRIGGERDBMC:2179,51,207' \
--geometryVersion 'default:ATLAS-R2-2016-01-00-01' \
--numberOfCavernBkg '0' \
--inputHITSFile HITS.pool.root \
--outputHISTFile HIST.root \
--maxEvents 10 \
--digiSeedOffset1 1 \
--digiSeedOffset2 1 \
--jobNumber 1 \
--runNumber 300000

jobNumberは1以上でなければいけない。

SCTLorentzMonitoringToolの変更

https://gitlab.cern.ch/oda/athena/commits/my-21.0.53

リンク

やりたいこと

https://cds.cern.ch/record/1596937/files/ATL-INDET-INT-2013-007.pdf

SCTのdigitizationは場所には依存するようにはしていないはずなので、HVも温度も一定にして、pionのsingle particle MCを走らせて、ローレンツ角が場所に依るか見てみる。

依存性があるなら、トラック再構成の位置依存性があって、バイアスしているということになる。

そうしたら、Geant4での通過位置と再構成された通過位置を比べてみる。

結果


トップ   編集 凍結解除 差分 バックアップ 添付 複製 名前変更 リロード   新規 一覧 単語検索 最終更新   ヘルプ   最終更新のRSS
Last-modified: 2020-01-06 (月) 12:13:17 (214d)