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

/* $Id$ */

//-------------------------------------------------------------------------
//     AOD event base class
//     Author: Markus Oldenburg, CERN
//-------------------------------------------------------------------------

#include "AliAODHeader.h"

ClassImp(AliAODHeader)

//______________________________________________________________________________
AliAODHeader::AliAODHeader() : 
  AliVHeader(),
  fMagneticField(-999.),
  fMuonMagFieldScale(-999.),
  fCentrality(-999.),
  fZDCN1Energy(-999.),
  fZDCP1Energy(-999.),
  fZDCN2Energy(-999.),
  fZDCP2Energy(-999.),
  fNQTheta(0),
  fQTheta(0x0),
  fTriggerMask(0),
  fRunNumber(-999),  
  fRefMult(-999),
  fRefMultPos(-999),
  fRefMultNeg(-999),
  fEventType(0),
  fOrbitNumber(0),
  fPeriodNumber(0),
  fBunchCrossNumber(0),
  fTriggerCluster(0)
{
  // default constructor

  SetName("header");
  for(int j=0; j<2; j++) fZDCEMEnergy[j] = -999.;
}

//______________________________________________________________________________
AliAODHeader::AliAODHeader(Int_t nRun, 
			   UShort_t nBunchX,
			   UInt_t nOrbit,
			   UInt_t nPeriod,
			   Char_t *title) :
  AliVHeader(),
  fMagneticField(-999.),
  fMuonMagFieldScale(-999.),
  fCentrality(-999.),
  fZDCN1Energy(-999.),
  fZDCP1Energy(-999.),
  fZDCN2Energy(-999.),
  fZDCP2Energy(-999.),
  fNQTheta(0),
  fQTheta(0x0),
  fTriggerMask(0),
  fRunNumber(nRun),
  fRefMult(-999),
  fRefMultPos(-999),
  fRefMultNeg(-999),
  fEventType(0),
  fOrbitNumber(nOrbit),
  fPeriodNumber(nPeriod),
  fBunchCrossNumber(nBunchX),
  fTriggerCluster(0)
{
  // constructor

  SetName("header");
  SetTitle(title);
  for(int j=0; j<2; j++) fZDCEMEnergy[j] = -999.;
}

//______________________________________________________________________________
AliAODHeader::AliAODHeader(Int_t nRun, 
			   UShort_t nBunchX,
			   UInt_t nOrbit,
			   UInt_t nPeriod,
			   Int_t refMult,
			   Int_t refMultPos,
			   Int_t refMultNeg,
			   Double_t magField,
			   Double_t muonMagFieldScale,
			   Double_t cent,
			   Double_t n1Energy,
			   Double_t p1Energy,
			   Double_t n2Energy,
			   Double_t p2Energy,
			   Double_t *emEnergy,
 			   ULong64_t trigMask,
			   UChar_t trigClus,
			   UInt_t evttype,
			   Char_t *title) :
  AliVHeader(),
  fMagneticField(magField),
  fMuonMagFieldScale(muonMagFieldScale),
  fCentrality(cent),
  fZDCN1Energy(n1Energy),
  fZDCP1Energy(p1Energy),
  fZDCN2Energy(n2Energy),
  fZDCP2Energy(p2Energy),
  fNQTheta(0),
  fQTheta(0x0),
  fTriggerMask(trigMask),
  fRunNumber(nRun),  
  fRefMult(refMult),
  fRefMultPos(refMultPos),
  fRefMultNeg(refMultNeg),
  fEventType(evttype),
  fOrbitNumber(nOrbit),
  fPeriodNumber(nPeriod),
  fBunchCrossNumber(nBunchX),
  fTriggerCluster(trigClus)
{
  // constructor

  SetName("header");
  SetTitle(title);
  for(int j=0; j<2; j++) fZDCEMEnergy[j] = emEnergy[j];
}

//______________________________________________________________________________
AliAODHeader::~AliAODHeader() 
{
  // destructor
  
  RemoveQTheta();
}

//______________________________________________________________________________
AliAODHeader::AliAODHeader(const AliAODHeader& hdr) :
  AliVHeader(hdr),
  fMagneticField(hdr.fMagneticField),
  fMuonMagFieldScale(hdr.fMuonMagFieldScale),
  fCentrality(hdr.fCentrality),
  fZDCN1Energy(hdr.fZDCN1Energy),
  fZDCP1Energy(hdr.fZDCP1Energy),
  fZDCN2Energy(hdr.fZDCN2Energy),
  fZDCP2Energy(hdr.fZDCP2Energy),
  fNQTheta(0),
  fQTheta(0x0),
  fTriggerMask(hdr.fTriggerMask),
  fRunNumber(hdr.fRunNumber),  
  fRefMult(hdr.fRefMult), 
  fRefMultPos(hdr.fRefMultPos), 
  fRefMultNeg(hdr.fRefMultNeg),
  fEventType(hdr.fEventType),
  fOrbitNumber(hdr.fOrbitNumber),
  fPeriodNumber(hdr.fPeriodNumber),
  fBunchCrossNumber(hdr.fBunchCrossNumber),
  fTriggerCluster(hdr.fTriggerCluster)
{
  // Copy constructor.
  
  SetName(hdr.fName);
  SetTitle(hdr.fTitle);
  SetQTheta(hdr.fQTheta, hdr.fNQTheta);
  SetZDCEMEnergy(hdr.fZDCEMEnergy[0], hdr.fZDCEMEnergy[1]);
}

//______________________________________________________________________________
AliAODHeader& AliAODHeader::operator=(const AliAODHeader& hdr)
{
  // Assignment operator
  if(this!=&hdr) {
    
     AliVHeader::operator=(hdr);
    
    fMagneticField    = hdr.fMagneticField;
    fMuonMagFieldScale= hdr.fMuonMagFieldScale;
    fCentrality       = hdr.fCentrality;
    fZDCN1Energy      = hdr.fZDCN1Energy;
    fZDCP1Energy      = hdr.fZDCP1Energy;
    fZDCN2Energy      = hdr.fZDCN2Energy;
    fZDCP2Energy      = hdr.fZDCP2Energy;
    fTriggerMask      = hdr.fTriggerMask;
    fRunNumber        = hdr.fRunNumber;
    fRefMult          = hdr.fRefMult;
    fRefMultPos       = hdr.fRefMultPos;
    fRefMultNeg       = hdr.fRefMultNeg;
    fEventType        = hdr.fEventType;
    fOrbitNumber      = hdr.fOrbitNumber;
    fPeriodNumber     = hdr.fPeriodNumber;
    fBunchCrossNumber = hdr.fBunchCrossNumber;
    fTriggerCluster   = hdr.fTriggerCluster;

    SetName(hdr.fName);
    SetTitle(hdr.fTitle);
    SetQTheta(hdr.fQTheta, hdr.fNQTheta);
    SetZDCEMEnergy(hdr.fZDCEMEnergy[0], hdr.fZDCEMEnergy[1]);
  }


  return *this;
}

//______________________________________________________________________________
void AliAODHeader::SetQTheta(Double_t *QTheta, UInt_t size) 
{
  if (QTheta && size>0) {
    if (size != (UInt_t)fNQTheta) {
      RemoveQTheta();
      fNQTheta = size;
      fQTheta = new Double_t[fNQTheta];
    }
    
    for (Int_t i = 0; i < fNQTheta; i++) {
      fQTheta[i] = QTheta[i];
    }
  } else {
    RemoveQTheta();
  }

  return;
}

//______________________________________________________________________________
Double_t AliAODHeader::GetQTheta(UInt_t i) const
{
  if (fQTheta && i < (UInt_t)fNQTheta) {
    return fQTheta[i];
  } else {
    return -999.;
  }
}

//______________________________________________________________________________
void AliAODHeader::RemoveQTheta()
{
  delete[] fQTheta;
  fQTheta = 0x0;
  fNQTheta = 0;

  return;
}

//______________________________________________________________________________
void AliAODHeader::Print(Option_t* /*option*/) const 
{
  // prints event information

  printf("Run #                   : %d\n", fRunNumber);
  printf("Bunch Crossing  #       : %d\n", fBunchCrossNumber);
  printf("Orbit Number #          : %d\n", fOrbitNumber);
  printf("Period Number #         : %d\n", fPeriodNumber);
  printf("Trigger mask            : %lld\n", fTriggerMask);
  printf("Trigger cluster         : %d\n", fTriggerCluster);
  printf("Event Type              : %d\n", fEventType);
  printf("Magnetic field          : %f\n", fMagneticField);
  printf("Muon mag. field scale   : %f\n", fMuonMagFieldScale);
  
  printf("Centrality              : %f\n", fCentrality);
  printf("ZDC N1 Energy           : %f\n", fZDCN1Energy);
  printf("ZDC P1 Energy           : %f\n", fZDCP1Energy);
  printf("ZDC N2 Energy           : %f\n", fZDCN2Energy);
  printf("ZDC P2 Energy           : %f\n", fZDCP2Energy);
  printf("ZDC EM1 Energy          : %f\n", fZDCEMEnergy[0]);
  printf("ZDC EM2 Energy          : %f\n", fZDCEMEnergy[1]);
  printf("ref. Multiplicity       : %d\n", fRefMult);
  printf("ref. Multiplicity (pos) : %d\n", fRefMultPos);
  printf("ref. Multiplicity (neg) : %d\n", fRefMultNeg);

  if (fQTheta) {
    for (UInt_t i = 0; i<(UInt_t)fNQTheta; i++) {
      printf("QTheta[%d]              : %d\n", i, GetQTheta(i));
    }
  }

  return;
}
Commit	Line	Data
df9db588	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	/* $Id$ */
	17
	18	//-------------------------------------------------------------------------
	19	// AOD event base class
	20	// Author: Markus Oldenburg, CERN
	21	//-------------------------------------------------------------------------
	22
	23	#include "AliAODHeader.h"
	24
	25	ClassImp(AliAODHeader)
	26
	27	//______________________________________________________________________________
	28	AliAODHeader::AliAODHeader() :
9ae2e5e6	29	AliVHeader(),
df9db588	30	fMagneticField(-999.),
6b6f8d32	31	fMuonMagFieldScale(-999.),
df9db588	32	fCentrality(-999.),
6b6f8d32	33	fZDCN1Energy(-999.),
	34	fZDCP1Energy(-999.),
	35	fZDCN2Energy(-999.),
	36	fZDCP2Energy(-999.),
ff254193	37	fNQTheta(0),
ff254193	38	fQTheta(0x0),
df9db588	39	fTriggerMask(0),
df9db588	40	fRunNumber(-999),
	41	fRefMult(-999),
	42	fRefMultPos(-999),
	43	fRefMultNeg(-999),
9333290e	44	fEventType(0),
	45	fOrbitNumber(0),
	46	fPeriodNumber(0),
	47	fBunchCrossNumber(0),
df9db588	48	fTriggerCluster(0)
	49	{
	50	// default constructor
9ae2e5e6	51
9ae2e5e6	52	SetName("header");
a85132e7	53	for(int j=0; j<2; j++) fZDCEMEnergy[j] = -999.;
df9db588	54	}
	55
	56	//______________________________________________________________________________
31fd97b2	57	AliAODHeader::AliAODHeader(Int_t nRun,
	58	UShort_t nBunchX,
	59	UInt_t nOrbit,
89cf15db	60	UInt_t nPeriod,
df9db588	61	Char_t *title) :
9ae2e5e6	62	AliVHeader(),
df9db588	63	fMagneticField(-999.),
6b6f8d32	64	fMuonMagFieldScale(-999.),
df9db588	65	fCentrality(-999.),
6b6f8d32	66	fZDCN1Energy(-999.),
	67	fZDCP1Energy(-999.),
	68	fZDCN2Energy(-999.),
	69	fZDCP2Energy(-999.),
ff254193	70	fNQTheta(0),
ff254193	71	fQTheta(0x0),
df9db588	72	fTriggerMask(0),
ff254193	73	fRunNumber(nRun),
df9db588	74	fRefMult(-999),
	75	fRefMultPos(-999),
	76	fRefMultNeg(-999),
9333290e	77	fEventType(0),
	78	fOrbitNumber(nOrbit),
	79	fPeriodNumber(nPeriod),
	80	fBunchCrossNumber(nBunchX),
df9db588	81	fTriggerCluster(0)
	82	{
	83	// constructor
9ae2e5e6	84
	85	SetName("header");
	86	SetTitle(title);
a85132e7	87	for(int j=0; j<2; j++) fZDCEMEnergy[j] = -999.;
df9db588	88	}
	89
	90	//______________________________________________________________________________
31fd97b2	91	AliAODHeader::AliAODHeader(Int_t nRun,
	92	UShort_t nBunchX,
	93	UInt_t nOrbit,
89cf15db	94	UInt_t nPeriod,
df9db588	95	Int_t refMult,
	96	Int_t refMultPos,
	97	Int_t refMultNeg,
	98	Double_t magField,
6b6f8d32	99	Double_t muonMagFieldScale,
df9db588	100	Double_t cent,
6b6f8d32	101	Double_t n1Energy,
	102	Double_t p1Energy,
	103	Double_t n2Energy,
	104	Double_t p2Energy,
a85132e7	105	Double_t *emEnergy,
6b6f8d32	106	ULong64_t trigMask,
df9db588	107	UChar_t trigClus,
	108	UInt_t evttype,
	109	Char_t *title) :
9ae2e5e6	110	AliVHeader(),
df9db588	111	fMagneticField(magField),
6b6f8d32	112	fMuonMagFieldScale(muonMagFieldScale),
df9db588	113	fCentrality(cent),
6b6f8d32	114	fZDCN1Energy(n1Energy),
	115	fZDCP1Energy(p1Energy),
	116	fZDCN2Energy(n2Energy),
	117	fZDCP2Energy(p2Energy),
ff254193	118	fNQTheta(0),
ff254193	119	fQTheta(0x0),
df9db588	120	fTriggerMask(trigMask),
df9db588	121	fRunNumber(nRun),
	122	fRefMult(refMult),
	123	fRefMultPos(refMultPos),
	124	fRefMultNeg(refMultNeg),
9333290e	125	fEventType(evttype),
	126	fOrbitNumber(nOrbit),
	127	fPeriodNumber(nPeriod),
	128	fBunchCrossNumber(nBunchX),
df9db588	129	fTriggerCluster(trigClus)
	130	{
	131	// constructor
9ae2e5e6	132
	133	SetName("header");
	134	SetTitle(title);
a85132e7	135	for(int j=0; j<2; j++) fZDCEMEnergy[j] = emEnergy[j];
df9db588	136	}
	137
	138	//______________________________________________________________________________
	139	AliAODHeader::~AliAODHeader()
	140	{
	141	// destructor
ff254193	142
ff254193	143	RemoveQTheta();
df9db588	144	}
	145
	146	//______________________________________________________________________________
	147	AliAODHeader::AliAODHeader(const AliAODHeader& hdr) :
9ae2e5e6	148	AliVHeader(hdr),
df9db588	149	fMagneticField(hdr.fMagneticField),
6b6f8d32	150	fMuonMagFieldScale(hdr.fMuonMagFieldScale),
df9db588	151	fCentrality(hdr.fCentrality),
6b6f8d32	152	fZDCN1Energy(hdr.fZDCN1Energy),
	153	fZDCP1Energy(hdr.fZDCP1Energy),
	154	fZDCN2Energy(hdr.fZDCN2Energy),
	155	fZDCP2Energy(hdr.fZDCP2Energy),
ff254193	156	fNQTheta(0),
ff254193	157	fQTheta(0x0),
df9db588	158	fTriggerMask(hdr.fTriggerMask),
df9db588	159	fRunNumber(hdr.fRunNumber),
	160	fRefMult(hdr.fRefMult),
	161	fRefMultPos(hdr.fRefMultPos),
	162	fRefMultNeg(hdr.fRefMultNeg),
9333290e	163	fEventType(hdr.fEventType),
	164	fOrbitNumber(hdr.fOrbitNumber),
	165	fPeriodNumber(hdr.fPeriodNumber),
	166	fBunchCrossNumber(hdr.fBunchCrossNumber),
df9db588	167	fTriggerCluster(hdr.fTriggerCluster)
	168	{
	169	// Copy constructor.
9ae2e5e6	170
	171	SetName(hdr.fName);
	172	SetTitle(hdr.fTitle);
ff254193	173	SetQTheta(hdr.fQTheta, hdr.fNQTheta);
a85132e7	174	SetZDCEMEnergy(hdr.fZDCEMEnergy[0], hdr.fZDCEMEnergy[1]);
df9db588	175	}
	176
	177	//______________________________________________________________________________
	178	AliAODHeader& AliAODHeader::operator=(const AliAODHeader& hdr)
	179	{
	180	// Assignment operator
	181	if(this!=&hdr) {
6b6f8d32	182
9ae2e5e6	183	AliVHeader::operator=(hdr);
df9db588	184
31fd97b2	185	fMagneticField = hdr.fMagneticField;
6b6f8d32	186	fMuonMagFieldScale= hdr.fMuonMagFieldScale;
31fd97b2	187	fCentrality = hdr.fCentrality;
6b6f8d32	188	fZDCN1Energy = hdr.fZDCN1Energy;
	189	fZDCP1Energy = hdr.fZDCP1Energy;
	190	fZDCN2Energy = hdr.fZDCN2Energy;
	191	fZDCP2Energy = hdr.fZDCP2Energy;
31fd97b2	192	fTriggerMask = hdr.fTriggerMask;
31fd97b2	193	fRunNumber = hdr.fRunNumber;
	194	fRefMult = hdr.fRefMult;
	195	fRefMultPos = hdr.fRefMultPos;
	196	fRefMultNeg = hdr.fRefMultNeg;
9333290e	197	fEventType = hdr.fEventType;
	198	fOrbitNumber = hdr.fOrbitNumber;
	199	fPeriodNumber = hdr.fPeriodNumber;
	200	fBunchCrossNumber = hdr.fBunchCrossNumber;
31fd97b2	201	fTriggerCluster = hdr.fTriggerCluster;
ff254193	202
	203	SetName(hdr.fName);
	204	SetTitle(hdr.fTitle);
	205	SetQTheta(hdr.fQTheta, hdr.fNQTheta);
a85132e7	206	SetZDCEMEnergy(hdr.fZDCEMEnergy[0], hdr.fZDCEMEnergy[1]);
df9db588	207	}
df9db588	208
9ae2e5e6	209
df9db588	210	return *this;
	211	}
	212
ff254193	213	//______________________________________________________________________________
	214	void AliAODHeader::SetQTheta(Double_t *QTheta, UInt_t size)
	215	{
	216	if (QTheta && size>0) {
	217	if (size != (UInt_t)fNQTheta) {
	218	RemoveQTheta();
	219	fNQTheta = size;
	220	fQTheta = new Double_t[fNQTheta];
	221	}
	222
	223	for (Int_t i = 0; i < fNQTheta; i++) {
	224	fQTheta[i] = QTheta[i];
	225	}
	226	} else {
	227	RemoveQTheta();
	228	}
	229
	230	return;
	231	}
	232
	233	//______________________________________________________________________________
	234	Double_t AliAODHeader::GetQTheta(UInt_t i) const
	235	{
	236	if (fQTheta && i < (UInt_t)fNQTheta) {
	237	return fQTheta[i];
	238	} else {
	239	return -999.;
	240	}
	241	}
	242
	243	//______________________________________________________________________________
	244	void AliAODHeader::RemoveQTheta()
	245	{
	246	delete[] fQTheta;
	247	fQTheta = 0x0;
	248	fNQTheta = 0;
	249
	250	return;
	251	}
	252
df9db588	253	//______________________________________________________________________________
	254	void AliAODHeader::Print(Option_t* /option/) const
	255	{
	256	// prints event information
	257
df9db588	258	printf("Run # : %d\n", fRunNumber);
31fd97b2	259	printf("Bunch Crossing # : %d\n", fBunchCrossNumber);
31fd97b2	260	printf("Orbit Number # : %d\n", fOrbitNumber);
89cf15db	261	printf("Period Number # : %d\n", fPeriodNumber);
df9db588	262	printf("Trigger mask : %lld\n", fTriggerMask);
	263	printf("Trigger cluster : %d\n", fTriggerCluster);
	264	printf("Event Type : %d\n", fEventType);
	265	printf("Magnetic field : %f\n", fMagneticField);
6b6f8d32	266	printf("Muon mag. field scale : %f\n", fMuonMagFieldScale);
df9db588	267
df9db588	268	printf("Centrality : %f\n", fCentrality);
6b6f8d32	269	printf("ZDC N1 Energy : %f\n", fZDCN1Energy);
	270	printf("ZDC P1 Energy : %f\n", fZDCP1Energy);
	271	printf("ZDC N2 Energy : %f\n", fZDCN2Energy);
	272	printf("ZDC P2 Energy : %f\n", fZDCP2Energy);
a85132e7	273	printf("ZDC EM1 Energy : %f\n", fZDCEMEnergy[0]);
a85132e7	274	printf("ZDC EM2 Energy : %f\n", fZDCEMEnergy[1]);
df9db588	275	printf("ref. Multiplicity : %d\n", fRefMult);
	276	printf("ref. Multiplicity (pos) : %d\n", fRefMultPos);
	277	printf("ref. Multiplicity (neg) : %d\n", fRefMultNeg);
	278
ff254193	279	if (fQTheta) {
	280	for (UInt_t i = 0; i<(UInt_t)fNQTheta; i++) {
	281	printf("QTheta[%d] : %d\n", i, GetQTheta(i));
	282	}
	283	}
	284
	285	return;
df9db588	286	}