Fork me on GitHub

Changeset e9c0d73 in git


Ignore:
Timestamp:
Dec 5, 2017, 1:11:19 PM (7 years ago)
Author:
Ulrike Schnoor <schnooru@…>
Branches:
ImprovedOutputFile, Timing, dual_readout, llp, master
Children:
0879ed1
Parents:
197fed7
Message:

added exclusive_ymerge for n=2,3,4,5

Files:
4 edited

Legend:

Unmodified
Added
Removed
  • classes/DelphesClasses.cc

    r197fed7 re9c0d73  
    157157  NSubJetsPruned(0),
    158158  NSubJetsSoftDropped(0),
     159  ExclYmerge23(0),
     160  ExclYmerge34(0),
     161  ExclYmerge45(0),
     162  ExclYmerge56(0),
    159163  fFactory(0),
    160164  fArray(0)
  • classes/DelphesClasses.h

    r197fed7 re9c0d73  
    392392  TLorentzVector SoftDroppedP4[5]; // first entry (i = 0) is the total SoftDropped Jet 4-momenta and from i = 1 to 4 are the pruned subjets 4-momenta
    393393
     394
    394395  Int_t NSubJetsTrimmed; // number of subjets trimmed
    395396  Int_t NSubJetsPruned; // number of subjets pruned
    396397  Int_t NSubJetsSoftDropped; // number of subjets soft-dropped
    397398
     399  Double_t ExclYmerge23;
     400  Double_t ExclYmerge34;
     401  Double_t ExclYmerge45;
     402  Double_t ExclYmerge56;
     403 
    398404  TRefArray Constituents; // references to constituents
    399405  TRefArray Particles; // references to generated particles
     
    633639  Int_t NSubJetsSoftDropped; // number of subjets soft-dropped
    634640
    635 
     641  // Exclusive clustering variables
     642  Double_t ExclYmerge23;
     643  Double_t ExclYmerge34;
     644  Double_t ExclYmerge45;
     645  Double_t ExclYmerge56;
     646     
    636647  static CompBase *fgCompare; //!
    637648  const CompBase *GetCompare() const { return fgCompare; }
  • modules/FastJetFinder.cc

    r197fed7 re9c0d73  
    8080
    8181FastJetFinder::FastJetFinder() :
    82   fPlugin(0), fRecomb(0), fAxesDef(0), fMeasureDef(0), fNjettinessPlugin(0),
    83   fDefinition(0), fAreaDefinition(0), fItInputArray(0), fValenciaPlugin(0)
     82  fPlugin(0), fRecomb(0), fAxesDef(0), fMeasureDef(0), fNjettinessPlugin(0), fValenciaPlugin(0),
     83  fDefinition(0), fAreaDefinition(0), fItInputArray(0)
    8484{
    8585
     
    332332  vector< PseudoJet >::iterator itInputList, itOutputList;
    333333  vector< TEstimatorStruct >::iterator itEstimators;
    334 
     334  Double_t excl_ymerge23 = 0.0;
     335  Double_t excl_ymerge34 = 0.0;
     336  Double_t excl_ymerge45 = 0.0;
     337  Double_t excl_ymerge56 = 0.0;
     338 
    335339  DelphesFactory *factory = GetFactory();
    336340
     
    377381  outputList.clear();
    378382
     383 
    379384  if(fExclusiveClustering)
    380385    {
    381       //not neccessary to sort
    382       outputList = sequence->exclusive_jets( fNJets );
     386      outputList = sorted_by_pt(sequence->exclusive_jets( fNJets ));
     387
     388      excl_ymerge23 = sequence->exclusive_ymerge( 2 );
     389      excl_ymerge34 = sequence->exclusive_ymerge( 3 );
     390      excl_ymerge45 = sequence->exclusive_ymerge( 4 );
     391      excl_ymerge56 = sequence->exclusive_ymerge( 5 );
    383392    }
    384393  else
     
    401410    if(fAreaDefinition) area = itOutputList->area_4vector();
    402411
     412
     413   
    403414    candidate = factory->NewCandidate();
    404415
     
    444455    candidate->NNeutrals = nneutrals;
    445456    candidate->NCharged = ncharged;
     457
     458
     459    //for exclusive clustering, access y_n,n+1 as exclusive_ymerge (fNJets);
     460    candidate->ExclYmerge23 = excl_ymerge23;
     461    candidate->ExclYmerge34 = excl_ymerge34;
     462    candidate->ExclYmerge45 = excl_ymerge45;
     463    candidate->ExclYmerge56 = excl_ymerge56;
    446464   
    447465    //------------------------------------
  • modules/TreeWriter.cc

    r197fed7 re9c0d73  
    690690    }
    691691
     692    //--- exclusive clustering variables ---
     693    entry->ExclYmerge23 = candidate->ExclYmerge23;
     694    entry->ExclYmerge34 = candidate->ExclYmerge34;
     695    entry->ExclYmerge45 = candidate->ExclYmerge45;
     696    entry->ExclYmerge56 = candidate->ExclYmerge56;   
     697
     698
    692699    FillParticles(candidate, &entry->Particles);
    693700  }
Note: See TracChangeset for help on using the changeset viewer.