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.10 by tschum, Fri Feb 12 15:21:17 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;
7	<	p.addEventInfo = 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");
11	<	p.init("run123596_lumi69-139.txt","","Events","run123596_lumi69-130");
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	>	//   p.varBlockNames += "Clusters";
14	>	//   p.varBlockNames += "EcalHits";
15	>
16	>	p.varBlockNames += "ALL";
17	>	p.recreateTree = true;
18	>
19	>
20	>	p.init("goldenRun.txt","dcap://dcache-ses-cms.desy.de//pnfs/desy.de/cms/tier2","Events","Jan29ReReco");
21
22		vector<string> vars;
23		vector<string> cuts;
24
25	<	string saveDir="/afs/desy.de/user/s/schum/cms/run_123596_Plots/";
25	>	string saveDir="/afs/naf.desy.de/group/cms/scratch/schum/";
26		int numEv=-1;
27
18	–	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	–	vars.push_back("run");
27	–	vars.push_back("event");
28	–	vars.push_back("lumi");
29	–
30	–	cuts.push_back("");
31	–	p.globalCuts="lumi>68";
32	–
33	–	//  p.samePad_trees   = false;
34	–	//  p.sameCanv_vars   = true;
35	–
36	–
37	–	p.loop(vars,cuts,numEv);
38	–	p.saveCanvases("EventAuxiliary",saveDir);
39	–	p.clearCanvases();
40	–
41	–	p.setVariables("offlinePrimaryVertices");
42	–	p.loop(p.autoVars,cuts,numEv);
43	–	p.saveCanvases("offlinePrimaryVertices",saveDir);
44	–	p.clearCanvases();
45	–
46	–	p.setVariables("offlineBeamSpot");
47	–	p.loop(p.autoVars,cuts,numEv);
48	–	p.saveCanvases("offlineBeamSpot",saveDir);
49	–	p.clearCanvases();
50	–
51	–	p.setVariables("recoCaloMETs_metNoHFHO");
52	–	p.loop(p.autoVars,cuts,numEv);
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	–	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();
28
29	<	cuts.clear();
29	>	//   vars.push_back("run");
30	>	//   vars.push_back("event");
31	>	//   vars.push_back("lumi");
32	>	//   cuts.push_back("");
33	>	//   p.loop(vars,cuts,numEv);
34	>
35	>
36	>	//   p.samePad_trees   = false;
37	>	//   p.sameCanv_vars   = true;
38	>
39	>
40	>	p.globalCuts="Trigger1";
41	>	p.logY = 0;
42	>
43	>
44	>	p.loop("TowPhiBinsE:PhiBins","TowTotalE>0",numEv,"","prof");
45	>	p.loop("TowEtaBinsE:EtaBins","TowTotalE>0",numEv,"","prof");
46	>	p.loop("TowRBinsE:RBins"    ,"TowTotalE>0",numEv,"","prof");
47	>
48	>	p.loop("TrackPhiBinsE:PhiBins","TrackTotalE>0",numEv,"","prof");
49	>	p.loop("TrackEtaBinsE:EtaBins","TrackTotalE>0",numEv,"","prof");
50	>	p.loop("TrackRBinsE:RBins"    ,"TrackTotalE>0",numEv,"","prof");
51	>
52
53
54	+	p.saveCanvases("TEST_Jan29ReReco",saveDir);
55
56		}

Diff Legend

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