Difference: SWGuideTkAlignmentValidationEnergyOverMomentum (1 vs. 12)

Revision 122013-05-10 - PalHidas

Line: 1 to 1
Changed:
<
<
META TOPICPARENT name="SWGuideTrackAlignment"
>
>
META TOPICPARENT name="TkAlignmentValidationTools"
 

Energy Over Momentum (E/p) Alignment Validation Tool

Complete: 2

Revision 112012-03-30 - ChristopheGoetzmann

Line: 1 to 1
 
META TOPICPARENT name="SWGuideTrackAlignment"

Energy Over Momentum (E/p) Alignment Validation Tool

Line: 144 to 144
 
Changed:
<
<
* variable (string): "eta" for twist validation, "phi", "phi1", "phi2", and "phi3" for validation of phi-modes for all tracks, eta<-0.9, |eta|<0.9, eta>0.9, respectively.
>
>
* variable (string): "eta" for twist validation, "phi", "phi1", "phi2", and "phi3" for validation of phi-modes for all tracks, eta<-0.9, |eta|<0.9, eta>0.9, respectively.
  * path (string): path to the produced tree files int the form: "/myRep/rep/.../"

Revision 102012-03-01 - ChristopheGoetzmann

Line: 1 to 1
 
META TOPICPARENT name="SWGuideTrackAlignment"

Energy Over Momentum (E/p) Alignment Validation Tool

Line: 184 to 184
 
  • 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.
Added:
>
>
  • Lower threshold of 50 GeV on energy deposit in 3x3 Hadronic Calorimeter tower cluster.
 

Outputs

Revision 92012-03-01 - ChristopheGoetzmann

Line: 1 to 1
 
META TOPICPARENT name="SWGuideTrackAlignment"

Energy Over Momentum (E/p) Alignment Validation Tool

Line: 187 to 187
 

Outputs

Changed:
<
<
Results obtained with global tag "GR_R_44_V13", and tracker geometry "TrackerAlignment_GR10_v4_offline"
>
>
Results obtained with global tag "GR_R_44_V13", and tracker geometry "TrackerAlignment_GR10_v4_offline", with dataset: /Commissioning/Run2011A-HcalCalIsoTrk-08Nov2011-v1/ALCARECO.
 
  • Sagitta validation:
    sagitta_validation_all.gif

Revision 82012-02-29 - ChristopheGoetzmann

Line: 1 to 1
 
META TOPICPARENT name="SWGuideTrackAlignment"

Energy Over Momentum (E/p) Alignment Validation Tool

Line: 190 to 190
 Results obtained with global tag "GR_R_44_V13", and tracker geometry "TrackerAlignment_GR10_v4_offline"

  • Sagitta validation:
Changed:
<
<
sagitta_validation_all.gif
>
>
sagitta_validation_all.gif
 
  • Sagitta validation for low eta:
Changed:
<
<
sagitta_validation_lowEta.gif
>
>
sagitta_validation_lowEta.gif
 
  • Sagitta validation for central eta:
Changed:
<
<
sagitta_validation_centralEta.gif
>
>
sagitta_validation_centralEta.gif
 
  • Sagitta validation for high eta:
Changed:
<
<
sagitta_validation_highEta.gif
>
>
sagitta_validation_highEta.gif
 
  • Twist validation:
Changed:
<
<
twist_validation.gif
>
>
twist_validation.gif
 

-- HolgerEnderle - 14-Nov-2011

Revision 72012-02-29 - ChristopheGoetzmann

Line: 1 to 1
 
META TOPICPARENT name="SWGuideTrackAlignment"

Energy Over Momentum (E/p) Alignment Validation Tool

Line: 187 to 187
 

Outputs

Added:
>
>
Results obtained with global tag "GR_R_44_V13", and tracker geometry "TrackerAlignment_GR10_v4_offline"
 
  • Sagitta validation:
    sagitta_validation_all.gif
Line: 204 to 206
 

-- HolgerEnderle - 14-Nov-2011

Changed:
<
<
-- ChristopheGoetzmann - 23-Feb-2012
>
>
-- ChristopheGoetzmann - 29-Feb-2012
 
META FILEATTACHMENT attachment="sagitta_validation_all.gif" attr="" comment="" date="1330530304" name="sagitta_validation_all.gif" path="sagitta_validation_all.gif" size="18166" stream="sagitta_validation_all.gif" tmpFilename="/usr/tmp/CGItemp52980" user="cgoetzma" version="1"
META FILEATTACHMENT attachment="sagitta_validation_centralEta.gif" attr="" comment="" date="1330530317" name="sagitta_validation_centralEta.gif" path="sagitta_validation_centralEta.gif" size="19214" stream="sagitta_validation_centralEta.gif" tmpFilename="/usr/tmp/CGItemp52946" user="cgoetzma" version="1"

Revision 62012-02-29 - ChristopheGoetzmann

Line: 1 to 1
 
META TOPICPARENT name="SWGuideTrackAlignment"

Energy Over Momentum (E/p) Alignment Validation Tool

Line: 136 to 136
 

Running the Validation Macro

Added:
>
>
4 variables are mandatory:


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

* variable (string): "eta" for twist validation, "phi", "phi1", "phi2", and "phi3" for validation of phi-modes for all tracks, eta<-0.9, |eta|<0.9, eta>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:

 
Changed:
<
<
root -l -b -q 'momentumBiasValidation.C+(variable,path,files+legend label)'
>
>
root -l -b -q 'momentumBiasValidation.C+(variable, path, files+legend label, outputType, verbose, givenMin, givenMax)'
 
Changed:
<
<
* variable: "eta" for twist validation, "phi", "phi1", "phi2", and "phi3" for validation of phi-modes for all tracks, eta<-0.9, |eta|<0.9, eta>0.9, respectively.
>
>
* verbose (bool): set true for verbose mode (default: false)
 
Changed:
<
<
* path: path to the produced tree files
>
>
* 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\\..."
Line: 166 to 185
 
  • 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.
Added:
>
>

Outputs

  • Sagitta validation:
    sagitta_validation_all.gif

  • Sagitta validation for low eta:
    sagitta_validation_lowEta.gif

  • Sagitta validation for central eta:
    sagitta_validation_centralEta.gif

  • Sagitta validation for high eta:
    sagitta_validation_highEta.gif

  • Twist validation:
    twist_validation.gif
 -- HolgerEnderle - 14-Nov-2011 -- ChristopheGoetzmann - 23-Feb-2012
Added:
>
>
META FILEATTACHMENT attachment="sagitta_validation_all.gif" attr="" comment="" date="1330530304" name="sagitta_validation_all.gif" path="sagitta_validation_all.gif" size="18166" stream="sagitta_validation_all.gif" tmpFilename="/usr/tmp/CGItemp52980" user="cgoetzma" version="1"
META FILEATTACHMENT attachment="sagitta_validation_centralEta.gif" attr="" comment="" date="1330530317" name="sagitta_validation_centralEta.gif" path="sagitta_validation_centralEta.gif" size="19214" stream="sagitta_validation_centralEta.gif" tmpFilename="/usr/tmp/CGItemp52946" user="cgoetzma" version="1"
META FILEATTACHMENT attachment="sagitta_validation_highEta.gif" attr="" comment="" date="1330530328" name="sagitta_validation_highEta.gif" path="sagitta_validation_highEta.gif" size="19218" stream="sagitta_validation_highEta.gif" tmpFilename="/usr/tmp/CGItemp53034" user="cgoetzma" version="1"
META FILEATTACHMENT attachment="sagitta_validation_lowEta.gif" attr="" comment="" date="1330530338" name="sagitta_validation_lowEta.gif" path="sagitta_validation_lowEta.gif" size="19262" stream="sagitta_validation_lowEta.gif" tmpFilename="/usr/tmp/CGItemp53028" user="cgoetzma" version="1"
META FILEATTACHMENT attachment="twist_validation.gif" attr="" comment="" date="1330530355" name="twist_validation.gif" path="twist_validation.gif" size="20422" stream="twist_validation.gif" tmpFilename="/usr/tmp/CGItemp53045" user="cgoetzma" version="1"

Revision 52012-02-29 - EricConte

Line: 1 to 1
 
META TOPICPARENT name="SWGuideTrackAlignment"

Energy Over Momentum (E/p) Alignment Validation Tool

Line: 15 to 15
 

Checking out the Code

Added:
>
>

Stable release

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

* interface/EopVariables.h (co)

Line: 37 to 39
  This version is tested in CMSSW_4_2_X and CMSSW_4_4_X.
Added:
>
>

release in development

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

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

 

Configuring the Tree Writer

1. define the input files:

Revision 42012-02-23 - ChristopheGoetzmann

Line: 1 to 1
 
META TOPICPARENT name="SWGuideTrackAlignment"

Energy Over Momentum (E/p) Alignment Validation Tool

Line: 140 to 140
 
  • 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:

Changed:
<
<
  • Trigger selection: high pT QVD events (as they frequently contain isolated tracks). Lower threshold of 38 GeV on track momentum magnitude.
>
>
  • 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.

Revision 32012-02-23 - ChristopheGoetzmann

Line: 1 to 1
 
META TOPICPARENT name="SWGuideTrackAlignment"

Energy Over Momentum (E/p) Alignment Validation Tool

Line: 130 to 130
  * 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\\..."
Added:
>
>

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 QVD 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.
 -- HolgerEnderle - 14-Nov-2011
Deleted:
<
<
-- ChristopheGoetzmann - 01-Feb-2012
 \ No newline at end of file
Added:
>
>
-- ChristopheGoetzmann - 23-Feb-2012

Revision 22012-02-01 - ChristopheGoetzmann

Line: 1 to 1
 
META TOPICPARENT name="SWGuideTrackAlignment"

Energy Over Momentum (E/p) Alignment Validation Tool

Line: 91 to 91
 
Added:
>
>

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

Line: 106 to 131
  * 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\\..."

-- HolgerEnderle - 14-Nov-2011

Added:
>
>
-- ChristopheGoetzmann - 01-Feb-2012
 \ No newline at end of file

Revision 12011-11-14 - HolgerEnderle

Line: 1 to 1
Added:
>
>
META TOPICPARENT name="SWGuideTrackAlignment"

Energy Over Momentum (E/p) Alignment Validation Tool

Complete: 2

Scope

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

Running

Checking out the Code

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.

Configuring the Tree Writer

1. define the input files:

<!-- SyntaxHighlightingPlugin -->
process.source = cms.Source( "PoolSource",
    fileNames = cms.untracked.vstring(
    'yourInputFiles.root'
    )
)
<!-- end SyntaxHighlightingPlugin -->

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

<!-- SyntaxHighlightingPlugin -->
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")
<!-- end SyntaxHighlightingPlugin -->

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

<!-- SyntaxHighlightingPlugin -->
process.TrackRefitter.src = cms.InputTag( 'yourInputTag')
<!-- end SyntaxHighlightingPlugin -->

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

<!-- SyntaxHighlightingPlugin -->
process.energyOverMomentumTree.src = cms.InputTag('yourInputTag') # 'TkAlIsoProd:' for HcalCalIsoTrk files
process.p = cms.Path(process.offlineBeamSpot*process.energyOverMomentumTree)
<!-- end SyntaxHighlightingPlugin -->

Running the Tree Writer


cmsRun eopTreeWriter_cfg.py

Running the Validation Macro


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

* variable: "eta" for twist validation, "phi", "phi1", "phi2", and "phi3" for validation of phi-modes for all tracks, eta<-0.9, |eta|<0.9, eta>0.9, respectively.

* path: path to the produced tree files

* 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\\..."

-- HolgerEnderle - 14-Nov-2011

 
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