Fork me on GitHub

Changeset 40 in svn


Ignore:
Timestamp:
Nov 18, 2008, 10:31:32 AM (16 years ago)
Author:
severine ovyn
Message:

ETmis using CaloTowers

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/Delphes.cpp

    r38 r40  
    338338            // for calorimetric towers and mission PT
    339339            if(genMomentum.E()!=0) {
    340               PTmis = PTmis + genMomentum;//ptmis
    341340              if(pid !=pMU) {
    342341                PhysicsTower CaloTower = PhysicsTower(LorentzVector(genMomentum.Px(),genMomentum.Py(),genMomentum.Pz(), genMomentum.E()));
     
    345344                // back of the input_particles vector
    346345                input_particles.push_back(fastjet::PseudoJet(genMomentum.Px(),genMomentum.Py(),genMomentum.Pz(), genMomentum.E()));
     346                //genMomentumCalo.SetPtEtaPhiE(CaloTower.Et(),CaloTower.iEta(),CaloTower.iPhi(),CaloTower.fourVector.E);
    347347                genMomentumCalo.SetPxPyPzE(CaloTower.fourVector.px,CaloTower.fourVector.py,CaloTower.fourVector.pz,CaloTower.fourVector.E);
    348348                elementCalo = (TRootCalo*) branchCalo->NewEntry();
     
    472472     
    473473      // computes the Missing Transverse Momentum
     474      TLorentzVector Att(0.,0.,0.,0.);
     475      for(unsigned int i=0; i < towers.size(); i++)
     476        {
     477          Att.SetPxPyPzE(towers[i].fourVector.px,towers[i].fourVector.py,towers[i].fourVector.pz,towers[i].fourVector.E);
     478          PTmis = PTmis + Att;
     479        }
    474480      elementEtmis = (TRootETmis*) branchETmis->NewEntry();
    475       elementEtmis->ET = (-PTmis).Pt();
     481      elementEtmis->ET = (PTmis).Pt();
    476482      elementEtmis->Phi = (-PTmis).Phi();
    477483      elementEtmis->Px = (-PTmis).Px();
Note: See TracChangeset for help on using the changeset viewer.