ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/L1RpcTriggerAnalysis/plugins/L1RpcTreeAnalysis.cc
(Generate patch)

Comparing UserCode/L1RpcTriggerAnalysis/plugins/L1RpcTreeAnalysis.cc (file contents):
Revision 1.3 by konec, Mon Oct 8 07:09:24 2012 UTC vs.
Revision 1.12 by konec, Sun Dec 16 22:56:15 2012 UTC

# Line 16 | Line 16
16  
17   #include "UserCode/L1RpcTriggerAnalysis/interface/L1Obj.h"
18   #include "UserCode/L1RpcTriggerAnalysis/interface/EventObj.h"
19 + #include "UserCode/L1RpcTriggerAnalysis/interface/EventObjBXExtra.h"
20   #include "UserCode/L1RpcTriggerAnalysis/interface/TrackObj.h"
21   #include "UserCode/L1RpcTriggerAnalysis/interface/MuonObj.h"
22   #include "UserCode/L1RpcTriggerAnalysis/interface/L1ObjColl.h"
23   #include "UserCode/L1RpcTriggerAnalysis/interface/SynchroCountsObj.h"
24 + #include "UserCode/L1RpcTriggerAnalysis/interface/DetCluDigiObj.h"
25 + #include "UserCode/L1RpcTriggerAnalysis/interface/TriggerMenuResultObj.h"
26 +
27  
28   #include "UserCode/L1RpcTriggerAnalysis/interface/ConverterRPCRawSynchroSynchroCountsObj.h"
29  
30 + L1RpcTreeAnalysis::L1RpcTreeAnalysis(const edm::ParameterSet & cfg)
31 +  : theConfig(cfg),
32 +    theAnaMuonDistribution( cfg.getParameter<edm::ParameterSet>("anaMuonDistribution") ),
33 +    theAnaTimingL1 (cfg.getParameter<edm::ParameterSet>("anaTimingL1") ),
34 +    theAnaMenu(cfg.getParameter<edm::ParameterSet>("anaMenu") ),
35 +    theAnaEvent(cfg.getParameter<edm::ParameterSet>("anaEvent") )
36 + { }
37  
38   void L1RpcTreeAnalysis::beginJob()
39   {
# Line 36 | Line 47 | void L1RpcTreeAnalysis::beginJob()
47    theAnaEmu.init(theHistos);
48    theAnaSynch.init(theHistos);
49    theAnaClu.init(theHistos);
50 <
50 >  theAnaTimingL1.init(theHistos);
51 >  theAnaEvent.init(theHistos);
52 >  theAnaMenu.init(theHistos);
53   }
54  
55   void L1RpcTreeAnalysis::beginRun(const edm::Run& ru, const edm::EventSetup& es)
# Line 70 | Line 83 | void L1RpcTreeAnalysis::analyze(const ed
83    std::vector<SynchroCountsObj> *counts= 0;
84    std::vector<uint32_t> *detsCrossedByMuon = 0;
85    std::vector<uint32_t> *detsCrossedByMuonDeepInside = 0;
86 <  std::vector<uint32_t> *detsHitsCompatibleWithMuon = 0;
86 >  std::vector<DetCluDigiObj> *detsHitsCompatibleWithMuon = 0;
87    std::vector<uint32_t> *detsSIMU = 0;
75 –  std::vector<uint32_t> *nDigisCompDets = 0;
76 –  std::vector<uint32_t> *clSizeCompDets = 0;
88  
89    EventObj * event = 0;
90    MuonObj * muon = 0;
91  
92 +  TriggerMenuResultObj *bitsL1  = 0;
93 +  TriggerMenuResultObj *bitsHLT = 0;
94 +
95    TBranch *bcounts=0;
96    TBranch *bdetsCrossedByMuon =0;
97    TBranch *bdetsCrossedByMuonDeepInside =0;
98    TBranch *bdetsHitsCompatibleWithMuon = 0;
99    TBranch *bdetsSIMU =0;
86 –  TBranch *bnDigisCompDets =0;
87 –  TBranch *bclSizeCompDets =0;
100  
101 <  L1ObjColl* l1RpcColl = 0;
90 <  L1ObjColl* l1OtherColl = 0;
91 <  L1ObjColl* l1RpcCollEmu = 0;
101 >  L1ObjColl* l1ObjColl = 0;
102  
103    chain.SetBranchAddress("event",&event);
104    chain.SetBranchAddress("muon",&muon);
105  
106 +  chain.SetBranchAddress("bitsL1",&bitsL1);
107 +  chain.SetBranchAddress("bitsHLT",&bitsHLT);
108 +  
109    chain.SetBranchAddress("counts",&counts,&bcounts);
110    chain.SetBranchAddress("detsCrossedByMuon",&detsCrossedByMuon,&bdetsCrossedByMuon);
111    chain.SetBranchAddress("detsCrossedByMuonDeepInside",&detsCrossedByMuonDeepInside,&bdetsCrossedByMuonDeepInside);
112    chain.SetBranchAddress("detsHitsCompatibleWithMuon",&detsHitsCompatibleWithMuon,&bdetsHitsCompatibleWithMuon);
113    chain.SetBranchAddress("detsSIMU",&detsSIMU,&bdetsSIMU);
101 –  chain.SetBranchAddress("nDigisCompDets",&nDigisCompDets, &bnDigisCompDets);
102 –  chain.SetBranchAddress("clSizeCompDets",&clSizeCompDets, &bclSizeCompDets);
114  
115 <  chain.SetBranchAddress("l1RpcColl",&l1RpcColl);
105 <  chain.SetBranchAddress("l1OtherColl",&l1OtherColl);
106 <  chain.SetBranchAddress("l1RpcCollEmu",&l1RpcCollEmu);
115 >  chain.SetBranchAddress("l1ObjColl",&l1ObjColl);
116  
117  
118    //
119 <  // number of efents
119 >  // number of events
120    //
121    Int_t nentries = (Int_t) chain.GetEntries();
122    std::cout <<" ENTRIES: " << nentries << std::endl;
# Line 119 | Line 128 | void L1RpcTreeAnalysis::analyze(const ed
128    unsigned int lastRun = 0;
129    for (int ev=0; ev<nentries; ev++) {
130      chain.GetEntry(ev);
131 +    theAnaMenu.updateMenu(bitsL1->names, bitsHLT->names);
132  
133      if (lastRun != (*event).run) {
134        lastRun = (*event).run;
135        std::cout <<"RUN:"    << std::setw(7) << (*event).run
136                  <<" event:" << std::setw(8) << ev
137 <                <<" done:"  << std::setw(6)<< std::setiosflags(std::ios::fixed) << std::setprecision(2) << ev*100./nentries<<"%"<<std::endl;
137 >                <<" done:"  << std::setw(6)<< std::setiosflags(std::ios::fixed) << std::setprecision(2) << ev*100./nentries<<"%";
138 >      std::cout<<std::endl;
139      }
140  
130 – //    if (ev > 10) break;
141  
142 < //    if (event->run != 178854) continue;
143 < //    if (lastLumi != (*event).lumi) { lastLumi = (*event).lumi; std::cout <<"lumi: " << (*event).lumi<<std::endl; }
144 < //   theAnaDet.debug = false;
145 < //    if ((*event).id==60422922)theAnaRpcMisc.debug = true;
146 <
147 <   theAnaMuonDistribution.run(muon);
148 <   theAnaRpcVsOth.run(muon,l1RpcColl,l1OtherColl);
149 <   theAnaEff.run(muon,l1RpcColl,l1OtherColl);
150 <   theAnaRpcMisc.run(event,muon,l1RpcColl,l1OtherColl);
151 <   theAnaDet.run( muon, *detsHitsCompatibleWithMuon,  *detsCrossedByMuon, *detsCrossedByMuonDeepInside);
152 <   theAnaEmu.run ( event, muon, l1RpcCollEmu, l1RpcColl);
153 <   theAnaSynch.run( event, muon, ConverterRPCRawSynchroSynchroCountsObj::toRawSynchro( *counts));
154 <   theAnaClu.run( event, muon, l1RpcColl, *detsHitsCompatibleWithMuon, *nDigisCompDets, *clSizeCompDets);
155 <
156 < //   theAnaEmu.debug =theAnaDet.debug;
157 < //    std::cout <<"----------"<<std::endl;
158 < //   theAnaDet.debug =theAnaEmu.debug;
159 < //    if (anaDet.debug) std::cout <<" Event: "<<(*event).id <<" Lumi: "<<(*event).lumi<< std::endl;
160 <  
142 >   // EVENT NUMBER, BX structure etc.
143 >   EventObjBXExtra eventBx(*event);
144 >   if ( !theAnaEvent.filter(&eventBx) && theConfig.getParameter<bool>("filterByAnaEvent") ) continue;
145 >   // ANALYSE AND FILTER KINEMCTICS
146 >   if ( !theAnaMuonDistribution.filter(muon) && theConfig.getParameter<bool>("filterByAnaMuonDistribution") ) continue;
147 >   // ANALYSE AND FILTER TRIGGER MENU
148 >   if ( !theAnaMenu.filter(event, muon, bitsL1, bitsHLT) && theConfig.getParameter<bool>("filterByAnaMenu") ) continue;
149 >
150 > /*
151 > //   theAnaRpcVsOth.run(muon,l1ObjColl);
152 > //   theAnaEff.run(muon, l1ObjColl);
153 > //   theAnaRpcMisc.run(event,muon,l1ObjColl);
154 > //   theAnaDet.run( muon, *detsHitsCompatibleWithMuon,  *detsCrossedByMuon, *detsCrossedByMuonDeepInside);
155 > //   theAnaEmu.run ( event, muon, l1ObjColl);
156 > //   theAnaSynch.run( event, muon, ConverterRPCRawSynchroSynchroCountsObj::toRawSynchro( *counts));
157 > //   theAnaClu.run( event, muon, l1ObjColl, *detsHitsCompatibleWithMuon);
158 > */
159 >   theAnaTimingL1.run( &eventBx, muon, l1ObjColl);
160 >
161    }
162   }
163  
# Line 157 | Line 167 | void L1RpcTreeAnalysis::endJob()
167    TGraph* hGraph_DetEff = theAnaDet.resume();
168    theAnaRpcMisc.resume(theHistos);
169    TGraph* hGraph_RunClu = theAnaClu.resume();
170 +  theAnaTimingL1.resume(theHistos);
171 +  theAnaMenu.resume(theHistos);
172 +  theAnaEvent.resume(theHistos);
173  
174    theAnaSynch.endJob();
175  

Diff Legend

– Removed lines
+ Added lines
< Changed lines
> Changed lines