41 |
|
theHistos.Add(new TH1F("hDistE","hDistE",50,0.,500.)); |
42 |
|
theHistos.Add(new TH1F("hPenetrationB","hPenetrationB",10,0.,10.)); |
43 |
|
theHistos.Add(new TH1F("hPenetrationE","hPenetrationE",10,0.,10.)); |
44 |
+ |
theHistos.Add(new TH1F("hPullX","hPullX",100,-10.,10.)); |
45 |
+ |
theHistos.Add(new TH1F("hDistX","hDistX",100,-100.,100.)); |
46 |
|
|
47 |
|
typedef std::vector<edm::ParameterSet> VPSet; |
48 |
|
VPSet selectorConfigs = theAnaConfig.getParameter<VPSet>("ORedSynchroFilters"); |
53 |
|
f.setHistoPhi(static_cast<TH1F* >(theHistos.FindObject("hDeltaPhi"))); |
54 |
|
f.setHistoEta(static_cast<TH1F* >(theHistos.FindObject("hDeltaEta"))); |
55 |
|
f.setHistoDxy(static_cast<TH1F* >(theHistos.FindObject("hDxy"))); |
56 |
+ |
f.setHistoPullX(static_cast<TH1F* >(theHistos.FindObject("hPullX"))); |
57 |
+ |
f.setHistoDistX(static_cast<TH1F* >(theHistos.FindObject("hDistX"))); |
58 |
+ |
} |
59 |
+ |
if (it->getParameter<std::string>("collection") == "globalMuons") { |
60 |
|
} |
61 |
|
theSynchroFilters.push_back(f); |
62 |
|
} |
107 |
|
if (takeIt) selected.push_back(*it); |
108 |
|
if (takeIt) wasTakeIt = true; |
109 |
|
} |
110 |
+ |
|
111 |
+ |
/* |
112 |
|
std::vector<GlobalPoint> pos = theSynchroFilters[0].positions(); |
113 |
|
std::sort(pos.begin(), pos.end(), LessGP()); |
114 |
|
GlobalPoint last(0.,0.,0.); |
132 |
|
if (barrel && endcap) std::cout << "PROBLEM" << std::endl; |
133 |
|
if (wasTakeIt && barrel ) static_cast<TH1F* >(theHistos.FindObject("hPenetrationB"))->Fill(points); |
134 |
|
if (wasTakeIt && endcap) static_cast<TH1F* >(theHistos.FindObject("hPenetrationE"))->Fill(points); |
135 |
+ |
|
136 |
+ |
*/ |
137 |
|
if (wasTakeIt) static_cast<TH1*>(theHistos.FindObject("hBX"))->Fill( ev.bunchCrossing()); |
138 |
|
if (wasTakeIt) theEventCounter++; |
139 |
|
return selected; |