QCDTruthMultiJetFilter.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2022 CERN for the benefit of the ATLAS collaboration
*/
 
#ifndef GENERATORFILTERS_QCDTRUTHMULTIJETFILTER_H
#define GENERATORFILTERS_QCDTRUTHMULTIJETFILTER_H
 
#include "GeneratorModules/GenFilter.h"
#include "GaudiKernel/ServiceHandle.h"
#include "AthenaKernel/IAthRNGSvc.h"
#include <string>
 
namespace CLHEP {
  class HepRandomEngine;
}
 
class QCDTruthMultiJetFilter : public GenFilter {
public:
 
  QCDTruthMultiJetFilter(const std::string& name, ISvcLocator* pSvcLocator);
  StatusCode filterInitialize();
  StatusCode filterFinalize();
  StatusCode filterEvent();
 
private:
 
  CLHEP::HepRandomEngine* getRandomEngine(const std::string& streamName,
                                          const EventContext& ctx) const;
 
  int m_Njet;  
  double m_NjetMinPt;  
  double m_MinLeadJetPt;  
  double m_MaxLeadJetPt;  
  double m_MaxEta; 
 
  std::string m_TruthJetContainerName;  
 
  ServiceHandle<IAthRNGSvc> m_rndmSvc{this, "RndmSvc", "AthRNGSvc"};  
 
  long m_total;    
  long m_passed;   
  long m_ptfailed; 
  long m_nJetsFailed; 
 
  double m_norm;   
  double m_high;   
 
  bool m_doShape;  
 
public:
 
  static double fitFnR(const double * x , const double *) {
    static const double p[7] = //{100000,-5,0.0004,0.0000001,1./3600.,23,0.8};
      { 1000000, -5.5, 0.0001, 0.00000012, 1./4150., 23, 0.4 };
    return p[0]*std::pow(x[0],p[1]+p[2]*x[0]+p[3]*pow(x[0],2))*pow(1-x[0]*p[4],p[5])*pow(x[0],p[6]);
  }
 
  static double fitFn(const double x ) {
    static const double p[7] = //{100000,-5,0.0004,0.0000001,1./3600.,23,0.8};
      { 1000000, -5.5, 0.0001, 0.00000012, 1./4150., 23, 0.4 };
    return p[0]*std::pow(x,p[1]+p[2]*x+p[3]*pow(x,2))*pow(1-x*p[4],p[5])*pow(x,p[6]);
  }
 
};
 
#endif