#include "AliLog.h"
+#include "AliESDtrack.h"
#include "AliTRDCalibraFillHisto.h"
#include "AliTRDCalibraMode.h"
#include "AliTRDCalibraVector.h"
,fNormalizeNbOfCluster(kFALSE)
,fMaxCluster(0)
,fNbMaxCluster(0)
+ ,fCutWithVdriftCalib(kFALSE)
+ ,fMinNbTRDtracklets(0)
+ ,fMinTRDMomentum(0.0)
,fFirstRunGain(0)
,fVersionGainUsed(0)
,fSubVersionGainUsed(0)
,fNormalizeNbOfCluster(c.fNormalizeNbOfCluster)
,fMaxCluster(c.fMaxCluster)
,fNbMaxCluster(c.fNbMaxCluster)
+ ,fCutWithVdriftCalib(c.fCutWithVdriftCalib)
+ ,fMinNbTRDtracklets(c.fMinNbTRDtracklets)
+ ,fMinTRDMomentum(c.fMinTRDMomentum)
,fFirstRunGain(c.fFirstRunGain)
,fVersionGainUsed(c.fVersionGainUsed)
,fSubVersionGainUsed(c.fSubVersionGainUsed)
}
//____________Offline tracking in the AliTRDtracker____________________________
-Bool_t AliTRDCalibraFillHisto::UpdateHistogramsV1(const AliTRDtrackV1 *t)
+Bool_t AliTRDCalibraFillHisto::UpdateHistogramsV1(const AliTRDtrackV1 *t,const AliESDtrack *esdtrack)
{
//
// Use AliTRDtrackV1 for the calibration
AliTRDcluster *cl = 0x0; // cluster attached now to the tracklet
AliTRDcluster *cls = 0x0; // shared cluster attached now to the tracklet
Bool_t newtr = kTRUE; // new track
-
- // Get cal
- // AliTRDcalibDB *cal = AliTRDcalibDB::Instance();
- /*
- if (!cal) {
- AliInfo("Could not get calibDB");
- return kFALSE;
- }
-*/
+
+
+ //
+ // Cut on the number of TRD tracklets
+ //
+ Int_t numberoftrdtracklets = t->GetNumberOfTracklets();
+ if(numberoftrdtracklets < fMinNbTRDtracklets) return kFALSE;
+
+ //
if (!fCalibDB) {
AliInfo("Could not get calibDB");
return kFALSE;
if(!tracklet->IsOK()) continue;
fNumberTrack++;
ResetfVariablestracklet();
+ Float_t momentum = t->GetMomentum(itr);
+ if(TMath::Abs(momentum) < fMinTRDMomentum) continue;
+
//////////////////////////////////////////
// localisation of the tracklet and dqdl
////////////////////////////
// loop over the clusters
////////////////////////////
+ Double_t chargeQ = 0.0;
Int_t nbclusters = 0;
for(int jc=0; jc<AliTRDseedV1::kNtb; jc++){
if(!(cl = tracklet->GetClusters(jc))) continue;
// Add the charge if shared cluster
cls = tracklet->GetClusters(jc+AliTRDseedV1::kNtb);
//
- StoreInfoCHPHtrack(cl, tracklet->GetQperTB(jc),group,row,col,cls); //tracklet->GetdQdl(jc)
+ chargeQ += StoreInfoCHPHtrack(cl, tracklet->GetQperTB(jc),group,row,col,cls); //tracklet->GetdQdl(jc)
}
////////////////////////////////////////
// Gain calibration
if (fCH2dOn) {
- FillTheInfoOfTheTrackCH(nbclusters);
+ if(fCutWithVdriftCalib) {
+ if(pass) FillTheInfoOfTheTrackCH(nbclusters);
+ } else {
+ FillTheInfoOfTheTrackCH(nbclusters);
+ }
}
// PH calibration
if (fPH2dOn) {
- FillTheInfoOfTheTrackPH();
+ if(fCutWithVdriftCalib) {
+ if(pass) FillTheInfoOfTheTrackPH();
+ }
+ else {
+ FillTheInfoOfTheTrackPH();
+ }
}
if(pass && fPRF2dOn) HandlePRFtrackletV1(tracklet,nbclusters);
-
+
+
+ /////////////////////////////////////////////////////////
+ // Debug
+ ////////////////////////////////////////////////////////
+ if(fDebugLevel > 0){
+ //printf("test\n");
+ if ( !fDebugStreamer ) {
+ //debug stream
+ TDirectory *backup = gDirectory;
+ fDebugStreamer = new TTreeSRedirector("TRDdebugCalibraFill.root");
+ if ( backup ) backup->cd(); //we don't want to be cd'd to the debug streamer
+ }
+
+ Int_t stacke = AliTRDgeometry::GetStack(detector);
+ Int_t sme = AliTRDgeometry::GetSector(detector);
+ Int_t layere = AliTRDgeometry::GetLayer(detector);
+ // Some variables
+ Float_t b[2] = {0.0,0.0};
+ Float_t bCov[3] = {0.0,0.0,0.0};
+ if(esdtrack) esdtrack->GetImpactParameters(b,bCov);
+ if (bCov[0]<=0 || bCov[2]<=0) {
+ bCov[0]=0; bCov[2]=0;
+ }
+ Float_t dcaxy = b[0];
+ Float_t dcaz = b[1];
+ Int_t tpcnbclusters = 0;
+ if(esdtrack) tpcnbclusters = esdtrack->GetTPCclusters(0);
+ Double_t tpcsignal = 0.0;
+ if(esdtrack) tpcsignal = esdtrack->GetTPCsignal();
+ Int_t cutvdriftlinear = 0;
+ if(!pass) cutvdriftlinear = 1;
+
+ (* fDebugStreamer) << "FillCharge"<<
+ "detector="<<detector<<
+ "stack="<<stacke<<
+ "sm="<<sme<<
+ "layere="<<layere<<
+ "dcaxy="<<dcaxy<<
+ "dcaz="<<dcaz<<
+ "nbtpccls="<<tpcnbclusters<<
+ "tpcsignal="<<tpcsignal<<
+ "cutvdriftlinear="<<cutvdriftlinear<<
+ "ptrd="<<momentum<<
+ "nbtrdtracklet="<<numberoftrdtracklets<<
+ "charge="<<chargeQ<<
+ "\n";
+ }
+
+
} // if a good tracklet
}
// Loop clusters
//////////////////////////////////
- Short_t sigArr[AliTRDfeeParam::GetNcol()];
+ Float_t sigArr[AliTRDfeeParam::GetNcol()];
memset(sigArr, 0, AliTRDfeeParam::GetNcol()*sizeof(sigArr[0]));
Int_t ncl=0, tbf=0, tbl=0;
for(int ip=-1, jp=2; jp<5; ip++, jp++){
Int_t idx=col+ip;
if(idx<0 || idx>=AliTRDfeeParam::GetNcol()) continue;
- sigArr[idx]+=cl->GetSignals()[jp];
+ sigArr[idx]+=((Float_t)cl->GetSignals()[jp]);
}
if((fLimitChargeIntegration) && (!cl->IsInChamber())) continue;
TVectorD pars;
fLinearFitterTracklet->Eval();
fLinearFitterTracklet->GetParameters(pars);
- pointError = TMath::Sqrt(fLinearFitterTracklet->GetChisquare()/(nbli-2));
+ pointError = TMath::Sqrt(TMath::Abs(fLinearFitterTracklet->GetChisquare()/(nbli-2)));
errorpar = fLinearFitterTracklet->GetParError(1)*pointError;
dydt = pars[1];
//printf("chis %f, nbli %d, pointError %f, parError %f, errorpar %f\n",fLinearFitterTracklet->GetChisquare(),nbli,pointError,fLinearFitterTracklet->GetParError(1),errorpar);
///////////////////////////////////
Float_t signalSum(0.);
- Float_t mean = -1, rms = -1;
+ Float_t mean = 0.0, rms = 0.0;
Float_t dydx = tracklet->GetYref(1), tilt = tracklet->GetTilt(); // ,dzdx = tracklet->GetZref(1); (identical to the previous definition!)
Float_t dz = dzdx*(tbl-tbf)/10;
if(ncl>10){
signalSum+=sigArr[ip];
mean+=ip*sigArr[ip];
}
- mean/=signalSum;
+ if(signalSum > 0.0) mean/=signalSum;
for(Int_t ip = 0; ip<AliTRDfeeParam::GetNcol(); ip++)
rms+=sigArr[ip]*(ip-mean)*(ip-mean);
- rms = TMath::Sqrt(rms/signalSum);
+
+ if(signalSum > 0.0) rms = TMath::Sqrt(TMath::Abs(rms/signalSum));
rms -= TMath::Abs(dz*tilt);
dydx -= dzdx*tilt;
///////////////////////////////
- //if(fDebugLevel > 0){
+ if(fDebugLevel > 1){
if ( !fDebugStreamer ) {
//debug stream
TDirectory *backup = gDirectory;
"ncl="<<ncl<<
"\n";
- //}
+ }
/////////////////////////
// Cuts quality
// Debug stuff
////////////////////
- if(fDebugLevel > 0){
+ if(fDebugLevel > 1){
if ( !fDebugStreamer ) {
//debug stream
TDirectory *backup = gDirectory;
// Per tracklet: store or reset the info, fill the histos with the info
//////////////////////////////////////////////////////////////////////////////////////////
//_____________________________________________________________________________
-void AliTRDCalibraFillHisto::StoreInfoCHPHtrack(const AliTRDcluster *cl,const Double_t dqdl,const Int_t *group,const Int_t row,const Int_t col,const AliTRDcluster *cls)
+Float_t AliTRDCalibraFillHisto::StoreInfoCHPHtrack(const AliTRDcluster *cl,const Double_t dqdl,const Int_t *group,const Int_t row,const Int_t col,const AliTRDcluster *cls)
{
//
// Store the infos in fAmpTotal, fPHPlace and fPHValue
// Correct from the gain correction before
// cls is shared cluster if any
+ // Return the charge
+ //
//
//printf("StoreInfoCHPHtrack\n");
fPHPlace[time] = group[1];
fPHValue[time] = charge;
}
+
+ return correction;
}
//____________Offine tracking in the AliTRDtracker_____________________________
/////////////////////////////////////////////////////////
// Debug
////////////////////////////////////////////////////////
- if(fDebugLevel > 0){
+ if(fDebugLevel > 1){
if ( !fDebugStreamer ) {
//debug stream
TDirectory *backup = gDirectory;
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff