%CERTIFY%

Analysis of ATLAS Cosmic-Ray Data 2008/2009



This page is the entry point for the different analyses of the cosmic data that has been taken since autumn 2008.

Data and Monte Carlo samples

Data samples

A detailed description on how to access the real data from the different cosmic runs and reprocessing campaigns can be found on this TWiki page.

Monte Carlo samples

The following samples are available (you need to replace TYPE with ESD or AOD):
Filter Volume Magnet configuration Events Dataset
Pixel Both magnets ON 825k valid2.108865.CosSimPixVolSolOnTorOn.recon.TYPE.s531_d167_r676/
Pixel Both magnets OFF 995k valid2.108864.CosSimPixVolSolOffTorOff.recon.TYPE.s531_d167_r677/
Inner Detector Both magnets ON 328k valid2.108867.CosSimIDVolSolOnTorOn.recon.TYPE.s531_d167_r676/
Inner Detector Both magnets OFF 282k valid2.108866.CosSimIDVolSolOffTorOff.recon.TYPE.s531_d167_r677/
Muon Both magnets ON 465k valid2.108869.CosSimMuonVolSolOnTorOn.recon.TYPE.s531_d167_r676/
Muon Both magnets OFF 462k valid2.108868.CosSimMuonVolSolOffTorOff.recon.TYPE.s531_d167_r677/
A more description of the availabe Monte Carlo samples can be found on this TWiki page.

Example analysis jobs

Rerunning the monitoring on DPDs on a local machine

If you want to develop your algorithm you'll need to test it on a small sample before sending it to the Grid. The jobOptions below are based on the file RecExCommission/readesd.py and simply run the monitoring on the input DPD/ESD in release 15.0.0.3,AtlasProduction:

from RecExConfig.RecFlags  import rec
from AthenaCommon.BeamFlags import jobproperties
from AthenaCommon.AthenaCommonFlags import athenaCommonFlags
from AthenaCommon.GlobalFlags import globalflags

# input file(s)
athenaCommonFlags.PoolESDInput=['castor:/castor/cern.ch/grid/atlas/tzero/prod1/perm/data08_cosmag/physics_IDCosmic/0091900/data08_cosmag.00091900.physics_IDCosmic.recon.ESD.o4_f73/data08_cosmag.00091900.physics_IDCosmic.recon.ESD.o4_f73._lb0001._sfo01._0001.1']
athenaCommonFlags.EvtMax.set_Value_and_Lock(10)

#tell the job that we are reading ESD/DPD
globalflags.InputFormat.set_Value_and_Lock('pool')
rec.readESD.set_Value_and_Lock(True)
rec.readRDO.set_Value_and_Lock(False)

#autoconfiguration
rec.AutoConfiguration.append('FieldAndGeo')
rec.AutoConfiguration.append('BeamType')

# what shall be done in the job
rec.doHist.set_Value_and_Lock(True)
rec.doCBNT.set_Value_and_Lock(False)
rec.doWriteTAGCOM.set_Value_and_Lock(False)
rec.doWriteTAG.set_Value_and_Lock(False)   
rec.doWriteAOD.set_Value_and_Lock(False)
rec.doAOD.set_Value_and_Lock(False)
rec.doMonitoring.set_Value_and_Lock(True)
rec.doESD.set_Value_and_Lock(False)
rec.doWriteESD.set_Value_and_Lock(False)
rec.doPerfMon=False


include("RecExCommission/RecExCommission.py")                
from AthenaMonitoring.DQMonFlags import DQMonFlags
DQMonFlags.monType='BS_ESD'
DQMonFlags.histogramFile.set_Value_and_Lock( "myMonitoringESD.root" )
include ("RecExCommon/RecExCommon_topOptions.py")
   

As a first step make sure that the above jobOptions work in your setup before you start adding your own algorithms at the end of the above jobOptions. This helps in disentangling the errors from the general setup and from your algorithm.

Rerunning the monitoring on DPDs on the Grid

The following example shows you how to rerun the monitoring from a DPD or ESD on the Grid. You can use it as a starting point for your analyses by adding the corresponding joOptions at the end of this job and disabling the monitoring (if needed). It is doing the same job as the example above with the difference that this time the job is run on the Grid. .
  • Release: 15.0.0.3,AtlasProduction
  • jobOptions
    from RecExConfig.RecFlags  import rec
    from AthenaCommon.BeamFlags import jobproperties
    from AthenaCommon.AthenaCommonFlags import athenaCommonFlags
    from AthenaCommon.GlobalFlags import globalflags
    
    # various configuration special to commissioning, needed for all jobs
    jobproperties.Beam.beamType.set_Value_and_Lock("cosmics")
    globalflags.DetDescrVersion.set_Value_and_Lock('ATLAS-GEO-03-00-00')
    
    
    ##
    ## The following try... except is only needed for ganga!
    ##
    
    try:
        import string
        inp=open("input.py").readlines()
        print inp
        print
        for i in inp:
            print i
            if i.find("ServiceMgr.EventSelector.InputCollections")>=0:
                print "--> Found match"
                t=i.split("=")
                print t[1].strip()
                exec("athenaCommonFlags.PoolESDInput="+t[1].strip())
                print "Result"
                print athenaCommonFlags.PoolESDInput
    except Exception, inst:
        print inst
        pass
    
    ##
    ## End of hack for ganga
    ##
    
    
    athenaCommonFlags.EvtMax.set_Value_and_Lock(-1)
    
    
    from AthenaCommon.BFieldFlags import jobproperties
    jobproperties.BField.solenoidOn.set_Value_and_Lock(True)
    jobproperties.BField.barrelToroidOn.set_Value_and_Lock(True)
    jobproperties.BField.endcapToroidOn.set_Value_and_Lock(True)
    
    rec.doJetMissingETTag.set_Value_and_Lock(False)
    from JetRec.JetRecFlags import jetFlags
    jetFlags.Enabled.set_Value_and_Lock(False)
    from RecExConfig.RecAlgsFlags  import recAlgs
    recAlgs.doMissingET.set_Value_and_Lock(False)
    rec.doTau.set_Value_and_Lock(False)
    recAlgs.doMuonSpShower.set_Value_and_Lock(False)
    
    
    
    rec.doHist.set_Value_and_Lock(True)
    rec.doCBNT.set_Value_and_Lock(False)
    rec.doWriteTAGCOM.set_Value_and_Lock(False)
    rec.doWriteTAG.set_Value_and_Lock(False)   
    rec.doWriteAOD.set_Value_and_Lock(False)
    rec.doAOD.set_Value_and_Lock(False)
    rec.doMonitoring.set_Value_and_Lock(True)
    
    rec.readESD.set_Value_and_Lock(True)
    rec.doESD.set_Value_and_Lock(False)
    rec.doWriteESD.set_Value_and_Lock(False)
    globalflags.InputFormat.set_Value_and_Lock('pool')
    
    rec.doPerfMon=False
    
    # the main jobOpt
    include("RecExCommission/RecExCommission.py")                
    
    # Tracking ntuple: it only works if ID reconstruction is re-ran (FIXME)
    #ATLASCosmicFlags.doTrkNtuple = False
    
    from AthenaMonitoring.DQMonFlags import DQMonFlags
    DQMonFlags.monType='BS_ESD'
    DQMonFlags.histogramFile.set_Value_and_Lock( "myMonitoringESD.root" )
    include ("RecExCommon/RecExCommon_topOptions.py")
    include.block("input.py")
    
    # use closest DB replica
    
    from AthenaCommon.AppMgr import ServiceMgr
    from PoolSvc.PoolSvcConf import PoolSvc
    ServiceMgr+=PoolSvc(SortReplicas=True)
    from DBReplicaSvc.DBReplicaSvcConf import DBReplicaSvc
    ServiceMgr+=DBReplicaSvc(UseCOOLSQLite=False)
    
         
    The important difference to note are
    • The statements specifying the input files have been removed. In principle ganga should take care of adding those back, however at the moment this does not work and therefore a patch as shown above had to be put it
    • On the Grid you have to specify for now the geometry and the field settings as the autoconfiguration doesn't seem to work properly
    • It is also needed to change the job to talk to the closest DB replica as the sqlite files do not contain enough information to fully process real data (last lines of the jobOptions above)
  • Ganga commands: The following lines will submit the job to the Grid using Ganga
    j = Job() 
    j.application=Athena() 
    j.application.option_file=['readesd_grid.py'] 
    j.application.prepare(athena_compile=False) 
    j.splitter=DQ2JobSplitter() 
    j.splitter.numsubjobs=2 
    j.inputdata=DQ2Dataset() 
    j.inputdata.dataset="data08_cosmag.00091391.physics_IDCosmic.merge.DPD_SCTCOMM.o4_r653_p26/" 
    j.outputdata=DQ2OutputDataset() 
    j.outputdata.outputdata=['myMonitoringESD.root'] 
    j.backend=LCG() 
    j.backend.requirements.cloud='DE' 
    j.submit() 
         

Contact persons

For questions regarding the cosmic reconstruction and analyses carried out in the different systems/groups please contact the corresponding persons below:

Links to the analyses TWiki pages of the different systems/groups

Uncovered topics


Major updates:
-- ChristianSchmitt - 18 May 2009

%RESPONSIBLE% Christian Schmitt
%REVIEW% Never reviewed

Edit | Attach | Watch | Print version | History: r7 < r6 < r5 < r4 < r3 | Backlinks | Raw View | WYSIWYG | More topic actions
Topic revision: r7 - 2009-05-25 - ChristianSchmitt
 
    • Cern Search Icon Cern Search
    • TWiki Search Icon TWiki Search
    • Google Search Icon Google Search

    Sandbox All webs login

This site is powered by the TWiki collaboration platform Powered by PerlCopyright & 2008-2020 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
or Ideas, requests, problems regarding TWiki? use Discourse or Send feedback