class AliTPCFCVoltError3D: public AliTPCCorrection
AliTPCFCVoltError3D class
The class calculates the space point distortions due to residual voltage errors on the Field Cage (FC) boundaries (rods and strips) of the TPC in 3D. It uses the Poisson relaxation technique in three dimension as implemented in the parent class.Although the calculation is performed in 3D, the calculation time was significantly reduced by using certain symmetry conditions within the calculation.
The input parameters can be set via the functions (e.g.) SetRodVoltShift(rod,dV) where
rod is the number of the rod on which the voltage offset dV is set. The corresponding
shift in z direction would be $dz=dV/40$ with an opposite sign for the C side. The
rods are numbered in anti-clock-wise direction starting at $\phi=0$. Rods in the IFC
are from 0 to 17, rods on the OFC go from 18 to 35.
This convention is following the offline numbering scheme of the ROCs.
Different misalignment scenarios can be modeled:
- A rotated clip scenario is only possible at two positions (Rod 11 at IFC, rod 3(+18) at OFC) and can be set via SetRotatedClipVolt. The key feature is that at the mentioned positions, the strip-ends were combined. At these positions, a systematic offset of one strip-end in respect to the other is possible.
- A normal rod offset, where the strips follow the movement of the rod, can be set via SetRodVoltShift. It has a anti-mirrored signature in terms of distortions when compared to the rotated clip. This misalignment is possible at each single rod of the FC.
- A simple rod offset, where the strips do not follow the shift, results in an even more localized distortion close to the rod. The difference to a rod shift, where the strips follow, is more dominant on the OFC. This effect can be set via the function SetCopperRodShift.
{
gROOT->SetStyle("Plain"); gStyle->SetPalette(1);
TCanvas *c2 = new TCanvas("cAliTPCVoltError3D","cAliTPCVoltError3D",500,450);
AliTPCFCVoltError3D fc;
fc.SetOmegaTauT1T2(0,1,1);
fc.SetRotatedClipVoltA(0,40);
fc.SetRodVoltShiftA(3,40);
fc.SetCopperRodShiftA(7+18,40);
fc.SetRodVoltShiftA(15+18,40);
fc.CreateHistoDRPhiinXY(10)->Draw("cont4z");
return c2;
}
Date: 08/08/2010
Authors: Jim Thomas, Stefan Rossegger
Function Members (Methods)
public:
protected:
private:
| AliTPCFCVoltError3D(const AliTPCFCVoltError3D&) | |
| AliTPCFCVoltError3D& | operator=(const AliTPCFCVoltError3D&) |
Data Members
public:
| enum { | kRows | |
| kColumns | ||
| kPhiSlicesPerSector | ||
| kPhiSlices | ||
| kIterations | ||
| }; | ||
| enum AliTPCCorrection::CompositionType { | kParallel | |
| kQueue | ||
| }; | ||
| enum AliTPCCorrection::[unnamed] { | kNR | |
| kNPhi | ||
| kNZ | ||
| }; | ||
| enum TObject::EStatusBits { | kCanDelete | |
| kMustCleanup | ||
| kObjInCanvas | ||
| kIsReferenced | ||
| kHasUUID | ||
| kCannotPick | ||
| kNoContextMenu | ||
| kInvalidObject | ||
| }; | ||
| enum TObject::[unnamed] { | kIsOnHeap | |
| kNotDeleted | ||
| kZombie | ||
| kBitMask | ||
| kSingleKey | ||
| kOverwrite | ||
| kWriteDelete | ||
| }; |
protected:
| Int_t | AliTPCCorrection::fILow | |
| Bool_t | AliTPCCorrection::fIsLocal | switch to indicate that the distortion is a local vector drphi/dz, dr/dz |
| Int_t | AliTPCCorrection::fJLow | |
| Int_t | AliTPCCorrection::fKLow | variable to help in the interpolation |
| TString | TNamed::fName | object identifier |
| Double_t | AliTPCCorrection::fT1 | tensor term of wt - T1 |
| Double_t | AliTPCCorrection::fT2 | tensor term of wt - T2 |
| TString | TNamed::fTitle | object title |
| static TObjArray* | AliTPCCorrection::fgVisualCorrection | array of orrection for visualization |
| static const Double_t | AliTPCCorrection::fgkCathodeV | Cathode Voltage (volts) |
| static const Double_t | AliTPCCorrection::fgkEM | charge/mass in [C/kg] |
| static const Double_t | AliTPCCorrection::fgkGG | Gating Grid voltage (volts) |
| static const Double_t | AliTPCCorrection::fgkIFCRadius | Mean Radius of the Inner Field Cage ( 82.43 min, 83.70 max) (cm) |
| static const Double_t | AliTPCCorrection::fgkOFCRadius | Mean Radius of the Outer Field Cage (252.55 min, 256.45 max) (cm) |
| Double_t | AliTPCCorrection::fgkPhiList[181] | points in the phi direction (for the lookup table) |
| Double_t | AliTPCCorrection::fgkRList[72] | points in the radial direction (for the lookup table) |
| static const Double_t | AliTPCCorrection::fgkTPCZ0 | nominal gating grid position |
| Double_t | AliTPCCorrection::fgkZList[166] | points in the z direction (for the lookup table) |
| static const Double_t | AliTPCCorrection::fgkZOffSet | Offset from CE: calculate all distortions closer to CE as if at this point |
| static const Double_t | AliTPCCorrection::fgkdvdE | [cm/V] drift velocity dependency on the E field (from Magboltz for NeCO2N2 at standard environment) |
| static const Double_t | AliTPCCorrection::fgke0 | vacuum permittivity [A·s/(V·m)] |
private:
| Float_t | fC0 | coefficient C0 (compare Jim Thomas's notes for definitions) |
| Float_t | fC1 | coefficient C1 (compare Jim Thomas's notes for definitions) |
| Float_t | fCopperRodShiftA[36] | only Rod shift |
| Float_t | fCopperRodShiftC[36] | only Rod shift |
| Bool_t | fInitLookUp | flag to check if the Look Up table was created (SUM) |
| Bool_t | fInitLookUpBasic[6] | ! flag if the basic lookup was created (shifted Rod (IFC,OFC) or rotated clip (IFC,OFC)) |
| TMatrixD* | fLookUpBasic1DeltaEz[31] | ! Array to store electric field integral (int Ez) |
| TMatrixD* | fLookUpBasic1EphiOverEz[31] | ! Array to store electric field integral (int Ephi/Ez) |
| TMatrixD* | fLookUpBasic1ErOverEz[31] | ! Array to store electric field integral (int Er/Ez) |
| TMatrixD* | fLookUpBasic2DeltaEz[31] | ! Array to store electric field integral |
| TMatrixD* | fLookUpBasic2EphiOverEz[31] | ! Array to store electric field integral |
| TMatrixD* | fLookUpBasic2ErOverEz[31] | ! Array to store electric field integral |
| TMatrixD* | fLookUpBasic3DeltaEz[31] | ! Array to store electric field integral |
| TMatrixD* | fLookUpBasic3EphiOverEz[31] | ! Array to store electric field integral |
| TMatrixD* | fLookUpBasic3ErOverEz[31] | ! Array to store electric field integral |
| TMatrixD* | fLookUpBasic4DeltaEz[31] | ! Array to store electric field integral |
| TMatrixD* | fLookUpBasic4EphiOverEz[31] | ! Array to store electric field integral |
| TMatrixD* | fLookUpBasic4ErOverEz[31] | ! Array to store electric field integral |
| TMatrixD* | fLookUpBasic5DeltaEz[31] | ! Array to store electric field integral |
| TMatrixD* | fLookUpBasic5EphiOverEz[31] | ! Array to store electric field integral |
| TMatrixD* | fLookUpBasic5ErOverEz[31] | ! Array to store electric field integral |
| TMatrixD* | fLookUpBasic6DeltaEz[31] | ! Array to store electric field integral |
| TMatrixD* | fLookUpBasic6EphiOverEz[31] | ! Array to store electric field integral |
| TMatrixD* | fLookUpBasic6ErOverEz[31] | ! Array to store electric field integral |
| TMatrixF* | fLookUpDeltaEz[181] | Array to store electric field integral (int Er/Ez) |
| TMatrixF* | fLookUpEphiOverEz[181] | Array to store electric field integral (int Er/Ez) |
| TMatrixF* | fLookUpErOverEz[181] | Array to store electric field integral (int Er/Ez) |
| Float_t | fRodVoltShiftA[36] | Rod (&strips) shift in Volt (40V~1mm) |
| Float_t | fRodVoltShiftC[36] | Rod (&strips) shift in Volt (40V~1mm) |
| Float_t | fRotatedClipVoltA[2] | rotated clips at HV rod |
| Float_t | fRotatedClipVoltC[2] | rotated clips at HV rod |
Class Charts
| Inheritance Chart: | ||||||||||||
|
Function documentation
Bool_t AddCorrectionCompact(AliTPCCorrection* corr, Double_t weight)
Add correction and make them compact Assumptions: - origin of distortion/correction are additive - only correction ot the same type supported ()
void GetCorrection(const Float_t* x, const Short_t roc, Float_t* dx)
Calculates the correction due e.g. residual voltage errors on the TPC boundaries
void InitFCVoltError3D()
Initialization of the Lookup table which contains the solutions of the Dirichlet boundary problem Calculation of the single 3D-Poisson solver is done just if needed (see basic lookup tables in header file)
void Print(Option_t* option = "") const
Print function to check the settings of the Rod shifts and the rotated clips option=="a" prints the C0 and C1 coefficents for calibration purposes
void SetOmegaTauT1T2(Float_t omegaTau, Float_t t1, Float_t t2)
common setters and getters for tangled ExB effect
void SetRodVoltShiftA(Int_t rod, Float_t voltOffset, Bool_t doInit = kTRUE)
setters and getters Set rod shift in Voltage equivalents (40V ~ 1mm) rod numbers: 0-17 (IFC), 18-35 (OFC) note: strips move accordingly
{fRodVoltShiftA[rod]=voltOffset; fInitLookUp=doInit;}void SetRodVoltShiftC(Int_t rod, Float_t voltOffset, Bool_t doInit = kTRUE)
{fRodVoltShiftC[rod]=voltOffset; fInitLookUp=doInit;}void SetRotatedClipVoltA(Int_t rod, Float_t voltOffset, Bool_t doInit = kTRUE)
Set rotated clip (just at High Voltage RODs) in Voltage equivalents (40V ~ 1mm) rod number: 0 (IFC), 1 (OFC)
{fRotatedClipVoltA[rod]=voltOffset; fInitLookUp=doInit;}void SetRotatedClipVoltC(Int_t rod, Float_t voltOffset, Bool_t doInit = kTRUE)
{fRotatedClipVoltC[rod]=voltOffset; fInitLookUp=doInit;}void SetCopperRodShiftA(Int_t rod, Float_t voltOffset, Bool_t doInit = kTRUE)
Set rod shift in Voltage equivalents (40V ~ 1mm) rod numbers: 0-17 (OFC) note: strips DO NOT move, only the copper rods do ...
{fCopperRodShiftA[rod]=voltOffset; fInitLookUp=doInit;}void SetCopperRodShiftC(Int_t rod, Float_t voltOffset, Bool_t doInit = kTRUE)
{fCopperRodShiftC[rod]=voltOffset; fInitLookUp=doInit;}AliTPCFCVoltError3D & operator=(const AliTPCFCVoltError3D& )