Fork me on GitHub

Changeset f382b7e in git


Ignore:
Timestamp:
Sep 23, 2020, 1:44:18 PM (4 years ago)
Author:
Michele Selvaggi <michele.selvaggi@…>
Branches:
master
Children:
2ca681a
Parents:
55f821a
Message:

now filling neutral hadrons in DR towers

File:
1 edited

Legend:

Unmodified
Added
Removed
  • modules/DualReadoutCalorimeter.cc

    r55f821a rf382b7e  
    645645    //cout<<"Creating tower with Pt, Eta, Phi, Energy: "<<pt<<","<<eta<<","<<phi<<","<<neutralEnergy<<endl;
    646646    tower->Momentum.SetPtEtaPhiE(pt, eta, phi, neutralEnergy);
    647     tower->Eem = neutralEnergy;
    648     tower->Ehad = 0.0;
    649     tower->PID = 22;
    650 
     647
     648    // if no hadronic energy, use ECAL resolution
     649    if (fHCalTowerEnergy <= fHCalEnergyMin)
     650    {
     651      tower->Eem = neutralEnergy;
     652      tower->Ehad = 0.0;
     653      tower->PID = 22;
     654    }
     655
     656    // if hadronic fraction > 0, use HCAL resolution
     657    else
     658    {
     659      tower->Eem = 0;
     660      tower->Ehad = neutralEnergy;
     661      tower->PID = 130;
     662    }
     663
     664    fEFlowPhotonOutputArray->Add(tower);
    651665   
    652 
    653     fEFlowPhotonOutputArray->Add(tower);
    654 
    655666
    656667    //clone tracks
Note: See TracChangeset for help on using the changeset viewer.