- Double_t MeanIdxRad ()const {return 1.29204;} //<--TEMPORAR--> to be removed in future. Mean ref index C6F14
- Double_t MeanIdxWin ()const {return 1.57819;} //<--TEMPORAR--> to be removed in future. Mean ref index quartz
- Float_t DistCut ()const {return 1.0;} //<--TEMPORAR--> to be removed in future. Cut for MIP-TRACK residual
- Float_t QCut ()const {return 100;} //<--TEMPORAR--> to be removed in future. Separation PHOTON-MIP charge
- Float_t MultCut ()const {return 200;} //<--TEMPORAR--> to be removed in future. Multiplicity cut to activate WEIGHT procedure
+ //For optical properties
+ static Double_t EPhotMin() {return 5.5;} //
+ static Double_t EPhotMax() {return 8.5;} //Photon energy range,[eV]
+ static Double_t NIdxRad(Double_t eV,Double_t temp) {return TMath::Sqrt(1+0.554*(1239.84/eV)*(1239.84/eV)/((1239.84/eV)*(1239.84/eV)-5769)-0.0005*(temp-20));}
+ static Double_t NIdxWin(Double_t eV) {return TMath::Sqrt(1+46.411/(10.666*10.666-eV*eV)+228.71/(18.125*18.125-eV*eV));}
+ static Double_t NMgF2Idx(Double_t eV) {return 1.7744 - 2.866e-3*(1239.842609/eV) + 5.5564e-6*(1239.842609/eV)*(1239.842609/eV);} // MgF2 idx of trasparency system
+ static Double_t NIdxGap(Double_t eV) {return 1+0.12489e-6/(2.62e-4 - eV*eV/1239.84/1239.84);}
+ static Double_t LAbsRad(Double_t eV) {return (eV<7.8)*(GausPar(eV,3.20491e16,-0.00917890,0.742402)+GausPar(eV,3035.37,4.81171,0.626309))+(eV>=7.8)*0.0001;}
+ static Double_t LAbsWin(Double_t eV) {return (eV<8.2)*(818.8638-301.0436*eV+36.89642*eV*eV-1.507555*eV*eV*eV)+(eV>=8.2)*0.0001;}//fit from DiMauro data 28.10.03
+ static Double_t LAbsGap(Double_t eV) {return (eV<7.75)*6512.399+(eV>=7.75)*3.90743e-2/(-1.655279e-1+6.307392e-2*eV-8.011441e-3*eV*eV+3.392126e-4*eV*eV*eV);}
+ static Double_t QEffCSI(Double_t eV) {return (eV>6.07267)*0.344811*(1-exp(-1.29730*(eV-6.07267)));}//fit from DiMauro data 28.10.03
+ static Double_t GausPar(Double_t x,Double_t a1,Double_t a2,Double_t a3) {return a1*TMath::Exp(-0.5*((x-a2)/a3)*((x-a2)/a3));}
+ inline static Double_t FindTemp(Double_t tLow,Double_t tUp,Double_t y); //find the temperature of the C6F14 in a given point with coord. y (in x is uniform)
+
+
+ Double_t GetEPhotMean ()const {return fPhotEMean;}
+ Double_t GetRefIdx ()const {return fRefIdx;} //running refractive index
+
+ Double_t MeanIdxRad ()const {return NIdxRad(fPhotEMean,fTemp);}
+ Double_t MeanIdxWin ()const {return NIdxWin(fPhotEMean);}
+ //
+ Float_t DistCut ()const {return 1.0;} //<--TEMPORAR--> to be removed in future. Cut for MIP-TRACK residual
+ Float_t QCut ()const {return 100;} //<--TEMPORAR--> to be removed in future. Separation PHOTON-MIP charge
+ Float_t MultCut ()const {return 200;} //<--TEMPORAR--> to be removed in future. Multiplicity cut to activate WEIGHT procedure