| 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 |
> |
// p.varBlockNames += "Clusters"; |
| 14 |
> |
// p.varBlockNames += "EcalHits"; |
| 15 |
> |
// p.varBlockNames += "HLT"; |
| 16 |
> |
|
| 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"); |
| 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"); |
| 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; |
| 24 |
> |
vector<string> cuts; |
| 25 |
|
|
| 26 |
< |
vars.push_back("EventAuxiliary.run()"); |
| 27 |
< |
vars.push_back("EventAuxiliary.isRealData()"); |
| 16 |
< |
vars.push_back("EventAuxiliary.bunchCrossing()"); |
| 17 |
< |
vars.push_back("EventAuxiliary.orbitNumber()"); |
| 18 |
< |
vars.push_back("EventAuxiliary.storeNumber()"); |
| 19 |
< |
vars.push_back("EventAuxiliary.event()"); |
| 20 |
< |
vars.push_back("EventAuxiliary.experimentType()"); |
| 21 |
< |
vars.push_back("EventAuxiliary.luminosityBlock()"); |
| 22 |
< |
|
| 23 |
< |
// cuts.push_back(""); |
| 24 |
< |
// cuts.push_back("EventAuxiliary.run()==122294"); |
| 25 |
< |
// cuts.push_back("EventAuxiliary.run()==122314"); |
| 26 |
< |
|
| 27 |
< |
// p.samePad_trees = false; |
| 28 |
< |
// p.sameCanv_vars = true; |
| 29 |
< |
|
| 30 |
< |
|
| 31 |
< |
p.loop(vars); |
| 32 |
< |
p.saveCanvases("EventAuxiliary","/afs/desy.de/user/s/schum/public/run_123592_Plots/"); |
| 33 |
< |
p.clearCanvases(); |
| 34 |
< |
|
| 35 |
< |
p.setVariables("Lumi"); |
| 36 |
< |
p.loop(p.autoVars); |
| 37 |
< |
p.saveCanvases("Lumi","/afs/desy.de/user/s/schum/public/run_123592_Plots/"); |
| 38 |
< |
p.clearCanvases(); |
| 39 |
< |
|
| 40 |
< |
p.setVariables("Level1"); |
| 41 |
< |
p.loop(p.autoVars); |
| 42 |
< |
p.saveCanvases("Level1","/afs/desy.de/user/s/schum/public/run_123592_Plots/"); |
| 43 |
< |
p.clearCanvases(); |
| 44 |
< |
|
| 45 |
< |
p.setVariables("offlinePrimaryVertices"); |
| 46 |
< |
p.loop(p.autoVars); |
| 47 |
< |
p.saveCanvases("offlinePrimaryVertices","/afs/desy.de/user/s/schum/public/run_123592_Plots/"); |
| 48 |
< |
p.clearCanvases(); |
| 49 |
< |
|
| 50 |
< |
p.setVariables("recoCaloMETs_metNoHFHO"); |
| 51 |
< |
p.loop(p.autoVars); |
| 52 |
< |
p.saveCanvases("metNoHFHO","/afs/desy.de/user/s/schum/public/run_123592_Plots/"); |
| 53 |
< |
p.clearCanvases(); |
| 54 |
< |
|
| 55 |
< |
p.setVariables("recoCaloJets_ak5CaloJets"); |
| 56 |
< |
p.loop(p.autoVars); |
| 57 |
< |
p.saveCanvases("ak5CaloJets","/afs/desy.de/user/s/schum/public/run_123592_Plots/"); |
| 58 |
< |
p.clearCanvases(); |
| 59 |
< |
|
| 60 |
< |
|
| 61 |
< |
p.setVariables("recoTracks_generalTracks"); |
| 62 |
< |
p.loop(p.autoVars); |
| 63 |
< |
p.saveCanvases("recoTracks_generalTracks","/afs/desy.de/user/s/schum/public/run_123592_Plots/"); |
| 64 |
< |
p.clearCanvases(); |
| 26 |
> |
string saveDir="/afs/naf.desy.de/group/cms/scratch/schum/"; |
| 27 |
> |
int numEv=-1; |
| 28 |
|
|
| 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 |
|
} |
