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

Timestamp:

Mar 9, 2009, 9:47:21 PM (16 years ago)

Author:

severine ovyn

Message:

LHCO advanced

File:

• trunk/src/LHCOConverter.cc (modified) (4 diffs)

trunk/src/LHCOConverter.cc

@@ -66 +66 @@
         }
         else {
-                outputfilename_.replace(r_find,5,".lhco");
+                outputfilename_.replace(r_find,5,"_events.lhco");
                 ofstream outfile( outputfilename_.c_str());
                 outfile.close();
@@ -138 +138 @@
 
     outfile << setw(3) << line++  << setw(4) << " " << setw(7) << event << setw(7) << triginfo << endl;
-    outfile.close();
 
     // 0 photon data
-    //outfile << BranchReader(branchPhoton,line,lhcoPhotonID);
     BranchReader(branchPhoton,line,lhcoPhotonID);
 
@@ -161 +159 @@
   } // event loop
 
+  outfile.close();
   delete triggerTree;
   delete analysisTree;
@@ -180 +179 @@
         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; }
         else if (lhcoID == lhcoTauJetID) { TRootTauJet taujet(*((TRootTauJet*) branch->At(i))); ntrk = taujet.Charge; }
-        float btag =0;
-        double ratioE = 0;
+        else if (lhcoID == lhcoJetID) { TRootJet jet(*((TRootJet*) branch->At(i))); ntrk = jet.NTracks; btag = jet.Btag;}
         if(lhcoID != lhcoETmisID) {
         outfile << fixed << setprecision(3)