20 |
|
#include "UserCode/L1RpcTriggerAnalysis/interface/TrackToL1ObjMatcher.h" |
21 |
|
#include "DataFormats/RPCDigi/interface/RPCRawSynchro.h" |
22 |
|
#include "UserCode/L1RpcTriggerAnalysis/interface/ConverterRPCRawSynchroSynchroCountsObj.h" |
23 |
< |
|
23 |
> |
#include "UserCode/L1RpcTriggerAnalysis/interface/TriggerMenuResultObj.h" |
24 |
|
|
25 |
|
#include "TFile.h" |
26 |
|
#include "TTree.h" |
29 |
|
|
30 |
|
|
31 |
|
L1RpcTreeMaker::L1RpcTreeMaker(const edm::ParameterSet& cfg) |
32 |
< |
: theConfig(cfg), theTree(0), event(0), muon(0), track(0), counts(0), |
32 |
> |
: theConfig(cfg), theTree(0), event(0), muon(0), track(0), |
33 |
> |
bitsL1(0), bitsHLT(0), |
34 |
> |
counts(0), |
35 |
|
l1RpcColl(0) , l1OtherColl(0), l1RpcCollEmu(0), |
36 |
|
theCounter(0), |
37 |
|
theBestMuonFinder(cfg.getParameter<edm::ParameterSet>("bestMuonFinder")), |
52 |
|
theTree->Branch("muon","MuonObj",&muon,32000,99); |
53 |
|
theTree->Branch("track", "TrackObj",&track,32000,99); |
54 |
|
|
55 |
+ |
theTree->Branch("bitsL1" ,"TriggerMenuResultObj",&bitsL1 ,32000,99); |
56 |
+ |
theTree->Branch("bitsHLT","TriggerMenuResultObj",&bitsHLT,32000,99); |
57 |
+ |
|
58 |
|
theTree->Branch("counts",&counts); |
59 |
|
theTree->Branch("detsCrossedByMuon",&detsCrossedByMuon); |
60 |
|
theTree->Branch("detsCrossedByMuonDeepInside",&detsCrossedByMuonDeepInside); |
61 |
|
theTree->Branch("detsHitsCompatibleWithMuon",&detsHitsCompatibleWithMuon); |
62 |
|
theTree->Branch("detsSIMU",&detsSIMU); |
63 |
|
theTree->Branch("nDigisCompDets",&nDigisCompDets); |
64 |
+ |
theTree->Branch("clSizeCompDets",&clSizeCompDets); |
65 |
|
|
66 |
|
theTree->Branch("l1RpcColl","L1ObjColl",&l1RpcColl,32000,99); |
67 |
|
theTree->Branch("l1RpcCollEmu","L1ObjColl",&l1RpcCollEmu,32000,99); |
116 |
|
// |
117 |
|
muon = new MuonObj(); |
118 |
|
track = new TrackObj(); |
119 |
+ |
|
120 |
+ |
bitsL1 = new TriggerMenuResultObj(); |
121 |
+ |
bitsHLT = new TriggerMenuResultObj(); |
122 |
+ |
|
123 |
|
// counts = new SynchroCountsObjVect; |
124 |
|
counts = std::vector<SynchroCountsObj>(); |
125 |
|
detsCrossedByMuon = std::vector<uint32_t>(); |
127 |
|
detsHitsCompatibleWithMuon = std::vector<uint32_t>(); |
128 |
|
detsSIMU = std::vector<uint32_t>(); |
129 |
|
nDigisCompDets = std::vector<uint32_t>(); |
130 |
+ |
clSizeCompDets = std::vector<uint32_t>(); |
131 |
|
l1RpcColl = new L1ObjColl; |
132 |
|
l1OtherColl = new L1ObjColl; |
133 |
|
l1RpcCollEmu = new L1ObjColl; |
156 |
|
detsCrossedByMuonDeepInside = theDetHitCollector.compatibleDets( theMuon, ev, es, true); |
157 |
|
if (theConfig.getParameter<bool>("checkDestSIMU")) detsSIMU = theDetHitCollector.compatibleSIMU( theMuon, ev, es); |
158 |
|
nDigisCompDets = theDetHitCollector.nDigisCompDets(detsHitsCompatibleWithMuon, ev, es); |
159 |
+ |
clSizeCompDets = theDetHitCollector.clSizeCompDets(detsHitsCompatibleWithMuon, ev, es); |
160 |
+ |
|
161 |
+ |
for (uint32_t i=0; i< nDigisCompDets.size(); i++) { |
162 |
+ |
if (clSizeCompDets[i] > nDigisCompDets[i]) std::cout <<" PROBLEM, event: "<<theCounter<<" cl:"<<clSizeCompDets[i]<<" nDigis:"<<nDigisCompDets[i]<<std::endl; |
163 |
+ |
} |
164 |
|
} |
165 |
|
|
166 |
|
|
212 |
|
delete event; event = 0; |
213 |
|
delete muon; muon = 0; |
214 |
|
delete track; track = 0; |
215 |
+ |
delete bitsL1; bitsL1= 0; |
216 |
+ |
delete bitsHLT; bitsHLT= 0; |
217 |
|
delete l1RpcColl; l1RpcColl = 0; |
218 |
|
delete l1OtherColl; l1OtherColl = 0; |
219 |
|
delete l1RpcCollEmu; l1RpcCollEmu = 0; |