|
[Rivet-svn] r3396 - trunk/src/Analysesblackhole at projects.hepforge.org blackhole at projects.hepforge.orgWed Sep 28 15:08:08 BST 2011
Author: hoeth Date: Wed Sep 28 15:08:08 2011 New Revision: 3396 Log: cosmetics Modified: trunk/src/Analyses/ATLAS_2011_S9120807.cc Modified: trunk/src/Analyses/ATLAS_2011_S9120807.cc ============================================================================== --- trunk/src/Analyses/ATLAS_2011_S9120807.cc Wed Sep 28 14:58:59 2011 (r3395) +++ trunk/src/Analyses/ATLAS_2011_S9120807.cc Wed Sep 28 15:08:08 2011 (r3396) @@ -67,7 +67,7 @@ double eta = fabs(eta_w); int v_iter=0; - for (v_iter=0; v_iter < (int)_eta_bins_areaoffset.size()-1; v_iter++){ + for (v_iter=0; v_iter < (int)_eta_bins_areaoffset.size()-1; v_iter++) { if (inRange(eta, _eta_bins_areaoffset[v_iter], _eta_bins_areaoffset[v_iter+1])) break; } @@ -84,7 +84,7 @@ /// require at least 2 photons in final state /// ParticleVector photons = applyProjection<IdentifiedFinalState>(event, "Photon").particlesByPt(); - if (photons.size() < 2){ + if (photons.size() < 2) { vetoEvent; } @@ -104,19 +104,18 @@ /// get the cluster sequence double area = clust_seq_area->area(jet); - if(area > 10e-4 && fabs(eta)<_eta_bins_areaoffset[_eta_bins_areaoffset.size()-1]){ + if(area > 10e-4 && fabs(eta) < _eta_bins_areaoffset[_eta_bins_areaoffset.size()-1]) { ptDensities.at(getEtaBin(fabs(eta))).push_back(pt/area); } } - for(int b=0; b<(int)_eta_bins_areaoffset.size()-1;b++){ + for(int b=0; b<(int)_eta_bins_areaoffset.size()-1; b++) { double median = 0.0; - if(ptDensities[b].size() > 0) - { + if(ptDensities[b].size() > 0) { std::sort(ptDensities[b].begin(), ptDensities[b].end()); int nDens = ptDensities[b].size(); if( nDens%2 == 0 ) - median = (ptDensities[b][nDens/2]+ptDensities[b][(nDens-2)/2])/2; + median = (ptDensities[b][nDens/2] + ptDensities[b][(nDens-2)/2]) / 2; else median = ptDensities[b][(nDens-1)/2]; } @@ -134,7 +133,7 @@ /// remove photons in crack /// double eta_P = photon.momentum().eta(); - if(fabs(eta_P)>=1.37 && fabs(eta_P)<1.52) continue; + if (fabs(eta_P)>=1.37 && fabs(eta_P)<1.52) continue; double phi_P = photon.momentum().phi(); @@ -162,7 +161,7 @@ /// shouldn't need to subtract photon /// note: using expected cut at hadron/particle level, not cut at reco level - if(mom_in_EtCone.Et()-correction > 4.0*GeV){ + if (mom_in_EtCone.Et()-correction > 4.0*GeV) { continue; } @@ -191,10 +190,10 @@ vetoEvent; } - FourMomentum yy=y1+y2; + FourMomentum yy = y1+y2; double Myy = yy.mass()/GeV; double pTyy = yy.pT()/GeV; - double dPhiyy = mapAngle0ToPi(y1.phi()-y2.phi()); + double dPhiyy = deltaPhi(y1.phi(), y2.phi()); _h_M->fill(Myy, weight); _h_pT->fill(pTyy, weight);
More information about the Rivet-svn mailing list |