#include "PndFTSCAGBTracker.h"
#include "PndFTSCAGBHit.h"
#include "PndFTSCAGBTrack.h"
#include "PndFTSArray.h"
#include "PndFTSCAMath.h"
#include "PndFTSCATrackLinearisationVector.h"
#include "PndFTSCAPerformance.h"
#include "TStopwatch.h"
#include "L1Timer.h"
#include "FTSCATarget.h"
#include "FTSCAHits.h"
#include "FTSCAHitsV.h"
#include "FTSCATracks.h"
#include "Math/SMatrix.h"
#include "TMatrixD.h"
#include "TVectorD.h"
#include <algorithm>
#include <fstream>
#include <iostream>
#include <sstream>
#include "PndFTSCAMCPoint.h"
#include "PndFTSPerformanceBase.h"

Classes
struct	TrackHitRecord

Functions
t	SetQPt (float_v(1.e-8f), CAMath::Abs(t.QPt())< 1.e-8f)

	for (unsigned char i=0;i< 15;i++) ok &

bool	IsRightNeighbour (const FTSCANPlet &a, const FTSCANPlet &b, float pick, float &chi2)

Variables
bool	SINGLE_THREADED = false

float_m	ok = active0

const int	StartStationShift = 1

vector< TrackHitRecord >	gTrackHitRecords

Function Documentation

for ( )

Definition at line 156 of file anal_hit_digi_cluster_fwendcap.C.

References b, c, PndMvdAdvancedPidAlgo::CalcLikelihood(), cemc, cluster_energy, CAMath::Cos(), dE, digi, digi_ene, PndEmcCluster::DigiList(), Double_t, dx, dxAbs, dyAbs, dySci, dzAbs, dzSci, PndEmcCluster::energy(), evt, f, fTrackParFinal, fTrackParGeane, PndMCTrack::Get4Momentum(), GFTrack::getCardinalRep(), PndEmcHit::GetCopy(), PndEmcDigi::GetCopy(), PndEmcHit::GetCrystal(), PndEmcDigi::GetCrystal(), PndEmcDigi::GetDetectorId(), PndEmcHit::GetEnergy(), PndEmcDigi::GetEnergy(), PndMvdPidCand::GetLikelihood(), PndEmcHit::GetModule(), PndEmcDigi::GetModule(), GFAbsTrackRep::getMom(), PndMCTrack::GetMomentum(), PndMvdPidCand::GetMvdHitdE(), PndMvdPidCand::GetMvdHitdx(), PndMvdPidCand::GetMvdHitMomentum(), PndMvdPidCand::GetMvdHits(), PndEmcDigi::GetPhi(), PndEmcDigi::GetPhiInt(), PndSdsMCPoint::GetPosition(), PndSdsMCPoint::GetPositionOut(), PndEmcHit::GetRow(), PndEmcDigi::GetRow(), GFAbsTrackRep::getStatusFlag(), PndEmcDigi::GetTheta(), PndEmcDigi::GetThetaInt(), GFTrack::getTrackRep(), PndEmcHit::GetXPad(), PndEmcDigi::GetXPad(), PndEmcHit::GetYPad(), PndEmcDigi::GetYPad(), hit, hit_ene, hit_pos(), i, kSpaceInSub, mctrack, name, nChFE, nFEBox, nhits, nsteps, offset(), offset1(), out, partGen, Pi, pos, PndMvdConvertApv::ReadNext(), rotNamech, CAMath::Sin(), thickness, TString, PndEmcCluster::where(), PndEmcDigi::where(), Y, PndEmcCluster::z(), and Z.

Referenced by PndFilteredPrimaryGenerator::AddFilter(), operator<<(), PndFsmCmpDet::readParameters(), and PndFsmCombiDet::readParameters().

                                    {
    //for (Int_t j=1;j< 3; j++){
      Float_t sum_hit_ene,hit_ene_cp1,hit_ene_cp2,hit_ene_cp3,hit_ene_cp4;
      Float_t sum_digi_ene,digi_ene_cp1,digi_ene_cp2,digi_ene_cp3,digi_ene_cp4;
      Float_t sum_clus_ene=0;
      Float_t cemc_sum=0;
 
      sum_hit_ene = hit_ene_cp1 = hit_ene_cp2 = hit_ene_cp3 = hit_ene_cp4 = 0;
      sum_digi_ene = digi_ene_cp1 = digi_ene_cp2 = digi_ene_cp3 = digi_ene_cp4 =0;
 
      tsim->GetEntry(j);
      tful->GetEntry(j);
 
         PndMCTrack *mctrack=(PndMCTrack *) mctrack_array->At(0);
         photon_momentum=mctrack->GetMomentum();
         theta_mc=photon_momentum.Theta()*(180./TMath::Pi());
         phi_mc=photon_momentum.Phi()*(180./TMath::Pi());
 
         p4mom=mctrack->Get4Momentum();
 
         Int_t nhits = hit_array->GetEntries();
         for (Int_t i=0; i<nhits; i++)
           {
             PndEmcHit *hit=(PndEmcHit*)hit_array->At(i);
 
             module_h = hit->GetModule();
 
             if (module_h==3){
 
               hene_mc->Fill(p4mom.E());
 
               hit_ene=hit->GetEnergy();
 
               hene_h->Fill(hit_ene);
 
               sum_hit_ene+=hit_ene;
 
               crystal_h = hit->GetCrystal();
               row_h = hit->GetRow();
               copy_h = hit->GetCopy();
 
               if (copy_h==1) hit_ene_cp1+=hit->GetEnergy();
               if (copy_h==2) hit_ene_cp2+=hit->GetEnergy();
               if (copy_h==3) hit_ene_cp3+=hit->GetEnergy();
               if (copy_h==4) hit_ene_cp4+=hit->GetEnergy();
 
               TVector3 hit_pos(hit->GetX(),hit->GetY(),hit->GetZ());
               hit_theta=hit_pos.Theta()*(180./TMath::Pi());
               hit_phi=hit_pos.Phi()*(180./TMath::Pi());
 
               //cout << "MC: theta= "<<theta_mc<<", phi= "<<phi_mc<<endl;
               h10->Fill(theta_mc);
               h20->Fill(phi_mc);
 
               h10d->Fill(hit_theta);
               h20d->Fill(hit_phi);
 
               h1->Fill(hit_theta-theta_mc);
               h2->Fill(hit_phi-phi_mc);
 
               if (copy_h==1) hz_cp1_h->Fill(hit->GetZ());
               if (copy_h==2) hz_cp2_h->Fill(hit->GetZ());
               if (copy_h==3) hz_cp3_h->Fill(hit->GetZ());
               if (copy_h==4) hz_cp4_h->Fill(hit->GetZ());
 
               hz_h->Fill(hit->GetZ());
               hcp_z_h->Fill(copy_h,hit->GetZ());
 
               hxpad_x_h->Fill(hit->GetXPad(),hit->GetX());
               hypad_y_h->Fill(hit->GetYPad(),hit->GetY());
 
               hx_row_h->Fill(hit->GetRow(),hit->GetX());
 
               id_h = hit->GetDetectorID();
 
               crystal_hid = (id_h%10000);
               row_hid = ((id_h/1000000)%100);
 
               hmodule_h->Fill(module_h);
               hrow_crystal_h->Fill(crystal_hid,row_hid);
               hx_y_h->Fill(hit->GetX(),hit->GetY());
 
               hid_h->Fill(id_h);
             }
           }
         hhit_ene->Fill(sum_hit_ene);
         hhit_ene_cp1->Fill(hit_ene_cp1);
         hhit_ene_cp2->Fill(hit_ene_cp2);
         hhit_ene_cp3->Fill(hit_ene_cp3);
         hhit_ene_cp4->Fill(hit_ene_cp4);
 
         // digits
         Int_t ndigits = digi_array->GetEntries();
         //cout << "DIGI array ==  " << ndigits<< endl;
         for (Int_t i=0; i<ndigits; i++)
           {
             PndEmcDigi *digi=(PndEmcDigi*)digi_array->At(i);
 
             TVector3 digi_pos=digi->where();
 
             module_d = digi->GetModule();
 
             if (module_d==3){
 
               dene_mc->Fill(p4mom.E());
 
               digi_ene=digi->GetEnergy();
               hene_d->Fill(digi_ene);
 
               sum_digi_ene+=digi_ene;
 
               crystal_d = digi->GetCrystal();
               row_d = digi->GetRow();
               copy_d = digi->GetCopy();
 
               Double_t digi_z=digi_pos.Z();
               hz_d->Fill(digi_z);
 
               if (copy_d==1){
                 digi_ene_cp1+=digi->GetEnergy();
                 hz_cp1_d->Fill(digi_z);
               }
               if (copy_d==2){
                 digi_ene_cp2+=digi->GetEnergy();
                 hz_cp2_d->Fill(digi_z);
               }
               if (copy_d==3){
                 digi_ene_cp3+=digi->GetEnergy();
                 hz_cp3_d->Fill(digi_z);
               }
               if (copy_d==4){
                 digi_ene_cp4+=digi->GetEnergy();
                 hz_cp4_d->Fill(digi_z);
               }
 
               digi_theta=(digi->GetTheta())*(180./TMath::Pi());
               digi_phi=(digi->GetPhi())*(180./TMath::Pi());
 
               h10digi->Fill(theta_mc);
               h20digi->Fill(phi_mc);
 
               h10ddigi->Fill(digi_theta);
               h20ddigi->Fill(digi_phi);
 
               h1d->Fill(digi_theta-theta_mc);
               h2d->Fill(digi_phi-phi_mc);
 
               hxpad_x_d->Fill(digi->GetXPad(),digi->GetThetaInt());
               hypad_y_d->Fill(digi->GetYPad(),digi->GetPhiInt());
 
               hx_row_d->Fill(digi->GetRow(),digi->GetThetaInt());
 
               id_d = digi->GetDetectorId();
 
               crystal_did = (id_d%10000);
               row_did = ((id_d/1000000)%100);
 
               hmodule_d->Fill(module_d);
               hrow_crystal_d->Fill(crystal_did,row_did);
               hx_y_d->Fill(digi->GetThetaInt(),digi->GetPhiInt());
 
               hid_d->Fill(id_d);
             }
           }
         hdigi_ene->Fill(sum_digi_ene);
         hdigi_ene_cp1->Fill(digi_ene_cp1);
         hdigi_ene_cp2->Fill(digi_ene_cp2);
         hdigi_ene_cp3->Fill(digi_ene_cp3);
         hdigi_ene_cp4->Fill(digi_ene_cp4);
 
         // reconstructed clusters
         Int_t ClustNoOrig,ClustNo,ClustNo_cp1,ClustNo_cp2,ClustNo_cp3,ClustNo_cp4;
         Int_t detectorID, module, mod3;
         if (cluster_array->GetEntriesFast()>0)
           {
             Int_t    idWithHighestEnergy = 0;
             Double_t highestEnergy = -1.;
 
             // First find the cluster with the highest energy
 
             if (module_h==3){
 
               ClustNoOrig=ClustNo=
                 ClustNo_cp1=ClustNo_cp2=ClustNo_cp3=ClustNo_cp4=0;
 
               ClustNoOrig=cluster_array->GetEntriesFast();
               for (Int_t i=0; i<ClustNoOrig; i++)
                 {
                   PndEmcCluster *cluster=(PndEmcCluster*)cluster_array->At(i);
                   cluster_energy=cluster->energy();
 
                   std::vector<PndEmcDigi*> digiList=cluster->DigiList();
                   ndigi=digiList.size();
 
                   if(ndigi != 0)
                     {
                       for(int k = 0; k < digiList.size(); k++)
                         {
                           detectorID = digiList[k]->GetDetectorId();
                           module = digiList[k]->GetModule();
                           row_c = digiList[k]->GetRow();
                           crystal_c = digiList[k]->GetCrystal();
                           copy_c = digiList[k]->GetCopy();
 
                           crystal_cid = (detectorID%10000);
                           row_cid = ((detectorID/1000000)%100);
 
                           hrow_crystal_c->Fill(crystal_cid,row_cid);
                           //cluster_energy=cluster->energy();
                         }
                     }
                   mod3 = module;
 
                   if (cluster_energy>highestEnergy)
                     {
                       idWithHighestEnergy = i;
                       highestEnergy = cluster_energy;
                       ClustNo++;
                     }
                 } // number of cluster
 
               // Lets analyze that cluster!
               if (mod3 !=3 )cout << "mod3 = "<< mod3 << endl;
               PndEmcCluster *cluster=(PndEmcCluster*)cluster_array->At(idWithHighestEnergy);
 
               TVector3 cluster_pos=cluster->where();
               cluster_theta=cluster_pos.Theta()*(180./TMath::Pi());
               cluster_phi=cluster_pos.Phi()*(180./TMath::Pi());//+360;
 
               cluster_energy=cluster->energy();
               //ClustNo++;
 
               if (copy_c==1){
                 hclus_ene_cp1->Fill(cluster->energy());
                 hz_cp1_c->Fill(cluster->z());
                 ClustNo_cp1++;
               }
               if (copy_c==2){
                 hclus_ene_cp2->Fill(cluster->energy());
                 hz_cp2_c->Fill(cluster->z());
                 ClustNo_cp2++;
               }
               if (copy_c==3){
                 hclus_ene_cp3->Fill(cluster->energy());
                 hz_cp3_c->Fill(cluster->z());
                 ClustNo_cp3++;
               }
               if (copy_c==4){
                 hclus_ene_cp4->Fill(cluster->energy());
                 hz_cp4_c->Fill(cluster->z());
                 ClustNo_cp4++;
               }
 
               hz_clust->Fill(cluster->z());
 
               //cout<<"CLUSTER: th="<<cluster_theta<<", ph="<<cluster_phi<<", energy="<<cluster_energy<<endl;
 
               sum_clus_ene+=cluster_energy;
 
               cemc=p4mom.E();
               cemc_sum+=cemc;
 
               hclus_ene_mc->Fill(cemc);
               hclus_ene->Fill(cluster_energy);
 
               hthc_mc->Fill(theta_mc); // Monte Carlo
               hphc_mc->Fill(phi_mc);   // Monte Carlo
 
               hthc->Fill(cluster_theta); // reconstructed
               hphc->Fill(cluster_phi);   // reconstructed
 
               hthc_diff->Fill(cluster_theta-theta_mc);
               hphc_diff->Fill(cluster_phi-phi_mc);
             }
           }
         hclus_ene_sum_mc->Fill(cemc_sum);
 
         hclus_ene_sum->Fill(sum_clus_ene);
 
         hClustNoOrig->Fill(ClustNoOrig);
         hClustNo->Fill(ClustNo);
 
         hClNo_cp1->Fill(ClustNo_cp1);
         hClNo_cp2->Fill(ClustNo_cp2);
         hClNo_cp3->Fill(ClustNo_cp3);
         hClNo_cp4->Fill(ClustNo_cp4);
    }

bool IsRightNeighbour	(	const FTSCANPlet &	a,
		const FTSCANPlet &	b,
		float	pick,
		float &	chi2
	)

inline

Definition at line 1863 of file PndFTSCAGBTracker.cxx.

References fabs(), i, FTSCANPlet::IHit(), FTSCANPlet::ISta(), FTSCANPlet::N(), FTSCANPlet::Param(), FTSCANPlet::QMomentum(), FTSCANPlet::QMomentumErr2(), sqrt(), and StartStationShift.

 {
   int start = (a.N() < b.N() ) ? 0 : a.N() - b.N();
   for( int i = start; i<a.N()-StartStationShift; i++)
   {
     if ( a.IHit(i+StartStationShift) != b.IHit(i) )
     {
       return false;
     }
   }
 
   //chi2 = fabs(a.QMomentum() - b.QMomentum())/sqrt(a.QMomentumErr2() + b.QMomentumErr2());
 
   chi2= fabs(a.Param().Tx() - b.Param().Tx())/sqrt(a.Param().Err2Tx() + b.Param().Err2Tx());
   //cout<<"b.ISta(b.N()-1) "<<b.ISta(b.N()-1)<<endl;
   if (b.ISta(b.N()-1)==26)
   {
     chi2 = fabs(a.QMomentum() - b.QMomentum())/sqrt(a.QMomentumErr2() + b.QMomentumErr2());
   }
   /*
   float chi2_qp = fabs(a.QMomentum() - b.QMomentum())/sqrt(a.Param().Err2QMomentum() + b.Param().Err2QMomentum());
   float chi2_x = fabs(a.Param().X() - b.Param().X())/sqrt(a.Param().Err2X() + b.Param().Err2X());
 
   cout<<"a.QP "<<a.QMomentum()<<" b.QP "<<b.QMomentum()<<endl;
   cout<<"a.C44 "<<a.Param().Err2QMomentum()<<" b.C44 "<<b.Param().Err2QMomentum()<<endl;
   cout<<"a.Tx() "<<a.Param().Tx()<<" b.Tx() "<<b.Param().Tx()<<endl;
   cout<<"a.C22 "<<a.Param().Err2Tx()<<" b.C22 "<<b.Param().Err2Tx()<<endl;
   cout<<"a.X() "<<a.Param().X()<<" b.X() "<<b.Param().X()<<endl;
   cout<<"a.C00 "<<a.Param().Err2X()<<" b.C00 "<<b.Param().Err2X()<<endl;
   cout<<"dif qp "<<a.QMomentum()-b.QMomentum()<<" tx "<<a.Param().Tx()-b.Param().Tx()<<" x "<<a.Param().X()-b.Param().X()<<endl;
   cout<<"chi2_tx "<<chi2<<" chi2_qp "<<chi2_qp<<" chi2_x "<<chi2_x<<endl;
   */
   if ( chi2 > pick ) return false; // neighbours must have same qp (or some other criteria)
 
   chi2 *= chi2;
   //cout<<"picked neigh \n";
   //cout<<endl;
   return true;
 }

t SetQPt ( float_v(1.e-8f) )

Variable Documentation

vector<TrackHitRecord> gTrackHitRecords

Definition at line 1556 of file PndFTSCAGBTracker.cxx.

ok &return ok = active0

Definition at line 820 of file PndFTSCAGBTracker.cxx.

Referenced by PndKFParticleFinder::Exec(), KFParticleFinder::FindParticles(), KFParticleFinder::FindTrackV0Decay(), PndFTSCAGBTracker::FitTrack(), KFParticleBaseSIMD::GetDistanceToVertexLine(), PndFTSCATrackParam::GetXYZPxPyPzQ(), PndFTSTopoReconstructor::Init(), init(), PndRadMapBoxMesh::IsInside(), KFParticleSIMD::KFParticleSIMD(), and PndSoftTriggerTask::ReadConfiguration().

bool SINGLE_THREADED = false

Definition at line 75 of file PndFTSCAGBTracker.cxx.

const int StartStationShift = 1

Definition at line 1412 of file PndFTSCAGBTracker.cxx.

Classes

Functions

Variables

Function Documentation

Variable Documentation