/*
$Log$
+ Revision 1.6 2000/10/02 21:28:12 fca
+ Removal of useless dependecies via forward declarations
+
Revision 1.5 2000/06/30 16:30:28 dibari
Disabled writing to rechits.
//printf("Area de uma elipse com teta 0 e Omega 45:%f",Area(0,45));
+ Int_t track;
- for (Int_t track=0; track<ntracks;track++) {
+ for (track=0; track<ntracks;track++) {
gAlice->ResetHits();
gAlice->TreeH()->GetEvent(track);
TClonesArray *pHits = pRICH->Hits();
//printf ("track %d, theta %f, phi %f, omega %f\n\n\n",track,rechit[0],rechit[1],rechit[2]);
// fill rechits
- //pRICH->AddRecHit(nch-1,rechit);
+ pRICH->AddRecHit3D(nch-1,rechit);
+ printf("Chamber:%d",nch);
}
//printf("\n\n\n\n");
gAlice->TreeR()->Fill();
//Stat_t ndig=TR->GetEntries();
TClonesArray *fRec;
for (i=0;i<kNCH;i++) {
- fRec=pRICH->RecHitsAddress(i);
+ fRec=pRICH->RecHitsAddress3D(i);
int ndig=fRec->GetEntriesFast();
printf ("Chamber %d, rings %d\n",i,ndig);
}
//printf("Number of rec. hits: %d",ndig);
- pRICH->ResetRecHits();
+ pRICH->ResetRecHits3D();
//char hname[30];
//sprintf(hname,"TreeR%d",track);
//gAlice->TreeR()->Write(hname);
/*
$Log$
+ Revision 1.7 2000/10/03 21:44:09 morsch
+ Use AliSegmentation and AliHit abstract base classes.
+
Revision 1.6 2000/10/02 21:28:12 fca
Removal of useless dependecies via forward declarations
#include "AliRICHPadHit.h"
#include "AliRICHDigit.h"
#include "AliRICHRawCluster.h"
-#include "AliRICHRecHit.h"
+#include "AliRICHRecHit1D.h"
#include "AliRun.h"
#include "AliDetector.h"
#include "AliRICH.h"
//printf("Center coordinates:%f %f\n",rechit[3],rechit[4]);
- pRICH->AddRecHit(ich,rechit,fEtaPhotons,padsUsedX,padsUsedY);
+ pRICH->AddRecHit1D(ich,rechit,fEtaPhotons,padsUsedX,padsUsedY);
}
gAlice->TreeR()->Fill();
TClonesArray *fRec;
for (i=0;i<kNCH;i++) {
- fRec=pRICH->RecHitsAddress(i);
+ fRec=pRICH->RecHitsAddress1D(i);
int ndig=fRec->GetEntriesFast();
printf ("Chamber %d, rings %d\n",i,ndig);
}
- pRICH->ResetRecHits();
+ pRICH->ResetRecHits1D();
}
Float_t part;
Float_t pX, pY, pZ;
- printf("Calling TrackParam\n");
+ //printf("Calling TrackParam\n");
gAlice->ResetHits();
TTree *treeH = gAlice->TreeH();
funBkg = tan(etaStepAvg)*TMath::Power((1.+TMath::Power(tan(etaStepAvg),2)),
5.52)-7.803 + 22.02*tan(etaStepAvg);
*/
- thetaSig = asin(nfreon/ngas*sin(etaStepAvg));
+
+ //printf("etaStepAvg: %f, etaStepMax: %f, etaStepMin: %f", etaStepAvg,etaStepMax,etaStepMin);
+
+ thetaSig = TMath::ASin(nfreon/ngas*TMath::Sin(etaStepAvg));
funBkg = tan(thetaSig)*(1.+TMath::Power(tan(thetaSig),2))*nfreon
/ngas*cos(etaStepAvg)/cos(thetaSig);
areaBkg += stepEta*funBkg;