TRadCor.h

Go to the documentation of this file.

#ifndef RAD_TRadCor
#define RAD_TRadCor
 
#include "TVCrossPart.h"
#include "TKinemCut.h"
 
class TRadCor{
protected:
  unsigned int fNPart;   // number of cross-section parts
  unsigned int fNEvents; // number events to calculate integral
  double *fMax;        //[fNPart] array with maximum of matrix elements
  double *fMin;        //[fNPart] array with minimum of matrix elements must
               // be greater than 0
  double *fSigma;      //[fNPart] array with fMax error
  double *fCross;      //[fNPart] array with part cross-section
  double *fCross2;     //[fNPart] array with part cross-section
  double *fECross;      //[fNPart] array with part cross-section
  double *fECross2;     //[fNPart] array with part cross-section
  double *fSM;         //[fNPart] array with part sum of matrix elements
  double *fSM2;        //[fNPart] array with part sum of matrix elements
  double *fPart;       //[fNPart]
  unsigned int *fNSim; //[fNPart] array with number of selected events
  unsigned int *fN; //[fNPart] array with number of selected events
  unsigned int *fWN; //[fNPart] array with number of selected events
  bool *fInList; // [fNPart] array of selected parts of cross section
  TVCrossPart *fMatrix; // reference to matrix elements 
  bool fNew;
  double fTheta_Aver_0; // special case of Theta for elactic
            // cross-section to increase efficiency of
            // bhabha generation
  bool fIsSmear;  // to smear or not to smear that is a question
  void SetSpecialCase(unsigned int i){
    if ( i == 0 && !fIsSmear )
      fMatrix->SetThetaMin(fTheta_Aver_0);
    else
      fMatrix->SetThetaMin(gCut->ThetaMin());
  }
public:
  TRadCor();
  TRadCor(TVCrossPart*);
  ~TRadCor();
  
  void MakeMaximum(const unsigned int&);
  void MakeIntegral(const unsigned int&);
  void SetCuts(){};
  void SetAuxiliaryParameters(){}
  double Norm(const unsigned int &);
  void SetNEvents(const unsigned int &n){fNEvents = n;}
  void MakeAllMaximums();
  void MakeAllIntegrals();
  void Init();
  unsigned int MakeEvent();
  void MakePart();
  unsigned int GetNPart(){return fNPart;}
  double GetMax(unsigned int i){return fMax[i];}
  double GetMin(unsigned int i){return fMin[i];}
  unsigned int GetNSim(unsigned int i){return fNSim[i];}
  unsigned int GetN(unsigned int i){return fN[i];}
  unsigned int GetWN(unsigned int i){return fWN[i];}
  double GetCross(unsigned int i){return fCross[i];}
  double GetECross(unsigned int i){return fECross[i];}
  double GetCross2(unsigned int i){return fCross2[i];}
  double GetECross2(unsigned int i){return fECross2[i];}
  inline void GetEvent(double &p0, double &t0, double &ph0, 
               double &p1, double &t1, double &ph){
    fMatrix->GetEvent(p0, t0, ph0, p1, t1, ph);
  }
  inline void MakeEvent(double &p0, double &t0, double &ph0, 
            double &p1, double &t1, double &ph, unsigned int &ipart){
    ipart = MakeEvent();
    fMatrix->GetEvent(p0, t0, ph0, p1, t1, ph);
  }
  inline void MakeEvent(double *m, int &n){
    MakeEvent();
    fMatrix->GetEvent(m,n);    
  }
  inline TVCrossPart *GetCrossPart(){ return fMatrix;}
  void SetMType(const bool& x){fNew = x;}
  void SetSmear(const bool& x){fIsSmear =x;}
  void MakeCrossSection();
  void SetPartList(const bool *InList){
    for(unsigned int i=0; i<fNPart; i++) fInList[i] = InList[i];
  }
};
 
#endif //#ifndef RAD_TRadCor