Changes in / [c1780a5:b4786d3] in git

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
@@ -41 +40 @@
     // Load selected branches with data from specified event
     treeReader->ReadEntry(entry);
-    
-    // main MC event weight
-    HepMCEvent *event = (HepMCEvent*) branchEvent -> At(0);
-    Double_t w = event->Weight;
-    
-    // read lhe event weights
+
+    Double_t w =1.0;
+    // read MC event weight
     if(branchWeight->GetEntries() > 0)
     {
       Weight *weight = (Weight*) branchWeight -> At(0);
-      Double_t lhe_weight = weight->Weight;
-      
-      //cout<<lhe_weight<<endl;
-      // do stuff ...
+      w = weight->Weight;
     }
-
+    
     // If event contains at least 1 jet
     if(branchJet->GetEntries() > 0)
@@ -73 +66 @@
         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")
-branchEvent    = treeReader.UseBranch("Event")
+branchWeight = treeReader.UseBranch("Weight")
 
 # Book histograms
@@ -47 +46 @@
   treeReader.ReadEntry(entry)
 
-  ## main MC event weight
-  w =  branchEvent[0].Weight
+  w = 1.0  
+  ## read MC event weight
+  if branchWeight.GetEntries() > 0:  
+    weight = branchWeight.At(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:
@@ -81 +76 @@
     eta = abs(electron.Eta)
 
-    ## looseCut = 0.3, mediumCut = 0.2, tightCut = 0.1
     IsoCut = 0.2
     IsoOk = electron.IsolationVar < IsoCut