Changeset 23 in svn for trunk/Resolutions.cpp

Timestamp:

Nov 7, 2008, 6:37:38 PM (16 years ago)

Author:

severine ovyn

Message:

Working on the resolution

File:

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

trunk/Resolutions.cpp

@@ -172 +172 @@
   vector<fastjet::PseudoJet> inclusive_jetsS;
   vector<fastjet::PseudoJet> sorted_jetsS;
+  vector<PhysicsTower> towers;
 
   fastjet::JetDefinition jet_def;
@@ -180 +181 @@
   // 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(0,DET->CONERADIUS,DET->C_ADJACENCYCUT,DET->C_MAXITERATIONS,DET->C_IRATCH,DET->C_OVERLAPTHRESHOLD);
   jet_def = fastjet::JetDefinition(plugins);
   #else
@@ -200 +201 @@
   for(entry = 0; entry < allEntries; ++entry)
     {
+      TLorentzVector PTmisS(0,0,0,0);
       TLorentzVector PTmis(0,0,0,0);
       treeReader->ReadEntry(entry);
@@ -216 +218 @@
       inclusive_jetsS.clear();
       sorted_jetsS.clear();
+      towers.clear();
 
       // Loop over all particles in event
@@ -238 +241 @@
               )
             ) {
+                if(pid != pMU)PTmis = PTmis + genMomentum;//ptmis
                 switch(pid) {
                     
@@ -264 +268 @@
                 // all final particles but muons and neutrinos  
                 // for calorimetric towers and mission PT
-                if(genMomentum.E()!=0) PTmis = PTmis + genMomentum;//ptmis
-
-                if(pid != pMU && genMomentum.Pt() > DET->PT_TRACKS_MIN)
+                //if(genMomentum.E()!=0) PTmis = PTmis + genMomentum;//ptmis
+
+                if(pid != pMU)
                   {
+                     PTmisS = PTmisS + genMomentum;
+                     towers.push_back(PhysicsTower(LorentzVector(genMomentum.Px(),genMomentum.Py(),genMomentum.Pz(), genMomentum.E())));
                      input_particlesS.push_back(fastjet::PseudoJet(genMomentum.Px(),genMomentum.Py(),genMomentum.Pz(), genMomentum.E()));
                   }
@@ -289 +295 @@
         } // while
 
+  TLorentzVector toWerS(0,0,0,0);
+  //calcul de ETmis au depart des tours calorimetriques.. 
+/*  for(unsigned int i=0; i < towers.size(); i++) {
+      if(towers[i].fourVector.pt() < 0.5) continue;
+          toWerS.SetPxPyPzE(towers[i].fourVector.px,towers[i].fourVector.py,towers[i].fourVector.pz,towers[i].fourVector.E);
+          PTmisS = PTmisS + toWerS;
+          }
+*/
+//PTmis = (0,0,0,0)-PTmis;
+//PTmisS = (0,0,0,0)-PTmisS;
+
+//float ET=PTmis.Et()<<endl;
+//float ETS=PTmisS.Et()<<endl;
+
+cout<<"Ptmis "<<PTmis.Et()<<" PTmis smeare "<<PTmisS.Et()<<endl;
+cout<<"Ptmis "<<(-PTmis).Pt()<<" PTmis smeare "<<(-PTmisS).Pt()<<endl;
+
+  elementEtmis= (ETMIS*) branchetmis->NewEntry();
+  elementEtmis->Et = (PTmis).Et();
+  elementEtmis->EtSmeare = (PTmisS).Et();
+
+
 
       //*****************************
@@ -314 +342 @@
         if(JETSm.Pt()>1)
           {
+            float NonSmeareEt=(JET.E()*sin(JET.Theta()));
+            float SmeareEt=(JETSm.E()*sin(JETSm.Theta()));
+
             elementJet= (RESOLJET*) branchjet->NewEntry();
-            elementJet->NonSmearePT = JET.Pt();
-            elementJet->SmearePT = (JETSm.Pt()/JET.Pt());
-            /*cout<<"valeur obtenue "<<JETSm.Pt()/JET.Pt()<<endl;
-            cout<<" pt non smeare "<<JET.Pt()<<" phi "<<JET.Phi()<<" eta "<<JET.Eta()<<endl;
-            cout<<"pt smeare "<<JETSm.Pt()<<" phi "<<JETSm.Phi()<<" eta "<<JETSm.Eta()<<endl;
-            cout<<"************"<<endl;
-*/
+            elementJet->NonSmearePT = JET.Et();
+            elementJet->SmearePT = JETSm.Et()/JET.Et();
+          
           }