Changeset 43 in svn for trunk/Resolutions.cpp

Timestamp:

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

Author:

severine ovyn

Message:

general jet variable

File:

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

trunk/Resolutions.cpp

@@ -224 +224 @@
   // set up a CDF midpoint jet definition
 #ifdef ENABLE_PLUGIN_CDFCONES
-  plugins = new fastjet::CDFJetCluPlugin(0,DET->CONERADIUS,DET->C_ADJACENCYCUT,DET->C_MAXITERATIONS,DET->C_IRATCH,DET->C_OVERLAPTHRESHOLD);
+  plugins = new fastjet::CDFJetCluPlugin(0,DET->CONERADIUS,DET->C_ADJACENCYCUT,DET->C_MAXITERATIONS,DET->C_IRATCH,DET->OVERLAPTHRESHOLD);
   jet_def = fastjet::JetDefinition(plugins);
 #else
@@ -232 +232 @@
   // set up a CDF midpoint jet definition
 #ifdef ENABLE_PLUGIN_CDFCONES
-  pluginsS = new fastjet::CDFJetCluPlugin(1,DET->CONERADIUS,DET->C_ADJACENCYCUT,DET->C_MAXITERATIONS,DET->C_IRATCH,DET->C_OVERLAPTHRESHOLD);
+  pluginsS = new fastjet::CDFJetCluPlugin(1,DET->CONERADIUS,DET->C_ADJACENCYCUT,DET->C_MAXITERATIONS,DET->C_IRATCH,DET->OVERLAPTHRESHOLD);
   jet_defS = fastjet::JetDefinition(pluginsS);
 #else
@@ -240 +240 @@
  // set up a CDF midpoint jet definition
  #ifdef ENABLE_PLUGIN_CDFCONES
- pluginsT = new fastjet::CDFJetCluPlugin(0,DET->CONERADIUS,DET->C_ADJACENCYCUT,DET->C_MAXITERATIONS,DET->C_IRATCH,DET->C_OVERLAPTHRESHOLD);
+ pluginsT = new fastjet::CDFJetCluPlugin(0,DET->CONERADIUS,DET->C_ADJACENCYCUT,DET->C_MAXITERATIONS,DET->C_IRATCH,DET->OVERLAPTHRESHOLD);
  jet_defT = fastjet::JetDefinition(pluginsT);
  #else
@@ -256 +256 @@
       treeReader->ReadEntry(entry);
       treeWriter->Clear();
-      
       if((entry % 100) == 0 && entry > 0 )  cout << "** Processing element # " << entry << endl;
       
@@ -329 +328 @@
               } // switch (pid)
               
-              if(pid != pMU)
+              if(pid != pMU && genMomentum.Et() !=0)
                 {
                   towers.push_back(PhysicsTower(LorentzVector(genMomentum.Px(),genMomentum.Py(),genMomentum.Pz(), genMomentum.E())));
@@ -354 +353 @@
     
      TLorentzVector Att(0.,0.,0.,0.);
-     float ScalarEt=0;
      for(unsigned int i=0; i < towers.size(); i++) 
        {
          Att.SetPxPyPzE(towers[i].fourVector.px,towers[i].fourVector.py,towers[i].fourVector.pz,towers[i].fourVector.E);
-         ScalarEt = ScalarEt + Att.Et();
          PTmisS = PTmisS + Att;
        }
-//cout<<"non smeare "<<PTmis.Pt()<<" "<<PTmisS.Pt()<<endl;
-//cout<<"smeare "<<PTmis.Px()<<" "<<PTmisS.Px()<<endl;
-//cout<<"**********"<<endl;
 
       elementEtmis= (ETMIS*) branchetmis->NewEntry();
       elementEtmis->Et = (PTmis).Pt();
-      elementEtmis->Ex = (-PTmis).Px();
-      elementEtmis->SEt = ScalarEt;
-
       elementEtmis->EtSmeare = (PTmisS).Pt()-(PTmis).Pt();
-      elementEtmis->ExSmeare = (-PTmisS).Px()-(PTmis).Px();
       
       //*****************************