[u/mrichter/AliRoot.git] / PWG / CaloTrackCorrBase / AliHistogramRanges.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.                  *
 **************************************************************************/

//_________________________________________________________________________
// Class containing histogram settings:
//    - Number of bins
//    - Min and max of range
//
//-- Author: Gustavo Conesa (LPSC-Grenoble)
//

#include "TMath.h"

#include "AliHistogramRanges.h"

ClassImp(AliHistogramRanges)


//_______________________________________________
AliHistogramRanges::AliHistogramRanges() : 
TObject(), 
fHistoPtBins(0),              fHistoPtMax(0.),             fHistoPtMin(0.),
fHistoPhiBins(0),             fHistoPhiMax(0.),            fHistoPhiMin(0.),
fHistoEtaBins(0),             fHistoEtaMax(0.),            fHistoEtaMin(0.),
fHistoDeltaPhiBins(0),        fHistoDeltaPhiMax(0.),       fHistoDeltaPhiMin(0.),
fHistoDeltaEtaBins(0),        fHistoDeltaEtaMax(0.),       fHistoDeltaEtaMin(0.),
fHistoMassBins(0),            fHistoMassMax(0.),           fHistoMassMin(0.),
fHistoAsymBins(0),            fHistoAsymMax(0.),           fHistoAsymMin(0.),
fHistoV0SBins(0),             fHistoV0SMax(0),             fHistoV0SMin(0),
fHistoV0MBins(0),             fHistoV0MMax(0),             fHistoV0MMin(0),
fHistoTrMBins(0),             fHistoTrMMax(0),             fHistoTrMMin(0),
fHistoFinePtBins(1000),       fHistoFinePtMax(5.),         fHistoFinePtMin(0.),
fHistoPOverEBins(100),        fHistoPOverEMax(100.),       fHistoPOverEMin(0.),
fHistodEdxBins(100),          fHistodEdxMax(100.),         fHistodEdxMin(0.),
fHistodRBins(100),            fHistodRMax(100.),           fHistodRMin(0.),
fHistoTimeBins(100),          fHistoTimeMax(100.),         fHistoTimeMin(0.),
fHistoNClusCellBins(100),     fHistoNClusCellMax(100),     fHistoNClusCellMin(0),
fHistoNCellsBins(100),        fHistoNCellsMax(100),        fHistoNCellsMin(0),
fHistoNClustersBins(100),     fHistoNClustersMax(100),     fHistoNClustersMin(0),
fHistoRatioBins(100),         fHistoRatioMax(100.),        fHistoRatioMin(0.),
fHistoVertexDistBins(100),    fHistoVertexDistMax(100.),   fHistoVertexDistMin(0.),
fHistoRBins(100),             fHistoRMax(1000),            fHistoRMin(-1000),
fHistoXBins(100),             fHistoXMax(1000),            fHistoXMin(-1000),
fHistoYBins(100),             fHistoYMax(1000),            fHistoYMin(-1000),
fHistoZBins(100),             fHistoZMax(1000),            fHistoZMin(-1000),
fHistoSSBins(0),              fHistoSSMax(0),              fHistoSSMin(0),
fHistoDiffTimeBins(0),        fHistoDiffTimeMax(0),        fHistoDiffTimeMin(0),
fHistoTrackResidualEtaBins(0),fHistoTrackResidualEtaMax(0),fHistoTrackResidualEtaMin(0),
fHistoTrackResidualPhiBins(0),fHistoTrackResidualPhiMax(0),fHistoTrackResidualPhiMin(0)
{
  //Default Ctor
  
  //Initialize parameters
  InitParameters();
  
}

//_______________________________________
void AliHistogramRanges::InitParameters()
{ 
  //Initialize histogram parameters
  
  fHistoPtBins           = 240 ;  fHistoPtMax           = 120   ; fHistoPtMin           = 0.  ;
  fHistoPhiBins          = 120 ;  fHistoPhiMax          = TMath::TwoPi(); fHistoPhiMin  = 0.  ;
  fHistoEtaBins          = 100 ;  fHistoEtaMax          = 1     ; fHistoEtaMin          = -1  ;
  fHistoDeltaPhiBins     = 310 ;  fHistoDeltaPhiMax     = 4.8   ; fHistoDeltaPhiMin     = -1.7 ;
  fHistoDeltaEtaBins     = 200 ;  fHistoDeltaEtaMax     = 2     ; fHistoDeltaEtaMin     = -2 ;
  fHistoMassBins         = 200 ;  fHistoMassMax         = 1.    ; fHistoMassMin         = 0.  ;
  fHistoAsymBins         = 10  ;  fHistoAsymMax         = 1.    ; fHistoAsymMin         = 0.  ;
  fHistoV0SBins          = 100 ;  fHistoV0SMax          = 10000 ; fHistoV0SMin          = 0   ;
  fHistoV0MBins          = 100;   fHistoV0MMax          = 10000 ; fHistoV0MMin          = 0   ;
  fHistoTrMBins          = 100 ;  fHistoTrMMax          = 2000  ; fHistoTrMMin          = 0   ;
  fHistoPOverEBins       = 100 ;  fHistoPOverEMax       = 10.   ; fHistoPOverEMin       = 0.  ;
  fHistodEdxBins         = 200 ;  fHistodEdxMax         = 400.  ; fHistodEdxMin         = 0.  ;  
  fHistodRBins           = 300 ;  fHistodRMax           = 3.15  ; fHistodRMin           = 0.  ;
  fHistoTimeBins         = 1000;  fHistoTimeMax         = 1.e3  ; fHistoTimeMin         = 0.  ;//ns
  fHistoNClusCellBins    = 200 ;  fHistoNClusCellMax    = 200   ; fHistoNClusCellMin    = 0   ;
  fHistoNCellsBins       = 300 ;  fHistoNCellsMax       = 300   ; fHistoNCellsMin       = 0   ;
  fHistoNClustersBins    = 50  ;  fHistoNClustersMax    = 50    ; fHistoNClustersMin    = 0   ;
  fHistoRatioBins        = 200 ;  fHistoRatioMax        = 2     ; fHistoRatioMin        = 0.  ;
  fHistoVertexDistBins   = 100 ;  fHistoVertexDistMax   = 500.  ; fHistoVertexDistMin   = 0.  ;
  fHistoRBins            = 100 ;  fHistoRMax            = 500   ; fHistoRMin            = -500;//cm
  fHistoXBins            = 100 ;  fHistoXMax            = 500   ; fHistoXMin            = -500;//cm
  fHistoYBins            = 100 ;  fHistoYMax            = 500   ; fHistoYMin            = -500;//cm
  fHistoZBins            = 100 ;  fHistoZMax            = 600   ; fHistoZMin            = -500;//cm
  fHistoSSBins           = 500 ;  fHistoSSMax           = 5     ; fHistoSSMin           = 0   ;  
  fHistoDiffTimeBins     = 800 ;  fHistoDiffTimeMax     = 400   ; fHistoDiffTimeMin     = -400;// ns
  
  fHistoTrackResidualEtaBins = 100 ; fHistoTrackResidualEtaMax = 0.15 ; fHistoTrackResidualEtaMin = -0.15;
  fHistoTrackResidualPhiBins = 100 ; fHistoTrackResidualPhiMax = 0.15 ; fHistoTrackResidualPhiMin = -0.15;

}

//________________________________________________________
void AliHistogramRanges::Print(const Option_t * /*opt*/) const
{
  //Print some relevant parameters set for the analysis
  	
  printf("Histograms: %3.1f < pT  < %3.1f,  Nbin = %d\n"             , fHistoPtMin,          fHistoPtMax,          fHistoPtBins);
  printf("Histograms: %3.1f < phi < %3.1f, Nbin = %d\n"              , fHistoPhiMin,         fHistoPhiMax,         fHistoPhiBins);
  printf("Histograms: %3.1f < eta < %3.1f, Nbin = %d\n"              , fHistoEtaMin,         fHistoEtaMax,         fHistoEtaBins);
  printf("Histograms: %3.1f < delta phi < %3.1f, Nbin = %d\n"        , fHistoDeltaPhiMin,    fHistoDeltaPhiMax,    fHistoDeltaPhiBins);
  printf("Histograms: %3.1f < delta eta < %3.1f, Nbin = %d\n"        , fHistoDeltaEtaMin,    fHistoDeltaEtaMax,    fHistoDeltaEtaBins);
  printf("Histograms: %3.1f < mass < %3.1f, Nbin = %d\n"             , fHistoMassMin,        fHistoMassMax,        fHistoMassBins);
  printf("Histograms: %3.1f < asymmetry < %3.1f, Nbin = %d\n"        , fHistoAsymMin,        fHistoAsymMax,        fHistoAsymBins);
  printf("Histograms: %d < V0 Signal < %d, Nbin = %d\n"              , fHistoV0SMin,         fHistoV0SMax,         fHistoV0SBins);
  printf("Histograms: %d < V0 Mult < %d, Nbin = %d\n"                , fHistoV0MMin,         fHistoV0MMax,         fHistoV0MBins);
  printf("Histograms: %d < Track Mult < %d, Nbin = %d\n"             , fHistoTrMMin,         fHistoTrMMax,         fHistoTrMBins);
  printf("Histograms: %3.1f < p/E  < %3.1f, Nbin = %d\n"             , fHistoPOverEMin,      fHistoPOverEMax,      fHistoPOverEBins);
  printf("Histograms: %3.1f < dEdx < %3.1f, Nbin = %d\n"             , fHistodEdxMin,        fHistodEdxMax,        fHistodEdxBins);
  printf("Histograms: %3.1f < dR (track cluster)< %3.1f, Nbin = %d\n", fHistodRMin,          fHistodRMax,          fHistodRBins);
  printf("Histograms: %3.1f < R=sqrt{x^2+y^2}   < %3.1f, Nbin = %d\n", fHistoRMin,           fHistoRMax,           fHistoRBins);
  printf("Histograms: %3.1f < X    < %3.1f, Nbin = %d\n"             , fHistoXMin,           fHistoXMax,           fHistoXBins);
  printf("Histograms: %3.1f < Y    < %3.1f, Nbin = %d\n"             , fHistoYMin,           fHistoYMax,           fHistoYBins);
  printf("Histograms: %3.1f < Z    < %3.1f, Nbin = %d\n"             , fHistoZMin,           fHistoZMax,           fHistoZBins);
  printf("Histograms: %g < Time < %g, Nbin = %d\n"                   , fHistoTimeMin,        fHistoTimeMax,        fHistoTimeBins);
  printf("Histograms: %d < N cells per cluster    < %d, Nbin = %d\n" , fHistoNClusCellMin,   fHistoNClusCellMax,   fHistoNClusCellBins);
  printf("Histograms: %d < N cells   < %d, Nbin = %d\n"              , fHistoNCellsMin,      fHistoNCellsMax,      fHistoNCellsBins);
  printf("Histograms: %d < N clusters   < %d, Nbin = %d\n"           , fHistoNClustersMin,   fHistoNClustersMax,   fHistoNClustersBins);
  printf("Histograms: %3.1f < Ratio< %3.1f, Nbin = %d\n"             , fHistoRatioMin,       fHistoRatioMax,       fHistoRatioBins);
  printf("Histograms: %3.1f < Vertex Distance < %3.1f,   Nbin = %d\n", fHistoVertexDistMin,  fHistoVertexDistMax,  fHistoVertexDistBins);
  printf("Histograms: %2.2f < Residual Eta(Z) < %2.2f,   Nbin = %d\n", fHistoTrackResidualEtaMin, fHistoTrackResidualEtaMax,fHistoTrackResidualEtaBins);
  printf("Histograms: %2.2f < Residual Phi(R,X) < %2.2f, Nbin = %d\n", fHistoTrackResidualPhiMin, fHistoTrackResidualPhiMax,fHistoTrackResidualPhiBins);

  printf("    \n") ;
  
}
Commit	Line	Data
745913ae	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
	16	//_________________________________________________________________________
	17	// Class containing histogram settings:
	18	// - Number of bins
	19	// - Min and max of range
	20	//
	21	//-- Author: Gustavo Conesa (LPSC-Grenoble)
	22	//
	23
	24	#include "TMath.h"
	25
	26	#include "AliHistogramRanges.h"
	27
	28	ClassImp(AliHistogramRanges)
	29
	30
	31	//_______________________________________________
	32	AliHistogramRanges::AliHistogramRanges() :
	33	TObject(),
	34	fHistoPtBins(0), fHistoPtMax(0.), fHistoPtMin(0.),
	35	fHistoPhiBins(0), fHistoPhiMax(0.), fHistoPhiMin(0.),
	36	fHistoEtaBins(0), fHistoEtaMax(0.), fHistoEtaMin(0.),
b1494997	37	fHistoDeltaPhiBins(0), fHistoDeltaPhiMax(0.), fHistoDeltaPhiMin(0.),
b1494997	38	fHistoDeltaEtaBins(0), fHistoDeltaEtaMax(0.), fHistoDeltaEtaMin(0.),
745913ae	39	fHistoMassBins(0), fHistoMassMax(0.), fHistoMassMin(0.),
	40	fHistoAsymBins(0), fHistoAsymMax(0.), fHistoAsymMin(0.),
	41	fHistoV0SBins(0), fHistoV0SMax(0), fHistoV0SMin(0),
	42	fHistoV0MBins(0), fHistoV0MMax(0), fHistoV0MMin(0),
	43	fHistoTrMBins(0), fHistoTrMMax(0), fHistoTrMMin(0),
	44	fHistoFinePtBins(1000), fHistoFinePtMax(5.), fHistoFinePtMin(0.),
	45	fHistoPOverEBins(100), fHistoPOverEMax(100.), fHistoPOverEMin(0.),
	46	fHistodEdxBins(100), fHistodEdxMax(100.), fHistodEdxMin(0.),
	47	fHistodRBins(100), fHistodRMax(100.), fHistodRMin(0.),
	48	fHistoTimeBins(100), fHistoTimeMax(100.), fHistoTimeMin(0.),
	49	fHistoNClusCellBins(100), fHistoNClusCellMax(100), fHistoNClusCellMin(0),
	50	fHistoNCellsBins(100), fHistoNCellsMax(100), fHistoNCellsMin(0),
	51	fHistoNClustersBins(100), fHistoNClustersMax(100), fHistoNClustersMin(0),
	52	fHistoRatioBins(100), fHistoRatioMax(100.), fHistoRatioMin(0.),
	53	fHistoVertexDistBins(100), fHistoVertexDistMax(100.), fHistoVertexDistMin(0.),
	54	fHistoRBins(100), fHistoRMax(1000), fHistoRMin(-1000),
	55	fHistoXBins(100), fHistoXMax(1000), fHistoXMin(-1000),
	56	fHistoYBins(100), fHistoYMax(1000), fHistoYMin(-1000),
	57	fHistoZBins(100), fHistoZMax(1000), fHistoZMin(-1000),
	58	fHistoSSBins(0), fHistoSSMax(0), fHistoSSMin(0),
09273901	59	fHistoDiffTimeBins(0), fHistoDiffTimeMax(0), fHistoDiffTimeMin(0),
	60	fHistoTrackResidualEtaBins(0),fHistoTrackResidualEtaMax(0),fHistoTrackResidualEtaMin(0),
	61	fHistoTrackResidualPhiBins(0),fHistoTrackResidualPhiMax(0),fHistoTrackResidualPhiMin(0)
745913ae	62	{
	63	//Default Ctor
	64
	65	//Initialize parameters
	66	InitParameters();
	67
	68	}
	69
	70	//_______________________________________
	71	void AliHistogramRanges::InitParameters()
	72	{
	73	//Initialize histogram parameters
	74
	75	fHistoPtBins = 240 ; fHistoPtMax = 120 ; fHistoPtMin = 0. ;
	76	fHistoPhiBins = 120 ; fHistoPhiMax = TMath::TwoPi(); fHistoPhiMin = 0. ;
69db1204	77	fHistoEtaBins = 100 ; fHistoEtaMax = 1 ; fHistoEtaMin = -1 ;
	78	fHistoDeltaPhiBins = 310 ; fHistoDeltaPhiMax = 4.8 ; fHistoDeltaPhiMin = -1.7 ;
	79	fHistoDeltaEtaBins = 200 ; fHistoDeltaEtaMax = 2 ; fHistoDeltaEtaMin = -2 ;
745913ae	80	fHistoMassBins = 200 ; fHistoMassMax = 1. ; fHistoMassMin = 0. ;
	81	fHistoAsymBins = 10 ; fHistoAsymMax = 1. ; fHistoAsymMin = 0. ;
	82	fHistoV0SBins = 100 ; fHistoV0SMax = 10000 ; fHistoV0SMin = 0 ;
	83	fHistoV0MBins = 100; fHistoV0MMax = 10000 ; fHistoV0MMin = 0 ;
	84	fHistoTrMBins = 100 ; fHistoTrMMax = 2000 ; fHistoTrMMin = 0 ;
	85	fHistoPOverEBins = 100 ; fHistoPOverEMax = 10. ; fHistoPOverEMin = 0. ;
	86	fHistodEdxBins = 200 ; fHistodEdxMax = 400. ; fHistodEdxMin = 0. ;
	87	fHistodRBins = 300 ; fHistodRMax = 3.15 ; fHistodRMin = 0. ;
	88	fHistoTimeBins = 1000; fHistoTimeMax = 1.e3 ; fHistoTimeMin = 0. ;//ns
	89	fHistoNClusCellBins = 200 ; fHistoNClusCellMax = 200 ; fHistoNClusCellMin = 0 ;
	90	fHistoNCellsBins = 300 ; fHistoNCellsMax = 300 ; fHistoNCellsMin = 0 ;
	91	fHistoNClustersBins = 50 ; fHistoNClustersMax = 50 ; fHistoNClustersMin = 0 ;
	92	fHistoRatioBins = 200 ; fHistoRatioMax = 2 ; fHistoRatioMin = 0. ;
	93	fHistoVertexDistBins = 100 ; fHistoVertexDistMax = 500. ; fHistoVertexDistMin = 0. ;
	94	fHistoRBins = 100 ; fHistoRMax = 500 ; fHistoRMin = -500;//cm
	95	fHistoXBins = 100 ; fHistoXMax = 500 ; fHistoXMin = -500;//cm
	96	fHistoYBins = 100 ; fHistoYMax = 500 ; fHistoYMin = -500;//cm
	97	fHistoZBins = 100 ; fHistoZMax = 600 ; fHistoZMin = -500;//cm
09273901	98	fHistoSSBins = 500 ; fHistoSSMax = 5 ; fHistoSSMin = 0 ;
745913ae	99	fHistoDiffTimeBins = 800 ; fHistoDiffTimeMax = 400 ; fHistoDiffTimeMin = -400;// ns
745913ae	100
09273901	101	fHistoTrackResidualEtaBins = 100 ; fHistoTrackResidualEtaMax = 0.15 ; fHistoTrackResidualEtaMin = -0.15;
	102	fHistoTrackResidualPhiBins = 100 ; fHistoTrackResidualPhiMax = 0.15 ; fHistoTrackResidualPhiMin = -0.15;
	103
745913ae	104	}
	105
	106	//________________________________________________________
	107	void AliHistogramRanges::Print(const Option_t * /opt/) const
	108	{
	109	//Print some relevant parameters set for the analysis
	110
09273901	111	printf("Histograms: %3.1f < pT < %3.1f, Nbin = %d\n" , fHistoPtMin, fHistoPtMax, fHistoPtBins);
	112	printf("Histograms: %3.1f < phi < %3.1f, Nbin = %d\n" , fHistoPhiMin, fHistoPhiMax, fHistoPhiBins);
	113	printf("Histograms: %3.1f < eta < %3.1f, Nbin = %d\n" , fHistoEtaMin, fHistoEtaMax, fHistoEtaBins);
b1494997	114	printf("Histograms: %3.1f < delta phi < %3.1f, Nbin = %d\n" , fHistoDeltaPhiMin, fHistoDeltaPhiMax, fHistoDeltaPhiBins);
b1494997	115	printf("Histograms: %3.1f < delta eta < %3.1f, Nbin = %d\n" , fHistoDeltaEtaMin, fHistoDeltaEtaMax, fHistoDeltaEtaBins);
09273901	116	printf("Histograms: %3.1f < mass < %3.1f, Nbin = %d\n" , fHistoMassMin, fHistoMassMax, fHistoMassBins);
	117	printf("Histograms: %3.1f < asymmetry < %3.1f, Nbin = %d\n" , fHistoAsymMin, fHistoAsymMax, fHistoAsymBins);
	118	printf("Histograms: %d < V0 Signal < %d, Nbin = %d\n" , fHistoV0SMin, fHistoV0SMax, fHistoV0SBins);
	119	printf("Histograms: %d < V0 Mult < %d, Nbin = %d\n" , fHistoV0MMin, fHistoV0MMax, fHistoV0MBins);
	120	printf("Histograms: %d < Track Mult < %d, Nbin = %d\n" , fHistoTrMMin, fHistoTrMMax, fHistoTrMBins);
	121	printf("Histograms: %3.1f < p/E < %3.1f, Nbin = %d\n" , fHistoPOverEMin, fHistoPOverEMax, fHistoPOverEBins);
	122	printf("Histograms: %3.1f < dEdx < %3.1f, Nbin = %d\n" , fHistodEdxMin, fHistodEdxMax, fHistodEdxBins);
	123	printf("Histograms: %3.1f < dR (track cluster)< %3.1f, Nbin = %d\n", fHistodRMin, fHistodRMax, fHistodRBins);
	124	printf("Histograms: %3.1f < R=sqrt{x^2+y^2} < %3.1f, Nbin = %d\n", fHistoRMin, fHistoRMax, fHistoRBins);
	125	printf("Histograms: %3.1f < X < %3.1f, Nbin = %d\n" , fHistoXMin, fHistoXMax, fHistoXBins);
	126	printf("Histograms: %3.1f < Y < %3.1f, Nbin = %d\n" , fHistoYMin, fHistoYMax, fHistoYBins);
	127	printf("Histograms: %3.1f < Z < %3.1f, Nbin = %d\n" , fHistoZMin, fHistoZMax, fHistoZBins);
	128	printf("Histograms: %g < Time < %g, Nbin = %d\n" , fHistoTimeMin, fHistoTimeMax, fHistoTimeBins);
	129	printf("Histograms: %d < N cells per cluster < %d, Nbin = %d\n" , fHistoNClusCellMin, fHistoNClusCellMax, fHistoNClusCellBins);
	130	printf("Histograms: %d < N cells < %d, Nbin = %d\n" , fHistoNCellsMin, fHistoNCellsMax, fHistoNCellsBins);
	131	printf("Histograms: %d < N clusters < %d, Nbin = %d\n" , fHistoNClustersMin, fHistoNClustersMax, fHistoNClustersBins);
	132	printf("Histograms: %3.1f < Ratio< %3.1f, Nbin = %d\n" , fHistoRatioMin, fHistoRatioMax, fHistoRatioBins);
	133	printf("Histograms: %3.1f < Vertex Distance < %3.1f, Nbin = %d\n", fHistoVertexDistMin, fHistoVertexDistMax, fHistoVertexDistBins);
	134	printf("Histograms: %2.2f < Residual Eta(Z) < %2.2f, Nbin = %d\n", fHistoTrackResidualEtaMin, fHistoTrackResidualEtaMax,fHistoTrackResidualEtaBins);
	135	printf("Histograms: %2.2f < Residual Phi(R,X) < %2.2f, Nbin = %d\n", fHistoTrackResidualPhiMin, fHistoTrackResidualPhiMax,fHistoTrackResidualPhiBins);
	136
745913ae	137	printf(" \n") ;
	138
	139	}
	140
	141
	142