Changeset 67 in svn for trunk/Delphes.cpp

Timestamp:

Dec 2, 2008, 12:14:45 PM (16 years ago)

Author:

severine ovyn

Message:

start trigger

File:

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

trunk/Delphes.cpp

@@ -30 +30 @@
 #include "interface/JetUtils.h"
 
+#include "Utilities/FROG/Examples/Sim_Delphes/FrogUtil.h"
 
 #include <vector>
@@ -353 +354 @@
   delete treeWriter;
 
-/*  if(DET->DOTRIGGER == 1) 
+  if(DET->DOTRIGGER == 1) 
     {
-//      TChain chainT("Analysis");
-//      chainT.Add(outputfilename.c_str());
-      ExRootTreeReader *treeReaderT = new ExRootTreeReader("Analysis");
-
+      TChain chainT("Analysis");
+      chainT.Add(outputfilename.c_str());
+      ExRootTreeReader *treeReaderT = new ExRootTreeReader(&chainT);
 
       const TClonesArray *branchElecTrig = treeReaderT->UseBranch("Electron");
@@ -366 +366 @@
       const TClonesArray *branchPhotonTrig = treeReaderT->UseBranch("Photon");
       const TClonesArray *branchETmisTrig = treeReaderT->UseBranch("ETmis");
-   
-     // Loop over all events
+  
+      ExRootTreeWriter *treeWriterT = new ExRootTreeWriter(outputfilename, "Trigger");
+      ExRootTreeBranch *branchTrigger = treeWriterT->NewBranch("Trigger", TRootTrigger::Class());
+
+      TRootTrigger *elementTrigger;
+
      Long64_t entryT, allEntriesT = treeReaderT->GetEntries();
-      cout << "** Chain contains " << allEntriesT << " events" << endl;
+     cout << "** Chain contains " << allEntriesT << " events" << endl;
      for(entryT = 0; entryT < allEntriesT; ++entryT)
         {
-//cout<<"entry "<<entryT<<endl;
+          treeWriterT->Clear();
+//          treeReaderT->ReadEntry(entry);
+          elementTrigger = (TRootTrigger*) branchTrigger->NewEntry();
+          elementTrigger->Accepted=true;
+          treeWriterT->Fill();
         }
+
+    treeWriterT->Write();
+    delete treeWriterT;
    }
-*/
 
   cout << "** Exiting..." << endl;