ViewVC Help

View File | Revision Log | Show Annotations | Root Listing

root/cvsroot/UserCode/tschum/FWlite_Analysis/start.C

Comparing UserCode/tschum/FWlite_Analysis/start.C (file contents):
Revision 1.7 by tschum, Wed Dec 9 13:02:16 2009 UTC vs.
Revision 1.8 by tschum, Thu Dec 10 18:29:27 2009 UTC

#	Line 2 | Line 2
2		gROOT->ProcessLine(".L PlotTool.cc+");
3		PlotTool p;
4
5	<	p.addTrackJets = 0;
5	>	p.addTrackJets = 1;
6		p.addEventInfo = 1;
7	+	p.addTower = 1;
8		//  p.init("","dcap://dcache-ses-cms.desy.de//pnfs/desy.de/cms/tier2/store/data/BeamCommissioning09/MinimumBias/RECO/rereco_FIRSTCOLL_v1/0083/","Events","FIRSTCOLL");
9		// p.init("goldenRun.txt","dcap://dcache-ses-cms.desy.de//pnfs/desy.de/cms/tier2","Events","FIRSTCOLL");
10		// p.init("","/scratch/hh/current/cms/user/schum/run123592/","Events","run123592");
#	Line 34 | Line 35
35		//  p.sameCanv_vars   = true;
36
37
38	<	p.loop(vars,cuts,numEv);
39	<	p.saveCanvases("EventAuxiliary",saveDir);
40	<	p.clearCanvases();
41	<
42	<	p.setVariables("offlinePrimaryVertices");
43	<	p.loop(p.autoVars,cuts,numEv);
44	<	p.saveCanvases("offlinePrimaryVertices",saveDir);
45	<	p.clearCanvases();
46	<
47	<	p.setVariables("offlineBeamSpot");
48	<	p.loop(p.autoVars,cuts,numEv);
49	<	p.saveCanvases("offlineBeamSpot",saveDir);
50	<	p.clearCanvases();
51	<
52	<	p.setVariables("recoCaloMETs_metNoHFHO");
53	<	p.loop(p.autoVars,cuts,numEv);
54	<	p.saveCanvases("metNoHFHO",saveDir);
55	<	p.clearCanvases();
38	>	//   p.loop(vars,cuts,numEv);
39	>	//   p.saveCanvases("EventAuxiliary",saveDir);
40	>	//   p.clearCanvases();
41	>
42	>	//   p.setVariables("offlinePrimaryVertices");
43	>	//   p.loop(p.autoVars,cuts,numEv);
44	>	//   p.saveCanvases("offlinePrimaryVertices",saveDir);
45	>	//   p.clearCanvases();
46	>
47	>	//   p.setVariables("offlineBeamSpot");
48	>	//   p.loop(p.autoVars,cuts,numEv);
49	>	//   p.saveCanvases("offlineBeamSpot",saveDir);
50	>	//   p.clearCanvases();
51	>
52	>	//   p.setVariables("recoCaloMETs_metNoHFHO");
53	>	//   p.loop(p.autoVars,cuts,numEv);
54	>	//   p.saveCanvases("metNoHFHO",saveDir);
55	>	//   p.clearCanvases();
56	>
57	>	//   cuts.clear();
58	>	//   cuts.push_back("");
59	>	//   cuts.push_back("abs(ak5CaloJets.eta())<2.5");
60	>	//   //  cuts.push_back(cuts.back()+string("&&ak5CaloJets.n90()>1"));
61	>	//   cuts.push_back(cuts.back()+string("&&ak5CaloJets.emEnergyFraction()>0.02&&ak5CaloJets.emEnergyFraction()<0.98"));
62	>	//   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"));
63	>
64	>	//   p.setVariables("recoCaloJets_ak5CaloJets");
65	>	//   p.loop(p.autoVars,cuts,numEv);
66	>	//   p.saveCanvases("ak5CaloJets",saveDir);
67	>	//   p.clearCanvases();
68	>
69	>	//   cuts.clear();
70	>	//   cuts.push_back("");
71	>	//   p.setVariables("recoTracks_generalTracks");
72	>	//   p.loop(p.autoVars,cuts,numEv);
73	>	//   p.saveCanvases("generalTracks",saveDir);
74	>	//   p.clearCanvases();
75	>
76	>	//   cuts.clear();
77	>	//   cuts.push_back("");
78	>	//   cuts.push_back("abs(recoTracks_generalTracks.eta())<2.5");
79	>	//   cuts.push_back(cuts.back()+string("&&recoTracks_generalTracks.numberOfValidHits()>8"));
80	>	//   cuts.push_back(cuts.back()+string("&&abs(recoTracks_generalTracks.d0())<2"));
81	>	//   cuts.push_back(cuts.back()+string("&&abs(recoTracks_generalTracks.dz())<20"));
82	>
83	>	//   p.logY=0;
84	>	//   p.loop("recoTracks_generalTracks.d0():recoTracks_generalTracks.phi()",cuts,numEv);
85	>	//   p.saveCanvases("generalTracks_2D",saveDir);
86	>	//   p.clearCanvases();
87
88		cuts.clear();
89	<	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"));
89	>	cuts.push_back("abs(ak5CaloJets.eta())<1.0&&ak5CaloJets.emEnergyFraction()>0.02&&ak5CaloJets.emEnergyFraction()<0.98");
90
91		p.setVariables("recoCaloJets_ak5CaloJets");
92		p.loop(p.autoVars,cuts,numEv);
93	<	p.saveCanvases("ak5CaloJets",saveDir);
93	>	p.saveCanvases("ak5CaloJets_eta1",saveDir);
94		p.clearCanvases();
95
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	–
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	–
96
97		}

Diff Legend

–	Removed lines
+	Added lines
<	Changed lines
>	Changed lines