- 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 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;}
+ static Double_t SectorSizeX() {return NpadsX()*PadSizeX()/2;} //sector size x in cm
+ static Double_t SectorSizeY() {return NpadsY()*PadSizeY()/3;} //sector size y in cm
+ static Double_t PcSizeX() {return NpadsX()*PadSizeX()+DeadZone();} //PC size x, cm
+ static Double_t PcSizeY() {return NpadsY()*PadSizeY()+2*DeadZone();} //PC size y, cm
+
+ static Double_t Zfreon() {return 1.5;} //freon thinkness, cm
+ static Double_t Zwin() {return 0.5;} //radiator quartz window, cm
+ static Double_t Pc2Win() {return 8.0;} //cm between CsI PC and radiator quartz window
+ static Double_t Pc2Coll() {return 7.0;} //cm between CsI PC and third wire grid (collection wires)
+ static Double_t Pc2Anod() {return 0.204;} //cm between CsI PC and first wire grid (anod wires)
+ static Double_t Pc2Cath() {return 0.445;} //cm between CsI PC and second wire grid (cathode wires)
+ static Double_t Freon2Pc() {return Zfreon()+Zwin()+Pc2Win();} //cm between CsI PC and entrance to freon
+ static Double_t PitchAnod() {return PadSizeY()/2;} //cm between anode wires
+ static Double_t PitchCath() {return PadSizeY()/4;} //cm between cathode wires
+ static Double_t PitchColl() {return 0.5;} //cm between collection wires
+
+ static Double_t IonisationPotential() {return 26.0e-9;} //for CH4 in GeV taken from ????
+ static TVector2 MathiesonDelta() {return TVector2(5*0.18,5*0.18);} //area of 5 sigmas of Mathieson distribution (cm)
+ static Int_t MaxQdc() {return 4095;} //QDC number of channels
+
+ static Bool_t IsResolveClusters() {return fgIsResolveClusters;} //go after resolved clusters?
+ static Bool_t IsWireSag() {return fgIsWireSag;} //take wire sagita in account?
+ static Bool_t IsRadioSrc() {return fgIsRadioSrc;} //add radioactive source inside CH4?
+ static Int_t HV(Int_t sector) {
+ if (sector>=1 && sector <=6)
+ return fgHV[sector-1];
+ else {
+ ::Error("HV","Wrong sector %d",sector);
+ return kBad;
+ }
+ } //high voltage for this sector
+ static void SetDeclustering(Bool_t a) {fgIsResolveClusters=a;}
+ static void SetRadioSrc(Bool_t a) {fgIsRadioSrc=a;}
+ static void SetWireSag(Bool_t status) {fgIsWireSag=status;}
+ static void SetHV(Int_t sector,Int_t hv){fgHV[sector-1]=hv;}
+ static void SetAngleRot(Double_t rot) {fgAngleRot =rot;}
+ static Double_t IndOfRefC6F14(Double_t eV) {return eV*0.0172+1.177;} // eV = photon energy in eV
+ static Double_t IndOfRefSiO2(Double_t eV) {Double_t e1=10.666,e2=18.125,f1=46.411,f2= 228.71;
+ return TMath::Sqrt(1.+f1/(e1*e1-TMath::Power(eV,2))+f2/(e2*e2-TMath::Power(eV,2)));}//TDR p.35
+ static Double_t IndOfRefCH4() {return 1.000444;}