if (!fLego) BookLego();
if (flag == 0) fLego->Reset();
- Int_t nbytes;
+ Int_t nbytes=0;
//
index, GetNPrimaries()) ;
return 0 ;
}
- else {
- Int_t dummy ;
- AliEMCALGetter * gime = AliEMCALGetter::Instance() ;
+ Int_t dummy ;
+ AliEMCALGetter * gime = AliEMCALGetter::Instance() ;
- Int_t ecaRPindex = dynamic_cast<AliEMCALTrackSegment*>(gime->TrackSegments()->At(GetEMCALTSIndex()))->GetECAIndex();
- Int_t primaryindex = dynamic_cast<AliEMCALTowerRecPoint*>(gime->ECARecPoints()->At(ecaRPindex))->GetPrimaries(dummy)[index] ;
- return gime->Primary(primaryindex) ;
- }
- return 0 ;
+ Int_t ecaRPindex = dynamic_cast<AliEMCALTrackSegment*>(gime->TrackSegments()->At(GetEMCALTSIndex()))->GetECAIndex();
+ Int_t primaryindex = dynamic_cast<AliEMCALTowerRecPoint*>(gime->ECARecPoints()->At(ecaRPindex))->GetPrimaries(dummy)[index] ;
+ return gime->Primary(primaryindex) ;
}
if(fNBinsToScale < 1)
{
- return 0.0;
Error("Scale","Number of bins for scaling is smaller thnan 1");
+ return 0.0;
}
UInt_t nbins = num->GetNbinsX();
if (fNBinsToScale > nbins)
if( (fNBinsToScaleX < 1) || (fNBinsToScaleY < 1) )
{
- return 0.0;
Error("Scale","Number of bins for scaling is smaller thnan 1");
+ return 0.0;
}
UInt_t nbinsX = fNumerator->GetNbinsX();
if (fNBinsToScaleX > nbinsX)
if( (fNBinsToScaleX < 1) || (fNBinsToScaleY < 1) || (fNBinsToScaleZ < 1))
{
- return 0.0;
Error("Scale","Number of bins for scaling is smaller thnan 1");
+ return 0.0;
}
UInt_t nbinsX = fNumerator->GetNbinsX();
if (fNBinsToScaleX > nbinsX)
return 100000000 ;
}
//____________________________________________________________________________
-const Int_t AliPHOSPIDv1::GetCPVBit(AliPHOSEmcRecPoint * emc,AliPHOSRecPoint * cpv,const Int_t effPur, const Float_t e) const
+const Int_t AliPHOSPIDv1::GetCPVBit(AliPHOSEmcRecPoint * emc,AliPHOSRecPoint * cpv,const Int_t effPur, Float_t e) const
{
if(effPur>2 || effPur<0)
Error("GetCPVBit","Invalid Efficiency-Purity choice %d",effPur);
}
//____________________________________________________________________________
-const Int_t AliPHOSPIDv1::GetPrincipalBit(TString particle, const Double_t* p,const Int_t effPur, const Float_t e)const
+const Int_t AliPHOSPIDv1::GetPrincipalBit(TString particle, const Double_t* p,const Int_t effPur, Float_t e)const
{
//Is the particle inside de PCA ellipse?
}
else return 0;//success
}
- else
- {//such an object already exists
- TFolder* fol = dynamic_cast<TFolder*>(obj);
- if (fol == 0x0)
- {
- Error("AddSubFodler(TFolder*, ....",
- "Object named %s already exists in folder %s AND IT IS NOT A FOLDER",newfoldname,infoler);
- return 3;
- }
- else
- {
- Warning("AddSubFodler(TFolder*, ....",
- "Folder named %s already exists in folder %s",newfoldname,infoler);
- return 0;
- }
- }
- return 0; //never reached
+ //such an object already exists
+ TFolder* fol = dynamic_cast<TFolder*>(obj);
+ if (fol == 0x0)
+ {
+ Error("AddSubFodler(TFolder*, ....",
+ "Object named %s already exists in folder %s AND IT IS NOT A FOLDER",newfoldname,infoler);
+ return 3;
+ }
+ Warning("AddSubFodler(TFolder*, ....",
+ "Folder named %s already exists in folder %s",newfoldname,infoler);
+ return 0;
}
TTree* TreeT(){return GetTracksDataLoader()->Tree();} //returns the tree from folder; shortcut method
TTree* TreeP(){return GetRecParticlesDataLoader()->Tree();} //returns the tree from folder; shortcut method
- Int_t LoadHits(Option_t* opt=""){return GetHitsDataLoader()->Load(opt);SetTAddrInDet();}
- Int_t LoadSDigits(Option_t* opt=""){return GetSDigitsDataLoader()->Load(opt);SetTAddrInDet();}
- Int_t LoadDigits(Option_t* opt=""){return GetDigitsDataLoader()->Load(opt);SetTAddrInDet();}
+ Int_t LoadHits(Option_t* opt=""){Int_t status = GetHitsDataLoader()->Load(opt);SetTAddrInDet();return status;}
+ Int_t LoadSDigits(Option_t* opt=""){Int_t status = GetSDigitsDataLoader()->Load(opt);SetTAddrInDet(); return status;}
+ Int_t LoadDigits(Option_t* opt=""){Int_t status = GetDigitsDataLoader()->Load(opt);SetTAddrInDet();return status;}
Int_t LoadRecPoints(Option_t* opt=""){return GetRecPointsDataLoader()->Load(opt);}
Int_t LoadTracks(Option_t* opt=""){return GetTracksDataLoader()->Load(opt);}
Int_t LoadRecParticles(Option_t* opt=""){return GetRecParticlesDataLoader()->Load(opt);}
Error("GetHbtPStatusCode","GetTrackEventIndex returned error");
return 0;
}
- else return fHbtPStatCodes[ev][idx];
+ return fHbtPStatCodes[ev][idx];
- return 0;
}
/*******************************************************************/
Error("GetTrack","GetTrackEventIndex returned error");
return 0x0;
}
- else
- {
- if (fgDebug > 5) Info("GetTrack","Number of Tracks in Event(%d) = %d",ev,cab->GetStack(ev)->GetNprimary());
- return cab->GetStack(ev)->Particle(idx); //safe - in case stack does not exist
- } //GetTrackEventIndex would have returned error
- return 0x0;
+ if (fgDebug > 5) Info("GetTrack","Number of Tracks in Event(%d) = %d",ev,cab->GetStack(ev)->GetNprimary());
+ return cab->GetStack(ev)->Particle(idx);//safe - in case stack does not exist
}
/*******************************************************************/
T.Write();
savedir->cd();
out->Close();
- return 0;
gBenchmark->Stop(name);
gBenchmark->Show(name);
+ return 0;
}
//__________________________________________________________________
fcharge = &(fChargeArray[i*fNPRF]);
return GetPRFActiv(xin);
}
- else{
- //make interpolation from more fore lines
- Int_t i= Int_t(y);
- Float_t res;
- if ((i<0) || (i>=fNYdiv) ) return 0;
- Float_t z0=0;
- Float_t z1=0;
- Float_t z2=0;
- Float_t z3=0;
- if (i>0) {
- fcharge =&(fChargeArray[(i-1)*fNPRF]);
- z0 = GetPRFActiv(xin);
- }
- fcharge =&(fChargeArray[i*fNPRF]);
- z1=GetPRFActiv(xin);
- if ((i+1)<fNYdiv){
- fcharge =&(fChargeArray[(i+1)*fNPRF]);
- z2 = GetPRFActiv(xin);
- }
- if ((i+2)<fNYdiv){
- fcharge =&(fChargeArray[(i+2)*fNPRF]);
- z3 = GetPRFActiv(xin);
- }
- Float_t a,b,c,d,k,l;
- a=z1;
- b=(z2-z0)/2.;
- k=z2-a-b;
- l=(z3-z1)/2.-b;
- d=l-2*k;
- c=k-d;
- Float_t dy=y-Float_t(i);
-
- res = a+b*dy+c*dy*dy+d*dy*dy*dy;
- return res;
- }
- return 0.;
+ //make interpolation from more fore lines
+ Int_t i= Int_t(y);
+ Float_t res;
+ if ((i<0) || (i>=fNYdiv) ) return 0;
+ Float_t z0=0;
+ Float_t z1=0;
+ Float_t z2=0;
+ Float_t z3=0;
+ if (i>0) {
+ fcharge =&(fChargeArray[(i-1)*fNPRF]);
+ z0 = GetPRFActiv(xin);
+ }
+ fcharge =&(fChargeArray[i*fNPRF]);
+ z1=GetPRFActiv(xin);
+ if ((i+1)<fNYdiv){
+ fcharge =&(fChargeArray[(i+1)*fNPRF]);
+ z2 = GetPRFActiv(xin);
+ }
+ if ((i+2)<fNYdiv){
+ fcharge =&(fChargeArray[(i+2)*fNPRF]);
+ z3 = GetPRFActiv(xin);
+ }
+ Float_t a,b,c,d,k,l;
+ a=z1;
+ b=(z2-z0)/2.;
+ k=z2-a-b;
+ l=(z3-z1)/2.-b;
+ d=l-2*k;
+ c=k-d;
+ Float_t dy=y-Float_t(i);
+
+ res = a+b*dy+c*dy*dy+d*dy*dy*dy;
+ return res;
}
fC30=c[6 ]; fC31=c[7 ]; fC32=c32; fC33=c[9 ];
fC40=c[10]; fC41=c[11]; fC42=c42; fC43=c[13]; fC44=c[14];
- if (t.GetStatus()&AliESDtrack::kTIME == 0) return;
+ if ((t.GetStatus()&AliESDtrack::kTIME) == 0) return;
StartTimeIntegral();
Double_t times[10]; t.GetIntegratedTimes(times); SetIntegratedTimes(times);
SetIntegratedLength(t.GetIntegratedLength());