ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/LJMet/MultivariateAnalysis/root/event_count.C
Revision: 1.3
Committed: Thu Mar 5 17:30:21 2009 UTC (16 years, 2 months ago) by kukartse
Content type: text/plain
Branch: MAIN
CVS Tags: V00-03-01, ZMorph_BASE_20100408, gak040610_morphing, V00-02-02, gak011410, gak010310, ejterm2010_25nov2009, V00-02-01, V00-02-00, gak112409, CMSSW_22X_branch_base, segala101609, V00-01-15, V00-01-14, V00-01-13, V00-01-12, V00-01-11, V00-01-10, gak031009, gak030509, HEAD
Branch point for: ZMorph-V00-03-01, CMSSW_22X_branch
Changes since 1.2: +9 -3 lines
Log Message: 
*** empty log message ***

File Contents

# User Rev Content
1 kukartse 1.1 {
2     TChain chain("ttljets");
3 kukartse 1.2 //chain . Add("/uscms_data/d1/lpcljm/MVA/Summer08/PAT_Layer1/InclusiveMu5Pt50_Summer08_IDEAL_V9_v1_GEN-SIM-RECO/secondary/InclusiveMu5Pt50_Summer08_IDEAL_V9_v1_GEN-SIM-RECO-secondary-22dec2008.root");
4 kukartse 1.3 chain . Add("/uscms/home/kukarzev/work/disk_ljmet/MVA/Summer08/Jan_23_2009/PAT_Layer1/InclusiveMuPt15_Summer08_IDEAL_V9_v1_GEN-SIM-RECO/secondary/ljmet_tmva_InclusiveMuPt15-12feb2009.root");
5 kukartse 1.2 //chain . Add("/uscms_data/d1/lpcljm/MVA/Summer08/Jan_23_2009/PAT_Layer1/TTJets-madgraph_Fall08_IDEAL_V9_v1_GEN-SIM-RECO/secondary/ljmet_tmva_ttjets-06feb2009.root");
6 kukartse 1.3 //chain . Add("/uscms_data/d1/lpcljm/MVA/Summer08/Jan_23_2009/PAT_Layer1/WJets-madgraph_Fall08_IDEAL_V9_v1_GEN-SIM-RECO/secondary/ljmet_tmva_wjets-06feb2009.root");
7 kukartse 1.2 //chain . Add("/uscms_data/d1/lpcljm/MVA/Summer08/Jan_23_2009/PAT_Layer1/ZJets-madgraph_Fall08_IDEAL_V9_v1_GEN-SIM-RECO/secondary/ljmet_tmva_zjets-06feb2009.root");
8 kukartse 1.1
9     //gROOT->Macro("./cuts.C");
10    
11 kukartse 1.3 // selection
12     #include "cuts.C"
13     TCut nom_cut = semilep_cut && reversible_cut;
14     TCut rev_cut = semilep_cut && !reversible_cut;
15     the_cut = nom_cut;
16     the_cut = semilep_cut && lepton_jet_dr05_cut && !reversible_cut;;
17 kukartse 1.1
18     TH1F h("h","Event weight",1000000,0, 1000);
19     TH1F h_weight("h_weight","Event weight",1000000,0, 1000);
20    
21     //TCut the_cut = "lepton_calo_iso<7.0 && muon_isGood_AllGlobalMuons==1";
22     //TCut the_cut = "(1.0/(1.0+(lepton_calo_iso+lepton_track_iso)/lepton_pt))>0.98 && muon_isGood_AllGlobalMuons==1";
23 kukartse 1.2 //TCut the_cut = "n_tagged_jets_trackCounting_loose>1";
24 kukartse 1.1 //TCut the_cut = "";
25 kukartse 1.3 //TCut the_cut = nominal_cut;
26 kukartse 1.1
27     // weighted count
28     //chain . Draw("event_weight>>h_weight", weight*the_cut);
29     //unweighted count
30     chain . Draw("event_weight>>h", the_cut);
31    
32     //cout << "Weighted events: " << h_weight . Integral() << endl;
33     cout << "Not weighted events: " << h . Integral() << endl;
34     }