//__________________________________________________________________________________________________
void AliRICHClusterFinder::FindRawClusters(Int_t iChamber)
{//finds neighbours and fill the tree with raw clusters
- Info("FindRawClusters","Start for Chamber %i.",iChamber);
-
Int_t nDigits=Rich()->Digits(iChamber)->GetEntriesFast();
+ Info("FindRawClusters","Start for Chamber %i with %i digits.",iChamber,nDigits);
if(nDigits==0)return;
fHitMap=new AliRICHMap(Rich()->Digits(iChamber));
AliRICHcluster *pRawCluster;
-
+
for(Int_t iDig=0;iDig<nDigits;iDig++){
- AliRICHdigit *dig=(AliRICHdigit*)Rich()->Digits(iChamber)->UncheckedAt(iDig);
+ AliRICHdigit *dig=(AliRICHdigit*)Rich()->Digits(iChamber)->At(iDig);
Int_t i=dig->X(); Int_t j=dig->Y();
if(fHitMap->TestHit(i,j)==kUsed) continue;
**************************************************************************/
#include "AliRICHMap.h"
-#include "AliRICHParam.h"
-#include "AliRICHDigit.h"
+#include "AliRICH.h"
ClassImp(AliRICHMap)
{//Fill hits into HitMap
if(!fNdigits) return;
for(Int_t iDigN=0;iDigN<fNdigits;iDigN++){
- AliRICHDigit *pDig= (AliRICHDigit*)fDigits->At(iDigN);
- SetHit(pDig->PadX(),pDig->PadY(),iDigN);
+ AliRICHdigit *pDig= (AliRICHdigit*)fDigits->At(iDigN);
+ SetHit(pDig->X(),pDig->Y(),iDigN);
}
}
//__________________________________________________________________________________________________
**************************************************************************/
#include "AliRICHParam.h"
-Bool_t AliRICHParam::fgIsWireSag =kTRUE;
-Double_t AliRICHParam::fgAngleRot =-60;
-Int_t AliRICHParam::fgHV =2150;
+Bool_t AliRICHParam::fgIsWireSag =kTRUE;
+Bool_t AliRICHParam::fgIsResolveClusters =kFALSE;
+Double_t AliRICHParam::fgAngleRot =-60;
+Int_t AliRICHParam::fgHV =2150;
ClassImp(AliRICHParam)
static Double_t QdcSlope(Int_t sec) {HV(sec);return 27;}
static Double_t AlphaFeedback(Int_t sec) {HV(sec);return 0.036;}
+ static Bool_t IsResolveClusters() {return fgIsResolveClusters;}
static Bool_t IsWireSag() {return fgIsWireSag;}
static Int_t HV(Int_t) {return fgHV;}
- static Double_t AngleRot() {return fgAngleRot*TMath::DegToRad();}
- static void SetWireSag(Bool_t status) {fgIsWireSag=status;}
- static void SetHV(Int_t hv) {fgHV =hv;}
- static void SetAngleRot(Double_t rot) {fgAngleRot =rot;}
+ static Double_t AngleRot() {return fgAngleRot*TMath::DegToRad();}
+ static void SetResolveClusters(Bool_t a){fgIsResolveClusters=a;}
+ static void SetWireSag(Bool_t status) {fgIsWireSag=status;}
+ static void SetHV(Int_t hv) {fgHV =hv;}
+ static void SetAngleRot(Double_t rot) {fgAngleRot =rot;}
inline static Double_t Mathienson(Double_t lx1,Double_t lx2,Double_t ly1,Double_t ly2);
inline static void Loc2Area(TVector3 hitX3,Int_t &padxMin,Int_t &padyMin,Int_t &padxMax,Int_t &padyMax);
Bool_t SigGenCond(Double_t,Double_t){return kFALSE;}
inline static Int_t Loc2Sec(Double_t &x,Double_t &y);
inline static Int_t Pad2Sec(Int_t &padx,Int_t &pady);
- inline static Bool_t IsResolveClusters() {return kTRUE;}
+
protected:
- static Bool_t fgIsWireSag; //is wire sagitta taken into account
- static Int_t fgHV; //HV applied to anod wires
- static Double_t fgAngleRot; //rotation of RICH from up postion (0,0,490)cm
+ static Bool_t fgIsWireSag; //is wire sagitta taken into account
+ static Bool_t fgIsResolveClusters; //performs declustering or not
+ static Int_t fgHV; //HV applied to anod wires
+ static Double_t fgAngleRot; //rotation of RICH from up postion (0,0,490)cm
ClassDef(AliRICHParam,4) //RICH main parameters
};