| 1 |
{
|
| 2 |
gROOT->ProcessLine(".L FWliteVariables.cc+");
|
| 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";
|
| 13 |
// p.varBlockNames += "Clusters";
|
| 14 |
// p.varBlockNames += "EcalHits";
|
| 15 |
// p.varBlockNames += "HLT";
|
| 16 |
|
| 17 |
p.varBlockNames += "ALL";
|
| 18 |
// p.recreateTree = true;
|
| 19 |
|
| 20 |
// p.init("runList900GeV.txt","dcap://dcache-ses-cms.desy.de//pnfs/desy.de/cms/tier2","Events","Jan29ReReco_900GeV");
|
| 21 |
// p.init("goldenRun.txt","dcap://dcache-ses-cms.desy.de//pnfs/desy.de/cms/tier2","Events","Jan29ReReco_2360GeV");
|
| 22 |
|
| 23 |
p.init("Feb9ReReco_2360GeV.txt","dcap://dcache-ses-cms.desy.de//pnfs/desy.de/cms/tier2","Events","Feb9ReReco_2360GeV");
|
| 24 |
// p.init("620C11D6-2F0D-DF11-964F-001A928116B8","dcap://dcache-ses-cms.desy.de//pnfs/desy.de/cms/tier2/store/data/BeamCommissioning09/MinimumBias/RECO/Jan29ReReco-v2/0017/","Events","TEST");
|
| 25 |
|
| 26 |
vector<string> vars;
|
| 27 |
vector<string> cuts;
|
| 28 |
|
| 29 |
string saveDir="/afs/naf.desy.de/group/cms/scratch/schum/";
|
| 30 |
int numEv=-1;
|
| 31 |
|
| 32 |
// vars.push_back("run");
|
| 33 |
// vars.push_back("event");
|
| 34 |
// vars.push_back("lumi");
|
| 35 |
// 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)
|
| 36 |
|
| 37 |
// cuts.push_back("Sum$(offlinePrimaryVertices.isFake()==0&&offlinePrimaryVertices.ndof()>4&&abs(offlinePrimaryVertices.z())<=15&&offlinePrimaryVertices.position().Rho()<=2) >= 1"); //event selection on primary vertices
|
| 38 |
// p.loop(vars,cuts,numEv);
|
| 39 |
|
| 40 |
|
| 41 |
// p.samePad_trees = false;
|
| 42 |
// p.sameCanv_vars = true;
|
| 43 |
|
| 44 |
|
| 45 |
p.globalCuts="TriggerDefault";
|
| 46 |
p.globalCuts+="&&(nTracks<10||nTracksHighPurity/nTracks>0.25)";
|
| 47 |
p.globalCuts+="&&Sum$(offlinePrimaryVertices.isFake()==0&&offlinePrimaryVertices.ndof()>4&&abs(offlinePrimaryVertices.z())<=15&&offlinePrimaryVertices.position().Rho()<=2) >= 1";
|
| 48 |
p.logY = 0;
|
| 49 |
|
| 50 |
// cuts.push_back("TowTotalE>0");
|
| 51 |
// cuts.push_back("NTowLeadJet>1");
|
| 52 |
cuts.push_back("NTowLeadJet>5");
|
| 53 |
// cuts.push_back("ak5CaloJets[0].pt()>5");
|
| 54 |
// cuts.push_back("ak5CaloJets[0].pt()>1&&abs(ak5CaloJets[0].eta())<1.3");
|
| 55 |
|
| 56 |
p.loop("TowPhiBinsE:PhiBins",cuts,numEv,"","prof");
|
| 57 |
p.loop("TowEtaBinsE:EtaBins",cuts,numEv,"","prof");
|
| 58 |
p.loop("TowRBinsE:RBins" ,cuts,numEv,"","prof");
|
| 59 |
|
| 60 |
p.loop("TowPhiBinsEm:PhiBins",cuts,numEv,"","prof");
|
| 61 |
p.loop("TowEtaBinsEm:EtaBins",cuts,numEv,"","prof");
|
| 62 |
p.loop("TowRBinsEm:RBins" ,cuts,numEv,"","prof");
|
| 63 |
|
| 64 |
p.loop("TowPhiBinsHad:PhiBins",cuts,numEv,"","prof");
|
| 65 |
p.loop("TowEtaBinsHad:EtaBins",cuts,numEv,"","prof");
|
| 66 |
p.loop("TowRBinsHad:RBins" ,cuts,numEv,"","prof");
|
| 67 |
|
| 68 |
p.loop("TrackPhiBinsE:PhiBins","TrackTotalE>0",numEv,"","prof");
|
| 69 |
p.loop("TrackEtaBinsE:EtaBins","TrackTotalE>0",numEv,"","prof");
|
| 70 |
p.loop("TrackRBinsE:RBins" ,"TrackTotalE>0",numEv,"","prof");
|
| 71 |
|
| 72 |
p.loop("TrackOuterPhiBinsE:PhiBins","TrackOuterTotalE>0",numEv,"","prof");
|
| 73 |
p.loop("TrackOuterEtaBinsE:EtaBins","TrackOuterTotalE>0",numEv,"","prof");
|
| 74 |
p.loop("TrackOuterRBinsE:RBins" ,"TrackOuterTotalE>0",numEv,"","prof");
|
| 75 |
|
| 76 |
|
| 77 |
|
| 78 |
|
| 79 |
p.saveCanvases("JetContent2_Feb9ReReco_2360GeV",saveDir);
|
| 80 |
|
| 81 |
// p.saveCanvases("TEST_Jan29ReReco_900GeV",saveDir);
|
| 82 |
}
|
