Fork me on GitHub

Changeset 43 in svn


Ignore:
Timestamp:
Nov 19, 2008, 2:06:15 PM (16 years ago)
Author:
severine ovyn
Message:

general jet variable

Location:
trunk
Files:
4 edited

Legend:

Unmodified
Added
Removed

  • trunk/Delphes.cpp

    r40 r43  
    168168  vector<TLorentzVector> TrackCentral; 
    169169  vector<PhysicsTower> towers;
    170  
    171170  vector<fastjet::PseudoJet> input_particles;//for FastJet algorithm
    172171  vector<fastjet::PseudoJet> inclusive_jets;
    … …  
    214213    // set up a CDF midpoint jet definition
    215214#ifdef ENABLE_PLUGIN_CDFCONES
    216     plugins = new fastjet::CDFJetCluPlugin(DET->C_SEEDTHRESHOLD,DET->CONERADIUS,DET->C_ADJACENCYCUT,DET->C_MAXITERATIONS,DET->C_IRATCH,DET->C_OVERLAPTHRESHOLD);
     215    plugins = new fastjet::CDFJetCluPlugin(DET->SEEDTHRESHOLD,DET->CONERADIUS,DET->C_ADJACENCYCUT,DET->C_MAXITERATIONS,DET->C_IRATCH,DET->OVERLAPTHRESHOLD);
    217216    jet_def = fastjet::JetDefinition(plugins);
    218217#else
    … …  
    226225    // set up a CDF midpoint jet definition
    227226#ifdef ENABLE_PLUGIN_CDFCONES
    228     plugins = new fastjet::CDFMidPointPlugin (DET->M_SEEDTHRESHOLD,DET->CONERADIUS,DET->M_CONEAREAFRACTION,DET->M_MAXPAIRSIZE,DET->M_MAXPAIRSIZE,DET->C_OVERLAPTHRESHOLD);
     227    plugins = new fastjet::CDFMidPointPlugin (DET->SEEDTHRESHOLD,DET->CONERADIUS,DET->M_CONEAREAFRACTION,DET->M_MAXPAIRSIZE,DET->M_MAXPAIRSIZE,DET->OVERLAPTHRESHOLD);
    229228    jet_def = fastjet::JetDefinition(plugins);
    230229#else
    … …  
    238237    int npass = 0;               // do infinite number of passes
    239238    double protojet_ptmin = 0.0; // use all protojets
    240     plugins = new fastjet::SISConePlugin (DET->CONERADIUS,DET->C_OVERLAPTHRESHOLD,npass, protojet_ptmin);
     239    plugins = new fastjet::SISConePlugin (DET->CONERADIUS,DET->OVERLAPTHRESHOLD,npass, protojet_ptmin);
    241240    jet_def = fastjet::JetDefinition(plugins);
    242241#else
    … …  
    262261      treeReader->ReadEntry(entry);
    263262      treeWriter->Clear();
    264      
    265263      if((entry % 100) == 0 && entry > 0 )  cout << "** Processing element # " << entry << endl;
    266264     
    … …  
    337335            // all final particles but muons and neutrinos     
    338336            // for calorimetric towers and mission PT
    339             if(genMomentum.E()!=0) {
     337           
     338            if(genMomentum.E() !=0) {
    340339              if(pid !=pMU) {
    341340                PhysicsTower CaloTower = PhysicsTower(LorentzVector(genMomentum.Px(),genMomentum.Py(),genMomentum.Pz(), genMomentum.E()));
    … …  
    371370                TrackCentral.push_back(genMomentum);
    372371              }
    373           } // switch
     372         
     373} // switch
    374374         
    375375          // Forward particles in CASTOR ?
    … …  
    470470      }
    471471     
    472      
    473472      // computes the Missing Transverse Momentum
    474473      TLorentzVector Att(0.,0.,0.,0.);
    … …  
    483482      elementEtmis->Px = (-PTmis).Px();
    484483      elementEtmis->Py = (-PTmis).Py();
    485      
    486484      //*****************************
    487485     
    488486      // run the jet clustering with the above jet definition
     487      cout<<input_particles.size()<<endl;
    489488      if(input_particles.size()!=0)
    490489        {
    491490          fastjet::ClusterSequence clust_seq(input_particles, jet_def);
    492          
    493          
    494491          // extract the inclusive jets with pt > 5 GeV
    495492          double ptmin = 5.0;
    496493          inclusive_jets = clust_seq.inclusive_jets(ptmin);
    497          
     494cout<<"inclusive_jets "<<inclusive_jets.size()<<endl;     
    498495          // sort jets into increasing pt
    499496          sorted_jets = sorted_by_pt(inclusive_jets);
    500497        }
     498cout<<"sorted_jets.size() "<<sorted_jets.size()<<endl;
    501499      for (unsigned int i = 0; i < sorted_jets.size(); i++) {
    502500        TLorentzVector JET;

  • trunk/Resolutions.cpp

    r39 r43  
    224224  // set up a CDF midpoint jet definition
    225225#ifdef ENABLE_PLUGIN_CDFCONES
    226   plugins = new fastjet::CDFJetCluPlugin(0,DET->CONERADIUS,DET->C_ADJACENCYCUT,DET->C_MAXITERATIONS,DET->C_IRATCH,DET->C_OVERLAPTHRESHOLD);
     226  plugins = new fastjet::CDFJetCluPlugin(0,DET->CONERADIUS,DET->C_ADJACENCYCUT,DET->C_MAXITERATIONS,DET->C_IRATCH,DET->OVERLAPTHRESHOLD);
    227227  jet_def = fastjet::JetDefinition(plugins);
    228228#else
    … …  
    232232  // set up a CDF midpoint jet definition
    233233#ifdef ENABLE_PLUGIN_CDFCONES
    234   pluginsS = new fastjet::CDFJetCluPlugin(1,DET->CONERADIUS,DET->C_ADJACENCYCUT,DET->C_MAXITERATIONS,DET->C_IRATCH,DET->C_OVERLAPTHRESHOLD);
     234  pluginsS = new fastjet::CDFJetCluPlugin(1,DET->CONERADIUS,DET->C_ADJACENCYCUT,DET->C_MAXITERATIONS,DET->C_IRATCH,DET->OVERLAPTHRESHOLD);
    235235  jet_defS = fastjet::JetDefinition(pluginsS);
    236236#else
    … …  
    240240 // set up a CDF midpoint jet definition
    241241 #ifdef ENABLE_PLUGIN_CDFCONES
    242  pluginsT = new fastjet::CDFJetCluPlugin(0,DET->CONERADIUS,DET->C_ADJACENCYCUT,DET->C_MAXITERATIONS,DET->C_IRATCH,DET->C_OVERLAPTHRESHOLD);
     242 pluginsT = new fastjet::CDFJetCluPlugin(0,DET->CONERADIUS,DET->C_ADJACENCYCUT,DET->C_MAXITERATIONS,DET->C_IRATCH,DET->OVERLAPTHRESHOLD);
    243243 jet_defT = fastjet::JetDefinition(pluginsT);
    244244 #else
    … …  
    256256      treeReader->ReadEntry(entry);
    257257      treeWriter->Clear();
    258      
    259258      if((entry % 100) == 0 && entry > 0 )  cout << "** Processing element # " << entry << endl;
    260259     
    … …  
    329328              } // switch (pid)
    330329             
    331               if(pid != pMU)
     330              if(pid != pMU && genMomentum.Et() !=0)
    332331                {
    333332                  towers.push_back(PhysicsTower(LorentzVector(genMomentum.Px(),genMomentum.Py(),genMomentum.Pz(), genMomentum.E())));
    … …  
    354353   
    355354     TLorentzVector Att(0.,0.,0.,0.);
    356      float ScalarEt=0;
    357355     for(unsigned int i=0; i < towers.size(); i++)
    358356       {
    359357         Att.SetPxPyPzE(towers[i].fourVector.px,towers[i].fourVector.py,towers[i].fourVector.pz,towers[i].fourVector.E);
    360          ScalarEt = ScalarEt + Att.Et();
    361358         PTmisS = PTmisS + Att;
    362359       }
    363 //cout<<"non smeare "<<PTmis.Pt()<<" "<<PTmisS.Pt()<<endl;
    364 //cout<<"smeare "<<PTmis.Px()<<" "<<PTmisS.Px()<<endl;
    365 //cout<<"**********"<<endl;
    366360
    367361      elementEtmis= (ETMIS*) branchetmis->NewEntry();
    368362      elementEtmis->Et = (PTmis).Pt();
    369       elementEtmis->Ex = (-PTmis).Px();
    370       elementEtmis->SEt = ScalarEt;
    371 
    372363      elementEtmis->EtSmeare = (PTmisS).Pt()-(PTmis).Pt();
    373       elementEtmis->ExSmeare = (-PTmisS).Px()-(PTmis).Px();
    374364     
    375365      //*****************************

  • trunk/interface/SmearUtil.h

    r33 r43  
    9393  double CONERADIUS;
    9494  int JETALGO;
     95
     96  //General jet variable
     97  double SEEDTHRESHOLD;   
     98  double OVERLAPTHRESHOLD;
     99
    95100  // MidPoint algorithm definition
    96   double M_SEEDTHRESHOLD;   
    97101  double M_CONEAREAFRACTION;
    98102  int    M_MAXPAIRSIZE;     
    99103  int    M_MAXITERATIONS;   
    100   double M_OVERLAPTHRESHOLD;
    101104 
    102105  // Define Cone algorithm.
    103   double C_SEEDTHRESHOLD;
    104106  int    C_ADJACENCYCUT;
    105107  int    C_MAXITERATIONS;
    106108  int    C_IRATCH;
    107   double C_OVERLAPTHRESHOLD;
    108109
    109110  /// Reads the data card for the initialisation of the parameters

  • trunk/src/SmearUtil.cc

    r33 r43  
    8080CONERADIUS   = 0.7;     // generic jet radius ; not for tau's !!!
    8181JETALGO      =   1;     // 1 for Cone algorithm, 2 for MidPoint algorithm, 3 for SIScone algorithm, 4 for kt algorithm
     82
     83//General jet parameters
     84SEEDTHRESHOLD     = 1.0;
     85OVERLAPTHRESHOLD = 0.75;
     86
    8287// Define Cone algorithm.
    83 C_SEEDTHRESHOLD     = 1.0;
    8488C_ADJACENCYCUT     = 2;
    8589C_MAXITERATIONS    = 100;
    8690C_IRATCH           = 1;
    87 C_OVERLAPTHRESHOLD = 0.75;
    8891
    8992//Define MidPoint algorithm.
    90 M_SEEDTHRESHOLD     =   1.0;
    9193M_CONEAREAFRACTION  =   0.25;
    9294M_MAXPAIRSIZE       =   2;
    9395M_MAXITERATIONS     =   100;
    94 M_OVERLAPTHRESHOLD  =  0.75;
    9596
    9697}
    … …  
    149150  }
    150151
     152// General jet variables
     153 SEEDTHRESHOLD     = 1.0;
     154 OVERLAPTHRESHOLD = 0.75;
     155
    151156// Define Cone algorithm.
    152  C_SEEDTHRESHOLD     = 1.0;
    153157 C_ADJACENCYCUT     = 2;
    154158 C_MAXITERATIONS    = 100;
    155159 C_IRATCH           = 1;
    156  C_OVERLAPTHRESHOLD = 0.75;
    157160
    158161//Define MidPoint algorithm.
    159  M_SEEDTHRESHOLD     =   1.0;
    160162 M_CONEAREAFRACTION  =   0.25;
    161163 M_MAXPAIRSIZE       =   2;
    162164 M_MAXITERATIONS     =   100;
    163  M_OVERLAPTHRESHOLD  =  0.75;
    164165
    165166}
    … …  
    250251  double Energie=0;
    251252  for(unsigned int i=0; i < towers.size(); i++) {
    252       if(towers[i].fourVector.pt() < M_SEEDTHRESHOLD) continue;
     253      if(towers[i].fourVector.pt() < SEEDTHRESHOLD) continue;
    253254      if((DeltaR(phi,eta,towers[i].fourVector.phi(),towers[i].fourVector.eta()) < TAU_CONE_ENERGY)) {
    254255          Energie += towers[i].fourVector.E;
Note: See TracChangeset for help on using the changeset viewer.