[u/mrichter/AliRoot.git] / PWGLF / FORWARD / analysis2 / AliForwardCreateResponseMatrices.h

#ifndef ALIFORWARDCREATERESPONSEMATRICES_H
#define ALIFORWARDCREATERESPONSEMATRICES_H
 
/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
 * See cxx source for full Copyright notice                               */
  
#include "AliAnalysisTaskSE.h"
#include "AliBasedNdetaTask.h"
#include <TList.h>
#include <iostream>

class TH2D;

/**
 * Task to do the multiplicity distibution
 * 
 */
class AliForwardCreateResponseMatrices : public AliBasedNdetaTask
{
public:
  /** 
   * 
   * Default Constructor 
   * 
   */
  AliForwardCreateResponseMatrices();
  /**
   *  Constructor
   *
   */
  AliForwardCreateResponseMatrices(const char* name);
  /**
   * Copy Constructor
   *
   */
  AliForwardCreateResponseMatrices(const AliForwardCreateResponseMatrices& o) : AliBasedNdetaTask(o), fTrigger(),fBins(), fOutput() { }
  /**
   * Assignment operator
   *
   */
  AliForwardCreateResponseMatrices& operator=(const AliForwardCreateResponseMatrices&){return *this;}
  /** 
   * 
   * Destructor
   * 
   */
  virtual ~AliForwardCreateResponseMatrices(){}
  /**
   *  Embedded Class begins here
   */
  struct Bin : public TNamed
  {
    /**
     * Default Constructor
     */
    Bin();
    /**
     * Constructor
     */
    Bin(Double_t etaLow, Double_t etaHigh);
    /**
     * Copy Constructor
     */ 
    Bin(const Bin&);
    /**
     * Assignment Operator
     */
    Bin&operator=(const Bin&){return*this;}
    /**
     * Destructor
     */
    ~Bin(){}
    /**
     * Define outputs of a single eta bin
     */
    virtual void DefineOutputs(TList* cont,  Int_t max);
    /**
     * Process a single eta bin
     */    
    virtual void Process(TH1D* dndetaForward, TH1D* dndetaCentral, TH1D* normForward,   TH1D* normCentral, TH1D* dndetaMC, Double_t VtxZ, Bool_t selectedTrigger,  Bool_t isMCNSDm, Bool_t isESDNSD, AliAODEvent* aodevent);
    Double_t fEtaLow;                  // low eta limit 
    Double_t fEtaHigh;                 // high eta limit 
    TH1D*    fHist;                    // multiplicity histogram 
    TH1D*    fHistMC;                  // multiplicity histogram MC truth primaries
    TH2D*    fAcceptance;              // histogram showing the 'holes' in acceptance. BinContent of 1 shows a hole, and BinContent of 10 shows data coverage
    TH2D*    fVtxZvsNdataBins;         // VtxZ vs. number of data acceptance bins (normalised to the eta range) 
    TH2D*    fResponseMatrix;          //Response matrix (MC truth vs. analysed multiplicity)
    TH2D*    fResponseMatrixPlus05;    //Response matrix with analysed multiplicity scaled up by 5%
    TH2D*    fResponseMatrixPlus075;   //Response matrix  with analysed multiplicity scaled up by 7.5%
    TH2D*    fResponseMatrixPlus10;    //Response matrix with analysed multiplicity scaled up by 10%
    TH2D*    fResponseMatrixMinus05;   //Response matrix with analysed multiplicity scaled down by 5%
    TH2D*    fResponseMatrixMinus075;  //Response matrix with analysed multiplicity scaled down by 7.55%
    TH2D*    fResponseMatrixMinus10;   //Response matrix with analysed multiplicity scaled down by 10%
    TH2D*    fResponseMatrixMinusSys;  //Response matrix with analysed multiplicity scaled up by event mult uncertainty
    TH2D*    fResponseMatrixPlusSys;   //Response matrix with analysed multiplicity scaled down by event mult uncertainty
    TH1D*    fESDNSD;                  //number of events found as NSD by the analysis vs. multiplicity
    TH1D*    fMCNSD;                   //number of events found as NSD by the MC truth vs. multiplicity
    TH1D*    fMCESDNSD;                //number of events found as NSD by both analysis and MC truth vs. multiplicity
    TH1D*    fTriggerBias;             // histogram for trigger vertex bias correction
   ClassDef(Bin,1); // Manager of data 
  };
  /**
   * Create Output Objects
   */
  virtual void UserCreateOutputObjects();
  /**
   * User Exec
   */
  void UserExec(Option_t *option);
  /**
   * Terminate
   */
  void Terminate(Option_t *option);
  /** 
   * implementation of pure virtual function, always returning 0
   */
  virtual TH2D* GetHistogram(const AliAODEvent* aod, Bool_t mc);
  /** 
   * Get single event forward and central @f$d^2N/d\eta d\phi@f$
   * histograms
   * 
   */  
  virtual void GetHistograms(const AliAODEvent* aod, TH2D& forward, TH2D& central); 
  /** 
   * Add another eta bin to the task
   */
  void AddBin(Double_t etaLow, Double_t etaHigh){fBins.Add(new Bin(etaLow, etaHigh)); }
 protected:
  TH1I* fTrigger;  //Trigger histogram
  TList  fBins;    // List of eta bins
  TList*  fOutput; // Output list
  ClassDef(AliForwardCreateResponseMatrices, 1); 
};

#endif
// Local Variables:
//   mode: C++
// End:
Commit	Line	Data
ad7be237	1	#ifndef ALIFORWARDCREATERESPONSEMATRICES_H
	2	#define ALIFORWARDCREATERESPONSEMATRICES_H
	3
	4	/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
	5	* See cxx source for full Copyright notice */
	6
	7	#include "AliAnalysisTaskSE.h"
	8	#include "AliBasedNdetaTask.h"
	9	#include <TList.h>
	10	#include <iostream>
	11
	12	class TH2D;
	13
	14	/**
	15	* Task to do the multiplicity distibution
	16	*
	17	*/
	18	class AliForwardCreateResponseMatrices : public AliBasedNdetaTask
	19	{
	20	public:
	21	/**
	22	*
	23	* Default Constructor
	24	*
	25	*/
	26	AliForwardCreateResponseMatrices();
	27	/**
	28	* Constructor
	29	*
	30	*/
	31	AliForwardCreateResponseMatrices(const char* name);
	32	/**
	33	* Copy Constructor
	34	*
	35	*/
	36	AliForwardCreateResponseMatrices(const AliForwardCreateResponseMatrices& o) : AliBasedNdetaTask(o), fTrigger(),fBins(), fOutput() { }
	37	/**
	38	* Assignment operator
	39	*
	40	*/
	41	AliForwardCreateResponseMatrices& operator=(const AliForwardCreateResponseMatrices&){return *this;}
	42	/**
	43	*
	44	* Destructor
	45	*
	46	*/
	47	virtual ~AliForwardCreateResponseMatrices(){}
	48	/**
	49	* Embedded Class begins here
	50	*/
	51	struct Bin : public TNamed
	52	{
	53	/**
	54	* Default Constructor
	55	*/
	56	Bin();
	57	/**
	58	* Constructor
	59	*/
	60	Bin(Double_t etaLow, Double_t etaHigh);
	61	/**
	62	* Copy Constructor
	63	*/
	64	Bin(const Bin&);
65	/**
66	* Assignment Operator
67	*/
68	Bin&operator=(const Bin&){return*this;}
69	/**
70	* Destructor
71	*/
72	~Bin(){}
73	/**
74	* Define outputs of a single eta bin
75	*/
76	virtual void DefineOutputs(TList* cont, Int_t max);
77	/**
78	* Process a single eta bin
79	*/
80	virtual void Process(TH1D* dndetaForward, TH1D* dndetaCentral, TH1D* normForward, TH1D* normCentral, TH1D* dndetaMC, Double_t VtxZ, Bool_t selectedTrigger, Bool_t isMCNSDm, Bool_t isESDNSD, AliAODEvent* aodevent);
81	Double_t fEtaLow; // low eta limit
82	Double_t fEtaHigh; // high eta limit
83	TH1D* fHist; // multiplicity histogram
84	TH1D* fHistMC; // multiplicity histogram MC truth primaries
85	TH2D* fAcceptance; // histogram showing the 'holes' in acceptance. BinContent of 1 shows a hole, and BinContent of 10 shows data coverage
86	TH2D* fVtxZvsNdataBins; // VtxZ vs. number of data acceptance bins (normalised to the eta range)
87	TH2D* fResponseMatrix; //Response matrix (MC truth vs. analysed multiplicity)
88	TH2D* fResponseMatrixPlus05; //Response matrix with analysed multiplicity scaled up by 5%
89	TH2D* fResponseMatrixPlus075; //Response matrix with analysed multiplicity scaled up by 7.5%
90	TH2D* fResponseMatrixPlus10; //Response matrix with analysed multiplicity scaled up by 10%
91	TH2D* fResponseMatrixMinus05; //Response matrix with analysed multiplicity scaled down by 5%
92	TH2D* fResponseMatrixMinus075; //Response matrix with analysed multiplicity scaled down by 7.55%
93	TH2D* fResponseMatrixMinus10; //Response matrix with analysed multiplicity scaled down by 10%
94	TH2D* fResponseMatrixMinusSys; //Response matrix with analysed multiplicity scaled up by event mult uncertainty
95	TH2D* fResponseMatrixPlusSys; //Response matrix with analysed multiplicity scaled down by event mult uncertainty
96	TH1D* fESDNSD; //number of events found as NSD by the analysis vs. multiplicity
97	TH1D* fMCNSD; //number of events found as NSD by the MC truth vs. multiplicity
98	TH1D* fMCESDNSD; //number of events found as NSD by both analysis and MC truth vs. multiplicity
99	TH1D* fTriggerBias; // histogram for trigger vertex bias correction
100	ClassDef(Bin,1); // Manager of data
101	};
102	/**
103	* Create Output Objects
104	*/
105	virtual void UserCreateOutputObjects();
106	/**
107	* User Exec
108	*/
109	void UserExec(Option_t *option);
110	/**
111	* Terminate
112	*/
113	void Terminate(Option_t *option);
114	/**
115	* implementation of pure virtual function, always returning 0
116	*/
117	virtual TH2D* GetHistogram(const AliAODEvent* aod, Bool_t mc);
118	/**
119	* Get single event forward and central @f$d^2N/d\eta d\phi@f$
120	* histograms
121	*
122	*/
123	virtual void GetHistograms(const AliAODEvent* aod, TH2D& forward, TH2D& central);
124	/**
125	* Add another eta bin to the task
126	*/
127	void AddBin(Double_t etaLow, Double_t etaHigh){fBins.Add(new Bin(etaLow, etaHigh)); }
128	protected:
129	TH1I* fTrigger; //Trigger histogram
130	TList fBins; // List of eta bins
131	TList* fOutput; // Output list
132	ClassDef(AliForwardCreateResponseMatrices, 1);
133	};
134
135	#endif
136	// Local Variables:
137	// mode: C++
138	// End: