(r33) CosmicSimulationSamplesRel14 < Main

Main Web>TWikiUsers>JamieBoyd>JamiesWiki>HowToSimulateCosmics>HowToSimulateCosmicIn1250>CosmicSimulationSamples>CosmicSimulationSamplesRel14 (revision 33) (raw view)~~EditAttachPDF~~
---+ Cosmic Simulation in rel 14: *Samples*
%TOC{title="Contents:"}%
---+ Overview
   * These samples are produced in 14.4.0 AtlasProduction on the batch system, although simulation step is also available with Simulation/SimuJobTransforms/scripts/csc_cosmics_sim_trf.py
   * What we used for these samples is now included in SimuJobTransforms-00-07-99, which needs tags G4CosmicFilter-00-00-12 and CosmicGenerator-00-00-34.
<!--- Tags for simulation: G4AtlasApps-00-02-87, CosmicGenerator-00-00-31,  InDetSimEvent-00-04-09, G4CosmicFilter-00-00-10, SimuJobTransforms-00-07-86 (and some later updates, see next bullet) -->
   * *The new simulation transform is included in 14.5.1.2,AtlasProduction,runtime*. Run without arguments for instructions - here is an example anyway:
<verbatim>  csc_cosmics_sim_trf.py outputHitsFile=sim.root maxEvents=10000 randomSeed=1234 geometryVersion=ATLAS-GEONF-04-00-00 CosmicFilterVolume=TRT_Barrel CosmicFilterVolume2=Pixel conditionsTag=OFLCOND-CSC-01-00-00 jobConfig=Cosmics.py
</verbatim> 
   * Tags for digitization: CommissionUtils-00-00-21, PixelDigitization-00-07-13-02, SiDigitization-03-05-04, SCT_Digitization-00-11-13 . jobOptions testdig.py attached.
   * The size of the ATLAS-GEO samples will be ~50k events when first stage finished.
   * To access the data use STAGE_SVCCLASS=atlasgroupdisk
   * Two filters: CosmicGenerator constrains the production vertex and direction of generator-level tracks, G4CosmicFilter requires simulated hits inside a certain volume - more details on the settings for the generator can be found in the attached jobOptions_ConfigCosmicProd.py.
   * Some limitations: 
      * Several questions lately: note that the mu+/mu- ratio in simulation (1.5) is not the same as in data. The simulation flux calculation uses this reference: A. Dar ,  Phys. Rev. Letters 51, 227 (1983)
      *  Timing is wrong for the muons, trigger will not work.

---++ List of geometries used:
   * *ATLAS-GEO-03-00-00* (aligned, Solenoid On, Toroids On)
   * *ATLAS-GEONF-04-00-00* (aligned, no magn field)
   * *ATLAS-GEO-04-00-00* (aligned, Solenoid off, Toroids On)

---++ List of simulation filter volumes used:
   * *Muon* (all of Atlas, geant volume is MUONQ2::MUONQ2)
   * *TRT barrel* (geant volume is TRT::BarrelOuterSupport)
   * *Pixel barrel* (has to have hits inside geant volumes TRT::BarrelOuterSupport _and_ Pixel::Pixel)
_Please note that all subsystems are switched on in all samples; the naming refers just to the filtering on the track direction_


---+ Samples
---++ Production with Muon Volume 
---+++++ ATLAS-GEO-03-00-00:
   * 5000 events at /castor/cern.ch/grid/atlas/atlasgroupdisk/proj-simcos/rel14/ATLAS-GEO-03-00-00/dig/dig.muon_B.70000x.root
where x=1-5

---+++++ ATLAS-GEONF-04-00-00:
   * so far only a few thousand test events. Pixel hit timing not correct for these. Can be found at ~gencomm/w0/cosmicSim_rel14/NoField/dig.muon*.root . 

---++ Production with TRT barrel Volume
---+++++ ATLAS-GEO-03-00-00:
   * 58666 events at /castor/cern.ch/grid/atlas/atlasgroupdisk/proj-simcos/rel14/ATLAS-GEO-03-00-00/dig/dig.GEO-03-TRTBarrel.*.root
   * The corresponding hits files can be found in the ATLAS-GEO-03-00-00/sim directory

---+++++ ATLAS-GEONF-04-00-00:
   * 50327 events at /castor/cern.ch/grid/atlas/atlasgroupdisk/proj-simcos/rel14/ATLAS-GEONF-04-00-00/dig.trtbarrel.200*.root
   * The corresponding hits files can be found in the ATLAS-GEONF-04-00-00/hits directory

---+++++ ATLAS-GEO-04-00-00:
   * 57756 events at /castor/cern.ch/grid/atlas/atlasgroupdisk/proj-simcos/rel14/ATLAS-GEO-04-00-00/dig
   * corresponding hits files can be seen in ATLAS-GEO-04-00-00/sim directory

---++ Production with pixel barrel Volume
---+++++ ATLAS-GEO-03-00-00:
   * Approx 46k events can be found at: /castor/cern.ch/grid/atlas/atlasgroupdisk/proj-simcos/rel14/ATLAS-GEO-03-00-00/dig/dig.pixelbarrel_B.400*.root
   * The corresponding hits files can be found in the ATLAS-GEO-03-00-00/sim directory

---+++++ ATLAS-GEONF-04-00-00:
   * Approx 48k events at /castor/cern.ch/grid/atlas/atlasgroupdisk/proj-simcos/rel14/ATLAS-GEONF-04-00-00/dig.pixelbarrel.500*.root
   * The corresponding hits files can be found in the ATLAS-GEONF-04-00-00/hits directory

---+++++ ATLAS-GEO-04-00-00:
   * 55368 events. Files are in /castor/cern.ch/grid/atlas/atlasgroupdisk/proj-simcos/rel14/ATLAS-GEO-04-00-00/dig/dig.pixelbarrel_B.600*.root
   * The corresponding hits files can be found in the ATLAS-GEO-04-00-00/sim directory

---++ How to run reconstruction:
   * Use release 14.5.0.2,AtlasProduction
   * Check out and compile <verbatim>cmt co -r RecExCommission-00-03-94 Reconstruction/RecExample/RecExCommission</verbatim>
   * get JO:
<verbatim>
get_files -jo Collisions_ATN.py
</verbatim>
   * edit this file with the following
<verbatim>
add the lines:
jobproperties.Beam.beamType.set_Value_and_Lock('cosmics')
rec.doTrigger.set_Value_and_Lock(False)
after the importing of the flags.
edit the lines:
<             globalflags.DetDescrVersion.set_Value_and_Lock('ATLAS-Comm-02-00-00')
<             athenaCommonFlags.PoolRDOInput.set_Value_and_Lock( ["/afs/cern.ch/atlas/maxidisk/d17/dig.trtbarrel_ma_B100.root"])
to point to the dataset you want and to set teh DetDescrVersion correctly for that datset.
(you may also have to edit the lines:
            jobproperties.BField.solenoidOn=True
            jobproperties.BField.barrelToroidOn=True
            jobproperties.BField.endcapToroidOn=True
to se the magnetic field configuration correctly).
for simulation made in rel 14 remove the line: 
           SpecialTileConfigHack=True
edit the number of events you want to run over:
athenaCommonFlags.EvtMax.set_Value_and_Lock(500)
</verbatim>
   * then run like:
<verbatim>
athena -c RTTsim=True Collisions_ATN.py
</verbatim>
   * This will run and produce ESD, CBNT. you can configure the job to also produce monitoring histograms, TAG_COMM, Jive xml by adding the appropriate rec flags.
(monitoring histograms fail if the run number changes during the job - which happens with old (<14) simulated cosmics)

---++ useful info
   * we can see how full our castor diskpool is by looking: [[http://atladcops.cern.ch:8000/mappings/groupdisk.log][here]] our limit is 5TB

-- Main.ElseLytken - 9 Dec 2008
   * [[%ATTACHURL%/testdig.py.txt][testdig.py.txt]]: example of digitization jobOptions


   * [[%ATTACHURL%/jobOptions_ConfigCosmicProd.py.txt][jobOptions_ConfigCosmicProd.py.txt]]: Cosmic generator settings