Fork me on GitHub

Changeset 49 in svn for trunk/src


Ignore:
Timestamp:
Nov 20, 2008, 11:28:51 AM (16 years ago)
Author:
severine ovyn
Message:

2 new jet algorithms

File:
1 edited

Legend:

Unmodified
Added
Removed

  • trunk/src/SmearUtil.cc

    r48 r49  
    353353  f_out<<"#***************                                                     *"<<"\n";
    354354  f_out<<"*                                                                    *"<<"\n";
     355  f_out<<"*  Six algorithms are currently available:                           *"<<"\n";
     356  f_out<<"*         - 1) CDF cone algorithm,                                   *"<<"\n";
     357  f_out<<"*         - 2) CDF MidPoint algorithm,                               *"<<"\n";
     358  f_out<<"*         - 3) SIScone algorithm,                                    *"<<"\n";
     359  f_out<<"*         - 4) kt algorithm,                                         *"<<"\n";
     360  f_out<<"*         - 5) Cambrigde/Aachen algorithm,                           *"<<"\n";
     361  f_out<<"*         - 6) Anti-kt algorithm.                                    *"<<"\n";
     362  f_out<<"*                                                                    *"<<"\n";
     363  f_out<<"* You have chosen                                                    *"<<"\n";
    355364  switch(JETALGO) {
    356365  default:
    357366  case 1: {
    358     f_out<<"* CDF JetClu jet algorithm with:                                     *"<<"\n";
    359     f_out << left << setw(40) <<"*        - Seed threshold:    "<<""
    360           << left << setw(10) <<SEEDTHRESHOLD                   <<""<< right << setw(20)<<"! not in datacard  *"<<"\n";
    361     f_out << left << setw(40) <<"*        - Cone radius:       "<<""
    362           << left << setw(10) <<CONERADIUS                      <<""<< right << setw(20)<<"*"<<"\n";
    363     f_out << left << setw(40) <<"*        - Adjacency cut:     "<<""
    364           << left << setw(10) <<C_ADJACENCYCUT                  <<""<< right << setw(20)<<"! not in datacard  *"<<"\n";
    365     f_out << left << setw(40) <<"*        - Max iterations:    "<<""
    366           << left << setw(10) <<C_MAXITERATIONS                 <<""<< right << setw(20)<<"! not in datacard  *"<<"\n";
    367     f_out << left << setw(40) <<"*        - Iratch:            "<<""
    368           << left << setw(10) <<C_IRATCH                        <<""<< right << setw(20)<<"! not in datacard  *"<<"\n";
    369     f_out << left << setw(40) <<"*        - Overlap threshold: "<<""
    370           << left << setw(10) <<OVERLAPTHRESHOLD                <<""<< right << setw(20)<<"! not in datacard  *"<<"\n";
     367  f_out<<"* CDF JetClu jet algorithm with parameters:                          *"<<"\n";
     368  f_out << left << setw(40) <<"*        - Seed threshold:    "<<""
     369        << left << setw(10) <<SEEDTHRESHOLD                   <<""<< right << setw(20)<<"! not in datacard  *"<<"\n";
     370  f_out << left << setw(40) <<"*        - Cone radius:       "<<""
     371        << left << setw(10) <<CONERADIUS                      <<""<< right << setw(20)<<"*"<<"\n";
     372  f_out << left << setw(40) <<"*        - Adjacency cut:     "<<""
     373        << left << setw(10) <<C_ADJACENCYCUT                  <<""<< right << setw(20)<<"! not in datacard  *"<<"\n";
     374  f_out << left << setw(40) <<"*        - Max iterations:    "<<""
     375        << left << setw(10) <<C_MAXITERATIONS                 <<""<< right << setw(20)<<"! not in datacard  *"<<"\n";
     376  f_out << left << setw(40) <<"*        - Iratch:            "<<""
     377        << left << setw(10) <<C_IRATCH                        <<""<< right << setw(20)<<"! not in datacard  *"<<"\n";
     378  f_out << left << setw(40) <<"*        - Overlap threshold: "<<""
     379        << left << setw(10) <<OVERLAPTHRESHOLD                <<""<< right << setw(20)<<"! not in datacard  *"<<"\n";
    371380  }
    372381    break;
    373382  case 2: {
    374     f_out<<"* CDF midpoint jet algorithm with:                                     *"<<"\n";
    375     f_out << left << setw(40) <<"*        - Seed threshold:    "<<""
    376           << left << setw(20) <<SEEDTHRESHOLD                   <<""<< right << setw(10)<<"! not in datacard  *"<<"\n";
    377     f_out << left << setw(40) <<"*        - Cone radius:       "<<""
    378           << left << setw(20) <<CONERADIUS                      <<""<< right << setw(10)<<"*"<<"\n";
    379     f_out << left << setw(40) <<"*        - Cone area fraction:"<<""
    380           << left << setw(20) <<M_CONEAREAFRACTION              <<""<< right << setw(10)<<"! not in datacard  *"<<"\n";
    381     f_out << left << setw(40) <<"*        - Maximum pair size: "<<""
    382           << left << setw(20) <<M_MAXPAIRSIZE                   <<""<< right << setw(10)<<"! not in datacard  *"<<"\n";
    383     f_out << left << setw(40) <<"*        - Max iterations:    "<<""
    384           << left << setw(20) <<M_MAXITERATIONS                 <<""<< right << setw(10)<<"! not in datacard  *"<<"\n";
    385     f_out << left << setw(40) <<"*        - Overlap threshold: "<<""
    386           << left << setw(20) <<OVERLAPTHRESHOLD                <<""<< right << setw(10)<<"! not in datacard  *"<<"\n";
     383  f_out<<"* CDF midpoint jet algorithm with parameters:                        *"<<"\n";
     384  f_out << left << setw(40) <<"*        - Seed threshold:    "<<""
     385        << left << setw(20) <<SEEDTHRESHOLD                   <<""<< right << setw(10)<<"! not in datacard  *"<<"\n";
     386  f_out << left << setw(40) <<"*        - Cone radius:       "<<""
     387        << left << setw(20) <<CONERADIUS                      <<""<< right << setw(10)<<"*"<<"\n";
     388  f_out << left << setw(40) <<"*        - Cone area fraction:"<<""
     389        << left << setw(20) <<M_CONEAREAFRACTION              <<""<< right << setw(10)<<"! not in datacard  *"<<"\n";
     390  f_out << left << setw(40) <<"*        - Maximum pair size: "<<""
     391        << left << setw(20) <<M_MAXPAIRSIZE                   <<""<< right << setw(10)<<"! not in datacard  *"<<"\n";
     392  f_out << left << setw(40) <<"*        - Max iterations:    "<<""
     393        << left << setw(20) <<M_MAXITERATIONS                 <<""<< right << setw(10)<<"! not in datacard  *"<<"\n";
     394  f_out << left << setw(40) <<"*        - Overlap threshold: "<<""
     395        << left << setw(20) <<OVERLAPTHRESHOLD                <<""<< right << setw(10)<<"! not in datacard  *"<<"\n";
    387396  }
    388397    break;
    389398  case 3: {
    390     f_out<<"* SISCone jet algorithm with:                                     *"<<"\n";
    391     f_out << left << setw(40) <<"*        - Cone radius:             "<<""
    392           << left << setw(20) <<CONERADIUS                            <<""<< right << setw(10)<<"*"<<"\n";
    393     f_out << left << setw(40) <<"*        - Overlap threshold:       "<<""
    394           << left << setw(20) <<OVERLAPTHRESHOLD                      <<""<< right << setw(10)<<"! not in datacard  *"<<"\n";
    395     f_out << left << setw(40) <<"*        - Number pass max:         "<<""
    396           << left << setw(20) <<NPASS                                 <<""<< right << setw(10)<<"! not in datacard  *"<<"\n";
    397     f_out << left << setw(40) <<"*        - Minimum pT for protojet: "<<""
    398           << left << setw(20) <<PROTOJET_PTMIN                        <<""<< right << setw(10)<<"! not in datacard  *"<<"\n";
     399  f_out <<"* SISCone jet algorithm with parameters:                            *"<<"\n";
     400  f_out << left << setw(40) <<"*        - Cone radius:             "<<""
     401        << left << setw(20) <<CONERADIUS                            <<""<< right << setw(10)<<"*"<<"\n";
     402  f_out << left << setw(40) <<"*        - Overlap threshold:       "<<""
     403        << left << setw(20) <<OVERLAPTHRESHOLD                      <<""<< right << setw(10)<<"! not in datacard  *"<<"\n";
     404  f_out << left << setw(40) <<"*        - Number pass max:         "<<""
     405        << left << setw(20) <<NPASS                                 <<""<< right << setw(10)<<"! not in datacard  *"<<"\n";
     406  f_out << left << setw(40) <<"*        - Minimum pT for protojet: "<<""
     407        << left << setw(20) <<PROTOJET_PTMIN                        <<""<< right << setw(10)<<"! not in datacard  *"<<"\n";
    399408  }
    400409    break;
    401410  case 4: {
    402     f_out<<"* KT jet algorithm with:                                          *"<<"\n";
    403     f_out << left << setw(40) <<"*        - Cone radius: "<<""
    404           << left << setw(20) <<CONERADIUS                <<""<< right << setw(10)<<"*"<<"\n";
     411  f_out <<"* KT jet algorithm with parameters:                                 *"<<"\n";
     412  f_out << left << setw(40) <<"*        - Cone radius: "<<""
     413        << left << setw(20) <<CONERADIUS                <<""<< right << setw(10)<<"*"<<"\n";
    405414  }
    406415    break;
     416  case 5: {
     417  f_out <<"* Cambridge/Aachen jet algorithm with parameters:                   *"<<"\n";
     418  f_out << left << setw(40) <<"*        - Cone radius: "<<""
     419        << left << setw(20) <<CONERADIUS                <<""<< right << setw(10)<<"*"<<"\n";
     420  }
     421    break;
     422  case 6: {
     423  f_out <<"* Anti-kt jet algorithm with parameters:                            *"<<"\n";
     424  f_out << left << setw(40) <<"*        - Cone radius: "<<""
     425        << left << setw(20) <<CONERADIUS                <<""<< right << setw(10)<<"*"<<"\n";
     426  }
     427    break;
     428
     429
    407430  }
    408431  f_out<<"*                                                                    *"<<"\n";
Note: See TracChangeset for help on using the changeset viewer.