← Back to branch summary

~ma5/madanalysis5/madanalysis-development

« back to all changes in this revision

Viewing changes to tools/SampleAnalyzer/Commons/Service/TransverseVariables.cpp

  • Committer: Benjamin Fuks
  • Date: 2018-05-04 10:44:49 UTC
  • mfrom: (115.1.81 v1.6beta)
  • Revision ID: fuks@cern.ch-20180504104449-60h8a00loxgr8zg0
Releasing v1.6

Show diffs side-by-side

added added

removed removed

Lines of Context:
tools/SampleAnalyzer/Commons/Service/TransverseVariables.cpp
23
23
 
24
24
 
25
25
// STL headers
 
26
#include <cmath>
 
27
 
26
28
// SampleAnalyzer headers
27
29
#include "SampleAnalyzer/Commons/Service/TransverseVariables.h"
28
30
#include "SampleAnalyzer/Commons/Service/LogService.h"
524
526
     double THT1 = 0; double THT2 = 0;
525
527
     for (unsigned int i=0;i<jets1.size();i++) THT1+=jets1[i]->et();
526
528
     for (unsigned int i=0;i<jets2.size();i++) THT2+=jets2[i]->et();
527
 
     double DeltaHT = fabs(THT1-THT2);
 
529
     double DeltaHT = std::abs(THT1-THT2);
528
530
     if (DeltaHT<MinDHT) MinDHT=DeltaHT;
529
531
 
530
532
    // Exit
560
562
     double THT1 = 0; double THT2 = 0;
561
563
     for (unsigned int i=0;i<jets1.size();i++) THT1+=jets1[i].et();
562
564
     for (unsigned int i=0;i<jets2.size();i++) THT2+=jets2[i].et();
563
 
     double DeltaHT = fabs(THT1-THT2);
 
565
     double DeltaHT = std::abs(THT1-THT2);
564
566
     if (DeltaHT<MinDHT) MinDHT=DeltaHT;
565
567
 
566
568
    // Exit
595
597
 
596
598
  // dijet event
597
599
  if (jets.size()==2) return std::min(jets[0]->et(),jets[1]->et()) /
598
 
                            ( jets[0]->momentum()+jets[1]->momentum() ).Mt();
 
600
                             ( jets[0]->momentum()+jets[1]->momentum() ).Mt();
599
601
 
600
602
  double MinDeltaHT = 1e6;
601
603