Changeset 72 in svn for trunk/Delphes.cpp

Timestamp:

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

Author:

severine ovyn

Message:

final trigger

File:

• trunk/Delphes.cpp (modified) (3 diffs)

trunk/Delphes.cpp

@@ -95 +95 @@
 
   //Trigger information
-  //Trigger *TRIG = new Trigger();
-  //TRIG->TriggerReader("data/trigger.dat");
+  TriggerTable *TRIGT = new TriggerTable();
+  TRIGT->TriggerCardReader("data/trigger.dat");
+  TRIGT->PrintTriggerTable(LogName);
 
   //Propagation of tracks in the B field
@@ -336 +337 @@
 
       //*****************************
-      treeWriter->Fill();
       
       sorted_jets=JETRUN->RunJets(input_particles);
@@ -344 +344 @@
       // Add here the trigger
       // Should test all the trigger table on the event, based on reconstructed objects
-     
+      treeWriter->Fill();
+    
     } // Loop over all events
 
   treeWriter->Write();
-  
+  delete treeWriter;
+
+  if(DET->DOTRIGGER == 1)
+    {
+      TChain chainT("Analysis");
+      chainT.Add(outputfilename.c_str());
+      ExRootTreeReader *treeReaderT = new ExRootTreeReader(&chainT);
+
+      TClonesArray *branchElecTrig = treeReaderT->UseBranch("Electron");
+      TClonesArray *branchMuonTrig = treeReaderT->UseBranch("Muon");
+      TClonesArray *branchJetTrig = treeReaderT->UseBranch("Jet");
+      TClonesArray *branchTauJetTrig = treeReaderT->UseBranch("TauJet");
+      TClonesArray *branchPhotonTrig = treeReaderT->UseBranch("Photon");
+      TClonesArray *branchETmisTrig = treeReaderT->UseBranch("ETmis");
+
+      ExRootTreeWriter *treeWriterT = new ExRootTreeWriter(outputfilename, "Trigger");
+      ExRootTreeBranch *branchTrigger = treeWriterT->NewBranch("TrigResult", TRootTrigger::Class());
+
+      Long64_t entryT, allEntriesT = treeReaderT->GetEntries();
+      cout << "** Chain contains " << allEntriesT << " events" << endl;
+      for(entryT = 0; entryT < allEntriesT; ++entryT)
+        {
+          treeWriterT->Clear();
+          treeReaderT->ReadEntry(entryT);
+          TRIGT->GetGlobalResult(branchElecTrig, branchMuonTrig,branchJetTrig, branchTauJetTrig,branchPhotonTrig, branchETmisTrig,branchTrigger);
+          treeWriterT->Fill();
+        }
+
+      treeWriterT->Write();
+      delete treeWriterT;
+    }
+
   cout << "** Exiting..." << endl;
   
-  delete treeWriter;
   delete treeReader;
   delete DET;