CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
SingleParticleJetResponse.cc
Go to the documentation of this file.
1 #include "JetMETCorrections/Utilities/interface/SingleParticleJetResponse.h"
2 
3 using namespace std;
4 
5 SingleParticleJetResponse::SingleParticleJetResponse() {}
6 
7 vector<double> SingleParticleJetResponse::response(double echar, double energycluster, int algo) const
8 {
9  vector<double> response;
10  double recal=0,rhcal=0;
11 
12  // threshold on energy of track
13  double e_thEC_hint;
14  double e_thHC_hint;
15  double e_thHC1_mip;
16  double e_thHC2_mip;
17 
18  // parameters for Ecal responces with interaction in Ecal
19  double PAR1_EC_hint;
20  double PAR2_EC_hint;
21  double PAR3_EC_hint;
22  double PAR4_EC_hint;
23 
24  // parameters for Hcal responces with interaction in Ecal
25  double PAR1_HC_hint;
26  double PAR2_HC_hint;
27  double PAR3_HC_hint;
28  double PAR4_HC_hint;
29 
30  // parameters for Hcal responces without interaction in Ecal (MIP)
31  double PAR1_HC_mip;
32  double PAR2_HC_mip;
33  double PAR3_HC_mip;
34  double PAR4_HC_mip;
35  double PAR5_HC_mip;
36 
37  if (algo == 0)
38  {
39  if(energycluster>0.5)
40  {
41  double epiecal=1.6/(1+0.6*0.11*log(0.4*echar));
42  double epihcal=1.39/(1+0.39*0.11*log(0.6*echar));
43  recal=0.4*echar/epiecal;
44  rhcal=0.6*echar/epihcal;
45  }
46  else
47  {
48  double epihcal=1.39/(1+0.39*0.11*log(echar));
49  rhcal=echar/epihcal;
50  recal=energycluster;
51  } // interact or not
52  }
53 
54  else if (algo == 1) // ORCA_6
55  {
56  // threshold on energy of track
57  e_thHC_hint = 65.;
58  e_thHC1_mip = 14.;
59  e_thHC2_mip = 54.;
60 
61  // parameters for Ecal responses with interaction in Ecal
62  PAR1_EC_hint = 3.2181;
63  PAR2_EC_hint = 4.8399;
64 
65  // parameters for Hcal responses with interaction in Ecal
66  PAR1_HC_hint = 0.2496;
67  PAR2_HC_hint =-1.1673;
68  PAR3_HC_hint = 0.7506;
69 
70  // parameters for Hcal responses without interaction in Ecal (MIP)
71  PAR1_HC_mip = 0.4108;
72  PAR2_HC_mip =-0.3474;
73  PAR3_HC_mip = 0.0293;
74  PAR4_HC_mip = 0.8388;
75  PAR5_HC_mip = 1.0520;
76 
77  if(energycluster>0.5)
78  {
79  double fecal_hint = PAR1_EC_hint/(sqrt(echar) + PAR2_EC_hint);
80  recal = echar*fecal_hint;
81 
82  if(echar <= e_thHC_hint)
83  {
84  double fhcal_hint = PAR1_HC_hint * (log(echar) + PAR2_HC_hint);
85  rhcal=echar*fhcal_hint;
86  }
87  else
88  {
89  rhcal =echar* PAR3_HC_hint;
90  } // hcal response, interact=1
91 
92  // MIP, interact=0
93  }
94  else
95  {
96  recal=energycluster;
97  if(echar <= e_thHC1_mip)
98  {
99  double fhcal_mip = PAR1_HC_mip * (log(echar) + PAR2_HC_mip);
100  rhcal=echar*fhcal_mip;
101  }
102  else if (echar > e_thHC1_mip && echar <= e_thHC2_mip)
103  {
104  double fhcal_mip = PAR3_HC_mip * sqrt(echar) + PAR4_HC_mip;
105  rhcal=echar*fhcal_mip;
106  }
107  else
108  {
109  rhcal=echar*PAR5_HC_mip;
110  } // hcal response, MIP
111 
112  } // interact or not
113  }
114 
115  else if (algo == 2) // cmsim133
116  {
117  // threshold on energy of track
118  e_thEC_hint = 0.1;
119  e_thHC_hint = 1000.;
120  e_thHC1_mip = 15.5;
121  e_thHC2_mip = 58.;
122 
123  // parameters for Ecal responses with interaction in Ecal
124  PAR1_EC_hint = -0.99228E-01;
125  PAR2_EC_hint = -8.2298;
126  PAR3_EC_hint = -0.82878E-03;
127  PAR4_EC_hint = 0.41902;
128 
129  // parameters for Hcal responses with interaction in Ecal
130  PAR1_HC_hint = 0.90693E-01;
131  PAR2_HC_hint = 0.86861E-01;
132  PAR3_HC_hint = 0.089524;
133  PAR4_HC_hint = 0.65067;
134 
135  // parameters for Hcal responses without interaction in Ecal (MIP)
136  PAR1_HC_mip = 0.29538;
137  PAR2_HC_mip = -0.12289;
138  PAR3_HC_mip = 0.14940E-01;
139  PAR4_HC_mip = 0.73503;
140  PAR5_HC_mip = 0.84801;
141 
142 
143  if(energycluster>0.5)
144  {
145 
146  if(echar <= e_thEC_hint)
147  {
148  double fecal_hint = PAR1_EC_hint * (echar + PAR2_EC_hint);
149  recal = echar*fecal_hint;
150  }
151  else
152  {
153  double fecal_hint = PAR3_EC_hint * echar + PAR4_EC_hint;
154  recal = echar*fecal_hint;
155  }
156 
157  if(echar <= e_thHC_hint)
158  {
159  double fhcal_hint = PAR1_HC_hint * (log(echar) + PAR2_HC_hint);
160  rhcal=echar*fhcal_hint;
161  }
162  else
163  {
164  double fhcal_hint = PAR3_HC_hint * (log(echar) + PAR4_HC_hint);
165  rhcal=echar*fhcal_hint;
166  } // hcal response, interact=1
167 
168  // MIP, interact=0
169  }
170  else
171  {
172  recal=energycluster;
173  if(echar <= e_thHC1_mip)
174  {
175  double fhcal_mip = PAR1_HC_mip * (log(echar) + PAR2_HC_mip);
176  rhcal=echar*fhcal_mip;
177  }
178  else if (echar > e_thHC1_mip && echar <= e_thHC2_mip)
179  {
180  double fhcal_mip = PAR3_HC_mip * sqrt(echar) + PAR4_HC_mip;
181  rhcal=echar*fhcal_mip;
182  }
183  else
184  {
185  rhcal=echar*PAR5_HC_mip;
186  } // hcal response, MIP
187 
188  } // interact or not
189  }
190 
191  else if (algo == 3) // OSCAR_3_6_0
192  {
193  // threshold on energy of track
194  e_thEC_hint = 1000.;
195  e_thHC_hint = 0.5;
196  e_thHC1_mip = 12.;
197  e_thHC2_mip = 30.;
198 
199  // parameters for Ecal responses with interaction in Ecal
200  PAR1_EC_hint = -0.54951E-03;
201  PAR2_EC_hint = 0.42609;
202  PAR3_EC_hint = -2.8831;
203  PAR4_EC_hint = 0.33487;
204 
205  // parameters for Hcal responses with interaction in Ecal
206  PAR1_HC_hint = -0.028222;
207  PAR2_HC_hint = 0.43868;
208  PAR3_HC_hint = 0.12144;
209  PAR4_HC_hint = -0.63474;
210 
211  // parameters for Hcal responses without interaction in Ecal (MIP)
212  PAR1_HC_mip = 0.86455E-01;
213  PAR2_HC_mip = -0.36659;
214  PAR3_HC_mip = 0.60879E-01;
215  PAR4_HC_mip = 0.64219;
216  PAR5_HC_mip = 0.96800;
217 
218  if(energycluster>0.5)
219  {
220  if(echar <= e_thEC_hint)
221  {
222  double fecal_hint = PAR1_EC_hint * (echar) + PAR2_EC_hint;
223  recal = echar*fecal_hint;
224  }
225  else
226  {
227  double fecal_hint = PAR3_EC_hint/(-echar) + PAR4_EC_hint;
228  recal = echar*fecal_hint;
229  }
230 
231  if(echar <= e_thHC_hint)
232  {
233  double fhcal_hint = PAR1_HC_hint * echar + PAR2_HC_hint;
234  rhcal=echar*fhcal_hint;
235  }
236  else
237  {
238  double fhcal_hint = PAR3_HC_hint * (log(echar) + PAR4_HC_hint);
239  rhcal=echar*fhcal_hint;
240  } // hcal response, interact=1
241 
242  // MIP, interact=0
243  }
244  else
245  {
246  recal=energycluster;
247  if(echar <= e_thHC1_mip)
248  {
249  double fhcal_mip = PAR1_HC_mip * (echar + PAR2_HC_mip);
250  rhcal=echar*fhcal_mip;
251  }
252  else if (echar > e_thHC1_mip && echar <= e_thHC2_mip)
253  {
254  double fhcal_mip = PAR3_HC_mip * sqrt(echar) + PAR4_HC_mip;
255  rhcal=echar*fhcal_mip;
256  } else
257  {
258  rhcal=echar*PAR5_HC_mip;
259  } // hcal response, MIP
260 
261  } // interact or not
262  }
263 
264  response.push_back(recal);
265  response.push_back(rhcal);
266  return response;
267 }
response
std::pair< double, double > response
Definition: HCALResponse.h:20
ExpressReco_HICollisions_FallBack.algo
tuple algo
Definition: ExpressReco_HICollisions_FallBack.py:2684
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:28
funct::log
Log< T >::type log(const T &t)
Definition: Log.h:22
SingleParticleJetResponse::response
std::vector< double > response(double echar, double energycluster, int algo=0) const
Definition: SingleParticleJetResponse.cc:7
Generated for CMSSW Reference Manual by doxygen 1.8.5