Changeset 72 in svn for trunk/src

Timestamp:

Dec 3, 2008, 12:01:04 PM (16 years ago)

Author:

severine ovyn

Message:

final trigger

Location:

trunk/src

Files:

• SmearUtil.cc (modified) (2 diffs)
• TriggerUtil.cc (modified) (1 diff)

trunk/src/SmearUtil.cc

@@ -86 +86 @@
 MISTAGGING_L =    1;
 
+//Trigger flag
+DOTRIGGER    =     1;
 
 CONERADIUS   = 0.7;     // generic jet radius ; not for tau's !!!
@@ -178 +180 @@
     else if(strstr(temp_string.c_str(),"TAUJET_pt")){curstring >> varname >> value; TAUJET_pt = value;}
     else if(strstr(temp_string.c_str(),"NTOWERS")){curstring >> varname >> ivalue; NTOWERS = ivalue;}
+    else if(strstr(temp_string.c_str(),"DOTRIGGER")){curstring >> varname >> ivalue; DOTRIGGER = ivalue;}
     else if(strstr(temp_string.c_str(),"TOWER_ETA_EDGES")){
         curstring >> varname;   for(unsigned int i=0; i<NTOWERS+1; i++) {curstring >> value; TOWER_ETA_EDGES[i] = value;} }

trunk/src/TriggerUtil.cc

@@ -65 +65 @@
 }
 
+bool TriggerTable::GetGlobalResult(TClonesArray *branchElecTrig, TClonesArray *branchMuonTrig,
+      TClonesArray *branchJetTrig, TClonesArray *branchTauJetTrig,
+      TClonesArray *branchPhotonTrig,  TClonesArray *branchETmisTrig, ExRootTreeBranch *branchTrigger) {
+
+  TRootTrigger *elementTrigger;
+  GlobalResult=false;
+  for(unsigned int i=0; i < list_of_trigger_bits.size(); i++) {
+  if(list_of_trigger_bits[i].GetResult(branchElecTrig, branchMuonTrig,  branchJetTrig, branchTauJetTrig, branchPhotonTrig, branchETmisTrig))
+    {
+       GlobalResult=true;
+       int flag=i+1;
+       elementTrigger = (TRootTrigger*) branchTrigger->NewEntry();
+       elementTrigger->Accepted=flag;
+    }
+  }
+  if(GlobalResult)
+    {
+      elementTrigger = (TRootTrigger*) branchTrigger->NewEntry();
+      elementTrigger->Accepted=0;
+    }
+    return GlobalResult;
+
+}
+
 //*************************************************************************
-
 TriggerBit::TriggerBit() {
- has_been_evaluated=false;
+}
+
+bool TriggerBit::GetResult(TClonesArray *branchElecTrig, TClonesArray *branchMuonTrig,
+      TClonesArray *branchJetTrig, TClonesArray *branchTauJetTrig,
+      TClonesArray *branchPhotonTrig,  TClonesArray *branchETmisTrig)
+{
+  int elec_size = ElecValues.size();
+  int muon_size = MuonValues.size();
+  int jet_size = JetValues.size();
+  int taujet_size = TauJetValues.size();
+  int gamma_size = GammaValues.size();
+  int etmis_size = EtmisValues.size();
+
+  Result=false;
+  if(
+     (branchElecTrig->GetEntries() >= elec_size )
+      && (branchMuonTrig->GetEntries() >= muon_size )
+      && (branchJetTrig->GetEntries() >= jet_size )
+      && (branchTauJetTrig->GetEntries() >= taujet_size )
+      && (branchPhotonTrig->GetEntries() >= gamma_size )
+    )
+    {
+      Result=true;
+      if(elec_size!=0){
+        TRootElectron *electron;
+        for(int i=0;i<elec_size;i++){
+           electron = (TRootElectron*)branchElecTrig->At(i);
+           if(electron->PT < ElecValues[i])Result=false;}}
+
+      if(muon_size!=0){
+        TRootMuon *muon;
+        for(int i=0;i<muon_size;i++){
+           muon = (TRootMuon*)branchMuonTrig->At(i);
+           if(muon->PT < MuonValues[i])Result=false;}}
+
+      if(jet_size!=0){
+        TRootJet *jet;
+        for(int i=0;i<jet_size;i++){
+           jet = (TRootJet*)branchJetTrig->At(i);
+           if(jet->PT < JetValues[i])Result=false;}}
+
+      if(taujet_size!=0){
+        TRootTauJet *taujet;
+        for(int i=0;i<taujet_size;i++){
+           taujet = (TRootTauJet*)branchTauJetTrig->At(i);
+           if(taujet->PT < TauJetValues[i])Result=false;}}
+
+      if(gamma_size!=0){
+        TRootPhoton *gamma;
+        for(int i=0;i<gamma_size;i++){
+           gamma = (TRootPhoton*)branchPhotonTrig->At(i);
+           if(gamma->PT < GammaValues[i])Result=false;}}
+
+
+      if(etmis_size!=0){
+        TRootETmis *etmis = (TRootETmis*)branchETmisTrig->At(0);
+           if(etmis->ET < EtmisValues[0])Result=false;}
+    
+   }
+  return Result;
 }