Fork me on GitHub

Changeset 49 in svn


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

2 new jet algorithms

Location:
trunk
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • trunk/Delphes.cpp

    r47 r49  
    252252  case 4: {
    253253    jet_def = fastjet::JetDefinition(fastjet::kt_algorithm, DET->CONERADIUS);
    254     //jet_defs[4] = fastjet::JetDefinition(fastjet::cambridge_algorithm,jet_radius);
    255     //jet_defs[5] = fastjet::JetDefinition(fastjet::antikt_algorithm,jet_radius);
     254  }
     255    break;
     256   
     257  case 5: {
     258    jet_def = fastjet::JetDefinition(fastjet::cambridge_algorithm,DET->CONERADIUS);
     259  }
     260    break;
     261   
     262  case 6: {
     263    jet_def = fastjet::JetDefinition(fastjet::antikt_algorithm,DET->CONERADIUS);
    256264  }
    257265    break;
    … …  
    375383                TrackCentral.push_back(genMomentum);
    376384              }
    377          
    378 } // switch
     385           
     386          } // switch
    379387         
    380388          // Forward particles in CASTOR ?
    … …  
    504512        // Tau jet identification : 1! track and electromagnetic collimation
    505513        if(fabs(JET.Eta()) < (DET->MAX_TRACKER - DET->TAU_CONE_TRACKS)) {
    506         double Energie_tau_central = DET->EnergySmallCone(towers,JET.Eta(),JET.Phi());
    507         if(
    508            ( Energie_tau_central/JET.E() > DET->TAU_EM_COLLIMATION ) &&
    509            ( DET->NumTracks(TrackCentral,DET->PT_TRACK_TAU,JET.Eta(),JET.Phi()) == 1 ) &&
    510            ( JET.Pt() > DET->TAUJET_pt)
    511           ) {
    512                elementTauJet = (TRootTauJet*) branchTauJet->NewEntry();
    513                elementTauJet->Set(JET);
    514             } // if tau jet
    515            } // if JET.eta < tracker - tau_cone : Tau jet identification
    516  
     514          double Energie_tau_central = DET->EnergySmallCone(towers,JET.Eta(),JET.Phi());
     515          if(
     516             ( Energie_tau_central/JET.E() > DET->TAU_EM_COLLIMATION ) &&
     517             ( DET->NumTracks(TrackCentral,DET->PT_TRACK_TAU,JET.Eta(),JET.Phi()) == 1 ) &&
     518             ( JET.Pt() > DET->TAUJET_pt)
     519             ) {
     520            elementTauJet = (TRootTauJet*) branchTauJet->NewEntry();
     521            elementTauJet->Set(JET);
     522          } // if tau jet
     523        } // if JET.eta < tracker - tau_cone : Tau jet identification
     524       
    517525        if(JET.Pt() > DET->JET_pt)
    518526          {

  • 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.