Utils/src/PUReweighting.cc

// $Id: PUReweighting.cc,v 1.1 2011/05/26 12:00:56 mzanetti Exp $


#include "MitPhysics/Utils/interface/PUReweighting.h"
#include <TFile.h>

ClassImp(mithep::PUReweighting)

using namespace mithep;


PUReweighting::PUReweighting(const char* fileName, const char*histoName) {

  TFile *referenceFile = TFile::Open(fileName);
  if (!referenceFile) {
    std::cout<<"Error: PU reference file does not exist, reweighting will not be applied"<<std::endl;
  }
  else {
    referenceHisto = (TH1F*) referenceFile->Get(histoName);
  }

}

double PUReweighting::reweightOOT(int inTimePUMultiplicity, int outOfTimePUMultiplicity) {

    if (!referenceHisto) {
      std::cout<<"Error: PU reference histo does not exist in the file, reweighting will not be applied"<<std::endl;
      return 1;
    }

  // residual correction to match lumi spectrum
  static double Correct_Weights2011[25] = { 
    5.30031,
    2.07903,
    1.40729,
    1.27687,
    1.0702,
    0.902094,
    0.902345,
    0.931449,
    0.78202,
    0.824686,
    0.837735,
    0.910261,
    1.01394,
    1.1599,
    1.12778,
    1.58423,
    1.78868,
    1.58296,
    2.3291,
    3.86641,
    0,
    0,
    0,
    0,
    0
  };                        


  static double weight_24[25] = {
    0,
    0,
    0,
    0,
    2.46277e-06,
    2.95532e-05,
    0.000104668,
    0.000401431,
    0.00130034,
    0.00342202,
    0.00818132,
    0.0175534,
    0.035784,
    0.0650836,
    0.112232,
    0.178699,
    0.268934,
    0.380868,
    0.507505,
    0.640922,
    0.768551,
    0.877829,
    0.958624,
    0.99939,
    1
  };

  static double weight_23[25] = {
    0,
    1.20628e-06,
    1.20628e-06,
    2.41255e-06,
    1.20628e-05,
    6.39326e-05,
    0.000252112,
    0.000862487,
    0.00244995,
    0.00616527,
    0.0140821,
    0.0293342,
    0.0564501,
    0.100602,
    0.164479,
    0.252659,
    0.36268,
    0.491427,
    0.627979,
    0.75918,
    0.873185,
    0.957934,
    0.999381,
    1,
    0.957738
  };

  static double weight_22[25] = {
    0,
    0,
    0,
    5.88636e-06,
    3.0609e-05,
    0.000143627,
    0.000561558,
    0.00173059,
    0.00460078,
    0.0110616,
    0.0238974,
    0.0475406,
    0.0875077,
    0.148682,
    0.235752,
    0.343591,
    0.473146,
    0.611897,
    0.748345,
    0.865978,
    0.953199,
    0.997848,
    1,
    0.954245,
    0.873688
  };

  static double weight_21[25] = {
    0,
    0,
    1.15381e-06,
    8.07665e-06,
    7.1536e-05,
    0.000280375,
    0.00107189,
    0.00327104,
    0.00809396,
    0.0190978,
    0.0401894,
    0.0761028,
    0.13472,
    0.216315,
    0.324649,
    0.455125,
    0.598241,
    0.739215,
    0.861866,
    0.953911,
    0.998918,
    1,
    0.956683,
    0.872272,
    0.76399
  };
 
 
  static double weight_20[25] = {
    0,
    0,
    1.12532e-06,
    2.58822e-05,
    0.000145166,
    0.000633552,
    0.00215048,
    0.00592816,
    0.0145605,
    0.0328367,
    0.0652649,
    0.11893,
    0.19803,
    0.305525,
    0.436588,
    0.581566,
    0.727048,
    0.8534,
    0.949419,
    0.999785,
    1,
    0.953008,
    0.865689,
    0.753288,
    0.62765
  }; 
  static double weight_19[25] = {
    0,
    0,
    1.20714e-05,
    5.92596e-05,
    0.000364337,
    0.00124994,
    0.00403953,
    0.0108149,
    0.025824,
    0.0544969,
    0.103567,
    0.17936,
    0.283532,
    0.416091,
    0.562078,
    0.714714,
    0.846523,
    0.947875,
    1,
    0.999448,
    0.951404,
    0.859717,
    0.742319,
    0.613601,
    0.48552
  };

  static double weight_18[25] = {
    0,
    3.20101e-06,
    2.88091e-05,
    0.000164319,
    0.000719161,
    0.00250106,
    0.00773685,
    0.0197513,
    0.0443693,
    0.0885998,
    0.159891,
    0.262607,
    0.392327,
    0.543125,
    0.69924,
    0.837474,
    0.943486,
    0.998029,
    1,
    0.945937,
    0.851807,
    0.729309,
    0.596332,
    0.467818,
    0.350434
  };

 
  static double weight_17[25] = {
    1.03634e-06,
    7.25437e-06,
    4.97443e-05,
    0.000340956,
    0.00148715,
    0.00501485,
    0.0143067,
    0.034679,
    0.0742009,
    0.140287,
    0.238288,
    0.369416,
    0.521637,
    0.682368,
    0.828634,
    0.939655,
    1,
    0.996829,
    0.94062,
    0.841575,
    0.716664,
    0.582053,
    0.449595,
    0.331336,
    0.234332
  };

 
  static double weight_16[25] = {
    4.03159e-06,
    2.41895e-05,
    0.000141106,
    0.00081942,
    0.00314565,
    0.00990662,
    0.026293,
    0.0603881,
    0.120973,
    0.214532,
    0.343708,
    0.501141,
    0.665978,
    0.820107,
    0.938149,
    1,
    0.99941,
    0.940768,
    0.837813,
    0.703086,
    0.564023,
    0.42928,
    0.312515,
    0.216251,
    0.14561
  };
 
 
  static double weight_15[25] = {
    9.76084e-07,
    5.07564e-05,
    0.000303562,
    0.00174036,
    0.00617959,
    0.0188579,
    0.047465,
    0.101656,
    0.189492,
    0.315673,
    0.474383,
    0.646828,
    0.809462,
    0.934107,
    0.998874,
    1,
    0.936163,
    0.827473,
    0.689675,
    0.544384,
    0.40907,
    0.290648,
    0.198861,
    0.12951,
    0.0808051
  };
 
 
  static double weight_14[25] = {
    1.13288e-05,
    0.000124617,
    0.000753365,
    0.00345056,
    0.0123909,
    0.0352712,
    0.0825463,
    0.16413,
    0.287213,
    0.44615,
    0.625826,
    0.796365,
    0.930624,
    0.999958,
    1,
    0.934414,
    0.816456,
    0.672939,
    0.523033,
    0.386068,
    0.269824,
    0.180342,
    0.114669,
    0.0698288,
    0.0406496
  };

 
  static double weight_13[25] = {
    2.54296e-05,
    0.000261561,
    0.00167018,
    0.00748083,
    0.0241308,
    0.0636801,
    0.138222,
    0.255814,
    0.414275,
    0.600244,
    0.779958,
    0.92256,
    0.999155,
    1,
    0.927126,
    0.804504,
    0.651803,
    0.497534,
    0.35976,
    0.245834,
    0.160904,
    0.0991589,
    0.0585434,
    0.0332437,
    0.0180159
  };

  static double weight_12[25] = {
    5.85742e-05,
    0.000627706,
    0.00386677,
    0.0154068,
    0.0465892,
    0.111683,
    0.222487,
    0.381677,
    0.5719,
    0.765001,
    0.915916,
    1,
    0.999717,
    0.921443,
    0.791958,
    0.632344,
    0.475195,
    0.334982,
    0.223666,
    0.141781,
    0.0851538,
    0.048433,
    0.0263287,
    0.0133969,
    0.00696683
  };

 
  static double weight_11[25] = {
    0.00015238,
    0.00156064,
    0.00846044,
    0.0310939,
    0.0856225,
    0.187589,
    0.343579,
    0.541892,
    0.74224,
    0.909269,
    0.998711,
    1,
    0.916889,
    0.77485,
    0.608819,
    0.447016,
    0.307375,
    0.198444,
    0.121208,
    0.070222,
    0.0386492,
    0.0201108,
    0.0100922,
    0.00484937,
    0.00222458
  };

  static double weight_10[25] = {
    0.000393044,
    0.00367001,
    0.0179474,
    0.060389,
    0.151477,
    0.302077,
    0.503113,
    0.720373,
    0.899568,
    1,
    0.997739,
    0.909409,
    0.75728,
    0.582031,
    0.415322,
    0.277663,
    0.174147,
    0.102154,
    0.0566719,
    0.0298642,
    0.0147751,
    0.00710995,
    0.00319628,
    0.00140601,
    0.000568796
  };

 
  static double weight_9[25] = {
    0.00093396,
    0.00854448,
    0.0380306,
    0.113181,
    0.256614,
    0.460894,
    0.690242,
    0.888781,
    1,
    0.998756,
    0.899872,
    0.735642,
    0.552532,
    0.382726,
    0.246114,
    0.147497,
    0.0825541,
    0.0441199,
    0.0218157,
    0.0103578,
    0.00462959,
    0.0019142,
    0.000771598,
    0.000295893,
    0.000111529
  };

 
  static double weight_8[25] = {
    0.00240233,
    0.0192688,
    0.0768653,
    0.205008,
    0.410958,
    0.65758,
    0.875657,
    0.999886,
    1,
    0.889476,
    0.711446,
    0.517781,
    0.345774,
    0.212028,
    0.121208,
    0.0644629,
    0.0324928,
    0.0152492,
    0.00673527,
    0.0028547,
    0.00117213,
    0.000440177,
    0.000168471,
    5.80689e-05,
    1.93563e-05
  };

  static double weight_7[25] = {
    0.00617233,
    0.0428714,
    0.150018,
    0.350317,
    0.612535,
    0.856525,
    0.999923,
    1,
    0.87544,
    0.679383,
    0.478345,
    0.303378,
    0.176923,
    0.0950103,
    0.0476253,
    0.0222211,
    0.00972738,
    0.00392962,
    0.0015258,
    0.000559168,
    0.000183928,
    6.77983e-05,
    1.67818e-05,
    7.38398e-06,
    6.71271e-07
  };
 
  static double weight_6[25] = {
    0.0154465,
    0.0923472,
    0.277322,
    0.55552,
    0.833099,
    0.999035,
    1,
    0.855183,
    0.641976,
    0.428277,
    0.256804,
    0.139798,
    0.0700072,
    0.0321586,
    0.0137971,
    0.00544756,
    0.00202316,
    0.000766228,
    0.000259348,
    8.45836e-05,
    1.80362e-05,
    8.70713e-06,
    3.73163e-06,
    6.21938e-07,
    0
  };
 
 
  static double weight_5[25] = {
    0.0382845,
    0.191122,
    0.478782,
    0.797314,
    1,
    0.997148,
    0.831144,
    0.59461,
    0.371293,
    0.205903,
    0.103102,
    0.0471424,
    0.0194997,
    0.00749415,
    0.00273709,
    0.000879189,
    0.000286049,
    0.000102364,
    1.70606e-05,
    3.98081e-06,
    2.27475e-06,
    0,
    0,
    0,
    0
  };
 
 
  static double weight_4[25] = {
    0.0941305,
    0.373824,
    0.750094,
    1,
    0.997698,
    0.800956,
    0.532306,
    0.304597,
    0.152207,
    0.0676275,
    0.0270646,
    0.00975365,
    0.00326077,
    0.00101071,
    0.000301781,
    7.41664e-05,
    1.58563e-05,
    3.58045e-06,
    1.02299e-06,
    0,
    5.11493e-07,
    0,
    0,
    0,
    0
  };
 
 
  static double weight_3[25] = {
    0.222714,
    0.667015,
    1,
    0.999208,
    0.750609,
    0.449854,
    0.224968,
    0.0965185,
    0.0361225,
    0.012084,
    0.00359618,
    0.000977166,
    0.000239269,
    6.29422e-05,
    1.16064e-05,
    1.78559e-06,
    0,
    4.46398e-07,
    0,
    0,
    0,
    0,
    0,
    0,
    0
  };
 
  static double weight_2[25] = {
    0.499541,
    0.999607,
    1,
    0.666607,
    0.333301,
    0.13279,
    0.0441871,
    0.0127455,
    0.00318434,
    0.00071752,
    0.000132204,
    2.69578e-05,
    5.16999e-06,
    2.21571e-06,
    0,
    0,
    0,
    0,
    0,
    0,
    0,
    0,
    0,
    0,
    0
  };
 
  static double weight_1[25] = {
    0.999165,
    1,
    0.499996,
    0.166868,
    0.0414266,
    0.00831053,
    0.00137472,
    0.000198911,
    2.66302e-05,
    2.44563e-06,
    2.71737e-07,
    2.71737e-07,
    0,
    0,
    0,
    0,
    0,
    0,
    0,
    0,
    0,
    0,
    0,
    0,
    0
  };
 
  static double weight_0[25] = {
    1,
    0,
    0,
    0,
    0,
    0,
    0,
    0,
    0,
    0,
    0,
    0,
    0,
    0,
    0,
    0,
    0,
    0,
    0,
    0,
    0,
    0,
    0,
    0,
    0
  };

  double WeightOOTPU[25][25];

  double* WeightPtr = 0;

  for(int iint = 0; iint<25; ++iint){
    if(iint ==0) WeightPtr = weight_0;
    if(iint ==1) WeightPtr = weight_1;
    if(iint ==2) WeightPtr = weight_2;
    if(iint ==3) WeightPtr = weight_3;
    if(iint ==4) WeightPtr = weight_4;
    if(iint ==5) WeightPtr = weight_5;
    if(iint ==6) WeightPtr = weight_6;
    if(iint ==7) WeightPtr = weight_7;
    if(iint ==8) WeightPtr = weight_8;
    if(iint ==9) WeightPtr = weight_9;
    if(iint ==10) WeightPtr = weight_10;
    if(iint ==11) WeightPtr = weight_11;
    if(iint ==12) WeightPtr = weight_12;
    if(iint ==13) WeightPtr = weight_13;
    if(iint ==14) WeightPtr = weight_14;
    if(iint ==15) WeightPtr = weight_15;
    if(iint ==16) WeightPtr = weight_16;
    if(iint ==17) WeightPtr = weight_17;
    if(iint ==18) WeightPtr = weight_18;
    if(iint ==19) WeightPtr = weight_19;
    if(iint ==20) WeightPtr = weight_20;
    if(iint ==21) WeightPtr = weight_21;
    if(iint ==22) WeightPtr = weight_22;
    if(iint ==23) WeightPtr = weight_23;
    if(iint ==24) WeightPtr = weight_24;

    for(int ibin = 0; ibin<25; ++ibin){
      WeightOOTPU[iint][ibin] = *(WeightPtr+ibin);
    }
  }


  double inTimeWeight = referenceHisto->GetBinContent( inTimePUMultiplicity );

  return inTimeWeight * WeightOOTPU[inTimePUMultiplicity-1][outOfTimePUMultiplicity-1] * Correct_Weights2011[inTimePUMultiplicity-1];

}

Revision:	1.1
Committed:	Wed Jun 1 12:36:59 2011 UTC (13 years, 11 months ago) by mzanetti
Content type:	text/plain
Branch:	MAIN
CVS Tags:	Mit_029c, Mit_029b, Mit_029a, Mit_028a, Mit_028, Mit_027, Mit_027a, Mit_025e, Mit_025d, Mit_025c, Mit_025b, Mit_025a, Mit_025, Mit_025pre2, Mit_024b, Mit_025pre1, Mit_024a, Mit_024, Mit_023, Mit_022a, Mit_022, HEAD
Log Message:	typo in the name of the previous file
#	Content
1	// $Id: PUReweighting.cc,v 1.1 2011/05/26 12:00:56 mzanetti Exp $
2
3
4	#include "MitPhysics/Utils/interface/PUReweighting.h"
5	#include <TFile.h>
6
7	ClassImp(mithep::PUReweighting)
8
9	using namespace mithep;
10
11
12	PUReweighting::PUReweighting(const char* fileName, const char*histoName) {
13
14	TFile *referenceFile = TFile::Open(fileName);
15	if (!referenceFile) {
16	std::cout<<"Error: PU reference file does not exist, reweighting will not be applied"<<std::endl;
17	}
18	else {
19	referenceHisto = (TH1F*) referenceFile->Get(histoName);
20	}
21
22	}
23
24	double PUReweighting::reweightOOT(int inTimePUMultiplicity, int outOfTimePUMultiplicity) {
25
26	if (!referenceHisto) {
27	std::cout<<"Error: PU reference histo does not exist in the file, reweighting will not be applied"<<std::endl;
28	return 1;
29	}
30
31	// residual correction to match lumi spectrum
32	static double Correct_Weights2011[25] = {
33	5.30031,
34	2.07903,
35	1.40729,
36	1.27687,
37	1.0702,
38	0.902094,
39	0.902345,
40	0.931449,
41	0.78202,
42	0.824686,
43	0.837735,
44	0.910261,
45	1.01394,
46	1.1599,
47	1.12778,
48	1.58423,
49	1.78868,
50	1.58296,
51	2.3291,
52	3.86641,
53	0,
54	0,
55	0,
56	0,
57	0
58	};
59
60
61	static double weight_24[25] = {
62	0,
63	0,
64	0,
65	0,
66	2.46277e-06,
67	2.95532e-05,
68	0.000104668,
69	0.000401431,
70	0.00130034,
71	0.00342202,
72	0.00818132,
73	0.0175534,
74	0.035784,
75	0.0650836,
76	0.112232,
77	0.178699,
78	0.268934,
79	0.380868,
80	0.507505,
81	0.640922,
82	0.768551,
83	0.877829,
84	0.958624,
85	0.99939,
86	1
87	};
88
89	static double weight_23[25] = {
90	0,
91	1.20628e-06,
92	1.20628e-06,
93	2.41255e-06,
94	1.20628e-05,
95	6.39326e-05,
96	0.000252112,
97	0.000862487,
98	0.00244995,
99	0.00616527,
100	0.0140821,
101	0.0293342,
102	0.0564501,
103	0.100602,
104	0.164479,
105	0.252659,
106	0.36268,
107	0.491427,
108	0.627979,
109	0.75918,
110	0.873185,
111	0.957934,
112	0.999381,
113	1,
114	0.957738
115	};
116
117	static double weight_22[25] = {
118	0,
119	0,
120	0,
121	5.88636e-06,
122	3.0609e-05,
123	0.000143627,
124	0.000561558,
125	0.00173059,
126	0.00460078,
127	0.0110616,
128	0.0238974,
129	0.0475406,
130	0.0875077,
131	0.148682,
132	0.235752,
133	0.343591,
134	0.473146,
135	0.611897,
136	0.748345,
137	0.865978,
138	0.953199,
139	0.997848,
140	1,
141	0.954245,
142	0.873688
143	};
144
145	static double weight_21[25] = {
146	0,
147	0,
148	1.15381e-06,
149	8.07665e-06,
150	7.1536e-05,
151	0.000280375,
152	0.00107189,
153	0.00327104,
154	0.00809396,
155	0.0190978,
156	0.0401894,
157	0.0761028,
158	0.13472,
159	0.216315,
160	0.324649,
161	0.455125,
162	0.598241,
163	0.739215,
164	0.861866,
165	0.953911,
166	0.998918,
167	1,
168	0.956683,
169	0.872272,
170	0.76399
171	};
172
173
174	static double weight_20[25] = {
175	0,
176	0,
177	1.12532e-06,
178	2.58822e-05,
179	0.000145166,
180	0.000633552,
181	0.00215048,
182	0.00592816,
183	0.0145605,
184	0.0328367,
185	0.0652649,
186	0.11893,
187	0.19803,
188	0.305525,
189	0.436588,
190	0.581566,
191	0.727048,
192	0.8534,
193	0.949419,
194	0.999785,
195	1,
196	0.953008,
197	0.865689,
198	0.753288,
199	0.62765
200	};
201	static double weight_19[25] = {
202	0,
203	0,
204	1.20714e-05,
205	5.92596e-05,
206	0.000364337,
207	0.00124994,
208	0.00403953,
209	0.0108149,
210	0.025824,
211	0.0544969,
212	0.103567,
213	0.17936,
214	0.283532,
215	0.416091,
216	0.562078,
217	0.714714,
218	0.846523,
219	0.947875,
220	1,
221	0.999448,
222	0.951404,
223	0.859717,
224	0.742319,
225	0.613601,
226	0.48552
227	};
228
229	static double weight_18[25] = {
230	0,
231	3.20101e-06,
232	2.88091e-05,
233	0.000164319,
234	0.000719161,
235	0.00250106,
236	0.00773685,
237	0.0197513,
238	0.0443693,
239	0.0885998,
240	0.159891,
241	0.262607,
242	0.392327,
243	0.543125,
244	0.69924,
245	0.837474,
246	0.943486,
247	0.998029,
248	1,
249	0.945937,
250	0.851807,
251	0.729309,
252	0.596332,
253	0.467818,
254	0.350434
255	};
256
257
258	static double weight_17[25] = {
259	1.03634e-06,
260	7.25437e-06,
261	4.97443e-05,
262	0.000340956,
263	0.00148715,
264	0.00501485,
265	0.0143067,
266	0.034679,
267	0.0742009,
268	0.140287,
269	0.238288,
270	0.369416,
271	0.521637,
272	0.682368,
273	0.828634,
274	0.939655,
275	1,
276	0.996829,
277	0.94062,
278	0.841575,
279	0.716664,
280	0.582053,
281	0.449595,
282	0.331336,
283	0.234332
284	};
285
286
287	static double weight_16[25] = {
288	4.03159e-06,
289	2.41895e-05,
290	0.000141106,
291	0.00081942,
292	0.00314565,
293	0.00990662,
294	0.026293,
295	0.0603881,
296	0.120973,
297	0.214532,
298	0.343708,
299	0.501141,
300	0.665978,
301	0.820107,
302	0.938149,
303	1,
304	0.99941,
305	0.940768,
306	0.837813,
307	0.703086,
308	0.564023,
309	0.42928,
310	0.312515,
311	0.216251,
312	0.14561
313	};
314
315
316	static double weight_15[25] = {
317	9.76084e-07,
318	5.07564e-05,
319	0.000303562,
320	0.00174036,
321	0.00617959,
322	0.0188579,
323	0.047465,
324	0.101656,
325	0.189492,
326	0.315673,
327	0.474383,
328	0.646828,
329	0.809462,
330	0.934107,
331	0.998874,
332	1,
333	0.936163,
334	0.827473,
335	0.689675,
336	0.544384,
337	0.40907,
338	0.290648,
339	0.198861,
340	0.12951,
341	0.0808051
342	};
343
344
345	static double weight_14[25] = {
346	1.13288e-05,
347	0.000124617,
348	0.000753365,
349	0.00345056,
350	0.0123909,
351	0.0352712,
352	0.0825463,
353	0.16413,
354	0.287213,
355	0.44615,
356	0.625826,
357	0.796365,
358	0.930624,
359	0.999958,
360	1,
361	0.934414,
362	0.816456,
363	0.672939,
364	0.523033,
365	0.386068,
366	0.269824,
367	0.180342,
368	0.114669,
369	0.0698288,
370	0.0406496
371	};
372
373
374	static double weight_13[25] = {
375	2.54296e-05,
376	0.000261561,
377	0.00167018,
378	0.00748083,
379	0.0241308,
380	0.0636801,
381	0.138222,
382	0.255814,
383	0.414275,
384	0.600244,
385	0.779958,
386	0.92256,
387	0.999155,
388	1,
389	0.927126,
390	0.804504,
391	0.651803,
392	0.497534,
393	0.35976,
394	0.245834,
395	0.160904,
396	0.0991589,
397	0.0585434,
398	0.0332437,
399	0.0180159
400	};
401
402	static double weight_12[25] = {
403	5.85742e-05,
404	0.000627706,
405	0.00386677,
406	0.0154068,
407	0.0465892,
408	0.111683,
409	0.222487,
410	0.381677,
411	0.5719,
412	0.765001,
413	0.915916,
414	1,
415	0.999717,
416	0.921443,
417	0.791958,
418	0.632344,
419	0.475195,
420	0.334982,
421	0.223666,
422	0.141781,
423	0.0851538,
424	0.048433,
425	0.0263287,
426	0.0133969,
427	0.00696683
428	};
429
430
431	static double weight_11[25] = {
432	0.00015238,
433	0.00156064,
434	0.00846044,
435	0.0310939,
436	0.0856225,
437	0.187589,
438	0.343579,
439	0.541892,
440	0.74224,
441	0.909269,
442	0.998711,
443	1,
444	0.916889,
445	0.77485,
446	0.608819,
447	0.447016,
448	0.307375,
449	0.198444,
450	0.121208,
451	0.070222,
452	0.0386492,
453	0.0201108,
454	0.0100922,
455	0.00484937,
456	0.00222458
457	};
458
459	static double weight_10[25] = {
460	0.000393044,
461	0.00367001,
462	0.0179474,
463	0.060389,
464	0.151477,
465	0.302077,
466	0.503113,
467	0.720373,
468	0.899568,
469	1,
470	0.997739,
471	0.909409,
472	0.75728,
473	0.582031,
474	0.415322,
475	0.277663,
476	0.174147,
477	0.102154,
478	0.0566719,
479	0.0298642,
480	0.0147751,
481	0.00710995,
482	0.00319628,
483	0.00140601,
484	0.000568796
485	};
486
487
488	static double weight_9[25] = {
489	0.00093396,
490	0.00854448,
491	0.0380306,
492	0.113181,
493	0.256614,
494	0.460894,
495	0.690242,
496	0.888781,
497	1,
498	0.998756,
499	0.899872,
500	0.735642,
501	0.552532,
502	0.382726,
503	0.246114,
504	0.147497,
505	0.0825541,
506	0.0441199,
507	0.0218157,
508	0.0103578,
509	0.00462959,
510	0.0019142,
511	0.000771598,
512	0.000295893,
513	0.000111529
514	};
515
516
517	static double weight_8[25] = {
518	0.00240233,
519	0.0192688,
520	0.0768653,
521	0.205008,
522	0.410958,
523	0.65758,
524	0.875657,
525	0.999886,
526	1,
527	0.889476,
528	0.711446,
529	0.517781,
530	0.345774,
531	0.212028,
532	0.121208,
533	0.0644629,
534	0.0324928,
535	0.0152492,
536	0.00673527,
537	0.0028547,
538	0.00117213,
539	0.000440177,
540	0.000168471,
541	5.80689e-05,
542	1.93563e-05
543	};
544
545	static double weight_7[25] = {
546	0.00617233,
547	0.0428714,
548	0.150018,
549	0.350317,
550	0.612535,
551	0.856525,
552	0.999923,
553	1,
554	0.87544,
555	0.679383,
556	0.478345,
557	0.303378,
558	0.176923,
559	0.0950103,
560	0.0476253,
561	0.0222211,
562	0.00972738,
563	0.00392962,
564	0.0015258,
565	0.000559168,
566	0.000183928,
567	6.77983e-05,
568	1.67818e-05,
569	7.38398e-06,
570	6.71271e-07
571	};
572
573	static double weight_6[25] = {
574	0.0154465,
575	0.0923472,
576	0.277322,
577	0.55552,
578	0.833099,
579	0.999035,
580	1,
581	0.855183,
582	0.641976,
583	0.428277,
584	0.256804,
585	0.139798,
586	0.0700072,
587	0.0321586,
588	0.0137971,
589	0.00544756,
590	0.00202316,
591	0.000766228,
592	0.000259348,
593	8.45836e-05,
594	1.80362e-05,
595	8.70713e-06,
596	3.73163e-06,
597	6.21938e-07,
598	0
599	};
600
601
602	static double weight_5[25] = {
603	0.0382845,
604	0.191122,
605	0.478782,
606	0.797314,
607	1,
608	0.997148,
609	0.831144,
610	0.59461,
611	0.371293,
612	0.205903,
613	0.103102,
614	0.0471424,
615	0.0194997,
616	0.00749415,
617	0.00273709,
618	0.000879189,
619	0.000286049,
620	0.000102364,
621	1.70606e-05,
622	3.98081e-06,
623	2.27475e-06,
624	0,
625	0,
626	0,
627	0
628	};
629
630
631	static double weight_4[25] = {
632	0.0941305,
633	0.373824,
634	0.750094,
635	1,
636	0.997698,
637	0.800956,
638	0.532306,
639	0.304597,
640	0.152207,
641	0.0676275,
642	0.0270646,
643	0.00975365,
644	0.00326077,
645	0.00101071,
646	0.000301781,
647	7.41664e-05,
648	1.58563e-05,
649	3.58045e-06,
650	1.02299e-06,
651	0,
652	5.11493e-07,
653	0,
654	0,
655	0,
656	0
657	};
658
659
660	static double weight_3[25] = {
661	0.222714,
662	0.667015,
663	1,
664	0.999208,
665	0.750609,
666	0.449854,
667	0.224968,
668	0.0965185,
669	0.0361225,
670	0.012084,
671	0.00359618,
672	0.000977166,
673	0.000239269,
674	6.29422e-05,
675	1.16064e-05,
676	1.78559e-06,
677	0,
678	4.46398e-07,
679	0,
680	0,
681	0,
682	0,
683	0,
684	0,
685	0
686	};
687
688	static double weight_2[25] = {
689	0.499541,
690	0.999607,
691	1,
692	0.666607,
693	0.333301,
694	0.13279,
695	0.0441871,
696	0.0127455,
697	0.00318434,
698	0.00071752,
699	0.000132204,
700	2.69578e-05,
701	5.16999e-06,
702	2.21571e-06,
703	0,
704	0,
705	0,
706	0,
707	0,
708	0,
709	0,
710	0,
711	0,
712	0,
713	0
714	};
715
716	static double weight_1[25] = {
717	0.999165,
718	1,
719	0.499996,
720	0.166868,
721	0.0414266,
722	0.00831053,
723	0.00137472,
724	0.000198911,
725	2.66302e-05,
726	2.44563e-06,
727	2.71737e-07,
728	2.71737e-07,
729	0,
730	0,
731	0,
732	0,
733	0,
734	0,
735	0,
736	0,
737	0,
738	0,
739	0,
740	0,
741	0
742	};
743
744	static double weight_0[25] = {
745	1,
746	0,
747	0,
748	0,
749	0,
750	0,
751	0,
752	0,
753	0,
754	0,
755	0,
756	0,
757	0,
758	0,
759	0,
760	0,
761	0,
762	0,
763	0,
764	0,
765	0,
766	0,
767	0,
768	0,
769	0
770	};
771
772	double WeightOOTPU[25][25];
773
774	double* WeightPtr = 0;
775
776	for(int iint = 0; iint<25; ++iint){
777	if(iint ==0) WeightPtr = weight_0;
778	if(iint ==1) WeightPtr = weight_1;
779	if(iint ==2) WeightPtr = weight_2;
780	if(iint ==3) WeightPtr = weight_3;
781	if(iint ==4) WeightPtr = weight_4;
782	if(iint ==5) WeightPtr = weight_5;
783	if(iint ==6) WeightPtr = weight_6;
784	if(iint ==7) WeightPtr = weight_7;
785	if(iint ==8) WeightPtr = weight_8;
786	if(iint ==9) WeightPtr = weight_9;
787	if(iint ==10) WeightPtr = weight_10;
788	if(iint ==11) WeightPtr = weight_11;
789	if(iint ==12) WeightPtr = weight_12;
790	if(iint ==13) WeightPtr = weight_13;
791	if(iint ==14) WeightPtr = weight_14;
792	if(iint ==15) WeightPtr = weight_15;
793	if(iint ==16) WeightPtr = weight_16;
794	if(iint ==17) WeightPtr = weight_17;
795	if(iint ==18) WeightPtr = weight_18;
796	if(iint ==19) WeightPtr = weight_19;
797	if(iint ==20) WeightPtr = weight_20;
798	if(iint ==21) WeightPtr = weight_21;
799	if(iint ==22) WeightPtr = weight_22;
800	if(iint ==23) WeightPtr = weight_23;
801	if(iint ==24) WeightPtr = weight_24;
802
803	for(int ibin = 0; ibin<25; ++ibin){
804	WeightOOTPU[iint][ibin] = *(WeightPtr+ibin);
805	}
806	}
807
808
809	double inTimeWeight = referenceHisto->GetBinContent( inTimePUMultiplicity );
810
811	return inTimeWeight * WeightOOTPU[inTimePUMultiplicity-1][outOfTimePUMultiplicity-1] * Correct_Weights2011[inTimePUMultiplicity-1];
812
813	}
814