TRadCor.h

Go to the documentation of this file.

#ifndef RAD_TRadCor
#define RAD_TRadCor
 
#include "TKinemCut.h"
#include "TVCrossPart.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