(r8) SWGuideTkAlignmentValidationEnergyOverMomentum < CMSPublic

CMSPublic Web>SWGuide>SWGuideCalAli>SWGuideTrackAlignment>SWGuideTkAlignmentValidationEnergyOverMomentum (revision 8) (raw view)~~EditAttachPDF~~

---+!! Energy Over Momentum (E/p) Alignment Validation Tool

%COMPLETE2%

%TOC{title="Contents:"}%


---++ Scope

This tool is developed to determine momentum changing weak modes in the tracker geometry.

---++ Running

---+++ Checking out the Code

---++++ Stable release

At the moment, you have to check out or update following files from the HEAD of Alignment/OfflineValidation

* interface/EopVariables.h (co)

* macros/momentumBiasValidation (co)

* plugins/EopTreeWriter.h (co)

* plugins/EopTreeWriter.cc (co)

* plugins/BuildFile.xml (update)

* python/eopTreeWriter_cfi.py (co)

* src/classes.h (update)

* src/classes_def.xml (update)

* test/eopTreeWriter_cfg.py (co)

This version is tested in CMSSW_4_2_X and CMSSW_4_4_X.

---++++ release in development

The code is under development and the last release is stored in a UserCode folder.

<verbatim>
scramv1 project CMSSW CMSSW_4_4_2
cd CMSSW_4_4_2
cmsenv
cd $CMSSW_BASE/src
cvs co UserCode/IPHC_Eop/install.sh
source UserCode/IPHC_Eop/install.sh
</verbatim>




---+++ Configuring the Tree Writer

1. define the input files:

%SYNTAX{ syntax="python"}%
process.source = cms.Source( "PoolSource",
    fileNames = cms.untracked.vstring(
    'yourInputFiles.root'
    )
)
%ENDSYNTAX%

2. choose the geometry you want to validate optionally with alignment errors and/or surface deformations:

%SYNTAX{ syntax="python"}%
from CondCore.DBCommon.CondDBSetup_cfi import *
process.trackerAlignment = cms.ESSource("PoolDBESSource",CondDBSetup,                                        
                                        connect = cms.string("yourDbFile"),
                                        toGet = cms.VPSet(cms.PSet(record = cms.string("TrackerAlignmentRcd"),
                                                                   tag = cms.string("yourAlignmentTag")
                                                                   )## ,
##                                                             cms.PSet(record = cms.string("TrackerAlignmentErrorRcd"),
##                                                                   tag = cms.string("yourAlignmentErrorTag")
##                                                                   ),
##                                                             cms.PSet(record = cms.string("TrackerSurfaceDeformationRcd"),
##                                                                   tag = cms.string("yourDeformationTag")
##                                                                   )
                                                            )
                                           )
process.prefer( "trackerAlignment")
%ENDSYNTAX%

(3.) if you run on files other than *HcalCalIsoTrk*, adapt the input tags:

%SYNTAX{ syntax="python"}%
process.TrackRefitter.src = cms.InputTag( 'yourInputTag')
%ENDSYNTAX%

(4.) if you want to validate the geometry used for the track reconstruction, switch of the track refitter:

%SYNTAX{ syntax="python"}%
process.energyOverMomentumTree.src = cms.InputTag('yourInputTag') # 'TkAlIsoProd:' for HcalCalIsoTrk files
process.p = cms.Path(process.offlineBeamSpot*process.energyOverMomentumTree)

%ENDSYNTAX%

---+++ Running the Tree Writer

<verbatim>

cmsRun eopTreeWriter_cfg.py

</verbatim>

---+++ Variables stored in Tree

| *Variable* | *comments* |
| track_outerRadius | reco::Track::outerRadius() |
| track_chi2 | reco::TrackBase::chi2() |
| track_normalizedChi2 | reco::TrackBase::normalizedChi2() |
| track_p | |
| track_pt | |
| track_ptError | |
| track_theta | |
| track_eta | |
| track_phi | |
| track_emc1 | energy deposit in Ecal considering 1x1 cell |
| track_emc3 | energy deposit in Ecal considering 3x3 cells |
| track_emc5 | energy deposit in Ecal considering 5x5 cells |
| track_hac1 | energy deposit in Hcal considering 1x1 cell |
| track_hac3 | energy deposit in Hcal considering 3x3 cells |
| track_hac5 | energy deposit in Hcal considering 5x5 cells |
| track_maxPNearby | for the considered track, gives the momentum magnitude of the "higher momentum track" in a cone whose differential radius is defined by the user. The informations used for calculation are the eta and phi parameters of propagated tracks in Ecal. Default value is -10 GeV (i.e. if there is no other tracks in the cone). An upper threshold on this value can be used as an isolation criteria.  |
| distofmax | Gives the radial distance in Ecal between the considered track and the eventual "higher momentum track" contained in the cone defined on the previous line. Default value is 50 cm (i.e. if there is no other tracks in the cone) |
| track_charge | |
| track_nHits | reco::TrackBase::numberOfValidHits() |
|track_nLostHits | reco::TrackBase::numberOfLostHits() |
| track_innerOK | reco::Track::innerOk() |

---+++ Running the Validation Macro

*4 variables are mandatory:*

<verbatim>

root -l -b -q 'momentumBiasValidation.C+(variable, path, files+legend label, outputType)'

</verbatim>

 * variable (string): "eta" for twist validation, "phi", "phi1", "phi2", and *"phi3"* for validation of phi-modes for all tracks, eta&lt;-0.9, |eta|&lt;0.9, eta&gt;0.9, respectively.

 * path (string): path to the produced tree files int the form: "/myRep/rep/.../"

 * files+legend label (string): list of the files with the corresponding legend entry in the following form: "firstFile.root=first legend entry\\secondFile.root=second legend entry\\..."

 * outputType (string): set the type of the produced plots, e.g. "gif", "root", "eps", etc...


*3 more optional variables:*

<verbatim>

root -l -b -q 'momentumBiasValidation.C+(variable, path, files+legend label, outputType, verbose, givenMin, givenMax)'

</verbatim>

 * verbose (bool): set true for verbose mode (default: false)

 * givenMin/givenMax (double): set the lowest and the highest value of the Y-axis of the output plot

 * files+legend label: list of the files with the corresponding legend entry in the following form: "firstFile.root=first legend entry\\secondFile.root=second legend entry\\..."

---+++ Track Selection

There are two levels were selections criteria are applied on tracks: in the tree writer ='eopTreeWriter_cfg.py'= and in the validation macro ='momentumBiasValidation.C'=.

The selection criteria applied for the Tree Writer are:
   * Lower threshold of 9 GeV on track momentum magnitude.
   * Lower threshold of 5 GeV on energy deposit in 3x3 Hadronic Calorimeter tower cluster.
   * For a given track, no track of more than 5 GeV is allowed in a dR of 40cm.

The selection criteria applied for the validation macro are:
   * Trigger selection: high pT QCD events (as they frequently contain isolated tracks). Lower threshold of 38 GeV on track momentum magnitude.
   * Number of valid hits greater or equal to 13.
   * No lost hits.
   * Radius of outermost hit greater than 99cm.
   * Normalized chi square smaller than 5.
   * Upper threshold of 1 GeV on energy deposit in ECAL (in a 9cm radius cylinder), to ensure that particles do not shower until they enter the HCAL.
   * Upper threshold of 8 GeV on energy deposit in ECAL between to concentric cylinders (respective radii: 15 and 35 cm) arround impact point, to ensure isolation against neutral particles.

---+++ Outputs

Results obtained with global tag "GR_R_44_V13", and tracker geometry "TrackerAlignment_GR10_v4_offline"

   * Sagitta validation: <br />
     <img src="%ATTACHURLPATH%/sagitta_validation_all.gif" alt="sagitta_validation_all.gif" width='917' height='618' />

   * Sagitta validation for low eta: <br />
     <img src="%ATTACHURLPATH%/sagitta_validation_lowEta.gif" alt="sagitta_validation_lowEta.gif" width='917' height='618' />

   * Sagitta validation for central eta: <br />
     <img src="%ATTACHURLPATH%/sagitta_validation_centralEta.gif" alt="sagitta_validation_centralEta.gif" width='917' height='618' />

   * Sagitta validation for high eta: <br />
     <img src="%ATTACHURLPATH%/sagitta_validation_highEta.gif" alt="sagitta_validation_highEta.gif" width='917' height='618' />

   * Twist validation: <br />
     <img src="%ATTACHURLPATH%/twist_validation.gif" alt="twist_validation.gif" width='917' height='618' />


-- Main.HolgerEnderle - 14-Nov-2011
-- Main.ChristopheGoetzmann - 29-Feb-2012

Attachments

Topic attachments
I	Attachment	History	Action	Size	Date	Who
gif	sagitta_validation_all.gif	r1	manage	17.7 K	2012-02-29 - 16:45	ChristopheGoetzmann
gif	sagitta_validation_centralEta.gif	r1	manage	18.8 K	2012-02-29 - 16:45	ChristopheGoetzmann
gif	sagitta_validation_highEta.gif	r1	manage	18.8 K	2012-02-29 - 16:45	ChristopheGoetzmann
gif	sagitta_validation_lowEta.gif	r1	manage	18.8 K	2012-02-29 - 16:45	ChristopheGoetzmann
gif	twist_validation.gif	r1	manage	19.9 K	2012-02-29 - 16:45	ChristopheGoetzmann

Topic revision: r8 - 2012-02-29 - ChristopheGoetzmann

CMSPublic

Create a LeftBar

Public webs

- Cern Search
- TWiki Search
- Google Search
CMSPublic All webs

Copyright &© 2008-2021 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