ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/tschum/FWlite_Analysis/start.C
(Generate patch)

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.7 by tschum, Wed Dec 9 13:02:16 2009 UTC

# Line 3 | Line 3
3    PlotTool p;
4  
5    p.addTrackJets = 0;
6 <
6 >  p.addEventInfo = 1;
7    //  p.init("","dcap://dcache-ses-cms.desy.de//pnfs/desy.de/cms/tier2/store/data/BeamCommissioning09/MinimumBias/RECO/rereco_FIRSTCOLL_v1/0083/","Events","FIRSTCOLL");
8    // p.init("goldenRun.txt","dcap://dcache-ses-cms.desy.de//pnfs/desy.de/cms/tier2","Events","FIRSTCOLL");
9    // p.init("","/scratch/hh/current/cms/user/schum/run123592/","Events","run123592");
# Line 12 | Line 12
12    vector<string> vars;
13    vector<string> cuts;
14  
15 <  string saveDir="/afs/desy.de/user/s/schum/public/run_123596_Plots/";
15 >  string saveDir="/afs/desy.de/user/s/schum/cms/run_123596_Plots/";
16    int numEv=-1;
17  
18    vars.push_back("EventAuxiliary.run()");
# Line 23 | Line 23
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="Level1TriggerScalerss_scalersRawToDigi.lumiSegmentNrLumiSeg()>=69&&Level1TriggerScalerss_scalersRawToDigi.lumiSegmentNrLumiSeg()<=130";
31 >  p.globalCuts="lumi>68";
32  
33    //  p.samePad_trees   = false;
34    //  p.sameCanv_vars   = true;
# Line 35 | Line 38
38    p.saveCanvases("EventAuxiliary",saveDir);
39    p.clearCanvases();
40  
41 <  p.setVariables("Lumi");
39 <  p.loop(p.autoVars,cuts,numEv);
40 <  p.saveCanvases("Lumi",saveDir);
41 <  p.clearCanvases();
42 <
43 <  p.setVariables("Level1");
41 >  p.setVariables("offlinePrimaryVertices");
42    p.loop(p.autoVars,cuts,numEv);
43 <  p.saveCanvases("Level1",saveDir);
43 >  p.saveCanvases("offlinePrimaryVertices",saveDir);
44    p.clearCanvases();
45  
46 <  p.setVariables("offlinePrimaryVertices");
46 >  p.setVariables("offlineBeamSpot");
47    p.loop(p.autoVars,cuts,numEv);
48 <  p.saveCanvases("offlinePrimaryVertices",saveDir);
48 >  p.saveCanvases("offlineBeamSpot",saveDir);
49    p.clearCanvases();
50  
51    p.setVariables("recoCaloMETs_metNoHFHO");
# Line 55 | Line 53
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 <
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  
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");
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  
90  
91   }

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines