Fork me on GitHub

Ignore:
File:
1 edited

Legend:

Unmodified
Added
Removed
  • examples/Example3.C

    r6e41190 r3b99510  
    144144
    145145      particle = (GenParticle*) photon->Particles.At(0);
    146 
    147       plots->fPhotonDeltaPT->Fill((particle->PT - photon->PT)/particle->PT);
    148       plots->fPhotonDeltaEta->Fill((particle->Eta - photon->Eta)/particle->Eta);
    149       plots->fPhotonDeltaE->Fill((particle->E - photon->E)/particle->E);
     146      if(object->IsA() == GenParticle::Class())
     147      {
     148        plots->fPhotonDeltaPT->Fill((particle->PT - photon->PT)/particle->PT);
     149        plots->fPhotonDeltaEta->Fill((particle->Eta - photon->Eta)/particle->Eta);
     150        plots->fPhotonDeltaE->Fill((particle->E - photon->E)/particle->E);
     151      }   
    150152    }
    151153
Note: See TracChangeset for help on using the changeset viewer.