Fork me on GitHub

Changeset 32 in svn for trunk/src/SmearUtil.cc


Ignore:
Timestamp:
Nov 17, 2008, 8:49:25 PM (16 years ago)
Author:
severine ovyn
Message:

remove HCALecal

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/src/SmearUtil.cc

    r31 r32  
    122122    else if(strstr(temp_string.c_str(),"ELG_Cfwd")){curstring >> varname >> value; ELG_Cfwd = value;}
    123123    else if(strstr(temp_string.c_str(),"ELG_Nfwd")){curstring >> varname >> value; ELG_Nfwd = value;}
    124     else if(strstr(temp_string.c_str(),"HAD_Secal")){curstring >> varname >> value; HAD_Secal = value;}
    125     else if(strstr(temp_string.c_str(),"HAD_Necal")){curstring >> varname >> value; HAD_Necal = value;}
    126     else if(strstr(temp_string.c_str(),"HAD_Cecal")){curstring >> varname >> value; HAD_Cecal = value;}
    127124    else if(strstr(temp_string.c_str(),"HAD_Shcal")){curstring >> varname >> value; HAD_Shcal = value;}
    128125    else if(strstr(temp_string.c_str(),"HAD_Nhcal")){curstring >> varname >> value; HAD_Nhcal = value;}
     
    227224     
    228225   energyS2 =  gRandom->Gaus(energy_ecal, sqrt(
    229                                       pow(HAD_Necal,2) +
    230                                       pow(HAD_Cecal*energy_ecal,2) +
    231                                       pow(HAD_Secal*sqrt(energy_ecal),2) ) );     
     226                                      pow(ELG_Ncen,2) +
     227                                      pow(ELG_Ccen*energy_ecal,2) +
     228                                      pow(ELG_Scen*sqrt(energy_ecal),2) ) );     
    232229
    233230  energyS = ((energyS1>0)?energyS1:0) + ((energyS2>0)?energyS2:0);
Note: See TracChangeset for help on using the changeset viewer.