SuperDijetMela Class Reference

#include <SuperDijetMela.h>

Public Member Functions
	SuperDijetMela (float sqrts_, TVar::VerbosityLevel verbosity_)
	SuperDijetMela (const SuperDijetMela &other)
	~SuperDijetMela ()
void	SetVerbosity (const TVar::VerbosityLevel verbosity_)
void	SetupResolutionModel (TVar::Production prod)
float	GetConvBW (TVar::Production prod, MELACandidate *cand, bool useTrueBW)

Protected Member Functions
void	Build ()

Protected Attributes
float	sqrts
TVar::VerbosityLevel	verbosity
std::unordered_map< int, MELADifermionResolutionModel * >	ResolutionModelMap

Detailed Description

Definition at line 8 of file SuperDijetMela.h.

Constructor & Destructor Documentation

SuperDijetMela() [1/2]

SuperDijetMela::SuperDijetMela	(	float	sqrts_,
		TVar::VerbosityLevel	verbosity_
	)

Definition at line 14 of file SuperDijetMela.cc.

                                                                          :
sqrts(sqrts_),
verbosity(verbosity_)
{
  Build();
}

SuperDijetMela() [2/2]

SuperDijetMela::SuperDijetMela

(

const SuperDijetMela &

other

)

Definition at line 20 of file SuperDijetMela.cc.

                                                          : 
sqrts(other.sqrts),
verbosity(other.verbosity)
{
  Build();
}

~SuperDijetMela()

SuperDijetMela::~SuperDijetMela

(

)

Definition at line 26 of file SuperDijetMela.cc.

                               {
  for (auto it=ResolutionModelMap.begin(); it!=ResolutionModelMap.end(); it++) delete it->second;
}

Member Function Documentation

Build()

void SuperDijetMela::Build ( )

protected

Definition at line 30 of file SuperDijetMela.cc.

                          {
  SetupResolutionModel(TVar::Had_ZH);
  SetupResolutionModel(TVar::Had_WH);
}

GetConvBW()

float SuperDijetMela::GetConvBW	(	TVar::Production	prod,
		MELACandidate *	cand,
		bool	useTrueBW
	)

Definition at line 64 of file SuperDijetMela.cc.

                                                                                       {
  float result=-1;
  int iprod = (int)prod;
  if (cand!=0){
    float mJJval=-1;
    bool modelExists = (ResolutionModelMap.find(iprod)!=ResolutionModelMap.end());
    int AssociationVCompatibility=0;
 
    if (modelExists || useTrueBW){
      // Get a simple event record, safest way to handle jet mass corrections
      int nRequested_AssociatedJets=0;
      int partIncCode=TVar::kNoAssociated; // Just to avoid warnings
      if (prod == TVar::Had_ZH || prod == TVar::Had_WH){ // Only use associated partons
        partIncCode=TVar::kUseAssociated_Jets;
        nRequested_AssociatedJets=2;
      }
      if (prod==TVar::Had_WH) AssociationVCompatibility=24;
      else if (prod==TVar::Had_ZH) AssociationVCompatibility=23;
      simple_event_record mela_event;
      mela_event.AssociationCode=partIncCode;
      mela_event.AssociationVCompatibility=AssociationVCompatibility;
      mela_event.nRequested_AssociatedJets=nRequested_AssociatedJets;
      TUtil::GetBoostedParticleVectors(
        cand,
        mela_event,
        verbosity
        );
 
      vector<TLorentzVector> pJets;
      const SimpleParticleCollection_t& pAssociated = mela_event.pAssociated;
      for (auto& part : pAssociated){
        if (PDGHelpers::isAJet(part.first)) pJets.push_back(part.second);
        if (pJets.size()==2) break;
      }
      if (pJets.size()==2) mJJval = (pJets[0] + pJets[1]).M();
    }
 
    if (useTrueBW){
      float truePoleMass=-1;
      float truePoleWidth=-1;
 
      if (AssociationVCompatibility!=0){
        truePoleMass=TUtil::GetMass(AssociationVCompatibility);
        truePoleWidth=TUtil::GetDecayWidth(AssociationVCompatibility);
        if (mJJval>=0.) result = 1./(pow(pow(mJJval, 2)-pow(truePoleMass, 2), 2) + pow(truePoleMass*truePoleWidth, 2));
      }
 
      if (verbosity>=TVar::DEBUG) MELAout
        << "SuperDijetMela::GetConvBW[idV=" << AssociationVCompatibility << "]::"
        << "sqrt(s) = " << mJJval
        << ", true (m,Gamma) = ( " << truePoleMass << " , " << truePoleWidth << " )"
        << ", ideal BW = " << result
        << endl;
    }
    else if (modelExists){
      result = ResolutionModelMap[iprod]->getVal(mJJval);
      if (verbosity>=TVar::DEBUG) MELAout
        << "SuperDijetMela::GetConvBW[" << TVar::ProductionName(prod) << "]::"
        << "sqrt(s) = " << mJJval
        << ", reco BW = " << result
        << endl;
    }
    else{
      if (verbosity>=TVar::DEBUG) MELAout
        << "SuperDijetMela::GetConvBW:: No method for " << TVar::ProductionName(prod) << " is known. "
        << "sqrt(s) = " << mJJval
        << ", return value = " << result
        << endl;
    }
  }
  return result;
}

SetupResolutionModel()

void SuperDijetMela::SetupResolutionModel

(

TVar::Production

prod

)

Definition at line 35 of file SuperDijetMela.cc.

                                                            {
  // Setup file path
  const string MELAPKGPATH = TVar::GetMELAPath();
  TString path = TString(MELAPKGPATH.c_str()) + "data/resolution_mJJ_recoVStrue_";
  TString prodName;
  switch (prod){
  case TVar::Had_ZH:
    prodName = "ZH";
    break;
  case TVar::Had_WH:
    prodName = "WH";
    break;
  default:
    MELAout << "SuperDijetMela::SetupResolutionModel: Production " << TVar::ProductionName(prod) << " is unknown." << endl;
    return;
  }
  sqrts = floor(sqrts); // interim solution until 13.6 TeV production is ready
  
  path += prodName;
  path += Form("_%.0fTeV%s", sqrts, ".root");
 
  TString appendName = Form("mJJReso_%.0fTeV", sqrts);
  MELADifermionResolutionModel* model = new MELADifermionResolutionModel(prod, sqrts, path, appendName);
  if (model->isValid()){
    int iprod = (int)prod;
    ResolutionModelMap[iprod] = model;
  }
  else MELAerr << "SuperDijetMela::SetupResolutionModel: Model MELADifermionResolutionModel(" << TVar::ProductionName(prod) << ", " << sqrts << ", " << path << ", " << appendName << ") cannot be built." << endl;
}

SetVerbosity()

void SuperDijetMela::SetVerbosity ( const TVar::VerbosityLevel  verbosity_ )

inline

Definition at line 23 of file SuperDijetMela.h.

{ verbosity=verbosity_; }

Member Data Documentation

ResolutionModelMap

std::unordered_map<int, MELADifermionResolutionModel*> SuperDijetMela::ResolutionModelMap

protected

Definition at line 13 of file SuperDijetMela.h.

sqrts

float SuperDijetMela::sqrts

protected

Definition at line 10 of file SuperDijetMela.h.

verbosity

TVar::VerbosityLevel SuperDijetMela::verbosity

protected

Definition at line 11 of file SuperDijetMela.h.

---

The documentation for this class was generated from the following files:

• MELA/interface/SuperDijetMela.h
• MELA/src/SuperDijetMela.cc