Changeset 63 in svn for trunk/Delphes.cpp

Timestamp:

Dec 2, 2008, 9:50:08 AM (16 years ago)

Author:

severine ovyn

Message:

Trigger log file

File:

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

trunk/Delphes.cpp

@@ -29 +29 @@
 #include "interface/VeryForward.h"
 #include "interface/JetUtils.h"
+
+#include "Utilities/FROG/Examples/Sim_Delphes/FrogUtil.h"
 
 #include <vector>
@@ -93 +95 @@
   DET->Logfile(LogName);
 
+
   //Trigger information
-  Trigger *TRIG = new Trigger();
-  TRIG->TriggerReader("data/trigger.dat");
+  TriggerTable *TRIGT = new TriggerTable();
+  TRIGT->TriggerCardReader("data/trigger.dat");
+  TRIGT->PrintTriggerTable(LogName);
+  //outputFile->Close();
+  
 
   //Propagation of tracks in the B field
@@ -227 +233 @@
             { 
             genMomentum.SetPxPyPzE(particle->Px, particle->Py, particle->Pz, particle->E);
+//cout<<"avant "<<genMomentum.Pt()<<endl;
             TRACP->Propagation(particle,genMomentum);
+//cout<<"apres "<<genMomentum.Pt()<<endl;
             float eta=fabs(genMomentum.Eta());
             switch(pid) {
@@ -234 +242 @@
               DET->SmearElectron(genMomentum);
               electron.push_back(genMomentum);
+  //            cout<<"apres le smearing "<<genMomentum.Pt()<<endl;
+
               elecPID.push_back(particle->PID);
               break; // case pE
@@ -267 +277 @@
                 // back of the input_particles vector
                 input_particles.push_back(fastjet::PseudoJet(genMomentum.Px(),genMomentum.Py(),genMomentum.Pz(), genMomentum.E())); 
-                
+            //   cout<<"on remplie avec "<<pid<<" Pt "<<genMomentum.Pt()<<" Eta "<<genMomentum.Eta()<<" Phi "<<genMomentum.Eta()<<endl;
+  
                 genMomentumCalo.SetPxPyPzE(CaloTower.fourVector.px,CaloTower.fourVector.py,CaloTower.fourVector.pz,CaloTower.fourVector.E);
                 elementCalo = (TRootCalo*) branchCalo->NewEntry();
@@ -300 +311 @@
           
         } // while
-      
+// cout<<"*************"<<endl;     
       for(unsigned int i=0; i < electron.size(); i++) {
         if(electron[i].E()!=0 && fabs(electron[i].Eta()) < DET->MAX_TRACKER && electron[i].Pt() > DET->ELEC_pt)