Fork me on GitHub

Ignore:
File:
1 edited

Legend:

Unmodified
Added
Removed
  • modules/DualReadoutCalorimeter.cc

    r5eda6767 r4afb18d  
    613613    // create new photon tower
    614614    tower = static_cast<Candidate*>(fTower->Clone());
    615     pt =  neutralEnergy / TMath::CosH(eta);
     615    pt = neutralEnergy / TMath::CosH(eta);
    616616    //cout<<"Creating tower with Pt, Eta, Phi, Energy: "<<pt<<","<<eta<<","<<phi<<","<<neutralEnergy<<endl;
    617617    tower->Momentum.SetPtEtaPhiE(pt, eta, phi, neutralEnergy);
     
    623623      tower->Ehad = 0.0;
    624624      tower->PID = 22;
    625     }
    626 
     625      fEFlowPhotonOutputArray->Add(tower);
     626    }
    627627    // if hadronic fraction > 0, use HCAL resolution
    628628    else
     
    631631      tower->Ehad = neutralEnergy;
    632632      tower->PID = 130;
    633     }
    634 
    635     fEFlowPhotonOutputArray->Add(tower);
    636 
     633      fEFlowNeutralHadronOutputArray->Add(tower);
     634    }
    637635
    638636    //clone tracks
Note: See TracChangeset for help on using the changeset viewer.