21 double eta = fabs(Jet.
p4().Eta());
32 if (EMF-HadF!=0.) EMvsHadF = (EMF+HadF)/(EMF-HadF);
43 double n90 = Jet.
n90();
53 double MaxEnergyTower = 0., SumTowPt=0., SumTowPtR=0.;
54 for(std::vector<CaloTowerPtr>::const_iterator tow = jetTowers.begin(),
55 towend = jetTowers.end(); tow != towend; ++tow){
57 SumTowPt += (*tow)->et();
58 SumTowPtR += (*tow)->et()*
deltaR( Jet.
p4().Eta(), Jet.
p4().Phi(),
59 (*tow)->eta(), (*tow)->phi() );
60 if ( (*tow)->et() > MaxEnergyTower )
61 MaxEnergyTower = (*tow)->et();
65 double EtTowerMaxOverEtJet = 0.;
66 if (Jet.
p4().Et()!=0.) EtTowerMaxOverEtJet = MaxEnergyTower /Jet.
p4().Et();
72 if (SumTowPt!=0.) RWidth = SumTowPtR /SumTowPt;
77 double PtJetoverArea = 0.;
83 double PtToweroverArea = 0.;
double PtTowerOverArea_max
Jets made from CaloTowers.
virtual std::vector< CaloTowerPtr > getCaloConstituents() const
get all constituents
1: Failed selection (without additional info)
bool filter(edm::Event &evt, const edm::EventSetup &es) override
process one event
const LorentzVector & p4() const final
four-momentum Lorentz vector
double HighestTowerOverJetmax
double HighestTowerOverJetmin
double PtTowerOverArea_min
ParticleStatus
Definition of particle status after selection.
float emEnergyFraction() const