Cherenkov ring finding package

The package RichRecTemplateRings is in the Rich/RichRecTemplateRings subdirectory directory of the REC packge which is part of the LHCb Brunel software. This RichRecTemplateRings is part of the Rich reconstruction packages and gets its input data from the RichRecSegments and RichRecHits available from the RICH software. It then creates the rings made by tracksegments in each of the the three radiators of the RICH system. The output is stored in to the standard RichRecRing class.

Using this package

One can activate the running of this package from the standard Brunel options file.
   The running of this package is configured in the  file 
    ../Rich/RichRecTemplateRings/python/RichRecTemplateRings/Configuration.py
 
In order to activate this package from Brunel
  In a file similar to ../Rec/Brunel/options/DC06-Files.py 
  use the database tags for 2009 data and use a data file similar to that indicated below.

 from Configurables import Brunel, LHCbApp, RichRecSysConf

LHCbApp().DDDBtag   = "head-20090108"
LHCbApp().CondDBtag = "sim-20090108"
from Configurables import DDDBConf
 The following line only if you have a private version of the geometry db.
 DDDBConf().DbRoot = "$SQLDDDBROOT/db/myDDDB/lhcb.xml"
 If you are using a digi file as an input, then some like the following would be useful

    datasetName = "Gauss_bbar_200ev_run1k"
    InpdatasetName ="Gauss_bbar_200ev_run1k"
    EventSelector().Input = ["DATAFILE='PFN:castor:/castor/cern.ch/user/s/seaso/BooleOutput/bevents/2009Geometry/" + 
     InpdatasetName + ".digi' TYP='POOL_ROOTTREE' OPT='READ'"]
    Brunel().DatasetName = datasetName

  If you want to run for example, 20 events

    nEvts = 20       
   ApplicationMgr( EvtMax     = nEvts )

  Now to activate the RichRecTemplateRings package

  RichRecSysConf().TracklessRingAlgs = ["Template"]
  from Configurables import RichTemplateRingFinderConf
 
# begin options local to RichRecTemplateRings package

  Now for the various configurations inside the !RichRecTemplateRings package 
    The options listed below are for specialed application of this package. One does nto need them for
  a default running of the package. 
 
  To activate the built in monitoring histograms
  RichTemplateRingFinderConf().MakeNtuple = False
  To activate finding the mass of the charged particle going through the radiator.
  from Configurables import Rich__Rec__TemplateRings__RichRingMasterAlg
   master = Rich__Rec__TemplateRings__RichRingMasterAlg() 
  master.RichTemplateRingMassFinderActivate=True

 To activatethe option to run upto a maximum number of
 track segments or all track segments.  Also to activate which all radiators to activate. By default the aerogel is switched off.
 By setting the MinRadiator=0, the aerogel is activated.
 (More info in the RichRingMasterAlg algorithm)
from Configurables import Rich__Rec__TemplateRings__RichRingReconParam
Rich__Rec__TemplateRings__RichRingReconParam().ActivateMaxNumberOfTrackSegmentsInRadiator=False
Rich__Rec__TemplateRings__RichRingReconParam().MinRadiator=0
#end Flags local to Template Rings

#begin options to activate comparing with MC Truth
RichRecQCConf().TemplateRingMonitoring =True
 This requires that in the python options of the RichRecQC package the following must be present.
   in ../Rich/RichRecQC/python/RichRecQC/Configuration.py
   from Configurables import RichTemplateRingFinderConf
    ## Steering options
    __slots__ = {
        "Context": "Offline"  # The context within which to run
       ,"DataType" : "2008" # Data type, can be ['DC06','2008']
       ,"RawMonitoring"             : True
       ,"PidMonitoring"             : True
       ,"PixelMonitoring"           : True
       ,"TrackMonitoring"           : True
       ,"PhotonMonitoring"          : True
       ,"TracklessRingMonitoring"   : False
       ,"TemplateRingMonitoring"    : False    (This is to be set True for activating the comparsion with MC Truth) 
       ,"AlignmentMonitoring"       : False
       ,"HPDIFBMonitoring"          : False
       ,"PidMomentumRanges": [ [2,100], [2,10], [10,70], [70,100] ]
       ,"PidTrackTypes":  [ ["All"], ["Forward","Match"] ] 
       .........................
          # Template rings
        if self.getProp("TemplateRingMonitoring") :
           from Configurables import RichTemplateRingFinderConf
           RichTemplateRingFinderConf().ringsTMoni(self.newSeq(sequence,"RichTracklessRingsMoni"))
           
#end options to activate comparing with MC Truth
In order to store the ntuples and histograms created in this package one may activate the following 
 from Configurables import NTupleSvc
NTupleSvc().Output = ["RICHTUPLE1 DATAFILE='myDirectory/Ntuple_RichRecTemplateRings.root' TYP='ROOT' OPT='NEW'"]

-- SajanEaso - 13 Aug 2009
Edit | Attach | Watch | Print version | History: r2 < r1 | Backlinks | Raw View | WYSIWYG | More topic actions
Topic revision: r2 - 2009-08-13 - SajanEaso
 
    • Cern Search Icon Cern Search
    • TWiki Search Icon TWiki Search
    • Google Search Icon Google Search

    LHCb All webs login

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