fTriggerClasses(kNTriggerClasses),
fDetInDAQ(0),
fDetInReco(0)
+
{
for (Int_t i=0; i<2; i++) fDiamondXY[i]=0.;
fDiamondCovXY[0]=fDiamondCovXY[2]=3.*3.;
fDiamondCovXY[1]=0.;
fTriggerClasses.SetOwner(kTRUE);
+ fMeanBeamInt[0][0]=fMeanBeamInt[0][1]=fMeanBeamInt[1][0]=fMeanBeamInt[1][1]=-1;
for (Int_t m=0; m<kNPHOSMatrix; m++) fPHOSMatrix[m]=NULL;
for (Int_t sm=0; sm<kNEMCALMatrix; sm++) fEMCALMatrix[sm]=NULL;
+ for (Int_t i=0; i<kT0spreadSize;i++) fT0spread[i]=0.;
+ for (Int_t it=0; it<8; it++) fCaloTriggerType[it]=0;
}
//______________________________________________________________________________
else
fPHOSMatrix[m]=NULL;
}
-
+
+ for (int ib=2;ib--;) for (int it=2;it--;) fMeanBeamInt[ib][it] = esd.fMeanBeamInt[ib][it];
+
for(Int_t sm=0; sm<kNEMCALMatrix; sm++){
if(esd.fEMCALMatrix[sm])
fEMCALMatrix[sm]=new TGeoHMatrix(*(esd.fEMCALMatrix[sm])) ;
else
fEMCALMatrix[sm]=NULL;
}
+ for (Int_t i=0; i<kT0spreadSize;i++) fT0spread[i]=esd.fT0spread[i];
+ for (Int_t it=0; it<8; it++) fCaloTriggerType[it]=esd.fCaloTriggerType[it];
+
}
//______________________________________________________________________________
fDetInDAQ = esd.fDetInDAQ;
fDetInReco = esd.fDetInReco;
+ for (int ib=2;ib--;) for (int it=2;it--;) fMeanBeamInt[ib][it] = esd.fMeanBeamInt[ib][it];
+
for(Int_t m=0; m<kNPHOSMatrix; m++){
if(esd.fPHOSMatrix[m])
fPHOSMatrix[m]=new TGeoHMatrix(*(esd.fPHOSMatrix[m])) ;
fEMCALMatrix[sm]=0;
}
}
+ for (Int_t i=0; i<kT0spreadSize;i++) fT0spread[i]=esd.fT0spread[i];
+ for (Int_t it=0; it<8; it++) fCaloTriggerType[it]=esd.fCaloTriggerType[it];
return *this;
}
void AliESDRun::SetDiamond(const AliESDVertex *vertex) {
// set the interaction diamond
- fDiamondXY[0]=vertex->GetXv();
- fDiamondXY[1]=vertex->GetYv();
- fDiamondZ=vertex->GetZv();
- Double32_t cov[6];
- vertex->GetCovMatrix(cov);
- fDiamondCovXY[0]=cov[0];
- fDiamondCovXY[1]=cov[1];
- fDiamondCovXY[2]=cov[2];
- fDiamondSig2Z=cov[5];
+ if (vertex) {
+ fDiamondXY[0]=vertex->GetXv();
+ fDiamondXY[1]=vertex->GetYv();
+ fDiamondZ=vertex->GetZv();
+ Double32_t cov[6];
+ vertex->GetCovMatrix(cov);
+ fDiamondCovXY[0]=cov[0];
+ fDiamondCovXY[1]=cov[1];
+ fDiamondCovXY[2]=cov[2];
+ fDiamondSig2Z=cov[5];
+ }
}
TNamed *str = (TNamed *)((fTriggerClasses).At(i));
if (str) printf("%s ",str->GetName());
}
+ printf("Mean intenstity for interacting : beam1:%+.3e beam2:%+.3e\n",fMeanBeamInt[0][0],fMeanBeamInt[1][0]);
+ printf("Mean intenstity for non-intecting : beam1:%+.3e beam2:%+.3e\n",fMeanBeamInt[0][1],fMeanBeamInt[1][1]);
printf("\n");
}
}
//
}
+
+//_____________________________________________________________________________
+void AliESDRun::SetT0spread(Int_t i,Float_t t)
+{
+ //
+ // Setting the T0 spread value at index i
+ //
+
+ if ( (i>=0) && (i<kT0spreadSize)) {
+ fT0spread[i]=t;
+ } else {
+ AliError(Form("Index %d out of bound",i));
+ }
+ return;
+}
+
+//_____________________________________________________________________________
+void AliESDRun::SetT0spread(Float_t *t)
+{
+ //
+ // Setting the T0 spread values
+ //
+ if (t == 0x0){
+ AliError(Form("Null pointer passed"));
+ }
+ else{
+ for (Int_t i=0;i<kT0spreadSize;i++) fT0spread[i]=t[i];
+ }
+ return;
+}
+