Skip to content

Commit

Permalink
apply code format
Browse files Browse the repository at this point in the history
  • Loading branch information
AndrewLevin committed Nov 10, 2019
1 parent 30645fe commit 5546d11
Showing 1 changed file with 29 additions and 32 deletions.
61 changes: 29 additions & 32 deletions GeneratorInterface/RivetInterface/interface/RivetAnalysis.h
Original file line number Diff line number Diff line change
Expand Up @@ -39,7 +39,7 @@ namespace Rivet {
double _jetConeSize, _jetMinPt, _jetMaxEta;
double _fatJetConeSize, _fatJetMinPt, _fatJetMaxEta;

double _phoMinPt, _phoMaxEta ,_phoIsoConeSize, _phoMaxRelIso;
double _phoMinPt, _phoMaxEta, _phoIsoConeSize, _phoMaxRelIso;

Particles _leptons, _photons, _neutrinos;
Jets _jets, _fatjets;
Expand Down Expand Up @@ -67,10 +67,10 @@ namespace Rivet {
_fatJetMinPt(pset.getParameter<double>("fatJetMinPt")),
_fatJetMaxEta(pset.getParameter<double>("fatJetMaxEta")),

_phoMinPt (pset.getParameter<double>("phoMinPt")),
_phoMaxEta (pset.getParameter<double>("phoMaxEta")),
_phoIsoConeSize (pset.getParameter<double>("phoIsoConeSize")),
_phoMaxRelIso (pset.getParameter<double>("phoMaxRelIso")) {}
_phoMinPt(pset.getParameter<double>("phoMinPt")),
_phoMaxEta(pset.getParameter<double>("phoMaxEta")),
_phoIsoConeSize(pset.getParameter<double>("phoIsoConeSize")),
_phoMaxRelIso(pset.getParameter<double>("phoMaxRelIso")) {}

// Initialize Rivet projections
void init() override {
Expand Down Expand Up @@ -167,34 +167,31 @@ namespace Rivet {
}

// Photons
Particles fsparticles = apply<FinalState>(event,"FS").particles();

for ( auto & photon : apply<FinalState>(event, "Photons").particlesByPt() ) {

if (photon.pt() < _phoMinPt)
continue;

if (abs(photon.eta()) > _phoMaxEta)
continue;

double photonptsum = 0;

for (auto &fsparticle : fsparticles) {

if (deltaR(fsparticle, photon) == 0)
continue;

if (deltaR(fsparticle, photon) > _phoIsoConeSize)
continue;

photonptsum += fsparticle.pt();
}
Particles fsparticles = apply<FinalState>(event, "FS").particles();

for (auto& photon : apply<FinalState>(event, "Photons").particlesByPt()) {
if (photon.pt() < _phoMinPt)
continue;

if (abs(photon.eta()) > _phoMaxEta)
continue;

double photonptsum = 0;

for (auto& fsparticle : fsparticles) {
if (deltaR(fsparticle, photon) == 0)
continue;

if (deltaR(fsparticle, photon) > _phoIsoConeSize)
continue;

photonptsum += fsparticle.pt();
}

if (photonptsum / photon.pt() > _phoMaxRelIso)
continue;

if (photonptsum/photon.pt() > _phoMaxRelIso)
continue;

_photons.push_back(photon);

_photons.push_back(photon);
}

_jets = apply<FastJets>(event, "Jets").jetsByPt(jet_cut);
Expand Down

0 comments on commit 5546d11

Please sign in to comment.