Fork me on GitHub

Changeset 182 in svn for trunk


Ignore:
Timestamp:
Jan 14, 2009, 11:38:51 AM (16 years ago)
Author:
Xavier Rouby
Message:

pt_cut du bfield retire si bfield simule

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/Resolutions.cpp

    r152 r182  
    292292               int charge=Charge(pid);
    293293                if(recoMomentum.E() !=0 && pid != pMU) {
    294                 if(charge == 0 || (charge !=0 && recoMomentum.Pt() >= DET->TRACK_ptmin)){
     294                        // in case the Bfield is not simulated, checks that charged particles have enough pt to reach the calos
     295                  if ( !DET->FLAG_bfield && charge!=0 && genMomentum.Pt() <= DET->TRACK_ptmin ) { /* particules do not reach calos */ }
     296                  else { // particles reach calos
    295297                  DET->BinEtaPhi(recoMomentum.Phi(), recoMomentum.Eta(), iPhi, iEta);
    296298                  if(iEta != -100 && iPhi != -100)
     
    303305                }
    304306              }
    305               
     307             
    306308              // all final charged particles
    307309              if(
    308310                 (recoMomentum.E()!=0) &&
    309311                 (fabs(recoMomentum.Eta()) < DET->CEN_max_tracker) &&
    310                  (recoMomentum.Pt() > DET->TRACK_ptmin ) &&     // pt too small to be taken into account
     312                 (DET->FLAG_bfield || ( !DET->FLAG_bfield && genMomentum.Pt() > DET->TRACK_ptmin )) &&
     313                        // if bfield not simulated, pt should be high enough to be taken into account
    311314                 ((rand()%100) < DET->TRACK_eff)  &&
    312315                 (charge!=0)
Note: See TracChangeset for help on using the changeset viewer.