AliEmcalJetShapeProperties.h

Go to the documentation of this file.

#ifndef ALIEMCALJETSHAPEPROPERTIES_H
#define ALIEMCALJETSHAPEPROPERTIES_H

/* Copyright(c) 1998-2016, ALICE Experiment at CERN, All rights reserved. *
 * See cxx source for full Copyright notice                               */

#include <TArrayF.h>
#include <Rtypes.h>
#include <TString.h>

class AliEmcalJetShapeProperties {

public:

  AliEmcalJetShapeProperties();
  AliEmcalJetShapeProperties(const AliEmcalJetShapeProperties &jetshape);
  AliEmcalJetShapeProperties& operator=(const AliEmcalJetShapeProperties &jetshape);

  //jet shape derivatives
  //jet mass
  void              SetFirstDerivative(Double_t d)                             { fJetShapeMassFirstDer = d                       ; }
  void              SetSecondDerivative(Double_t d)                            { fJetShapeMassSecondDer = d                      ; }
  void              SetFirstOrderSubtracted(Double_t d)                        { fJetShapeMassFirstSub = d                       ; }
  void              SetSecondOrderSubtracted(Double_t d)                       { fJetShapeMassSecondSub = d                      ; }
  Double_t          GetFirstDerivative()                                 const { return fJetShapeMassFirstDer                    ; }
  Double_t          GetSecondDerivative()                                const { return fJetShapeMassSecondDer                   ; }
  Double_t          GetFirstOrderSubtracted()                            const { return fJetShapeMassFirstSub                    ; }
  Double_t          GetSecondOrderSubtracted()                           const { return fJetShapeMassSecondSub                   ; }

  //jet structure function
  TArrayF           GetGRNumerator()                                     const { return fGRNumerator                             ; }
  TArrayF           GetGRDenominator()                                   const { return fGRDenominator                           ; }
  TArrayF           GetGRNumeratorSub()                                  const { return fGRNumeratorSub                          ; }
  TArrayF           GetGRDenominatorSub()                                const { return fGRDenominatorSub                        ; }
  void              AddGRNumAt(Float_t num, Int_t idx)                         { fGRNumerator.AddAt(num, idx)                    ; }
  void              AddGRDenAt(Float_t den, Int_t idx)                         { fGRDenominator.AddAt(den, idx)                  ; }
  void              SetGRNumSize(UInt_t s)                                     { fGRNumerator.Set(s)                             ; }
  void              SetGRDenSize(UInt_t s)                                     { fGRDenominator.Set(s)                           ; }

  void              AddGRNumSubAt(Float_t num, Int_t idx)                      { fGRNumeratorSub.AddAt(num, idx)                 ; }
  void              AddGRDenSubAt(Float_t den, Int_t idx)                      { fGRDenominatorSub.AddAt(den, idx)               ; }
  void              SetGRNumSubSize(UInt_t s)                                  { fGRNumeratorSub.Set(s)                          ; }
  void              SetGRDenSubSize(UInt_t s)                                  { fGRDenominatorSub.Set(s)                        ; }
  void              PrintGR();

  //Angularity
  void              SetFirstDerivativeAngularity(Double_t d)                   { fJetShapeAngularityFirstDer = d                 ; }
  void              SetSecondDerivativeAngularity(Double_t d)                  { fJetShapeAngularitySecondDer = d                ; }
  void              SetFirstOrderSubtractedAngularity(Double_t d)              { fJetShapeAngularityFirstSub = d                 ; }
  void              SetSecondOrderSubtractedAngularity(Double_t d)             { fJetShapeAngularitySecondSub = d                ; }
  Double_t          GetFirstDerivativeAngularity()                       const { return fJetShapeAngularityFirstDer              ; }
  Double_t          GetSecondDerivativeAngularity()                      const { return fJetShapeAngularitySecondDer             ; }
  Double_t          GetFirstOrderSubtractedAngularity()                  const { return fJetShapeAngularityFirstSub              ; }
  Double_t          GetSecondOrderSubtractedAngularity()                 const { return fJetShapeAngularitySecondSub             ; }

  //pTD
  void              SetFirstDerivativepTD(Double_t d)                          { fJetShapepTDFirstDer = d                        ; }
  void              SetSecondDerivativepTD(Double_t d)                         { fJetShapepTDSecondDer = d                       ; }
  void              SetFirstOrderSubtractedpTD(Double_t d)                     { fJetShapepTDFirstSub = d                        ; }
  void              SetSecondOrderSubtractedpTD(Double_t d)                    { fJetShapepTDSecondSub = d                       ; }
  Double_t          GetFirstDerivativepTD()                              const { return fJetShapepTDFirstDer                     ; }
  Double_t          GetSecondDerivativepTD()                             const { return fJetShapepTDSecondDer                    ; }
  Double_t          GetFirstOrderSubtractedpTD()                         const { return fJetShapepTDFirstSub                     ; }
  Double_t          GetSecondOrderSubtractedpTD()                        const { return fJetShapepTDSecondSub                    ; }

  //Circularity
  void              SetFirstDerivativeCircularity(Double_t d)                  { fJetShapeCircularityFirstDer = d                ; }
  void              SetSecondDerivativeCircularity(Double_t d)                 { fJetShapeCircularitySecondDer = d               ; }
  void              SetFirstOrderSubtractedCircularity(Double_t d)             { fJetShapeCircularityFirstSub = d                ; }
  void              SetSecondOrderSubtractedCircularity(Double_t d)            { fJetShapeCircularitySecondSub = d               ; }
  Double_t          GetFirstDerivativeCircularity()                      const { return fJetShapeCircularityFirstDer             ; }
  Double_t          GetSecondDerivativeCircularity()                     const { return fJetShapeCircularitySecondDer            ; }
  Double_t          GetFirstOrderSubtractedCircularity()                 const { return fJetShapeCircularityFirstSub             ; }
  Double_t          GetSecondOrderSubtractedCircularity()                const { return fJetShapeCircularitySecondSub            ; }

  //Sigma2
  void              SetFirstDerivativeSigma2(Double_t d)                       { fJetShapeSigma2FirstDer = d                     ; }
  void              SetSecondDerivativeSigma2(Double_t d)                      { fJetShapeSigma2SecondDer = d                    ; }
  void              SetFirstOrderSubtractedSigma2(Double_t d)                  { fJetShapeSigma2FirstSub = d                     ; }
  void              SetSecondOrderSubtractedSigma2(Double_t d)                 { fJetShapeSigma2SecondSub = d                    ; }
  Double_t          GetFirstDerivativeSigma2()                           const { return fJetShapeSigma2FirstDer                  ; }
  Double_t          GetSecondDerivativeSigma2()                          const { return fJetShapeSigma2SecondDer                 ; }
  Double_t          GetFirstOrderSubtractedSigma2()                      const { return fJetShapeSigma2FirstSub                  ; }
  Double_t          GetSecondOrderSubtractedSigma2()                     const { return fJetShapeSigma2SecondSub                 ; }


  //number of contituents
  void              SetFirstDerivativeConstituent(Double_t d)                  { fJetShapeConstituentFirstDer = d                ; }
  void              SetSecondDerivativeConstituent(Double_t d)                 { fJetShapeConstituentSecondDer = d               ; }
  void              SetFirstOrderSubtractedConstituent(Double_t d)             { fJetShapeConstituentFirstSub = d                ; }
  void              SetSecondOrderSubtractedConstituent(Double_t d)            { fJetShapeConstituentSecondSub = d               ; }
  Double_t          GetFirstDerivativeConstituent()                      const { return fJetShapeConstituentFirstDer             ; }
  Double_t          GetSecondDerivativeConstituent()                     const { return fJetShapeConstituentSecondDer            ; }
  Double_t          GetFirstOrderSubtractedConstituent()                 const { return fJetShapeConstituentFirstSub             ; }
  Double_t          GetSecondOrderSubtractedConstituent()                const { return fJetShapeConstituentSecondSub            ; }

  //leading minus subleading constituent
  void              SetFirstDerivativeLeSub(Double_t d)                        { fJetShapeLeSubFirstDer = d                      ; }
  void              SetSecondDerivativeLeSub(Double_t d)                       { fJetShapeLeSubSecondDer = d                     ; }
  void              SetFirstOrderSubtractedLeSub(Double_t d)                   { fJetShapeLeSubFirstSub = d                      ; }
  void              SetSecondOrderSubtractedLeSub(Double_t d)                  { fJetShapeLeSubSecondSub = d                     ; }
  Double_t          GetFirstDerivativeLeSub()                            const { return fJetShapeLeSubFirstDer                   ; }
  Double_t          GetSecondDerivativeLeSub()                           const { return fJetShapeLeSubSecondDer                  ; }
  Double_t          GetFirstOrderSubtractedLeSub()                       const { return fJetShapeLeSubFirstSub                   ; }
  Double_t          GetSecondOrderSubtractedLeSub()                      const { return fJetShapeLeSubSecondSub                  ; }

  //1subjettiness_kt
  void              SetFirstDerivative1subjettiness_kt(Double_t d)             { fJetShape1subjettinessktFirstDer = d                      ; }
  void              SetSecondDerivative1subjettiness_kt(Double_t d)            { fJetShape1subjettinessktSecondDer = d                     ; }
  void              SetFirstOrderSubtracted1subjettiness_kt(Double_t d)        { fJetShape1subjettinessktFirstSub = d                      ; }
  void              SetSecondOrderSubtracted1subjettiness_kt(Double_t d)       { fJetShape1subjettinessktSecondSub = d                     ; }
  Double_t          GetFirstDerivative1subjettiness_kt()                 const { return fJetShape1subjettinessktFirstDer                   ; }
  Double_t          GetSecondDerivative1subjettiness_kt()                const { return fJetShape1subjettinessktSecondDer                  ; }
  Double_t          GetFirstOrderSubtracted1subjettiness_kt()            const { return fJetShape1subjettinessktFirstSub                   ; }
  Double_t          GetSecondOrderSubtracted1subjettiness_kt()           const { return fJetShape1subjettinessktSecondSub                  ; }

  //2subjettiness_kt
  void              SetFirstDerivative2subjettiness_kt(Double_t d)             { fJetShape2subjettinessktFirstDer = d                      ; }
  void              SetSecondDerivative2subjettiness_kt(Double_t d)            { fJetShape2subjettinessktSecondDer = d                     ; }
  void              SetFirstOrderSubtracted2subjettiness_kt(Double_t d)        { fJetShape2subjettinessktFirstSub = d                      ; }
  void              SetSecondOrderSubtracted2subjettiness_kt(Double_t d)       { fJetShape2subjettinessktSecondSub = d                     ; }
  Double_t          GetFirstDerivative2subjettiness_kt()                 const { return fJetShape2subjettinessktFirstDer                   ; }
  Double_t          GetSecondDerivative2subjettiness_kt()                const { return fJetShape2subjettinessktSecondDer                  ; }
  Double_t          GetFirstOrderSubtracted2subjettiness_kt()            const { return fJetShape2subjettinessktFirstSub                   ; }
  Double_t          GetSecondOrderSubtracted2subjettiness_kt()           const { return fJetShape2subjettinessktSecondSub                  ; }

  //3subjettiness_kt
  void              SetFirstDerivative3subjettiness_kt(Double_t d)             { fJetShape3subjettinessktFirstDer = d                      ; }
  void              SetSecondDerivative3subjettiness_kt(Double_t d)            { fJetShape3subjettinessktSecondDer = d                     ; }
  void              SetFirstOrderSubtracted3subjettiness_kt(Double_t d)        { fJetShape3subjettinessktFirstSub = d                      ; }
  void              SetSecondOrderSubtracted3subjettiness_kt(Double_t d)       { fJetShape3subjettinessktSecondSub = d                     ; }
  Double_t          GetFirstDerivative3subjettiness_kt()                 const { return fJetShape3subjettinessktFirstDer                   ; }
  Double_t          GetSecondDerivative3subjettiness_kt()                const { return fJetShape3subjettinessktSecondDer                  ; }
  Double_t          GetFirstOrderSubtracted3subjettiness_kt()            const { return fJetShape3subjettinessktFirstSub                   ; }
  Double_t          GetSecondOrderSubtracted3subjettiness_kt()           const { return fJetShape3subjettinessktSecondSub                  ; }

  //OpeningAngle_kt
  void              SetFirstDerivativeOpeningAngle_kt(Double_t d)             { fJetShapeOpeningAnglektFirstDer = d                      ; }
  void              SetSecondDerivativeOpeningAngle_kt(Double_t d)            { fJetShapeOpeningAnglektSecondDer = d                     ; }
  void              SetFirstOrderSubtractedOpeningAngle_kt(Double_t d)        { fJetShapeOpeningAnglektFirstSub = d                      ; }
  void              SetSecondOrderSubtractedOpeningAngle_kt(Double_t d)       { fJetShapeOpeningAnglektSecondSub = d                     ; }
  Double_t          GetFirstDerivativeOpeningAngle_kt()                 const { return fJetShapeOpeningAnglektFirstDer                   ; }
  Double_t          GetSecondDerivativeOpeningAngle_kt()                const { return fJetShapeOpeningAnglektSecondDer                  ; }
  Double_t          GetFirstOrderSubtractedOpeningAngle_kt()            const { return fJetShapeOpeningAnglektFirstSub                   ; }
  Double_t          GetSecondOrderSubtractedOpeningAngle_kt()           const { return fJetShapeOpeningAnglektSecondSub                  ; }

  //SoftDrop
  void              SetSoftDropZg(Double_t d)                                 { fSoftDropZg = d                      ; }
  void              SetSoftDropdR(Double_t d)                                 { fSoftDropdR = d                      ; }
  void              SetSoftDropPtfrac(Double_t d)                             { fSoftDropPtfrac = d                  ; }
  void              SetSoftDropDropCount(Int_t d)                             { fSoftDropDropCount = d               ; }
  Double_t          GetSoftDropZg()                                     const { return fSoftDropZg                   ; }
  Double_t          GetSoftDropdR()                                     const { return fSoftDropdR                   ; }
  Double_t          GetSoftDropPtfrac()                                 const { return fSoftDropPtfrac               ; }
  Int_t             GetSoftDropDropCount()                              const { return fSoftDropDropCount;           ; }

  void              PrintGR() const;

protected:

  Double_t          fJetShapeMassFirstDer;                   
  Double_t          fJetShapeMassSecondDer;                  
  Double_t          fJetShapeMassFirstSub;                   
  Double_t          fJetShapeMassSecondSub;                  

  TArrayF           fGRNumerator;                            
  TArrayF           fGRDenominator;                          
  TArrayF           fGRNumeratorSub;                         
  TArrayF           fGRDenominatorSub;                       

  Double_t          fJetShapeAngularityFirstDer;             
  Double_t          fJetShapeAngularitySecondDer;            
  Double_t          fJetShapeAngularityFirstSub;             
  Double_t          fJetShapeAngularitySecondSub;            

  Double_t          fJetShapepTDFirstDer;                    
  Double_t          fJetShapepTDSecondDer;                   
  Double_t          fJetShapepTDFirstSub;                    
  Double_t          fJetShapepTDSecondSub;                   

  Double_t          fJetShapeCircularityFirstDer;            
  Double_t          fJetShapeCircularitySecondDer;           
  Double_t          fJetShapeCircularityFirstSub;            
  Double_t          fJetShapeCircularitySecondSub;           

  Double_t          fJetShapeSigma2FirstDer;                 
  Double_t          fJetShapeSigma2SecondDer;                
  Double_t          fJetShapeSigma2FirstSub;                 
  Double_t          fJetShapeSigma2SecondSub;                

  Double_t          fJetShapeConstituentFirstDer;            
  Double_t          fJetShapeConstituentSecondDer;           
  Double_t          fJetShapeConstituentFirstSub;            
  Double_t          fJetShapeConstituentSecondSub;           

  Double_t          fJetShapeLeSubFirstDer;                  
  Double_t          fJetShapeLeSubSecondDer;                 
  Double_t          fJetShapeLeSubFirstSub;                  
  Double_t          fJetShapeLeSubSecondSub;                 

  Double_t          fJetShape1subjettinessktFirstDer;        
  Double_t          fJetShape1subjettinessktSecondDer;       
  Double_t          fJetShape1subjettinessktFirstSub;        
  Double_t          fJetShape1subjettinessktSecondSub;       
  
  Double_t          fJetShape2subjettinessktFirstDer;        
  Double_t          fJetShape2subjettinessktSecondDer;       
  Double_t          fJetShape2subjettinessktFirstSub;        
  Double_t          fJetShape2subjettinessktSecondSub;       

  Double_t          fJetShape3subjettinessktFirstDer;        
  Double_t          fJetShape3subjettinessktSecondDer;       
  Double_t          fJetShape3subjettinessktFirstSub;        
  Double_t          fJetShape3subjettinessktSecondSub;       

  Double_t          fJetShapeOpeningAnglektFirstDer;         
  Double_t          fJetShapeOpeningAnglektSecondDer;        
  Double_t          fJetShapeOpeningAnglektFirstSub;         
  Double_t          fJetShapeOpeningAnglektSecondSub;        

  Double_t          fSoftDropZg;                             
  Double_t          fSoftDropdR;                             
  Double_t          fSoftDropPtfrac;                         
  Int_t             fSoftDropDropCount;                      
};

#endif