[u/mrichter/AliRoot.git] / PWGLF / totEt / AliAnalysisEtCuts.cxx

//_________________________________________________________________________
//  Utility Class for transverse energy studies
//  Base class for ESD & MC analysis
//  - cuts for reconstruction and MonteCarlo 
// implementation file
//  
//*-- Authors: Oystein Djuvsland (Bergen), David Silvermyr (ORNL)
//_________________________________________________________________________

#include "AliAnalysisEtCuts.h"

#include <iostream>
#include "TNamed.h"

using namespace std;

ClassImp(AliAnalysisEtCuts);


AliAnalysisEtCuts::AliAnalysisEtCuts() : 
  TNamed()
				    //
  ,fCommonEtaCut(0.5)
  ,fCommonClusterEnergyCut(0.15)
  ,fCommonTrackPtCut(0.0)
  ,fCommonSingleCell(1)
  ,fEmcalTrackDistanceCut(15.0)
  ,fEmcalTrackDxCut(5.0)
  ,fEmcalTrackDzCut(5.0)
  ,fPhosTrackDistanceCut(10.0)  
  ,fPhosTrackDxCut(8.0)
  ,fPhosTrackDzCut(3.0)
  ,fPhosTrackRCut(5.0)
  ,fPhosBadDistanceCut(3.0)
  
  ,fGeometryPhosEtaAccCut(0.1)//gives wiggle room
  ,fGeometryPhosPhiAccMinCut(-100)//260.0)
  ,fGeometryPhosPhiAccMaxCut(-40)//320.0)
  ,fGeometryPhosDetectorRadius(460.0)
				    //
  ,fGeometryEmcalEtaAccCut(0.6)
  ,fGeometryEmcalPhiAccMinCut(80.0) // rad 1.4
  ,fGeometryEmcalPhiAccMaxCut(120.0) // rad 2.1
				    //,fGeometryEmcalPhiAccMaxCut(180.0) // rad 3.14
  ,fGeometryEmcalDetectorRadius(440.0)
				    //
  ,fReconstructedVertexXCut(0.5)
  ,fReconstructedVertexYCut(0.5)
  ,fReconstructedVertexZCut(12.0)
  ,fReconstructedIPxyCut(1.5)
  ,fReconstructedIPzCut(1.5)
  ,fReconstructedNTpcClustersCut(30)
  ,fReconstructedNItsClustersCut(3)
  ,fReconstructedPidCut(0.0)
				    //
  ,fReconstructedPhosClusterType(-1)
  ,fReconstructedPhosClusterEnergyCut(0.25)// GeV
  ,fReconstructedPhosSingleCellEnergyCut(0.5)
  ,fReconstructedPhosTrackDistanceTightCut(3.0)
  ,fReconstructedPhosTrackDistanceMediumCut(5.0)
  ,fReconstructedPhosTrackDistanceLooseCut(15.0)
				    //
  ,fReconstructedEmcalClusterType(1)
  ,fReconstructedEmcalClusterEnergyCut(0.30) // GeV
  ,fReconstructedEmcalSingleCellEnergyCut(0.5)
  ,fReconstructedEmcalTrackDistanceTightCut(5.0)
  ,fReconstructedEmcalTrackDistanceMediumCut(10.0)
  ,fReconstructedEmcalTrackDistanceLooseCut(15.0)
  
  ,fMonteCarloSingleChargedParticle(3)
  ,fMonteCarloNeutralParticle(0)

  ,fHistMakeTree(kFALSE)
  ,fHistMakeTreeDeposit(kFALSE)
  ,fHistNbinsMult(2000)
  ,fHistMinMult(-0.5)
  ,fHistMaxMult(1999.5)
  ,fHistNbinsTotEt(10000)
  ,fHistMinTotEt(0.000)
  ,fHistMaxTotEt(1000)
  ,fHistNbinsParticleEt(5000)
  ,fHistMinParticleEt(0)
  ,fHistMaxParticleEt(500)
  ,fHistNbinsParticlePt(200) 
  ,fHistMinParticlePt(0)
  ,fHistMaxParticlePt(20)
  
  ,fPrimaryVertexCutXY(4.0)
  ,fPrimaryVertexCutZ(20.0)
{ // ctor
}

AliAnalysisEtCuts::~AliAnalysisEtCuts()
{ // dtor
}

void AliAnalysisEtCuts::SetPbPbDefaults()
{ // just increase seom histogram max values for now
  // enough to multiply conservative p+p defaults by a factor 100?
  fHistMaxMult = 20000;
  fHistMaxTotEt = 10000;
}
Commit	Line	Data
4998becf	1	//_________________________________________________________________________
	2	// Utility Class for transverse energy studies
	3	// Base class for ESD & MC analysis
	4	// - cuts for reconstruction and MonteCarlo
	5	// implementation file
	6	//
	7	//*-- Authors: Oystein Djuvsland (Bergen), David Silvermyr (ORNL)
	8	//_________________________________________________________________________
	9
	10	#include "AliAnalysisEtCuts.h"
	11
	12	#include <iostream>
	13	#include "TNamed.h"
	14
	15	using namespace std;
	16
	17	ClassImp(AliAnalysisEtCuts);
	18
	19
	20	AliAnalysisEtCuts::AliAnalysisEtCuts() :
	21	TNamed()
	22	//
393b8149	23	,fCommonEtaCut(0.5)
ef647350	24	,fCommonClusterEnergyCut(0.15)
4998becf	25	,fCommonTrackPtCut(0.0)
4998becf	26	,fCommonSingleCell(1)
393b8149	27	,fEmcalTrackDistanceCut(15.0)
	28	,fEmcalTrackDxCut(5.0)
	29	,fEmcalTrackDzCut(5.0)
	30	,fPhosTrackDistanceCut(10.0)
	31	,fPhosTrackDxCut(8.0)
	32	,fPhosTrackDzCut(3.0)
b2c10007	33	,fPhosTrackRCut(5.0)
ef647350	34	,fPhosBadDistanceCut(3.0)
393b8149	35
cc8322de	36	,fGeometryPhosEtaAccCut(0.1)//gives wiggle room
43dd5a38	37	,fGeometryPhosPhiAccMinCut(-100)//260.0)
43dd5a38	38	,fGeometryPhosPhiAccMaxCut(-40)//320.0)
4998becf	39	,fGeometryPhosDetectorRadius(460.0)
4998becf	40	//
76edf1b5	41	,fGeometryEmcalEtaAccCut(0.6)
4998becf	42	,fGeometryEmcalPhiAccMinCut(80.0) // rad 1.4
6ab109c8	43	,fGeometryEmcalPhiAccMaxCut(120.0) // rad 2.1
6ab109c8	44	//,fGeometryEmcalPhiAccMaxCut(180.0) // rad 3.14
4998becf	45	,fGeometryEmcalDetectorRadius(440.0)
	46	//
	47	,fReconstructedVertexXCut(0.5)
	48	,fReconstructedVertexYCut(0.5)
	49	,fReconstructedVertexZCut(12.0)
	50	,fReconstructedIPxyCut(1.5)
	51	,fReconstructedIPzCut(1.5)
	52	,fReconstructedNTpcClustersCut(30)
	53	,fReconstructedNItsClustersCut(3)
8e6681b9	54	,fReconstructedPidCut(0.0)
4998becf	55	//
4998becf	56	,fReconstructedPhosClusterType(-1)
b813e06f	57	,fReconstructedPhosClusterEnergyCut(0.25)// GeV
4998becf	58	,fReconstructedPhosSingleCellEnergyCut(0.5)
393b8149	59	,fReconstructedPhosTrackDistanceTightCut(3.0)
	60	,fReconstructedPhosTrackDistanceMediumCut(5.0)
	61	,fReconstructedPhosTrackDistanceLooseCut(15.0)
4998becf	62	//
e16aac77	63	,fReconstructedEmcalClusterType(1)
b813e06f	64	,fReconstructedEmcalClusterEnergyCut(0.30) // GeV
4998becf	65	,fReconstructedEmcalSingleCellEnergyCut(0.5)
393b8149	66	,fReconstructedEmcalTrackDistanceTightCut(5.0)
	67	,fReconstructedEmcalTrackDistanceMediumCut(10.0)
	68	,fReconstructedEmcalTrackDistanceLooseCut(15.0)
4998becf	69
	70	,fMonteCarloSingleChargedParticle(3)
	71	,fMonteCarloNeutralParticle(0)
ce546038	72
3f510c32	73	,fHistMakeTree(kFALSE)
393b8149	74	,fHistMakeTreeDeposit(kFALSE)
393b8149	75	,fHistNbinsMult(2000)
0fa8c632	76	,fHistMinMult(-0.5)
393b8149	77	,fHistMaxMult(1999.5)
	78	,fHistNbinsTotEt(10000)
	79	,fHistMinTotEt(0.000)
	80	,fHistMaxTotEt(1000)
	81	,fHistNbinsParticleEt(5000)
0fa8c632	82	,fHistMinParticleEt(0)
393b8149	83	,fHistMaxParticleEt(500)
0fa8c632	84	,fHistNbinsParticlePt(200)
	85	,fHistMinParticlePt(0)
	86	,fHistMaxParticlePt(20)
b2c10007	87
	88	,fPrimaryVertexCutXY(4.0)
	89	,fPrimaryVertexCutZ(20.0)
4998becf	90	{ // ctor
	91	}
	92
	93	AliAnalysisEtCuts::~AliAnalysisEtCuts()
	94	{ // dtor
	95	}
	96
0fa8c632	97	void AliAnalysisEtCuts::SetPbPbDefaults()
	98	{ // just increase seom histogram max values for now
	99	// enough to multiply conservative p+p defaults by a factor 100?
	100	fHistMaxMult = 20000;
	101	fHistMaxTotEt = 10000;
	102	}
	103
4998becf	104