[u/mrichter/AliRoot.git] / PWG3 / dielectron / AliDielectronCF.h

#ifndef ALIDIELECTRONCF_H
#define ALIDIELECTRONCF_H
/* Copyright(c) 1998-2009, ALICE Experiment at CERN, All rights reserved. *
 * See cxx source for full Copyright notice                               */

//#############################################################
//#                                                           #
//#             Class AliDielectronCF                         #
//#       Dielectron Correction Framework Manager             #
//#                                                           #
//#  Authors:                                                 #
//#   Anton     Andronic, GSI / A.Andronic@gsi.de             #
//#   Ionut C.  Arsene,   GSI / I.C.Arsene@gsi.de             #
//#   Julian    Book,     Uni Ffm / Julian.Book@cern.ch       #
//#   Frederick Kramer,   Uni Ffm, / Frederick.Kramer@cern.ch #
//#   Magnus    Mager,    CERN / Magnus.Mager@cern.ch         #
//#   WooJin J. Park,     GSI / W.J.Park@gsi.de               #
//#   Jens      Wiechula, Uni HD / Jens.Wiechula@cern.ch      #
//#                                                           #
//#############################################################


#include <TNamed.h>
#include <TVectorDfwd.h>
#include "AliDielectronVarManager.h"

class AliAnalysisCuts;
class AliAnalysisFilter;
class AliCFContainer;
class AliDielectronPair;
class TObjArray;

class AliDielectronCF : public TNamed {
public:
  enum {kNmaxAddSteps=50};
  
  AliDielectronCF();
  AliDielectronCF(const char* name, const char* title);
  virtual ~AliDielectronCF();

  void SetStepForMCtruth(Bool_t steps=kTRUE)           { fStepForMCtruth=steps;           }
  void SetStepForNoCutsMCmotherPid(Bool_t steps=kTRUE) { fStepForNoCutsMCmotherPid=steps; }
  void SetStepForAfterAllCuts(Bool_t steps=kTRUE)      { fStepForAfterAllCuts=steps;      }
  void SetStepsForEachCut(Bool_t steps=kTRUE)          { fStepsForEachCut=steps;          }
  void SetStepsForCutsIncreasing(Bool_t steps=kTRUE)   { fStepsForCutsIncreasing=steps;   }
  void SetStepsForSignal(Bool_t steps=kTRUE)           { fStepsForSignal=steps;           }
  void SetStepsForBackground(Bool_t steps=kTRUE)       { fStepsForBackground=steps;       }
  
  void SetPdgMother(Int_t pdg) { fPdgMother=pdg; }
  
  void AddStepMask(UInt_t mask)                  { fStepMasks[fNStepMasks++]=mask; }
  
  void AddVariable(AliDielectronVarManager::ValueTypes type, Int_t nbins, Double_t min, Double_t max, Bool_t leg=kFALSE);
  void AddVariable(AliDielectronVarManager::ValueTypes type, TVectorD * const binLimits);
  
  void InitialiseContainer(const AliAnalysisFilter& filter);
  
//   void Fill(UInt_t mask, const TObject *particle);
  void Fill(UInt_t mask, const AliDielectronPair *particle);
  void FillMC(const TObject *particle);
  
  AliCFContainer* GetContainer() const { return fCfContainer; }
  
private:
  UInt_t          fVariables[AliDielectronVarManager::kNMaxValues]; //configured variables
  UInt_t          fVariablesLeg[AliDielectronVarManager::kNMaxValues]; //configured variables for the legs
  
  Int_t           fNSteps;                     // number of selection steps
  
  Int_t           fNVars;                      // number of variables
  Int_t           fNBins[kNmaxAddSteps];       // array of numbers ob bins of the vars
  Double_t        fVarLoLimit[kNmaxAddSteps];  // array of the lower limits of the vars
  Double_t        fVarUpLimit[kNmaxAddSteps];  // array of the upper limits of the vars
  TObjArray      *fVarBinLimits;               // alternative array of bin limits

  Int_t           fNVarsLeg;                      // number of variables for the legs
  Int_t           fNBinsLeg[kNmaxAddSteps];       // array of numbers ob bins of the vars for the legs
  Double_t        fVarLoLimitLeg[kNmaxAddSteps];  // array of the lower limits of the vars for the legs
  Double_t        fVarUpLimitLeg[kNmaxAddSteps];  // array of the upper limits of the vars for the legs
  
  Int_t           fNCuts;                         // Number of cuts in the filter concerned

  Double_t        *fValues;                       //! Value array for filling the container
  
  Bool_t fStepForMCtruth;               //create a step for the MC truth
  Bool_t fStepForNoCutsMCmotherPid;     //create a step for before cuts, but with MC truth of the mother
  Bool_t fStepForAfterAllCuts;          //create a step for before cuts, but with MC truth of the mother
  Bool_t fStepsForEachCut;              //create steps for each cut?
  Bool_t fStepsForCutsIncreasing;       //create steps for increasing cut combinatons?
                                        //e.g. cut1&cut2, cut1&cut2&cut3 ...
  Bool_t fStepsForSignal;               //steps for pure signal
  Bool_t fStepsForBackground;           //steps for pure background
  
  UInt_t fStepMasks[kNmaxAddSteps];      //steps for additional cut combinatons
  UInt_t fNStepMasks;                    //number of configured step masks

  Int_t fPdgMother;                      //Pdg code of MCtruth validation
  AliCFContainer* fCfContainer;          //the CF container

  Bool_t fHasMC;                         //if MC info is available
  Int_t  fNAddSteps;                     //number of additional MC related steps per cut step
  
  AliDielectronCF(const AliDielectronCF &c);
  AliDielectronCF &operator=(const AliDielectronCF &c);
  
  ClassDef(AliDielectronCF,3)  //Dielectron Correction Framework handler
};

#endif
Commit	Line	Data
b2a297fa	1	#ifndef ALIDIELECTRONCF_H
	2	#define ALIDIELECTRONCF_H
	3	/* Copyright(c) 1998-2009, ALICE Experiment at CERN, All rights reserved. *
	4	* See cxx source for full Copyright notice */
	5
	6	//#############################################################
	7	//# #
	8	//# Class AliDielectronCF #
	9	//# Dielectron Correction Framework Manager #
	10	//# #
	11	//# Authors: #
	12	//# Anton Andronic, GSI / A.Andronic@gsi.de #
	13	//# Ionut C. Arsene, GSI / I.C.Arsene@gsi.de #
	14	//# Julian Book, Uni Ffm / Julian.Book@cern.ch #
	15	//# Frederick Kramer, Uni Ffm, / Frederick.Kramer@cern.ch #
	16	//# Magnus Mager, CERN / Magnus.Mager@cern.ch #
	17	//# WooJin J. Park, GSI / W.J.Park@gsi.de #
	18	//# Jens Wiechula, Uni HD / Jens.Wiechula@cern.ch #
	19	//# #
	20	//#############################################################
	21
	22
	23
	24	#include <TNamed.h>
61d106d3	25	#include <TVectorDfwd.h>
b2a297fa	26	#include "AliDielectronVarManager.h"
	27
	28	class AliAnalysisCuts;
	29	class AliAnalysisFilter;
	30	class AliCFContainer;
6551594b	31	class AliDielectronPair;
61d106d3	32	class TObjArray;
b2a297fa	33
	34	class AliDielectronCF : public TNamed {
	35	public:
	36	enum {kNmaxAddSteps=50};
	37
	38	AliDielectronCF();
	39	AliDielectronCF(const char* name, const char* title);
	40	virtual ~AliDielectronCF();
	41
e123f993	42	void SetStepForMCtruth(Bool_t steps=kTRUE) { fStepForMCtruth=steps; }
	43	void SetStepForNoCutsMCmotherPid(Bool_t steps=kTRUE) { fStepForNoCutsMCmotherPid=steps; }
	44	void SetStepForAfterAllCuts(Bool_t steps=kTRUE) { fStepForAfterAllCuts=steps; }
	45	void SetStepsForEachCut(Bool_t steps=kTRUE) { fStepsForEachCut=steps; }
	46	void SetStepsForCutsIncreasing(Bool_t steps=kTRUE) { fStepsForCutsIncreasing=steps; }
572b0139	47	void SetStepsForSignal(Bool_t steps=kTRUE) { fStepsForSignal=steps; }
572b0139	48	void SetStepsForBackground(Bool_t steps=kTRUE) { fStepsForBackground=steps; }
b2a297fa	49
	50	void SetPdgMother(Int_t pdg) { fPdgMother=pdg; }
	51
	52	void AddStepMask(UInt_t mask) { fStepMasks[fNStepMasks++]=mask; }
	53
6551594b	54	void AddVariable(AliDielectronVarManager::ValueTypes type, Int_t nbins, Double_t min, Double_t max, Bool_t leg=kFALSE);
61d106d3	55	void AddVariable(AliDielectronVarManager::ValueTypes type, TVectorD * const binLimits);
61d106d3	56
b2a297fa	57	void InitialiseContainer(const AliAnalysisFilter& filter);
b2a297fa	58
6551594b	59	// void Fill(UInt_t mask, const TObject *particle);
6551594b	60	void Fill(UInt_t mask, const AliDielectronPair *particle);
b2a297fa	61	void FillMC(const TObject *particle);
	62
	63	AliCFContainer* GetContainer() const { return fCfContainer; }
	64
6551594b	65	private:
b2a297fa	66	UInt_t fVariables[AliDielectronVarManager::kNMaxValues]; //configured variables
6551594b	67	UInt_t fVariablesLeg[AliDielectronVarManager::kNMaxValues]; //configured variables for the legs
b2a297fa	68
b2a297fa	69	Int_t fNSteps; // number of selection steps
6551594b	70
b2a297fa	71	Int_t fNVars; // number of variables
	72	Int_t fNBins[kNmaxAddSteps]; // array of numbers ob bins of the vars
	73	Double_t fVarLoLimit[kNmaxAddSteps]; // array of the lower limits of the vars
	74	Double_t fVarUpLimit[kNmaxAddSteps]; // array of the upper limits of the vars
61d106d3	75	TObjArray *fVarBinLimits; // alternative array of bin limits
b2a297fa	76
6551594b	77	Int_t fNVarsLeg; // number of variables for the legs
	78	Int_t fNBinsLeg[kNmaxAddSteps]; // array of numbers ob bins of the vars for the legs
	79	Double_t fVarLoLimitLeg[kNmaxAddSteps]; // array of the lower limits of the vars for the legs
	80	Double_t fVarUpLimitLeg[kNmaxAddSteps]; // array of the upper limits of the vars for the legs
	81
	82	Int_t fNCuts; // Number of cuts in the filter concerned
b2a297fa	83
572b0139	84	Double_t *fValues; //! Value array for filling the container
6551594b	85
e123f993	86	Bool_t fStepForMCtruth; //create a step for the MC truth
	87	Bool_t fStepForNoCutsMCmotherPid; //create a step for before cuts, but with MC truth of the mother
	88	Bool_t fStepForAfterAllCuts; //create a step for before cuts, but with MC truth of the mother
	89	Bool_t fStepsForEachCut; //create steps for each cut?
	90	Bool_t fStepsForCutsIncreasing; //create steps for increasing cut combinatons?
	91	//e.g. cut1&cut2, cut1&cut2&cut3 ...
572b0139	92	Bool_t fStepsForSignal; //steps for pure signal
	93	Bool_t fStepsForBackground; //steps for pure background
	94
b2a297fa	95	UInt_t fStepMasks[kNmaxAddSteps]; //steps for additional cut combinatons
	96	UInt_t fNStepMasks; //number of configured step masks
	97
	98	Int_t fPdgMother; //Pdg code of MCtruth validation
	99	AliCFContainer* fCfContainer; //the CF container
572b0139	100
	101	Bool_t fHasMC; //if MC info is available
	102	Int_t fNAddSteps; //number of additional MC related steps per cut step
b2a297fa	103
	104	AliDielectronCF(const AliDielectronCF &c);
	105	AliDielectronCF &operator=(const AliDielectronCF &c);
	106
6551594b	107	ClassDef(AliDielectronCF,3) //Dielectron Correction Framework handler
b2a297fa	108	};
	109
	110	#endif