AliAnalysisTaskEmcalJetHPerformance.h

Go to the documentation of this file.

#ifndef ALIANALYSISTASKEMCALJETHPERFORMANCE_H
#define ALIANALYSISTASKEMCALJETHPERFORMANCE_H

#include <string>

class AliJetContainer;
class AliEmcalJet;
#include "THistManager.h"
#include "AliYAMLConfiguration.h"
#include "AliAnalysisTaskEmcalJet.h"
#include "AliEmcalEmbeddingQA.h"
#include "AliAnalysisTaskEmcalJetHUtils.h"

// operator<< has to be forward declared carefully to stay in the global namespace so that it works with CINT.
// For generally how to keep the operator in the global namespace, See: https://stackoverflow.com/a/38801633
namespace PWGJE { namespace EMCALJetTasks { class AliAnalysisTaskEmcalJetHPerformance; } }
std::ostream & operator<< (std::ostream &in, const PWGJE::EMCALJetTasks::AliAnalysisTaskEmcalJetHPerformance &myTask);
void swap(PWGJE::EMCALJetTasks::AliAnalysisTaskEmcalJetHPerformance & first, PWGJE::EMCALJetTasks::AliAnalysisTaskEmcalJetHPerformance & second); 

namespace PWGJE {
namespace EMCALJetTasks {

class AliAnalysisTaskEmcalJetHPerformance : public AliAnalysisTaskEmcalJet {
 public:
  struct ResponseMatrixFillWrapper {
    double fPt;                 
    double fArea;               
    double fPhi;                
    double fRelativeEPAngle;    
    double fLeadingHadronPt;    
    double fDistance;           
    double fCentrality;         
  };

  AliAnalysisTaskEmcalJetHPerformance();
  AliAnalysisTaskEmcalJetHPerformance(const char * name);
  // Additional constructors
  AliAnalysisTaskEmcalJetHPerformance(const AliAnalysisTaskEmcalJetHPerformance & other);
  AliAnalysisTaskEmcalJetHPerformance& operator=(AliAnalysisTaskEmcalJetHPerformance other);
  friend void ::swap(AliAnalysisTaskEmcalJetHPerformance & first, AliAnalysisTaskEmcalJetHPerformance & second); 
  // Avoid implementing move since c++11 is not allowed in the header

  void UserCreateOutputObjects();

  // Initialize the task
  // Configuration is handled via the YAML configuration file
  bool Initialize();
  void SetConfigurationPath(const std::string & configurationPath) { fConfigurationPath = configurationPath; }

  // Utility functions
  // AddTask
  static AliAnalysisTaskEmcalJetHPerformance * AddTaskEmcalJetHPerformance(const char * suffix = "");
  // Printing
  std::string toString() const;
  friend std::ostream & ::operator<<(std::ostream &in, const AliAnalysisTaskEmcalJetHPerformance &myTask);
  void Print(Option_t* opt = "") const;
  std::ostream & Print(std::ostream &in) const;

 private:

  Bool_t Run();

  // Configuration
  void RetrieveAndSetTaskPropertiesFromYAMLConfig();
  void SetupJetContainersFromYAMLConfig();

  // Response matrix functions
  void SetupResponseMatrixHists();
  void CreateResponseMatrix();
  void FillResponseMatrix(AliEmcalJet * jet1, AliEmcalJet * jet2);
  ResponseMatrixFillWrapper CreateResponseMatrixFillWrapper(AliEmcalJet * jet) const;

  // Basic configuration
  PWG::Tools::AliYAMLConfiguration fYAMLConfig; 
  std::string fConfigurationPath;     
  bool fConfigurationInitialized;     

  // Histograms
  THistManager fHistManager;          
  AliEmcalEmbeddingQA fEmbeddingQA;   

  // Configuration options
  bool fCreateResponseMatrix;         

  // Response matrix variables
  // Response matrix fill map
  // For whatever reason, it appears that CINT can't handle this definition
  #if !(defined(__CINT__) || defined(__MAKECINT__))
  // For some reason, root can't handle this typedef
  //typedef double ResponseMatrixFillWrapper::*MP;
  //std::map <std::string, std::pair<int, MP>> fResponseMatrixFillMap;   //!<!
  std::map <std::string, std::pair<int, double ResponseMatrixFillWrapper::*>> fResponseMatrixFillMap;   
  #endif
  bool fResponseFromThreeJetCollections; 
  double fMinFractionShared;             
  AliAnalysisTaskEmcalJetHUtils::ELeadingHadronBiasType_t fLeadingHadronBiasType; 

  ClassDef(AliAnalysisTaskEmcalJetHPerformance, 1);
};

} /* namespace EMCALJetTasks */
} /* namespace PWGJE */

#endif /* AliAnalysisTaskEmcalJetHPerformance.h */