|
[Rivet-svn] r3746 - trunk/src/Analysesblackhole at projects.hepforge.org blackhole at projects.hepforge.orgFri Jun 8 15:04:50 BST 2012
Author: hoeth Date: Fri Jun 8 15:04:50 2012 New Revision: 3746 Log: whitespace cleanup to follow Rivet conventions Modified: trunk/src/Analyses/ATLAS_2012_I1082009.cc Modified: trunk/src/Analyses/ATLAS_2012_I1082009.cc ============================================================================== --- trunk/src/Analyses/ATLAS_2012_I1082009.cc Fri Jun 8 15:04:07 2012 (r3745) +++ trunk/src/Analyses/ATLAS_2012_I1082009.cc Fri Jun 8 15:04:50 2012 (r3746) @@ -22,8 +22,8 @@ /// Constructor ATLAS_2012_I1082009() : Analysis("ATLAS_2012_I1082009"), - _weight25_30(0.),_weight30_40(0.),_weight40_50(0.), - _weight50_60(0.),_weight60_70(0.),_weight25_70(0.) + _weight25_30(0.),_weight30_40(0.),_weight40_50(0.), + _weight50_60(0.),_weight60_70(0.),_weight25_70(0.) { } //@} @@ -62,54 +62,54 @@ // get the jets Jets jets; foreach (const Jet& jet, applyProjection<FastJets>(event, "jets").jetsByPt(25.0*GeV)) { - if ( fabs(jet.eta()) < 2.5 ) jets.push_back(jet); + if ( fabs(jet.eta()) < 2.5 ) jets.push_back(jet); } // get the D* mesons const UnstableFinalState& ufs = applyProjection<UnstableFinalState>(event, "UFS"); ParticleVector Dstar; foreach (const Particle& p, ufs.particles()) { const int id = abs(p.pdgId()); - if(id==413) Dstar.push_back(p); + if(id==413) Dstar.push_back(p); } // loop over the jobs foreach (const Jet& jet, jets ) { - double perp = jet.momentum().perp(); - bool found = false; - double z(0.); - if(perp<25.||perp>70.) continue; - foreach(const Particle & p, Dstar) { - if(p.momentum().perp()<7.5) continue; - if(deltaR(p, jet.momentum())<0.6) { - Vector3 axis = jet.momentum().vector3().unit(); - z = axis.dot(p.momentum().vector3())/jet.momentum().E(); - if(z<0.3) continue; - found = true; - break; - } - } - _weight25_70 += weight; - if(found) _h_pt25_70->fill(z,weight); - if(perp>=25.&&perp<30.) { - _weight25_30 += weight; - if(found) _h_pt25_30->fill(z,weight); - } - else if(perp>=30.&&perp<40.) { - _weight30_40 += weight; - if(found) _h_pt30_40->fill(z,weight); - } - else if(perp>=40.&&perp<50.) { - _weight40_50 += weight; - if(found) _h_pt40_50->fill(z,weight); - } - else if(perp>=50.&&perp<60.) { - _weight50_60 += weight; - if(found) _h_pt50_60->fill(z,weight); - } - else if(perp>=60.&&perp<70.) { - _weight60_70 += weight; - if(found) _h_pt60_70->fill(z,weight); - } + double perp = jet.momentum().perp(); + bool found = false; + double z(0.); + if(perp<25.||perp>70.) continue; + foreach(const Particle & p, Dstar) { + if(p.momentum().perp()<7.5) continue; + if(deltaR(p, jet.momentum())<0.6) { + Vector3 axis = jet.momentum().vector3().unit(); + z = axis.dot(p.momentum().vector3())/jet.momentum().E(); + if(z<0.3) continue; + found = true; + break; + } + } + _weight25_70 += weight; + if(found) _h_pt25_70->fill(z,weight); + if(perp>=25.&&perp<30.) { + _weight25_30 += weight; + if(found) _h_pt25_30->fill(z,weight); + } + else if(perp>=30.&&perp<40.) { + _weight30_40 += weight; + if(found) _h_pt30_40->fill(z,weight); + } + else if(perp>=40.&&perp<50.) { + _weight40_50 += weight; + if(found) _h_pt40_50->fill(z,weight); + } + else if(perp>=50.&&perp<60.) { + _weight50_60 += weight; + if(found) _h_pt50_60->fill(z,weight); + } + else if(perp>=60.&&perp<70.) { + _weight60_70 += weight; + if(found) _h_pt60_70->fill(z,weight); + } } } @@ -144,8 +144,6 @@ }; - - // The hook for the plugin system DECLARE_RIVET_PLUGIN(ATLAS_2012_I1082009);
More information about the Rivet-svn mailing list |