StevenGoldfarbWireSag < Main

Main Web>TWikiUsers>StevenGoldfarb>StevenGoldfarbNotes>StevenGoldfarbWireSag (2010-06-10, StevenGoldfarb) (raw view)
<table width="99%">
  <tr>
    <td valign="bottom" width="60">
      <img src="https://twiki.cern.ch/twiki/pub/Atlas/WebPreferences/atlaslogo.png" alt="atlaslogo"/>
    </td>
    <td valign="bottom">
---+!! Steve's Notes
---+!! _Wire Sag Validation (Release 15.6.9.13 - SLC5)_
    </td>
  </tr>
</table>

---
[[StevenGoldfarbAthena][Athena]] | [[StevenGoldfarbVP1][VP1]] | [[StevenGoldfarbDQ2][DQ2]] | [[StevenGoldfarbWireSag][Wire Sag]] | [[StevenGoldfarbCluster][UM Cluster]]
---

%TOC%

---++ Introduction

This is an informal summary of notes on installing and running the MDT wire sag validation software. The validation environment was set up by Niels van Eldik and the Rt correction function was developed by Dan Levin.  I am only trying to put things in place.  Please send comments and corrections to [[mailto:Steven.Goldfarb@cern.ch][Steven.Goldfarb@cern.ch]].

---++ Set Up the CMT Environment

Here we define and export variables for CMT.  The correct version of CMT to use for any particular release is found on the [[http://cern.ch/atlas-computing/projects/releases/status/][Offline Release Status Page]].  Here, we use [[http://cern.ch/atlas-computing/projects/releases/status/15.6.9.13/AtlasProduction][Release 15.6.9.13]].  The corresponding version of CMT is *v1r21*.

   1. Login to the appropriate machine _(the "-Y" option passes the necessary X-Window information)_
   <verbatim>
CERN> ssh -Y myusername@lx64slc5.cern.ch *OR* ssh -Y myusername@atint01.cern.ch (can also use atint02 or atint03)
UM> ssh -Y myusername@umt3int01.physics.lsa.umich.edu (can also use umt3int02 or umt3int03)
   </verbatim>
   1. Setup the cmt environment corresponding to the release _(do whenever there is a new version of CMT)_
   <verbatim>
CERN> source /afs/cern.ch/sw/contrib/CMT/v1r21/mgr/setup.sh
UMT3> source /atlas/data08/OSG/APP/atlas_app/atlas_rel/15.6.9/CMT/v1r21/mgr/setup.sh
   </verbatim>

---++ Set Up the Athena Environment

This sets up the Athena working environment.  More details at Atlas.WorkBookSetAccount .

   1. Create a directory structure  _(do this once for the release)_
   <verbatim>
mkdir -p ~/TestArea/15.6.9.13/cmthome
cd ~/TestArea/15.6.9.13/cmthome
   </verbatim>
   1. At *CERN (on lx64slc5 or atint)*, create a file called =requirements= with the following content:
   <verbatim>
#---------------------------------------------------------------------
set CMTSITE STANDALONE
set SITEROOT /afs/cern.ch/atlas/software/releases/15.6.9
macro ATLAS_TEST_AREA ${HOME}/TestArea
macro ATLAS_DIST_AREA ${SITEROOT}
macro SITE_PROJECT_AREA ${SITEROOT}
macro EXTERNAL_PROJECT_AREA ${SITEROOT}
apply_tag projectArea
apply_tag opt
apply_tag setup
apply_tag simpleTest
apply_tag AtlasProduction
apply_tag 15.6.9.13
use AtlasLogin AtlasLogin-* $(ATLAS_DIST_AREA)
set CMTCONFIG i686-slc5-gcc43-opt
#---------------------------------------------------------------------
   </verbatim>
   1. On *UMT3*, create a file called =requirements= with the following content:
   <verbatim>
#---------------------------------------------------------------------
set CMTSITE STANDALONE
set SITEROOT /atlas/data08/OSG/APP/atlas_app/atlas_rel/15.6.9
macro ATLAS_TEST_AREA ${HOME}/TestArea
macro ATLAS_DIST_AREA ${SITEROOT}
macro SITE_PROJECT_AREA ${SITEROOT}
macro EXTERNAL_PROJECT_AREA ${SITEROOT}
apply_tag projectArea
apply_tag opt
apply_tag setup
apply_tag simpleTest
apply_tag AtlasProduction
apply_tag 15.6.9.13
use AtlasLogin AtlasLogin-* $(ATLAS_DIST_AREA)
set CMTCONFIG i686-slc5-gcc43-opt
#---------------------------------------------------------------------   </verbatim>
   1. Configure cmt (creates the scripts for setting up the environment):
   <verbatim>
cmt config
   </verbatim>
   1. Setup the Athena environmental variables: _(every time you login)_
   <verbatim>
source ~/TestArea/15.6.9.13/cmthome/setup.sh
   </verbatim>
   1. Check the CMT variables:
   <verbatim>
echo $CMTCONFIG
CERN> i686-slc5-gcc43-opt
UMT3> i686-slc5-gcc43-opt
echo $CMTPATH
CERN> ~/TestArea/15.6.9.13:/afs/cern.ch/atlas/software/releases/15.6.9/AtlasProduction/15.6.9.13:(etc.)
UMT3> ~/TestArea/15.6.9.13:/atlas/data08/OSG/APP/atlas_app/atlas_rel/15.6.9/AtlasTier0/15.6.9.13:(etc)
   </verbatim>

---++ Check Out the !MuonRecExample Package

We will run standalone muon reconstruction from within this package

   1. If not done already, set up the environmental variables.
   <verbatim>
cd ~/TestArea/15.6.9.13
source cmthome/setup.sh
   </verbatim>
   1. Find out which version of !UserAnalysis is used (in this case !MuonRecExample-01-10-08):
   <verbatim>
cmt show versions MuonSpectrometer/MuonReconstruction/MuonRecExample
   </verbatim>
   1. On UMT3, get kerberos permission (using CERN afs password):
   <verbatim>
UMT3> /usr/kerberos/bin/kinit myusername@CERN.CH
   </verbatim>
   1. Checkout, configure and build the package
   <verbatim>
cmt co -r MuonRecExample-01-10-08 MuonSpectrometer/MuonReconstruction/MuonRecExample
cd MuonSpectrometer/MuonReconstruction/MuonRecExample/cmt
source setup.sh
gmake
   </verbatim>

---++ Test the !MuonRecExample Package

See more detailed instructions at Atlas.WorkBookReconstruction.

*CURRENTLY FAILS DUE TO LACK OF INPUT FILE*

   1. Set up Athena and !UserAnalysis environments: _(once per login)_
   <verbatim>
cd ~/TestArea/15.6.9.13
source cmthome/setup.sh
cd MuonSpectrometer/MuonReconstruction/MuonRecExample/run
source ../cmt/setup.sh
   </verbatim>
   1. See what options you can run with:
   <verbatim>
athena -h
   </verbatim>
   1. Try running !MuonRec_myTopOptions (as is):
   <verbatim>
get_files -jo MuonRecExample/MuonRec_myTopOptions.py
athena MuonRecExample/MuonRec_myTopOptions.py > myTopOutput.txt &
tail -f myTopOutput.txt
   </verbatim>
   1. Browse the n-tuple to see if something was produced:
   <verbatim>
root CalibrationNtuple.root 
   </verbatim>

---++ Check out the Validation Package for testing PRD

   1. If not done already, set up the environmental variables:
   <verbatim>
cd ~/TestArea/15.6.9.13
source cmthome/setup.sh
   </verbatim>
   1. Check out the appropriate package (MuonPRDTest-00-00-15 has Steve's fixes)
   <verbatim>
cmt co -r MuonPRDTest-00-00-15 MuonSpectrometer/MuonValidation/MuonPRDTest
   </verbatim>
   1. Set up the environment
   <verbatim>
cd MuonSpectrometer/MuonValidation/MuonPRDTest/cmt
source setup.sh
gmake
   </verbatim>

---++ Check out the TrackDistortedSurfaces package (geometrical corrections)

   1. If not done already, set up the environmental variables:
   <verbatim>
cd ~/TestArea/15.6.9.13
source cmthome/setup.sh
   </verbatim>
   1. Check out the package (tag TrkDistortedSurfaces-00-01-02 has Steve's fixes)
   <verbatim>
cmt co -r TrkDistortedSurfaces-00-01-02 Tracking/TrkDetDescr/TrkDistortedSurfaces
   </verbatim>
   1. Build it
   <verbatim>
cd Tracking/TrkDetDescr/TrkDistortedSurfaces/cmt
source setup.sh
gmake
   </verbatim>

---++ Check out the MuonCalibData package (Rt correction function)

   1. If not done already, set up the environmental variables:
   <verbatim>
cd ~/TestArea/15.6.9.13
source cmthome/setup.sh
   </verbatim>
   1. Check out the package (tag MdtCalibData-01-03-25 has Steve's implementation)
   <verbatim>
cmt co -r MdtCalibData-01-03-25 MuonSpectrometer/MuonCalib/MdtCalib/MdtCalibData
   </verbatim>
   1. Build it
   <verbatim>
cd MuonSpectrometer/MuonCalib/MdtCalib/MdtCalibData/cmt
source setup.sh
gmake
   </verbatim>

---++ Check out the MuonCalibSvc package (call to Rt correction functions)

   1. If not done already, set up the environmental variables:
   <verbatim>
cd ~/TestArea/15.6.9.13
source cmthome/setup.sh
   </verbatim>
   1. Check out the package (tag MdtCalibSvc-02-02-12 has Steve's implementation)
   <verbatim>
cmt co -r MdtCalibSvc-02-02-12 MuonSpectrometer/MuonCalib/MdtCalib/MdtCalibSvc
   </verbatim>
   1. Build it
   <verbatim>
cd MuonSpectrometer/MuonCalib/MdtCalib/MdtCalibData/cmt
source setup.sh
gmake
   </verbatim>

---++ Wire Sag Data Samples Generated by Rob Harrington

Using release 15.2.0, Rob Harrington generated data samples of non-interacting muons with and without a magnetic field in the barrel and endcap, sectors 1 and 5.  The goal is to digitize these data samples with and without wire-sag to validate the effects.  Samples in sector 1 should have no sag.  Samples in sector 5 should have maximum sag.

   1. The data samples are replicated in the following directories:
   <verbatim>
afs: /afs/cern.ch/user/m/muonprod/roberth/CSC.singlepart_mu100_15.2.0_ATLAS-GEO-07-00-00
castor: /castor/cern.ch/user/g/goldfarb/CSC.singlepart_mu100_15.2.0_ATLAS-GEO-07-00-00
atint0X: /net/s3_datac/goldfarb/CSC.singlepart_mu100_15.2.0_ATLAS-GEO-07-00-00
   </verbatim>
   1. Here is a listing of directory contents:
   <verbatim>
> ls -l /afs/cern.ch/user/m/muonprod/roberth/CSC.singlepart_mu100_15.2.0_ATLAS-GEO-07-00-00
-rw-r--r--  1 23674 zp 4301 Jun 25 21:59 batchscripts.tar.gz
drwxr-xr-x  2 23674 zp 6144 Jun 25 22:47 CSC.008001.singlepart_mu100_sector1_barrel
drwxr-xr-x  2 23674 zp 6144 Jun 25 21:59 CSC.008002.singlepart_mu100_sector5_barrel
drwxr-xr-x  2 23674 zp 6144 Jun 25 21:59 CSC.008003.singlepart_mu100_sector1_endcap
drwxr-xr-x  2 23674 zp 6144 Jun 25 21:59 CSC.008004.singlepart_mu100_sector5_endcap
-rwxr-xr-x  1 23674 zp 1063 Jun 25 21:59 CSC.008233.singlepart_mu20_barrelEC.py
-rwxr-xr-x  1 23674 zp 1060 Jun 25 21:59 CSC.singlepart_mu100_sector1_barrel.py
-rwxr-xr-x  1 23674 zp 1060 Jun 25 21:59 CSC.singlepart_mu100_sector1_endcap.py
-rwxr-xr-x  1 23674 zp 1060 Jun 25 21:59 CSC.singlepart_mu100_sector5_barrel.py
-rwxr-xr-x  1 23674 zp 1060 Jun 25 21:59 CSC.singlepart_mu100_sector5_endcap.py
-rwxr-xr-x  1 23674 zp  133 Jun 25 21:59 IDCaloOffDigitConfig.py
-rw-r--r--  1 23674 zp 1149 Jun 25 22:01 README.txt
-rwxr-xr-x  1 23674 zp 1981 Jun 25 21:59 submit_mu100.sh
   </verbatim>
   1. Check the README.txt file to understand the contents of the files.  The geometry tags are
   <verbatim>
field: ATLAS-GEO-07-00-00
no field: ATLAS-GEONF-07-00-00
   </verbatim>

---++ Modify !jobOptions to Run !PRDTest Algorithm on One of the Data Sets

This algorithm creates histograms comparing !SimHits with ROT. It plots the residuals and other relevant values.

   1. Start with the example:
<verbatim>
   cp !MuonRec_myTopOptions.py !PRDTestField.py
</verbatim>
   1. Add the following lines to *!PRDTestField.py* to run on one of the data samples (endcap sector 1, with B field):
   <verbatim>
from MuonRecExample.DataFiles import DataInfo,DataFilesInfoMap
DataFilesInfoMap['wField'] = DataInfo(["/afs/cern.ch/user/m/muonprod/roberth/CSC.singlepart_mu100_15.2.0_ATLAS-GEO-07-00-00/CSC.008003.singlepart_mu100_sector1_endcap/CSC.008003.singlepart_mu100_sector1_endcap_magfield_[0-2].HITS.pool.root"],
                                    'HITS',
                                    'ATLAS-GEO-07-00-00' )
muonRecFlags.InputData = 'wField'
   </verbatim>
   1. Add the following lines to specify the number of events:
   <verbatim>
if not 'EvtMax' in dir():
    EvtMax = 250
   </verbatim>
   1. Add the following flags:
   <verbatim>
doMuonRecStatistics = True
doAuditors          = True
recFlags.doCBNT     = True
doTrkNtuple         = False
readMuonDigit       = False
   </verbatim>
   1. After =include("MuonRecExample/MuonRec_topOptions.py")= add the following flags to include the PRD testing algorithm:
   <verbatim>
    include("MuonPRDTest/MDTPRDTest.py")
    if muonRecFlags.InputData() == 'wField':
        topSequence.MDTPRDValAlg.Propagator = ToolSvc.MuonPropagator
   </verbatim>

---++ Example !jobOption Files

!jobOption files have been created to run over events in the data sets described above using the following variations: sector 1/sector 5, barrel/endcap, field/no field, geometric sag in digitization, rt sag effects in digitization, geometric sag corrections in reconstruction, rt sag corrections in calibration.  Here are some examples:

   * [[%ATTACHURL%/sector1_endcap_nofield_nogeosag_nortsag_nogeocor_nortcor_jobOptions.py.txt][sector1_endcap_nofield_nogeosag_nortsag_nogeocor_nortcor_jobOptions.py.txt]]
   * [[%ATTACHURL%/sector5_endcap_nofield_nogeosag_nortsag_nogeocor_nortcor_jobOptions.py.txt][sector5_endcap_nofield_nogeosag_nortsag_nogeocor_nortcor_jobOptions.py.txt]]
   * [[%ATTACHURL%/sector5_endcap_nofield_geosag_nortsag_nogeocor_nortcor_jobOptions.py.txt][sector5_endcap_nofield_geosag_nortsag_nogeocor_nortcor_jobOptions.py.txt]]
   * [[%ATTACHURL%/sector5_endcap_nofield_geosag_nortsag_geocor_nortcor_jobOptions.py.txt][sector5_endcap_nofield_geosag_nortsag_geocor_nortcor_jobOptions.py.txt]]
   * [[%ATTACHURL%/sector5_endcap_magfield_nogeosag_nortsag_nogeocor_nortcor_jobOptions.py.txt][sector5_endcap_magfield_nogeosag_nortsag_nogeocor_nortcor_jobOptions.py.txt]]

---++ Run the Test

The following output files will be created:
   * *config.txt* - Configuration used in job
   * *PRDTestFieldOutput.txt* Output from running the job
   * *MDTPRDTest.root* root file containing n-tuples, histograms

   1. Run the Test :
   <verbatim>
athena PRDTestField.py >! PRDTestFieldOutput.txt &
tail -f PRDTestFieldOutput.txt
   </verbatim>
   1. Copy the macro file from MuonValidation
   <verbatim>
cp ~/TestArea/15.6.3.6/MuonSpectrometer/MuonValidation/MuonPRDTest/macros/MDTPRDVal.C .
   </verbatim>
   1. Be sure to modify the macros to handle the correct file names, etc.
   1. Open the root file:
   <verbatim>
root MDTPRDTest.root
   </verbatim>
   1. Load the macros
   <verbatim>
.L MDTPRDVal.C
   </verbatim>
   1. Plot the residuals
   <verbatim>
draw()
   </verbatim>

---++ Example Plots

The attached talk contains plots generated with the above !jobOptions.  There are still many technical problems to be resolved (track propagation, applying geometrical effects, etc.)

   * Analysis (28/07/2009): [ [[%ATTACHURL%/WireSag-20090728.pdf][pdf]] ] [ [[%ATTACHURL%/WireSag-20090728.ppt][ppt]] ]
   * Analysis (14/09/2009): [ [[%ATTACHURL%/WireSag-20090914.pdf][pdf]] ] [ [[%ATTACHURL%/WireSag-20090914.ppt][ppt]] ]

---++ Interpreting the Plots

| *Variable Name*       | *Variable Meaning*  |
| MDT_SimRadius      | Local R of pure simulated drift circle  |
| MDT_LocalZ            | Local Z of pure simulated drift circle  |
| MDT_SdoRadius      | Local R of simulated drift circle with RT corrections  |
| MDT_ExNomWireR  | Local R from wire to extrapolated track  |
| MDT_ExNomWireZ   | Local Z intercept from wire to extrapolated track  |
| MDT_ExSagWireR   | Local R from corrected wire to extrapolated track  |
| MDT_ExSagWireZ   | Local Z intercept from corrected wire to extrapolated track  |
| MDT_RotRadius      | Local R of reconstructed drift circle  |
| MDT_ExSagRotR    | Local R from corrected ROT to extrapolated track  |

---++ References

   * Atlas.WorkBook - ATLAS Workbook for Computing
   * Atlas.WorkBookReconstruction - ATLAS Workbook Reconstruction Recipe
   * [[https://hep.pa.msu.edu/twiki/bin/view/AGLT2/MuonReco][AGLT2.MuonReco]] - Muon Cosmic Reconstruction Tutorial (Dec 2008)
   * Atlas.RecExCommonFlags - Flags to be set for Reconstruction
   * [[http://cern.ch/atlas-computing/links/nightlyDevDirectory/AtlasOffline/latest_doxygen/InstallArea/doc/MuonRecExample/html/MuonDataRec__myTopOptions_8py-source.html][MuonDataRec_myTopOptions.py]] - Latest Top !jobOptions for Muon Reconstruction
   * Atlas.AtlasGeomDBTags - Geometry Tags Corresponding to Releases

---
*Major updates*:%BR%
Main.StevenGoldfarb - 23 Apr 2009 (creation)
-- Main.StevenGoldfarb - 28 Apr 2009 (added WireSag instructions)
-- Main.StevenGoldfarb - 01 Jul 2009 (modified for 15.2.0)
-- Main.StevenGoldfarb - 13 Jul 2009 (modified for 15.5.2.2)
-- Main.StevenGoldfarb - 28 Jul 2009 (linked some example job, plots)
-- Main.StevenGoldfarb - 2009-09-14 (added plot interpretation, analysis update)
-- Main.StevenGoldfarb - 23-Oct-2009 (modified for 15.3.1.7)
-- Main.StevenGoldfarb - 04-Nov-2009 (modified for 15.5.1.8)
-- Main.StevenGoldfarb - 04-Nov-2009 (reverted to 15.5.2.2)
-- Main.StevenGoldfarb - 07-Nov-2009 (modified for 15.5.2.2)
-- Main.StevenGoldfarb - 12-Feb-2010 (modified for 15.6.3.6 - SLC5 - and MdtCalibSvc changes)
-- Main.StevenGoldfarb - 01-Mar-2010 (updated some tags)
-- Main.StevenGoldfarb - 10-Jun-2010 (modified for 15.6.9.13)
Attachments
Topic attachments
I	Attachment	History	Action	Size	Date	Who	Comment
pdf	WireSag-20090728.pdf	r2 r1	manage	1332.2 K	2009-07-29 - 11:14	StevenGoldfarb	First Plots from Analysis
ppt	WireSag-20090728.ppt	r2 r1	manage	890.0 K	2009-07-29 - 11:15	StevenGoldfarb	First Plots from Analysis
pdf	WireSag-20090914.pdf	r1	manage	1639.2 K	2009-09-14 - 17:41	StevenGoldfarb	Wire Sag Presentation
ppt	WireSag-20090914.ppt	r1	manage	1262.0 K	2009-09-14 - 17:40	StevenGoldfarb	Wire Sag Presentation
txt	sector1_endcap_nofield_nogeosag_nortsag_nogeocor_nortcor_jobOptions.py.txt	r1	manage	2.3 K	2009-07-28 - 17:47	StevenGoldfarb
txt	sector5_endcap_magfield_nogeosag_nortsag_nogeocor_nortcor_jobOptions.py.txt	r1	manage	2.4 K	2009-07-28 - 17:49	StevenGoldfarb
txt	sector5_endcap_nofield_geosag_nortsag_geocor_nortcor_jobOptions.py.txt	r1	manage	2.3 K	2009-07-28 - 17:49	StevenGoldfarb
txt	sector5_endcap_nofield_geosag_nortsag_nogeocor_nortcor_jobOptions.py.txt	r1	manage	2.3 K	2009-07-28 - 17:48	StevenGoldfarb
txt	sector5_endcap_nofield_nogeosag_nortsag_nogeocor_nortcor_jobOptions.py.txt	r1	manage	2.3 K	2009-07-28 - 17:47	StevenGoldfarb
Topic revision: r30 - 2010-06-10 - StevenGoldfarb
Main
Webs
Public webs
Welcome Guest
- Cern Search
- TWiki Search
- Google Search
Main All webs
Copyright &© 2008-2024 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