ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/tschum/FWlite_Analysis/start.C
Revision: 1.7
Committed: Wed Dec 9 13:02:16 2009 UTC (15 years, 5 months ago) by tschum
Content type: text/plain
Branch: MAIN
Changes since 1.6: +27 -17 lines
Log Message: 
add EventInfo to friendTree

File Contents

# User Rev Content
1 tschum 1.1 {
2     gROOT->ProcessLine(".L PlotTool.cc+");
3     PlotTool p;
4 tschum 1.2
5     p.addTrackJets = 0;
6 tschum 1.7 p.addEventInfo = 1;
7 tschum 1.4 // p.init("","dcap://dcache-ses-cms.desy.de//pnfs/desy.de/cms/tier2/store/data/BeamCommissioning09/MinimumBias/RECO/rereco_FIRSTCOLL_v1/0083/","Events","FIRSTCOLL");
8 tschum 1.5 // p.init("goldenRun.txt","dcap://dcache-ses-cms.desy.de//pnfs/desy.de/cms/tier2","Events","FIRSTCOLL");
9 tschum 1.6 // p.init("","/scratch/hh/current/cms/user/schum/run123592/","Events","run123592");
10     p.init("run123596_lumi69-139.txt","","Events","run123596_lumi69-130");
11 tschum 1.1
12     vector<string> vars;
13 tschum 1.6 vector<string> cuts;
14    
15 tschum 1.7 string saveDir="/afs/desy.de/user/s/schum/cms/run_123596_Plots/";
16 tschum 1.6 int numEv=-1;
17 tschum 1.1
18 tschum 1.2 vars.push_back("EventAuxiliary.run()");
19     vars.push_back("EventAuxiliary.isRealData()");
20     vars.push_back("EventAuxiliary.bunchCrossing()");
21     vars.push_back("EventAuxiliary.orbitNumber()");
22     vars.push_back("EventAuxiliary.storeNumber()");
23     vars.push_back("EventAuxiliary.event()");
24     vars.push_back("EventAuxiliary.experimentType()");
25     vars.push_back("EventAuxiliary.luminosityBlock()");
26 tschum 1.7 vars.push_back("run");
27     vars.push_back("event");
28     vars.push_back("lumi");
29 tschum 1.1
30 tschum 1.6 cuts.push_back("");
31 tschum 1.7 p.globalCuts="lumi>68";
32 tschum 1.3
33 tschum 1.5 // p.samePad_trees = false;
34     // p.sameCanv_vars = true;
35    
36    
37 tschum 1.6 p.loop(vars,cuts,numEv);
38     p.saveCanvases("EventAuxiliary",saveDir);
39 tschum 1.5 p.clearCanvases();
40    
41 tschum 1.7 p.setVariables("offlinePrimaryVertices");
42 tschum 1.6 p.loop(p.autoVars,cuts,numEv);
43 tschum 1.7 p.saveCanvases("offlinePrimaryVertices",saveDir);
44 tschum 1.5 p.clearCanvases();
45    
46 tschum 1.7 p.setVariables("offlineBeamSpot");
47 tschum 1.6 p.loop(p.autoVars,cuts,numEv);
48 tschum 1.7 p.saveCanvases("offlineBeamSpot",saveDir);
49 tschum 1.5 p.clearCanvases();
50    
51     p.setVariables("recoCaloMETs_metNoHFHO");
52 tschum 1.6 p.loop(p.autoVars,cuts,numEv);
53     p.saveCanvases("metNoHFHO",saveDir);
54 tschum 1.5 p.clearCanvases();
55    
56 tschum 1.7 cuts.clear();
57     cuts.push_back("");
58     cuts.push_back("abs(ak5CaloJets.eta())<2.5");
59     // cuts.push_back(cuts.back()+string("&&ak5CaloJets.n90()>1"));
60     cuts.push_back(cuts.back()+string("&&ak5CaloJets.emEnergyFraction()>0.02&&ak5CaloJets.emEnergyFraction()<0.98"));
61     cuts.push_back(cuts.back()+string("&&(( - ak5CaloJets.vx() * ak5CaloJets.py() + ak5CaloJets.vy() * ak5CaloJets.px() ) / ak5CaloJets.pt())<2&&(ak5CaloJets.vz() - (ak5CaloJets.vx()*ak5CaloJets.px()+ak5CaloJets.vy()*ak5CaloJets.py())/ak5CaloJets.pt() * (ak5CaloJets.pz()/ak5CaloJets.pt()))<20"));
62    
63 tschum 1.5 p.setVariables("recoCaloJets_ak5CaloJets");
64 tschum 1.6 p.loop(p.autoVars,cuts,numEv);
65     p.saveCanvases("ak5CaloJets",saveDir);
66 tschum 1.5 p.clearCanvases();
67    
68 tschum 1.7 cuts.clear();
69     cuts.push_back("");
70 tschum 1.5 p.setVariables("recoTracks_generalTracks");
71 tschum 1.6 p.loop(p.autoVars,cuts,numEv);
72     p.saveCanvases("generalTracks",saveDir);
73 tschum 1.5 p.clearCanvases();
74    
75 tschum 1.7 cuts.clear();
76     cuts.push_back("");
77     cuts.push_back("abs(recoTracks_generalTracks.eta())<2.5");
78     cuts.push_back(cuts.back()+string("&&recoTracks_generalTracks.numberOfValidHits()>8"));
79     cuts.push_back(cuts.back()+string("&&abs(recoTracks_generalTracks.d0())<2"));
80     cuts.push_back(cuts.back()+string("&&abs(recoTracks_generalTracks.dz())<20"));
81 tschum 1.1
82 tschum 1.6 p.logY=0;
83     p.loop("recoTracks_generalTracks.d0():recoTracks_generalTracks.phi()",cuts,numEv);
84     p.saveCanvases("generalTracks_2D",saveDir);
85     p.clearCanvases();
86    
87 tschum 1.7 cuts.clear();
88    
89 tschum 1.1
90    
91     }