//____________________________________________________________________________
AliPHOS::~AliPHOS()
{
+ // dtor
delete fHits ;
delete fDigits ;
}
gMC->Gstpar(idtmed[712], "LOSS",3) ;
gMC->Gstpar(idtmed[712], "DRAY",1) ;
// Tracking threshold for photons and electrons in the gas ArC02
- gMC->Gstpar(idtmed[715], "CUTGAM",5.E-7) ;
- gMC->Gstpar(idtmed[715], "CUTELE",5.E-7) ;
- gMC->Gstpar(idtmed[715], "CUTNEU",5.E-7) ;
- gMC->Gstpar(idtmed[715], "CUTHAD",5.E-7) ;
- gMC->Gstpar(idtmed[715], "CUTMUO",5.E-7) ;
- gMC->Gstpar(idtmed[715], "BCUTE",5.E-7) ;
- gMC->Gstpar(idtmed[715], "BCUTM",5.E-7) ;
- gMC->Gstpar(idtmed[715], "DCUTE",5.E-7) ;
- gMC->Gstpar(idtmed[715], "DCUTM",5.E-7) ;
- gMC->Gstpar(idtmed[715], "PPCUTM",5.E-7) ;
+ gMC->Gstpar(idtmed[715], "CUTGAM",1.E-5) ;
+ gMC->Gstpar(idtmed[715], "CUTELE",1.E-5) ;
+ gMC->Gstpar(idtmed[715], "CUTNEU",1.E-5) ;
+ gMC->Gstpar(idtmed[715], "CUTHAD",1.E-5) ;
+ gMC->Gstpar(idtmed[715], "CUTMUO",1.E-5) ;
+ gMC->Gstpar(idtmed[715], "BCUTE",1.E-5) ;
+ gMC->Gstpar(idtmed[715], "BCUTM",1.E-5) ;
+ gMC->Gstpar(idtmed[715], "DCUTE",1.E-5) ;
+ gMC->Gstpar(idtmed[715], "DCUTM",1.E-5) ;
+ gMC->Gstpar(idtmed[715], "PPCUTM",1.E-5) ;
gMC->Gstpar(idtmed[715], "LOSS",2.) ;
gMC->Gstpar(idtmed[715], "DRAY",0.) ;
gMC->Gstpar(idtmed[715], "STRA",2.) ;
+}
+
+//____________________________________________________________________________
+AliPHOSRecPoint::RecPointsList * AliPHOS::EmcRecPoints(Int_t evt)
+{
+ // returns the pointer to the EMCA RecPoints list
+ // if the list is empty, get it from TreeR on the disk file
+
+ AliPHOSRecPoint::RecPointsList * rv = 0 ;
+
+ if ( fEmcRecPoints )
+ rv = fEmcRecPoints ;
+
+ else {
+ fEmcRecPoints = new TClonesArray("AliPHOSEmcRecPoint", 100) ;
+ gAlice->GetEvent(evt) ;
+ TTree * fReconstruct = gAlice->TreeR() ;
+ fReconstruct->SetBranchAddress( "PHOSEmcRP", &fEmcRecPoints) ;
+ fReconstruct->GetEvent(0) ;
+ fEmcRecPoints->Expand( fEmcRecPoints->GetEntries() ) ;
+ rv = fEmcRecPoints ;
+ }
+
+ return rv ;
+
+}
+
+//____________________________________________________________________________
+AliPHOSRecParticle::RecParticlesList * AliPHOS::RecParticles(Int_t evt)
+{
+ // returns the pointer to the RecParticles list
+ // if the list is empty, get it from TreeR on the disk file
+
+ AliPHOSRecParticle::RecParticlesList * rv = 0 ;
+ if ( fRecParticles )
+ rv = fRecParticles ;
+
+ else {
+ fRecParticles = new TClonesArray("AliPHOSRecParticle", 100) ;
+ gAlice->GetEvent(evt) ;
+ TTree * fReconstruct = gAlice->TreeR() ;
+ fReconstruct->SetBranchAddress( "PHOSRP", &fRecParticles) ;
+ fReconstruct->GetEvent(0) ;
+ fRecParticles->Expand( fRecParticles->GetEntries() ) ;
+ rv = fRecParticles ;
+ }
+ return rv ;
+
+}
+
+//____________________________________________________________________________
+AliPHOSRecParticle::RecParticlesList * AliPHOS::TrackSegments(Int_t evt)
+{
+ // returns the pointer to the TrackSegments list
+ // if the list is empty, get it from TreeR on the disk file
+ AliPHOSTrackSegment::TrackSegmentsList * rv = 0 ;
+
+ if ( fTrackSegments )
+ rv = fTrackSegments ;
+
+ else {
+ fTrackSegments = new TClonesArray("AliPHOSTrackSegment", 100) ;
+ gAlice->GetEvent(evt) ;
+ TTree * fReconstruct = gAlice->TreeR() ;
+ fReconstruct->SetBranchAddress( "PHOSTS", &fTrackSegments) ;
+ fReconstruct->GetEvent(0) ;
+ fTrackSegments->Expand( fTrackSegments->GetEntries() ) ;
+ rv = fTrackSegments ;
+ }
+
+ return rv ;
+
}