[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"
#include "AliCentrality.h"
#include <TGeoMatrix.h>
#include <TObjString.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),
  fFiredTriggers(),
  fRunNumber(-999),  
  fRefMult(-999),
  fRefMultPos(-999),
  fRefMultNeg(-999),
  fNMuons(0),
  fNDimuons(0),
  fEventType(0),
  fOrbitNumber(0),
  fPeriodNumber(0),
  fBunchCrossNumber(0),
  fTriggerCluster(0), 
  fDiamondZ(0.), 
  fDiamondSig2Z(0.),
  fOfflineTrigger(0),
  fESDFileName(""),
  fEventNumberESDFile(-1),
  fL0TriggerInputs(0),
  fL1TriggerInputs(0),
  fL2TriggerInputs(0),
  fCentralityP(0)
{
  // default constructor

  SetName("header");
  for(int j=0; j<2; j++) fZDCEMEnergy[j] = -999.;
  for(Int_t i=0; i<2; i++) fDiamondXY[i]=0.;
  fDiamondCovXY[0]=fDiamondCovXY[2]=3.*3.;
  fDiamondCovXY[1]=0.;
  for (Int_t m=0; m<kNPHOSMatrix; m++) fPHOSMatrix[m]=NULL;
  for (Int_t sm=0; sm<kNEMCALMatrix; sm++) fEMCALMatrix[sm]=NULL;
}

//______________________________________________________________________________
AliAODHeader::AliAODHeader(Int_t nRun, 
			   UShort_t nBunchX,
			   UInt_t nOrbit,
			   UInt_t nPeriod,
			   const Char_t *title) :
  AliVHeader(),
  fMagneticField(-999.),
  fMuonMagFieldScale(-999.),
  fCentrality(-999.),
  fZDCN1Energy(-999.),
  fZDCP1Energy(-999.),
  fZDCN2Energy(-999.),
  fZDCP2Energy(-999.),
  fNQTheta(0),
  fQTheta(0x0),
  fTriggerMask(0),
  fFiredTriggers(),
  fRunNumber(nRun),
  fRefMult(-999),
  fRefMultPos(-999),
  fRefMultNeg(-999),
  fNMuons(0),
  fNDimuons(0),
  fEventType(0),
  fOrbitNumber(nOrbit),
  fPeriodNumber(nPeriod),
  fBunchCrossNumber(nBunchX),
  fTriggerCluster(0), 
  fDiamondZ(0.), 
  fDiamondSig2Z(0.),
  fOfflineTrigger(0),
  fESDFileName(""),
  fEventNumberESDFile(-1),
  fL0TriggerInputs(0),
  fL1TriggerInputs(0),
  fL2TriggerInputs(0),
  fCentralityP(0)
{
  // constructor

  SetName("header");
  SetTitle(title);
  for(int j=0; j<2; j++) fZDCEMEnergy[j] = -999.;
  for(Int_t i=0; i<2; i++) fDiamondXY[i]=0.;
  fDiamondCovXY[0]=fDiamondCovXY[2]=3.*3.;
  fDiamondCovXY[1]=0.;
  for (Int_t m=0; m<kNPHOSMatrix; m++) fPHOSMatrix[m]=NULL;
  for (Int_t sm=0; sm<kNEMCALMatrix; sm++) fEMCALMatrix[sm]=NULL;
}

//______________________________________________________________________________
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,
			   const Char_t *title,
			   Int_t nMuons,
			   Int_t nDimuons) :
  AliVHeader(),
  fMagneticField(magField),
  fMuonMagFieldScale(muonMagFieldScale),
  fCentrality(cent),
  fZDCN1Energy(n1Energy),
  fZDCP1Energy(p1Energy),
  fZDCN2Energy(n2Energy),
  fZDCP2Energy(p2Energy),
  fNQTheta(0),
  fQTheta(0x0),
  fTriggerMask(trigMask),
  fFiredTriggers(),
  fRunNumber(nRun),  
  fRefMult(refMult),
  fRefMultPos(refMultPos),
  fRefMultNeg(refMultNeg),
  fNMuons(nMuons),
  fNDimuons(nDimuons),
  fEventType(evttype),
  fOrbitNumber(nOrbit),
  fPeriodNumber(nPeriod),
  fBunchCrossNumber(nBunchX),
  fTriggerCluster(trigClus),
  fDiamondZ(0.), 
  fDiamondSig2Z(0.),
  fOfflineTrigger(0),
  fESDFileName(""),
  fEventNumberESDFile(-1),
  fL0TriggerInputs(0),
  fL1TriggerInputs(0),
  fL2TriggerInputs(0),
  fCentralityP(0)
{
  // constructor

  SetName("header");
  SetTitle(title);
  for(int j=0; j<2; j++) fZDCEMEnergy[j] = emEnergy[j];
  for(Int_t i=0; i<2; i++) fDiamondXY[i]=0.;
  fDiamondCovXY[0]=fDiamondCovXY[2]=3.*3.;
  fDiamondCovXY[1]=0.;
  for (Int_t m=0; m<kNPHOSMatrix; m++) fPHOSMatrix[m]=NULL;
  for (Int_t sm=0; sm<kNEMCALMatrix; sm++) fEMCALMatrix[sm]=NULL;
}

//______________________________________________________________________________
AliAODHeader::~AliAODHeader() 
{
  // destructor
  delete fCentralityP;
  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),
  fFiredTriggers(hdr.fFiredTriggers),
  fRunNumber(hdr.fRunNumber),  
  fRefMult(hdr.fRefMult), 
  fRefMultPos(hdr.fRefMultPos), 
  fRefMultNeg(hdr.fRefMultNeg),
  fNMuons(hdr.fNMuons),
  fNDimuons(hdr.fNDimuons),
  fEventType(hdr.fEventType),
  fOrbitNumber(hdr.fOrbitNumber),
  fPeriodNumber(hdr.fPeriodNumber),
  fBunchCrossNumber(hdr.fBunchCrossNumber),
  fTriggerCluster(hdr.fTriggerCluster), 
  fDiamondZ(hdr.fDiamondZ), 
  fDiamondSig2Z(hdr.fDiamondSig2Z),
  fOfflineTrigger(hdr.fOfflineTrigger),
  fESDFileName(hdr.fESDFileName),
  fEventNumberESDFile(hdr.fEventNumberESDFile),
  fL0TriggerInputs(hdr.fL0TriggerInputs),
  fL1TriggerInputs(hdr.fL1TriggerInputs),
  fL2TriggerInputs(hdr.fL2TriggerInputs),
  fCentralityP(hdr.fCentralityP)
{
  // Copy constructor.
  
  SetName(hdr.fName);
  SetTitle(hdr.fTitle);
  SetQTheta(hdr.fQTheta, hdr.fNQTheta);
  SetZDCEMEnergy(hdr.fZDCEMEnergy[0], hdr.fZDCEMEnergy[1]);
  for(Int_t i=0; i<2; i++) fDiamondXY[i]=hdr.fDiamondXY[i];
  for(Int_t i=0; i<3; i++) fDiamondCovXY[i]=hdr.fDiamondCovXY[i];


  for(Int_t m=0; m<kNPHOSMatrix; m++){
      if(hdr.fPHOSMatrix[m])
	  fPHOSMatrix[m]=new TGeoHMatrix(*(hdr.fPHOSMatrix[m])) ;
      else
	  fPHOSMatrix[m]=0;
  }
  
  for(Int_t sm=0; sm<kNEMCALMatrix; sm++){
      if(hdr.fEMCALMatrix[sm])
	  fEMCALMatrix[sm]=new TGeoHMatrix(*(hdr.fEMCALMatrix[sm])) ;
      else
	  fEMCALMatrix[sm]=0;
  }

}

//______________________________________________________________________________
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;
    fFiredTriggers    = hdr.fFiredTriggers;
    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;
    fNMuons           = hdr.fNMuons;
    fNDimuons         = hdr.fNDimuons;
    fDiamondZ         = hdr.fDiamondZ;
    fDiamondSig2Z     = hdr.fDiamondSig2Z;
    fOfflineTrigger   = hdr.fOfflineTrigger;
    fESDFileName      = hdr.fESDFileName;
    fEventNumberESDFile = hdr.fEventNumberESDFile;
    fL0TriggerInputs    = hdr.fL0TriggerInputs;
    fL1TriggerInputs    = hdr.fL1TriggerInputs;
    fL2TriggerInputs    = hdr.fL2TriggerInputs;
    fCentralityP        = hdr.fCentralityP;

    SetName(hdr.fName);
    SetTitle(hdr.fTitle);
    SetQTheta(hdr.fQTheta, hdr.fNQTheta);
    SetZDCEMEnergy(hdr.fZDCEMEnergy[0], hdr.fZDCEMEnergy[1]);
    for(Int_t i=0; i<2; i++) fDiamondXY[i]=hdr.fDiamondXY[i];
    for(Int_t i=0; i<3; i++) fDiamondCovXY[i]=hdr.fDiamondCovXY[i];

    for(Int_t m=0; m<kNPHOSMatrix; m++){
	if(hdr.fPHOSMatrix[m])
	    fPHOSMatrix[m]=new TGeoHMatrix(*(hdr.fPHOSMatrix[m])) ;
	else
	    fPHOSMatrix[m]=0;
    }
    
    for(Int_t sm=0; sm<kNEMCALMatrix; sm++){
	if(hdr.fEMCALMatrix[sm])
	    fEMCALMatrix[sm]=new TGeoHMatrix(*(hdr.fEMCALMatrix[sm])) ;
	else
	    fEMCALMatrix[sm]=0;
  }

  }


  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::Clear(Option_t* /*opt*/)
{
  RemoveQTheta();
  if (fCentralityP) delete fCentralityP;
  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);
  printf("number of muons         : %d\n", fNMuons);
  printf("number of dimuons       : %d\n", fNDimuons);
  printf("offline trigger         : %u\n", fOfflineTrigger);

  if (fQTheta) {
    for (UInt_t i = 0; i<(UInt_t)fNQTheta; i++) {
      printf("QTheta[%d]              : %13.3e\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"
1d7d8ba7	24	#include "AliCentrality.h"
1aa76f71	25	#include <TGeoMatrix.h>
27346f69	26	#include <TObjString.h>
df9db588	27
	28	ClassImp(AliAODHeader)
	29
	30	//______________________________________________________________________________
	31	AliAODHeader::AliAODHeader() :
9ae2e5e6	32	AliVHeader(),
df9db588	33	fMagneticField(-999.),
6b6f8d32	34	fMuonMagFieldScale(-999.),
df9db588	35	fCentrality(-999.),
6b6f8d32	36	fZDCN1Energy(-999.),
	37	fZDCP1Energy(-999.),
	38	fZDCN2Energy(-999.),
	39	fZDCP2Energy(-999.),
ff254193	40	fNQTheta(0),
ff254193	41	fQTheta(0x0),
df9db588	42	fTriggerMask(0),
27346f69	43	fFiredTriggers(),
df9db588	44	fRunNumber(-999),
	45	fRefMult(-999),
	46	fRefMultPos(-999),
	47	fRefMultNeg(-999),
dba4ebde	48	fNMuons(0),
dba4ebde	49	fNDimuons(0),
9333290e	50	fEventType(0),
	51	fOrbitNumber(0),
	52	fPeriodNumber(0),
	53	fBunchCrossNumber(0),
65b25288	54	fTriggerCluster(0),
65b25288	55	fDiamondZ(0.),
0c6c629b	56	fDiamondSig2Z(0.),
a5f7aba4	57	fOfflineTrigger(0),
	58	fESDFileName(""),
	59	fEventNumberESDFile(-1),
	60	fL0TriggerInputs(0),
	61	fL1TriggerInputs(0),
1d7d8ba7	62	fL2TriggerInputs(0),
1d7d8ba7	63	fCentralityP(0)
df9db588	64	{
df9db588	65	// default constructor
9ae2e5e6	66
9ae2e5e6	67	SetName("header");
a85132e7	68	for(int j=0; j<2; j++) fZDCEMEnergy[j] = -999.;
613fc341	69	for(Int_t i=0; i<2; i++) fDiamondXY[i]=0.;
	70	fDiamondCovXY[0]=fDiamondCovXY[2]=3.*3.;
	71	fDiamondCovXY[1]=0.;
1aa76f71	72	for (Int_t m=0; m<kNPHOSMatrix; m++) fPHOSMatrix[m]=NULL;
1aa76f71	73	for (Int_t sm=0; sm<kNEMCALMatrix; sm++) fEMCALMatrix[sm]=NULL;
df9db588	74	}
	75
	76	//______________________________________________________________________________
31fd97b2	77	AliAODHeader::AliAODHeader(Int_t nRun,
	78	UShort_t nBunchX,
	79	UInt_t nOrbit,
89cf15db	80	UInt_t nPeriod,
abfce367	81	const Char_t *title) :
9ae2e5e6	82	AliVHeader(),
df9db588	83	fMagneticField(-999.),
6b6f8d32	84	fMuonMagFieldScale(-999.),
df9db588	85	fCentrality(-999.),
6b6f8d32	86	fZDCN1Energy(-999.),
	87	fZDCP1Energy(-999.),
	88	fZDCN2Energy(-999.),
	89	fZDCP2Energy(-999.),
ff254193	90	fNQTheta(0),
ff254193	91	fQTheta(0x0),
df9db588	92	fTriggerMask(0),
27346f69	93	fFiredTriggers(),
ff254193	94	fRunNumber(nRun),
df9db588	95	fRefMult(-999),
	96	fRefMultPos(-999),
	97	fRefMultNeg(-999),
dba4ebde	98	fNMuons(0),
dba4ebde	99	fNDimuons(0),
9333290e	100	fEventType(0),
	101	fOrbitNumber(nOrbit),
	102	fPeriodNumber(nPeriod),
	103	fBunchCrossNumber(nBunchX),
65b25288	104	fTriggerCluster(0),
65b25288	105	fDiamondZ(0.),
0c6c629b	106	fDiamondSig2Z(0.),
a5f7aba4	107	fOfflineTrigger(0),
	108	fESDFileName(""),
	109	fEventNumberESDFile(-1),
	110	fL0TriggerInputs(0),
	111	fL1TriggerInputs(0),
1d7d8ba7	112	fL2TriggerInputs(0),
1d7d8ba7	113	fCentralityP(0)
df9db588	114	{
df9db588	115	// constructor
9ae2e5e6	116
	117	SetName("header");
	118	SetTitle(title);
a85132e7	119	for(int j=0; j<2; j++) fZDCEMEnergy[j] = -999.;
613fc341	120	for(Int_t i=0; i<2; i++) fDiamondXY[i]=0.;
	121	fDiamondCovXY[0]=fDiamondCovXY[2]=3.*3.;
	122	fDiamondCovXY[1]=0.;
1aa76f71	123	for (Int_t m=0; m<kNPHOSMatrix; m++) fPHOSMatrix[m]=NULL;
1aa76f71	124	for (Int_t sm=0; sm<kNEMCALMatrix; sm++) fEMCALMatrix[sm]=NULL;
df9db588	125	}
	126
	127	//______________________________________________________________________________
31fd97b2	128	AliAODHeader::AliAODHeader(Int_t nRun,
	129	UShort_t nBunchX,
	130	UInt_t nOrbit,
89cf15db	131	UInt_t nPeriod,
df9db588	132	Int_t refMult,
	133	Int_t refMultPos,
	134	Int_t refMultNeg,
	135	Double_t magField,
6b6f8d32	136	Double_t muonMagFieldScale,
df9db588	137	Double_t cent,
6b6f8d32	138	Double_t n1Energy,
	139	Double_t p1Energy,
	140	Double_t n2Energy,
	141	Double_t p2Energy,
a85132e7	142	Double_t *emEnergy,
6b6f8d32	143	ULong64_t trigMask,
df9db588	144	UChar_t trigClus,
df9db588	145	UInt_t evttype,
fa8b0e56	146	const Char_t *title,
	147	Int_t nMuons,
	148	Int_t nDimuons) :
9ae2e5e6	149	AliVHeader(),
df9db588	150	fMagneticField(magField),
6b6f8d32	151	fMuonMagFieldScale(muonMagFieldScale),
df9db588	152	fCentrality(cent),
6b6f8d32	153	fZDCN1Energy(n1Energy),
	154	fZDCP1Energy(p1Energy),
	155	fZDCN2Energy(n2Energy),
	156	fZDCP2Energy(p2Energy),
ff254193	157	fNQTheta(0),
ff254193	158	fQTheta(0x0),
df9db588	159	fTriggerMask(trigMask),
27346f69	160	fFiredTriggers(),
df9db588	161	fRunNumber(nRun),
	162	fRefMult(refMult),
	163	fRefMultPos(refMultPos),
	164	fRefMultNeg(refMultNeg),
dba4ebde	165	fNMuons(nMuons),
dba4ebde	166	fNDimuons(nDimuons),
9333290e	167	fEventType(evttype),
	168	fOrbitNumber(nOrbit),
	169	fPeriodNumber(nPeriod),
	170	fBunchCrossNumber(nBunchX),
0c6c629b	171	fTriggerCluster(trigClus),
65b25288	172	fDiamondZ(0.),
0c6c629b	173	fDiamondSig2Z(0.),
a5f7aba4	174	fOfflineTrigger(0),
	175	fESDFileName(""),
	176	fEventNumberESDFile(-1),
	177	fL0TriggerInputs(0),
	178	fL1TriggerInputs(0),
1d7d8ba7	179	fL2TriggerInputs(0),
1d7d8ba7	180	fCentralityP(0)
df9db588	181	{
df9db588	182	// constructor
9ae2e5e6	183
	184	SetName("header");
	185	SetTitle(title);
a85132e7	186	for(int j=0; j<2; j++) fZDCEMEnergy[j] = emEnergy[j];
613fc341	187	for(Int_t i=0; i<2; i++) fDiamondXY[i]=0.;
	188	fDiamondCovXY[0]=fDiamondCovXY[2]=3.*3.;
	189	fDiamondCovXY[1]=0.;
1aa76f71	190	for (Int_t m=0; m<kNPHOSMatrix; m++) fPHOSMatrix[m]=NULL;
1aa76f71	191	for (Int_t sm=0; sm<kNEMCALMatrix; sm++) fEMCALMatrix[sm]=NULL;
df9db588	192	}
	193
	194	//______________________________________________________________________________
	195	AliAODHeader::~AliAODHeader()
	196	{
	197	// destructor
1d7d8ba7	198	delete fCentralityP;
ff254193	199	RemoveQTheta();
df9db588	200	}
	201
	202	//______________________________________________________________________________
	203	AliAODHeader::AliAODHeader(const AliAODHeader& hdr) :
9ae2e5e6	204	AliVHeader(hdr),
df9db588	205	fMagneticField(hdr.fMagneticField),
6b6f8d32	206	fMuonMagFieldScale(hdr.fMuonMagFieldScale),
df9db588	207	fCentrality(hdr.fCentrality),
6b6f8d32	208	fZDCN1Energy(hdr.fZDCN1Energy),
	209	fZDCP1Energy(hdr.fZDCP1Energy),
	210	fZDCN2Energy(hdr.fZDCN2Energy),
	211	fZDCP2Energy(hdr.fZDCP2Energy),
ff254193	212	fNQTheta(0),
ff254193	213	fQTheta(0x0),
df9db588	214	fTriggerMask(hdr.fTriggerMask),
27346f69	215	fFiredTriggers(hdr.fFiredTriggers),
df9db588	216	fRunNumber(hdr.fRunNumber),
	217	fRefMult(hdr.fRefMult),
	218	fRefMultPos(hdr.fRefMultPos),
	219	fRefMultNeg(hdr.fRefMultNeg),
dba4ebde	220	fNMuons(hdr.fNMuons),
dba4ebde	221	fNDimuons(hdr.fNDimuons),
9333290e	222	fEventType(hdr.fEventType),
	223	fOrbitNumber(hdr.fOrbitNumber),
	224	fPeriodNumber(hdr.fPeriodNumber),
	225	fBunchCrossNumber(hdr.fBunchCrossNumber),
80a3f83f	226	fTriggerCluster(hdr.fTriggerCluster),
80a3f83f	227	fDiamondZ(hdr.fDiamondZ),
0c6c629b	228	fDiamondSig2Z(hdr.fDiamondSig2Z),
a5f7aba4	229	fOfflineTrigger(hdr.fOfflineTrigger),
	230	fESDFileName(hdr.fESDFileName),
	231	fEventNumberESDFile(hdr.fEventNumberESDFile),
	232	fL0TriggerInputs(hdr.fL0TriggerInputs),
	233	fL1TriggerInputs(hdr.fL1TriggerInputs),
1d7d8ba7	234	fL2TriggerInputs(hdr.fL2TriggerInputs),
1d7d8ba7	235	fCentralityP(hdr.fCentralityP)
df9db588	236	{
df9db588	237	// Copy constructor.
9ae2e5e6	238
	239	SetName(hdr.fName);
	240	SetTitle(hdr.fTitle);
ff254193	241	SetQTheta(hdr.fQTheta, hdr.fNQTheta);
a85132e7	242	SetZDCEMEnergy(hdr.fZDCEMEnergy[0], hdr.fZDCEMEnergy[1]);
613fc341	243	for(Int_t i=0; i<2; i++) fDiamondXY[i]=hdr.fDiamondXY[i];
613fc341	244	for(Int_t i=0; i<3; i++) fDiamondCovXY[i]=hdr.fDiamondCovXY[i];
1aa76f71	245
	246
	247	for(Int_t m=0; m<kNPHOSMatrix; m++){
	248	if(hdr.fPHOSMatrix[m])
	249	fPHOSMatrix[m]=new TGeoHMatrix(*(hdr.fPHOSMatrix[m])) ;
	250	else
	251	fPHOSMatrix[m]=0;
	252	}
	253
	254	for(Int_t sm=0; sm<kNEMCALMatrix; sm++){
	255	if(hdr.fEMCALMatrix[sm])
	256	fEMCALMatrix[sm]=new TGeoHMatrix(*(hdr.fEMCALMatrix[sm])) ;
	257	else
	258	fEMCALMatrix[sm]=0;
	259	}
	260
df9db588	261	}
	262
	263	//______________________________________________________________________________
	264	AliAODHeader& AliAODHeader::operator=(const AliAODHeader& hdr)
	265	{
	266	// Assignment operator
	267	if(this!=&hdr) {
6b6f8d32	268
9ae2e5e6	269	AliVHeader::operator=(hdr);
df9db588	270
31fd97b2	271	fMagneticField = hdr.fMagneticField;
6b6f8d32	272	fMuonMagFieldScale= hdr.fMuonMagFieldScale;
31fd97b2	273	fCentrality = hdr.fCentrality;
6b6f8d32	274	fZDCN1Energy = hdr.fZDCN1Energy;
	275	fZDCP1Energy = hdr.fZDCP1Energy;
	276	fZDCN2Energy = hdr.fZDCN2Energy;
	277	fZDCP2Energy = hdr.fZDCP2Energy;
31fd97b2	278	fTriggerMask = hdr.fTriggerMask;
27346f69	279	fFiredTriggers = hdr.fFiredTriggers;
31fd97b2	280	fRunNumber = hdr.fRunNumber;
	281	fRefMult = hdr.fRefMult;
	282	fRefMultPos = hdr.fRefMultPos;
	283	fRefMultNeg = hdr.fRefMultNeg;
9333290e	284	fEventType = hdr.fEventType;
	285	fOrbitNumber = hdr.fOrbitNumber;
	286	fPeriodNumber = hdr.fPeriodNumber;
	287	fBunchCrossNumber = hdr.fBunchCrossNumber;
31fd97b2	288	fTriggerCluster = hdr.fTriggerCluster;
fa8b0e56	289	fNMuons = hdr.fNMuons;
fa8b0e56	290	fNDimuons = hdr.fNDimuons;
0c6c629b	291	fDiamondZ = hdr.fDiamondZ;
	292	fDiamondSig2Z = hdr.fDiamondSig2Z;
	293	fOfflineTrigger = hdr.fOfflineTrigger;
a5f7aba4	294	fESDFileName = hdr.fESDFileName;
	295	fEventNumberESDFile = hdr.fEventNumberESDFile;
	296	fL0TriggerInputs = hdr.fL0TriggerInputs;
	297	fL1TriggerInputs = hdr.fL1TriggerInputs;
	298	fL2TriggerInputs = hdr.fL2TriggerInputs;
1d7d8ba7	299	fCentralityP = hdr.fCentralityP;
ff254193	300
	301	SetName(hdr.fName);
	302	SetTitle(hdr.fTitle);
	303	SetQTheta(hdr.fQTheta, hdr.fNQTheta);
a85132e7	304	SetZDCEMEnergy(hdr.fZDCEMEnergy[0], hdr.fZDCEMEnergy[1]);
613fc341	305	for(Int_t i=0; i<2; i++) fDiamondXY[i]=hdr.fDiamondXY[i];
613fc341	306	for(Int_t i=0; i<3; i++) fDiamondCovXY[i]=hdr.fDiamondCovXY[i];
1aa76f71	307
	308	for(Int_t m=0; m<kNPHOSMatrix; m++){
	309	if(hdr.fPHOSMatrix[m])
	310	fPHOSMatrix[m]=new TGeoHMatrix(*(hdr.fPHOSMatrix[m])) ;
	311	else
	312	fPHOSMatrix[m]=0;
	313	}
	314
	315	for(Int_t sm=0; sm<kNEMCALMatrix; sm++){
	316	if(hdr.fEMCALMatrix[sm])
	317	fEMCALMatrix[sm]=new TGeoHMatrix(*(hdr.fEMCALMatrix[sm])) ;
	318	else
	319	fEMCALMatrix[sm]=0;
	320	}
	321
df9db588	322	}
df9db588	323
9ae2e5e6	324
df9db588	325	return *this;
	326	}
	327
ff254193	328	//______________________________________________________________________________
	329	void AliAODHeader::SetQTheta(Double_t *QTheta, UInt_t size)
	330	{
	331	if (QTheta && size>0) {
	332	if (size != (UInt_t)fNQTheta) {
	333	RemoveQTheta();
	334	fNQTheta = size;
	335	fQTheta = new Double_t[fNQTheta];
	336	}
	337
	338	for (Int_t i = 0; i < fNQTheta; i++) {
	339	fQTheta[i] = QTheta[i];
	340	}
	341	} else {
	342	RemoveQTheta();
	343	}
	344
	345	return;
	346	}
	347
	348	//______________________________________________________________________________
	349	Double_t AliAODHeader::GetQTheta(UInt_t i) const
	350	{
	351	if (fQTheta && i < (UInt_t)fNQTheta) {
	352	return fQTheta[i];
	353	} else {
	354	return -999.;
	355	}
	356	}
	357
	358	//______________________________________________________________________________
	359	void AliAODHeader::RemoveQTheta()
	360	{
	361	delete[] fQTheta;
	362	fQTheta = 0x0;
	363	fNQTheta = 0;
	364
	365	return;
	366	}
	367
99dbf027	368	void AliAODHeader::Clear(Option_t* /opt/)
	369	{
	370	RemoveQTheta();
	371	if (fCentralityP) delete fCentralityP;
	372	return;
	373	}
	374
df9db588	375	//______________________________________________________________________________
	376	void AliAODHeader::Print(Option_t* /option/) const
	377	{
	378	// prints event information
	379
df9db588	380	printf("Run # : %d\n", fRunNumber);
31fd97b2	381	printf("Bunch Crossing # : %d\n", fBunchCrossNumber);
31fd97b2	382	printf("Orbit Number # : %d\n", fOrbitNumber);
89cf15db	383	printf("Period Number # : %d\n", fPeriodNumber);
df9db588	384	printf("Trigger mask : %lld\n", fTriggerMask);
	385	printf("Trigger cluster : %d\n", fTriggerCluster);
	386	printf("Event Type : %d\n", fEventType);
	387	printf("Magnetic field : %f\n", fMagneticField);
6b6f8d32	388	printf("Muon mag. field scale : %f\n", fMuonMagFieldScale);
df9db588	389
df9db588	390	printf("Centrality : %f\n", fCentrality);
6b6f8d32	391	printf("ZDC N1 Energy : %f\n", fZDCN1Energy);
	392	printf("ZDC P1 Energy : %f\n", fZDCP1Energy);
	393	printf("ZDC N2 Energy : %f\n", fZDCN2Energy);
	394	printf("ZDC P2 Energy : %f\n", fZDCP2Energy);
a85132e7	395	printf("ZDC EM1 Energy : %f\n", fZDCEMEnergy[0]);
a85132e7	396	printf("ZDC EM2 Energy : %f\n", fZDCEMEnergy[1]);
df9db588	397	printf("ref. Multiplicity : %d\n", fRefMult);
	398	printf("ref. Multiplicity (pos) : %d\n", fRefMultPos);
	399	printf("ref. Multiplicity (neg) : %d\n", fRefMultNeg);
fa8b0e56	400	printf("number of muons : %d\n", fNMuons);
fa8b0e56	401	printf("number of dimuons : %d\n", fNDimuons);
0c6c629b	402	printf("offline trigger : %u\n", fOfflineTrigger);
df9db588	403
ff254193	404	if (fQTheta) {
ff254193	405	for (UInt_t i = 0; i<(UInt_t)fNQTheta; i++) {
d067f357	406	printf("QTheta[%d] : %13.3e\n", i, GetQTheta(i));
ff254193	407	}
	408	}
	409
	410	return;
df9db588	411	}