ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/MitAna/TreeMod/src/HLTMod.cc
(Generate patch)

Comparing UserCode/MitAna/TreeMod/src/HLTMod.cc (file contents):
Revision 1.11 by loizides, Mon Jun 15 15:00:17 2009 UTC vs.
Revision 1.20 by paus, Mon Mar 21 15:58:37 2011 UTC

# Line 8 | Line 8
8   #include "MitAna/DataTree/interface/TriggerMask.h"
9   #include "MitAna/DataTree/interface/TriggerObject.h"
10   #include "MitAna/DataTree/interface/TriggerObjectCol.h"
11 + #include "MitAna/DataTree/interface/TriggerObjectsTable.h"
12  
13   using namespace mithep;
14  
# Line 18 | Line 19 | HLTMod::HLTMod(const char *name, const c
19    BaseMod(name,title),
20    fAbort(kTRUE),
21    fPrintTable(kFALSE),
22 +  fIgnoreBits(kFALSE),
23 +  fObjMode(kHlt),
24    fBitsName(Names::gkHltBitBrn),
25    fMyObjsNamePub(Form("%sTrigObjs", name)),
26    fBits(0),
# Line 25 | Line 28 | HLTMod::HLTMod(const char *name, const c
28    fTriggers(0),
29    fTrigObjs(0),
30    fNEvents(0),
31 <  fNAcceped(0),
31 >  fNAccepted(0),
32    fNFailed(0)
33   {
34    // Constructor.
# Line 38 | Line 41 | HLTMod::~HLTMod()
41   }
42  
43   //--------------------------------------------------------------------------------------------------
44 < void HLTMod::AddTrigger(const char *expr)
44 > void HLTMod::AddTrigger(const char *expr, UInt_t firstRun, UInt_t lastRun)
45   {
46    // Add trigger search pattern to the list of patters. Each element of the list is logically
47    // "ored". The given expression can contain several trigger names logically "anded" (using "&").
# Line 46 | Line 49 | void HLTMod::AddTrigger(const char *expr
49    // "A", "!A", "A&B", "A&!B" or "A&B&C"  
50  
51    string tname(expr);
52 <  fTrigNames.push_back(tname);
52 >  std::pair<std::string,std::pair<UInt_t,UInt_t> > element(tname,std::pair<UInt_t,UInt_t>(firstRun,lastRun));
53 >  fTrigNames.push_back(element);
54   }
55  
56   //--------------------------------------------------------------------------------------------------
# Line 57 | Line 61 | void HLTMod::AddTrigObjs(UInt_t tid)
61    const BitMask256 &ba = fTrigBitsAnd.at(tid);
62    const BitMask256 &bm = fTrigBitsCmp.at(tid);
63    for (UInt_t i=0; i<bm.Size(); ++i) {
64 <    if (ba.TestBit(i)==0)
64 >    if (ba.TestBit(i)==0 && !fIgnoreBits)
65        continue; // not an active trigger bit
66      if (fBitsDone.TestBit(i))
67        continue; // objects for this bit are already obtained
# Line 65 | Line 69 | void HLTMod::AddTrigObjs(UInt_t tid)
69        continue; // excluded trigger bit (ie a !trgname)
70  
71      const TList *list = fTrigObjs->GetList(i);
72 <    fMyTrgObjs->Add(list);    
72 >    if (list) {
73 >      TIter iter(list->MakeIterator());
74 >      const TriggerObject *to = dynamic_cast<const TriggerObject*>(iter.Next());
75 >      while (to) {
76 >        if ( (fObjMode == kAll) ||
77 >             ((fObjMode==kHlt) && (to->IsHLT())) ||
78 >             ((fObjMode==kL1) && (to->IsL1())) )
79 >          fMyTrgObjs->Add(to);    
80 >        to = dynamic_cast<const TriggerObject*>(iter.Next());
81 >      }
82 >    }
83      fBitsDone.SetBit(i);
84    }
85   }
# Line 78 | Line 92 | void HLTMod::BeginRun()
92    fTrigBitsAnd.clear();
93    fTrigBitsCmp.clear();
94  
95 <  if (fPrintTable)
95 >  if (fPrintTable) {
96 >    Info("BeginRun", "Get trigger table for run %d", GetEventHeader()->RunNum());
97      fTriggers->Print();
98 +  }
99 +  
100 +  UInt_t runNumber = GetEventHeader()->RunNum();
101  
102    for (UInt_t i=0; i<fTrigNames.size(); ++i) {
103 +    
104 +    UInt_t firstRun = fTrigNames.at(i).second.first;
105 +    UInt_t lastRun  = fTrigNames.at(i).second.second;
106 +    
107 +    // implement run dependence of the triggers
108 +    if ( (!(firstRun==0 && lastRun==0)) &&
109 +         ( runNumber<firstRun || runNumber>lastRun ) )
110 +      continue;
111 +    
112      BitMask256 tmask; //trigger mask
113      BitMask256 amask; //bitand mask
114 <    TString names(fTrigNames.at(i).c_str());
114 >    TString names(fTrigNames.at(i).first.c_str());
115  
116      TObjArray *arr = names.Tokenize("&");
117      if (arr) {
# Line 103 | Line 130 | void HLTMod::BeginRun()
130            Warning("BeginRun", "Trigger %s not found.", sptr);
131            continue;
132          }
133 +
134          UShort_t bit = tn->Id();
135 <        amask.SetBit(bit); //always set and-mask bit
136 <        if (!invert)
137 <          tmask.SetBit(bit); //set trigger bit
135 >        if (amask.TestBit(bit)) {
136 >          if (tmask.TestBit(bit)==invert) {
137 >            amask.ClearBit(bit);
138 >            tmask.ClearBit(bit);
139 >            Warning("BeginRun", "Trigger expression %s always false.", names.Data());
140 >            break;
141 >          }
142 >        } else { //always set and-mask bit
143 >          amask.SetBit(bit);
144 >          if (!invert)
145 >            tmask.SetBit(bit); //set trigger bit
146 >        }
147        }
148        delete arr;
149      }
150 <    fTrigBitsAnd.push_back(amask);
151 <    fTrigBitsCmp.push_back(tmask);
150 >    if (amask.NBitsSet()) {
151 >      fTrigBitsAnd.push_back(amask);
152 >      fTrigBitsCmp.push_back(tmask);
153 >    }
154    }
155   }
156  
# Line 151 | Line 190 | void HLTMod::Process()
190    }
191  
192    // take action if accepted
193 <  ++fNAcceped;
193 >  ++fNAccepted;
194    IncNEventsProcessed();
195    OnAccepted();
196    if (!AddObjThisEvt(fMyTrgObjs)) {
# Line 187 | Line 226 | void HLTMod::SlaveTerminate()
226   {
227    // Save number of accepted events.
228  
229 +  printf(" %s - Accepted events: %d (/%d)\n",GetName(),fNAccepted,fNEvents);
230    SaveNEventsProcessed("hDHLTEvents");
231   }

Diff Legend

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