3 |
|
PlotTool p; |
4 |
|
|
5 |
|
p.addTrackJets = 0; |
6 |
< |
|
6 |
> |
p.addEventInfo = 1; |
7 |
|
// 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 |
|
// p.init("goldenRun.txt","dcap://dcache-ses-cms.desy.de//pnfs/desy.de/cms/tier2","Events","FIRSTCOLL"); |
9 |
|
// p.init("","/scratch/hh/current/cms/user/schum/run123592/","Events","run123592"); |
12 |
|
vector<string> vars; |
13 |
|
vector<string> cuts; |
14 |
|
|
15 |
< |
string saveDir="/afs/desy.de/user/s/schum/public/run_123596_Plots/"; |
15 |
> |
string saveDir="/afs/desy.de/user/s/schum/cms/run_123596_Plots/"; |
16 |
|
int numEv=-1; |
17 |
|
|
18 |
|
vars.push_back("EventAuxiliary.run()"); |
23 |
|
vars.push_back("EventAuxiliary.event()"); |
24 |
|
vars.push_back("EventAuxiliary.experimentType()"); |
25 |
|
vars.push_back("EventAuxiliary.luminosityBlock()"); |
26 |
+ |
vars.push_back("run"); |
27 |
+ |
vars.push_back("event"); |
28 |
+ |
vars.push_back("lumi"); |
29 |
|
|
30 |
|
cuts.push_back(""); |
31 |
< |
p.globalCuts="Level1TriggerScalerss_scalersRawToDigi.lumiSegmentNrLumiSeg()>=69&&Level1TriggerScalerss_scalersRawToDigi.lumiSegmentNrLumiSeg()<=130"; |
31 |
> |
p.globalCuts="lumi>68"; |
32 |
|
|
33 |
|
// p.samePad_trees = false; |
34 |
|
// p.sameCanv_vars = true; |
38 |
|
p.saveCanvases("EventAuxiliary",saveDir); |
39 |
|
p.clearCanvases(); |
40 |
|
|
41 |
< |
p.setVariables("Lumi"); |
39 |
< |
p.loop(p.autoVars,cuts,numEv); |
40 |
< |
p.saveCanvases("Lumi",saveDir); |
41 |
< |
p.clearCanvases(); |
42 |
< |
|
43 |
< |
p.setVariables("Level1"); |
41 |
> |
p.setVariables("offlinePrimaryVertices"); |
42 |
|
p.loop(p.autoVars,cuts,numEv); |
43 |
< |
p.saveCanvases("Level1",saveDir); |
43 |
> |
p.saveCanvases("offlinePrimaryVertices",saveDir); |
44 |
|
p.clearCanvases(); |
45 |
|
|
46 |
< |
p.setVariables("offlinePrimaryVertices"); |
46 |
> |
p.setVariables("offlineBeamSpot"); |
47 |
|
p.loop(p.autoVars,cuts,numEv); |
48 |
< |
p.saveCanvases("offlinePrimaryVertices",saveDir); |
48 |
> |
p.saveCanvases("offlineBeamSpot",saveDir); |
49 |
|
p.clearCanvases(); |
50 |
|
|
51 |
|
p.setVariables("recoCaloMETs_metNoHFHO"); |
53 |
|
p.saveCanvases("metNoHFHO",saveDir); |
54 |
|
p.clearCanvases(); |
55 |
|
|
56 |
+ |
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 |
|
p.setVariables("recoCaloJets_ak5CaloJets"); |
64 |
|
p.loop(p.autoVars,cuts,numEv); |
65 |
|
p.saveCanvases("ak5CaloJets",saveDir); |
66 |
|
p.clearCanvases(); |
67 |
|
|
68 |
< |
|
68 |
> |
cuts.clear(); |
69 |
> |
cuts.push_back(""); |
70 |
|
p.setVariables("recoTracks_generalTracks"); |
71 |
|
p.loop(p.autoVars,cuts,numEv); |
72 |
|
p.saveCanvases("generalTracks",saveDir); |
73 |
|
p.clearCanvases(); |
74 |
|
|
75 |
+ |
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 |
|
|
70 |
– |
p.setVariables("generalTracks"); |
71 |
– |
cuts.push_back("abs(recoTracks_generalTracks.d0())<1"); |
72 |
– |
cuts.push_back("abs(recoTracks_generalTracks.d0())<1&&recoTracks_generalTracks.numberOfValidHits()>7"); |
73 |
– |
cuts.push_back("abs(recoTracks_generalTracks.d0())<1&&recoTracks_generalTracks.numberOfValidHits()>7&&abs(recoTracks_generalTracks.dz())<10"); |
82 |
|
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 |
+ |
cuts.clear(); |
88 |
+ |
|
89 |
|
|
90 |
|
|
91 |
|
} |