PndGemMagneticFieldVsTrackParameters.h

Go to the documentation of this file.

//* $Id: */

// -------------------------------------------------------------------------
// -----                    PndGemMagneticFieldVsTrackParameters header file                 -----
// -----                  Created 12/02/2009 by R. Karabowicz          -----
// -------------------------------------------------------------------------


#ifndef PNDGEMMAGNETICFIELDVSTRACKPARAMETERS_H
#define PNDGEMMAGNETICFIELDVSTRACKPARAMETERS_H 1

#include "FairTask.h"

#include "TStopwatch.h"

#include "TH1F.h"
#include "TH2F.h"
#include "TH3F.h"
#include "TF1.h"
#include "TF2.h"

#include <list>
#include <map>

class TClonesArray;
class PndGemDigiPar;



class PndGemMagneticFieldVsTrackParameters : public FairTask
{

 public:

  PndGemMagneticFieldVsTrackParameters();


  PndGemMagneticFieldVsTrackParameters(Int_t iVerbose);


  PndGemMagneticFieldVsTrackParameters(const char* name, Int_t iVerbose);


  virtual ~PndGemMagneticFieldVsTrackParameters();


  virtual void Exec(Option_t* opt);



 private:

  //  FairGeoStsPar*     fGeoPar;       /** Geometry parameter container **/
  //  PndGemDigiPar*    fDigiPar;      /** Digitisation parameter container **/
  //  PndGemDigiScheme* fDigiScheme;   /** Digitisation scheme **/
  PndGemDigiPar* fDigiPar;
  
  Double_t fTrackFinderOnHits_ParThetaA;
  Double_t fTrackFinderOnHits_ParThetaB;
  
  Double_t fTrackFinderOnHits_ParTheta0;
  Double_t fTrackFinderOnHits_ParTheta1;
  Double_t fTrackFinderOnHits_ParTheta2;
  Double_t fTrackFinderOnHits_ParTheta3;
  
  Double_t fTrackFinderOnHits_ParRadPhi0;
  Double_t fTrackFinderOnHits_ParRadPhi2;
  
  Double_t fTrackFinderOnHits_ParMat0[3];
  Double_t fTrackFinderOnHits_ParMat1[3];

  TClonesArray* fMCTrackArray;
  TClonesArray* fMCPointArray;

  TStopwatch        fTimer;

  Int_t fNofEvents;  // event counter


  TList* fHistoList;
  // theta angle vs radius vs momentum
  TH3F* fhThetaVsRadiusVsMomentumAll;
  TH3F* fhThetaVsRadiusVsMomentum[10];
  TF2*  ffThetaVsRadiusVsMomentum[10];
  // radius vs angle
  TH2F*  fhRadiusVsAngle;
  TH2F*  fhRadiusVsAnglePair[20];
  TF1*   ffRadiusVsAngle[3];
  // momentum vs phi angle difference
  TH3F* fhMomentumVsPhiDiffVsRadius[10];
  TH2F* fhMomentumVsPhiDiffAll[10];
  TH2F* fhMomentumVsPhiDiff   [10][20];
  TF2*  ffMomentumVsPhiDiffVsRadius[10];
  TF1*  ffMomentumVsPhiDiffAll[10];
  TF1*  ffMomentumVsPhiDiff   [10][20];
  // transverse momentum vs trajectory helix radius
  TH2F* fhMomTransVsHRadius[10];
  // track phi versus hits phis
  TH3F* fhTrackPhiVsHitPhis   [10];
  TH2F* fhTrackPhiVsCalcPhi   [10];

  Int_t Fill1StationHistograms();
  Int_t Fill2StationsHistograms();

  void CreateHistos();

  virtual void SetParContainers();

  virtual void Finish();

  virtual InitStatus Init();


  virtual InitStatus ReInit();


  void Reset();



  ClassDef(PndGemMagneticFieldVsTrackParameters,1);

};

#endif