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.6 by tschum, Tue Dec 8 10:04:48 2009 UTC vs.
Revision 1.9 by tschum, Wed Feb 10 11:41:58 2010 UTC

#	Line 1 | Line 1
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");
20	>	// p.init("run123596_lumi69-139.txt","","Events","run123596_lumi69-130");
21
22		vector<string> vars;
23		vector<string> cuts;
24
25	<	string saveDir="/afs/desy.de/user/s/schum/public/run_123596_Plots/";
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()");
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	<	p.globalCuts="Level1TriggerScalerss_scalersRawToDigi.lumiSegmentNrLumiSeg()>=69&&Level1TriggerScalerss_scalersRawToDigi.lumiSegmentNrLumiSeg()<=130";
41	>	//  p.globalCuts="lumi>20";
42
30	–	//  p.samePad_trees   = false;
31	–	//  p.sameCanv_vars   = true;
43
44
45	<	p.loop(vars,cuts,numEv);
46	<	p.saveCanvases("EventAuxiliary",saveDir);
36	<	p.clearCanvases();
37	<
38	<	p.setVariables("Lumi");
39	<	p.loop(p.autoVars,cuts,numEv);
40	<	p.saveCanvases("Lumi",saveDir);
41	<	p.clearCanvases();
42	<
43	<	p.setVariables("Level1");
44	<	p.loop(p.autoVars,cuts,numEv);
45	<	p.saveCanvases("Level1",saveDir);
46	<	p.clearCanvases();
47	<
48	<	p.setVariables("offlinePrimaryVertices");
49	<	p.loop(p.autoVars,cuts,numEv);
50	<	p.saveCanvases("offlinePrimaryVertices",saveDir);
51	<	p.clearCanvases();
52	<
53	<	p.setVariables("recoCaloMETs_metNoHFHO");
54	<	p.loop(p.autoVars,cuts,numEv);
55	<	p.saveCanvases("metNoHFHO",saveDir);
56	<	p.clearCanvases();
57	<
58	<	p.setVariables("recoCaloJets_ak5CaloJets");
59	<	p.loop(p.autoVars,cuts,numEv);
60	<	p.saveCanvases("ak5CaloJets",saveDir);
61	<	p.clearCanvases();
62	<
63	<
64	<	p.setVariables("recoTracks_generalTracks");
65	<	p.loop(p.autoVars,cuts,numEv);
66	<	p.saveCanvases("generalTracks",saveDir);
67	<	p.clearCanvases();
68	<
69	<
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");
74	<	p.logY=0;
75	<	p.loop("recoTracks_generalTracks.d0():recoTracks_generalTracks.phi()",cuts,numEv);
76	<	p.saveCanvases("generalTracks_2D",saveDir);
77	<	p.clearCanvases();
45	>	//  p.samePad_trees   = false;
46	>	//  p.sameCanv_vars   = true;
47
48
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		}

Diff Legend

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