[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.),
  fMaxNonBendingSlope(0.),
  fNonBendingVertexDispersion(0.),
  fBendingVertexDispersion(0.),
  fMaxNonBendingDistanceToTrack(0.),
  fMaxBendingDistanceToTrack(0.),
  fSigmaCutForTracking(0.),
  fSigmaCutForImprovement(0.),
  fSigmaCutForTrigger(0.),
  fMaxNormChi2MatchTrigger(0.),
  fPercentOfFullClusterInESD(10.),
  fCombinedClusterTrackReco(kFALSE),
  fTrackAllTracks(kFALSE),
  fRecoverTracks(kFALSE),
  fMakeTrackCandidatesFast(kFALSE),
  fMakeMoreTrackCandidates(kFALSE),
  fComplementTracks(kFALSE),
  fImproveTracks(kFALSE),
  fUseSmoother(kFALSE),
  fSaveFullClusterInESD(kTRUE),
  fCalibrationMode("NOGAIN"),
  fBypassSt45(kFALSE)
{
  /// Constructor
  SetNameTitle("MUON","MUON");
  
  // use the default parameters for low flux environment
  SetLowFluxParam();
}

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

//_____________________________________________________________________________
Option_t*
AliMUONRecoParam::GetCalibrationMode() const
{
  /// Return the calibration mode. Can be : 
  /// NOGAIN : only do pedestal subtraction
  /// GAIN : do pedestal subtraction, and apply gain correction, but with a
  ///        single capacitance value for all channels
  /// GAINCONSTANTCAPA : as GAIN, but with a channel-dependent capacitance value
  
  return fCalibrationMode.Data();
}

//_____________________________________________________________________________
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 = 1.;
  fMaxBendingMomentum = 3000.;
  fMaxNonBendingSlope = 0.3;
  fNonBendingVertexDispersion = 10.;
  fBendingVertexDispersion = 10.;
  fMaxNonBendingDistanceToTrack = 1.;
  fMaxBendingDistanceToTrack = 1.;
  fSigmaCutForTracking = 6.;
  fSigmaCutForImprovement = 5.;
  fSigmaCutForTrigger = 8.;
  fMaxNormChi2MatchTrigger = 16.;
  fTrackAllTracks = kTRUE;
  fRecoverTracks = kTRUE;
  fMakeTrackCandidatesFast = kFALSE;
  fMakeMoreTrackCandidates = kFALSE;
  fComplementTracks = kTRUE;
  fImproveTracks = kTRUE;
  fUseSmoother = kTRUE;
  for (Int_t iCh = 0; iCh < 10; iCh++) fUseChamber[iCh] = kTRUE;
  for (Int_t iSt = 0; iSt < 5; iSt++) fRequestStation[iSt] = kTRUE;
  
}

//_____________________________________________________________________________
void AliMUONRecoParam::SetHighFluxParam() 
{
  /// Set reconstruction parameters for high flux environment
  
  fMinBendingMomentum = 1.;
  fMaxBendingMomentum = 3000.;
  fMaxNonBendingSlope = 0.3;
  fNonBendingVertexDispersion = 10.;
  fBendingVertexDispersion = 10.;
  fMaxNonBendingDistanceToTrack = 1.;
  fMaxBendingDistanceToTrack = 1.;
  fSigmaCutForTracking = 6.;
  fSigmaCutForImprovement = 5.;
  fSigmaCutForTrigger = 8.;
  fMaxNormChi2MatchTrigger = 16.;
  fTrackAllTracks = kTRUE;
  fRecoverTracks = kTRUE;
  fMakeTrackCandidatesFast = kFALSE;
  fMakeMoreTrackCandidates = kFALSE;
  fComplementTracks = kTRUE;
  fImproveTracks = kTRUE;
  fUseSmoother = kTRUE;
  for (Int_t iCh = 0; iCh < 10; iCh++) fUseChamber[iCh] = kTRUE;
  for (Int_t iSt = 0; iSt < 5; iSt++) fRequestStation[iSt] = 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("Calibration mode = %s",fCalibrationMode.Data())<<endl;
  cout<<Form("Clustering mode = %s",fClusteringMode.Data())<<endl;
  cout<<Form("Tracking mode = %s",fTrackingMode.Data())<<endl;

  if (BypassSt45()) cout << "Will bypass St45 (replacing their clusters by generated ones from trigger tracks)" << endl;
  
  if (fCombinedClusterTrackReco) cout<<"Combined cluster/track reconstruction: ON"<<endl;
  else cout<<"Combined cluster/track reconstruction: OFF"<<endl;
  
  if (fSaveFullClusterInESD) cout<<Form("Save all cluster info in ESD for %5.2f %% of events",fPercentOfFullClusterInESD)<<endl;
  else cout<<"Save partial cluster info in ESD"<<endl;
  
  cout<<Form("Bending momentum range = [%5.2f,%5.2f]",fMinBendingMomentum,fMaxBendingMomentum)<<endl;
  
  cout<<Form("Maximum non bending slope = %5.2f",fMaxNonBendingSlope)<<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<<"Make track candidates starting from 1 cluster in each of the stations 4 and 5: ";
    if (fMakeMoreTrackCandidates) cout<<"ON"<<endl;
    else cout<<"OFF"<<endl;
  } else if (fMakeMoreTrackCandidates)
    cout<<"Make track candidates starting from 1 cluster in each of the 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;
  
  Bool_t discardedCh = kFALSE;
  Int_t ch = 0;
  do {
    if (!UseChamber(ch)) {
      if (!discardedCh) {
	cout<<"Discarded chambers(1..): "<<ch+1;
	discardedCh = kTRUE;
      }
      else cout<<" "<<ch+1;
    }
  } while (++ch < 10);
  if (discardedCh) cout<<endl;
  
  Bool_t discardedSt = kFALSE;
  Int_t st = 0;
  do {
    if (!RequestStation(st)) {
      if (!discardedSt) {
	cout<<"Not requested stations(1..): "<<st+1;
	discardedSt = kTRUE;
      }
      else cout<<" "<<st+1;
    }
  } while (++st < 5);
  if (discardedSt) cout<<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.),
9bf6860b	43	fMaxNonBendingSlope(0.),
3304fa09	44	fNonBendingVertexDispersion(0.),
	45	fBendingVertexDispersion(0.),
	46	fMaxNonBendingDistanceToTrack(0.),
	47	fMaxBendingDistanceToTrack(0.),
	48	fSigmaCutForTracking(0.),
	49	fSigmaCutForImprovement(0.),
fda59e58	50	fSigmaCutForTrigger(0.),
3304fa09	51	fMaxNormChi2MatchTrigger(0.),
0a18ba02	52	fPercentOfFullClusterInESD(10.),
38bcf0ef	53	fCombinedClusterTrackReco(kFALSE),
3304fa09	54	fTrackAllTracks(kFALSE),
	55	fRecoverTracks(kFALSE),
	56	fMakeTrackCandidatesFast(kFALSE),
9bf6860b	57	fMakeMoreTrackCandidates(kFALSE),
3304fa09	58	fComplementTracks(kFALSE),
3304fa09	59	fImproveTracks(kFALSE),
0a18ba02	60	fUseSmoother(kFALSE),
de98fdc9	61	fSaveFullClusterInESD(kTRUE),
9bf6860b	62	fCalibrationMode("NOGAIN"),
9bf6860b	63	fBypassSt45(kFALSE)
3304fa09	64	{
3304fa09	65	/// Constructor
15d30ed4	66	SetNameTitle("MUON","MUON");
3304fa09	67
	68	// use the default parameters for low flux environment
	69	SetLowFluxParam();
	70	}
	71
	72	//_____________________________________________________________________________
	73	AliMUONRecoParam::~AliMUONRecoParam()
	74	{
	75	/// Destructor
	76	}
	77
9bf6860b	78	//_____________________________________________________________________________
	79	Option_t*
	80	AliMUONRecoParam::GetCalibrationMode() const
	81	{
	82	/// Return the calibration mode. Can be :
	83	/// NOGAIN : only do pedestal subtraction
	84	/// GAIN : do pedestal subtraction, and apply gain correction, but with a
	85	/// single capacitance value for all channels
	86	/// GAINCONSTANTCAPA : as GAIN, but with a channel-dependent capacitance value
	87
	88	return fCalibrationMode.Data();
	89	}
	90
3304fa09	91	//_____________________________________________________________________________
	92	AliMUONRecoParam *AliMUONRecoParam::GetLowFluxParam()
	93	{
	94	/// Return default reconstruction parameters for low flux environment
	95
	96	AliMUONRecoParam *param = new AliMUONRecoParam();
	97	param->SetLowFluxParam();
	98
	99	return param;
	100	}
	101
	102	//_____________________________________________________________________________
	103	AliMUONRecoParam *AliMUONRecoParam::GetHighFluxParam()
	104	{
	105	/// Return default reconstruction parameters for high flux environment
	106
	107	AliMUONRecoParam *param = new AliMUONRecoParam();
	108	param->SetHighFluxParam();
	109
	110	return param;
	111	}
	112
	113	//_____________________________________________________________________________
	114	void AliMUONRecoParam::SetLowFluxParam()
	115	{
	116	/// Set reconstruction parameters for low flux environment
	117
9bf6860b	118	fMinBendingMomentum = 1.;
3304fa09	119	fMaxBendingMomentum = 3000.;
9bf6860b	120	fMaxNonBendingSlope = 0.3;
3304fa09	121	fNonBendingVertexDispersion = 10.;
3304fa09	122	fBendingVertexDispersion = 10.;
9bf6860b	123	fMaxNonBendingDistanceToTrack = 1.;
9bf6860b	124	fMaxBendingDistanceToTrack = 1.;
3304fa09	125	fSigmaCutForTracking = 6.;
3304fa09	126	fSigmaCutForImprovement = 5.;
fda59e58	127	fSigmaCutForTrigger = 8.;
3304fa09	128	fMaxNormChi2MatchTrigger = 16.;
	129	fTrackAllTracks = kTRUE;
	130	fRecoverTracks = kTRUE;
	131	fMakeTrackCandidatesFast = kFALSE;
9bf6860b	132	fMakeMoreTrackCandidates = kFALSE;
3304fa09	133	fComplementTracks = kTRUE;
	134	fImproveTracks = kTRUE;
	135	fUseSmoother = kTRUE;
9bf6860b	136	for (Int_t iCh = 0; iCh < 10; iCh++) fUseChamber[iCh] = kTRUE;
9bf6860b	137	for (Int_t iSt = 0; iSt < 5; iSt++) fRequestStation[iSt] = kTRUE;
3304fa09	138
	139	}
	140
	141	//_____________________________________________________________________________
	142	void AliMUONRecoParam::SetHighFluxParam()
	143	{
	144	/// Set reconstruction parameters for high flux environment
	145
9bf6860b	146	fMinBendingMomentum = 1.;
3304fa09	147	fMaxBendingMomentum = 3000.;
9bf6860b	148	fMaxNonBendingSlope = 0.3;
3304fa09	149	fNonBendingVertexDispersion = 10.;
3304fa09	150	fBendingVertexDispersion = 10.;
9bf6860b	151	fMaxNonBendingDistanceToTrack = 1.;
9bf6860b	152	fMaxBendingDistanceToTrack = 1.;
3304fa09	153	fSigmaCutForTracking = 6.;
3304fa09	154	fSigmaCutForImprovement = 5.;
fda59e58	155	fSigmaCutForTrigger = 8.;
3304fa09	156	fMaxNormChi2MatchTrigger = 16.;
	157	fTrackAllTracks = kTRUE;
	158	fRecoverTracks = kTRUE;
	159	fMakeTrackCandidatesFast = kFALSE;
9bf6860b	160	fMakeMoreTrackCandidates = kFALSE;
3304fa09	161	fComplementTracks = kTRUE;
	162	fImproveTracks = kTRUE;
	163	fUseSmoother = kTRUE;
9bf6860b	164	for (Int_t iCh = 0; iCh < 10; iCh++) fUseChamber[iCh] = kTRUE;
9bf6860b	165	for (Int_t iSt = 0; iSt < 5; iSt++) fRequestStation[iSt] = kTRUE;
3304fa09	166
	167	}
	168
	169	//_____________________________________________________________________________
	170	void AliMUONRecoParam::Print(Option_t *option) const
	171	{
	172	/// print reconstruction parameters
	173	/// if option = FULL then print also unused parameters
	174
	175	cout<<endl<<"\t------Reconstruction parameters------"<<endl;
	176
de98fdc9	177	cout<<Form("Calibration mode = %s",fCalibrationMode.Data())<<endl;
3304fa09	178	cout<<Form("Clustering mode = %s",fClusteringMode.Data())<<endl;
3304fa09	179	cout<<Form("Tracking mode = %s",fTrackingMode.Data())<<endl;
9bf6860b	180
9bf6860b	181	if (BypassSt45()) cout << "Will bypass St45 (replacing their clusters by generated ones from trigger tracks)" << endl;
3304fa09	182
15d30ed4	183	if (fCombinedClusterTrackReco) cout<<"Combined cluster/track reconstruction: ON"<<endl;
15d30ed4	184	else cout<<"Combined cluster/track reconstruction: OFF"<<endl;
38bcf0ef	185
0a18ba02	186	if (fSaveFullClusterInESD) cout<<Form("Save all cluster info in ESD for %5.2f %% of events",fPercentOfFullClusterInESD)<<endl;
	187	else cout<<"Save partial cluster info in ESD"<<endl;
	188
3304fa09	189	cout<<Form("Bending momentum range = [%5.2f,%5.2f]",fMinBendingMomentum,fMaxBendingMomentum)<<endl;
3304fa09	190
9bf6860b	191	cout<<Form("Maximum non bending slope = %5.2f",fMaxNonBendingSlope)<<endl;
9bf6860b	192
3304fa09	193	if (strstr(fTrackingMode,"ORIGINAL"))
	194	cout<<Form("Vertex dispertion = (%5.2f,%5.2f)",fNonBendingVertexDispersion,fBendingVertexDispersion)<<endl;
	195	else if (strstr(option,"FULL"))
	196	cout<<Form("Vertex dispertion (used for original tracking only) = (%5.2f,%5.2f)",fNonBendingVertexDispersion,fBendingVertexDispersion)<<endl;
	197
	198	cout<<Form("Maximum distance to track = (%5.2f,%5.2f)",fMaxNonBendingDistanceToTrack,fMaxBendingDistanceToTrack)<<endl;
	199
	200	cout<<Form("Sigma cut for tracking = %5.2f",fSigmaCutForTracking)<<endl;
fda59e58	201
fda59e58	202	cout<<Form("Sigma cut for trigger hit pattern = %5.2f",fSigmaCutForTrigger)<<endl;
3304fa09	203
	204	if (fTrackAllTracks) cout<<"Track all the possible candidates"<<endl;
	205	else cout<<"Track only the best candidates"<<endl;
	206
	207	if (strstr(option,"FULL")) {
	208	cout<<"Make track candidates assuming linear propagation between stations 4 and 5: ";
	209	if (fMakeTrackCandidatesFast) cout<<"ON"<<endl;
	210	else cout<<"OFF"<<endl;
	211	} else if (fMakeTrackCandidatesFast)
	212	cout<<"Make track candidates assuming linear propagation between stations 4 and 5"<<endl;
	213
9bf6860b	214	if (strstr(option,"FULL")) {
	215	cout<<"Make track candidates starting from 1 cluster in each of the stations 4 and 5: ";
	216	if (fMakeMoreTrackCandidates) cout<<"ON"<<endl;
	217	else cout<<"OFF"<<endl;
	218	} else if (fMakeMoreTrackCandidates)
	219	cout<<"Make track candidates starting from 1 cluster in each of the stations 4 and 5"<<endl;
	220
3304fa09	221	if (strstr(option,"FULL")) {
	222	cout<<"Try to recover tracks getting lost during tracking: ";
	223	if (fRecoverTracks) cout<<"ON"<<endl;
	224	else cout<<"OFF"<<endl;
	225	} else if (fRecoverTracks)
	226	cout<<"Try to recover tracks getting lost during tracking"<<endl;
	227
	228	if (strstr(option,"FULL")) {
	229	cout<<"Try to complete the reconstructed tracks by adding missing clusters: ";
	230	if (fComplementTracks) cout<<"ON"<<endl;
	231	else cout<<"OFF"<<endl;
	232	} else if (fComplementTracks)
	233	cout<<"Try to complete the reconstructed tracks by adding missing clusters"<<endl;
	234
	235	if (strstr(option,"FULL")) {
	236	cout<<"Try to improve the reconstructed tracks by removing bad clusters: ";
	237	if (fImproveTracks) cout<<Form("ON (sigma cut = %5.2f)",fSigmaCutForImprovement)<<endl;
	238	else cout<<"OFF"<<endl;
	239	} else if (fImproveTracks)
	240	cout<<Form("Try to improve the reconstructed tracks by removing bad clusters (sigma cut = %5.2f)",fSigmaCutForImprovement)<<endl;
	241
	242	if (strstr(option,"FULL")) {
	243	cout<<"Use smoother to compute final track parameters, etc, at each cluster (used for Kalman tracking only): ";
	244	if (fUseSmoother) cout<<"ON"<<endl;
	245	else cout<<"OFF"<<endl;
	246	} else if (fUseSmoother)
	247	cout<<"Use smoother to compute final track parameters, etc, at each cluster"<<endl;
	248
	249	cout<<Form("Maximum normalized chi2 of tracking/trigger track matching = %5.2f",fMaxNormChi2MatchTrigger)<<endl;
	250
9bf6860b	251	Bool_t discardedCh = kFALSE;
	252	Int_t ch = 0;
	253	do {
	254	if (!UseChamber(ch)) {
	255	if (!discardedCh) {
	256	cout<<"Discarded chambers(1..): "<<ch+1;
	257	discardedCh = kTRUE;
	258	}
	259	else cout<<" "<<ch+1;
	260	}
	261	} while (++ch < 10);
	262	if (discardedCh) cout<<endl;
	263
	264	Bool_t discardedSt = kFALSE;
	265	Int_t st = 0;
	266	do {
	267	if (!RequestStation(st)) {
	268	if (!discardedSt) {
	269	cout<<"Not requested stations(1..): "<<st+1;
	270	discardedSt = kTRUE;
	271	}
	272	else cout<<" "<<st+1;
	273	}
	274	} while (++st < 5);
	275	if (discardedSt) cout<<endl;
	276
	277
3304fa09	278	cout<<"\t-------------------------------------"<<endl<<endl;
	279
	280	}
	281