| 23 |
|
#include "UserCode/L1RpcTriggerAnalysis/interface/SynchroCountsObj.h" |
| 24 |
|
#include "UserCode/L1RpcTriggerAnalysis/interface/DetCluDigiObj.h" |
| 25 |
|
#include "UserCode/L1RpcTriggerAnalysis/interface/TriggerMenuResultObj.h" |
| 26 |
+ |
#include "UserCode/L1RpcTriggerAnalysis/interface/HitSpecObj.h" |
| 27 |
|
|
| 28 |
|
|
| 29 |
|
#include "UserCode/L1RpcTriggerAnalysis/interface/ConverterRPCRawSynchroSynchroCountsObj.h" |
| 30 |
|
|
| 31 |
|
L1RpcTreeAnalysis::L1RpcTreeAnalysis(const edm::ParameterSet & cfg) |
| 32 |
|
: theConfig(cfg), |
| 33 |
< |
theAnaMuonDistribution( cfg.getParameter<edm::ParameterSet>("anaMuonDistribution") ), |
| 34 |
< |
theAnaTimingL1 (cfg.getParameter<edm::ParameterSet>("anaTimingL1") ), |
| 35 |
< |
theAnaMenu(cfg.getParameter<edm::ParameterSet>("anaMenu") ), |
| 36 |
< |
theAnaEvent(cfg.getParameter<edm::ParameterSet>("anaEvent") ) |
| 37 |
< |
{ } |
| 33 |
> |
theAnaMuonDistribution(0), |
| 34 |
> |
theAnaRpcVsOth(0), |
| 35 |
> |
theAnaRpcMisc(0), |
| 36 |
> |
theAnaEff(0), |
| 37 |
> |
theAnaDet(0), |
| 38 |
> |
theAnaEmu(0), |
| 39 |
> |
theAnaSynch(0), |
| 40 |
> |
theAnaClu(0), |
| 41 |
> |
theAnaTimingL1(0), |
| 42 |
> |
theAnaMenu(0), |
| 43 |
> |
theAnaEvent(0), |
| 44 |
> |
theAnaDigiSpec(0), |
| 45 |
> |
theAnaHitSpec(0), |
| 46 |
> |
thePatternProducer(0), |
| 47 |
> |
thePatternProvider(0) |
| 48 |
> |
{ |
| 49 |
> |
if (theConfig.exists("anaMuonDistribution")) theAnaMuonDistribution = new AnaMuonDistribution( cfg.getParameter<edm::ParameterSet>("anaMuonDistribution") ); |
| 50 |
> |
if (theConfig.exists("anaMenu")) theAnaMenu = new AnaMenu(theConfig.getParameter<edm::ParameterSet>("anaMenu")); |
| 51 |
> |
if (theConfig.exists("anaTimingL1")) theAnaTimingL1 = new AnaTimingL1( theConfig.getParameter<edm::ParameterSet>("anaTimingL1") ); |
| 52 |
> |
if (theConfig.exists("anaEvent")) theAnaEvent = new AnaEvent(cfg.getParameter<edm::ParameterSet>("anaEvent") ); |
| 53 |
> |
if (theConfig.exists("anaDigiSpec")) theAnaDigiSpec = new AnaDigiSpec(cfg.getParameter<edm::ParameterSet>("anaDigiSpec")); |
| 54 |
> |
if (theConfig.exists("anaHitSpec")) theAnaHitSpec = new AnaHitSpec(cfg.getParameter<edm::ParameterSet>("anaHitSpec")); |
| 55 |
> |
if (theConfig.exists("patternProducer")) thePatternProducer = new PatternManager(cfg.getParameter<edm::ParameterSet>("patternProducer")); |
| 56 |
> |
if (theConfig.exists("patternProvider")) thePatternProvider = new PatternManager(cfg.getParameter<edm::ParameterSet>("patternProvider")); |
| 57 |
> |
|
| 58 |
> |
} |
| 59 |
|
|
| 60 |
|
void L1RpcTreeAnalysis::beginJob() |
| 61 |
|
{ |
| 62 |
|
theHistos.SetOwner(); |
| 63 |
|
|
| 64 |
< |
theAnaMuonDistribution.init(theHistos); |
| 65 |
< |
theAnaRpcVsOth.init(theHistos); |
| 66 |
< |
theAnaRpcMisc.init(theHistos); |
| 67 |
< |
theAnaEff.init(theHistos); |
| 68 |
< |
theAnaDet.init(theHistos); |
| 69 |
< |
theAnaEmu.init(theHistos); |
| 70 |
< |
theAnaSynch.init(theHistos); |
| 71 |
< |
theAnaClu.init(theHistos); |
| 72 |
< |
theAnaTimingL1.init(theHistos); |
| 73 |
< |
theAnaEvent.init(theHistos); |
| 74 |
< |
theAnaMenu.init(theHistos); |
| 64 |
> |
if (theAnaMuonDistribution) theAnaMuonDistribution->init(theHistos); |
| 65 |
> |
if (theAnaRpcVsOth) theAnaRpcVsOth->init(theHistos); |
| 66 |
> |
if (theAnaRpcMisc) theAnaRpcMisc->init(theHistos); |
| 67 |
> |
if (theAnaEff) theAnaEff->init(theHistos); |
| 68 |
> |
if (theAnaDet) theAnaDet->init(theHistos); |
| 69 |
> |
if (theAnaEmu) theAnaEmu->init(theHistos); |
| 70 |
> |
if (theAnaSynch) theAnaSynch->init(theHistos); |
| 71 |
> |
if (theAnaClu) theAnaClu->init(theHistos); |
| 72 |
> |
if (theAnaTimingL1) theAnaTimingL1->init(theHistos); |
| 73 |
> |
if (theAnaEvent) theAnaEvent->init(theHistos); |
| 74 |
> |
if (theAnaMenu) theAnaMenu->init(theHistos); |
| 75 |
> |
if (theAnaDigiSpec) theAnaDigiSpec->init(theHistos); |
| 76 |
> |
if (theAnaHitSpec) theAnaHitSpec->init(theHistos); |
| 77 |
> |
|
| 78 |
> |
if (thePatternProvider) thePatternProvider->beginJob(); |
| 79 |
|
} |
| 80 |
|
|
| 81 |
|
void L1RpcTreeAnalysis::beginRun(const edm::Run& ru, const edm::EventSetup& es) |
| 82 |
|
{ |
| 83 |
< |
theAnaSynch.beginRun(ru,es); |
| 83 |
> |
if (theAnaSynch) theAnaSynch->beginRun(ru,es); |
| 84 |
|
} |
| 85 |
|
|
| 86 |
|
void L1RpcTreeAnalysis::analyze(const edm::Event&, const edm::EventSetup&) |
| 111 |
|
std::vector<uint32_t> *detsCrossedByMuonDeepInside = 0; |
| 112 |
|
std::vector<DetCluDigiObj> *detsHitsCompatibleWithMuon = 0; |
| 113 |
|
std::vector<uint32_t> *detsSIMU = 0; |
| 114 |
+ |
std::vector<std::pair<uint32_t, uint32_t> > *digSpec = 0; |
| 115 |
|
|
| 116 |
|
EventObj * event = 0; |
| 117 |
|
MuonObj * muon = 0; |
| 118 |
+ |
TrackObj * simu = 0; |
| 119 |
|
|
| 120 |
|
TriggerMenuResultObj *bitsL1 = 0; |
| 121 |
|
TriggerMenuResultObj *bitsHLT = 0; |
| 125 |
|
TBranch *bdetsCrossedByMuonDeepInside =0; |
| 126 |
|
TBranch *bdetsHitsCompatibleWithMuon = 0; |
| 127 |
|
TBranch *bdetsSIMU =0; |
| 128 |
+ |
TBranch *bdigSpec = 0; |
| 129 |
|
|
| 130 |
|
L1ObjColl* l1ObjColl = 0; |
| 131 |
+ |
HitSpecObj* hitSpec = 0; |
| 132 |
+ |
|
| 133 |
|
|
| 134 |
|
chain.SetBranchAddress("event",&event); |
| 135 |
|
chain.SetBranchAddress("muon",&muon); |
| 136 |
+ |
chain.SetBranchAddress("simu",&simu); |
| 137 |
|
|
| 138 |
|
chain.SetBranchAddress("bitsL1",&bitsL1); |
| 139 |
|
chain.SetBranchAddress("bitsHLT",&bitsHLT); |
| 143 |
|
chain.SetBranchAddress("detsCrossedByMuonDeepInside",&detsCrossedByMuonDeepInside,&bdetsCrossedByMuonDeepInside); |
| 144 |
|
chain.SetBranchAddress("detsHitsCompatibleWithMuon",&detsHitsCompatibleWithMuon,&bdetsHitsCompatibleWithMuon); |
| 145 |
|
chain.SetBranchAddress("detsSIMU",&detsSIMU,&bdetsSIMU); |
| 146 |
+ |
chain.SetBranchAddress("digSpec",&digSpec,&bdigSpec); |
| 147 |
|
|
| 148 |
|
chain.SetBranchAddress("l1ObjColl",&l1ObjColl); |
| 149 |
+ |
chain.SetBranchAddress("hitSpec",&hitSpec); |
| 150 |
|
|
| 151 |
|
|
| 152 |
|
// |
| 162 |
|
unsigned int lastRun = 0; |
| 163 |
|
for (int ev=0; ev<nentries; ev++) { |
| 164 |
|
chain.GetEntry(ev); |
| 165 |
< |
theAnaMenu.updateMenu(bitsL1->names, bitsHLT->names); |
| 165 |
> |
if (theAnaMenu) theAnaMenu->updateMenu(bitsL1->names, bitsHLT->names); |
| 166 |
|
|
| 167 |
< |
if (lastRun != (*event).run) { |
| 167 |
> |
if ( (lastRun != (*event).run) || (ev/10000*10000==ev) ) { |
| 168 |
|
// if (true) { |
| 169 |
+ |
// if (! ((*event).run==204601 && (*event).id ==109463402)) { continue; |
| 170 |
|
lastRun = (*event).run; |
| 171 |
|
std::cout <<"RUN:" << std::setw(7) << (*event).run |
| 172 |
|
<<" event:" << std::setw(8) << ev |
| 196 |
|
|
| 197 |
|
// EVENT NUMBER, BX structure etc. |
| 198 |
|
EventObjBXExtra eventBx(*event); |
| 199 |
< |
if ( !theAnaEvent.filter(&eventBx) && theConfig.getParameter<bool>("filterByAnaEvent") ) continue; |
| 199 |
> |
if ( theAnaEvent && !theAnaEvent->filter(&eventBx) && theConfig.getParameter<bool>("filterByAnaEvent") ) continue; |
| 200 |
|
// ANALYSE AND FILTER KINEMCTICS |
| 201 |
< |
if ( !theAnaMuonDistribution.filter(muon) && theConfig.getParameter<bool>("filterByAnaMuonDistribution") ) continue; |
| 201 |
> |
if ( theAnaMuonDistribution && !theAnaMuonDistribution->filter(muon) && theConfig.getParameter<bool>("filterByAnaMuonDistribution") ) continue; |
| 202 |
|
// ANALYSE AND FILTER TRIGGER MENU |
| 203 |
< |
if ( !theAnaMenu.filter(event, muon, bitsL1, bitsHLT) && theConfig.getParameter<bool>("filterByAnaMenu") ) continue; |
| 169 |
< |
|
| 170 |
< |
// theAnaRpcVsOth.run(muon,l1ObjColl); |
| 171 |
< |
// theAnaEff.run(muon, l1ObjColl); |
| 172 |
< |
// theAnaRpcMisc.run(event,muon,l1ObjColl); |
| 173 |
< |
// theAnaDet.run( muon, *detsHitsCompatibleWithMuon, *detsCrossedByMuon, *detsCrossedByMuonDeepInside); |
| 174 |
< |
// theAnaEmu.run ( event, muon, l1ObjColl); |
| 175 |
< |
// theAnaSynch.run( event, muon, ConverterRPCRawSynchroSynchroCountsObj::toRawSynchro( *counts)); |
| 176 |
< |
// theAnaClu.run( event, muon, l1ObjColl, *detsHitsCompatibleWithMuon); |
| 177 |
< |
theAnaTimingL1.run( &eventBx, muon, l1ObjColl); |
| 203 |
> |
if ( theAnaMenu && !theAnaMenu->filter(event, muon, bitsL1, bitsHLT) && theConfig.getParameter<bool>("filterByAnaMenu") ) continue; |
| 204 |
|
|
| 205 |
+ |
// ANALYSES |
| 206 |
+ |
if (theAnaRpcVsOth) theAnaRpcVsOth->run(muon,l1ObjColl); |
| 207 |
+ |
if (theAnaEff) theAnaEff->run(muon, l1ObjColl); |
| 208 |
+ |
if (theAnaRpcMisc) theAnaRpcMisc->run(event,muon,l1ObjColl); |
| 209 |
+ |
if (theAnaDet) theAnaDet->run( muon, *detsHitsCompatibleWithMuon, *detsCrossedByMuon, *detsCrossedByMuonDeepInside); |
| 210 |
+ |
if (theAnaEmu) theAnaEmu->run ( event, muon, l1ObjColl); |
| 211 |
+ |
if (theAnaSynch) theAnaSynch->run( event, muon, ConverterRPCRawSynchroSynchroCountsObj::toRawSynchro( *counts)); |
| 212 |
+ |
if (theAnaClu) theAnaClu->run( event, muon, l1ObjColl, *detsHitsCompatibleWithMuon); |
| 213 |
+ |
if (theAnaTimingL1) theAnaTimingL1->run( &eventBx, muon, l1ObjColl); |
| 214 |
+ |
|
| 215 |
+ |
// HITPATTERN ANALYSES & MC EFFICIENCY |
| 216 |
+ |
if (theAnaHitSpec) theAnaHitSpec->run(event, simu, hitSpec); |
| 217 |
+ |
if (theAnaDigiSpec) theAnaDigiSpec->run(event, simu, hitSpec, *digSpec); |
| 218 |
+ |
if (thePatternProducer) thePatternProducer->run(event, simu, hitSpec, *digSpec); |
| 219 |
+ |
if (thePatternProvider) thePatternProvider->check(event, simu, hitSpec, *digSpec); |
| 220 |
|
} |
| 221 |
|
} |
| 222 |
|
|
| 223 |
|
void L1RpcTreeAnalysis::endJob() |
| 224 |
|
{ |
| 225 |
|
std::cout <<"ENDJOB, summaries:"<<std::endl; |
| 226 |
< |
TGraph* hGraph_DetEff = theAnaDet.resume(); |
| 227 |
< |
theAnaRpcMisc.resume(theHistos); |
| 228 |
< |
TGraph* hGraph_RunClu = theAnaClu.resume(); |
| 229 |
< |
theAnaTimingL1.resume(theHistos); |
| 230 |
< |
theAnaMenu.resume(theHistos); |
| 231 |
< |
theAnaEvent.resume(theHistos); |
| 232 |
< |
|
| 233 |
< |
theAnaSynch.endJob(); |
| 226 |
> |
TGraph* hGraph_DetEff = (theAnaDet) ? theAnaDet->resume() : 0; |
| 227 |
> |
TGraph* hGraph_RunClu = (theAnaClu) ? theAnaClu->resume(): 0; |
| 228 |
> |
if (theAnaRpcMisc) theAnaRpcMisc->resume(theHistos); |
| 229 |
> |
if (theAnaTimingL1) theAnaTimingL1->resume(theHistos); |
| 230 |
> |
if (theAnaMenu) theAnaMenu->resume(theHistos); |
| 231 |
> |
if (theAnaMenu) theAnaEvent->resume(theHistos); |
| 232 |
> |
if (theAnaHitSpec) theAnaHitSpec->resume(theHistos); |
| 233 |
> |
if (theAnaDigiSpec) theAnaDigiSpec->resume(theHistos); |
| 234 |
> |
|
| 235 |
> |
if (theAnaSynch) theAnaSynch->endJob(); |
| 236 |
> |
if (thePatternProducer) thePatternProducer->endJob(); |
| 237 |
> |
if (thePatternProvider) thePatternProvider->endJob(); |
| 238 |
|
|
| 239 |
|
std::string histoFile = theConfig.getParameter<std::string>("histoFileName"); |
| 240 |
|
TFile f(histoFile.c_str(),"RECREATE"); |
| 241 |
|
theHistos.Write(); |
| 242 |
< |
hGraph_DetEff->Write("hGraph_DetEff"); |
| 243 |
< |
hGraph_RunClu->Write("hGraph_RunClu"); |
| 242 |
> |
if (hGraph_DetEff) hGraph_DetEff->Write("hGraph_DetEff"); |
| 243 |
> |
if (hGraph_RunClu) hGraph_RunClu->Write("hGraph_RunClu"); |
| 244 |
|
f.Close(); |
| 245 |
|
std::cout <<"END"<<std::endl; |
| 246 |
+ |
|
| 247 |
+ |
delete theAnaMuonDistribution; |
| 248 |
+ |
delete theAnaRpcVsOth; |
| 249 |
+ |
delete theAnaRpcMisc; |
| 250 |
+ |
delete theAnaEff; |
| 251 |
+ |
delete theAnaDet; |
| 252 |
+ |
delete theAnaEmu; |
| 253 |
+ |
delete theAnaSynch; |
| 254 |
+ |
delete theAnaClu; |
| 255 |
+ |
delete theAnaTimingL1; |
| 256 |
+ |
delete theAnaMenu; |
| 257 |
+ |
delete theAnaEvent; |
| 258 |
+ |
delete theAnaDigiSpec; |
| 259 |
+ |
delete theAnaHitSpec; |
| 260 |
+ |
delete thePatternProducer; |
| 261 |
+ |
delete thePatternProvider; |
| 262 |
|
} |
