Changeset b03c0d0 in git for examples/Example7.C

Timestamp:

May 26, 2021, 4:46:34 PM (3 years ago)

Author:

Michele Selvaggi <michele.selvaggi@…>

Branches:

master

Children:

c1780a5

Parents:

db316cd

Message:

updated examples to include gen/LHE weights

File:

• examples/Example7.C (modified) (3 diffs)

examples/Example7.C

@@ -30 +30 @@
   TClonesArray *branchElectron = treeReader->UseBranch("ElectronMedium");
   TClonesArray *branchWeight = treeReader->UseBranch("Weight");
+  TClonesArray *branchEvent = treeReader->UseBranch("Event");
 
   // Book histograms
@@ -40 +41 @@
     // Load selected branches with data from specified event
     treeReader->ReadEntry(entry);
-
-    Double_t w =1.0;
-    // read MC event weight
+    
+    // main MC event weight
+    HepMCEvent *event = (HepMCEvent*) branchEvent -> At(0);
+    Double_t w = event->Weight;
+    
+    // read lhe event weights
     if(branchWeight->GetEntries() > 0)
     {
       Weight *weight = (Weight*) branchWeight -> At(0);
-      w = weight->Weight;
+      Double_t lhe_weight = weight->Weight;
+      
+      //cout<<lhe_weight<<endl;
+      // do stuff ...
     }
-    
+
     // If event contains at least 1 jet
     if(branchJet->GetEntries() > 0)
@@ -66 +73 @@
         histJetPT->Fill(jet->PT, w);
     }
-
 
     // If event contains at least 1 jet