Changeset 260 in svn for trunk/Resolutions.cpp

Timestamp:

Feb 9, 2009, 3:37:19 PM (16 years ago)

Author:

severine ovyn

Message:

add header

File:

• trunk/Resolutions.cpp (modified) (4 diffs)

trunk/Resolutions.cpp

@@ -1 @@
-/*
-  ---- FastSim ----
-  A Fast Simulator for general purpose LHC detector
-  S. Ovyn ~~~~ severine.ovyn@uclouvain.be
-
-  Center for Particle Physics and Phenomenology (CP3)
-  Universite Catholique de Louvain (UCL)
-  Louvain-la-Neuve, Belgium
-*/
-
-/// \file Smearing.cpp
-/// \brief executable for the FastSim
+/***********************************************************************
+**                                                                    **
+**          /----------------------------------------------\          **
+**         |  Delphes, a framework for the fast simulation  |         **
+**         |       of a  generic collider experiment        |         **
+**          \----------------------------------------------/          **
+**                                                                    **
+**                                                                    **
+**   This package uses:                                               **
+**   ------------------                                               **
+**   FastJet algorithm: Phys. Lett. B641 (2006) [hep-ph/0512210]      ** 
+**   Hector: JINST 2:P09005 (2007) [physics.acc-ph:0707.1198v2]       ** 
+**   FROG: [hep-ex/0901.2718v1]                                       **
+**                                                                    **
+** ------------------------------------------------------------------ **
+**                                                                    **
+**   Main authors:                                                    **
+**   -------------                                                    **
+**                                                                    **
+**                Severine Ovyn                Xavier Rouby           ** 
+**          severine.ovyn@uclouvain.be      xavier.rouby@cern         ** 
+**                                                                    ** 
+**         Center for Particle Physics and Phenomenology (CP3)        ** 
+**                Universite catholique de Louvain (UCL)              **       
+**                     Louvain-la-Neuve, Belgium                      **            
+**                                                                    **
+**                      Copyright (C) 2008-2009,                      **
+**                        All rights reserved.                        **  
+**                                                                    **
+***********************************************************************/
 
 #include "TChain.h"
@@ -29 +47 @@
 #include "JetsUtil.h"
 #include "BFieldProp.h"
-
-//#include "PseudoJet.hh"
-//#include "ClusterSequence.hh"
 
 #include<vector>
@@ -367 +382 @@
       TRootJet *jet;
       itJet.Reset();
-//cout<<"a"<<endl;
       while( (jet = (TRootJet*) itJet.Next()) )
         {
-//cout<<"b"<<endl;
           TLorentzVector JETT(0,0,0,0);
           JETT.SetPxPyPzE(jet->Px,jet->Py,jet->Pz,jet->E);
           if(fabs(JETT.Eta()) < (DET->CEN_max_tracker - DET->TAU_track_scone)) 
            {
-//cout<<"c"<<endl;
             for(Int_t i=0; i<TausHadr.GetEntries();i++)
               {
-//cout<<"d"<<endl;
                 if(DeltaR(TausHadr[i]->Phi,TausHadr[i]->Eta,JETT.Phi(),JETT.Eta())<0.1)
                   {
-//cout<<"e"<<endl;
                     elementTaujet= (TAUHAD*) branchtaujet->NewEntry();
                     elementTaujet->EnergieCen = (EnergySmallCone(towers,JETT.Eta(),JETT.Phi(),DET->TAU_energy_scone,DET->JET_seed)/JETT.E());
@@ -387 +397 @@
                     if( (EnergySmallCone(towers,JETT.Eta(),JETT.Phi(),DET->TAU_energy_scone,DET->JET_seed)/JETT.E()) > 0.95 
                      &&  (NumTracks(branchTracks,DET->TAU_track_pt,JETT.Eta(),JETT.Phi(),DET->TAU_track_scone))==1)numTauRec++;
-//cout<<"f"<<endl;
                     
                  }
-//cout<<"g"<<endl;
               }
            }
-//cout<<"h"<<endl;
          
         
       } // for itJet : loop on all jets
-//cout<<"i"<<endl;
       
       treeWriter->Fill();
     } // Loop over all events
   treeWriter->Write();
-float frac = numTauRec/numTau;  
-cout<<numTauRec<<endl;
-cout<<numTau<<endl;
+  float frac = numTauRec/numTau;  
+  cout<<numTauRec<<endl;
+  cout<<numTau<<endl;
 
   cout << "** Exiting..." << endl;