Fork me on GitHub

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


Ignore:
Timestamp:
May 10, 2009, 8:23:32 PM (15 years ago)
Author:
Xavier Rouby
Message:

added Resolution terms for energy and timing for ZDC

File:
1 edited

Legend:

Unmodified
Added
Removed

  • trunk/src/SmearUtil.cc

    r350 r374  
    6262  ELG_Nfwd          = 0.0;               // N term for FCAL
    6363  ELG_Cfwd          = 0.107;             // C term for FCAL
     64  ELG_Szdc          = 0.70;              // S term for ZDC
     65  ELG_Nzdc          = 0.0;               // N term for ZDC
     66  ELG_Czdc          = 0.08;              // C term for ZDC
    6467
    6568  // Energy resolution for hadrons in ecal/hcal/hf
    … …  
    7174  HAD_Nhf           = 0.;                // N term for FCAL
    7275  HAD_Chf           = 0.13;              // C term for FCAL
     76  HAD_Szdc          = 1.38;              // S term for ZDC
     77  HAD_Nzdc          = 0.;                // N term for ZDC
     78  HAD_Czdc          = 0.13;              // C term for ZDC
    7379
    7480  // Muon smearing
    7581  MU_SmearPt        = 0.01;
     82
     83  // time resolution
     84  ZDC_T_resolution   = 0;       // resolution for time measurement [s]
     85  RP220_T_resolution = 0;
     86  RP420_T_resolution = 0;
    7687
    7788  // Tracking efficiencies
    … …  
    95106
    96107
    97   // Thresholds for reconstructed objetcs
     108  // Thresholds for reconstructed objetcs (GeV)
    98109  PTCUT_elec      = 10.0;
    99110  PTCUT_muon      = 10.0;
    … …  
    101112  PTCUT_gamma     = 10.0;
    102113  PTCUT_taujet    = 10.0;
     114
     115  ZDC_gamma_E     = 20; // GeV
     116  ZDC_n_E         = 50; // GeV
    103117
    104118  // Isolation
    … …  
    199213  ELG_Sfwd          = DET.ELG_Sfwd;
    200214  ELG_Nfwd          = DET.ELG_Nfwd;
    201 
    202   // Energy resolution for hadrons in ecal/hcal/hf
     215  ELG_Czdc          = DET.ELG_Czdc;
     216  ELG_Szdc          = DET.ELG_Szdc;
     217  ELG_Nzdc          = DET.ELG_Nzdc;
     218
     219  // Energy resolution for hadrons in ecal/hcal/hf/zdc
    203220  HAD_Shcal         = DET.HAD_Shcal;
    204221  HAD_Nhcal         = DET.HAD_Nhcal;
    … …  
    207224  HAD_Nhf           = DET.HAD_Nhf;
    208225  HAD_Chf           = DET.HAD_Chf;
     226  HAD_Szdc          = DET.HAD_Szdc;
     227  HAD_Nzdc          = DET.HAD_Nzdc;
     228  HAD_Czdc          = DET.HAD_Czdc;
     229
     230  // time resolution
     231  ZDC_T_resolution   = DET.ZDC_T_resolution;       // resolution for time measurement [s]
     232  RP220_T_resolution = DET.RP220_T_resolution;
     233  RP420_T_resolution = DET.RP420_T_resolution;
    209234
    210235  // Muon smearing
    … …  
    229254  PTCUT_gamma     = DET.PTCUT_gamma;
    230255  PTCUT_taujet    = DET.PTCUT_taujet;
     256
     257  ZDC_gamma_E     = DET.ZDC_gamma_E;
     258  ZDC_n_E         = DET.ZDC_n_E;
    231259
    232260  // Isolation
    … …  
    321349  ELG_Sfwd          = DET.ELG_Sfwd;
    322350  ELG_Nfwd          = DET.ELG_Nfwd;
     351  ELG_Czdc          = DET.ELG_Czdc;
     352  ELG_Szdc          = DET.ELG_Szdc;
     353  ELG_Nzdc          = DET.ELG_Nzdc;
    323354
    324355  // Energy resolution for hadrons in ecal/hcal/hf
    … …  
    329360  HAD_Nhf           = DET.HAD_Nhf;
    330361  HAD_Chf           = DET.HAD_Chf;
     362  HAD_Szdc          = DET.HAD_Szdc;
     363  HAD_Nzdc          = DET.HAD_Nzdc;
     364  HAD_Czdc          = DET.HAD_Czdc;
     365
     366  // time resolution
     367  ZDC_T_resolution   = DET.ZDC_T_resolution;       // resolution for time measurement [s]
     368  RP220_T_resolution = DET.RP220_T_resolution;
     369  RP420_T_resolution = DET.RP420_T_resolution;
    331370
    332371  // Muon smearing
    … …  
    351390  PTCUT_gamma     = DET.PTCUT_gamma;
    352391  PTCUT_taujet    = DET.PTCUT_taujet;
     392
     393  ZDC_gamma_E     = DET.ZDC_gamma_E;
     394  ZDC_n_E         = DET.ZDC_n_E;
    353395
    354396  // Isolation
    … …  
    475517    else if(strstr(temp_string.c_str(),"ELG_Cfwd"))         {curstring >> varname >> value; ELG_Cfwd          = value;}
    476518    else if(strstr(temp_string.c_str(),"ELG_Nfwd"))         {curstring >> varname >> value; ELG_Nfwd          = value;}
     519    else if(strstr(temp_string.c_str(),"ELG_Szdc"))         {curstring >> varname >> value; ELG_Szdc          = value;}
     520    else if(strstr(temp_string.c_str(),"ELG_Czdc"))         {curstring >> varname >> value; ELG_Czdc          = value;}
     521    else if(strstr(temp_string.c_str(),"ELG_Nzdc"))         {curstring >> varname >> value; ELG_Nzdc          = value;}
     522
    477523    else if(strstr(temp_string.c_str(),"HAD_Shcal"))        {curstring >> varname >> value; HAD_Shcal         = value;}
    478524    else if(strstr(temp_string.c_str(),"HAD_Nhcal"))        {curstring >> varname >> value; HAD_Nhcal         = value;}
    … …  
    481527    else if(strstr(temp_string.c_str(),"HAD_Nhf"))          {curstring >> varname >> value; HAD_Nhf           = value;}
    482528    else if(strstr(temp_string.c_str(),"HAD_Chf"))          {curstring >> varname >> value; HAD_Chf           = value;}
     529    else if(strstr(temp_string.c_str(),"HAD_Szdc"))         {curstring >> varname >> value; HAD_Szdc          = value;}
     530    else if(strstr(temp_string.c_str(),"HAD_Nzdc"))         {curstring >> varname >> value; HAD_Nzdc          = value;}
     531    else if(strstr(temp_string.c_str(),"HAD_Czdc"))         {curstring >> varname >> value; HAD_Czdc          = value;}
     532    else if(strstr(temp_string.c_str(),"ZDC_T_resolution")) {curstring >> varname >> value; ZDC_T_resolution  = value;}
     533    else if(strstr(temp_string.c_str(),"RP220_T_resolution")) {curstring >> varname >> value; RP220_T_resolution  = value;}
     534    else if(strstr(temp_string.c_str(),"RP420_T_resolution")) {curstring >> varname >> value; RP420_T_resolution  = value;}
    483535    else if(strstr(temp_string.c_str(),"MU_SmearPt"))       {curstring >> varname >> value; MU_SmearPt        = value;}
    484536   
    … …  
    503555    else if(strstr(temp_string.c_str(),"PTCUT_gamma"))      {curstring >> varname >> value; PTCUT_gamma       = value;}
    504556    else if(strstr(temp_string.c_str(),"PTCUT_taujet"))     {curstring >> varname >> value; PTCUT_taujet      = value;}
     557    else if(strstr(temp_string.c_str(),"ZDC_gamma_E"))      {curstring >> varname >> value; ZDC_gamma_E      = value;}
     558    else if(strstr(temp_string.c_str(),"ZDC_n_E"))          {curstring >> varname >> value; ZDC_n_E      = value;}
    505559
    506560    else if(strstr(temp_string.c_str(),"ISOL_PT"))          {curstring >> varname >> value; ISOL_PT           = value;}
    … …  
    694748  f_out << left << setw(30) <<"* C term for FCAL: "<<""
    695749        << left << setw(30) <<ELG_Cfwd       <<""<< right << setw(10)<<"*"<<"\n";
     750  f_out << left << setw(30) <<"* S term for ZDC: "<<""
     751        << left << setw(30) <<ELG_Szdc       <<""<< right << setw(10)<<"*"<<"\n";
     752  f_out << left << setw(30) <<"* N term for ZDC: "<<""
     753        << left << setw(30) <<ELG_Nzdc       <<""<< right << setw(10)<<"*"<<"\n";
     754  f_out << left << setw(30) <<"* C term for ZDC: "<<""
     755        << left << setw(30) <<ELG_Czdc       <<""<< right << setw(10)<<"*"<<"\n";
     756
    696757  f_out<<"*                                                                    *"<<"\n";
    697758  f_out<<"#*****************************                                       *"<<"\n";
    … …  
    711772  f_out << left << setw(30) <<"* C term for FCAL: "<<""
    712773        << left << setw(30) <<HAD_Chf         <<""<< right << setw(10)<<"*"<<"\n";
     774  f_out << left << setw(30) <<"* S term for ZDC: "<<""
     775        << left << setw(30) <<HAD_Szdc        <<""<< right << setw(10)<<"*"<<"\n";
     776  f_out << left << setw(30) <<"* N term for ZDC: "<<""
     777        << left << setw(30) <<HAD_Nzdc        <<""<< right << setw(10)<<"*"<<"\n";
     778  f_out << left << setw(30) <<"* C term for ZDC: "<<""
     779        << left << setw(30) <<HAD_Czdc        <<""<< right << setw(10)<<"*"<<"\n";
     780
     781  f_out<<"*                                                                    *"<<"\n";
     782  f_out<<"#*************************                                           *"<<"\n";
     783  f_out<<"# Time smearing parameters                                           *"<<"\n";
     784  f_out<<"#*************************                                           *"<<"\n";
     785  f_out<<"*                                                                    *"<<"\n";
     786  f_out << left << setw(55) <<"* Time resolution for ZDC        :              "<<""
     787        << left << setw(5) <<ZDC_T_resolution       <<""<< right << setw(10)<<"*"<<"\n";
     788  f_out << left << setw(55) <<"* Time resolution for RP220      :              "<<""
     789        << left << setw(5) <<RP220_T_resolution     <<""<< right << setw(10)<<"*"<<"\n";
     790  f_out << left << setw(55) <<"* Time resolution for RP420      :              "<<""
     791        << left << setw(5) <<RP420_T_resolution     <<""<< right << setw(10)<<"*"<<"\n";
     792  f_out<<"*                                                                    *"<<"\n";
     793
    713794  f_out<<"*                                                                    *"<<"\n";
    714795  f_out<<"#*************************                                           *"<<"\n";
    … …  
    786867  f_out << left << setw(40) <<"* Minimum pT for photons: "<<""
    787868        << left << setw(20) <<PTCUT_gamma        <<""<< right << setw(10)<<"*"<<"\n";
     869  f_out << left << setw(40) <<"* Minimum E for photons in ZDC: "<<""
     870        << left << setw(20) <<ZDC_gamma_E        <<""<< right << setw(10)<<"*"<<"\n";
     871  f_out << left << setw(40) <<"* Minimum E for neutrons in ZDC: "<<""
     872        << left << setw(20) <<ZDC_n_E        <<""<< right << setw(10)<<"*"<<"\n";
     873
    788874  f_out<<"*                                                                    *"<<"\n";
    789875  f_out<<"#*******************                                                 *"<<"\n";
Note: See TracChangeset for help on using the changeset viewer.