Changeset 43 in svn for trunk/Delphes.cpp

Timestamp:

Nov 19, 2008, 2:06:15 PM (16 years ago)

Author:

severine ovyn

Message:

general jet variable

File:

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

trunk/Delphes.cpp

@@ -168 +168 @@
   vector<TLorentzVector> TrackCentral;  
   vector<PhysicsTower> towers;
-  
   vector<fastjet::PseudoJet> input_particles;//for FastJet algorithm
   vector<fastjet::PseudoJet> inclusive_jets;
@@ -214 +213 @@
     // set up a CDF midpoint jet definition
 #ifdef ENABLE_PLUGIN_CDFCONES
-    plugins = new fastjet::CDFJetCluPlugin(DET->C_SEEDTHRESHOLD,DET->CONERADIUS,DET->C_ADJACENCYCUT,DET->C_MAXITERATIONS,DET->C_IRATCH,DET->C_OVERLAPTHRESHOLD);
+    plugins = new fastjet::CDFJetCluPlugin(DET->SEEDTHRESHOLD,DET->CONERADIUS,DET->C_ADJACENCYCUT,DET->C_MAXITERATIONS,DET->C_IRATCH,DET->OVERLAPTHRESHOLD);
     jet_def = fastjet::JetDefinition(plugins);
 #else
@@ -226 +225 @@
     // set up a CDF midpoint jet definition
 #ifdef ENABLE_PLUGIN_CDFCONES
-    plugins = new fastjet::CDFMidPointPlugin (DET->M_SEEDTHRESHOLD,DET->CONERADIUS,DET->M_CONEAREAFRACTION,DET->M_MAXPAIRSIZE,DET->M_MAXPAIRSIZE,DET->C_OVERLAPTHRESHOLD);
+    plugins = new fastjet::CDFMidPointPlugin (DET->SEEDTHRESHOLD,DET->CONERADIUS,DET->M_CONEAREAFRACTION,DET->M_MAXPAIRSIZE,DET->M_MAXPAIRSIZE,DET->OVERLAPTHRESHOLD);
     jet_def = fastjet::JetDefinition(plugins);
 #else
@@ -238 +237 @@
     int npass = 0;               // do infinite number of passes
     double protojet_ptmin = 0.0; // use all protojets
-    plugins = new fastjet::SISConePlugin (DET->CONERADIUS,DET->C_OVERLAPTHRESHOLD,npass, protojet_ptmin);
+    plugins = new fastjet::SISConePlugin (DET->CONERADIUS,DET->OVERLAPTHRESHOLD,npass, protojet_ptmin);
     jet_def = fastjet::JetDefinition(plugins);
 #else
@@ -262 +261 @@
       treeReader->ReadEntry(entry);
       treeWriter->Clear();
-      
       if((entry % 100) == 0 && entry > 0 )  cout << "** Processing element # " << entry << endl;
       
@@ -337 +335 @@
             // all final particles but muons and neutrinos      
             // for calorimetric towers and mission PT
-            if(genMomentum.E()!=0) {
+            
+            if(genMomentum.E() !=0) {
               if(pid !=pMU) {
                 PhysicsTower CaloTower = PhysicsTower(LorentzVector(genMomentum.Px(),genMomentum.Py(),genMomentum.Pz(), genMomentum.E()));
@@ -371 +370 @@
                 TrackCentral.push_back(genMomentum);
               }
-          } // switch
+          
+} // switch
           
           // Forward particles in CASTOR ?
@@ -470 +470 @@
       }
       
-      
       // computes the Missing Transverse Momentum
       TLorentzVector Att(0.,0.,0.,0.);
@@ -483 +482 @@
       elementEtmis->Px = (-PTmis).Px();
       elementEtmis->Py = (-PTmis).Py();
-      
       //*****************************
       
       // run the jet clustering with the above jet definition
+      cout<<input_particles.size()<<endl;
       if(input_particles.size()!=0)
         {
           fastjet::ClusterSequence clust_seq(input_particles, jet_def);
-          
-          
           // extract the inclusive jets with pt > 5 GeV
           double ptmin = 5.0;
           inclusive_jets = clust_seq.inclusive_jets(ptmin);
-          
+cout<<"inclusive_jets "<<inclusive_jets.size()<<endl;     
           // sort jets into increasing pt
           sorted_jets = sorted_by_pt(inclusive_jets);
         }
+cout<<"sorted_jets.size() "<<sorted_jets.size()<<endl;
       for (unsigned int i = 0; i < sorted_jets.size(); i++) {
         TLorentzVector JET;