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


//-------------------------------------------------------------------------
//                      Implementation of   Class AliESDZDC
//   This is a class that summarizes the ZDC data
//   for the ESD   
//   Origin: Christian Klein-Boesing, CERN, Christian.Klein-Boesing@cern.ch 
//-------------------------------------------------------------------------

#include <TMath.h>

#include "AliESDZDC.h"

ClassImp(AliESDZDC)

//______________________________________________________________________________
AliESDZDC::AliESDZDC() :
  TObject(),
  fZDCN1Energy(0),
  fZDCP1Energy(0),
  fZDCN2Energy(0),
  fZDCP2Energy(0),
  fZDCEMEnergy(0),
  fZDCEMEnergy1(0),
  fZDCParticipants(0),
  fZDCPartSideA(0),
  fZDCPartSideC(0),
  fImpactParameter(0),
  fImpactParamSideA(0),
  fImpactParamSideC(0),
  fESDQuality(0)
{
  for(Int_t i=0; i<5; i++){
    fZN1TowerEnergy[i] = fZN2TowerEnergy[i] = 0.;
    fZP1TowerEnergy[i] = fZP2TowerEnergy[i] = 0.;
    fZN1TowerEnergyLR[i] = fZN2TowerEnergyLR[i] = 0.;
    fZP1TowerEnergyLR[i] = fZP2TowerEnergyLR[i] = 0.;
  }
  for(Int_t i=0; i<2; i++){
    fZNACentrCoord[i] = fZNCCentrCoord[i] = 0.;
  }
  for(Int_t i=0; i<32; i++){
    fVMEScaler[i]=0;
    for(Int_t y=0; y<4; y++) fZDCTDCData[i][y]=0;
  }
}

//______________________________________________________________________________
AliESDZDC::AliESDZDC(const AliESDZDC& zdc) :
  TObject(zdc),
  fZDCN1Energy(zdc.fZDCN1Energy),
  fZDCP1Energy(zdc.fZDCP1Energy),
  fZDCN2Energy(zdc.fZDCN2Energy),
  fZDCP2Energy(zdc.fZDCP2Energy),
  fZDCEMEnergy(zdc.fZDCEMEnergy),
  fZDCEMEnergy1(zdc.fZDCEMEnergy1),
  fZDCParticipants(zdc.fZDCParticipants),
  fZDCPartSideA(zdc.fZDCPartSideA),
  fZDCPartSideC(zdc.fZDCPartSideC),
  fImpactParameter(zdc.fImpactParameter),
  fImpactParamSideA(zdc.fImpactParamSideA),
  fImpactParamSideC(zdc.fImpactParamSideC),
  fESDQuality(zdc.fESDQuality)
{
  // copy constructor
  for(Int_t i=0; i<5; i++){
     fZN1TowerEnergy[i] = zdc.fZN1TowerEnergy[i];
     fZN2TowerEnergy[i] = zdc.fZN2TowerEnergy[i];
     fZP1TowerEnergy[i] = zdc.fZP1TowerEnergy[i];
     fZP2TowerEnergy[i] = zdc.fZP2TowerEnergy[i];
     fZN1TowerEnergyLR[i] = zdc.fZN1TowerEnergyLR[i];
     fZN2TowerEnergyLR[i] = zdc.fZN2TowerEnergyLR[i];
     fZP1TowerEnergyLR[i] = zdc.fZP1TowerEnergyLR[i];
     fZP2TowerEnergyLR[i] = zdc.fZP2TowerEnergyLR[i];
  }
  for(Int_t i=0; i<2; i++){
    fZNACentrCoord[i] = zdc.fZNACentrCoord[i];
    fZNCCentrCoord[i] = zdc.fZNCCentrCoord[i];
  }
  for(Int_t i=0; i<32; i++){
    fVMEScaler[i] = zdc.fVMEScaler[i];
    for(Int_t y=0; y<4; y++) fZDCTDCData[i][y] = zdc.fZDCTDCData[i][y];
  }
}

//______________________________________________________________________________
AliESDZDC& AliESDZDC::operator=(const AliESDZDC&zdc)
{
  // assigment operator
  if(this!=&zdc) {
    TObject::operator=(zdc);
    fZDCN1Energy = zdc.fZDCN1Energy;
    fZDCP1Energy = zdc.fZDCP1Energy;
    fZDCN2Energy = zdc.fZDCN2Energy;
    fZDCP2Energy = zdc.fZDCP2Energy;
    fZDCEMEnergy = zdc.fZDCEMEnergy;
    fZDCEMEnergy1 = zdc.fZDCEMEnergy1;
    for(Int_t i=0; i<5; i++){
       fZN1TowerEnergy[i] = zdc.fZN1TowerEnergy[i];
       fZN2TowerEnergy[i] = zdc.fZN2TowerEnergy[i];
       fZP1TowerEnergy[i] = zdc.fZP1TowerEnergy[i];
       fZP2TowerEnergy[i] = zdc.fZP2TowerEnergy[i];
       fZN1TowerEnergyLR[i] = zdc.fZN1TowerEnergyLR[i];
       fZN2TowerEnergyLR[i] = zdc.fZN2TowerEnergyLR[i];
       fZP1TowerEnergyLR[i] = zdc.fZP1TowerEnergyLR[i];
       fZP2TowerEnergyLR[i] = zdc.fZP2TowerEnergyLR[i];
    }
    //
    fZDCParticipants = zdc.fZDCParticipants;
    fZDCPartSideA = zdc.fZDCPartSideA;
    fZDCPartSideC = zdc.fZDCPartSideC;
    fImpactParameter = zdc.fImpactParameter;
    fImpactParamSideA = zdc.fImpactParamSideA;
    fImpactParamSideC = zdc.fImpactParamSideC;
    //
    for(Int_t i=0; i<2; i++){
         fZNACentrCoord[i] = zdc.fZNACentrCoord[i];
         fZNCCentrCoord[i] = zdc.fZNCCentrCoord[i];
    }
    //
    fESDQuality = zdc.fESDQuality;
    for(Int_t i=0; i<32; i++){
      fVMEScaler[i] = zdc.fVMEScaler[i];
      for(Int_t y=0; y<4; y++) fZDCTDCData[i][y] = zdc.fZDCTDCData[i][y];
    }
  } 
  return *this;
}

//______________________________________________________________________________
void AliESDZDC::Copy(TObject &obj) const {
  
  // this overwrites the virtual TOBject::Copy()
  // to allow run time copying without casting
  // in AliESDEvent

  if(this==&obj)return;
  AliESDZDC *robj = dynamic_cast<AliESDZDC*>(&obj);
  if(!robj)return; // not an AliESDZDC
  *robj = *this;

}


//______________________________________________________________________________
void AliESDZDC::Reset()
{
  // reset all data members
  fZDCN1Energy=0;
  fZDCP1Energy=0;
  fZDCN2Energy=0;
  fZDCP2Energy=0;
  fZDCEMEnergy=0;
  fZDCEMEnergy1=0;
  for(Int_t i=0; i<5; i++){
    fZN1TowerEnergy[i] = fZN2TowerEnergy[i] = 0.;
    fZP1TowerEnergy[i] = fZP2TowerEnergy[i] = 0.;
    fZN1TowerEnergyLR[i] = fZN2TowerEnergyLR[i] = 0.;
    fZP1TowerEnergyLR[i] = fZP2TowerEnergyLR[i] = 0.;
  }
  fZDCParticipants=0;  
  fZDCPartSideA=0;  
  fZDCPartSideC=0;  
  fImpactParameter=0;
  fImpactParamSideA=0;
  fImpactParamSideC=0;
  for(Int_t i=0; i<2; i++){
       fZNACentrCoord[i] = fZNCCentrCoord[i] = 0.;
  }
  fESDQuality=0;
  for(Int_t i=0; i<32; i++){
     fVMEScaler[i] = 0;
     for(Int_t y=0; y<4; y++) fZDCTDCData[i][y] = 0;
  }
}

//______________________________________________________________________________
void AliESDZDC::Print(const Option_t *) const
{
  //  Print ESD for the ZDC
  printf("\n \t E_{ZNC} = %f TeV, E_{ZPC} = %f TeV, E_{ZNA} = %f TeV, E_{ZPA} = %f TeV,"
  " E_{ZEM} = %f GeV, Npart = %d, b = %1.2f fm\n",
  fZDCN1Energy/1000.,fZDCP1Energy/1000.,fZDCN2Energy/1000.,fZDCP2Energy/1000.,
  fZDCEMEnergy+fZDCEMEnergy1, fZDCParticipants,fImpactParameter);
  //
  printf(" ### VMEScaler (!=0): \n");
  for(Int_t i=0; i<32; i++) if(fVMEScaler[i]!=0) printf("\t %d \n",fVMEScaler[i]);
  printf("\n");
  //
  printf(" ### TDCData (!=0): \n");
  for(Int_t i=0; i<32; i++){
    for(Int_t j=0; j<4; j++)
      if(fZDCTDCData[i][j]!=0) printf("\t %d \n",fZDCTDCData[i][j]);
  }
  printf("\n");
}

//______________________________________________________________________________
Bool_t AliESDZDC::GetZNCentroidInPbPb(Float_t beamEne, Double_t centrZNC[2], Double_t centrZNA[2]) 
{
  // Provide coordinates of centroid over ZN (side C) front face
  if(beamEne==0){
    printf(" ZDC centroid in PbPb can't be calculated with E_beam = 0 !!!\n");
    for(Int_t jj=0; jj<2; jj++) fZNCCentrCoord[jj] = 999.;
    return kFALSE;
  }

  const Float_t x[4] = {-1.75, 1.75, -1.75, 1.75};
  const Float_t y[4] = {-1.75, -1.75, 1.75, 1.75};
  const Float_t alpha=0.395;
  Float_t numXZNC=0., numYZNC=0., denZNC=0., cZNC, wZNC; 
  Float_t numXZNA=0., numYZNA=0., denZNA=0., cZNA, wZNA; 
  //
  for(Int_t i=0; i<4; i++){
    if(fZN1TowerEnergy[i+1]>0.) {
      wZNC = TMath::Power(fZN1TowerEnergy[i+1], alpha);
      numXZNC += x[i]*wZNC;
      numYZNC += y[i]*wZNC;
      denZNC += wZNC;
    }
    if(fZN2TowerEnergy[i+1]>0.) {
      wZNA = TMath::Power(fZN2TowerEnergy[i+1], alpha);
      numXZNA += x[i]*wZNA;
      numYZNA += y[i]*wZNA;
      denZNA += wZNA;
    }
  }
  //
  if(denZNC!=0){
    Float_t nSpecnC = fZDCN1Energy/beamEne;
    cZNC = 1.89358-0.71262/(nSpecnC+0.71789);
    fZNCCentrCoord[0] = cZNC*numXZNC/denZNC;
    fZNCCentrCoord[1] = cZNC*numYZNC/denZNC;
  } 
  else{
    fZNCCentrCoord[0] = fZNCCentrCoord[1] = 999.;
  }
  if(denZNA!=0){
    Float_t nSpecnA = fZDCN1Energy/beamEne;
    cZNA = 1.89358-0.71262/(nSpecnA+0.71789);
    fZNACentrCoord[0] = cZNA*numXZNA/denZNA;
    fZNACentrCoord[1] = cZNA*numYZNA/denZNA;
  } 
  else{
    fZNACentrCoord[0] = fZNACentrCoord[1] = 999.;
  }
  //
  for(Int_t il=0; il<2; il++){
    centrZNC[il] = fZNCCentrCoord[il];
    centrZNA[il] = fZNACentrCoord[il];
  }
  
  return kTRUE;
}

//______________________________________________________________________________
Bool_t AliESDZDC::GetZNCentroidInpp(Double_t centrZNC[2], Double_t centrZNA[2]) 
{
  // Provide coordinates of centroid over ZN (side C) front face
  const Float_t x[4] = {-1.75, 1.75, -1.75, 1.75};
  const Float_t y[4] = {-1.75, -1.75, 1.75, 1.75};
  const Float_t alpha=0.5;
  Float_t numXZNC=0., numYZNC=0., denZNC=0., wZNC; 
  Float_t numXZNA=0., numYZNA=0., denZNA=0., wZNA; 
  //
  for(Int_t i=0; i<4; i++){
    if(fZN1TowerEnergy[i+1]>0.) {
      wZNC = TMath::Power(fZN1TowerEnergy[i+1], alpha);
      numXZNC += x[i]*wZNC;
      numYZNC += y[i]*wZNC;
      denZNC += wZNC;
    }
    if(fZN2TowerEnergy[i+1]>0.) {
      wZNA = TMath::Power(fZN2TowerEnergy[i+1], alpha);
      numXZNA += x[i]*wZNA;
      numYZNA += y[i]*wZNA;
      denZNA += wZNA;
    }
  }
  //
  if(denZNC!=0){
    fZNCCentrCoord[0] = numXZNC/denZNC;
    fZNCCentrCoord[1] = numYZNC/denZNC;
  } 
  else{
    fZNCCentrCoord[0] = fZNCCentrCoord[1] = 999.;
  }
  if(denZNA!=0){
    fZNACentrCoord[0] = numXZNA/denZNA;
    fZNACentrCoord[1] = numYZNA/denZNA;
  } 
  else{
    fZNACentrCoord[0] = fZNACentrCoord[1] = 999.;
  }
  //
  for(Int_t il=0; il<2; il++){
    centrZNC[il] = fZNCCentrCoord[il];
    centrZNA[il] = fZNACentrCoord[il];
  }
  
  return kTRUE;
}
Commit	Line	Data
d5ebf00e	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	//-------------------------------------------------------------------------
	18	// Implementation of Class AliESDZDC
	19	// This is a class that summarizes the ZDC data
	20	// for the ESD
	21	// Origin: Christian Klein-Boesing, CERN, Christian.Klein-Boesing@cern.ch
	22	//-------------------------------------------------------------------------
	23
cb6a0388	24	#include <TMath.h>
d5ebf00e	25
	26	#include "AliESDZDC.h"
	27
	28	ClassImp(AliESDZDC)
	29
	30	//______________________________________________________________________________
	31	AliESDZDC::AliESDZDC() :
	32	TObject(),
	33	fZDCN1Energy(0),
	34	fZDCP1Energy(0),
	35	fZDCN2Energy(0),
	36	fZDCP2Energy(0),
f6536de8	37	fZDCEMEnergy(0),
f6536de8	38	fZDCEMEnergy1(0),
2b6cdc43	39	fZDCParticipants(0),
5f0ec064	40	fZDCPartSideA(0),
	41	fZDCPartSideC(0),
	42	fImpactParameter(0),
	43	fImpactParamSideA(0),
	44	fImpactParamSideC(0),
	45	fESDQuality(0)
d5ebf00e	46	{
89291556	47	for(Int_t i=0; i<5; i++){
a85132e7	48	fZN1TowerEnergy[i] = fZN2TowerEnergy[i] = 0.;
	49	fZP1TowerEnergy[i] = fZP2TowerEnergy[i] = 0.;
	50	fZN1TowerEnergyLR[i] = fZN2TowerEnergyLR[i] = 0.;
	51	fZP1TowerEnergyLR[i] = fZP2TowerEnergyLR[i] = 0.;
	52	}
89291556	53	for(Int_t i=0; i<2; i++){
	54	fZNACentrCoord[i] = fZNCCentrCoord[i] = 0.;
	55	}
82dffa48	56	for(Int_t i=0; i<32; i++){
82dffa48	57	fVMEScaler[i]=0;
32e2fda5	58	for(Int_t y=0; y<4; y++) fZDCTDCData[i][y]=0;
82dffa48	59	}
d5ebf00e	60	}
d5ebf00e	61
2b6cdc43	62	//______________________________________________________________________________
d5ebf00e	63	AliESDZDC::AliESDZDC(const AliESDZDC& zdc) :
	64	TObject(zdc),
	65	fZDCN1Energy(zdc.fZDCN1Energy),
	66	fZDCP1Energy(zdc.fZDCP1Energy),
	67	fZDCN2Energy(zdc.fZDCN2Energy),
	68	fZDCP2Energy(zdc.fZDCP2Energy),
f6536de8	69	fZDCEMEnergy(zdc.fZDCEMEnergy),
f6536de8	70	fZDCEMEnergy1(zdc.fZDCEMEnergy1),
2b6cdc43	71	fZDCParticipants(zdc.fZDCParticipants),
5f0ec064	72	fZDCPartSideA(zdc.fZDCPartSideA),
	73	fZDCPartSideC(zdc.fZDCPartSideC),
	74	fImpactParameter(zdc.fImpactParameter),
8c755c03	75	fImpactParamSideA(zdc.fImpactParamSideA),
8c755c03	76	fImpactParamSideC(zdc.fImpactParamSideC),
5f0ec064	77	fESDQuality(zdc.fESDQuality)
d5ebf00e	78	{
d5ebf00e	79	// copy constructor
89291556	80	for(Int_t i=0; i<5; i++){
a85132e7	81	fZN1TowerEnergy[i] = zdc.fZN1TowerEnergy[i];
	82	fZN2TowerEnergy[i] = zdc.fZN2TowerEnergy[i];
	83	fZP1TowerEnergy[i] = zdc.fZP1TowerEnergy[i];
	84	fZP2TowerEnergy[i] = zdc.fZP2TowerEnergy[i];
	85	fZN1TowerEnergyLR[i] = zdc.fZN1TowerEnergyLR[i];
	86	fZN2TowerEnergyLR[i] = zdc.fZN2TowerEnergyLR[i];
	87	fZP1TowerEnergyLR[i] = zdc.fZP1TowerEnergyLR[i];
	88	fZP2TowerEnergyLR[i] = zdc.fZP2TowerEnergyLR[i];
	89	}
89291556	90	for(Int_t i=0; i<2; i++){
	91	fZNACentrCoord[i] = zdc.fZNACentrCoord[i];
	92	fZNCCentrCoord[i] = zdc.fZNCCentrCoord[i];
	93	}
82dffa48	94	for(Int_t i=0; i<32; i++){
82dffa48	95	fVMEScaler[i] = zdc.fVMEScaler[i];
32e2fda5	96	for(Int_t y=0; y<4; y++) fZDCTDCData[i][y] = zdc.fZDCTDCData[i][y];
82dffa48	97	}
d5ebf00e	98	}
d5ebf00e	99
2b6cdc43	100	//______________________________________________________________________________
d5ebf00e	101	AliESDZDC& AliESDZDC::operator=(const AliESDZDC&zdc)
	102	{
	103	// assigment operator
	104	if(this!=&zdc) {
	105	TObject::operator=(zdc);
	106	fZDCN1Energy = zdc.fZDCN1Energy;
	107	fZDCP1Energy = zdc.fZDCP1Energy;
	108	fZDCN2Energy = zdc.fZDCN2Energy;
	109	fZDCP2Energy = zdc.fZDCP2Energy;
f6536de8	110	fZDCEMEnergy = zdc.fZDCEMEnergy;
f6536de8	111	fZDCEMEnergy1 = zdc.fZDCEMEnergy1;
a85132e7	112	for(Int_t i=0; i<5; i++){
7b4458a2	113	fZN1TowerEnergy[i] = zdc.fZN1TowerEnergy[i];
7b4458a2	114	fZN2TowerEnergy[i] = zdc.fZN2TowerEnergy[i];
a85132e7	115	fZP1TowerEnergy[i] = zdc.fZP1TowerEnergy[i];
	116	fZP2TowerEnergy[i] = zdc.fZP2TowerEnergy[i];
	117	fZN1TowerEnergyLR[i] = zdc.fZN1TowerEnergyLR[i];
	118	fZN2TowerEnergyLR[i] = zdc.fZN2TowerEnergyLR[i];
	119	fZP1TowerEnergyLR[i] = zdc.fZP1TowerEnergyLR[i];
	120	fZP2TowerEnergyLR[i] = zdc.fZP2TowerEnergyLR[i];
7b4458a2	121	}
5f0ec064	122	//
	123	fZDCParticipants = zdc.fZDCParticipants;
	124	fZDCPartSideA = zdc.fZDCPartSideA;
	125	fZDCPartSideC = zdc.fZDCPartSideC;
	126	fImpactParameter = zdc.fImpactParameter;
	127	fImpactParamSideA = zdc.fImpactParamSideA;
	128	fImpactParamSideC = zdc.fImpactParamSideC;
	129	//
89291556	130	for(Int_t i=0; i<2; i++){
	131	fZNACentrCoord[i] = zdc.fZNACentrCoord[i];
	132	fZNCCentrCoord[i] = zdc.fZNCCentrCoord[i];
	133	}
5f0ec064	134	//
5f0ec064	135	fESDQuality = zdc.fESDQuality;
82dffa48	136	for(Int_t i=0; i<32; i++){
82dffa48	137	fVMEScaler[i] = zdc.fVMEScaler[i];
32e2fda5	138	for(Int_t y=0; y<4; y++) fZDCTDCData[i][y] = zdc.fZDCTDCData[i][y];
82dffa48	139	}
d5ebf00e	140	}
	141	return *this;
	142	}
	143
2b6cdc43	144	//______________________________________________________________________________
732a24fe	145	void AliESDZDC::Copy(TObject &obj) const {
	146
	147	// this overwrites the virtual TOBject::Copy()
	148	// to allow run time copying without casting
	149	// in AliESDEvent
	150
	151	if(this==&obj)return;
	152	AliESDZDC robj = dynamic_cast<AliESDZDC>(&obj);
	153	if(!robj)return; // not an AliESDZDC
	154	robj = this;
	155
	156	}
	157
d5ebf00e	158
	159	//______________________________________________________________________________
	160	void AliESDZDC::Reset()
	161	{
	162	// reset all data members
	163	fZDCN1Energy=0;
	164	fZDCP1Energy=0;
	165	fZDCN2Energy=0;
	166	fZDCP2Energy=0;
f6536de8	167	fZDCEMEnergy=0;
f6536de8	168	fZDCEMEnergy1=0;
89291556	169	for(Int_t i=0; i<5; i++){
a85132e7	170	fZN1TowerEnergy[i] = fZN2TowerEnergy[i] = 0.;
	171	fZP1TowerEnergy[i] = fZP2TowerEnergy[i] = 0.;
	172	fZN1TowerEnergyLR[i] = fZN2TowerEnergyLR[i] = 0.;
	173	fZP1TowerEnergyLR[i] = fZP2TowerEnergyLR[i] = 0.;
	174	}
5f0ec064	175	fZDCParticipants=0;
	176	fZDCPartSideA=0;
	177	fZDCPartSideC=0;
	178	fImpactParameter=0;
	179	fImpactParamSideA=0;
	180	fImpactParamSideC=0;
89291556	181	for(Int_t i=0; i<2; i++){
	182	fZNACentrCoord[i] = fZNCCentrCoord[i] = 0.;
	183	}
5f0ec064	184	fESDQuality=0;
82dffa48	185	for(Int_t i=0; i<32; i++){
82dffa48	186	fVMEScaler[i] = 0;
32e2fda5	187	for(Int_t y=0; y<4; y++) fZDCTDCData[i][y] = 0;
82dffa48	188	}
d5ebf00e	189	}
	190
	191	//______________________________________________________________________________
	192	void AliESDZDC::Print(const Option_t *) const
	193	{
a85132e7	194	// Print ESD for the ZDC
5f0ec064	195	printf("\n \t E_{ZNC} = %f TeV, E_{ZPC} = %f TeV, E_{ZNA} = %f TeV, E_{ZPA} = %f TeV,"
	196	" E_{ZEM} = %f GeV, Npart = %d, b = %1.2f fm\n",
	197	fZDCN1Energy/1000.,fZDCP1Energy/1000.,fZDCN2Energy/1000.,fZDCP2Energy/1000.,
	198	fZDCEMEnergy+fZDCEMEnergy1, fZDCParticipants,fImpactParameter);
54c6a92c	199	//
82dffa48	200	printf(" ### VMEScaler (!=0): \n");
	201	for(Int_t i=0; i<32; i++) if(fVMEScaler[i]!=0) printf("\t %d \n",fVMEScaler[i]);
	202	printf("\n");
	203	//
	204	printf(" ### TDCData (!=0): \n");
32e2fda5	205	for(Int_t i=0; i<32; i++){
	206	for(Int_t j=0; j<4; j++)
	207	if(fZDCTDCData[i][j]!=0) printf("\t %d \n",fZDCTDCData[i][j]);
	208	}
54c6a92c	209	printf("\n");
d5ebf00e	210	}
d5ebf00e	211
d1d3e36b	212	//______________________________________________________________________________
789b8a8d	213	Bool_t AliESDZDC::GetZNCentroidInPbPb(Float_t beamEne, Double_t centrZNC[2], Double_t centrZNA[2])
d1d3e36b	214	{
d1d3e36b	215	// Provide coordinates of centroid over ZN (side C) front face
789b8a8d	216	if(beamEne==0){
	217	printf(" ZDC centroid in PbPb can't be calculated with E_beam = 0 !!!\n");
	218	for(Int_t jj=0; jj<2; jj++) fZNCCentrCoord[jj] = 999.;
	219	return kFALSE;
	220	}
	221
b19cc8d0	222	const Float_t x[4] = {-1.75, 1.75, -1.75, 1.75};
b19cc8d0	223	const Float_t y[4] = {-1.75, -1.75, 1.75, 1.75};
b19cc8d0	224	const Float_t alpha=0.395;
789b8a8d	225	Float_t numXZNC=0., numYZNC=0., denZNC=0., cZNC, wZNC;
789b8a8d	226	Float_t numXZNA=0., numYZNA=0., denZNA=0., cZNA, wZNA;
5f0ec064	227	//
789b8a8d	228	for(Int_t i=0; i<4; i++){
b19cc8d0	229	if(fZN1TowerEnergy[i+1]>0.) {
789b8a8d	230	wZNC = TMath::Power(fZN1TowerEnergy[i+1], alpha);
	231	numXZNC += x[i]*wZNC;
	232	numYZNC += y[i]*wZNC;
	233	denZNC += wZNC;
b19cc8d0	234	}
b19cc8d0	235	if(fZN2TowerEnergy[i+1]>0.) {
6dec52cf	236	wZNA = TMath::Power(fZN2TowerEnergy[i+1], alpha);
789b8a8d	237	numXZNA += x[i]*wZNA;
	238	numYZNA += y[i]*wZNA;
	239	denZNA += wZNA;
b19cc8d0	240	}
789b8a8d	241	}
ce099621	242	//
789b8a8d	243	if(denZNC!=0){
	244	Float_t nSpecnC = fZDCN1Energy/beamEne;
	245	cZNC = 1.89358-0.71262/(nSpecnC+0.71789);
	246	fZNCCentrCoord[0] = cZNC*numXZNC/denZNC;
	247	fZNCCentrCoord[1] = cZNC*numYZNC/denZNC;
	248	}
	249	else{
	250	fZNCCentrCoord[0] = fZNCCentrCoord[1] = 999.;
b19cc8d0	251	}
789b8a8d	252	if(denZNA!=0){
	253	Float_t nSpecnA = fZDCN1Energy/beamEne;
	254	cZNA = 1.89358-0.71262/(nSpecnA+0.71789);
db3c6990	255	fZNACentrCoord[0] = cZNA*numXZNA/denZNA;
db3c6990	256	fZNACentrCoord[1] = cZNA*numYZNA/denZNA;
789b8a8d	257	}
	258	else{
	259	fZNACentrCoord[0] = fZNACentrCoord[1] = 999.;
	260	}
	261	//
	262	for(Int_t il=0; il<2; il++){
	263	centrZNC[il] = fZNCCentrCoord[il];
	264	centrZNA[il] = fZNACentrCoord[il];
	265	}
	266
	267	return kTRUE;
b19cc8d0	268	}
	269
	270	//______________________________________________________________________________
789b8a8d	271	Bool_t AliESDZDC::GetZNCentroidInpp(Double_t centrZNC[2], Double_t centrZNA[2])
b19cc8d0	272	{
	273	// Provide coordinates of centroid over ZN (side C) front face
	274	const Float_t x[4] = {-1.75, 1.75, -1.75, 1.75};
	275	const Float_t y[4] = {-1.75, -1.75, 1.75, 1.75};
b19cc8d0	276	const Float_t alpha=0.5;
789b8a8d	277	Float_t numXZNC=0., numYZNC=0., denZNC=0., wZNC;
789b8a8d	278	Float_t numXZNA=0., numYZNA=0., denZNA=0., wZNA;
b19cc8d0	279	//
789b8a8d	280	for(Int_t i=0; i<4; i++){
b19cc8d0	281	if(fZN1TowerEnergy[i+1]>0.) {
789b8a8d	282	wZNC = TMath::Power(fZN1TowerEnergy[i+1], alpha);
	283	numXZNC += x[i]*wZNC;
	284	numYZNC += y[i]*wZNC;
	285	denZNC += wZNC;
b19cc8d0	286	}
b19cc8d0	287	if(fZN2TowerEnergy[i+1]>0.) {
789b8a8d	288	wZNA = TMath::Power(fZN2TowerEnergy[i+1], alpha);
	289	numXZNA += x[i]*wZNA;
	290	numYZNA += y[i]*wZNA;
	291	denZNA += wZNA;
b19cc8d0	292	}
789b8a8d	293	}
	294	//
	295	if(denZNC!=0){
	296	fZNCCentrCoord[0] = numXZNC/denZNC;
	297	fZNCCentrCoord[1] = numYZNC/denZNC;
	298	}
	299	else{
	300	fZNCCentrCoord[0] = fZNCCentrCoord[1] = 999.;
	301	}
	302	if(denZNA!=0){
	303	fZNACentrCoord[0] = numXZNA/denZNA;
	304	fZNACentrCoord[1] = numYZNA/denZNA;
	305	}
	306	else{
	307	fZNACentrCoord[0] = fZNACentrCoord[1] = 999.;
	308	}
b19cc8d0	309	//
789b8a8d	310	for(Int_t il=0; il<2; il++){
	311	centrZNC[il] = fZNCCentrCoord[il];
	312	centrZNA[il] = fZNACentrCoord[il];
ce099621	313	}
789b8a8d	314
789b8a8d	315	return kTRUE;
d1d3e36b	316	}