CMS 3D CMS Logo
List of all members | Public Member Functions | Private Attributes
reweight::PoissonMeanShifter Class Reference

#include <LumiReWeighting.h>

Public Member Functions

 PoissonMeanShifter ()
 
 PoissonMeanShifter (float Shift)
 
 PoissonMeanShifter ()
 
 PoissonMeanShifter (float Shift)
 
double ShiftWeight (int ibin)
 
double ShiftWeight (float pvnum)
 
double ShiftWeight (int ibin)
 
double ShiftWeight (float pvnum)
 

Private Attributes

double Pweight_ [25]
 

Detailed Description

Definition at line 32 of file LumiReWeighting.h.

Constructor & Destructor Documentation

reweight::PoissonMeanShifter::PoissonMeanShifter ( )
inline

Definition at line 36 of file LumiReWeighting.h.

36 { };
reweight::PoissonMeanShifter::PoissonMeanShifter ( float  Shift)
inline

Definition at line 38 of file LumiReWeighting.h.

References constexpr, JetChargeProducer_cfi::exp, and Pweight_.

38  {
39 
40  // these are the polynomial or exponential coefficients for each bin of a 25-bin sequence that
41  // convert the Distribution of the 2011 luminosity to something with a lower or higher peak luminosity.
42  // The distributions aren't quite poisson because they model luminosity decreasing during a fill. This implies that
43  // they do get wider as the mean increases, so the weights are not linear with increasing mean.
44 
45  constexpr double p0_minus[20] = { 1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1. };
46  constexpr double p1_minus[20] = {
47  -0.677786,
48  -0.619614,
49  -0.49465,
50  -0.357963,
51  -0.238359,
52  -0.110002,
53  0.0348629,
54  0.191263,
55  0.347648,
56  0.516615,
57  0.679646,
58  0.836673,
59  0.97764,
60  1.135,
61  1.29922,
62  1.42467,
63  1.55901,
64  1.61762,
65  1.67275,
66  1.96008
67  };
68  constexpr double p2_minus[20] = {
69  0.526164,
70  0.251816,
71  0.11049,
72  0.026917,
73  -0.0464692,
74  -0.087022,
75  -0.0931581,
76  -0.0714295,
77  -0.0331772,
78  0.0347473,
79  0.108658,
80  0.193048,
81  0.272314,
82  0.376357,
83  0.4964,
84  0.58854,
85  0.684959,
86  0.731063,
87  0.760044,
88  1.02386
89  };
90 
91  constexpr double p1_expoM[5] = {
92  1.63363e-03,
93  6.79290e-04,
94  3.69900e-04,
95  2.24349e-04,
96  9.87156e-06
97  };
98 
99  constexpr double p2_expoM[5] = {
100  2.64692,
101  3.26585,
102  3.53229,
103  4.18035,
104  5.64027
105  };
106 
107 
108  constexpr double p0_plus[20] = { 1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1. };
109  constexpr double p1_plus[20] = {
110  -0.739059,
111  -0.594445,
112  -0.477276,
113  -0.359707,
114  -0.233573,
115  -0.103458,
116  0.0373401,
117  0.176571,
118  0.337617,
119  0.499074,
120  0.675126,
121  0.840522,
122  1.00917,
123  1.15847,
124  1.23816,
125  1.44271,
126  1.52982,
127  1.46385,
128  1.5802,
129  0.988689
130  };
131  constexpr double p2_plus[20] = {
132  0.208068,
133  0.130033,
134  0.0850356,
135  0.0448344,
136  0.000749832,
137  -0.0331347,
138  -0.0653281,
139  -0.0746009,
140  -0.0800667,
141  -0.0527636,
142  -0.00402649,
143  0.103338,
144  0.261261,
145  0.491084,
146  0.857966,
147  1.19495,
148  1.75071,
149  2.65559,
150  3.35433,
151  5.48835
152  };
153 
154  constexpr double p1_expoP[5] = {
155  1.42463e-01,
156  4.18966e-02,
157  1.12697e-01,
158  1.66197e-01,
159  1.50768e-01
160  };
161 
162  constexpr double p2_expoP[5] = {
163  1.98758,
164  2.27217,
165  2.26799,
166  2.38455,
167  2.52428
168  };
169 
170  // initialize weights based on desired Shift
171 
172 
173 
174  for (int ibin=0;ibin<20;ibin++) {
175 
176  if( Shift < .0) {
177  Pweight_[ibin] = p0_minus[ibin] + p1_minus[ibin]*Shift + p2_minus[ibin]*Shift*Shift;
178  }
179  else {
180  Pweight_[ibin] = p0_plus[ibin] + p1_plus[ibin]*Shift + p2_plus[ibin]*Shift*Shift;
181  }
182  }
183 
184  // last few bins fit better to an exponential...
185 
186  for (int ibin=20;ibin<25;ibin++) {
187  if( Shift < 0.) {
188  Pweight_[ibin] = p1_expoM[ibin-20]*exp(p2_expoM[ibin-20]*Shift);
189  }
190  else {
191  Pweight_[ibin] = p1_expoP[ibin-20]*exp(p2_expoP[ibin-20]*Shift);
192  }
193  }
194 
195  };
constexpr
#define constexpr
JetChargeProducer_cfi.exp
exp
Definition: JetChargeProducer_cfi.py:6
reweight::PoissonMeanShifter::PoissonMeanShifter ( )
inline

Definition at line 43 of file LumiReweightingStandAlone.h.

43 { };
reweight::PoissonMeanShifter::PoissonMeanShifter ( float  Shift)
inline

Definition at line 45 of file LumiReweightingStandAlone.h.

References JetChargeProducer_cfi::exp, and Pweight_.

45  {
46 
47  // these are the polynomial or exponential coefficients for each bin of a 25-bin sequence that
48  // convert the Distribution of the 2011 luminosity to something with a lower or higher peak luminosity.
49  // The distributions aren't quite poisson because they model luminosity decreasing during a fill. This implies that
50  // they do get wider as the mean increases, so the weights are not linear with increasing mean.
51 
52  static const double p0_minus[20] = { 1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1. };
53  static const double p1_minus[20] = {
54  -0.677786,
55  -0.619614,
56  -0.49465,
57  -0.357963,
58  -0.238359,
59  -0.110002,
60  0.0348629,
61  0.191263,
62  0.347648,
63  0.516615,
64  0.679646,
65  0.836673,
66  0.97764,
67  1.135,
68  1.29922,
69  1.42467,
70  1.55901,
71  1.61762,
72  1.67275,
73  1.96008
74  };
75  static const double p2_minus[20] = {
76  0.526164,
77  0.251816,
78  0.11049,
79  0.026917,
80  -0.0464692,
81  -0.087022,
82  -0.0931581,
83  -0.0714295,
84  -0.0331772,
85  0.0347473,
86  0.108658,
87  0.193048,
88  0.272314,
89  0.376357,
90  0.4964,
91  0.58854,
92  0.684959,
93  0.731063,
94  0.760044,
95  1.02386
96  };
97 
98  static const double p1_expoM[5] = {
99  1.63363e-03,
100  6.79290e-04,
101  3.69900e-04,
102  2.24349e-04,
103  9.87156e-06
104  };
105 
106  static const double p2_expoM[5] = {
107  2.64692,
108  3.26585,
109  3.53229,
110  4.18035,
111  5.64027
112  };
113 
114 
115  static const double p0_plus[20] = { 1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1.,1. };
116  static const double p1_plus[20] = {
117  -0.739059,
118  -0.594445,
119  -0.477276,
120  -0.359707,
121  -0.233573,
122  -0.103458,
123  0.0373401,
124  0.176571,
125  0.337617,
126  0.499074,
127  0.675126,
128  0.840522,
129  1.00917,
130  1.15847,
131  1.23816,
132  1.44271,
133  1.52982,
134  1.46385,
135  1.5802,
136  0.988689
137  };
138  static const double p2_plus[20] = {
139  0.208068,
140  0.130033,
141  0.0850356,
142  0.0448344,
143  0.000749832,
144  -0.0331347,
145  -0.0653281,
146  -0.0746009,
147  -0.0800667,
148  -0.0527636,
149  -0.00402649,
150  0.103338,
151  0.261261,
152  0.491084,
153  0.857966,
154  1.19495,
155  1.75071,
156  2.65559,
157  3.35433,
158  5.48835
159  };
160 
161  static const double p1_expoP[5] = {
162  1.42463e-01,
163  4.18966e-02,
164  1.12697e-01,
165  1.66197e-01,
166  1.50768e-01
167  };
168 
169  static const double p2_expoP[5] = {
170  1.98758,
171  2.27217,
172  2.26799,
173  2.38455,
174  2.52428
175  };
176 
177  // initialize weights based on desired Shift
178 
179 
180 
181  for (int ibin=0;ibin<20;ibin++) {
182 
183  if( Shift < .0) {
184  Pweight_[ibin] = p0_minus[ibin] + p1_minus[ibin]*Shift + p2_minus[ibin]*Shift*Shift;
185  }
186  else {
187  Pweight_[ibin] = p0_plus[ibin] + p1_plus[ibin]*Shift + p2_plus[ibin]*Shift*Shift;
188  }
189  }
190 
191  // last few bins fit better to an exponential...
192 
193  for (int ibin=20;ibin<25;ibin++) {
194  if( Shift < 0.) {
195  Pweight_[ibin] = p1_expoM[ibin-20]*exp(p2_expoM[ibin-20]*Shift);
196  }
197  else {
198  Pweight_[ibin] = p1_expoP[ibin-20]*exp(p2_expoP[ibin-20]*Shift);
199  }
200  }
201 
202  };
JetChargeProducer_cfi.exp
exp
Definition: JetChargeProducer_cfi.py:6

Member Function Documentation

double reweight::PoissonMeanShifter::ShiftWeight ( int  ibin)
inline

Definition at line 197 of file LumiReWeighting.h.

References Pweight_.

Referenced by ShiftWeight().

197  {
198  if(ibin<25 && ibin>=0) { return Pweight_[ibin]; }
199  else { return 0;}
200 
201  };
double reweight::PoissonMeanShifter::ShiftWeight ( float  pvnum)
inline

Definition at line 203 of file LumiReWeighting.h.

References ShiftWeight().

203  {
204  return ShiftWeight(int(pvnum));
205  };
double reweight::PoissonMeanShifter::ShiftWeight ( int  ibin)
inline

Definition at line 204 of file LumiReweightingStandAlone.h.

References Pweight_.

204  {
205 
206  if(ibin<25 && ibin>=0) { return Pweight_[ibin]; }
207  else { return 0;}
208 
209  };
double reweight::PoissonMeanShifter::ShiftWeight ( float  pvnum)
inline

Definition at line 211 of file LumiReweightingStandAlone.h.

References createfilelist::int, and Pweight_.

211  {
212 
213  int ibin = int(pvnum);
214 
215  if(ibin<25 && ibin>=0) { return Pweight_[ibin]; }
216  else { return 0;}
217 
218  };

Member Data Documentation

double reweight::PoissonMeanShifter::Pweight_
private

Definition at line 205 of file LumiReWeighting.h.

Referenced by PoissonMeanShifter(), and ShiftWeight().

Generated for CMSSW Reference Manual by doxygen 1.8.11