Changeset 307 in svn for trunk/src/LHCOConverter.cc
- Timestamp:
- Mar 9, 2009, 9:47:21 PM (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/src/LHCOConverter.cc
r288 r307 66 66 } 67 67 else { 68 outputfilename_.replace(r_find,5," .lhco");68 outputfilename_.replace(r_find,5,"_events.lhco"); 69 69 ofstream outfile( outputfilename_.c_str()); 70 70 outfile.close(); … … 138 138 139 139 outfile << setw(3) << line++ << setw(4) << " " << setw(7) << event << setw(7) << triginfo << endl; 140 outfile.close();141 140 142 141 // 0 photon data 143 //outfile << BranchReader(branchPhoton,line,lhcoPhotonID);144 142 BranchReader(branchPhoton,line,lhcoPhotonID); 145 143 … … 161 159 } // event loop 162 160 161 outfile.close(); 163 162 delete triggerTree; 164 163 delete analysisTree; … … 180 179 jmass = pmu.M(); 181 180 float ntrk = 0.0; 181 float btag =0; 182 double ratioE = 0; 183 182 184 if(lhcoID == lhcoElectronID) { TRootElectron elec(*((TRootElectron*) branch->At(i))); ntrk = elec.Charge; } 183 185 else if (lhcoID == lhcoMuonID) { TRootMuon muon(*((TRootMuon*) branch->At(i))); ntrk = muon.Charge; } 184 186 else if (lhcoID == lhcoTauJetID) { TRootTauJet taujet(*((TRootTauJet*) branch->At(i))); ntrk = taujet.Charge; } 185 float btag =0; 186 double ratioE = 0; 187 else if (lhcoID == lhcoJetID) { TRootJet jet(*((TRootJet*) branch->At(i))); ntrk = jet.NTracks; btag = jet.Btag;} 187 188 if(lhcoID != lhcoETmisID) { 188 189 outfile << fixed << setprecision(3)
Note:
See TracChangeset
for help on using the changeset viewer.