Changeset 33 in svn for trunk/src/SmearUtil.cc

Timestamp:

Nov 17, 2008, 9:01:32 PM (16 years ago)

Author:

severine ovyn

Message:

add offline Thresholds

File:

• trunk/src/SmearUtil.cc (modified) (4 diffs)

trunk/src/SmearUtil.cc

@@ -47 +47 @@
 ELG_Nfwd =        0.0   ;                 // N term for central ECAL
 
-HAD_Secal =       0.05 ;                    // S term for central ECAL // electromagnetic calorimeter
-HAD_Necal =       0.25 ;                    // N term for central ECAL
-HAD_Cecal =       0.0055 ;                  // C term for central ECAL
 HAD_Shcal =       0.91 ;                    // S term for central HCAL // hadronic calorimeter
 HAD_Nhcal  =      0.   ;                    // N term for central HCAL
@@ -58 +55 @@
 
 MU_SmearPt =      0.01 ;
+
+ELEC_pt =         10.0;
+MUON_pt =         10.0;
+JET_pt =          20.0;
+TAUJET_pt =       10.0;
+
 
 TAU_CONE_ENERGY = 0.15  ;                  // Delta R = radius of the cone // for "electromagnetic collimation"
@@ -139 +142 @@
     else if(strstr(temp_string.c_str(),"JETALGO")){curstring >> varname >> value; JETALGO = (int)value;}
     else if(strstr(temp_string.c_str(),"TRACKING_EFF")){curstring >> varname >> value; TRACKING_EFF = (int)value;}
+    else if(strstr(temp_string.c_str(),"ELEC_pt")){curstring >> varname >> value; ELEC_pt = value;}
+    else if(strstr(temp_string.c_str(),"MUON_pt")){curstring >> varname >> value; MUON_pt = value;}
+    else if(strstr(temp_string.c_str(),"JET_pt")){curstring >> varname >> value; JET_pt = value;}
+    else if(strstr(temp_string.c_str(),"TAUJET_pt")){curstring >> varname >> value; TAUJET_pt = value;}
+
   }
 
@@ -177 +185 @@
                                          pow(ELG_Sfwd*sqrt(energy),2) ) ); 
   }
-float theta=electron.Theta();
-float phi=electron.Phi();
-  float px = energyS * sin(theta) * cos(phi);
-  float py = energyS * sin(theta) * sin(phi);
-  float pz = energyS * cos(theta);
-  //cout<<" px "<<px<<" py "<<py<<" pz "<<pz<<endl;
   electron.SetPtEtaPhiE(energyS/cosh(electron.Eta()), electron.Eta(), electron.Phi(), energyS);
-  //cout<<"nous "<<electron.Px()<<"  "<<electron.Py()<<"  "<<electron.Pz()<<endl;
   if(electron.E() < 0)electron.SetPxPyPzE(0,0,0,0); // no negative values after smearing !
 }