[u/mrichter/AliRoot.git] / MUON / AliMUONRecoParam.cxx

/**************************************************************************
* Copyright(c) 1998-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.                  *
**************************************************************************/


//-----------------------------------------------------------------------------
/// \class AliMUONRecoParam
///
/// Class with MUON reconstruction parameters
///
///  \author Philippe Pillot
//-----------------------------------------------------------------------------


#include "AliMUONRecoParam.h"

#include "AliLog.h"

#include <Riostream.h>

ClassImp(AliMUONRecoParam)


//_____________________________________________________________________________
AliMUONRecoParam::AliMUONRecoParam()
: AliDetectorRecoParam(),
  fClusteringMode("MLEM"),
  fTrackingMode("KALMAN"),
  fMinBendingMomentum(0.),
  fMaxBendingMomentum(0.),
  fNonBendingVertexDispersion(0.),
  fBendingVertexDispersion(0.),
  fMaxNonBendingDistanceToTrack(0.),
  fMaxBendingDistanceToTrack(0.),
  fSigmaCutForTracking(0.),
  fSigmaCutForImprovement(0.),
  fSigmaCutForTrigger(0.),
  fMaxNormChi2MatchTrigger(0.),
  fCombinedClusterTrackReco(kFALSE),
  fTrackAllTracks(kFALSE),
  fRecoverTracks(kFALSE),
  fMakeTrackCandidatesFast(kFALSE),
  fComplementTracks(kFALSE),
  fImproveTracks(kFALSE),
  fUseSmoother(kFALSE)
{
  /// Constructor
  SetNameTitle("MUON","MUON");
  
  // use the default parameters for low flux environment
  SetLowFluxParam();
}

//_____________________________________________________________________________
AliMUONRecoParam::~AliMUONRecoParam() 
{
  /// Destructor
}

//_____________________________________________________________________________
AliMUONRecoParam *AliMUONRecoParam::GetLowFluxParam() 
{
  /// Return default reconstruction parameters for low flux environment
  
  AliMUONRecoParam *param = new AliMUONRecoParam();
  param->SetLowFluxParam();
  
  return param;
}

//_____________________________________________________________________________
AliMUONRecoParam *AliMUONRecoParam::GetHighFluxParam() 
{
  /// Return default reconstruction parameters for high flux environment
  
  AliMUONRecoParam *param = new AliMUONRecoParam();
  param->SetHighFluxParam();
 
  return param;
}

//_____________________________________________________________________________
void AliMUONRecoParam::SetLowFluxParam() 
{
  /// Set reconstruction parameters for low flux environment
  
  fMinBendingMomentum = 0.5;
  fMaxBendingMomentum = 3000.;
  fNonBendingVertexDispersion = 10.;
  fBendingVertexDispersion = 10.;
  fMaxNonBendingDistanceToTrack = 2.;
  fMaxBendingDistanceToTrack = 2.;
  fSigmaCutForTracking = 6.;
  fSigmaCutForImprovement = 5.;
  fSigmaCutForTrigger = 8.;
  fMaxNormChi2MatchTrigger = 16.;
  fTrackAllTracks = kTRUE;
  fRecoverTracks = kTRUE;
  fMakeTrackCandidatesFast = kFALSE;
  fComplementTracks = kTRUE;
  fImproveTracks = kTRUE;
  fUseSmoother = kTRUE;
  
}

//_____________________________________________________________________________
void AliMUONRecoParam::SetHighFluxParam() 
{
  /// Set reconstruction parameters for high flux environment
  
  fMinBendingMomentum = 0.5;
  fMaxBendingMomentum = 3000.;
  fNonBendingVertexDispersion = 10.;
  fBendingVertexDispersion = 10.;
  fMaxNonBendingDistanceToTrack = 2.;
  fMaxBendingDistanceToTrack = 2.;
  fSigmaCutForTracking = 6.;
  fSigmaCutForImprovement = 5.;
  fSigmaCutForTrigger = 8.;
  fMaxNormChi2MatchTrigger = 16.;
  fTrackAllTracks = kTRUE;
  fRecoverTracks = kTRUE;
  fMakeTrackCandidatesFast = kFALSE;
  fComplementTracks = kTRUE;
  fImproveTracks = kTRUE;
  fUseSmoother = kTRUE;
  
}

//_____________________________________________________________________________
void AliMUONRecoParam::Print(Option_t *option) const
{
  /// print reconstruction parameters
  /// if option = FULL then print also unused parameters
  
  cout<<endl<<"\t------Reconstruction parameters------"<<endl;
  
  cout<<Form("Clustering mode = %s",fClusteringMode.Data())<<endl;
  
  cout<<Form("Tracking mode = %s",fTrackingMode.Data())<<endl;
  
  if (fCombinedClusterTrackReco) cout<<"Combined cluster/track reconstruction: ON"<<endl;
  else cout<<"Combined cluster/track reconstruction: OFF"<<endl;
  
  cout<<Form("Bending momentum range = [%5.2f,%5.2f]",fMinBendingMomentum,fMaxBendingMomentum)<<endl;
  
  if (strstr(fTrackingMode,"ORIGINAL"))
    cout<<Form("Vertex dispertion = (%5.2f,%5.2f)",fNonBendingVertexDispersion,fBendingVertexDispersion)<<endl;
  else if (strstr(option,"FULL"))
    cout<<Form("Vertex dispertion (used for original tracking only) = (%5.2f,%5.2f)",fNonBendingVertexDispersion,fBendingVertexDispersion)<<endl;
  
  cout<<Form("Maximum distance to track = (%5.2f,%5.2f)",fMaxNonBendingDistanceToTrack,fMaxBendingDistanceToTrack)<<endl;
  
  cout<<Form("Sigma cut for tracking = %5.2f",fSigmaCutForTracking)<<endl;

  cout<<Form("Sigma cut for trigger hit pattern = %5.2f",fSigmaCutForTrigger)<<endl;
  
  if (fTrackAllTracks) cout<<"Track all the possible candidates"<<endl;
  else cout<<"Track only the best candidates"<<endl;
  
  if (strstr(option,"FULL")) {
    cout<<"Make track candidates assuming linear propagation between stations 4 and 5: ";
    if (fMakeTrackCandidatesFast) cout<<"ON"<<endl;
    else cout<<"OFF"<<endl;
  } else if (fMakeTrackCandidatesFast)
    cout<<"Make track candidates assuming linear propagation between stations 4 and 5"<<endl;
  
  if (strstr(option,"FULL")) {
    cout<<"Try to recover tracks getting lost during tracking: ";
    if (fRecoverTracks) cout<<"ON"<<endl;
    else cout<<"OFF"<<endl;
  } else if (fRecoverTracks)
    cout<<"Try to recover tracks getting lost during tracking"<<endl;
  
  if (strstr(option,"FULL")) {
    cout<<"Try to complete the reconstructed tracks by adding missing clusters: ";
    if (fComplementTracks) cout<<"ON"<<endl;
    else cout<<"OFF"<<endl;
  } else if (fComplementTracks)
    cout<<"Try to complete the reconstructed tracks by adding missing clusters"<<endl;
  
  if (strstr(option,"FULL")) {
    cout<<"Try to improve the reconstructed tracks by removing bad clusters: ";
    if (fImproveTracks) cout<<Form("ON (sigma cut = %5.2f)",fSigmaCutForImprovement)<<endl;
    else cout<<"OFF"<<endl;
  } else if (fImproveTracks)
    cout<<Form("Try to improve the reconstructed tracks by removing bad clusters (sigma cut = %5.2f)",fSigmaCutForImprovement)<<endl;
  
  if (strstr(option,"FULL")) {
    cout<<"Use smoother to compute final track parameters, etc, at each cluster (used for Kalman tracking only): ";
    if (fUseSmoother) cout<<"ON"<<endl;
    else cout<<"OFF"<<endl;
  } else if (fUseSmoother)
    cout<<"Use smoother to compute final track parameters, etc, at each cluster"<<endl;
  
  cout<<Form("Maximum normalized chi2 of tracking/trigger track matching = %5.2f",fMaxNormChi2MatchTrigger)<<endl;
  
  cout<<"\t-------------------------------------"<<endl<<endl;
  
}
Commit	Line	Data
3304fa09	1	/**************************************************************************
	2	* Copyright(c) 1998-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
	17	//-----------------------------------------------------------------------------
	18	/// \class AliMUONRecoParam
	19	///
	20	/// Class with MUON reconstruction parameters
	21	///
	22	/// \author Philippe Pillot
	23	//-----------------------------------------------------------------------------
	24
	25
	26
	27	#include "AliMUONRecoParam.h"
	28
	29	#include "AliLog.h"
	30
	31	#include <Riostream.h>
	32
	33	ClassImp(AliMUONRecoParam)
	34
	35
	36	//_____________________________________________________________________________
	37	AliMUONRecoParam::AliMUONRecoParam()
15d30ed4	38	: AliDetectorRecoParam(),
3304fa09	39	fClusteringMode("MLEM"),
	40	fTrackingMode("KALMAN"),
	41	fMinBendingMomentum(0.),
	42	fMaxBendingMomentum(0.),
	43	fNonBendingVertexDispersion(0.),
	44	fBendingVertexDispersion(0.),
	45	fMaxNonBendingDistanceToTrack(0.),
	46	fMaxBendingDistanceToTrack(0.),
	47	fSigmaCutForTracking(0.),
	48	fSigmaCutForImprovement(0.),
fda59e58	49	fSigmaCutForTrigger(0.),
3304fa09	50	fMaxNormChi2MatchTrigger(0.),
38bcf0ef	51	fCombinedClusterTrackReco(kFALSE),
3304fa09	52	fTrackAllTracks(kFALSE),
	53	fRecoverTracks(kFALSE),
	54	fMakeTrackCandidatesFast(kFALSE),
	55	fComplementTracks(kFALSE),
	56	fImproveTracks(kFALSE),
	57	fUseSmoother(kFALSE)
	58	{
	59	/// Constructor
15d30ed4	60	SetNameTitle("MUON","MUON");
3304fa09	61
	62	// use the default parameters for low flux environment
	63	SetLowFluxParam();
	64	}
	65
	66	//_____________________________________________________________________________
	67	AliMUONRecoParam::~AliMUONRecoParam()
	68	{
	69	/// Destructor
	70	}
	71
	72	//_____________________________________________________________________________
	73	AliMUONRecoParam *AliMUONRecoParam::GetLowFluxParam()
	74	{
	75	/// Return default reconstruction parameters for low flux environment
	76
	77	AliMUONRecoParam *param = new AliMUONRecoParam();
	78	param->SetLowFluxParam();
	79
	80	return param;
	81	}
	82
	83	//_____________________________________________________________________________
	84	AliMUONRecoParam *AliMUONRecoParam::GetHighFluxParam()
	85	{
	86	/// Return default reconstruction parameters for high flux environment
	87
	88	AliMUONRecoParam *param = new AliMUONRecoParam();
	89	param->SetHighFluxParam();
	90
	91	return param;
	92	}
	93
	94	//_____________________________________________________________________________
	95	void AliMUONRecoParam::SetLowFluxParam()
	96	{
	97	/// Set reconstruction parameters for low flux environment
	98
96ebe67e	99	fMinBendingMomentum = 0.5;
3304fa09	100	fMaxBendingMomentum = 3000.;
	101	fNonBendingVertexDispersion = 10.;
	102	fBendingVertexDispersion = 10.;
	103	fMaxNonBendingDistanceToTrack = 2.;
	104	fMaxBendingDistanceToTrack = 2.;
	105	fSigmaCutForTracking = 6.;
	106	fSigmaCutForImprovement = 5.;
fda59e58	107	fSigmaCutForTrigger = 8.;
3304fa09	108	fMaxNormChi2MatchTrigger = 16.;
	109	fTrackAllTracks = kTRUE;
	110	fRecoverTracks = kTRUE;
	111	fMakeTrackCandidatesFast = kFALSE;
	112	fComplementTracks = kTRUE;
	113	fImproveTracks = kTRUE;
	114	fUseSmoother = kTRUE;
	115
	116	}
	117
	118	//_____________________________________________________________________________
	119	void AliMUONRecoParam::SetHighFluxParam()
	120	{
	121	/// Set reconstruction parameters for high flux environment
	122
96ebe67e	123	fMinBendingMomentum = 0.5;
3304fa09	124	fMaxBendingMomentum = 3000.;
	125	fNonBendingVertexDispersion = 10.;
	126	fBendingVertexDispersion = 10.;
	127	fMaxNonBendingDistanceToTrack = 2.;
	128	fMaxBendingDistanceToTrack = 2.;
	129	fSigmaCutForTracking = 6.;
	130	fSigmaCutForImprovement = 5.;
fda59e58	131	fSigmaCutForTrigger = 8.;
3304fa09	132	fMaxNormChi2MatchTrigger = 16.;
	133	fTrackAllTracks = kTRUE;
	134	fRecoverTracks = kTRUE;
	135	fMakeTrackCandidatesFast = kFALSE;
	136	fComplementTracks = kTRUE;
	137	fImproveTracks = kTRUE;
	138	fUseSmoother = kTRUE;
	139
	140	}
	141
	142	//_____________________________________________________________________________
	143	void AliMUONRecoParam::Print(Option_t *option) const
	144	{
	145	/// print reconstruction parameters
	146	/// if option = FULL then print also unused parameters
	147
	148	cout<<endl<<"\t------Reconstruction parameters------"<<endl;
	149
	150	cout<<Form("Clustering mode = %s",fClusteringMode.Data())<<endl;
	151
	152	cout<<Form("Tracking mode = %s",fTrackingMode.Data())<<endl;
	153
15d30ed4	154	if (fCombinedClusterTrackReco) cout<<"Combined cluster/track reconstruction: ON"<<endl;
15d30ed4	155	else cout<<"Combined cluster/track reconstruction: OFF"<<endl;
38bcf0ef	156
3304fa09	157	cout<<Form("Bending momentum range = [%5.2f,%5.2f]",fMinBendingMomentum,fMaxBendingMomentum)<<endl;
	158
	159	if (strstr(fTrackingMode,"ORIGINAL"))
	160	cout<<Form("Vertex dispertion = (%5.2f,%5.2f)",fNonBendingVertexDispersion,fBendingVertexDispersion)<<endl;
	161	else if (strstr(option,"FULL"))
	162	cout<<Form("Vertex dispertion (used for original tracking only) = (%5.2f,%5.2f)",fNonBendingVertexDispersion,fBendingVertexDispersion)<<endl;
	163
	164	cout<<Form("Maximum distance to track = (%5.2f,%5.2f)",fMaxNonBendingDistanceToTrack,fMaxBendingDistanceToTrack)<<endl;
	165
	166	cout<<Form("Sigma cut for tracking = %5.2f",fSigmaCutForTracking)<<endl;
fda59e58	167
fda59e58	168	cout<<Form("Sigma cut for trigger hit pattern = %5.2f",fSigmaCutForTrigger)<<endl;
3304fa09	169
	170	if (fTrackAllTracks) cout<<"Track all the possible candidates"<<endl;
	171	else cout<<"Track only the best candidates"<<endl;
	172
	173	if (strstr(option,"FULL")) {
	174	cout<<"Make track candidates assuming linear propagation between stations 4 and 5: ";
	175	if (fMakeTrackCandidatesFast) cout<<"ON"<<endl;
	176	else cout<<"OFF"<<endl;
	177	} else if (fMakeTrackCandidatesFast)
	178	cout<<"Make track candidates assuming linear propagation between stations 4 and 5"<<endl;
	179
	180	if (strstr(option,"FULL")) {
	181	cout<<"Try to recover tracks getting lost during tracking: ";
	182	if (fRecoverTracks) cout<<"ON"<<endl;
	183	else cout<<"OFF"<<endl;
	184	} else if (fRecoverTracks)
	185	cout<<"Try to recover tracks getting lost during tracking"<<endl;
	186
	187	if (strstr(option,"FULL")) {
	188	cout<<"Try to complete the reconstructed tracks by adding missing clusters: ";
	189	if (fComplementTracks) cout<<"ON"<<endl;
	190	else cout<<"OFF"<<endl;
	191	} else if (fComplementTracks)
	192	cout<<"Try to complete the reconstructed tracks by adding missing clusters"<<endl;
	193
	194	if (strstr(option,"FULL")) {
	195	cout<<"Try to improve the reconstructed tracks by removing bad clusters: ";
	196	if (fImproveTracks) cout<<Form("ON (sigma cut = %5.2f)",fSigmaCutForImprovement)<<endl;
	197	else cout<<"OFF"<<endl;
	198	} else if (fImproveTracks)
	199	cout<<Form("Try to improve the reconstructed tracks by removing bad clusters (sigma cut = %5.2f)",fSigmaCutForImprovement)<<endl;
	200
	201	if (strstr(option,"FULL")) {
	202	cout<<"Use smoother to compute final track parameters, etc, at each cluster (used for Kalman tracking only): ";
	203	if (fUseSmoother) cout<<"ON"<<endl;
	204	else cout<<"OFF"<<endl;
	205	} else if (fUseSmoother)
	206	cout<<"Use smoother to compute final track parameters, etc, at each cluster"<<endl;
	207
	208	cout<<Form("Maximum normalized chi2 of tracking/trigger track matching = %5.2f",fMaxNormChi2MatchTrigger)<<endl;
	209
	210	cout<<"\t-------------------------------------"<<endl<<endl;
	211
	212	}
	213