Changeset 310 in svn for trunk/src/LHCOConverter.cc

Timestamp:

Mar 10, 2009, 11:52:16 AM (16 years ago)

Author:

severine ovyn

Message:

add Ehad/Eem

File:

• trunk/src/LHCOConverter.cc (modified) (1 diff)

trunk/src/LHCOConverter.cc

@@ -174 +174 @@
         particle = (TRootParticle*) branch->At(i);
         double temp =  sqrt(particle->PT*particle->PT + particle->Pz * particle->Pz);
-        double jmass = (particle->E - temp) * (particle->E + temp); // prevents some numerical sensitivity
         TLorentzVector pmu;
         pmu.SetPtEtaPhiE(particle->PT, particle->Eta, particle->Phi, particle->E);
-        jmass = pmu.M();
+        float jmass = pmu.M();
         float ntrk = 0.0;
         float btag =0;
         double ratioE = 0;
 
-        if(lhcoID == lhcoElectronID) { TRootElectron elec(*((TRootElectron*) branch->At(i))); ntrk = elec.Charge; }
-        else if (lhcoID == lhcoMuonID) { TRootMuon muon(*((TRootMuon*) branch->At(i))); ntrk = muon.Charge; }
+        if(lhcoID == lhcoPhotonID){TRootPhoton gam(*((TRootPhoton*) branch->At(i))); jmass=0; ratioE = gam.EHoverEE;}
+        if(lhcoID == lhcoElectronID) { TRootElectron elec(*((TRootElectron*) branch->At(i))); ntrk = elec.Charge; ratioE = elec.EHoverEE;}
+        else if (lhcoID == lhcoMuonID) { TRootMuon muon(*((TRootMuon*) branch->At(i))); ntrk = muon.Charge;  ratioE = muon.EHoverEE;}
         else if (lhcoID == lhcoTauJetID) { TRootTauJet taujet(*((TRootTauJet*) branch->At(i))); ntrk = taujet.Charge; }
-        else if (lhcoID == lhcoJetID) { TRootJet jet(*((TRootJet*) branch->At(i))); ntrk = jet.NTracks; btag = jet.Btag;}
+        else if (lhcoID == lhcoJetID) { TRootJet jet(*((TRootJet*) branch->At(i))); ntrk = jet.NTracks; btag = jet.Btag; ratioE = jet.EHoverEE;}
         if(lhcoID != lhcoETmisID) {
         outfile << fixed << setprecision(3)