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

Timestamp:

Mar 11, 2009, 12:53:50 PM (16 years ago)

Author:

Xavier Rouby

Message:

Isolation updated. ptiso implemented. etrat prepared but not finished

File:

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

trunk/src/LHCOConverter.cc

@@ -205 +205 @@
           case lhcoMuonID: {
            cout << "ERROR: LHCOConverter: for muons use LHCOConverter::BranchReaderMuon!\n";
-                TRootMuon muon(*((TRootMuon*) branch->At(i))); 
-                ntrk = muon.Charge;  ratioE = muon.EHoverEE; 
           } break;
           case lhcoETmisID: {
            cout << "ERROR: LHCOConverter: for ETmis use LHCOConverter::BranchReaderETmis!\n";
-                TRootMuon muon(*((TRootMuon*) branch->At(i))); 
-                ntrk = muon.Charge;  ratioE = muon.EHoverEE; 
           } break;
 
@@ -284 +280 @@
         } // loop on jets
 
+        int ratio = 0; char ratio_string[3];
+        if(muon->EtRatio != UNDEFINED) ratio =  static_cast<int>(muon->EtRatio);
+        if(ratio<0) sprintf(ratio_string,"00");
+        else if(ratio<10) sprintf(ratio_string,"0%d",ratio);
+        else if(ratio<100)sprintf(ratio_string,"%d",ratio);
+        else { cout << "Error: muon EtRatio > 99" << endl; sprintf(ratio_string,"99"); }
         // output
         outfile << fixed << setprecision(3)
@@ -293 +295 @@
                 << setw(7) << jmass <<" "       // invariant mass
                 << setw(7) << (float)muon->Charge  <<" "        // -1 for mu-, +1 for mu+
-                << setw(7) << (float)closest_jet_ID << " "
-                << setw(7) << muon->EHoverEE <<" "      // E_had/E_em
+                << setw(7) << (float)closest_jet_ID << " ";
+        outfile << fixed << setprecision(0)
+                << setw(4) << floor(muon->IsolPt) << "." << ratio_string<<" ";  // E_had/E_em
+        outfile << fixed << setprecision(3)
                 << setw(7) << 0.0 << setw(7) << 0.0  // dummy/dummy
                 << endl;