[u/mrichter/AliRoot.git] / PHOS / AliPHOSRecoParam.cxx

/**************************************************************************
 * Copyright(c) 2007, 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$ */

// Base class for the PHOS reconstruction parameters.
// Do not use in the reconstruction; use derivative classes instead.
// Author: Boris Polichtchouk.

// --- AliRoot header files ---
#include "AliCDBManager.h"
#include "AliCDBEntry.h"
#include "AliLog.h"
#include "AliPHOSRecoParam.h"

ClassImp(AliPHOSRecoParam)

TObjArray* AliPHOSRecoParam::fgkMaps =0; //ALTRO mappings
TObject*   AliPHOSRecoParam::fgkTrigParams =0; // Trigger parameters

//-----------------------------------------------------------------------------
AliPHOSRecoParam::AliPHOSRecoParam() :
  AliDetectorRecoParam(),
  fNonLinearityParams(3),
  fEMCClusteringThreshold(0.2),
  fEMCLocMaxCut(0.03),
  fEMCRawDigitThreshold(2),
  fEMCMinE(0.012),
  fEMCW0(4.5),
  fEMCSampleQualityCut(1.),
  fTimeGateAmpThresh(10.),
  fTimeGateLow(1.e-6),
  fTimeGateHigh(1.e-6),
  fEMCEcoreRadius(3.),
  fEMCEcore2ESD(kFALSE),
  fEMCSubtractPedestals(kTRUE),
  fEMCUnfold(kTRUE),
  fEMCEnergyCorrectionOn(kTRUE),
  fEMCFitterVersion(""),
  fNonlinearityCorrVersion(""),
  fGlobalAltroOffset(0),
  fGlobalAltroThreshold(0),
  fCPVClusteringThreshold(0.0),
  fCPVLocMaxCut(0.03),
  fCPVMinE(0.0),
  fCPVW0(4.0),
  fCPVUnfold(kTRUE)
{
  //Default constructor.
  SetNonlinearityCorrectionVersion("Gustavo2005") ;
}

//-----------------------------------------------------------------------------
AliPHOSRecoParam::AliPHOSRecoParam(const AliPHOSRecoParam& ):
  AliDetectorRecoParam(),
  fNonLinearityParams(3),
  fEMCClusteringThreshold(0.2),
  fEMCLocMaxCut(0.03),
  fEMCRawDigitThreshold(2),
  fEMCMinE(0.012),
  fEMCW0(4.5),
  fEMCSampleQualityCut(1.),
  fTimeGateAmpThresh(10.),
  fTimeGateLow(1.e-6),
  fTimeGateHigh(1.e-6),
  fEMCEcoreRadius(3.),
  fEMCEcore2ESD(kFALSE),
  fEMCSubtractPedestals(kTRUE),
  fEMCUnfold(kTRUE),
  fEMCEnergyCorrectionOn(kTRUE),
  fEMCFitterVersion(""),
  fNonlinearityCorrVersion(""),
  fGlobalAltroOffset(0),
  fGlobalAltroThreshold(0),
  fCPVClusteringThreshold(0.0),
  fCPVLocMaxCut(0.03),
  fCPVMinE(0.0),
  fCPVW0(4.0),
  fCPVUnfold(kTRUE)
{
  //Copy constructor.
}

//-----------------------------------------------------------------------------
AliPHOSRecoParam& AliPHOSRecoParam::operator = (const AliPHOSRecoParam& recoParam)
{
  //Assignment operator.

  if(this != &recoParam) {
    fNonLinearityParams     = recoParam.fNonLinearityParams;
    fEMCClusteringThreshold = recoParam.fEMCClusteringThreshold;
    fEMCLocMaxCut           = recoParam.fEMCLocMaxCut;
    fEMCRawDigitThreshold   = recoParam.fEMCRawDigitThreshold;
    fEMCMinE                = recoParam.fEMCMinE;
    fEMCW0                  = recoParam.fEMCW0;
    fEMCSampleQualityCut    = recoParam.fEMCSampleQualityCut;
    fTimeGateAmpThresh      = recoParam.fTimeGateAmpThresh ;
    fTimeGateLow            = recoParam.fTimeGateLow ;
    fTimeGateHigh           = recoParam.fTimeGateHigh ;
    fEMCEcoreRadius         = recoParam.fEMCEcoreRadius;
    fEMCEcore2ESD           = recoParam.fEMCEcore2ESD;
    fEMCSubtractPedestals   = recoParam.fEMCSubtractPedestals;
    fEMCUnfold              = recoParam.fEMCUnfold;
    fEMCEnergyCorrectionOn  = recoParam.fEMCEnergyCorrectionOn;
    fEMCFitterVersion       = recoParam.fEMCFitterVersion;
    fNonlinearityCorrVersion= recoParam.fNonlinearityCorrVersion;
    fGlobalAltroOffset      = recoParam.fGlobalAltroOffset;
    fGlobalAltroThreshold   = recoParam.fGlobalAltroThreshold;
    fCPVClusteringThreshold = recoParam.fCPVClusteringThreshold;
    fCPVLocMaxCut           = recoParam.fCPVLocMaxCut;
    fCPVMinE                = recoParam.fCPVMinE;
    fCPVW0                  = recoParam.fCPVW0;
    fCPVUnfold              = recoParam.fCPVUnfold;
    fNonLinearityParams     = recoParam.fNonLinearityParams ;
  }

  return *this;
}

//-----------------------------------------------------------------------------
void AliPHOSRecoParam::Print(Option_t * /*option*/) const
{
  AliDebug(2,Form("PHOS reconstruction parameters:\n"
		  "\tEMCClusteringThreshold = %f\n"
		  "\tEMCLocMaxCut           = %f\n"
		  "\tEMCRawDigitThreshold   = %f\n"
		  "\tEMCMinE                = %f\n"
		  "\tEMCW0                  = %f\n"
		  "\tEMCSampleQualityCut    = %f\n"
		  "\tEMCEcoreRadius         = %f\n"
		  "\tEMCEcore2ESD           = %d\n"
		  "\tEMCSubtractPedestals   = %d\n"
		  "\tEMCUnfold              = %d\n"
		  "\tEMCEnergyCorrectionOn  = %d\n"
		  "\tEMCFitterVersion       = \"%s\"\n"
		  "\tGlobalAltroOffset      = %d\n"
		  "\tGlobalAltroThreshold   = %d",
		  fEMCClusteringThreshold,
		  fEMCLocMaxCut,
		  fEMCRawDigitThreshold,
		  fEMCMinE,
		  fEMCW0,
		  fEMCSampleQualityCut,
		  fEMCEcoreRadius,
		  fEMCEcore2ESD,
		  fEMCSubtractPedestals,
		  fEMCUnfold,
		  fEMCEnergyCorrectionOn,
		  fEMCFitterVersion.Data(),
		  fGlobalAltroOffset,
		  fGlobalAltroThreshold));

}

//-----------------------------------------------------------------------------
AliPHOSRecoParam* AliPHOSRecoParam::GetDefaultParameters()
{
  //Default parameters for the reconstruction

  AliPHOSRecoParam* params = new AliPHOSRecoParam();
  return params;
}

//-----------------------------------------------------------------------------
const TObjArray* AliPHOSRecoParam::GetMappings()
{
  //Returns array of AliAltroMappings for RCU0..RCU3.
  //If not found, read it from OCDB.

  //Quick check as follows:
  //  root [0] AliCDBManager::Instance()->SetDefaultStorage("local://$ALICE_ROOT/OCDB");
  //  root [1] AliCDBManager::Instance()->SetRun(1);
  //  root [2] TObjArray* maps = AliPHOSRecoParam::GetMappings();
  //  root [3] maps->Print();
  
  if(fgkMaps) return fgkMaps;
  
  AliCDBEntry* entry = AliCDBManager::Instance()->Get("PHOS/Calib/Mapping");
  if(entry)
    fgkMaps = (TObjArray*)entry->GetObject();
  
  return fgkMaps;
  
}
//-----------------------------------------------------------------------------
void AliPHOSRecoParam::SetNonlinearityCorrectionVersion(const char * ver){
  //Set non-linearity correction version and adjust parameters

  if(strcmp(ver,"NoCorrection")==0){
    fNonLinearityParams.Set(1) ;
    fNonLinearityParams.AddAt(0.,0) ;
    fNonlinearityCorrVersion=ver ;
    return ;
  }
  if(strcmp(ver,"Henrik2010")==0){
    Float_t pars[7]={1.051,2.54e-003,-1.737e-002,5.69e-002,3.3e-001,1.42e-001,1.50e-002} ;
    fNonLinearityParams.Set(7,pars) ;
    fNonlinearityCorrVersion=ver ;
    return ;
  }
  if(strcmp(ver,"Gustavo2005")==0){
    Float_t pars[3]={0.0241, 1.0504, 0.000249} ;
    fNonLinearityParams.Set(3,pars) ;
    fNonlinearityCorrVersion=ver ;
    return ;
  }
  AliError(Form("Non known correction version: %s, still using default \n",ver)) ;


}

//-----------------------------------------------------------------------------
const TObject* AliPHOSRecoParam::GetTriggerParameters()
{
  //Returns trigger parameters.

  //Quick check as follows:
  //  root [0] AliCDBManager::Instance()->SetDefaultStorage("local://$ALICE_ROOT/OCDB");
  //  root [1] AliCDBManager::Instance()->SetRun(158171);
  //  root [2] TObject* parameters = AliPHOSRecoParam::GetTriggerParameters();
  //  root [3] parameters->Print();
  
  if(fgkTrigParams) return fgkTrigParams;
  
  AliCDBEntry* entry = AliCDBManager::Instance()->Get("PHOS/Trigger/Parameters");
  if(entry)
    fgkTrigParams = entry->GetObject();
  
  return fgkTrigParams;
  
}
Commit	Line	Data
	1	/**************************************************************************
	2	* Copyright(c) 2007, 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
	16	/* $Id$ */
	17
	18	// Base class for the PHOS reconstruction parameters.
	19	// Do not use in the reconstruction; use derivative classes instead.
	20	// Author: Boris Polichtchouk.
	21
	22	// --- AliRoot header files ---
	23	#include "AliCDBManager.h"
	24	#include "AliCDBEntry.h"
	25	#include "AliLog.h"
	26	#include "AliPHOSRecoParam.h"
	27
	28	ClassImp(AliPHOSRecoParam)
	29
	30	TObjArray* AliPHOSRecoParam::fgkMaps =0; //ALTRO mappings
	31	TObject* AliPHOSRecoParam::fgkTrigParams =0; // Trigger parameters
	32
	33	//-----------------------------------------------------------------------------
	34	AliPHOSRecoParam::AliPHOSRecoParam() :
	35	AliDetectorRecoParam(),
	36	fNonLinearityParams(3),
	37	fEMCClusteringThreshold(0.2),
	38	fEMCLocMaxCut(0.03),
	39	fEMCRawDigitThreshold(2),
	40	fEMCMinE(0.012),
	41	fEMCW0(4.5),
	42	fEMCSampleQualityCut(1.),
	43	fTimeGateAmpThresh(10.),
	44	fTimeGateLow(1.e-6),
	45	fTimeGateHigh(1.e-6),
	46	fEMCEcoreRadius(3.),
	47	fEMCEcore2ESD(kFALSE),
	48	fEMCSubtractPedestals(kTRUE),
	49	fEMCUnfold(kTRUE),
	50	fEMCEnergyCorrectionOn(kTRUE),
	51	fEMCFitterVersion(""),
	52	fNonlinearityCorrVersion(""),
	53	fGlobalAltroOffset(0),
	54	fGlobalAltroThreshold(0),
	55	fCPVClusteringThreshold(0.0),
	56	fCPVLocMaxCut(0.03),
	57	fCPVMinE(0.0),
	58	fCPVW0(4.0),
	59	fCPVUnfold(kTRUE)
	60	{
	61	//Default constructor.
	62	SetNonlinearityCorrectionVersion("Gustavo2005") ;
	63	}
	64
	65	//-----------------------------------------------------------------------------
	66	AliPHOSRecoParam::AliPHOSRecoParam(const AliPHOSRecoParam& ):
	67	AliDetectorRecoParam(),
	68	fNonLinearityParams(3),
	69	fEMCClusteringThreshold(0.2),
	70	fEMCLocMaxCut(0.03),
	71	fEMCRawDigitThreshold(2),
	72	fEMCMinE(0.012),
	73	fEMCW0(4.5),
	74	fEMCSampleQualityCut(1.),
	75	fTimeGateAmpThresh(10.),
	76	fTimeGateLow(1.e-6),
	77	fTimeGateHigh(1.e-6),
	78	fEMCEcoreRadius(3.),
	79	fEMCEcore2ESD(kFALSE),
	80	fEMCSubtractPedestals(kTRUE),
	81	fEMCUnfold(kTRUE),
	82	fEMCEnergyCorrectionOn(kTRUE),
	83	fEMCFitterVersion(""),
	84	fNonlinearityCorrVersion(""),
	85	fGlobalAltroOffset(0),
	86	fGlobalAltroThreshold(0),
	87	fCPVClusteringThreshold(0.0),
	88	fCPVLocMaxCut(0.03),
	89	fCPVMinE(0.0),
	90	fCPVW0(4.0),
	91	fCPVUnfold(kTRUE)
	92	{
	93	//Copy constructor.
	94	}
	95
	96	//-----------------------------------------------------------------------------
	97	AliPHOSRecoParam& AliPHOSRecoParam::operator = (const AliPHOSRecoParam& recoParam)
	98	{
	99	//Assignment operator.
	100
	101	if(this != &recoParam) {
	102	fNonLinearityParams = recoParam.fNonLinearityParams;
	103	fEMCClusteringThreshold = recoParam.fEMCClusteringThreshold;
	104	fEMCLocMaxCut = recoParam.fEMCLocMaxCut;
	105	fEMCRawDigitThreshold = recoParam.fEMCRawDigitThreshold;
	106	fEMCMinE = recoParam.fEMCMinE;
	107	fEMCW0 = recoParam.fEMCW0;
	108	fEMCSampleQualityCut = recoParam.fEMCSampleQualityCut;
	109	fTimeGateAmpThresh = recoParam.fTimeGateAmpThresh ;
	110	fTimeGateLow = recoParam.fTimeGateLow ;
	111	fTimeGateHigh = recoParam.fTimeGateHigh ;
	112	fEMCEcoreRadius = recoParam.fEMCEcoreRadius;
	113	fEMCEcore2ESD = recoParam.fEMCEcore2ESD;
	114	fEMCSubtractPedestals = recoParam.fEMCSubtractPedestals;
	115	fEMCUnfold = recoParam.fEMCUnfold;
	116	fEMCEnergyCorrectionOn = recoParam.fEMCEnergyCorrectionOn;
	117	fEMCFitterVersion = recoParam.fEMCFitterVersion;
	118	fNonlinearityCorrVersion= recoParam.fNonlinearityCorrVersion;
	119	fGlobalAltroOffset = recoParam.fGlobalAltroOffset;
	120	fGlobalAltroThreshold = recoParam.fGlobalAltroThreshold;
	121	fCPVClusteringThreshold = recoParam.fCPVClusteringThreshold;
	122	fCPVLocMaxCut = recoParam.fCPVLocMaxCut;
	123	fCPVMinE = recoParam.fCPVMinE;
	124	fCPVW0 = recoParam.fCPVW0;
	125	fCPVUnfold = recoParam.fCPVUnfold;
	126	fNonLinearityParams = recoParam.fNonLinearityParams ;
	127	}
	128
	129	return *this;
	130	}
	131
	132	//-----------------------------------------------------------------------------
	133	void AliPHOSRecoParam::Print(Option_t * /option/) const
	134	{
	135	AliDebug(2,Form("PHOS reconstruction parameters:\n"
	136	"\tEMCClusteringThreshold = %f\n"
	137	"\tEMCLocMaxCut = %f\n"
	138	"\tEMCRawDigitThreshold = %f\n"
	139	"\tEMCMinE = %f\n"
	140	"\tEMCW0 = %f\n"
	141	"\tEMCSampleQualityCut = %f\n"
	142	"\tEMCEcoreRadius = %f\n"
	143	"\tEMCEcore2ESD = %d\n"
	144	"\tEMCSubtractPedestals = %d\n"
	145	"\tEMCUnfold = %d\n"
	146	"\tEMCEnergyCorrectionOn = %d\n"
	147	"\tEMCFitterVersion = \"%s\"\n"
	148	"\tGlobalAltroOffset = %d\n"
	149	"\tGlobalAltroThreshold = %d",
	150	fEMCClusteringThreshold,
	151	fEMCLocMaxCut,
	152	fEMCRawDigitThreshold,
	153	fEMCMinE,
	154	fEMCW0,
	155	fEMCSampleQualityCut,
	156	fEMCEcoreRadius,
	157	fEMCEcore2ESD,
	158	fEMCSubtractPedestals,
	159	fEMCUnfold,
	160	fEMCEnergyCorrectionOn,
	161	fEMCFitterVersion.Data(),
	162	fGlobalAltroOffset,
	163	fGlobalAltroThreshold));
	164
	165	}
	166
	167	//-----------------------------------------------------------------------------
	168	AliPHOSRecoParam* AliPHOSRecoParam::GetDefaultParameters()
	169	{
	170	//Default parameters for the reconstruction
	171
	172	AliPHOSRecoParam* params = new AliPHOSRecoParam();
	173	return params;
	174	}
	175
	176	//-----------------------------------------------------------------------------
	177	const TObjArray* AliPHOSRecoParam::GetMappings()
	178	{
	179	//Returns array of AliAltroMappings for RCU0..RCU3.
	180	//If not found, read it from OCDB.
	181
	182	//Quick check as follows:
	183	// root [0] AliCDBManager::Instance()->SetDefaultStorage("local://$ALICE_ROOT/OCDB");
	184	// root [1] AliCDBManager::Instance()->SetRun(1);
	185	// root [2] TObjArray* maps = AliPHOSRecoParam::GetMappings();
	186	// root [3] maps->Print();
	187
	188	if(fgkMaps) return fgkMaps;
	189
	190	AliCDBEntry* entry = AliCDBManager::Instance()->Get("PHOS/Calib/Mapping");
	191	if(entry)
	192	fgkMaps = (TObjArray*)entry->GetObject();
	193
	194	return fgkMaps;
	195
	196	}
	197	//-----------------------------------------------------------------------------
	198	void AliPHOSRecoParam::SetNonlinearityCorrectionVersion(const char * ver){
	199	//Set non-linearity correction version and adjust parameters
	200
	201	if(strcmp(ver,"NoCorrection")==0){
	202	fNonLinearityParams.Set(1) ;
	203	fNonLinearityParams.AddAt(0.,0) ;
	204	fNonlinearityCorrVersion=ver ;
	205	return ;
	206	}
	207	if(strcmp(ver,"Henrik2010")==0){
	208	Float_t pars[7]={1.051,2.54e-003,-1.737e-002,5.69e-002,3.3e-001,1.42e-001,1.50e-002} ;
	209	fNonLinearityParams.Set(7,pars) ;
	210	fNonlinearityCorrVersion=ver ;
	211	return ;
	212	}
	213	if(strcmp(ver,"Gustavo2005")==0){
	214	Float_t pars[3]={0.0241, 1.0504, 0.000249} ;
	215	fNonLinearityParams.Set(3,pars) ;
	216	fNonlinearityCorrVersion=ver ;
	217	return ;
	218	}
	219	AliError(Form("Non known correction version: %s, still using default \n",ver)) ;
	220
	221
	222	}
	223
	224	//-----------------------------------------------------------------------------
	225	const TObject* AliPHOSRecoParam::GetTriggerParameters()
	226	{
	227	//Returns trigger parameters.
	228
	229	//Quick check as follows:
	230	// root [0] AliCDBManager::Instance()->SetDefaultStorage("local://$ALICE_ROOT/OCDB");
	231	// root [1] AliCDBManager::Instance()->SetRun(158171);
	232	// root [2] TObject* parameters = AliPHOSRecoParam::GetTriggerParameters();
	233	// root [3] parameters->Print();
	234
	235	if(fgkTrigParams) return fgkTrigParams;
	236
	237	AliCDBEntry* entry = AliCDBManager::Instance()->Get("PHOS/Trigger/Parameters");
	238	if(entry)
	239	fgkTrigParams = entry->GetObject();
	240
	241	return fgkTrigParams;
	242
	243	}