// --------------------------------------------------------------------- // ----- TRootKRATEvent header file ----- // ----- Created 10/11/13 by J. Lukasik ----- // ----- P. Pawlowski, S. Kupny ----- // --------------------------------------------------------------------- #ifndef __TROOT_KRATTA_CALIB_CLASS__ #define __TROOT_KRATTA_CALIB_CLASS__ #include #include #include #include #include "Rtypes.h" #include "TClass.h" #include "TClonesArray.h" #include "TCutG.h" #include "TFile.h" #include "TObject.h" #include "TRandom3.h" #include "TString.h" //______________________________________________________________________ class TRootKRATParticle : public TNamed{ public: /// Swiatowid members: Int_t Module; // module No. [0..34] int Code; // Quality code int Spect; // Number of the identification map float Xident; // X value on the identification map float Yident; // Y value on the identification map int PID; // Particle identifier int Z; // Charge int A; // Mass number float RealZ; // Real charge float RealA; // Real mass number float Energy; // in the middle of the target from the sum of DE_i // new calibration using "anti-quenching" for Z=1 float EnergyOld; // Parlog for all /// DeltaE in consecutive layers: float de_csi2; // float de_myl1; // float de_csi1; //from R-E tables float de_csi1_cal; //from calibration float de_myl0; // float de_sil1; //from R-E tables float de_sil1_cal; //from calibration float de_ddl1; // float de_sil0; //from R-E tables float de_sil0_cal; //from calibration !!!!!!!!!!!!!!!!!!!!!!!!!!!! float de_ddl0; // float de_win; // float de_air; // float de_tgt; // float WeightCol; // efficiency correction // from nuclear collision length CsI 22.30 cm (Sigma_TOT) float WeightInt; // efficiency correction // from nuclear interaction length CsI 38.04 cm // (Sigma_TOT - Sigma_elastic - Sigma_quasi_elastic) float WeightGeant; // efficiency correction from GEANT float RangeCsI; // range in CsI [cm] float DeltaECsI; // energy loss in CsI [MeV] float EnergyDE01; // from DeltaE in PD0+PD1 float EnergyDE0; // from DeltaE in PD0 float EneResDE01; // residual energy from DeltaE in PD0+PD1 float EneResDE0; // residual energy from DeltaE in PD0 float Theta; // polar angle of the module center float Phi; // azimuthal angle of the module center float ThetaRndm; // random angle (uniform within active area) float PhiRndm; // random angle (uniform within active area) bool Pextra; // proton punching through CSI1 w/o signal in CSI2 bool Gamma; // gamma bool Gamma2[2]; // gamma in [CSI1, CSI2] bool Star; // in "star" cut bool Star2[2]; // in "star" cut [CSI1, CSI2] bool Stopped[4]; // stoped in [PD0, PD1, CSI1, CSI2] bool PunchThrough; // punch through CSI2 or background bool InTime; // in time window bool InTime3[3]; // in time window [PD0, PD1, PD2] float DiffTimeSigma[3]; // time- in the units of sigma [PD0, PD1, PD2] float RapidityLab; // float RapidityCM; // float Pt; // float Px; // float Py; // float Pz; // float Xt; // float Winv; // float Gam; // gamma=1/sqrt(1-b*b) /// Raw channels float ch_csi2; // CsI2 total light float ch_csi1; // CsI1 total light float ch_sil1; // Si1 amplitude float ch_sil0; // Si0 amplitude float ch_slo2; // Slow for CsI2 float ch_fas2; // Fast for CsI2 float ch_sof2; // Ratio Slow/Fast for CsI2 float ch_mod0; // Mode for Si0 float ch_sum1; // CsI1 sum float ch_dif1; // CsI1 difference float ch_sum2; // CsI2 sum float ch_dif2; // CsI2 difference float ch_log01; // log(Si0+Si1) float ch_slo1; // Slow for CsI1 float ch_fas1; // Fast for CsI1 float ch_sof1; // Ratio Slow/Fast for CsI1 float ch_log1; // log(ch_csi1) float ch_log2; // log(ch_csi2) float ch_tm0; // time0 PD0 float ch_tm1; // time0 PD1 float ch_tm2; // time0 PD2 int pthr_flag; // non-zero if outside punch through channel /// Methods: TRootKRATParticle(); virtual ~TRootKRATParticle(); void Clear(Option_t *option=""); void Copy(const TRootKRATParticle& ds); TRootKRATParticle& operator=(const TRootKRATParticle& ds); virtual void ls(Option_t* opt = 0)const; friend std::ostream& operator<<(std::ostream&, const TRootKRATParticle& ae); ClassDef(TRootKRATParticle,1) }; //______________________________________________________________________ class TRootKRATEvent : public TNamed{ public: TClonesArray *fHit; Int_t fNhit; //multiplicity of hits (incremented automatically) Int_t Run; // run No. Int_t Evt; // event No. Int_t TgtThck; // target thickness [um] char TgtMat[20]; // target material ///--> tgt|air|win|ddl0|sil0|ddl1|sil1|myl0|csi1|myl1|csi2|myl2|sil2|ddl2 char Beam[20]; // beam float BeamEnergy; // beam kinetic energy in the middle of the target Int_t BeamA; // beam A Int_t BeamZ; // beam Z float BeamRapidity;// beam rapidity float CMRapidity; // CM rapidity TRootKRATEvent(); virtual ~TRootKRATEvent(); ///TRootKRATParticle* AddEmptyHit() {return new((*fHit)[fNhit++]) TRootKRATParticle();} // v1 TRootKRATParticle* GetPtrToAnotherOneHit(); void Clear(Option_t* opt=""); TRootKRATEvent& operator=(const TRootKRATEvent& kc); friend std::ostream& operator<<(std::ostream&, const TRootKRATEvent& kc); ClassDef(TRootKRATEvent,1) }; #endif ///__TROOT_KRATTA_CALIB_CLASS__