52 active_area = passive_area = 0.0;
67 active_area = passive_area = 0.0;
126 int _n_pass_max, Esplit_merge_scale _split_merge_scale,
129 vector<Cmomentum> all_particles;
135 SM_var2_hardest_cut_off = pt_soft_min*pt_soft_min;
142 int n_hard = _particles.size();
143 all_particles = _particles;
147 double eta_g,phi_g,pt_g;
149 for (i=0;i<grid_size;i++){
150 for (j=0;j<grid_size;j++){
151 eta_g = grid_eta_max*(-1.0+2.0*(i+0.5+grid_shift*(-1.0+2.0*(rand()/(RAND_MAX+1.0))))/grid_size);
152 phi_g = M_PI *(-1.0+2.0*(j+0.5+grid_shift*(-1.0+2.0*(rand()/(RAND_MAX+1.0))))/grid_size);
153 pt_g = pt_soft*(1.0+pt_shift*(-1.0+2.0*(rand()/(RAND_MAX+1.0))));
154 all_particles.push_back(
Cmomentum(pt_g*cos(phi_g),pt_g*sin(phi_g),pt_g*sinh(eta_g),pt_g*cosh(eta_g)));
161 int n_jets = compute_jets(all_particles, _radius, _f, _n_pass_max, 0.0, _split_merge_scale);
168 double area_factor = (2.0*grid_eta_max/grid_size)*(
twopi/grid_size);
170 for (i=0;i<(int) jets.size();i++){
171 jet_areas.push_back(jets[i]);
173 while ((j<jets[i].n) && (jets[i].contents[j]<n_hard)) j++;
174 jet_areas[i].active_area = (jets[i].n-j)*area_factor;
180 recompute_jets(_f, pt_soft_min);
183 for (i=0;i<(int) jets.size();i++){
185 while ((j<jets[i].n) && (jets[i].contents[j]<n_hard)) j++;
186 jet_areas[i].passive_area = (jets[i].n-j)*area_factor;
218 int _n_pass_max, Esplit_merge_scale _split_merge_scale){
220 vector<Cmomentum> all_particles;
226 stable_cone_soft_pt2_cutoff = pt_soft_min*pt_soft_min;
232 int n_hard = _particles.size();
233 all_particles = _particles;
237 double eta_g,phi_g,pt_g;
239 for (i=0;i<grid_size;i++){
240 for (j=0;j<grid_size;j++){
241 eta_g = grid_eta_max*(-1.0+2.0*(i+0.5+grid_shift*(-1.0+2.0*(rand()/(RAND_MAX+1.0))))/grid_size);
242 phi_g = M_PI *(-1.0+2.0*(j+0.5+grid_shift*(-1.0+2.0*(rand()/(RAND_MAX+1.0))))/grid_size);
243 pt_g = pt_soft*(1.0+pt_shift*(-1.0+2.0*(rand()/(RAND_MAX+1.0))));
244 all_particles.push_back(
Cmomentum(pt_g*cos(phi_g),pt_g*sin(phi_g),pt_g*sinh(eta_g),pt_g*cosh(eta_g)));
251 int n_jets = compute_jets(all_particles, _radius, _f, _n_pass_max, pt_soft_min, _split_merge_scale);
258 double area_factor = (2.0*grid_eta_max/grid_size)*(
twopi/grid_size);
259 for (i=0;i<(int) jets.size();i++){
261 while ((j<jets[i].n) && (jets[i].contents[j]<n_hard)) j++;
262 jet_areas[i].passive_area = (jets[i].n-j)*area_factor;
292 int _n_pass_max, Esplit_merge_scale _split_merge_scale,
295 vector<Cmomentum> all_particles;
301 SM_var2_hardest_cut_off = pt_soft_min*pt_soft_min;
308 int n_hard = _particles.size();
309 all_particles = _particles;
313 double eta_g,phi_g,pt_g;
315 for (i=0;i<grid_size;i++){
316 for (j=0;j<grid_size;j++){
317 eta_g = grid_eta_max*(-1.0+2.0*(i+0.5+grid_shift*(-1.0+2.0*(rand()/(RAND_MAX+1.0))))/grid_size);
318 phi_g = M_PI *(-1.0+2.0*(j+0.5+grid_shift*(-1.0+2.0*(rand()/(RAND_MAX+1.0))))/grid_size);
319 pt_g = pt_soft*(1.0+pt_shift*(-1.0+2.0*(rand()/(RAND_MAX+1.0))));
320 all_particles.push_back(
Cmomentum(pt_g*cos(phi_g),pt_g*sin(phi_g),pt_g*sinh(eta_g),pt_g*cosh(eta_g)));
327 int n_jets = compute_jets(all_particles, _radius, _f, _n_pass_max, 0.0, _split_merge_scale);
334 double area_factor = (2.0*grid_eta_max/grid_size)*(
twopi/grid_size);
336 for (i=0;i<(int) jets.size();i++){
337 jet_areas.push_back(jets[i]);
339 while ((j<jets[i].n) && (jets[i].contents[j]<n_hard)) j++;
340 jet_areas[i].active_area = (jets[i].n-j)*area_factor;