Difference: HIPMC15 (1 vs. 3)

Revision 32016-03-21 - AnthonyEricLionti

Line: 1 to 1
 
META TOPICPARENT name="AnthonyEricLionti"

Highly Ionizing Particles - MC15b full chain recipe

Line: 25 to 25
 
  • setupATLAS
  • asetup 19.2.3.7,AtlasProduction,here
  • pkgco.py InDetCondFolders -00-00-07
Changed:
<
<
>
>
 If you want to simulate samples with delta-rays fixed you should also modify the G4AtlasApps package :
  • sed -i '153s/0.05/30.0/' Simulation/G4Atlas/G4AtlasApps/python/atlas_idet.py
  • sed -i '154s/0.05/30.0/' Simulation/G4Atlas/G4AtlasApps/python/atlas_idet.py

Revision 22016-03-14 - AnthonyEricLionti

Line: 1 to 1
 
META TOPICPARENT name="AnthonyEricLionti"

Highly Ionizing Particles - MC15b full chain recipe

Line: 6 to 6
 

Event Generation

The first step is to generate some EVNT files. There are two different kind of events. The single particles ones generated with ParticleGun and the Drell-Yan ones generated with MadGraph+Pythia. In a new folder setup an athena release by doing :

Added:
>
>
For the single particles samples :
 
  • setupATLAS
  • asetup 19.2.3.7,AtlasProduction,here
Changed:
<
<
Then simply run athena with the corresponding jobOptions file here :
  • athena jobOption.ParticleGun_corrhist_el.py
Or here :
  • athena jobOption.ParticleGun_corrhist_mu.py
Make sure that the corresponding eta-pT histograms attached to this page are in the same folder.
>
>
For the Drell-Yan samples :
  • setupATLAS
  • asetup 19.2.4.10.2,MCProd,here

Then download the file here and remove the .txt extension. Finaly run athena with the corresponding jobOptions file here that you can edit before to change the sample mass and charge:

  • Generate_tf.py --firstEvent 1 --maxEvents 10 --runNumber 000001 --randomSeed $RANDOM --ecmEnergy 13000 --jobConfig MC15.000001.MGPy8EG_A14N23LO_DYMonopole.py --outputEVNTFile drellYan.EVNT.pool.root
Or for the single particles download the file here and remove the .txt extension, then run athena with the corresponding jobOptions file here :
  • Generate_tf.py --firstEvent 1 --maxEvents 10 --runNumber 000001 --randomSeed $RANDOM --ecmEnergy 13000 --jobConfig MC15.000001.ParticleGun_single_magneticMonopole.py --outputEVNTFile singleParticle.EVNT.pool.root
 

Simulation

From the EVNT file produced previously we will run geant 4. In a new folder :

  • setupATLAS
Changed:
<
<
  • asetup 19.2.1.5,AtlasProduction,here
>
>
  • asetup 19.2.3.7,AtlasProduction,here
  If you want to simulate samples with delta-rays fixed you should also modify the G4AtlasApps package :
Line: 31 to 35
 
  • cmt br cmt config
  • source setup.sh
  • cmt br cmt make
Changed:
<
<
Finally you can either submit the simulation to the GRID thanks to the script attached here that you would have to adapt a bit or locally by doing :
>
>
Then run the simulation :
 
  • cd ../run
Changed:
<
<
  • AtlasG4 _tf.py --inputEvgenFile SingleElectrons.2013.pool.root --outputHitsFile SingleElectrons.2013.HITS.pool.root --maxEvents 30000 --randomSeed 12345 --geometryVersion ATLAS-R2-2015-02-01-00_VALIDATION --conditionsTag OFLCOND-RUN12-SDR-21 --DataRunNumber 222222
>
>
  • Sim_tf.py --inputEvgenFile singleParticle.EVNT.pool.root --outputHitsFile singleParticle.HITS.pool.root --maxEvents 10 --DBRelease current --DataRunNumber 222525 --conditionsTag OFLCOND-RUN12-SDR-19 --geometryVersion ATLAS-R2-2015-03-01-00_VALIDATION --physicsList FTFP_BERT --postInclude PyJobTransforms /UseFrontier.py --preInclude "EVNTtoHITS:SimulationJobOptions/preInclude.BeamPipeKill.py,SimulationJobOptions/preInclude.FrozenShowersFCalOnly.py" --simulator MC12G4 --truthStrategy MC12
 In each case you would need to be sure that the previously generated EVNT files are in the run directory.
Changed:
<
<

Digitization

>
>

Digitization an reconstruction

  In a clean environment in an empty directory do the following:
  • setupATLAS
Changed:
<
<
  • asetup 20.1.5.12,AtlasProduction,here
  • pkgco.py TRT_Digitization-01-01-03
>
>
  • asetup 20.1.5.10.1,AtlasProd1,here
  • pkgco.py TRT_Digitization-01-01-07
 
  • pkgco.py TRT_PAI_Process-00-01-01
  • setupWorkArea.py
  • cd WorkArea /cmt
Line: 49 to 53
 
  • source setup.sh
  • cmt br cmt make
  • cd ../run
Changed:
<
<
Now you are ready to run the digitization either locally or on the GRID with the script here :
  • Reco_tf.py --ignoreErrors True --maxEvents 30000 --inputHITSFile SingleElectrons.2013.HITS.pool.root --outputRDOFile SingleElectrons.2013.RDO.pool.root --conditionsTag OFLCOND-RUN12-SDR-21 --geometryVersion ATLAS-R2-2015-02-01-00_VALIDATION --DataRunNumber 222222 --postInclude 'TRT_Digitization/postInclude.OverrideTRTparameters.py,RecJobTransforms/UseFrontier.py' --autoConfiguration 'everything'
>
>
Now you are ready to run the digitization :
  • Reco_tf.py --maxEvents 10 --inputHITSFile singleParticle.HITS.pool.root --outputESDFile singleParticle.ESD.pool.root --conditionsTag OFLCOND-RUN12-SDR-31 --geometryVersion ATLAS-R2-2015-03-01-00 --inputHighPtMinbiasHitsFile highPt.pool.root --inputLowPtMinbiasHitsFile lowPt.pool.root --numberOfCavernBkg 0 --numberOfHighPtMinBias 0.12268057 --numberOfLowPtMinBias 39.8773194 --pileupFinalBunch 6 --postExec "all:CfgMgr.MessageSvc().setError+=[\"HepMcParticleLink\"];svcMgr.IOVDbSvc.CacheTime=100000" "HITtoRDO:job.StandardPileUpToolsAlg.PileUpTools[\"MergeMcEventCollTool\"].OnlySaveSignalTruth=True;job.StandardPileUpToolsAlg.PileUpTools[\"MdtDigitizationTool\"].LastXing=150" "RDOtoRDOTrigger:from AlgSequence import AlgSequence;AlgSequence().LVL1TGCTrigger.TILEMU=True;from AppMgr import ServiceMgr;import MuonRPC _CablingConfig;ServiceMgr.MuonRPC_CablingSvc.RPCTriggerRoadsfromCool=False" --postInclude 'TRT_Digitization/postInclude.OverrideTRTparameters.py,RecJobTransforms/UseFrontier.py' --preExec "all:rec.Commissioning.set_Value_and_Lock(True);from BeamFlags import jobproperties;jobproperties.Beam.numberOfCollisions.set_Value_and_Lock(20.0);from 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)" "RAWtoESD:from CaloCellFlags import jobproperties;jobproperties.CaloCellFlags.doLArCellEmMisCalib=False" "ESDtoAOD:TriggerFlags.AODEDMSet=\"AODSLIM\"" "HITtoRDO:RunDMC_testing_configuration=True" --preInclude "HITtoRDO:Digitization/ForceUseOfPileUpTools.py,SimulationJobOptions/preInclude.PileUpBunchTrainsMC15_2015_25ns_Config1.py,RunDependentSimData/configLumi_run222526_v1.py,RunDependentSimData/configCommon.py" --autoConfiguration 'everything' --steering "RAWtoESD:in-RDO,in+RDO_TRIG,in-BS" --triggerConfig "MCRECO:DBF:TRIGGERDBMC:2013,7,11"
 If you want to modify the digitization parameters like the high-threshold values or the TR efficiencies, this should be done at this step by modifying the file TRT_Digitization/postInclude.OverrideTRTparameters.py before running the digitization.
Changed:
<
<
The default values of this digitization package correspond to the delta-rays suppressed tune.
>
>
The default values of this digitization package correspond to the delta-rays non suppressed tune.
 
Changed:
<
<
It is highly adviced when you have to generate samples with different threshold and TR efficiencies for both muons and electrons with or without the delta-rays fix to modify the script to automatically submit all the jobs.
>
>
If you want to use delta-rays suppressed HITS files, then you have to run the digitization without pile-up since there is no validated min bias samples yet :
  • Reco_tf.py --maxEvents 10 --inputHITSFile singleParticle.HITS.pool.root --outputESDFile singleParticle.ESD.pool.root --conditionsTag OFLCOND-RUN12-SDR-31 --geometryVersion ATLAS-R2-2015-03-01-00 --postInclude 'TRT_Digitization/postInclude.OverrideTRTparameters.py,RecJobTransforms/UseFrontier.py' --autoConfiguration 'everything'
 

Reconstruction and Ntuples

Changed:
<
<
In a clean environment in a new directory download the file here, then do :
  • tar -xvzf ntup.tar.gz
  • setupATLAS
  • asetup 20.1.5.12,AtlasProduction,here
  • setupWorkArea.py
  • cd WorkArea /cmt
  • cmt br cmt config
  • source setup.sh
  • cmt br cmt make
  • cd ../..
Now after modifying the input and path in either the job options or the script you can run either locally :
  • athena jobOptions_RecExCommon.py
Or on the GRID :
  • ./pathena.sh
>
>
In a clean environment in a new directory download the ntuple maker code here : co svn+ssh://svn.cern.ch/reps/atlasphys-exo/Physics/Exotic/UEH/HighlyIonising/Run2/NtupleCode/
  -- AnthonyEricLionti - 2016-03-14

Revision 12016-03-14 - AnthonyEricLionti

Line: 1 to 1
Added:
>
>
META TOPICPARENT name="AnthonyEricLionti"

Highly Ionizing Particles - MC15b full chain recipe

This TWiki page is designed to produce HIP MC15b private samples.

Event Generation

The first step is to generate some EVNT files. There are two different kind of events. The single particles ones generated with ParticleGun and the Drell-Yan ones generated with MadGraph+Pythia. In a new folder setup an athena release by doing :

  • setupATLAS
  • asetup 19.2.3.7,AtlasProduction,here

Then simply run athena with the corresponding jobOptions file here :

  • athena jobOption.ParticleGun_corrhist_el.py
Or here :
  • athena jobOption.ParticleGun_corrhist_mu.py
Make sure that the corresponding eta-pT histograms attached to this page are in the same folder.

Simulation

From the EVNT file produced previously we will run geant 4. In a new folder :

If you want to simulate samples with delta-rays fixed you should also modify the G4AtlasApps package :
  • sed -i '153s/0.05/30.0/' Simulation/G4Atlas/G4AtlasApps/python/atlas_idet.py
  • sed -i '154s/0.05/30.0/' Simulation/G4Atlas/G4AtlasApps/python/atlas_idet.py
Then you need to compile :
  • setupWorkArea.py
  • cd WorkArea /cmt
  • cmt br cmt config
  • source setup.sh
  • cmt br cmt make
Finally you can either submit the simulation to the GRID thanks to the script attached here that you would have to adapt a bit or locally by doing :
  • cd ../run
  • AtlasG4 _tf.py --inputEvgenFile SingleElectrons.2013.pool.root --outputHitsFile SingleElectrons.2013.HITS.pool.root --maxEvents 30000 --randomSeed 12345 --geometryVersion ATLAS-R2-2015-02-01-00_VALIDATION --conditionsTag OFLCOND-RUN12-SDR-21 --DataRunNumber 222222
In each case you would need to be sure that the previously generated EVNT files are in the run directory.

Digitization

In a clean environment in an empty directory do the following:

  • setupATLAS
  • asetup 20.1.5.12,AtlasProduction,here
  • pkgco.py TRT_Digitization-01-01-03
  • pkgco.py TRT_PAI_Process-00-01-01
  • setupWorkArea.py
  • cd WorkArea /cmt
  • cmt br cmt config
  • source setup.sh
  • cmt br cmt make
  • cd ../run
Now you are ready to run the digitization either locally or on the GRID with the script here :
  • Reco_tf.py --ignoreErrors True --maxEvents 30000 --inputHITSFile SingleElectrons.2013.HITS.pool.root --outputRDOFile SingleElectrons.2013.RDO.pool.root --conditionsTag OFLCOND-RUN12-SDR-21 --geometryVersion ATLAS-R2-2015-02-01-00_VALIDATION --DataRunNumber 222222 --postInclude 'TRT_Digitization/postInclude.OverrideTRTparameters.py,RecJobTransforms/UseFrontier.py' --autoConfiguration 'everything'
If you want to modify the digitization parameters like the high-threshold values or the TR efficiencies, this should be done at this step by modifying the file TRT_Digitization/postInclude.OverrideTRTparameters.py before running the digitization. The default values of this digitization package correspond to the delta-rays suppressed tune.

It is highly adviced when you have to generate samples with different threshold and TR efficiencies for both muons and electrons with or without the delta-rays fix to modify the script to automatically submit all the jobs.

Reconstruction and Ntuples

In a clean environment in a new directory download the file here, then do :

  • tar -xvzf ntup.tar.gz
  • setupATLAS
  • asetup 20.1.5.12,AtlasProduction,here
  • setupWorkArea.py
  • cd WorkArea /cmt
  • cmt br cmt config
  • source setup.sh
  • cmt br cmt make
  • cd ../..
Now after modifying the input and path in either the job options or the script you can run either locally :
  • athena jobOptions_RecExCommon.py
Or on the GRID :
  • ./pathena.sh

-- AnthonyEricLionti - 2016-03-14

META FILEATTACHMENT attachment="MadGraphControl_MGPy8EG_A14N23LO_DYMonopole.py.txt" attr="" comment="" date="1457956038" name="MadGraphControl_MGPy8EG_A14N23LO_DYMonopole.py.txt" path="MadGraphControl_MGPy8EG_A14N23LO_DYMonopole.py.txt" size="4532" user="anlionti" version="1"
META FILEATTACHMENT attachment="MC15.000001.MGPy8EG_A14N23LO_DYMonopole.py.txt" attr="" comment="" date="1457956038" name="MC15.000001.MGPy8EG_A14N23LO_DYMonopole.py.txt" path="MC15.000001.MGPy8EG_A14N23LO_DYMonopole.py.txt" size="119" user="anlionti" version="1"
META FILEATTACHMENT attachment="MC15.000001.ParticleGun_single_magneticMonopole.py.txt" attr="" comment="" date="1457956038" name="MC15.000001.ParticleGun_single_magneticMonopole.py.txt" path="MC15.000001.ParticleGun_single_magneticMonopole.py.txt" size="70" user="anlionti" version="1"
META FILEATTACHMENT attachment="ParticleGun_SingleMonopole.py.txt" attr="" comment="" date="1457956038" name="ParticleGun_SingleMonopole.py.txt" path="ParticleGun_SingleMonopole.py.txt" size="2426" user="anlionti" version="1"
 
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