PndGemMagneticFieldVsTrackParameters Class Reference

#include <PndGemMagneticFieldVsTrackParameters.h>

Inheritance diagram for PndGemMagneticFieldVsTrackParameters:

Public Member Functions
	PndGemMagneticFieldVsTrackParameters ()
	PndGemMagneticFieldVsTrackParameters (Int_t iVerbose)
	PndGemMagneticFieldVsTrackParameters (const char *name, Int_t iVerbose)
virtual	~PndGemMagneticFieldVsTrackParameters ()
virtual void	Exec (Option_t *opt)

Private Member Functions
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)

Private Attributes
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
TList *	fHistoList
TH3F *	fhThetaVsRadiusVsMomentumAll
TH3F *	fhThetaVsRadiusVsMomentum [10]
TF2 *	ffThetaVsRadiusVsMomentum [10]
TH2F *	fhRadiusVsAngle
TH2F *	fhRadiusVsAnglePair [20]
TF1 *	ffRadiusVsAngle [3]
TH3F *	fhMomentumVsPhiDiffVsRadius [10]
TH2F *	fhMomentumVsPhiDiffAll [10]
TH2F *	fhMomentumVsPhiDiff [10][20]
TF2 *	ffMomentumVsPhiDiffVsRadius [10]
TF1 *	ffMomentumVsPhiDiffAll [10]
TF1 *	ffMomentumVsPhiDiff [10][20]
TH2F *	fhMomTransVsHRadius [10]
TH3F *	fhTrackPhiVsHitPhis [10]
TH2F *	fhTrackPhiVsCalcPhi [10]

Detailed Description

Definition at line 41 of file PndGemMagneticFieldVsTrackParameters.h.

Constructor & Destructor Documentation

PndGemMagneticFieldVsTrackParameters::PndGemMagneticFieldVsTrackParameters

(

)

Default constructor

Definition at line 45 of file PndGemMagneticFieldVsTrackParameters.cxx.

References fDigiPar, fMCPointArray, fMCTrackArray, and Reset().

                                                                           : FairTask("GEM MagneticFieldVsTrackParametersr", 1) {
  fDigiPar     = NULL;
  fMCPointArray      = NULL;
  fMCTrackArray       = NULL;
  Reset();
}

PndGemMagneticFieldVsTrackParameters::PndGemMagneticFieldVsTrackParameters

(

Int_t

iVerbose

)

Standard constructor

Definition at line 56 of file PndGemMagneticFieldVsTrackParameters.cxx.

References fDigiPar, fMCPointArray, fMCTrackArray, and Reset().

  : FairTask("GEM MagneticFieldVsTrackParametersr", iVerbose) { 
  fDigiPar     = NULL;
  fMCPointArray      = NULL;
  fMCTrackArray       = NULL;
  Reset();
}

PndGemMagneticFieldVsTrackParameters::PndGemMagneticFieldVsTrackParameters	(	const char *	name,
		Int_t	iVerbose
	)

Constructor with name

Definition at line 68 of file PndGemMagneticFieldVsTrackParameters.cxx.

References fDigiPar, fMCPointArray, fMCTrackArray, and Reset().

  : FairTask(name, iVerbose) { 
  fDigiPar     = NULL;
  fMCPointArray      = NULL;
  fMCTrackArray       = NULL;
  Reset();
}

PndGemMagneticFieldVsTrackParameters::~PndGemMagneticFieldVsTrackParameters ( )

virtual

Destructor

Definition at line 80 of file PndGemMagneticFieldVsTrackParameters.cxx.

References fDigiPar, and Reset().

                                                                            { 
  if ( fDigiPar)   delete fDigiPar;
  Reset();
}

Member Function Documentation

PndGemMagneticFieldVsTrackParameters::ClassDef ( PndGemMagneticFieldVsTrackParameters  , 1    )

private

void PndGemMagneticFieldVsTrackParameters::CreateHistos ( )

private

Definition at line 298 of file PndGemMagneticFieldVsTrackParameters.cxx.

References Double_t, fDigiPar, ffMomentumVsPhiDiff, ffMomentumVsPhiDiffAll, ffMomentumVsPhiDiffVsRadius, ffRadiusVsAngle, ffThetaVsRadiusVsMomentum, fHistoList, fhMomentumVsPhiDiff, fhMomentumVsPhiDiffAll, fhMomentumVsPhiDiffVsRadius, fhMomTransVsHRadius, fhRadiusVsAngle, fhRadiusVsAnglePair, fhThetaVsRadiusVsMomentum, fhThetaVsRadiusVsMomentumAll, fhTrackPhiVsCalcPhi, fhTrackPhiVsHitPhis, and PndGemDigiPar::GetNStations().

Referenced by Init().

                                                        {
  fHistoList = new TList();
  // number of mc tracks, reco tracks, efficiency as function of MOMENTUM

  const Int_t nofMomBins = 19;
  Double_t momBins[nofMomBins+1] = {0.05,0.15,0.25,0.35,0.45,
                                    0.55,0.65,0.75,0.85,0.95,
                                    1.05,1.15,1.25,1.75,2.25,
                                    2.75,3.25,4.75,7.25,10.75};
  const Int_t nofRadiusBins = 60; 
  Double_t minRadius =  0.;
  Double_t maxRadius = 60.; 
  Double_t minRRadius =  0.; // radius/zdist
  Double_t maxRRadius =  0.4; // radius/zdist
  Double_t radBins[nofRadiusBins+1];
  for ( Int_t itemp = 0 ; itemp < nofRadiusBins+1 ; itemp++ ) radBins[itemp] = minRadius + (Double_t)itemp*(maxRadius-minRadius)/(Double_t(nofRadiusBins));
  Double_t rRadBins[nofRadiusBins+1];
  for ( Int_t itemp = 0 ; itemp < nofRadiusBins+1 ; itemp++ ) rRadBins[itemp] = minRRadius + (Double_t)itemp*(maxRRadius-minRRadius)/(Double_t(nofRadiusBins));
  const Int_t nofThetaBins = 60;
  Double_t minTheta =  0.;
  Double_t maxTheta = 30.;
  Double_t theBins[nofThetaBins+1];
  for ( Int_t itemp = 0 ; itemp < nofThetaBins+1 ; itemp++ ) theBins[itemp] = minTheta + (Double_t)itemp*(maxTheta-minTheta)/(Double_t(nofThetaBins));
  const Int_t nofPhiBins = 160;
  Double_t minPhi = 0.;
  Double_t maxPhi = 40.;
  Double_t phiBins[nofPhiBins+1];
  for ( Int_t itemp = 0 ; itemp < nofPhiBins+1 ; itemp++ ) phiBins[itemp] = minPhi + (Double_t)itemp*(maxPhi-minPhi)/(Double_t(nofPhiBins));
  const Int_t nofDetPhiBins = 400;
  const Int_t nofDetMomBins = 100;
  Double_t minMom =  0.;
  Double_t maxMom = 10.;
  const Int_t nofPhi1Bins = 360;
  Double_t phi1Bins[nofPhi1Bins+1];
  for ( Int_t itemp = 0 ; itemp < nofPhi1Bins+1 ; itemp++ ) phi1Bins[itemp] = itemp;
  const Int_t nofPhi2Bins = 120;
  Double_t minPhi2 = -60.;
  Double_t maxPhi2 =  60.;
  Double_t phi2Bins[nofPhi2Bins+1];
  for ( Int_t itemp = 0 ; itemp < nofPhi2Bins+1 ; itemp++ ) phi2Bins[itemp] = minPhi2 + (Double_t)itemp*(maxPhi2-minPhi2)/(Double_t(nofPhi2Bins));

  fhThetaVsRadiusVsMomentumAll = new TH3F("fhThetaVsRadiusVsMomentumAll",
                                          "Track Theta vs Point Radius/Z vs Track Momentum, all;p [GeV/c];r [cm];#theta [#circ]",
                                          nofMomBins,momBins,
                                          nofRadiusBins,rRadBins,
                                          nofThetaBins,theBins);
  for ( Int_t istat = 0 ; istat < fDigiPar->GetNStations() ; istat++ ) {
    fhThetaVsRadiusVsMomentum[istat] = new TH3F(Form("fhThetaVsRadiusVsMomentum_s%d",istat+1),
                                                Form("Track Theta vs Point Radius vs Track Momentum on station %d;p [GeV/c];r [cm];#theta [#circ]",istat+1),
                                                nofMomBins,momBins,
                                                nofRadiusBins,radBins,
                                                nofThetaBins,theBins);
    ffThetaVsRadiusVsMomentum[istat] = new TF2 (Form("ffThetaVsRadiusVsMomentum_s%d",istat+1),
                                                "([0]+1/(x+[1]*[4]+[2]))*y/[4]+[3]",0.,20.,minRRadius,maxRRadius);
    fHistoList->Add(fhThetaVsRadiusVsMomentum[istat]);
    fHistoList->Add(ffThetaVsRadiusVsMomentum[istat]);
  }
  fHistoList->Add(fhThetaVsRadiusVsMomentumAll);



  fhRadiusVsAngle = new TH2F("fhRadiusVsAngle","Radius ratio/z ratio vs phi angle difference;#phi_{2}-#phi_{1} [#circ];#frac{r_{2}}{r_{1}}#frac{z_{1}}{z_{2}}",800,-40.,40.,200,0.,2.);
  fHistoList->Add(fhRadiusVsAngle);
  Int_t crHist = 0; 
  for ( Int_t istat1 = 0 ; istat1 < fDigiPar->GetNStations() ; istat1++ ) {
    for ( Int_t istat2 = istat1+1 ; istat2 < fDigiPar->GetNStations() ; istat2++ ) {
      fhRadiusVsAnglePair[crHist] = new TH2F(Form("fhRadiusVsAnglePair_s%d_s%d",istat1+1,istat2+1),
                                             Form("Radius ratio/z ratio vs phi angle difference, stat%d vs stat%d;#phi_{2} - #phi_{1} [#circ];#frac{r_{2}}{r_{1}}#frac{z_{1}}{z_{2}}",istat1+1,istat2+1),
                                             800,-40.,40.,200,0.,2.);
      fHistoList->Add(fhRadiusVsAnglePair[crHist]);

      fhMomentumVsPhiDiffVsRadius[crHist] = new TH3F(Form("fhMomentumVsPhiDiffVsRadius_s%d_s%d",istat1+1,istat2+1),
                                                     Form("Momentum vs phi angle difference vs radius, stat%d vs stat%d;r [cm];#phi_{2} - #phi_{1} [#circ];p [GeV/c]",istat1+1,istat2+1),
                                                     nofRadiusBins,radBins,
                                                     nofPhiBins,phiBins,
                                                     nofMomBins,momBins);
      fHistoList->Add(fhMomentumVsPhiDiffVsRadius[crHist]);
                      

      ffMomentumVsPhiDiffVsRadius[crHist] = new TF2(Form("ffMomentumVsPhiDiffVsRadius_s%d_s%d",istat1+1,istat2+1),
                                                    "([0]+[1]/x)/y",0.,100.,0.,40.);
      fHistoList->Add(ffMomentumVsPhiDiffVsRadius[crHist]);

      fhMomentumVsPhiDiffAll[crHist] = new TH2F(Form("fhMomentumVsPhiDiffAll_s%d_s%d",istat1+1,istat2+1),
                                                Form("Momentum vs phi angle difference, stat%d vs stat%d;#phi_{2} - #phi_{1} [#circ];p [GeV/c]",istat1+1,istat2+1),
                                                nofDetPhiBins,minPhi,maxPhi,
                                                nofDetMomBins,minMom,maxMom);
      fHistoList->Add(fhMomentumVsPhiDiffAll[crHist]);

      ffMomentumVsPhiDiffAll[crHist] = new TF1(Form("ffMomentumVsPhiDiffAll_s%d_s%d",istat1+1,istat2+1),
                                                    "[0]/x",0.,40.);
      fHistoList->Add(ffMomentumVsPhiDiffAll[crHist]);

      for ( Int_t irad = 0 ; irad < 20 ; irad++ ) {
        fhMomentumVsPhiDiff[crHist][irad] = new TH2F(Form("fhMomentumVsPhiDiff_s%d_s%d_r%d",istat1+1,istat2+1,irad),
                                                     Form("Momentum vs phi angle difference, stat%d vs stat%d, rad %d to %d;#phi_{2} - #phi_{1} [#circ];p [GeV/c]",istat1+1,istat2+1,irad*5,(irad+1)*5),
                                                     nofDetPhiBins,minPhi,maxPhi,
                                                     nofDetMomBins,minMom,maxMom);
        fHistoList->Add(fhMomentumVsPhiDiff[crHist][irad]);

        ffMomentumVsPhiDiff[crHist][irad] = new TF1(Form("ffMomentumVsPhiDiff_s%d_s%d_r%d",istat1+1,istat2+1,irad),
                                                 "[0]/x",0.,40.);
        fHistoList->Add(ffMomentumVsPhiDiff[crHist][irad]);
      }
      
      fhTrackPhiVsHitPhis   [crHist] = new TH3F(Form("fhTrackPhiVsHitPhis_s%d_s%d",istat1+1,istat2+1),
                                                Form("Track phi vs hit1 phi vs phi difference/z ratio, stat%d vs stat%d;#phi_{1} [#circ];(#phi_{2} - #phi_{1})*z__{2}/z_{1} [#circ];#phi_{T} [#circ]",istat1+1,istat2+1),
                                                nofPhi1Bins,phi1Bins,
                                                nofPhi2Bins,phi2Bins,
                                                nofPhi1Bins,phi1Bins);
      fhTrackPhiVsCalcPhi   [crHist] = new TH2F(Form("fhTrackPhiVsCalcPhi_s%d_s%d",istat1+1,istat2+1),
                                                Form("Track MC phi vs track calculated phi, stat%d vs stat%d;#phi_{calc} [#circ];#phi_{MC} [#circ]",istat1+1,istat2+1),
                                                nofPhi1Bins,phi1Bins,
                                                nofPhi1Bins,phi1Bins);
      fHistoList->Add(fhTrackPhiVsHitPhis[crHist]);
      fHistoList->Add(fhTrackPhiVsCalcPhi[crHist]);


      fhMomTransVsHRadius[crHist] = new TH2F(Form("fhMomTransVsHRadius_s%d_s%d",istat1+1,istat2+1),
                                             Form("Track MC momentum trans vs track calculated helix radius;radius_{radius} [cm];pt_{MC} [GeV/c]"),
                                             nofDetMomBins,0.,800.,
                                             nofDetMomBins,0.,maxMom/3.);
      fHistoList->Add(fhMomTransVsHRadius[crHist]);

      crHist++;
    }
  }
  for ( Int_t itemp = 0 ; itemp < 3 ; itemp++ ) {
    ffRadiusVsAngle[itemp] = new TF1(Form("ffRadiusVsAngle%d",itemp),"pol2",-40.,40.);
    fHistoList->Add(ffRadiusVsAngle[itemp]); 
  }


}

void PndGemMagneticFieldVsTrackParameters::Exec ( Option_t *  opt )

virtual

Execution

Definition at line 87 of file PndGemMagneticFieldVsTrackParameters.cxx.

References Fill1StationHistograms(), Fill2StationsHistograms(), and Reset().

                                                         {

  // Size of fMCTrackArray
  //cout <<"# MC Tracks: "<< fMCTrackArray->GetEntriesFast() << endl;
  // Size of fMCTrackArray
  //cout <<"# MC Points: "<< fMCPointArray->GetEntriesFast() << endl;
  
  Fill1StationHistograms();
  Fill2StationsHistograms();
  
  Reset();
  
}

Int_t PndGemMagneticFieldVsTrackParameters::Fill1StationHistograms ( )

private

Definition at line 188 of file PndGemMagneticFieldVsTrackParameters.cxx.

References Double_t, fDigiPar, fhThetaVsRadiusVsMomentum, fhThetaVsRadiusVsMomentumAll, fMCPointArray, fMCTrackArray, PndMCTrack::GetMomentum(), PndGemMCPoint::GetSensorId(), PndGemDigiPar::GetStationNr(), and CAMath::Sqrt().

Referenced by Exec().

                                                                   {
  Int_t nofGemPoints = fMCPointArray->GetEntriesFast();
  PndGemMCPoint* mcPoint;
  PndMCTrack*    mcTrack;

  for ( Int_t imcp = 0 ; imcp < nofGemPoints ; imcp++ ) {
    mcPoint = (PndGemMCPoint*)fMCPointArray->At(imcp);

    Int_t stationNumber = fDigiPar->GetStationNr(mcPoint->GetSensorId());

    mcTrack = (PndMCTrack*)fMCTrackArray->At(mcPoint->GetTrackID());
    TVector3 mcMomVec = mcTrack->GetMomentum();

    Double_t pRadius = TMath::Sqrt(mcPoint->GetX()*mcPoint->GetX()+mcPoint->GetY()*mcPoint->GetY());
    Double_t pZ      = mcPoint->GetZ();
    Double_t tMom    = mcMomVec.Mag();
    Double_t tTheta  = mcMomVec.Theta()*TMath::RadToDeg();

    fhThetaVsRadiusVsMomentum[stationNumber-1]->Fill(tMom,pRadius,   tTheta);
    fhThetaVsRadiusVsMomentumAll              ->Fill(tMom,pRadius/pZ,tTheta);

  }
  return 0;
}

Int_t PndGemMagneticFieldVsTrackParameters::Fill2StationsHistograms ( )

private

Definition at line 215 of file PndGemMagneticFieldVsTrackParameters.cxx.

References CAMath::Abs(), Double_t, fDigiPar, fhMomentumVsPhiDiff, fhMomentumVsPhiDiffAll, fhMomentumVsPhiDiffVsRadius, fhMomTransVsHRadius, fhRadiusVsAngle, fhRadiusVsAnglePair, fhTrackPhiVsCalcPhi, fhTrackPhiVsHitPhis, fMCPointArray, fMCTrackArray, PndMCTrack::GetMomentum(), PndMCTrack::GetMotherID(), PndGemDigiPar::GetNStations(), PndGemDigiPar::GetStation(), PndGemStation::GetZ(), Pi, CAMath::Sin(), and CAMath::Sqrt().

Referenced by Exec().

                                                                    {
  Int_t nofGemPoints = fMCPointArray->GetEntriesFast();
  PndGemMCPoint* mcPoint1;
  PndGemMCPoint* mcPoint2;
  PndMCTrack* mcTrack;

  for ( Int_t imcp1 = 0 ; imcp1 < nofGemPoints ; imcp1++ ) {
    mcPoint1 = (PndGemMCPoint*)fMCPointArray->At(imcp1);

    Double_t p1Radius = TMath::Sqrt(mcPoint1->GetX()*mcPoint1->GetX()+mcPoint1->GetY()*mcPoint1->GetY());
    Double_t p1Z      = mcPoint1->GetZ();
    Double_t p1PhiAng = TMath::ACos(mcPoint1->GetX()/p1Radius);
    if ( mcPoint1->GetY() < 0 )
      p1PhiAng = 2.*TMath::Pi() - p1PhiAng;
    p1PhiAng *= TMath::RadToDeg();

    Int_t mcTrackId = mcPoint1->GetTrackID();

    mcTrack = (PndMCTrack*)fMCTrackArray->At(mcTrackId);

    if ( mcTrack->GetMotherID() != -1 ) continue;

    TVector3 mcMomVec = mcTrack->GetMomentum();
    Double_t tMom    = mcMomVec.Mag();
    Double_t  tPt    = mcMomVec.Pt();
    Double_t pMom    = TMath::RadToDeg()*mcMomVec.Phi() + (mcMomVec.Phi()>=0?0.:360.);

    for ( Int_t imcp2 = imcp1+1 ; imcp2 < nofGemPoints ; imcp2++ ) {
      mcPoint2 = (PndGemMCPoint*)fMCPointArray->At(imcp2);
      Double_t p2Z      = mcPoint2->GetZ();
      if ( TMath::Abs(p2Z-p1Z) < 5. ) continue;

      Double_t p2Radius = TMath::Sqrt(mcPoint2->GetX()*mcPoint2->GetX()+mcPoint2->GetY()*mcPoint2->GetY());

      Double_t p2PhiAng = TMath::ACos(mcPoint2->GetX()/p2Radius);
      if ( mcPoint2->GetY() < 0 )
        p2PhiAng = 2.*TMath::Pi() - p2PhiAng;
      p2PhiAng *= TMath::RadToDeg();

      if ( p1PhiAng <  90. && p2PhiAng > 270. ) p2PhiAng -= 360.;
      if ( p1PhiAng > 270. && p2PhiAng <  90. ) p2PhiAng += 360.;

      fhRadiusVsAngle->Fill((p2PhiAng-p1PhiAng)*p1Z/p2Z,p2Radius*p1Z/(p1Radius*p2Z));

      Int_t histNo = -1;
      Int_t ihist  =  0;
      for ( Int_t istat1 = 0 ; istat1 < fDigiPar->GetNStations() ; istat1++ ) {
        PndGemStation* stat1 = (PndGemStation*)fDigiPar->GetStation(istat1);
        Double_t zStation1 = stat1->GetZ();
        for ( Int_t istat2 = istat1+1 ; istat2 < fDigiPar->GetNStations() ; istat2++ ) {
          PndGemStation* stat2 = (PndGemStation*)fDigiPar->GetStation(istat2);
          Double_t zStation2 = stat2->GetZ();
          if ( TMath::Abs(p1Z-zStation1) < 2. && TMath::Abs(p2Z-zStation2) < 2. ) histNo = ihist;      
          ihist++;
        }
      }
      if ( histNo == -1 ) continue;
      fhRadiusVsAnglePair[histNo]->Fill((p2PhiAng-p1PhiAng)*p1Z/p2Z,p2Radius*p1Z/(p1Radius*p2Z));

      //      Double_t circRad = TMath::Sqrt( (gemHit2X-gemHit1X)*(gemHit2X-gemHit1X) + (gemHit2Y-gemHit1Y)*(gemHit2Y-gemHit1Y) ) / TMath::Sin(pangle-pangle2) / 2.;
      fhMomTransVsHRadius[histNo]->Fill(TMath::Abs(TMath::Sqrt((mcPoint2->GetX()-mcPoint1->GetX())*(mcPoint2->GetX()-mcPoint1->GetX())+
                                                               (mcPoint2->GetY()-mcPoint1->GetY())*(mcPoint2->GetY()-mcPoint1->GetY())) /
                                                   TMath::Sin(TMath::DegToRad()*(p2PhiAng-p1PhiAng)) / 2.),
                                        tPt);

      if ( mcTrackId != mcPoint2->GetTrackID() ) 
        continue;

      fhMomentumVsPhiDiffVsRadius[histNo]->Fill(p1Radius,TMath::Abs(p2PhiAng-p1PhiAng),tMom);

      fhMomentumVsPhiDiffAll[histNo]                    ->Fill(TMath::Abs(p2PhiAng-p1PhiAng),tMom);
      fhMomentumVsPhiDiff   [histNo][Int_t(p1Radius/5.)]->Fill(TMath::Abs(p2PhiAng-p1PhiAng),tMom);

      fhTrackPhiVsHitPhis   [histNo]                    ->Fill(p1PhiAng,(p1PhiAng-p2PhiAng)*p1Z/(p2Z-p1Z),pMom);
      fhTrackPhiVsCalcPhi   [histNo]                    ->Fill(p1PhiAng+(p1PhiAng-p2PhiAng)*p1Z/(p2Z-p1Z),pMom);

    }
  }
  return 0;
}

void PndGemMagneticFieldVsTrackParameters::Finish ( )

privatevirtual

Finish

Definition at line 435 of file PndGemMagneticFieldVsTrackParameters.cxx.

References Double_t, fDigiPar, ffMomentumVsPhiDiff, ffMomentumVsPhiDiffAll, ffMomentumVsPhiDiffVsRadius, ffRadiusVsAngle, ffThetaVsRadiusVsMomentum, fHistoList, fTrackFinderOnHits_ParMat0, fTrackFinderOnHits_ParMat1, fTrackFinderOnHits_ParRadPhi0, fTrackFinderOnHits_ParRadPhi2, fTrackFinderOnHits_ParTheta0, fTrackFinderOnHits_ParTheta1, fTrackFinderOnHits_ParTheta2, fTrackFinderOnHits_ParTheta3, PndGemDigiPar::GetNStations(), PndGemDigiPar::GetStation(), PndGemStation::GetZ(), next, and par.

                                                  {
 
  for ( Int_t istat = 0 ; istat < fDigiPar->GetNStations() ; istat++ ) {
    PndGemStation* station = (PndGemStation*)fDigiPar->GetStation(istat);
    Double_t zStation = station->GetZ();
    
    ffThetaVsRadiusVsMomentum[istat]->SetParameters(fTrackFinderOnHits_ParTheta0,
                                                    fTrackFinderOnHits_ParTheta1,
                                                    fTrackFinderOnHits_ParTheta2,
                                                    fTrackFinderOnHits_ParTheta3,
                                                    zStation);
  }
  for ( Int_t itemp = 0 ; itemp < 3 ; itemp++ ) {
    ffRadiusVsAngle[itemp]->SetParameters(fTrackFinderOnHits_ParRadPhi0+(-1+itemp)*0.05,0.,fTrackFinderOnHits_ParRadPhi2);
  }
  Int_t crHist = 0; 
  for ( Int_t istat1 = 0 ; istat1 < fDigiPar->GetNStations() ; istat1++ ) {
    PndGemStation* stat1 = (PndGemStation*)fDigiPar->GetStation(istat1);
    Double_t zStation1 = stat1->GetZ();
    for ( Int_t istat2 = istat1+1 ; istat2 < fDigiPar->GetNStations() ; istat2++ ) {
      PndGemStation* stat2 = (PndGemStation*)fDigiPar->GetStation(istat2);
      Double_t zStation2 = stat2->GetZ();
  
      Double_t zCuDiff = zStation2*zStation2*zStation2-zStation1*zStation1*zStation1;
      Double_t zSqDiff = zStation2*zStation2-zStation1*zStation1;
      Double_t zDiff   = zStation2-zStation1;
      
      Double_t par0_mom = fTrackFinderOnHits_ParMat0[0]*zCuDiff+fTrackFinderOnHits_ParMat0[1]*zSqDiff+fTrackFinderOnHits_ParMat0[2]*zDiff;
      Double_t par1_mom = fTrackFinderOnHits_ParMat1[0]*zCuDiff+fTrackFinderOnHits_ParMat1[1]*zSqDiff+fTrackFinderOnHits_ParMat1[2]*zDiff;

      ffMomentumVsPhiDiffVsRadius[crHist]->SetParameters(par0_mom,par1_mom);
      ffMomentumVsPhiDiffAll     [crHist]->SetParameter(0,par0_mom);
      for ( Int_t irad = 0 ; irad < 20 ; irad++ ) {
        Double_t par = par0_mom + par1_mom/(Double_t(irad)*5.+2.5);
        ffMomentumVsPhiDiff   [crHist][irad]->SetParameter(0,par);
      }
      
      crHist++;
    }
  }
  
  gDirectory->mkdir("GemMFVTP");
  gDirectory->cd("GemMFVTP");
  TIter next(fHistoList);
  while ( TH1* histo = ((TH1*)next()) ) histo->Write();
  gDirectory->cd("..");
}

InitStatus PndGemMagneticFieldVsTrackParameters::Init ( )

privatevirtual

Intialisation

Definition at line 123 of file PndGemMagneticFieldVsTrackParameters.cxx.

References CreateHistos(), fDigiPar, fMCPointArray, fMCTrackArray, fTrackFinderOnHits_ParMat0, fTrackFinderOnHits_ParMat1, fTrackFinderOnHits_ParRadPhi0, fTrackFinderOnHits_ParRadPhi2, fTrackFinderOnHits_ParTheta0, fTrackFinderOnHits_ParTheta1, fTrackFinderOnHits_ParTheta2, fTrackFinderOnHits_ParTheta3, fTrackFinderOnHits_ParThetaA, fTrackFinderOnHits_ParThetaB, PndGemDigiPar::GetNStations(), PndGemDigiPar::GetTrackFinderOnHits_ParMat0(), PndGemDigiPar::GetTrackFinderOnHits_ParMat1(), PndGemDigiPar::GetTrackFinderOnHits_ParRadPhi0(), PndGemDigiPar::GetTrackFinderOnHits_ParRadPhi2(), PndGemDigiPar::GetTrackFinderOnHits_ParTheta0(), PndGemDigiPar::GetTrackFinderOnHits_ParTheta1(), PndGemDigiPar::GetTrackFinderOnHits_ParTheta2(), PndGemDigiPar::GetTrackFinderOnHits_ParTheta3(), PndGemDigiPar::GetTrackFinderOnHits_ParThetaA(), and PndGemDigiPar::GetTrackFinderOnHits_ParThetaB().

                                                      {

  // Get input array 
  FairRootManager* ioman = FairRootManager::Instance();
  if ( ! ioman ) Fatal("Init", "No FairRootManager");

  // Get MCPoint array
  fMCPointArray = (TClonesArray*) ioman->GetObject("GEMPoint");
  if( !fMCPointArray ) {
    cout << "-E- "<< GetName() <<"::Init: No MCPoint array!"
         << endl;
    return kERROR;
  }

  // Get MCTrack array
  fMCTrackArray  = (TClonesArray*) ioman->GetObject("MCTrack");  
  if( !fMCTrackArray ) {
    cout << "-E- "<< GetName() <<"::Init: No MCTrack array!"
         << endl;
    return kERROR;
  }

  fTrackFinderOnHits_ParThetaA = fDigiPar->GetTrackFinderOnHits_ParThetaA();
  fTrackFinderOnHits_ParThetaB = fDigiPar->GetTrackFinderOnHits_ParThetaB();
  
  fTrackFinderOnHits_ParTheta0 = fDigiPar->GetTrackFinderOnHits_ParTheta0();
  fTrackFinderOnHits_ParTheta1 = fDigiPar->GetTrackFinderOnHits_ParTheta1();
  fTrackFinderOnHits_ParTheta2 = fDigiPar->GetTrackFinderOnHits_ParTheta2();
  fTrackFinderOnHits_ParTheta3 = fDigiPar->GetTrackFinderOnHits_ParTheta3();
  
  fTrackFinderOnHits_ParRadPhi0 = fDigiPar->GetTrackFinderOnHits_ParRadPhi0();
  fTrackFinderOnHits_ParRadPhi2 = fDigiPar->GetTrackFinderOnHits_ParRadPhi2();
  for ( Int_t in = 0 ; in < 3 ; in++ ) {
    fTrackFinderOnHits_ParMat0[in] = fDigiPar->GetTrackFinderOnHits_ParMat0(in);
    fTrackFinderOnHits_ParMat1[in] = fDigiPar->GetTrackFinderOnHits_ParMat1(in);
  }

  CreateHistos();

  cout << "-I- " << fName.Data() << "::Init(). There are " << fDigiPar->GetNStations() << " GEM stations." << endl;
  cout << "-I- " << fName.Data() << "::Init(). Initialization succesfull." << endl;

  return kSUCCESS;

}

InitStatus PndGemMagneticFieldVsTrackParameters::ReInit ( )

privatevirtual

Reinitialisation

Definition at line 173 of file PndGemMagneticFieldVsTrackParameters.cxx.

                                                        {
  return kSUCCESS;
}

void PndGemMagneticFieldVsTrackParameters::Reset ( )

private

Reset eventwise counters

Definition at line 181 of file PndGemMagneticFieldVsTrackParameters.cxx.

References fNofEvents.

Referenced by Exec(), PndGemMagneticFieldVsTrackParameters(), and ~PndGemMagneticFieldVsTrackParameters().

                                                 {
  fNofEvents = 0;
}

void PndGemMagneticFieldVsTrackParameters::SetParContainers ( )

privatevirtual

Get parameter containers

Definition at line 105 of file PndGemMagneticFieldVsTrackParameters.cxx.

References fDigiPar, and run.

                                                            {
  
  // Get run and runtime database
  FairRunAna* run = FairRunAna::Instance();
  if ( ! run ) Fatal("SetParContainers", "No analysis run");

  FairRuntimeDb* db = run->GetRuntimeDb();
  if ( ! db ) Fatal("SetParContainers", "No runtime database");

  // Get GEM digitisation parameter container
  fDigiPar = (PndGemDigiPar*)(db->getContainer("PndGemDetectors"));

}

Member Data Documentation

PndGemDigiPar* PndGemMagneticFieldVsTrackParameters::fDigiPar

private

Definition at line 72 of file PndGemMagneticFieldVsTrackParameters.h.

Referenced by CreateHistos(), Fill1StationHistograms(), Fill2StationsHistograms(), Finish(), Init(), PndGemMagneticFieldVsTrackParameters(), SetParContainers(), and ~PndGemMagneticFieldVsTrackParameters().

TF1* PndGemMagneticFieldVsTrackParameters::ffMomentumVsPhiDiff[10][20]

private

Definition at line 112 of file PndGemMagneticFieldVsTrackParameters.h.

Referenced by CreateHistos(), and Finish().

TF1* PndGemMagneticFieldVsTrackParameters::ffMomentumVsPhiDiffAll[10]

private

Definition at line 111 of file PndGemMagneticFieldVsTrackParameters.h.

Referenced by CreateHistos(), and Finish().

TF2* PndGemMagneticFieldVsTrackParameters::ffMomentumVsPhiDiffVsRadius[10]

private

Definition at line 110 of file PndGemMagneticFieldVsTrackParameters.h.

Referenced by CreateHistos(), and Finish().

TF1* PndGemMagneticFieldVsTrackParameters::ffRadiusVsAngle[3]

private

Definition at line 105 of file PndGemMagneticFieldVsTrackParameters.h.

Referenced by CreateHistos(), and Finish().

TF2* PndGemMagneticFieldVsTrackParameters::ffThetaVsRadiusVsMomentum[10]

private

Definition at line 101 of file PndGemMagneticFieldVsTrackParameters.h.

Referenced by CreateHistos(), and Finish().

TList* PndGemMagneticFieldVsTrackParameters::fHistoList

private

Definition at line 97 of file PndGemMagneticFieldVsTrackParameters.h.

Referenced by CreateHistos(), and Finish().

TH2F* PndGemMagneticFieldVsTrackParameters::fhMomentumVsPhiDiff[10][20]

private

Definition at line 109 of file PndGemMagneticFieldVsTrackParameters.h.

Referenced by CreateHistos(), and Fill2StationsHistograms().

TH2F* PndGemMagneticFieldVsTrackParameters::fhMomentumVsPhiDiffAll[10]

private

Definition at line 108 of file PndGemMagneticFieldVsTrackParameters.h.

Referenced by CreateHistos(), and Fill2StationsHistograms().

TH3F* PndGemMagneticFieldVsTrackParameters::fhMomentumVsPhiDiffVsRadius[10]

private

Definition at line 107 of file PndGemMagneticFieldVsTrackParameters.h.

Referenced by CreateHistos(), and Fill2StationsHistograms().

TH2F* PndGemMagneticFieldVsTrackParameters::fhMomTransVsHRadius[10]

private

Definition at line 114 of file PndGemMagneticFieldVsTrackParameters.h.

Referenced by CreateHistos(), and Fill2StationsHistograms().

TH2F* PndGemMagneticFieldVsTrackParameters::fhRadiusVsAngle

private

Definition at line 103 of file PndGemMagneticFieldVsTrackParameters.h.

Referenced by CreateHistos(), and Fill2StationsHistograms().

TH2F* PndGemMagneticFieldVsTrackParameters::fhRadiusVsAnglePair[20]

private

Definition at line 104 of file PndGemMagneticFieldVsTrackParameters.h.

Referenced by CreateHistos(), and Fill2StationsHistograms().

TH3F* PndGemMagneticFieldVsTrackParameters::fhThetaVsRadiusVsMomentum[10]

private

Definition at line 100 of file PndGemMagneticFieldVsTrackParameters.h.

Referenced by CreateHistos(), and Fill1StationHistograms().

TH3F* PndGemMagneticFieldVsTrackParameters::fhThetaVsRadiusVsMomentumAll

private

Definition at line 99 of file PndGemMagneticFieldVsTrackParameters.h.

Referenced by CreateHistos(), and Fill1StationHistograms().

TH2F* PndGemMagneticFieldVsTrackParameters::fhTrackPhiVsCalcPhi[10]

private

Definition at line 117 of file PndGemMagneticFieldVsTrackParameters.h.

Referenced by CreateHistos(), and Fill2StationsHistograms().

TH3F* PndGemMagneticFieldVsTrackParameters::fhTrackPhiVsHitPhis[10]

private

Definition at line 116 of file PndGemMagneticFieldVsTrackParameters.h.

Referenced by CreateHistos(), and Fill2StationsHistograms().

TClonesArray* PndGemMagneticFieldVsTrackParameters::fMCPointArray

private

Definition at line 89 of file PndGemMagneticFieldVsTrackParameters.h.

Referenced by Fill1StationHistograms(), Fill2StationsHistograms(), Init(), and PndGemMagneticFieldVsTrackParameters().

TClonesArray* PndGemMagneticFieldVsTrackParameters::fMCTrackArray

private

Definition at line 88 of file PndGemMagneticFieldVsTrackParameters.h.

Referenced by Fill1StationHistograms(), Fill2StationsHistograms(), Init(), and PndGemMagneticFieldVsTrackParameters().

Int_t PndGemMagneticFieldVsTrackParameters::fNofEvents

private

Event counter

Definition at line 94 of file PndGemMagneticFieldVsTrackParameters.h.

Referenced by Reset().

TStopwatch PndGemMagneticFieldVsTrackParameters::fTimer

private

Definition at line 91 of file PndGemMagneticFieldVsTrackParameters.h.

Double_t PndGemMagneticFieldVsTrackParameters::fTrackFinderOnHits_ParMat0[3]

private

Definition at line 85 of file PndGemMagneticFieldVsTrackParameters.h.

Referenced by Finish(), and Init().

Double_t PndGemMagneticFieldVsTrackParameters::fTrackFinderOnHits_ParMat1[3]

private

Definition at line 86 of file PndGemMagneticFieldVsTrackParameters.h.

Referenced by Finish(), and Init().

Double_t PndGemMagneticFieldVsTrackParameters::fTrackFinderOnHits_ParRadPhi0

private

Definition at line 82 of file PndGemMagneticFieldVsTrackParameters.h.

Referenced by Finish(), and Init().

Double_t PndGemMagneticFieldVsTrackParameters::fTrackFinderOnHits_ParRadPhi2

private

Definition at line 83 of file PndGemMagneticFieldVsTrackParameters.h.

Referenced by Finish(), and Init().

Double_t PndGemMagneticFieldVsTrackParameters::fTrackFinderOnHits_ParTheta0

private

Definition at line 77 of file PndGemMagneticFieldVsTrackParameters.h.

Referenced by Finish(), and Init().

Double_t PndGemMagneticFieldVsTrackParameters::fTrackFinderOnHits_ParTheta1

private

Definition at line 78 of file PndGemMagneticFieldVsTrackParameters.h.

Referenced by Finish(), and Init().

Double_t PndGemMagneticFieldVsTrackParameters::fTrackFinderOnHits_ParTheta2

private

Definition at line 79 of file PndGemMagneticFieldVsTrackParameters.h.

Referenced by Finish(), and Init().

Double_t PndGemMagneticFieldVsTrackParameters::fTrackFinderOnHits_ParTheta3

private

Definition at line 80 of file PndGemMagneticFieldVsTrackParameters.h.

Referenced by Finish(), and Init().

Double_t PndGemMagneticFieldVsTrackParameters::fTrackFinderOnHits_ParThetaA

private

Definition at line 74 of file PndGemMagneticFieldVsTrackParameters.h.

Referenced by Init().

Double_t PndGemMagneticFieldVsTrackParameters::fTrackFinderOnHits_ParThetaB

private

Definition at line 75 of file PndGemMagneticFieldVsTrackParameters.h.

Referenced by Init().

---

The documentation for this class was generated from the following files:

• PndGemMagneticFieldVsTrackParameters.h
• PndGemMagneticFieldVsTrackParameters.cxx