Line Code
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64
import FWCore.ParameterSet.Config as cms


from Alignment.APEEstimation.ApeEstimator_cfi import *
from Alignment.APEEstimation.SectorBuilder_cff import *


ApeEstimator = ApeEstimatorTemplate.clone(
  maxTracksPerEvent = 0,
  #applyTrackCuts = False,
  minGoodHitsPerTrack = 1,
  residualErrorBinning = [0.0005,0.0010,0.0015,0.0020,0.0025,0.0030,0.0035,0.0040,0.0050,0.0070,0.0100], # 5-100um
  #zoomHists = False,
  vErrHists = [1],
  #Sectors = SubdetSectors,
  #Sectors = TIBTOBQuarters,
  #Sectors = TIBTOBQuarters2DSeparation,
  #Sectors = TIBTOBPitchAnd2DSeparation,
  #Sectors = TIBTOBLayerAndOrientationSeparation,
  #Sectors = TIDTECSideAndRingAndOrientationSeparation,
  Sectors = RecentSectors,
  tjTkAssociationMapTag = "TrackRefitterHighPurityForApeEstimator",
)
ApeEstimator.HitSelector.width = [3,3]
ApeEstimator.HitSelector.maxIndex = [1,1]
#ApeEstimator.HitSelector.edgeStrips = [2,800]   # exclude first (and so also last) strip
ApeEstimator.HitSelector.sOverN = [20.,50.]

ApeEstimator.HitSelector.chargePixel = [10000., 2000000.]
ApeEstimator.HitSelector.widthX = [2,1000]
ApeEstimator.HitSelector.widthY = [2,1000]
ApeEstimator.HitSelector.logClusterProbability = [-5.,1.]
ApeEstimator.HitSelector.isOnEdge = [0,0]
ApeEstimator.HitSelector.qBin = [1,3]



# Why is charge and maxCharge double, not int?
#ApeEstimator.HitSelector.maxCharge = [0.,250.]
ApeEstimator.HitSelector.chargeOnEdges = [0.,0.5]


#ApeEstimator.HitSelector.phiSensX = [-1.0472,1.0472]  # [-60,60] degree
#ApeEstimator.HitSelector.phiSensY = [-1.0472,1.0472]  # [-60,60] degree


#ApeEstimator.HitSelector.errXHit = cms.vdouble(0.,0.0060)  # 60um, to exclude very large clusters






ApeAnalyzer = ApeEstimator.clone(
    Sectors = ValidationSectors,
    analyzerMode = True,
    calculateApe = True,
)