[u/mrichter/AliRoot.git] / ZDC / AliZDCReco.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.                  *
 **************************************************************************/

/* $Id$ */

////////////////////////////////////////////////
//  RecPoints classes for set ZDC             //
//  This class reconstructs the space         //
//  points from digits                        //
//  for the ZDC calorimeter                   //
////////////////////////////////////////////////


#include "AliZDCReco.h"

ClassImp(AliZDCReco)
  

//_____________________________________________________________________________
AliZDCReco::AliZDCReco() :
	
  TObject(),
  //
  fNDetSpecNSideA(0),
  fNDetSpecPSideA(0),
  fNDetSpecNSideC(0),
  fNDetSpecPSideC(0),
  fNTrueSpectators(0),
  fNTrueSpecSideA(0),
  fNTrueSpecSideC(0),
  fNParticipants(0),
  fNPartSideA(0),
  fNPartSideC(0),
  fImpParameter(0),
  fImpParSideA(0),
  fImpParSideC(0),
  fRecoFlag(0x0),
  fIsScalerOn(kFALSE)
{ 
  //
  // Default constructor
  //
  for(Int_t i=0; i<10; i++){
     fZN1EnTow[i] = fZP1EnTow[i] = fZN2EnTow[i] = fZP2EnTow[i] = 0.;
     if(i<2){
       fZN1Energy[i] = fZP1Energy[i] = fZN2Energy[i] = fZP2Energy[i] =  0.;
       fZEM1signal[i] = fZEM2signal[i] = 0.;
       fPMRef1[i] = fPMRef2[i] = 0.;
     }
  }
  
  for(Int_t i=0; i<32; i++){
    fZDCScaler[i] = 0;
    for(Int_t ij=0; ij<4; ij++) fZDCTDCData[i][ij] = 0;
  }
}
  

//_____________________________________________________________________________
AliZDCReco::AliZDCReco(
     Float_t* ezn1, Float_t* ezp1, Float_t* ezn2, Float_t* ezp2,  
     Float_t* ezn1tow, Float_t* ezp1tow, Float_t* ezn2tow, Float_t* ezp2tow, 
     Float_t* ezem1, Float_t* ezem2, Float_t* ref1, Float_t* ref2, 
     //    
     Int_t detspnSideA,  Int_t detsppSideA, Int_t detspnSideC, Int_t detsppSideC,  
     Int_t trsp, Int_t trspSideA,Int_t trspSideC,
     Int_t npart, Int_t npartSideA, Int_t npartSideC, 
     Float_t b, Float_t bSideA, Float_t bSideC,
     UInt_t recoFlag, Bool_t scalerOn, UInt_t* scaler, Int_t tdcData[32][4]) :
	
  TObject(),
  //
  fNDetSpecNSideA(detspnSideA),
  fNDetSpecPSideA(detsppSideA),
  fNDetSpecNSideC(detspnSideC),
  fNDetSpecPSideC(detsppSideC),
  fNTrueSpectators(trsp),
  fNTrueSpecSideA(trspSideA),
  fNTrueSpecSideC(trspSideC),
  fNParticipants(npart),
  fNPartSideA(npartSideA),
  fNPartSideC(npartSideC),
  fImpParameter(b),
  fImpParSideA(bSideA),
  fImpParSideC(bSideC),
  fRecoFlag(recoFlag),
  fIsScalerOn(scalerOn)
{ 
  //
  // Constructor
  //
  for(Int_t j=0; j<10; j++){
     fZN1EnTow[j] =  ezn1tow[j];
     fZP1EnTow[j] =  ezp1tow[j];
     fZN2EnTow[j] =  ezn2tow[j];
     fZP2EnTow[j] =  ezp2tow[j];
     if(j<2){
       fZN1Energy[j] = ezn1[j];
       fZP1Energy[j] = ezp1[j];
       fZN2Energy[j] = ezn2[j];
       fZP2Energy[j] = ezp2[j];
       fZEM1signal[j] = ezem1[j];
       fZEM2signal[j] = ezem2[j];
       fPMRef1[j] = ref1[j];
       fPMRef2[j] = ref2[j];
     }
  }
  for(Int_t j=0; j<32; j++){
    fZDCScaler[j] = scaler[j];
    for(Int_t y=0; y<4; y++) fZDCTDCData[j][y] = tdcData[j][y];
  }
}

//______________________________________________________________________________
AliZDCReco::AliZDCReco(const AliZDCReco &oldreco) :
TObject(),
fNDetSpecNSideA(oldreco.GetNDetSpecNSideA()),
fNDetSpecPSideA(oldreco.GetNDetSpecPSideA()),
fNDetSpecNSideC(oldreco.GetNDetSpecNSideC()),        
fNDetSpecPSideC(oldreco.GetNDetSpecPSideC()),       
fNTrueSpectators(oldreco.GetNTrueSpectators()),
fNTrueSpecSideA(oldreco.GetNTrueSpecSideA()),
fNTrueSpecSideC(oldreco.GetNTrueSpecSideC()),  	
fNParticipants(oldreco.GetNParticipants()),		       
fNPartSideA(oldreco.GetNPartSideA()),		       
fNPartSideC(oldreco.GetNPartSideC()),  		       
fImpParameter(oldreco.GetImpParameter()),      
fImpParSideA(oldreco.GetImpParSideA()),      
fImpParSideC(oldreco.GetImpParSideC()),
fRecoFlag(oldreco.GetRecoFlag()),
fIsScalerOn(oldreco.IsScalerOn())    
{
  // Copy constructor

  fZN1Energy[0]  = oldreco.GetZN1HREnergy();
  fZP1Energy[0]  = oldreco.GetZP1HREnergy();		
  fZN2Energy[0]  = oldreco.GetZN2HREnergy();	     
  fZP2Energy[0]  = oldreco.GetZP2HREnergy();	 
  //    
  fZN1Energy[1]  = oldreco.GetZN1LREnergy();
  fZP1Energy[1]  = oldreco.GetZP1LREnergy();	       
  fZN2Energy[1]  = oldreco.GetZN2LREnergy();	    
  fZP2Energy[1]  = oldreco.GetZP2LREnergy();	    
  //
  for(Int_t i=0; i<5; i++){	  
     fZN1EnTow[i]  = oldreco.GetZN1HREnTow(i);
     fZP1EnTow[i]  = oldreco.GetZP1HREnTow(i);
     fZN2EnTow[i]  = oldreco.GetZN2HREnTow(i);
     fZP2EnTow[i]  = oldreco.GetZP2HREnTow(i);
     fZN1EnTow[i+5]  = oldreco.GetZN1LREnTow(i);
     fZP1EnTow[i+5]  = oldreco.GetZP1LREnTow(i);
     fZN2EnTow[i+5]  = oldreco.GetZN2LREnTow(i);
     fZP2EnTow[i+5]  = oldreco.GetZP2LREnTow(i);
  }
  fZEM1signal[0] = oldreco.GetZEM1HRsignal();
  fZEM1signal[1] = oldreco.GetZEM1LRsignal();
  fZEM2signal[0] = oldreco.GetZEM2HRsignal();
  fZEM2signal[1] = oldreco.GetZEM2LRsignal();
  fPMRef1[0] = oldreco.GetPMRef1HRsignal();
  fPMRef1[1] = oldreco.GetPMRef1LRsignal();
  fPMRef2[0] = oldreco.GetPMRef2HRsignal();
  fPMRef2[1] = oldreco.GetPMRef2LRsignal();
  for(Int_t j=0; j<32; j++){
    fZDCScaler[j] = oldreco.GetZDCScaler(j);
    for(Int_t y=0; y<4; y++) fZDCTDCData[j][y] = oldreco.GetZDCTDCData(j, y);
  }
}

//______________________________________________________________________________
void AliZDCReco::Print(Option_t *) const {
  //
  // Printing Reconstruction Parameters
  //
  printf(" ****************** AliZDCReco object ******************\n"
  	 "       ---------------   side A ---------------\n"
	 " E_ZN = %1.2f TeV, E_ZP = %1.2f TeV, "
	 " E_ZEM1 =  %1.2f TeV,  E_ZEM2 = %1.2f TeV\n "
	 " N_spec_n = %d, N_spec_p = %d,"
	 " N_part = %d, b = %1.4f fm\n"
  	 "       ---------------   side C ---------------\n"
	 " E_ZN = %1.2f TeV, E_ZP = %1.2f TeV, "
	 " N_spec_n = %d, N_spec_p = %d,"
	 " N_part = %d, b = %1.4f fm\n"
         " *******************************************************\n",
	 fZN2Energy[0]/1000., fZP2Energy[0]/1000.,
	 fZEM1signal[0]/1000.,fZEM2signal[0]/1000.,
	 fNDetSpecNSideA,fNDetSpecPSideA, fNPartSideA,fImpParSideA,
	 fZN1Energy[0]/1000.,fZP1Energy[0]/1000.,	
         fNDetSpecNSideC,fNDetSpecPSideC,fNPartSideC,fImpParSideC);
	 
}
Commit	Line	Data
1450a7cd	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
803d1ab0	16	/* $Id$ */
1450a7cd	17
	18	////////////////////////////////////////////////
	19	// RecPoints classes for set ZDC //
93f80a6f	20	// This class reconstructs the space //
	21	// points from digits //
	22	// for the ZDC calorimeter //
1450a7cd	23	////////////////////////////////////////////////
	24
	25
	26	#include "AliZDCReco.h"
	27
	28	ClassImp(AliZDCReco)
	29
cc2abffd	30
	31	//_____________________________________________________________________________
	32	AliZDCReco::AliZDCReco() :
	33
	34	TObject(),
646f1679	35	//
73bc3a3f	36	fNDetSpecNSideA(0),
	37	fNDetSpecPSideA(0),
	38	fNDetSpecNSideC(0),
	39	fNDetSpecPSideC(0),
	40	fNTrueSpectators(0),
	41	fNTrueSpecSideA(0),
	42	fNTrueSpecSideC(0),
	43	fNParticipants(0),
	44	fNPartSideA(0),
	45	fNPartSideC(0),
	46	fImpParameter(0),
	47	fImpParSideA(0),
1e319f71	48	fImpParSideC(0),
81f09162	49	fRecoFlag(0x0),
81f09162	50	fIsScalerOn(kFALSE)
cc2abffd	51	{
	52	//
	53	// Default constructor
	54	//
c35ed519	55	for(Int_t i=0; i<10; i++){
	56	fZN1EnTow[i] = fZP1EnTow[i] = fZN2EnTow[i] = fZP2EnTow[i] = 0.;
	57	if(i<2){
	58	fZN1Energy[i] = fZP1Energy[i] = fZN2Energy[i] = fZP2Energy[i] = 0.;
	59	fZEM1signal[i] = fZEM2signal[i] = 0.;
	60	fPMRef1[i] = fPMRef2[i] = 0.;
	61	}
84d6255e	62	}
81f09162	63
32e2fda5	64	for(Int_t i=0; i<32; i++){
32e2fda5	65	fZDCScaler[i] = 0;
f53e5ecb	66	for(Int_t ij=0; ij<4; ij++) fZDCTDCData[i][ij] = 0;
32e2fda5	67	}
cc2abffd	68	}
	69
	70
1450a7cd	71	//_____________________________________________________________________________
73bc3a3f	72	AliZDCReco::AliZDCReco(
	73	Float_t* ezn1, Float_t* ezp1, Float_t* ezn2, Float_t* ezp2,
	74	Float_t* ezn1tow, Float_t* ezp1tow, Float_t* ezn2tow, Float_t* ezp2tow,
	75	Float_t* ezem1, Float_t* ezem2, Float_t* ref1, Float_t* ref2,
	76	//
	77	Int_t detspnSideA, Int_t detsppSideA, Int_t detspnSideC, Int_t detsppSideC,
	78	Int_t trsp, Int_t trspSideA,Int_t trspSideC,
	79	Int_t npart, Int_t npartSideA, Int_t npartSideC,
1e319f71	80	Float_t b, Float_t bSideA, Float_t bSideC,
f53e5ecb	81	UInt_t recoFlag, Bool_t scalerOn, UInt_t* scaler, Int_t tdcData[32][4]) :
cc2abffd	82
cc2abffd	83	TObject(),
646f1679	84	//
73bc3a3f	85	fNDetSpecNSideA(detspnSideA),
	86	fNDetSpecPSideA(detsppSideA),
	87	fNDetSpecNSideC(detspnSideC),
	88	fNDetSpecPSideC(detsppSideC),
	89	fNTrueSpectators(trsp),
	90	fNTrueSpecSideA(trspSideA),
	91	fNTrueSpecSideC(trspSideC),
	92	fNParticipants(npart),
	93	fNPartSideA(npartSideA),
	94	fNPartSideC(npartSideC),
	95	fImpParameter(b),
	96	fImpParSideA(bSideA),
1e319f71	97	fImpParSideC(bSideC),
81f09162	98	fRecoFlag(recoFlag),
81f09162	99	fIsScalerOn(scalerOn)
1450a7cd	100	{
93f80a6f	101	//
cc2abffd	102	// Constructor
93f80a6f	103	//
c35ed519	104	for(Int_t j=0; j<10; j++){
646f1679	105	fZN1EnTow[j] = ezn1tow[j];
	106	fZP1EnTow[j] = ezp1tow[j];
	107	fZN2EnTow[j] = ezn2tow[j];
	108	fZP2EnTow[j] = ezp2tow[j];
c35ed519	109	if(j<2){
	110	fZN1Energy[j] = ezn1[j];
	111	fZP1Energy[j] = ezp1[j];
	112	fZN2Energy[j] = ezn2[j];
	113	fZP2Energy[j] = ezp2[j];
	114	fZEM1signal[j] = ezem1[j];
	115	fZEM2signal[j] = ezem2[j];
	116	fPMRef1[j] = ref1[j];
	117	fPMRef2[j] = ref2[j];
	118	}
646f1679	119	}
82dffa48	120	for(Int_t j=0; j<32; j++){
82dffa48	121	fZDCScaler[j] = scaler[j];
32e2fda5	122	for(Int_t y=0; y<4; y++) fZDCTDCData[j][y] = tdcData[j][y];
82dffa48	123	}
1450a7cd	124	}
93f80a6f	125
cc2abffd	126	//______________________________________________________________________________
cc2abffd	127	AliZDCReco::AliZDCReco(const AliZDCReco &oldreco) :
a718c993	128	TObject(),
73bc3a3f	129	fNDetSpecNSideA(oldreco.GetNDetSpecNSideA()),
	130	fNDetSpecPSideA(oldreco.GetNDetSpecPSideA()),
	131	fNDetSpecNSideC(oldreco.GetNDetSpecNSideC()),
	132	fNDetSpecPSideC(oldreco.GetNDetSpecPSideC()),
	133	fNTrueSpectators(oldreco.GetNTrueSpectators()),
	134	fNTrueSpecSideA(oldreco.GetNTrueSpecSideA()),
	135	fNTrueSpecSideC(oldreco.GetNTrueSpecSideC()),
	136	fNParticipants(oldreco.GetNParticipants()),
	137	fNPartSideA(oldreco.GetNPartSideA()),
	138	fNPartSideC(oldreco.GetNPartSideC()),
	139	fImpParameter(oldreco.GetImpParameter()),
	140	fImpParSideA(oldreco.GetImpParSideA()),
1e319f71	141	fImpParSideC(oldreco.GetImpParSideC()),
81f09162	142	fRecoFlag(oldreco.GetRecoFlag()),
81f09162	143	fIsScalerOn(oldreco.IsScalerOn())
cc2abffd	144	{
	145	// Copy constructor
	146
c35ed519	147	fZN1Energy[0] = oldreco.GetZN1HREnergy();
	148	fZP1Energy[0] = oldreco.GetZP1HREnergy();
	149	fZN2Energy[0] = oldreco.GetZN2HREnergy();
	150	fZP2Energy[0] = oldreco.GetZP2HREnergy();
	151	//
	152	fZN1Energy[1] = oldreco.GetZN1LREnergy();
	153	fZP1Energy[1] = oldreco.GetZP1LREnergy();
	154	fZN2Energy[1] = oldreco.GetZN2LREnergy();
	155	fZP2Energy[1] = oldreco.GetZP2LREnergy();
646f1679	156	//
646f1679	157	for(Int_t i=0; i<5; i++){
c35ed519	158	fZN1EnTow[i] = oldreco.GetZN1HREnTow(i);
	159	fZP1EnTow[i] = oldreco.GetZP1HREnTow(i);
	160	fZN2EnTow[i] = oldreco.GetZN2HREnTow(i);
	161	fZP2EnTow[i] = oldreco.GetZP2HREnTow(i);
	162	fZN1EnTow[i+5] = oldreco.GetZN1LREnTow(i);
	163	fZP1EnTow[i+5] = oldreco.GetZP1LREnTow(i);
	164	fZN2EnTow[i+5] = oldreco.GetZN2LREnTow(i);
	165	fZP2EnTow[i+5] = oldreco.GetZP2LREnTow(i);
646f1679	166	}
c35ed519	167	fZEM1signal[0] = oldreco.GetZEM1HRsignal();
	168	fZEM1signal[1] = oldreco.GetZEM1LRsignal();
	169	fZEM2signal[0] = oldreco.GetZEM2HRsignal();
	170	fZEM2signal[1] = oldreco.GetZEM2LRsignal();
	171	fPMRef1[0] = oldreco.GetPMRef1HRsignal();
	172	fPMRef1[1] = oldreco.GetPMRef1LRsignal();
	173	fPMRef2[0] = oldreco.GetPMRef2HRsignal();
	174	fPMRef2[1] = oldreco.GetPMRef2LRsignal();
82dffa48	175	for(Int_t j=0; j<32; j++){
82dffa48	176	fZDCScaler[j] = oldreco.GetZDCScaler(j);
32e2fda5	177	for(Int_t y=0; y<4; y++) fZDCTDCData[j][y] = oldreco.GetZDCTDCData(j, y);
82dffa48	178	}
cc2abffd	179	}
cc2abffd	180
93f80a6f	181	//______________________________________________________________________________
	182	void AliZDCReco::Print(Option_t *) const {
	183	//
	184	// Printing Reconstruction Parameters
	185	//
73bc3a3f	186	printf(" **************** AliZDCReco object ****************\n"
73bc3a3f	187	" --------------- side A ---------------\n"
c2bb665a	188	" E_ZN = %1.2f TeV, E_ZP = %1.2f TeV, "
	189	" E_ZEM1 = %1.2f TeV, E_ZEM2 = %1.2f TeV\n "
	190	" N_spec_n = %d, N_spec_p = %d,"
	191	" N_part = %d, b = %1.4f fm\n"
73bc3a3f	192	" --------------- side C ---------------\n"
c2bb665a	193	" E_ZN = %1.2f TeV, E_ZP = %1.2f TeV, "
	194	" N_spec_n = %d, N_spec_p = %d,"
	195	" N_part = %d, b = %1.4f fm\n"
73bc3a3f	196	" *******************************************************\n",
4e1fc320	197	fZN2Energy[0]/1000., fZP2Energy[0]/1000.,
73bc3a3f	198	fZEM1signal[0]/1000.,fZEM2signal[0]/1000.,
	199	fNDetSpecNSideA,fNDetSpecPSideA, fNPartSideA,fImpParSideA,
	200	fZN1Energy[0]/1000.,fZP1Energy[0]/1000.,
	201	fNDetSpecNSideC,fNDetSpecPSideC,fNPartSideC,fImpParSideC);
	202
93f80a6f	203	}