Abstract base class for particle identfication using functionals developed by L. Tassan-Got (IPN Orsay)
These functionals are defined in the KVTGIDFunctions namespace, and the functional used by each KVTGID object is defined by giving its name to the constructor. These objects are persistent, i.e. can be retrieved from a ROOT file and used for identification.
GetIdentification() method performs an identification. Status code can be retrieved afterwards using GetStatus();
A KVIDGrid identification grid can be generated from the functional in order to visualise the corresponding identification lines. (MakeIDGrid)
The following methods must be defined by child-classes:
SetIdent() -- determine how to set the identification
for each line in the grid
AddLine() -- define the type of IDLine added to the grid
NewGrid() -- define the type of IDGrid to generate
virtual KVIDLine * AddLine(KVIDGrid *)=0
virtual void SetIdent(KVIDLine *, Double_t ID)=0
Definition at line 44 of file KVTGID.h.
|
| | KVTGID () |
| | Default ctor. More...
|
| |
| | KVTGID (const Char_t *name, const Char_t *function, Double_t xmin, Double_t xmax, Int_t npar, Int_t par_x, Int_t par_y) |
| |
| | KVTGID (const KVTGID &) |
| | copy constructor More...
|
| |
| virtual | ~ KVTGID () |
| |
| void | AddIDTelescope (KVBase *tel) |
| |
| virtual void | AddLineToGrid (KVIDGrid *g, Int_t ID, Int_t npoints, Double_t xmin, Double_t xmax, Bool_t log_scale=kFALSE) |
| |
| void | ClearIDTelescopes () |
| |
| Int_t | Compare (const TObject *obj) const |
| |
| virtual void | Copy (TObject &tgid) const |
| | Copy this KVTGID function into the KVTGID object referenced by tgid. More...
|
| |
| Double_t | GetAlpha () const |
| |
| Double_t | GetBeta () const |
| |
| Double_t | GetDistanceToLine (Double_t x, Double_t y, Int_t id, Double_t *params=0) |
| |
| Double_t | GetEta () const |
| |
| Int_t | GetFunctionalType () const |
| |
| const Char_t * | GetFunctionName () const |
| |
| Double_t | GetG () const |
| |
| virtual Double_t | GetIdentification (Double_t ID_min, Double_t ID_max, Double_t &ID_quality, Double_t *par=0) |
| |
| Double_t | GetIDmax () const |
| |
| Double_t | GetIDmin () const |
| |
| const KVString & | GetIDTelescopes () const |
| |
| Double_t | GetLambda () const |
| |
| Int_t | GetLightEnergyDependence () const |
| |
| Int_t | GetMassFormula () const |
| |
| Double_t | GetMu () const |
| |
| Double_t | GetNu () const |
| |
| Double_t | GetPdx () const |
| |
| Double_t | GetPdy () const |
| |
| virtual Int_t | GetStatus () const |
| |
| virtual const Char_t * | GetStatusString () const |
| | Returns explanatory message for value of GetStatus() More...
|
| |
| const KVNumberList & | GetValidRuns () const |
| |
| const Char_t * | GetVarX () const |
| |
| const Char_t * | GetVarY () const |
| |
| Double_t | GetXi () const |
| |
| Int_t | GetZorA () const |
| |
| Bool_t | IsValidForRun (Int_t run) const |
| |
| Bool_t | IsValidForTelescope (KVBase *tel) const |
| |
| virtual void | Print (Option_t *option="") const |
| | Print info on functional and grid. More...
|
| |
| void | SetAlpha (Double_t val) |
| |
| void | SetBeta (Double_t val) |
| |
| void | SetEta (Double_t val) |
| |
| void | SetG (Double_t val) |
| |
| void | SetIDmax (Double_t x) |
| |
| void | SetIDmin (Double_t x) |
| |
| void | SetIDTelescopes (const TCollection *) |
| |
| void | SetLambda (Double_t val) |
| |
| void | SetLTGParameterNames () |
| |
| void | SetLTGParameters (Double_t *par) |
| |
| void | SetLTGParameters (Float_t *par) |
| |
| void | SetMassformula (Int_t val) |
| |
| void | SetMu (Double_t val) |
| |
| void | SetNu (Double_t val) |
| |
| void | SetPdx (Double_t val) |
| |
| void | SetPdy (Double_t val) |
| |
| void | SetValidRuns (const KVNumberList &r) |
| |
| void | SetVarX (const Char_t *x) |
| |
| void | SetVarY (const Char_t *x) |
| |
| void | SetXi (Double_t val) |
| |
| void | WriteToAsciiFile (std::ofstream &) const |
| | Write parameters of LTG fit in file. More...
|
| |
| KVTGID::KVTGID |
( |
const Char_t * |
name, |
|
|
const Char_t * |
function, |
|
|
Double_t |
xmin, |
|
|
Double_t |
xmax, |
|
|
Int_t |
npar, |
|
|
Int_t |
par_x, |
|
|
Int_t |
par_y |
|
) |
| |
Create TF1 named "name" using "function", range xmin->xmax and number of parameters npar.
The given function name must correspond to one of the functions defined in the KVTGIDFunctions namespace (without the "KVTGIDFunctions::" scope, i.e. function="tassangot_Z" is a valid argument for this constructor). The resulting object can be written to and read back from a ROOT file and the function pointer will be reinitialised automatically when the object is read back from the file.
par_x and par_y are the indices of the parameters of the functional which correspond to the 'X' and 'Y' coordinates of the identification map, respectively.
NOTE: the function pointer passed to the TF1 constructor is just a dummy to make sure that the right constructor is called, and that the TF1 is set up correctly. The actual function used is set using the name given as argument 'function'.
Definition at line 88 of file KVTGID.cpp.
| void KVTGID::AddLineToGrid |
( |
KVIDGrid * |
g, |
|
|
Int_t |
ID, |
|
|
Int_t |
npoints, |
|
|
Double_t |
xmin, |
|
|
Double_t |
xmax, |
|
|
Bool_t |
log_scale = kFALSE |
|
) |
| |
|
virtual |
Add a line to the grid 'g' for identification label 'ID' with 'npoints' points calculated between X-coordinates xmin and xmax. Points are omitted if the resulting value of the functional is not a number (TMath::IsNan = kTRUE).
If log_scale=kTRUE (default is kFALSE), more points are used at the beginning of the line than at the end, with a logarithmic dependence. add new line to grid
Reimplemented in KVTGIDZA.
Definition at line 241 of file KVTGID.cpp.
| Double_t KVTGID::GetDistanceToLine |
( |
Double_t |
x, |
|
|
Double_t |
y, |
|
|
Int_t |
id, |
|
|
Double_t * |
params = 0 |
|
) |
| |
Given a point (x,y) (which could, for example, be a point in a KVIDGrid line which we want to fit) we give the (vertical) distance to the functional identification line 'id'. If the point is below the line the distance is positive, if the point is above the line it is negative. If the parameter array 'params' is not given, we use the current values of the parameters.
Definition at line 364 of file KVTGID.cpp.
| Double_t KVTGID::GetIdentification |
( |
Double_t |
ID_min, |
|
|
Double_t |
ID_max, |
|
|
Double_t & |
ID_quality, |
|
|
Double_t * |
par = 0 |
|
) |
| |
|
virtual |
Use the functional with the current parameter values in order to perform an identification. The value returned is the estimated identification value. If the identification is not possible, -1 is returned (ID_quality=-1 also).
For status code use GetStatus().
ID_min and ID_max are the lower & upper limits for the identification. The functional must change sign between these two limits in order for the identification to be possible. ID_quality is the (absolute) value of the functional corresponding to the estimated identification: it is the distance from the identified point to the nearest identification line. Ideally it should be zero or as small as possible.
Optional argument "par" allows to replace the current parameters.
Definition at line 188 of file KVTGID.cpp.
| Int_t KVTGID::GetNumberOfLTGParameters |
( |
Int_t |
type, |
|
|
Int_t |
light |
|
) |
| |
|
static |
Static function returning number of parameters used by Tassan-Got functional. Depends on type of functional (type) and whether non-linear light response is calculated or not (light). According to the (type, light) combination the numbers and orders of parameters are : type=0 light=0 5 parameters: lambda, mu, g, pdx, pdy type=0 light<>0 6 parameters: lambda, mu, g, pdx, pdy, eta type<>0 light=0 9 parameters: lambda, alpha, beta, mu, nu, xi, g, pdx, pdy type<>0 light<>0 10 parameters: lambda, alpha, beta, mu, nu, xi, g, pdx, pdy, eta}
Definition at line 700 of file KVTGID.cpp.
| KVTGID * KVTGID::MakeTGID |
( |
const Char_t * |
name, |
|
|
Int_t |
type, |
|
|
Int_t |
light, |
|
|
Int_t |
ZorA, |
|
|
Int_t |
mass |
|
) |
| |
|
static |
Static function used to create TGID objects
type --—> type of functional
- type : =0->basic functional <>0->extended functional
- * For the basic formula :
- yy = ((g*E)**(mu+1)+lambda**(mu+1)*Z**2*A**mu)**(1/(mu+1))-g*E + pdy
- * For the extended formula :
- yy = ((g*E)**(mu+nu+1)+(lambda*Z**alpha*A**beta)**(mu+nu+1)+
- xi*A**mu*(g*E)**nu)**(1/(mu+mu+1))-g*E + pdy
light --—> treatment of CsI total light output
- light : =0->no non-linear light response <>0->non-linear light response included
- * If ih=0 no non-linear light response : E=xx-pdx
- * If ih<>0 non-linear light response included :
- E = sqrt(h**2+2*rho*h*(1+log(1+h/rho)))
- rho=eta*Z**2*A and h=xx-pdx
ZorA --—> functional used to find Z or A
Definition at line 429 of file KVTGID.cpp.
| void KVTGID::SetLTGParameters |
( |
Double_t * |
par | ) |
|
'par' is an array containing the LTG functional parameters, whose number and order depends on fType & fLight:
fType=0, fLight=0: Double_t par[5];
par[0] = lambda par[1] = mu par[2] = g par[3] = pdx par[4] = pdy
fType=0, fLight<>0: Double_t par[6];
par[0] = lambda par[1] = mu par[2] = g par[3] = pdx par[4] = pdy par[5] = eta
fType<>0, fLight=0: Double_t par[9];
par[0] = lambda par[1] = alpha par[2] = beta par[3] = mu par[4] = nu par[5] = xi par[6] = g par[7] = pdx par[8] = pdy
fType<>0, fLight<>0: Double_t par[10];
par[0] = lambda par[1] = alpha par[2] = beta par[3] = mu par[4] = nu par[5] = xi par[6] = g par[7] = pdx par[8] = pdy par[9] = eta
Definition at line 527 of file KVTGID.cpp.
| void KVTGID::SetLTGParameters |
( |
Float_t * |
par | ) |
|
'par' is an array containing the LTG functional parameters, whose number and order depends on fType & fLight:
fType=0, fLight=0: Double_t par[5];
par[0] = lambda par[1] = mu par[2] = g par[3] = pdx par[4] = pdy
fType=0, fLight<>0: Double_t par[6];
par[0] = lambda par[1] = mu par[2] = g par[3] = pdx par[4] = pdy par[5] = eta
fType<>0, fLight=0: Double_t par[9];
par[0] = lambda par[1] = alpha par[2] = beta par[3] = mu par[4] = nu par[5] = xi par[6] = g par[7] = pdx par[8] = pdy
fType<>0, fLight<>0: Double_t par[10];
par[0] = lambda par[1] = alpha par[2] = beta par[3] = mu par[4] = nu par[5] = xi par[6] = g par[7] = pdx par[8] = pdy par[9] = eta
Definition at line 630 of file KVTGID.cpp.
1.9.1