Changes between Version 13 and Version 14 of WorkBook/QuickTour

Timestamp:

Feb 25, 2013, 1:08:55 AM (12 years ago)

Author:

Pavel Demin

Comment:

--

WorkBook/QuickTour

@@ -98 +98 @@
 void Example1(const char *inputFile)
 {
+  gSystem->Load("libDelphes");
+
   // Create chain of root trees
   TChain chain("Delphes");
@@ -117 +119 @@
   for(Int_t entry = 0; entry < numberOfEntries; ++entry)
   {
-
     // Load selected branches with data from specified event
     treeReader->ReadEntry(entry);
@@ -124 +125 @@
     if(branchJet->GetEntries() > 0)
     {
-
       // Take first jet
-      TRootJet *jet = (TRootJet*) branchJet->At(0);
+      Jet *jet = (Jet*) branchJet->At(0);
       
       // Plot jet transverse momentum
@@ -135 +135 @@
     }
 
-    TRootElectron *elec1, *elec2;
-    TLorentzVector vec1, vec2;
+    Electron *elec1, *elec2;
 
     // If event contains at least 2 electrons
     if(branchElectron->GetEntries() > 1)
     {
-
       // Take first two electrons
-      elec1 = (TRootElectron *) branchElectron->At(0);
-      elec2 = (TRootElectron *) branchElectron->At(1);
-
-      // Create two 4-vectors for the electrons
-      vec1.SetPtEtaPhiM(elec1->PT, elec1->Eta, elec1->Phi, 0.0);
-      vec2.SetPtEtaPhiM(elec2->PT, elec2->Eta, elec2->Phi, 0.0);
+      elec1 = (Electron *) branchElectron->At(0);
+      elec2 = (Electron *) branchElectron->At(1);
 
       // Plot their invariant mass
-      histMass->Fill((vec1 + vec2).M());
+      histMass->Fill(((elec1->P4()) + (elec2->P4())).M());
     }
   }