Changeset fc4e460 in git for examples/Example3.C

Timestamp:

Sep 9, 2014, 1:26:52 PM (10 years ago)

Author:

Michele <michele.selvaggi@…>

Branches:

ImprovedOutputFile, Timing, dual_readout, llp, master

Children:

fcc959c

Parents:

492aef9

Message:

removed muon check in jet constituent loop. Added comments to example macros

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\"\)
 */
@@ -36 +40 @@
 
   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);
@@ -179 +183 @@
     }
 
-    cout<<"--  New event -- "<<endl;
+  //  cout<<"--  New event -- "<<endl;
 
     // Loop over all jets in event
@@ -188 +192 @@
       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
@@ -201 +205 @@
         {
           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();
         }
@@ -207 +211 @@
         {
           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();
         }
@@ -213 +217 @@
         {
           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();
         }
       }