|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
|
// Original Author: Gero Flucke
// last change : $Date: 2012/03/29 08:47:43 $
// by : $Author: flucke $
#ifndef COMPAREMILLEPEDE_H
#define COMPAREMILLEPEDE_H
#include <TString.h>
class PlotMillePede;
class GFHistManager;
class CompareMillePede
{
public:
CompareMillePede(const char *fileName1, const char *fileName2, Int_t iov1 = 1, Int_t iov2 = 1,
Int_t hieraLevel = 0);// iov1/2: which IOV ; hieraLev: -1 ignore, 0 lowest level, etc.
virtual ~CompareMillePede();
void DrawPedeParam(Option_t *option = "", unsigned int nNonRigidParam = 12);//"add": keep old canvas, "free1/2": if free param in file 1/2, "line": draw line in 2D plot
void DrawPedeParamVsLocation(Option_t *option = "", unsigned int nNonRigidParam = 12);//"add": keep old canvas, "free1/2": if free param in file 1/2
void DrawParam(Option_t *option="");//"add": keep old canvas, "free1/2": if free param in file 1/2
void DrawParamVsLocation(Option_t *option="");//"add": keep old canvas, "free1/2": if free param in file 1/2
void DrawParamDeltaMis(Option_t *option="");//"add": keep old canvas, "free1/2": if free param in file 1/2
void DrawParamDeltaMisVsLoc(Option_t *option="");//"add": keep old canvas, "free1/2": if free param in file 1/2
void DrawNumHits(Option_t *opt="");//"add": keep old canvas
void DrawAbsPos(Option_t *opt="");//"start": at start (else end), "add": keep old canvas
void DrawSurfaceDeformations(Option_t *option="", UInt_t firstPar=0, UInt_t lastPar=11,
const TString &whichOne = "result"
); //"add": keep old canvases, "limit": use GetMaxDev(), "noVs", noDiff"; "result"||"start"||"diff"
bool IsConsistent(); // check correct order of alignables, if false draw some hists
TString DeltaPar(UInt_t iPar) const; // par_2 - par_1
TString DeltaParBySigma(UInt_t iPar, const PlotMillePede *sigmaSource) const;
TString DeltaMisPar(UInt_t iPar) const; // abs(misalignment_2)-abs(misalignment_1)
TString DeltaMisParBySigma(UInt_t iPar, const PlotMillePede *sigmaSource) const;
TString DeltaPos(UInt_t iPos) const;
void AddIsFreeSel(TString &sel, const TString &option, UInt_t iPar) const;
void SetSubDetId(Int_t subDetId); // 1-6 are TPB, TPE, TIB, TID, TOB, TEC, -1 means: take all
void AddSubDetId(Int_t subDetId); // 1-6 are TPB, TPE, TIB, TID, TOB, TEC
void SetAlignableTypeId(Int_t alignableTypeId);//detunit=1,det=2,rod=3,etc. from AlignableObjectIdType (-1: all)
void SetHieraLevel(Int_t hieraLevel); // select hierarchical level (-1: all)
void AddAdditionalSel(const char *selection);// special select; StripDoubleOr1D,StripRphi,StripStereo
void AddAdditionalSel(const TString &xyzrPhiNhit, Float_t min, Float_t max); // x,y,z,r,phi,Nhit
// const TString GetAdditionalSel () const { return fAdditionalSel;}
void ClearAdditionalSel ();
void SetSurfDefDeltaBows(bool deltaBows); // take care: problems for false if drawing 1-sensor modules!
TString TitleAdd() const;
PlotMillePede* GetPlotMillePede1() {return fPlotMp1;}
PlotMillePede* GetPlotMillePede2() {return fPlotMp2;}
GFHistManager* GetHistManager() { return fHistManager;}
static const unsigned int kNpar; // number of parameters we have...
private:
CompareMillePede() : fPlotMp1(0), fPlotMp2(0), fHistManager(0) {}
Int_t PrepareAdd(bool addPlots);
PlotMillePede *fPlotMp1;
PlotMillePede *fPlotMp2;
GFHistManager *fHistManager;
};
#endif
|