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.9 by tschum, Wed Feb 10 11:41:58 2010 UTC vs.
Revision 1.11 by tschum, Wed Feb 17 15:42:23 2010 UTC

#	Line 3 | Line 3
3		gROOT->ProcessLine(".L PlotTool.cc+");
4		PlotTool p;
5
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";
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	>	//   p.varBlockNames += "HLT";
16
17	<	p.recreateTree = true;
17	>	p.varBlockNames += "ALL";
18	>	//   p.recreateTree = true;
19
20
21	<	//  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","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	>	p.init("runList900GeV.txt","dcap://dcache-ses-cms.desy.de//pnfs/desy.de/cms/tier2","Events","Jan29ReReco_900GeV");
22
23		vector<string> vars;
24		vector<string> cuts;
#	Line 25 | Line 26 | p.init("goldenRun.txt","dcap://dcache-se
26		string saveDir="/afs/naf.desy.de/group/cms/scratch/schum/";
27		int numEv=-1;
28
29	<	//   vars.push_back("EventAuxiliary.run()");
30	<	//   vars.push_back("EventAuxiliary.isRealData()");
31	<	//   vars.push_back("EventAuxiliary.bunchCrossing()");
32	<	//   vars.push_back("EventAuxiliary.orbitNumber()");
33	<	//   vars.push_back("EventAuxiliary.storeNumber()");
34	<	//   vars.push_back("EventAuxiliary.event()");
35	<	//   vars.push_back("EventAuxiliary.experimentType()");
36	<	//   vars.push_back("EventAuxiliary.luminosityBlock()");
37	<	vars.push_back("run");
38	<	vars.push_back("event");
39	<	vars.push_back("lumi");
40	<
41	<	cuts.push_back("");
42	<	//  p.globalCuts="lumi>20";
43	<
44	<
45	<
46	<	//  p.samePad_trees   = false;
47	<	//  p.sameCanv_vars   = true;
48	<
49	<
50	<	p.loop(vars,cuts,numEv);
51	<	//   p.saveCanvases("EventAuxiliary",saveDir);
52	<	//   p.clearCanvases();
53	<
54	<	//   p.setVariables("offlinePrimaryVertices");
55	<	//   p.loop(p.autoVars,cuts,numEv);
56	<	//   p.saveCanvases("offlinePrimaryVertices",saveDir);
57	<	//   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);
29	>
30	>	//   vars.push_back("run");
31	>	//   vars.push_back("event");
32	>	//   vars.push_back("lumi");
33	>	cuts.push_back("recoTracks_generalTracks@.size()<10||Sum$(recoTracks_generalTracks.quality(reco::Track::highPurity))/recoTracks_generalTracks@.size()>0.25"); //rejection of beam scraping events (MONSTER Track Events)
34	>
35	>	cuts.push_back("Sum$(offlinePrimaryVertices.isFake()==0&&offlinePrimaryVertices.ndof()>4&&abs(offlinePrimaryVertices.z())<=15&&offlinePrimaryVertices.position().Rho()<=2) >= 1"); //event selection on primary vertices
36	>	//   p.loop(vars,cuts,numEv);
37	>
38	>
39	>	//   p.samePad_trees   = false;
40	>	//   p.sameCanv_vars   = true;
41	>
42	>
43	>	p.globalCuts="TriggerDefault";
44	>	p.logY = 0;
45	>
46	>
47	>	p.loop("TowPhiBinsE:PhiBins","TowTotalE>0",numEv,"","prof");
48	>	p.loop("TowEtaBinsE:EtaBins","TowTotalE>0",numEv,"","prof");
49	>	p.loop("TowRBinsE:RBins"    ,"TowTotalE>0",numEv,"","prof");
50	>
51	>	p.loop("TrackPhiBinsE:PhiBins","TrackTotalE>0",numEv,"","prof");
52	>	p.loop("TrackEtaBinsE:EtaBins","TrackTotalE>0",numEv,"","prof");
53	>	p.loop("TrackRBinsE:RBins"    ,"TrackTotalE>0",numEv,"","prof");
54	>
55	>
56	>
57	>	p.saveCanvases("TEST_Jan29ReReco_900GeV",saveDir);
58
59		}

Diff Legend

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