Changeset 182 in svn for trunk/Resolutions.cpp

Timestamp:

Jan 14, 2009, 11:38:51 AM (16 years ago)

Author:

Xavier Rouby

Message:

pt_cut du bfield retire si bfield simule

File:

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

trunk/Resolutions.cpp

@@ -292 +292 @@
                int charge=Charge(pid);
                 if(recoMomentum.E() !=0 && pid != pMU) {
-                if(charge == 0 || (charge !=0 && recoMomentum.Pt() >= DET->TRACK_ptmin)){
+                        // in case the Bfield is not simulated, checks that charged particles have enough pt to reach the calos
+                  if ( !DET->FLAG_bfield && charge!=0 && genMomentum.Pt() <= DET->TRACK_ptmin ) { /* particules do not reach calos */ }
+                  else { // particles reach calos
                   DET->BinEtaPhi(recoMomentum.Phi(), recoMomentum.Eta(), iPhi, iEta);
                   if(iEta != -100 && iPhi != -100)
@@ -303 +305 @@
                 }
               }
-              
+             
               // all final charged particles
               if(
                  (recoMomentum.E()!=0) &&
                  (fabs(recoMomentum.Eta()) < DET->CEN_max_tracker) &&
-                 (recoMomentum.Pt() > DET->TRACK_ptmin ) &&     // pt too small to be taken into account
+                 (DET->FLAG_bfield || ( !DET->FLAG_bfield && genMomentum.Pt() > DET->TRACK_ptmin )) &&
+                        // if bfield not simulated, pt should be high enough to be taken into account
                  ((rand()%100) < DET->TRACK_eff)  &&
                  (charge!=0)