Changes in examples/Example3.C [fc4e460:21f3c04] in git

File:

• examples/Example3.C (modified) (7 diffs)

examples/Example3.C

@@ -1 +1 @@
 /*
-
-This macro shows how to access the particle-level reference for reconstructed objects.
-It is also shown how to loop over the jet constituents. 
-
-root -l examples/Example3.C\(\"delphes_output.root\"\)
+root -l examples/Example3.C'("delphes_output.root")'
 */
 
@@ -40 +36 @@
 
   plots->fElectronDeltaPT = result->AddHist1D(
-    "electron_delta_pt", "(p_{T}^{particle} - p_{T}^{electron})/p_{T}^{particle}",
+    "electron delta pt", "(p_{T}^{particle} - p_{T}^{electron})/p_{T}^{particle}",
     "(p_{T}^{particle} - p_{T}^{electron})/p_{T}^{particle}", "number of electrons",
     100, -0.1, 0.1);
 
   plots->fElectronDeltaEta = result->AddHist1D(
-    "electron_delta_eta", "(#eta^{particle} - #eta^{electron})/#eta^{particle}",
+    "electron delta eta", "(#eta^{particle} - #eta^{electron})/#eta^{particle}",
     "(#eta^{particle} - #eta^{electron})/#eta^{particle}", "number of electrons",
     100, -0.1, 0.1);
 
   plots->fPhotonDeltaPT = result->AddHist1D(
-    "photon_delta_pt", "(p_{T}^{particle} - p_{T}^{photon})/p_{T}^{particle}",
+    "photon delta pt", "(p_{T}^{particle} - p_{T}^{photon})/p_{T}^{particle}",
     "(p_{T}^{particle} - p_{T}^{photon})/p_{T}^{particle}", "number of photons",
     100, -0.1, 0.1);
 
   plots->fPhotonDeltaEta = result->AddHist1D(
-    "photon_delta_eta", "(#eta^{particle} - #eta^{photon})/#eta^{particle}",
+    "photon delta eta", "(#eta^{particle} - #eta^{photon})/#eta^{particle}",
     "(#eta^{particle} - #eta^{photon})/#eta^{particle}", "number of photons",
     100, -0.1, 0.1);
 
   plots->fPhotonDeltaE = result->AddHist1D(
-    "photon_delta_energy", "(E^{particle} - E^{photon})/E^{particle}",
+    "photon delta energy", "(E^{particle} - E^{photon})/E^{particle}",
     "(E^{particle} - E^{photon})/E^{particle}", "number of photons",
     100, -0.1, 0.1);
 
   plots->fMuonDeltaPT = result->AddHist1D(
-    "muon_delta_pt", "(p_{T}^{particle} - p_{T}^{muon})/p_{T}^{particle}",
+    "muon delta pt", "(p_{T}^{particle} - p_{T}^{muon})/p_{T}^{particle}",
     "(p_{T}^{particle} - p_{T}^{muon})/p_{T}^{particle}", "number of muons",
     100, -0.1, 0.1);
 
   plots->fMuonDeltaEta = result->AddHist1D(
-    "muon_delta_eta", "(#eta^{particle} - #eta^{muon})/#eta^{particle}",
+    "muon delta eta", "(#eta^{particle} - #eta^{muon})/#eta^{particle}",
     "(#eta^{particle} - #eta^{muon})/#eta^{particle}", "number of muons",
     100, -0.1, 0.1);
 
   plots->fTrackDeltaPT = result->AddHist1D(
-    "track_delta_pt", "(p_{T}^{particle} - p_{T}^{track})/p_{T}^{particle}",
+    "track delta pt", "(p_{T}^{particle} - p_{T}^{track})/p_{T}^{particle}",
     "(p_{T}^{particle} - p_{T}^{track})/p_{T}^{particle}", "number of tracks",
     100, -0.1, 0.1);
 
   plots->fTrackDeltaEta = result->AddHist1D(
-    "track_delta_eta", "(#eta^{particle} - #eta^{track})/#eta^{particle}",
+    "track delta eta", "(#eta^{particle} - #eta^{track})/#eta^{particle}",
     "(#eta^{particle} - #eta^{track})/#eta^{particle}", "number of tracks",
     100, -0.1, 0.1);
 
   plots->fJetDeltaPT = result->AddHist1D(
-    "jet_delta_pt", "(p_{T}^{jet} - p_{T}^{constituents})/p_{T}^{jet}",
+    "jet delta pt", "(p_{T}^{jet} - p_{T}^{constituents})/p_{T}^{jet}",
     "(p_{T}^{jet} - p_{T}^{constituents})/p_{T}^{jet}", "number of jets",
     100, -1.0e-1, 1.0e-1);
@@ -183 +179 @@
     }
 
-  //  cout<<"--  New event -- "<<endl;
+    cout<<"--  New event -- "<<endl;
 
     // Loop over all jets in event
@@ -192 +188 @@
       momentum.SetPxPyPzE(0.0, 0.0, 0.0, 0.0);
 
-      //cout<<"Looping over jet constituents. Jet pt: "<<jet->PT<<", eta: "<<jet->Eta<<", phi: "<<jet->Phi<<endl;      
+      cout<<"Looping over jet constituents. Jet pt: "<<jet->PT<<", eta: "<<jet->Eta<<", phi: "<<jet->Phi<<endl;      
 
       // Loop over all jet's constituents
@@ -205 +201 @@
         {
           particle = (GenParticle*) object;
-          //cout << "    GenPart pt: " << particle->PT << ", eta: " << particle->Eta << ", phi: " << particle->Phi << endl;
+          cout << "    GenPart pt: " << particle->PT << ", eta: " << particle->Eta << ", phi: " << particle->Phi << endl;
           momentum += particle->P4();
         }
@@ -211 +207 @@
         {
           track = (Track*) object;
-          //cout << "    Track pt: " << track->PT << ", eta: " << track->Eta << ", phi: " << track->Phi << endl;
+          cout << "    Track pt: " << track->PT << ", eta: " << track->Eta << ", phi: " << track->Phi << endl;
           momentum += track->P4();
         }
@@ -217 +213 @@
         {
           tower = (Tower*) object;
-          //cout << "    Tower pt: " << tower->ET << ", eta: " << tower->Eta << ", phi: " << tower->Phi << endl;
+          cout << "    Tower pt: " << tower->ET << ", eta: " << tower->Eta << ", phi: " << tower->Phi << endl;
           momentum += tower->P4();
+        }
+        else if(object->IsA() == Muon::Class())
+        {
+          muon = (Muon*) object;
+          cout << "    Muon pt: " << muon->PT << ", eta: " << muon->Eta << ", phi: " << muon->Phi << endl;
+          momentum += muon->P4();
         }
       }