SiStripLorentzAngle/python/MeasureLA_cff.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 byLayer = cms.int32(0)
0004 byModule = cms.int32(1)
0005 byModuleSummary = cms.int32(2)
0006
0007 METHOD_WIDTH =  cms.int32(  1 )
0008 METHOD_PROB1 =  cms.int32(  2 )
0009 METHOD_AVGV2 =  cms.int32(  4 )
0010 METHOD_AVGV3 =  cms.int32(  8 )
0011 METHOD_RMSV2 =  cms.int32( 16 )
0012 METHOD_RMSV3 =  cms.int32( 32 )
0013
0014 def LA_Report(method, granularity, filename) :
0015     return cms.PSet( Method = method,
0016                      Granularity = granularity,
0017                      ReportName = cms.string(filename) )
0018
0019 def LA_Measurement(method, granularity, minimumEntries, maxChi2ndof) :
0020     return cms.PSet( Method = method,
0021                      Granularity = granularity,
0022                      MinEntries = cms.uint32(minimumEntries),
0023                      MaxChi2ndof = cms.double(maxChi2ndof) )
0024
0025 uncalSlopes  = cms.vdouble(1,1,1,1,1,1,1,1,1,1,1,1,1,1)
0026 uncalOffsets = cms.vdouble(0,0,0,0,0,0,0,0,0,0,0,0,0,0)
0027 uncalPulls   = cms.vdouble(1,1,1,1,1,1,1,1,1,1,1,1,1,1)
0028
0029 LorentzAngleCalibrations_Uncalibrated = cms.VPSet(
0030     cms.PSet( Method = METHOD_WIDTH, Slopes = uncalSlopes , Offsets = uncalOffsets, Pulls = uncalPulls),
0031     cms.PSet( Method = METHOD_PROB1, Slopes = uncalSlopes , Offsets = uncalOffsets, Pulls = uncalPulls),
0032     cms.PSet( Method = METHOD_AVGV2, Slopes = uncalSlopes , Offsets = uncalOffsets, Pulls = uncalPulls),
0033     cms.PSet( Method = METHOD_AVGV3, Slopes = uncalSlopes , Offsets = uncalOffsets, Pulls = uncalPulls),
0034     cms.PSet( Method = METHOD_RMSV2, Slopes = uncalSlopes , Offsets = uncalOffsets, Pulls = uncalPulls),
0035     cms.PSet( Method = METHOD_RMSV3, Slopes = uncalSlopes , Offsets = uncalOffsets, Pulls = uncalPulls)
0036     )
0037
0038 LorentzAngleCalibrations_PeakModeBEAM = cms.VPSet(
0039     cms.PSet( Method = METHOD_WIDTH, Slopes = uncalSlopes , Offsets = uncalOffsets, Pulls = uncalPulls),
0040     cms.PSet( Method = METHOD_PROB1, Slopes = uncalSlopes , Offsets = uncalOffsets, Pulls = uncalPulls),
0041     cms.PSet( Method = METHOD_AVGV2, Slopes = uncalSlopes , Offsets = uncalOffsets, Pulls = uncalPulls),
0042     cms.PSet( Method = METHOD_AVGV3, Slopes = uncalSlopes , Offsets = uncalOffsets, Pulls = uncalPulls),
0043     cms.PSet( Method = METHOD_RMSV2, Slopes = uncalSlopes , Offsets = uncalOffsets, Pulls = uncalPulls),
0044     cms.PSet( Method = METHOD_RMSV3, Slopes = uncalSlopes , Offsets = uncalOffsets, Pulls = uncalPulls)
0045     )
0046
0047 LorentzAngleCalibrations_DeconvolutionModeBEAM = cms.VPSet(
0048     cms.PSet( Method = METHOD_WIDTH, Slopes = uncalSlopes , Offsets = uncalOffsets, Pulls = uncalPulls),
0049     cms.PSet( Method = METHOD_PROB1, Slopes = uncalSlopes , Offsets = uncalOffsets, Pulls = uncalPulls),
0050     cms.PSet( Method = METHOD_AVGV2, Slopes = uncalSlopes , Offsets = uncalOffsets, Pulls = uncalPulls),
0051     cms.PSet( Method = METHOD_AVGV3, Slopes = uncalSlopes , Offsets = uncalOffsets, Pulls = uncalPulls),
0052     cms.PSet( Method = METHOD_RMSV2, Slopes = uncalSlopes , Offsets = uncalOffsets, Pulls = uncalPulls),
0053     cms.PSet( Method = METHOD_RMSV3, Slopes = uncalSlopes , Offsets = uncalOffsets, Pulls = uncalPulls)
0054     )
0055
0056 # Calibration found from ensemble of 30 pseudo-datasets at each of 6 points, tanLA = {-0.15,-0.12,-0.09,-0.06,-0.03,0.005}
0057 # Sample was 2/3 of /TkCosmics38T/Summer09-STARTUP31X_V3-v1/GEN-SIM-DIGI-RAW at each tanLA point
0058 # redigitized in peak mode with peak mode noises, CMSSW_3_3_2
0059 LorentzAngleCalibrations_PeakModeCOSMIC = cms.VPSet(
0060     cms.PSet( Method = METHOD_WIDTH,
0061               Slopes = cms.vdouble(0.985407,0.986015,0.982955,0.992045,0.988984,0.987343,0.954721,0.95576,0.954863,0.960566,0.959319,0.956621,0.977128,0.976453),
0062               Offsets = cms.vdouble(0.000363377,0.000347239,-8.73609e-05,0.000633789,0.000785546,0.00063365,0.000723714,0.000720315,0.000644327,0.000814116,0.00164101,0.00123942,0.00180739,0.00171959),
0063               Pulls = cms.vdouble(1.05217,1.06731,0.97208,1.14943,1.07177,0.988835,0.942119,1.025,1.11031,0.930538,1.08191,1.04613,1.03081,1.17641) ),
0064     cms.PSet( Method = METHOD_PROB1,
0065               Slopes = cms.vdouble(0.982061,0.98068,0.976091,0.978934,0.971368,0.971559,0.964595,0.959935,0.964505,0.961342,0.956929,0.972855,0.974695,0.972261),
0066               Offsets = cms.vdouble(-8.24181e-06,-1.06929e-05,-6.53133e-05,-0.000149074,7.90199e-05,-0.000329831,0.00013365,-0.000472845,-7.92843e-05,-0.000348687,-0.000571311,7.07575e-05,3.87507e-05,-0.000166651),
0067               Pulls = cms.vdouble(1.07221,0.909626,0.861862,0.990638,1.01422,1.01225,1.05827,1.00192,1.10314,0.916937,0.901404,0.932825,0.922407,0.947059) ),
0068     cms.PSet( Method = METHOD_AVGV2,
0069               Slopes = cms.vdouble(0.991046,0.988139,0.992464,0.989286,0.978063,0.970409,0.977013,0.981871,0.976017,0.980569,0.978308,0.977915,0.99055,0.989702),
0070               Offsets = cms.vdouble(0.000179484,-0.000315592,0.000286534,8.59587e-05,0.000123831,-0.000232847,-0.000386919,9.07647e-05,-0.000112936,-1.02089e-05,0.000157863,0.000102052,0.000240302,-8.68151e-05),
0071               Pulls = cms.vdouble(0.99978,1.00856,1.0204,1.17001,1.07818,1.01825,1.01474,0.958657,1.00109,1.09488,0.928755,1.00347,1.10553,1.07716) ),
0072     cms.PSet( Method = METHOD_AVGV3,
0073               Slopes = cms.vdouble(1.0234,1.01731,1.01839,1.01912,1.03934,1.03615,1.01273,1.00858,1.01181,1.01492,1.01467,1.01318,1.01644,1.00849),
0074               Offsets = cms.vdouble(0.000702506,0.000216634,3.31415e-05,0.000210135,0.000578267,0.000550148,0.000451231,-1.4985e-05,0.000224097,0.000587055,0.000725274,0.000455932,0.00118975,0.00021313),
0075               Pulls = cms.vdouble(1.05618,0.912573,1.08878,1.07971,0.957048,0.982786,1.00972,1.01669,0.993294,0.928439,0.980928,0.988378,1.08049,1.16658) ),
0076     cms.PSet( Method = METHOD_RMSV2,
0077               Slopes = cms.vdouble(0.974224,0.974741,0.98057,0.987761,0.952585,0.949676,0.955757,0.959137,0.953516,0.954911,0.959072,0.955053,0.963291,0.96448),
0078               Offsets = cms.vdouble(-0.000242094,-0.000455473,0.000236359,0.000771006,0.000200029,-0.000187993,0.000151929,0.000249159,0.000239323,0.000289015,0.000260618,0.000420144,-0.00024393,-0.000411812),
0079               Pulls = cms.vdouble(1.30967,1.23259,1.14261,1.24064,1.24455,1.28206,1.37354,1.39888,1.26128,1.30592,1.36302,1.23031,1.40948,1.32819) ),
0080     cms.PSet( Method = METHOD_RMSV3,
0081               Slopes = cms.vdouble(1.05309,1.0322,1.02728,1.04035,1.07448,1.0385,1.04547,1.05374,1.04194,1.04169,1.05094,1.04238,1.0321,1.03053),
0082               Offsets = cms.vdouble(0.00105634,-0.000423326,-0.000977181,-9.60822e-05,-0.000328321,-0.00243066,-0.00163576,-0.000720268,-0.00164131,-0.00167024,-0.000967582,-0.00190743,-0.00107932,-0.00120297),
0083               Pulls = cms.vdouble(1.24586,1.20579,1.11604,1.22117,1.42447,1.2723,1.31163,1.40946,1.30212,1.4361,1.32736,1.31096,1.35414,1.40196) )
0084     )
0085
0086 LorentzAngleCalibrations_DeconvolutionModeCOSMIC = cms.VPSet(
0087     cms.PSet( Method = METHOD_WIDTH, Slopes = uncalSlopes , Offsets = uncalOffsets, Pulls = uncalPulls),
0088     cms.PSet( Method = METHOD_PROB1, Slopes = uncalSlopes , Offsets = uncalOffsets, Pulls = uncalPulls),
0089     cms.PSet( Method = METHOD_AVGV2, Slopes = uncalSlopes , Offsets = uncalOffsets, Pulls = uncalPulls),
0090     cms.PSet( Method = METHOD_AVGV3, Slopes = uncalSlopes , Offsets = uncalOffsets, Pulls = uncalPulls),
0091     cms.PSet( Method = METHOD_RMSV2, Slopes = uncalSlopes , Offsets = uncalOffsets, Pulls = uncalPulls),
0092     cms.PSet( Method = METHOD_RMSV3, Slopes = uncalSlopes , Offsets = uncalOffsets, Pulls = uncalPulls)
0093     )