Changeset b03c0d0 in git for examples

Timestamp:

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

Author:

Michele Selvaggi <michele.selvaggi@…>

Branches:

master

Children:

c1780a5

Parents:

db316cd

Message:

updated examples to include gen/LHE weights

Location:

examples

Files:

• Example7.C (modified) (3 diffs)
• Example7.py (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

examples/Example7.py

@@ -33 +33 @@
 
 # Get pointers to branches used in this analysis
-branchJet = treeReader.UseBranch("JetPUPPITight")
+branchJet      = treeReader.UseBranch("JetPUPPITight")
 branchElectron = treeReader.UseBranch("ElectronMedium")
-branchWeight = treeReader.UseBranch("Weight")
+branchWeight   = treeReader.UseBranch("Weight")
+branchEvent    = treeReader.UseBranch("Event")
 
 # Book histograms
@@ -46 +47 @@
   treeReader.ReadEntry(entry)
 
-  w = 1.0  
-  ## read MC event weight
-  if branchWeight.GetEntries() > 0:  
-    weight = branchWeight.At(0).Weight
+  ## main MC event weight
+  w =  branchEvent[0].Weight
 
+  ## read lhe event weight
+  for weight in branchWeight:  
+    lhe_weight = weight.Weight 
+    ## do stuff ... 
+    ## print lhe_weight
+    
   # If event contains at least 1 jet
   if branchJet.GetEntries() > 0:
@@ -76 +81 @@
     eta = abs(electron.Eta)
 
+    ## looseCut = 0.3, mediumCut = 0.2, tightCut = 0.1
     IsoCut = 0.2
     IsoOk = electron.IsolationVar < IsoCut