Report problems to ATLAS LXR Team (with time and IP address indicated)

The LXR Cross Referencer
[ source navigation ]
[ diff markup ]
[ identifier search ]
[ general search ]
 
 
Architecture: [ linux ]
Version: [ head ] [ nightly ] [ GaudiDev ]
  Links to LXR source navigation pages for stable releases [ 12.*.* ]   [ 13.*.* ]   [ 14.*.* ]  
001 ###############################################################
002 #
003 #  iPatRecExample jobOptions file:
004 #  runs everything needed for iPatRec (Truth/InnerDetector etc)
005 #
006 #==============================================================
007     
008 #--------------------------------------------------------------
009 # Needed for Configurables
010 #--------------------------------------------------------------
011 from AthenaCommon.AppMgr import theApp
012 from AthenaCommon.AppMgr import ServiceMgr
013 from AthenaCommon.AppMgr import ToolSvc
014 from AthenaCommon.AlgSequence import AlgSequence
015 topSequence = AlgSequence()
016 
017 #--------------------------------------------------------------
018 # use auditors to write out a summary of the time spent
019 #--------------------------------------------------------------
020 from AthenaCommon.ConfigurableDb import getConfigurable
021 ServiceMgr.AuditorSvc   += getConfigurable("ChronoAuditor")()
022 #ServiceMgr.AuditorSvc   += getConfigurable("NameAuditor")()
023 theApp.AuditAlgorithms  =  True # FIXME necessary even to get CPU time
024 theApp.AuditServices    =  True
025 theApp.CreateSvc += ["FPEControlSvc"]
026 
027 #--------------------------------------------------------------
028 # set the number of events to be processed
029 #--------------------------------------------------------------
030 theApp.EvtMax   = 1000
031  
032 #--------------------------------------------------------------
033 # basic configuration: flags (a la RecExCommon)
034 #--------------------------------------------------------------
035 from RecExConfig.RecFlags import rec
036 from RecExConfig.RecAlgsFlags import recAlgs
037 doTruth         = True
038 makeESD         = True
039 makeNtuple      = False
040             
041 #--------------------------------------------------------------
042 # select DC3 data 
043 #--------------------------------------------------------------
044 test5       = False
045 test10      = False
046 test100     = False
047 test500     = False
048 testJ7      = False
049 testZmumu   = False
050 testDC7     = False
051 dataMuon10  = False
052 dataMuon100 = False
053 dataMuon300 = False
054 dataElec25  = False
055 dataH120bb  = False
056 dataTop500  = False
057 dataZmumu   = False
058 dataMbias   = False
059 
060 #test5       = True
061 #test10      = True
062 #test100     = True
063 #test500     = True
064 #testJ7      = True
065 testZmumu   = True
066 #testDC7     = True
067 
068 #--------------------------------------------------------------
069 # select layout 
070 #--------------------------------------------------------------
071 if test5:
072     DetDescrVersion = "ATLAS-CSC-01-00-00"
073 elif test10 or test100 or test500 or testJ7 or testZmumu:
074     DetDescrVersion = "ATLAS-CSC-01-02-00"     
075 elif testDC7:
076     DetDescrVersion = "ATLAS-DC3-07"
077 else:
078     DetDescrVersion = "ATLAS-DC3-02"
079 from AthenaCommon.GlobalFlags import globalflags
080 globalflags.DetDescrVersion = DetDescrVersion
081 print globalflags.DetDescrVersion()
082 
083 #--------------------------------------------------------------
084 # select detectors to be used for reconstruction
085 # uncomment appropriate line to turn pixel, SCT and/or TRT off)
086 #--------------------------------------------------------------
087 from AthenaCommon.DetFlags import DetFlags
088 DetFlags.readRDOPool.ID_setOn()
089 DetFlags.makeRIO.ID_setOn()
090 #DetFlags.pixel_setOff()
091 #DetFlags.SCT_setOff()
092 #DetFlags.TRT_setOff()
093 
094 #--------------------------------------------------------------
095 # can switch off HistorySvc when it's a bottleneck
096 #--------------------------------------------------------------
097 ServiceMgr.StoreGateSvc.ActivateHistory = False
098 #ServiceMgr.HistorySvc.Activate = False
099 
100 #--------------------------------------------------------------
101 # set appropriate configuration into AthenaCommon.globalflags
102 #--------------------------------------------------------------
103 globalflags.DetGeo      = 'atlas'
104 globalflags.DataSource  = 'geant4'
105 globalflags.InputFormat = 'pool'
106 
107 # setup POOL access in ATHENA
108 import AthenaPoolCnvSvc.ReadAthenaPool
109 #include( "AthenaPoolCnvSvc/ReadAthenaPool_jobOptions.py" )
110 #ServiceMgr.EventSelector.SkipEvents = 31
111 
112 #--------------------------------------------------------------
113 # Geomodel and MagneticField setup
114 #--------------------------------------------------------------
115 DetFlags.detdescr.ID_setOn()
116 from AtlasGeoModel import SetGeometryVersion
117 from AtlasGeoModel import GeoModelInit
118 include( "BFieldAth/BFieldAth_jobOptions.py" )
119 
120 # Beam Spot service
121 include( "InDetBeamSpotService/BeamCondSvc.py" )
122 
123 #--------------------------------------------------------------
124 # print final selection 
125 #--------------------------------------------------------------
126 globalflags.print_JobProperties()
127 DetFlags.Print()
128 
129 #--------------------------------------------------------------
130 # particle property service
131 #--------------------------------------------------------------
132 include( "PartPropSvc/PartPropSvc.py" )
133 
134 #--------------------------------------------------------------
135 # InnerDetector jobOptions fragment
136 #--------------------------------------------------------------
137 include( "iPatRecExample/InnerDetector_jobOptions.py" )
138 
139 #--------------------------------------------------------------
140 # iPatRec jobOptions fragment(s)
141 #--------------------------------------------------------------
142 include( "iPatRecExample/iPatRec_jobOptions.py" )
143 
144 from iPatRecAlgs.iPatRecAlgsConf import iPatShortTracks
145 #iPatShortTracks = iPatShortTracks( DetectorOption = 'pixel')  # pixel, sct, or trt
146 #topSequence += iPatShortTracks
147 
148 if makeNtuple:
149     from GaudiSvc.GaudiSvcConf import THistSvc
150     ServiceMgr += THistSvc()
151 
152     # this makes the iPatRec ntuple
153     from iPatRecAlgs.iPatRecAlgsConf import iPatNtuple
154     iPatNtuple = iPatNtuple()
155     topSequence += iPatNtuple
156     SeviceMgr.THistSvc.Output = [ "FILE1 DATAFILE='tuple1.root' OPT='NEW'" ]
157      
158 #--------------------------------------------------------------
159 # Legacy Conversion
160 #--------------------------------------------------------------
161 if makeESD:
162     include( "iPatRecExample/MakeESD_jobOptions.py" )
163 
164 #--------------------------------------------------------------
165 # jobO selection for data from POOL
166 #--------------------------------------------------------------
167 if globalflags.InputFormat == 'pool':
168     if DetFlags.readRDOPool.ID_on():
169         include ( "InDetEventCnvTools/InDetEventCnvTools_jobOptions.py" )
170         include ( "TrkEventAthenaPool/TrkEventAthenaPool_joboptions.py" )
171         include ( "InDetEventAthenaPool/InDetEventAthenaPool_joboptions.py" )
172     if DetFlags.readRDOPool.LAr_on() :
173         include( "LArAthenaPool/LArAthenaPool_joboptions.py" )
174     if DetFlags.readRDOPool.Tile_on():
175         include( "TileEventAthenaPool/TileEventAthenaPool_joboptions.py" )
176     if DetFlags.readRDOPool.Muon_on() or DetFlags.readRIOPool.Muon_on()  :
177         include( "MuonEventAthenaPool/MuonEventAthenaPool_joboptions.py" )
178     
179     #include( "EventAthenaPool/EventAthenaPool_joboptions.py" )
180     if doTruth:
181         include( "GeneratorObjectsAthenaPool/GeneratorObjectsAthenaPool_joboptions.py" )
182         #needed to read MuonEntryRecord
183         include( "G4SimAthenaPOOL/G4SimAthenaPOOL_joboptions.py" )
184 
185 #--------------------------------------------------------------
186 # file selection for data from POOL
187 #--------------------------------------------------------------
188 if testDC7:
189     ServiceMgr.EventSelector.InputCollections = [ "/afs/cern.ch/atlas/maxidisk/d38/ID_ValidationData/testIdeal_07.005200.T1_McAtNlo_Jimmy.digit.RDO.v12000201_tid002760._00002.pool.root.1"]
190 else:
191     # careful when ignoring tag - this is dangerous !
192     #GeoModelSvc.IgnoreTagDifference = True
193     t = []
194     filseq = range(50)
195     if test5:
196         filseq = [1,2,3,4,5,6,7,8,9,10]
197     elif test10:
198         filseq = [1,2,3,4,5,6,7,8,10,11,12,13,14,15,16,17,18,19,20,21]
199     elif test100:
200         filseq = [121,166]
201     elif test500:
202         filseq = [13]    
203     elif testJ7:
204         filseq = [4382]
205     elif testZmumu:
206         filseq = [1,2,3,5,6]
207 
208     for i in filseq:
209         if test5:
210             s = "castor:/castor/cern.ch/user/m/muonprod/1206/digit.RDO/mc12.007206.digit_atlas-csc-100.mu_p5GeV/mc12.csc100.007206.digit.muminus_p5GeV._"+"%05d"%(i)+".pool.root"
211         elif test10:
212             s = "castor:/castor/cern.ch/user/r/rtt/RTTdata/misal1_mc12.007211.singlepart_mu10.digit.RDO.v12000604/misal1_mc12.007211.singlepart_mu10.digit.RDO.v12000604_tid007981._"+"%05d"%(i)+".pool.root.1"
213         elif test100:
214             #s = "castor:/castor/cern.ch/user/r/rtt/RTTdata/misal1_mc12.007233.singlepart_mu100.digit.RDO.v12000604/misal1_mc12.007233.singlepart_mu100.digit.RDO.v12000604_tid007982._"+"%05d"%(i)+".pool.root.1"
215             s = "/afs/cern.ch/atlas/project/muon/data/csc12/misal1_csc11.007233.singlepart_mu100.digit.RDO.v12003101_tid003587._"+"%05d"%(i)+".pool.root"
216         elif test500:
217             s = "castor:/castor/cern.ch/user/r/rtt/RTTdata/misal1_mc12.007235.singlepart_mu500.digit.RDO.v12003107/misal1_mc12.007235.singlepart_mu500.digit.RDO.v12003107_tid004498._"+"%05d"%(i)+".pool.root.1"
218         elif testJ7:
219             s = "/afs/cern.ch/atlas/project/muon/data/csc12/misal1_csc11.005016.J7_pythia_jetjet.digit.RDO.v12003105_tid004302._"+"%05d"%(i)+".pool.root.1"
220         elif testZmumu:
221             s = "/afs/cern.ch/atlas/project/muon/data/csc12/misal1_csc11.005145.PythiaZmumu.digit.RDO.v12003101_tid003501._"+"%05d"%(i)+".pool.root.1"        
222         elif dataMbias:
223             s = "castor:/castor/cern.ch/atlas/csc/valiprod/sampleA/csc11/005001.pythia_minbias/digit/v11004201/csc11.005001.pythia_minbias.digit.RDO.v11004201._"+"%05d"%(i+1)+".pool.root"
224 
225 
226         t.append(s)
227 
228     #print "%s"%t
229     ServiceMgr.EventSelector.InputCollections = t
230 
231 #--------------------------------------------------------------
232 # Set output level threshold (ALL, VERBOSE, DEBUG, INFO, WARNING, ERROR, FATAL) 
233 #--------------------------------------------------------------
234 ServiceMgr.MessageSvc.OutputLevel  = INFO
235 ServiceMgr.MessageSvc.defaultLimit = 999999
236 
237 #--------------------------------------------------------------
238 # (re-)configure the iPatRec algorithms
239 #--------------------------------------------------------------
240 from iPatRecAlgs.iPatRecAlgsConf import iPatRec
241 iPatRec.maxSecondaryImpact           = 100.0*mm
242 iPatRec.minPt                        = 0.4*GeV
243 iPatRec.printLevel                   = 2   # 2 gives parameters of found tracks
244 iPatRec.TracksLocation               = "Tracks"
245 
246 from iPatRecAlgs.iPatRecAlgsConf import iPatStatistics
247 iPatStatistics.minPt             = iPatRec.minPt
248 
249 #iPatTrackTruthAssociator.OutputLevel      = DEBUG
250 
251 
252 #--------------------------------------------------------------
253 # ... and tools
254 #--------------------------------------------------------------
255 from iPatTrackFitter.iPatTrackFitterConf import TrackFitter
256 #ToolSvc += TrackFitter(name = 'TrackFitter', OutputLevel  = DEBUG)
257 from TrkiPatFitter.TrkiPatFitterConf import Trk__iPatFitter
258 #ToolSvc += Trk__iPatFitter(name = 'iPatFitter', OutputLevel  = VERBOSE)
259 from iPatTrackFinder.iPatTrackFinderConf import TrackFinder
260 ToolSvc += TrackFinder(name = 'TrackFinder',
261                        OutputLevel           = INFO,
262                        searchForPrimaries    = True,
263                        searchForSecondaries  = True,
264                        searchForVertexRegion = True)
265     
266 #--------------------------------------------------------------
267 # in case you want to set the magnetic field off 
268 #--------------------------------------------------------------
269 #TrajectorySvc = Service( "TrajectorySvc" )
270 # ServiceMgr.TrajectorySvc.magneticFieldOff = True
271 # ServiceMgr.TrajectorySvc.parametrizeSolenoid = False
272 
273 #--------------------------------------------------------------
274 # set tolerances for misaligned data
275 #--------------------------------------------------------------
276 # ToolSvc.AlignmentTolerances.sigmaPixelBarrel  = [0.100*mm, 0.010*mm, 0.100*mm]  #  r, rphi, z
277 # ToolSvc.AlignmentTolerances.sigmaPixelEndcap  = [0.100*mm, 0.010*mm, 0.100*mm]  #  r, rphi, z
278 # ToolSvc.AlignmentTolerances.sigmaSCT_Barrel   = [0.100*mm, 0.010*mm, 0.100*mm]  #  r, rphi, z
279 # ToolSvc.AlignmentTolerances.sigmaSCT_Endcap   = [0.100*mm, 0.010*mm, 0.100*mm]  #  r, rphi, z
280 
281 #==============================================================
282 # End of iPatRecExample job options file
283 #
284 ###############################################################
285
[ source navigation ] [ diff markup ] [ identifier search ] [ general search ]
Due to the LXR bug, the updates fail sometimes to remove references to deleted files. The Saturday's full rebuilds fix these problems
This page was automatically generated by the LXR engine. Valid HTML 4.01!