[u/mrichter/AliRoot.git] / PWG4 / AliGammaReader.cxx

/**************************************************************************
 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
 *                                                                        *
 * Author: The ALICE Off-line Project.                                    *
 * Contributors are mentioned in the code where appropriate.              *
 *                                                                        *
 * Permission to use, copy, modify and distribute this software and its   *
 * documentation strictly for non-commercial purposes is hereby granted   *
 * without fee, provided that the above copyright notice appears in all   *
 * copies and that both the copyright notice and this permission notice   *
 * appear in the supporting documentation. The authors make no claims     *
 * about the suitability of this software for any purpose. It is          *
 * provided "as is" without express or implied warranty.                  *
 **************************************************************************/
/* $Id$ */

/* History of cvs commits:
 *
 * $Log$
 * Revision 1.1.2.1  2007/07/26 10:32:09  schutz
 * new analysis classes in the the new analysis framework
 *
 *
 */

//_________________________________________________________________________
// Base class for reading data in order to do prompt gamma correlations
//*-- Author: Gustavo Conesa (LNF-INFN) 
//////////////////////////////////////////////////////////////////////////////


// --- ROOT system ---

//---- ANALYSIS system ----
#include "Riostream.h"
#include "AliLog.h"
#include "AliGammaReader.h"

ClassImp(AliGammaReader)


//____________________________________________________________________________
AliGammaReader::AliGammaReader() : 
  TObject(), fDataType(0),
  fCTSEtaCut(0.), fEMCALEtaCut(0.), fPHOSEtaCut(0.),
  fNeutralPtCut(0.),
  fChargedPtCut(0.)
{
  //Ctor

  fPhiEMCALCut[0]=0.;
  fPhiEMCALCut[1]=0.;
  fPhiPHOSCut[0]=0.;
  fPhiPHOSCut[1]=0.;
 
  //Initialize parameters
  InitParameters();
}

//____________________________________________________________________________
AliGammaReader::AliGammaReader(const AliGammaReader & g) :   
  TObject(g), fDataType(g.fDataType),
  fCTSEtaCut(g.fCTSEtaCut),  fEMCALEtaCut(g.fEMCALEtaCut),  fPHOSEtaCut(g.fPHOSEtaCut),
  fNeutralPtCut(g.fNeutralPtCut),
  fChargedPtCut(g.fChargedPtCut)
{
  // cpy ctor

  fPhiEMCALCut[0]=g.fPhiEMCALCut[0];
  fPhiEMCALCut[1]=g.fPhiEMCALCut[1];
  fPhiPHOSCut[0]=g.fPhiPHOSCut[0];
  fPhiPHOSCut[1]=g.fPhiPHOSCut[1];
}

//_________________________________________________________________________
AliGammaReader & AliGammaReader::operator = (const AliGammaReader & source)
{
  // assignment operator

  if(&source == this) return *this;

  fDataType = source.fDataType ;
  fCTSEtaCut = source.fCTSEtaCut;  
  fEMCALEtaCut = source.fEMCALEtaCut;  
  fPHOSEtaCut = source.fPHOSEtaCut;
  fNeutralPtCut = source.fNeutralPtCut;
  fChargedPtCut = source.fChargedPtCut; 

  fPhiEMCALCut[0]=source.fPhiEMCALCut[0];
  fPhiEMCALCut[1]=source.fPhiEMCALCut[1];
  fPhiPHOSCut[0]=source.fPhiPHOSCut[0];
  fPhiPHOSCut[1]=source.fPhiPHOSCut[1];

  return *this;

}

//_______________________________________________________________
void AliGammaReader::InitParameters()
{
 
  //Initialize the parameters of the analysis.
  fDataType = kData ;
  fCTSEtaCut         = 0.7 ;  
  fEMCALEtaCut         = 0.7 ;  
  fPHOSEtaCut         = 0.12 ;
  fPhiEMCALCut[0] = 80 *TMath::DegToRad();
  fPhiEMCALCut[1] = 190*TMath::DegToRad();
  fPhiPHOSCut[0] = 220. *TMath::DegToRad();
  fPhiPHOSCut[1] = 320.*TMath::DegToRad();
  fNeutralPtCut   = 0.5 ;
  fChargedPtCut   = 0.5 ;

}


//________________________________________________________________
void AliGammaReader::Print(const Option_t * opt) const
{

  //Print some relevant parameters set for the analysis
  if(! opt)
    return;

  Info("Print", "%s %s", GetName(), GetTitle() ) ;
  printf("Data type           : %d\n", fDataType) ;
  printf("CTS Eta cut           : %f\n", fCTSEtaCut) ;
  printf("EMCAL Eta cut           : %f\n", fEMCALEtaCut) ;
  printf("PHOS Eta cut           : %f\n", fPHOSEtaCut) ;
  printf("Phi EMCAL cut           : [%f, %f]\n", fPhiEMCALCut[0],fPhiEMCALCut[1]) ;
  printf("Phi PHOS cut           : [%f, %f]\n", fPhiPHOSCut[0],fPhiPHOSCut[1]) ;
  printf("pT neutral cut           : %f GeV/c\n", fNeutralPtCut) ;
  printf("pT charged cut           : %f GeV/c\n", fChargedPtCut) ;

}
Commit	Line	Data
bdcfac30	1	/**************************************************************************
	2	* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
	3	* *
	4	* Author: The ALICE Off-line Project. *
	5	* Contributors are mentioned in the code where appropriate. *
	6	* *
	7	* Permission to use, copy, modify and distribute this software and its *
	8	* documentation strictly for non-commercial purposes is hereby granted *
	9	* without fee, provided that the above copyright notice appears in all *
	10	* copies and that both the copyright notice and this permission notice *
	11	* appear in the supporting documentation. The authors make no claims *
	12	* about the suitability of this software for any purpose. It is *
	13	* provided "as is" without express or implied warranty. *
	14	**************************************************************************/
	15	/* $Id$ */
	16
	17	/* History of cvs commits:
	18	*
	19	* $Log$
	20	* Revision 1.1.2.1 2007/07/26 10:32:09 schutz
	21	* new analysis classes in the the new analysis framework
	22	*
	23	*
	24	*/
	25
	26	//_________________________________________________________________________
	27	// Base class for reading data in order to do prompt gamma correlations
	28	//*-- Author: Gustavo Conesa (LNF-INFN)
	29	//////////////////////////////////////////////////////////////////////////////
	30
	31
	32	// --- ROOT system ---
	33
	34	//---- ANALYSIS system ----
	35	#include "Riostream.h"
	36	#include "AliLog.h"
	37	#include "AliGammaReader.h"
	38
	39	ClassImp(AliGammaReader)
	40
	41
	42	//____________________________________________________________________________
	43	AliGammaReader::AliGammaReader() :
	44	TObject(), fDataType(0),
	45	fCTSEtaCut(0.), fEMCALEtaCut(0.), fPHOSEtaCut(0.),
	46	fNeutralPtCut(0.),
	47	fChargedPtCut(0.)
	48	{
	49	//Ctor
	50
	51	fPhiEMCALCut[0]=0.;
	52	fPhiEMCALCut[1]=0.;
	53	fPhiPHOSCut[0]=0.;
	54	fPhiPHOSCut[1]=0.;
	55
	56	//Initialize parameters
	57	InitParameters();
	58	}
	59
	60	//____________________________________________________________________________
	61	AliGammaReader::AliGammaReader(const AliGammaReader & g) :
	62	TObject(g), fDataType(g.fDataType),
	63	fCTSEtaCut(g.fCTSEtaCut), fEMCALEtaCut(g.fEMCALEtaCut), fPHOSEtaCut(g.fPHOSEtaCut),
	64	fNeutralPtCut(g.fNeutralPtCut),
65	fChargedPtCut(g.fChargedPtCut)
66	{
67	// cpy ctor
68
69	fPhiEMCALCut[0]=g.fPhiEMCALCut[0];
70	fPhiEMCALCut[1]=g.fPhiEMCALCut[1];
71	fPhiPHOSCut[0]=g.fPhiPHOSCut[0];
72	fPhiPHOSCut[1]=g.fPhiPHOSCut[1];
73	}
74
75	//_________________________________________________________________________
76	AliGammaReader & AliGammaReader::operator = (const AliGammaReader & source)
77	{
78	// assignment operator
79
80	if(&source == this) return *this;
81
82	fDataType = source.fDataType ;
83	fCTSEtaCut = source.fCTSEtaCut;
84	fEMCALEtaCut = source.fEMCALEtaCut;
85	fPHOSEtaCut = source.fPHOSEtaCut;
86	fNeutralPtCut = source.fNeutralPtCut;
87	fChargedPtCut = source.fChargedPtCut;
88
89	fPhiEMCALCut[0]=source.fPhiEMCALCut[0];
90	fPhiEMCALCut[1]=source.fPhiEMCALCut[1];
91	fPhiPHOSCut[0]=source.fPhiPHOSCut[0];
92	fPhiPHOSCut[1]=source.fPhiPHOSCut[1];
93
94	return *this;
95
96	}
97
98	//_______________________________________________________________
99	void AliGammaReader::InitParameters()
100	{
101
102	//Initialize the parameters of the analysis.
103	fDataType = kData ;
104	fCTSEtaCut = 0.7 ;
105	fEMCALEtaCut = 0.7 ;
106	fPHOSEtaCut = 0.12 ;
107	fPhiEMCALCut[0] = 80 *TMath::DegToRad();
108	fPhiEMCALCut[1] = 190*TMath::DegToRad();
109	fPhiPHOSCut[0] = 220. *TMath::DegToRad();
110	fPhiPHOSCut[1] = 320.*TMath::DegToRad();
111	fNeutralPtCut = 0.5 ;
112	fChargedPtCut = 0.5 ;
113
114	}
115
116
117	//________________________________________________________________
118	void AliGammaReader::Print(const Option_t * opt) const
119	{
120
121	//Print some relevant parameters set for the analysis
122	if(! opt)
123	return;
124
125	Info("Print", "%s %s", GetName(), GetTitle() ) ;
126	printf("Data type : %d\n", fDataType) ;
127	printf("CTS Eta cut : %f\n", fCTSEtaCut) ;
128	printf("EMCAL Eta cut : %f\n", fEMCALEtaCut) ;
129	printf("PHOS Eta cut : %f\n", fPHOSEtaCut) ;
130	printf("Phi EMCAL cut : [%f, %f]\n", fPhiEMCALCut[0],fPhiEMCALCut[1]) ;
131	printf("Phi PHOS cut : [%f, %f]\n", fPhiPHOSCut[0],fPhiPHOSCut[1]) ;
132	printf("pT neutral cut : %f GeV/c\n", fNeutralPtCut) ;
133	printf("pT charged cut : %f GeV/c\n", fChargedPtCut) ;
134
135	}
136
137