Fork me on GitHub

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


Ignore:
Timestamp:
Feb 8, 2009, 11:28:45 PM (16 years ago)
Author:
severine ovyn
Message:

small changes

File:
1 edited

Legend:

Unmodified
Added
Removed

  • trunk/src/SmearUtil.cc

    r254 r257  
    5151  ELG_Ncen          = 0.25;              // N term for central ECAL
    5252  ELG_Ccen          = 0.005;             // C term for central ECAL
    53   ELG_Cfwd          = 0.107;             // S term for forward ECAL
    54   ELG_Sfwd          = 2.084;             // C term for forward ECAL
    55   ELG_Nfwd          = 0.0;               // N term for central ECAL
     53  ELG_Sfwd          = 2.084;             // S term for FCAL
     54  ELG_Nfwd          = 0.0;               // N term for FCAL
     55  ELG_Cfwd          = 0.107;             // C term for FCAL
    5656
    5757  // Energy resolution for hadrons in ecal/hcal/hf
    5858  // \sigma/E = C + N/E + S/\sqrt{E}
    59   HAD_Shcal         = 1.5;               // S term for central HCAL // hadronic calorimeter
     59  HAD_Shcal         = 1.5;               // S term for central HCAL
    6060  HAD_Nhcal         = 0.;                // N term for central HCAL
    6161  HAD_Chcal         = 0.05;              // C term for central HCAL
    62   HAD_Shf           = 2.7;               // S term for HF // forward calorimeter
    63   HAD_Nhf           = 0.;                // N term for HF
    64   HAD_Chf           = 0.13;              // C term for HF
     62  HAD_Shf           = 2.7;               // S term for FCAL
     63  HAD_Nhf           = 0.;                // N term for FCAL
     64  HAD_Chf           = 0.13;              // C term for FCAL
    6565
    6666  // Muon smearing
    … …  
    127127  RP_420_s          = 420;                   // distance of the RP to the IP, in meters
    128128  RP_420_x          = 0.004;                 // distance of the RP to the beam, in meters
     129  RP_IP_name        = "IP5";
    129130  RP_beam1Card      = "data/LHCB1IR5_v6.500.tfs";
    130131  RP_beam2Card      = "data/LHCB1IR5_v6.500.tfs";
    … …  
    255256  RP_cross_y        = DET.RP_cross_y;
    256257  RP_cross_ang      = DET.RP_cross_ang;
     258  RP_IP_name        = DET.RP_IP_name;
    257259
    258260
    … …  
    368370  RP_cross_y        = DET.RP_cross_y;
    369371  RP_cross_ang      = DET.RP_cross_ang;
    370 
     372  RP_IP_name        = DET.RP_IP_name;
    371373
    372374
    … …  
    433435    else if(strstr(temp_string.c_str(),"RP_420_s"))         {curstring >> varname >> value; RP_420_s          = value;}
    434436    else if(strstr(temp_string.c_str(),"RP_420_x"))         {curstring >> varname >> value; RP_420_x          = value;}
    435     else if(strstr(temp_string.c_str(),"Beam1Card"))        {curstring >> varname >> svalue;RP_beam1Card         = svalue;}
    436     else if(strstr(temp_string.c_str(),"Beam2Card"))        {curstring >> varname >> svalue;RP_beam2Card         = svalue;}
     437    else if(strstr(temp_string.c_str(),"RP_beam1Card"))     {curstring >> varname >> svalue;RP_beam1Card      = svalue;}
     438    else if(strstr(temp_string.c_str(),"RP_beam2Card"))     {curstring >> varname >> svalue;RP_beam2Card      = svalue;}
     439    else if(strstr(temp_string.c_str(),"RP_IP_name"))       {curstring >> varname >> svalue;RP_IP_name        = svalue;}
    437440   
    438441    else if(strstr(temp_string.c_str(),"ELG_Scen"))         {curstring >> varname >> value; ELG_Scen          = value;}
    … …  
    575578    f_out << left << setw(55) <<"* Distance of the RP to the beam, in meters: "<<""
    576579          << left << setw(5) <<RP_420_x             <<""<< right << setw(10)<<"*"<<"\n";
     580    f_out << left << setw(55) <<"* Interaction point at the LHC named: "<<""
     581          << left << setw(5) <<RP_IP_name             <<""<< right << setw(10)<<"*"<<"\n";
    577582    f_out << left << setw(35) <<"* Datacard for beam 1: "<<""
    578583          << left << setw(25) <<RP_beam1Card            <<""<< right << setw(10)<<"*"<<"\n";
    … …  
    609614  f_out << left << setw(30) <<"* C term for central ECAL: "<<""
    610615        << left << setw(30) <<ELG_Ccen       <<""<< right << setw(10)<<"*"<<"\n";
    611   f_out << left << setw(30) <<"* S term for forward ECAL: "<<""
     616  f_out << left << setw(30) <<"* S term for FCAL: "<<""
    612617        << left << setw(30) <<ELG_Sfwd       <<""<< right << setw(10)<<"*"<<"\n";
    613   f_out << left << setw(30) <<"* N term for forward ECAL: "<<""
     618  f_out << left << setw(30) <<"* N term for FCAL: "<<""
    614619        << left << setw(30) <<ELG_Nfwd       <<""<< right << setw(10)<<"*"<<"\n";
    615   f_out << left << setw(30) <<"* C term for forward ECAL: "<<""
     620  f_out << left << setw(30) <<"* C term for FCAL: "<<""
    616621        << left << setw(30) <<ELG_Cfwd       <<""<< right << setw(10)<<"*"<<"\n";
    617622  f_out<<"*                                                                    *"<<"\n";
    … …  
    626631  f_out << left << setw(30) <<"* C term for central HCAL: "<<""
    627632        << left << setw(30) <<HAD_Chcal       <<""<< right << setw(10)<<"*"<<"\n";
    628   f_out << left << setw(30) <<"* S term for forward HCAL: "<<""
     633  f_out << left << setw(30) <<"* S term for FCAL: "<<""
    629634        << left << setw(30) <<HAD_Shf         <<""<< right << setw(10)<<"*"<<"\n";
    630   f_out << left << setw(30) <<"* N term for forward HCAL: "<<""
     635  f_out << left << setw(30) <<"* N term for FCAL: "<<""
    631636        << left << setw(30) <<HAD_Nhf         <<""<< right << setw(10)<<"*"<<"\n";
    632   f_out << left << setw(30) <<"* C term for forward HCAL: "<<""
     637  f_out << left << setw(30) <<"* C term for FCAL: "<<""
    633638        << left << setw(30) <<HAD_Chf         <<""<< right << setw(10)<<"*"<<"\n";
    634639  f_out<<"*                                                                    *"<<"\n";
Note: See TracChangeset for help on using the changeset viewer.