1 |
|
{ |
2 |
+ |
gROOT->ProcessLine(".L FWliteVariables.cc+"); |
3 |
|
gROOT->ProcessLine(".L PlotTool.cc+"); |
4 |
|
PlotTool p; |
5 |
|
|
6 |
< |
p.addTrackJets = 0; |
6 |
> |
p.varBlockNames += "EventInfo"; |
7 |
> |
p.varBlockNames += "TrackJets"; |
8 |
> |
p.varBlockNames += "CaloTower"; |
9 |
> |
p.varBlockNames += "HitDetInfo"; |
10 |
> |
p.varBlockNames += "DiTrackMass"; |
11 |
> |
p.varBlockNames += "TechTrigger"; |
12 |
> |
p.varBlockNames += "JetContent"; |
13 |
> |
|
14 |
> |
p.recreateTree = true; |
15 |
> |
|
16 |
|
|
17 |
|
// p.init("","dcap://dcache-ses-cms.desy.de//pnfs/desy.de/cms/tier2/store/data/BeamCommissioning09/MinimumBias/RECO/rereco_FIRSTCOLL_v1/0083/","Events","FIRSTCOLL"); |
18 |
< |
p.init("goldenRun.txt","dcap://dcache-ses-cms.desy.de//pnfs/desy.de/cms/tier2","Events","FIRSTCOLL"); |
18 |
> |
p.init("goldenRun.txt","dcap://dcache-ses-cms.desy.de//pnfs/desy.de/cms/tier2","Events","Jan29ReReco"); |
19 |
> |
// p.init("","/scratch/hh/current/cms/user/schum/run123592/","Events","run123592"); |
20 |
> |
// p.init("run123596_lumi69-139.txt","","Events","run123596_lumi69-130"); |
21 |
|
|
22 |
|
vector<string> vars; |
23 |
|
vector<string> cuts; |
24 |
|
|
25 |
< |
vars.push_back("EventAuxiliary.run()"); |
26 |
< |
vars.push_back("EventAuxiliary.isRealData()"); |
27 |
< |
vars.push_back("EventAuxiliary.bunchCrossing()"); |
28 |
< |
vars.push_back("EventAuxiliary.orbitNumber()"); |
29 |
< |
vars.push_back("EventAuxiliary.storeNumber()"); |
30 |
< |
vars.push_back("EventAuxiliary.event()"); |
31 |
< |
vars.push_back("EventAuxiliary.experimentType()"); |
32 |
< |
vars.push_back("EventAuxiliary.luminosityBlock()"); |
25 |
> |
string saveDir="/afs/naf.desy.de/group/cms/scratch/schum/"; |
26 |
> |
int numEv=-1; |
27 |
> |
|
28 |
> |
// vars.push_back("EventAuxiliary.run()"); |
29 |
> |
// vars.push_back("EventAuxiliary.isRealData()"); |
30 |
> |
// vars.push_back("EventAuxiliary.bunchCrossing()"); |
31 |
> |
// vars.push_back("EventAuxiliary.orbitNumber()"); |
32 |
> |
// vars.push_back("EventAuxiliary.storeNumber()"); |
33 |
> |
// vars.push_back("EventAuxiliary.event()"); |
34 |
> |
// vars.push_back("EventAuxiliary.experimentType()"); |
35 |
> |
// vars.push_back("EventAuxiliary.luminosityBlock()"); |
36 |
> |
vars.push_back("run"); |
37 |
> |
vars.push_back("event"); |
38 |
> |
vars.push_back("lumi"); |
39 |
|
|
40 |
|
cuts.push_back(""); |
41 |
< |
// cuts.push_back("EventAuxiliary.run()==122294"); |
42 |
< |
// cuts.push_back("EventAuxiliary.run()==122314"); |
41 |
> |
// p.globalCuts="lumi>20"; |
42 |
> |
|
43 |
> |
|
44 |
|
|
45 |
< |
p.samePad_trees = false; |
46 |
< |
p.sameCanv_vars = true; |
45 |
> |
// p.samePad_trees = false; |
46 |
> |
// p.sameCanv_vars = true; |
47 |
|
|
29 |
– |
p.loop(vars,cuts); |
48 |
|
|
49 |
< |
p.saveCanvases(); |
49 |
> |
p.loop(vars,cuts,numEv); |
50 |
> |
// p.saveCanvases("EventAuxiliary",saveDir); |
51 |
> |
// p.clearCanvases(); |
52 |
> |
|
53 |
> |
// p.setVariables("offlinePrimaryVertices"); |
54 |
> |
// p.loop(p.autoVars,cuts,numEv); |
55 |
> |
// p.saveCanvases("offlinePrimaryVertices",saveDir); |
56 |
> |
// p.clearCanvases(); |
57 |
> |
|
58 |
> |
// p.setVariables("offlineBeamSpot"); |
59 |
> |
// p.loop(p.autoVars,cuts,numEv); |
60 |
> |
// p.saveCanvases("offlineBeamSpot",saveDir); |
61 |
> |
// p.clearCanvases(); |
62 |
> |
|
63 |
> |
// p.setVariables("recoCaloMETs_metNoHFHO"); |
64 |
> |
// p.loop(p.autoVars,cuts,numEv); |
65 |
> |
// p.saveCanvases("metNoHFHO",saveDir); |
66 |
> |
// p.clearCanvases(); |
67 |
> |
|
68 |
> |
// cuts.clear(); |
69 |
> |
// cuts.push_back(""); |
70 |
> |
// cuts.push_back("abs(ak5CaloJets.eta())<2.5"); |
71 |
> |
// // cuts.push_back(cuts.back()+string("&&ak5CaloJets.n90()>1")); |
72 |
> |
// cuts.push_back(cuts.back()+string("&&ak5CaloJets.emEnergyFraction()>0.02&&ak5CaloJets.emEnergyFraction()<0.98")); |
73 |
> |
// 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")); |
74 |
> |
|
75 |
> |
// p.setVariables("recoCaloJets_ak5CaloJets"); |
76 |
> |
// p.loop(p.autoVars,cuts,numEv); |
77 |
> |
// p.saveCanvases("ak5CaloJets",saveDir); |
78 |
> |
// p.clearCanvases(); |
79 |
> |
|
80 |
> |
// cuts.clear(); |
81 |
> |
// cuts.push_back(""); |
82 |
> |
// p.setVariables("recoTracks_generalTracks"); |
83 |
> |
// p.loop(p.autoVars,cuts,numEv); |
84 |
> |
// p.saveCanvases("generalTracks",saveDir); |
85 |
> |
// p.clearCanvases(); |
86 |
> |
|
87 |
> |
// cuts.clear(); |
88 |
> |
// cuts.push_back(""); |
89 |
> |
// cuts.push_back("abs(recoTracks_generalTracks.eta())<2.5"); |
90 |
> |
// cuts.push_back(cuts.back()+string("&&recoTracks_generalTracks.numberOfValidHits()>8")); |
91 |
> |
// cuts.push_back(cuts.back()+string("&&abs(recoTracks_generalTracks.d0())<2")); |
92 |
> |
// cuts.push_back(cuts.back()+string("&&abs(recoTracks_generalTracks.dz())<20")); |
93 |
> |
|
94 |
> |
// p.logY=0; |
95 |
> |
// p.loop("recoTracks_generalTracks.d0():recoTracks_generalTracks.phi()",cuts,numEv); |
96 |
> |
// p.saveCanvases("generalTracks_2D",saveDir); |
97 |
> |
// p.clearCanvases(); |
98 |
> |
|
99 |
> |
// cuts.clear(); |
100 |
> |
// cuts.push_back("abs(ak5CaloJets.eta())<1.0&&ak5CaloJets.emEnergyFraction()>0.02&&ak5CaloJets.emEnergyFraction()<0.98"); |
101 |
> |
|
102 |
> |
// p.setVariables("recoCaloJets_ak5CaloJets"); |
103 |
> |
// p.loop(p.autoVars,cuts,numEv); |
104 |
> |
p.saveCanvases("TEST_Jan29ReReco",saveDir); |
105 |
|
|
106 |
|
} |