EdbShowPVRQuality Class Reference

#include <EdbShowPVRQuality.h>

Inheritance diagram for EdbShowPVRQuality:

Collaboration diagram for EdbShowPVRQuality:

Public Member Functions
void	CheckEdbPVRec ()
void	CheckFilledXYSize ()
Bool_t	CheckSegmentQualityInPattern_ConstBTDens (EdbPVRec *ali, Int_t PatternAtNr, EdbSegP *seg)
Bool_t	CheckSegmentQualityInPattern_ConstQual (EdbPVRec *ali, Int_t PatternAtNr, EdbSegP *seg)
	ClassDef (EdbShowPVRQuality, 1)
void	CreateEdbPVRec ()
	EdbShowPVRQuality ()
	EdbShowPVRQuality (EdbPVRec *ali)
	EdbShowPVRQuality (EdbPVRec *ali, Float_t BTDensityTargetLevel)
void	Execute_ConstantBTDensity ()
void	Execute_ConstantQuality ()
Int_t	FindFirstBinAbove (TH1 *hist, Double_t threshold, Int_t axis)
Int_t	FindLastBinAbove (TH1 *hist, Double_t threshold, Int_t axis)
Float_t *	GetagreementChi2Cut ()
Float_t	GetagreementChi2Cut (Int_t patNR)
Float_t	GetagreementChi2CutMeanChi2 ()
Float_t	GetagreementChi2CutMeanW ()
Float_t	GetagreementChi2CutRMSChi2 ()
Float_t	GetagreementChi2CutRMSW ()
Float_t	GetBTDensityLevel ()
Bool_t	GetBTDensityLevelCalcMethodMC ()
Int_t	GetBTDensityLevelCalcMethodMCConfirmation ()
Int_t	GetCutMethod ()
Bool_t	GetCutMethodIsDone (Int_t type)
Float_t *	GetCutp0 ()
Float_t	GetCutp0 (Int_t patNR)
Float_t *	GetCutp1 ()
Float_t	GetCutp1 (Int_t patNR)
EdbPVRec *	GetEdbPVRec ()
EdbPVRec *	GetEdbPVRec (Bool_t NeedModified)
EdbPVRec *	GetEdbPVRec (Int_t EdbPVRecType)
EdbPVRec *	GetEdbPVRec_modified ()
EdbPVRec *	GetEdbPVRec_orig ()
TH2F *	GetHistChi2W ()
TH2F *	GetHistYX ()
TObjArray *	GetTracksFromLinkedTracksRootFile ()
void	Help ()
void	Print ()
void	PrintCutType ()
void	PrintCutType0 ()
void	PrintCutType1 ()
EdbPVRec *	Remove_DoubleBT (EdbPVRec *aliSource)
EdbPVRec *	Remove_Passing (EdbPVRec *aliSource)
EdbPVRec *	Remove_Segment (EdbSegP *seg, EdbPVRec *aliSource=NULL, Int_t Option=0)
EdbPVRec *	Remove_SegmentArray (TObjArray *segarray, EdbPVRec *aliSource=NULL, Int_t Option=0)
EdbPVRec *	Remove_Track (EdbTrackP *track, EdbPVRec *aliSource=NULL, Int_t Option=0)
EdbPVRec *	Remove_TrackArray (TObjArray *trackArray, EdbPVRec *aliSource=NULL, Int_t Option=0)
void	SetBTDensityLevel (Float_t BTDensityLevel)
void	SetBTDensityLevelCalcMethodMC (Bool_t BTDensityLevelCalcMethodMC)
void	SetBTDensityLevelCalcMethodMCConfirmation (Int_t BTDensityLevelCalcMethodMCConfirmationNumber)
void	SetCutMethod (Int_t CutMethod)
void	SetEdbPVRec (EdbPVRec *Ali_orig)
void	SetHistGeometry_MC ()
void	SetHistGeometry_OPERA ()
void	SetHistGeometry_OPERAandMC ()
TObjArray *	TrackArrayToSegmentArray (TObjArray *trackArray)
TObjArray *	TrackToSegmentArray (EdbTrackP *track)
virtual	~EdbShowPVRQuality ()

Protected Member Functions
void	Init ()
void	Set0 ()

Private Attributes
Float_t	eAgreementChi2CutMeanChi2
Float_t	eAgreementChi2CutMeanW
Float_t	eAgreementChi2CutRMSChi2
Float_t	eAgreementChi2CutRMSW
Float_t	eAgreementChi2WDistCut [114]
Int_t	eAli_maxNpatterns
EdbPVRec *	eAli_modified
EdbPVRec *	eAli_orig
Float_t	eBTDensityLevel
Bool_t	eBTDensityLevelCalcMethodMC
Int_t	eBTDensityLevelCalcMethodMCConfirmationNumber
Float_t	eCutDistChi2 [114]
Float_t	eCutDistW [114]
Int_t	eCutMethod
Bool_t	eCutMethodIsDone [2]
Float_t	eCutp0 [114]
Float_t	eCutp1 [114]
TH2F *	eHistChi2W
Int_t	eHistGeometry
TH2F *	eHistYX
Bool_t	eIsSource
Bool_t	eIsTarget
Bool_t	eNeedModified
Float_t	ePatternBTDensity_modified [114]
Float_t	ePatternBTDensity_orig [114]
TProfile *	eProfileBTdens_vs_PID_source
Float_t	eProfileBTdens_vs_PID_source_meanX
Float_t	eProfileBTdens_vs_PID_source_meanY
Float_t	eProfileBTdens_vs_PID_source_rmsX
Float_t	eProfileBTdens_vs_PID_source_rmsY
TProfile *	eProfileBTdens_vs_PID_target
Float_t	eProfileBTdens_vs_PID_target_meanX
Float_t	eProfileBTdens_vs_PID_target_meanY
Float_t	eProfileBTdens_vs_PID_target_rmsX
Float_t	eProfileBTdens_vs_PID_target_rmsY
Float_t	maxX
Float_t	maxY
Float_t	minX
Float_t	minY
Int_t	NbinsX
Int_t	NbinsY

Constructor & Destructor Documentation

EdbShowPVRQuality() [1/3]

EdbShowPVRQuality::EdbShowPVRQuality

(

)

EdbShowPVRQuality() [2/3]

EdbShowPVRQuality::EdbShowPVRQuality

(

EdbPVRec *

ali

)

{
    // Default Constructor with a EdbPVRec object.
    // (the EdbPVRec object corresponds to the collection of plates scanned with the collection
    //  of basetracks (and additionally linked tracks)).
    // This constructor automatically checks the original data for background level and
    // creates a new EdbPVRec object that contains only segments that fulfill the quality
    // cut in accordance with the desired (predefined) background level.
    // If general basetrack density is lower than 20BT/mm2 then no cleaning is done.
    cout << "EdbShowPVRQuality::EdbShowPVRQuality(EdbPVRec* ali)   Constructor (does automatically all in one...)"<<endl;
 
    Init();
    Set0();
 
    if (NULL == ali || ali->Npatterns()<1) {
        cout << "WARNING   EdbShowPVRQuality(EdbPVRec* ali)   ali is no good EdbPVRec object! Check!" << endl;
        return;
    }
 
    // Set ali as eAli_orig
    eAli_orig=ali;
    eIsSource=kTRUE;
 
    CheckEdbPVRec();
    Execute_ConstantBTDensity();
    CreateEdbPVRec();
    Print();
}

EdbShowPVRQuality() [3/3]

EdbShowPVRQuality::EdbShowPVRQuality	(	EdbPVRec *	ali,
		Float_t	BTDensityTargetLevel
	)

{
    // Default Constructor with a EdbPVRec object and the desired basetrack density target level.
    // Does same as constructor EdbShowPVRQuality::EdbShowPVRQuality(EdbPVRec* ali) but now with adjustable
    // background level.
    // If general basetrack density is lower than BTDensityTargetLevel BT/mm2 then no cleaning is done.
    cout << "EdbShowPVRQuality::EdbShowPVRQuality(EdbPVRec* ali)   Constructor (does automatically all in one...)"<<endl;
 
    Init();
    Set0();
 
    if (NULL == ali || ali->Npatterns()<1) {
        cout << "WARNING   EdbShowPVRQuality(EdbPVRec* ali)   ali is no good EdbPVRec object! Check!" << endl;
        return;
    }
 
    // Set ali as eAli_orig
    eAli_orig=ali;
    eIsSource=kTRUE;
 
    // Set BTDensityTargetLevel
    SetBTDensityLevel(BTDensityTargetLevel);
    cout << "EdbShowPVRQuality::EdbShowPVRQuality(EdbPVRec* ali)  GetBTDensityLevel() " << GetBTDensityLevel() << endl;
 
    CheckEdbPVRec();
    Execute_ConstantBTDensity();
    CreateEdbPVRec();
    Print();
}

~EdbShowPVRQuality()

EdbShowPVRQuality::~EdbShowPVRQuality ( )

virtual

{
    // Default Destructor
    cout << "EdbShowPVRQuality::~EdbShowPVRQuality()"<<endl;
    delete      eHistChi2W;
    delete      eHistYX;
    delete      eProfileBTdens_vs_PID_source;
}

Member Function Documentation

CheckEdbPVRec()

void EdbShowPVRQuality::CheckEdbPVRec

(

)

{
    // Main function to check if the EdbPVRec object of the scanned data is of low/high background.
    // Following steps are carried out:
    //  Get plate, count number of basetracks in the unit area (1x1cm^2).
    //  Fill (draw if desired (like in EDA display)) histogram with the entries of the unit area.
    //  Get mean of the histogram, compare this value with the reference value.
    // The histogram covers all the area of one emulsion. (for the record: the old ORFEO MC
    // simulation gives not the same position as data does. The area of the histogramm was largely
    // increased to cover both cases).
    // It gives a good estimation of the density. Spikes in some plates, or in some zones are not
    // checked for, this is on the todo list, but maybe not so important.
 
    Log(2,"EdbShowPVRQuality::CheckEdbPVRec","EdbShowPVRQuality::CheckEdbPVRec");
 
    if (!eIsSource) {
        cout << "EdbShowPVRQuality::CheckEdbPVRec  eIsSource=  " << eIsSource << ". This means no source set. Return!" << endl;
        return;
    }
    // Check the patterns of the EdbPVRec:
    eAli_maxNpatterns= eAli_orig->Npatterns();
//     cout << "EdbShowPVRQuality::CheckEdbPVRec  eAli_orig->Npatterns()=  " << eAli_maxNpatterns << endl;
    if (eAli_maxNpatterns>57) cout << " This tells us not yet if we do have one/two brick reconstruction done. A possibility could also be that the dataset was read with microtracks. Further investigation is needed! (On todo list)." << endl;
    if (eAli_maxNpatterns>114) {
        cout << "ERROR! EdbShowPVRQuality::CheckEdbPVRec  eAli_orig->Npatterns()=  " << eAli_maxNpatterns << " is greater than possible basetrack data two bricks. This class does (not yet) work with this large number of patterns. Set maximum patterns to 114!!!." << endl;
        eAli_maxNpatterns=114;
    }
 
    int Npat = eAli_maxNpatterns;
    TH1F* histPatternBTDensity = new TH1F("histPatternBTDensity","histPatternBTDensity",200,0,200);
 
    // Loop over the patterns of the EdbPVRec:
    for (Int_t i=0; i<Npat; i++) {
 
        if (i>56) {
            cout << "WARNING     EdbShowPVRQuality::CheckEdbPVRec()    Your EdbPVRec object has more than 57 patterns! " << endl;
            cout << "WARNING     EdbShowPVRQuality::CheckEdbPVRec()    Check it! " << endl;
        }
        if (gEDBDEBUGLEVEL>2) cout << "CheckEdbPVRec   Doing Pattern " << i << endl;
 
 
        eHistYX->Reset(); // important to clean the histogram
        eHistChi2W->Reset(); // important to clean the histogram
 
        EdbPattern* pat = (EdbPattern*)eAli_orig->GetPattern(i);
        Int_t npat=pat->N();
 
        EdbSegP* seg=0;
        // Loop over the segments of the pattern
        for (Int_t j=0; j<npat; j++) {
            seg=(EdbSegP*)pat->At(j);
            // Very important:
            // For the data case, we assume the following:
            // Data (MCEvt<0) will     be taken for BTdensity calculation
            // Sim (MCEvt>0)      will NOT be taken for BTdensity calculation
            // We take it ONLY and ONLY into account if it is especially wished
            // by the user!
            // Therefore (s)he needs to know how many Gauge Coupling Parameters
            // in the Standard Model exist (at least)...
            Bool_t result=kTRUE;
            if (seg->MCEvt()>0) {
                if (eBTDensityLevelCalcMethodMCConfirmationNumber==18&&eBTDensityLevelCalcMethodMC==kTRUE) {
                    result = kTRUE;
                    // cout << "result = kTRUE !! " << endl;
                }
                else {
                    result = kFALSE;
                }
            }
 
            if (gEDBDEBUGLEVEL>4)  cout << "Doing segment " << j << " result for bool query is: " << result << endl;
 
            // Main decision for segment to be kept or not  (seg is of MC or data type).
            if ( kFALSE == result ) continue;
 
            // For the check, fill the histograms in any case:
            eHistYX->Fill(seg->Y(),seg->X());
            eHistChi2W->Fill(seg->W(),seg->Chi2());
        }
 
        if (gEDBDEBUGLEVEL>2) cout << "I have filled the eHistYX Histogram. Entries = " << eHistYX->GetEntries() << endl;
 
        // Important to reset histogram before it is filled.
        histPatternBTDensity->Reset();
 
        // Search for empty bins, because they can spoil the overall calulation
        // of the mean value.
        Int_t nbins=eHistYX->GetNbinsX()*eHistYX->GetNbinsY();
        Int_t nemptybinsXY=0;
        Int_t bincontentXY=0;
        for (int k=1; k<nbins-1; k++) {
            if (eHistYX->GetBinContent(k)==0) {
                ++nemptybinsXY;
                continue;
            }
            bincontentXY=eHistYX->GetBinContent(k);
            histPatternBTDensity->Fill(bincontentXY);
            eProfileBTdens_vs_PID_source->Fill(i,bincontentXY);
        }
 
        // failsafe warning in case that there are many bins with zero content.
        // for now we print a error message:
        CheckFilledXYSize();
 
        // Save the density in the variable.
        ePatternBTDensity_orig[i]=histPatternBTDensity->GetMean();
 
    }
 
    eProfileBTdens_vs_PID_source_meanX=eProfileBTdens_vs_PID_source->GetMean(1);
    eProfileBTdens_vs_PID_source_meanY=eProfileBTdens_vs_PID_source->GetMean(2);
    eProfileBTdens_vs_PID_source_rmsX=eProfileBTdens_vs_PID_source->GetRMS(1);
    eProfileBTdens_vs_PID_source_rmsY=eProfileBTdens_vs_PID_source->GetRMS(2);
 
    // No assignment for the  eProfileBTdens_vs_PID_target  histogram yet.
    // This will be done in one of the two Execute_ functions.
 
 
    // This will be commented when using in batch mode...
    // For now its there for clarity reasons.
    TCanvas* c1 = new TCanvas();
    c1->Divide(2,2);
    c1->cd(1);
    eHistYX->DrawCopy("colz");
    c1->cd(2);
    eHistChi2W->DrawCopy("colz");
    c1->cd(3);
    histPatternBTDensity->DrawCopy("");
    c1->cd(4);
    eProfileBTdens_vs_PID_source->Draw("profileZ");
    eProfileBTdens_vs_PID_source->GetXaxis()->SetRangeUser(0,eAli_maxNpatterns+2);
    c1->cd();
 
    eHistYX->Reset();
    eHistChi2W->Reset();
 
    Log(2,"EdbShowPVRQuality::CheckEdbPVRec","EdbShowPVRQuality::CheckEdbPVRec...done");
    return;
}

CheckFilledXYSize()

void EdbShowPVRQuality::CheckFilledXYSize

(

)

{
    // Check the bins filled in the actual pattern.
    // The histogram of the XY distribution is analysed.
    // A warning is given if more than 10 percent of
    // the bins are empty.
    // In this case one should look closer at the specific
    // plate distribution.
 
    if (gEDBDEBUGLEVEL>3)   cout << "-----     void EdbShowPVRQuality::CheckFilledXYSize() return maximum/minimum entries of  eHistYX    -----" << endl;
    Int_t nbx,nby=0;
    nbx=eHistYX->GetNbinsX();
    nby=eHistYX->GetNbinsY();
 
    Int_t n1x= FindFirstBinAbove(eHistYX,0,1);
    Int_t n1y= FindFirstBinAbove(eHistYX,0,2);
    Int_t n2x= FindLastBinAbove(eHistYX,0,1);
    Int_t n2y= FindLastBinAbove(eHistYX,0,2);
    Int_t width_x=0;
    Int_t width_y=0;
    width_x=TMath::Abs(eHistYX->GetXaxis()->GetBinCenter(n1x)-eHistYX->GetXaxis()->GetBinCenter(n2x));
    width_y=TMath::Abs(eHistYX->GetYaxis()->GetBinCenter(n1y)-eHistYX->GetYaxis()->GetBinCenter(n2y));
 
    // Now check the number of empty bins between! the filled area
    // within (FindFirstBinAbove,FindLastBinAbove)
    // This function is NOT optimized for speed :-)
    Int_t NonEmptyBins=0;
    Int_t nBins=0;
    for (Int_t i=n1x; i<n2x; ++i) {
        for (Int_t j=n1y; j<n2y; ++j) {
            ++nBins;
            if (eHistYX->GetBinContent(i,j)==0) continue;
            if (eHistYX->GetBinContent(i,j)==0) continue;
            ++NonEmptyBins;
        }
    }
    Float_t FractionOfEmptyBins=1-(Float_t(NonEmptyBins)/Float_t(nBins));
    if (FractionOfEmptyBins>0.1) cout << "WARNING: void EdbShowPVRQuality::CheckFilledXYSize() FractionOfEmptyBins = " << FractionOfEmptyBins << endl;
 
    if (gEDBDEBUGLEVEL>3) {
        cout << "----      NonEmptyBins bins in covered area:  = " << NonEmptyBins << endl;
        cout << "----      nBins totale bins in covered area:  = " << nBins << endl;
        cout << "----      Extrem Center Endpoints of eHistYX --- " << endl;
        cout << "----      nbxMin= " << n1x << endl;
        cout << "----      nbxMax= " << n1y << endl;
        cout << "----      nbyMin= " << n2x << endl;
        cout << "----      nbyMax= " << n2y << endl;
        cout << "----      nbxMin= " << eHistYX->GetXaxis()->GetBinCenter(n1x) << endl;
        cout << "----      nbxMax= " << eHistYX->GetYaxis()->GetBinCenter(n1y) << endl;
        cout << "----      nbxMin= " << eHistYX->GetXaxis()->GetBinCenter(n2x) << endl;
        cout << "----      nbxMax= " << eHistYX->GetYaxis()->GetBinCenter(n2y) << endl;
        cout << "-----     void EdbShowPVRQuality::CheckFilledXYSize() ... done. " << endl;
    }
    return;
}

CheckSegmentQualityInPattern_ConstBTDens()

Bool_t EdbShowPVRQuality::CheckSegmentQualityInPattern_ConstBTDens	(	EdbPVRec *	ali,
		Int_t	PatternAtNr,
		EdbSegP *	seg
	)

{
    // Core function to check if a basetrack of the specific pattern matches the expectations
    // for the desired quality cut (calculated on the estimations in CheckEdbPVRec(). ).
    // Implementation for the Cuttype 0: Constant BT Density
    // ---
    // Note: since the eCutp1[i] values are calculated with
    // this pattern->At() scheme labelling,
    // its not necessarily guaranteed that seg->PID gives correct this
    // number back. Thats why we have to give the PatternAtNr here again.
    //
    // And: it is not checked here if seg is contained in this specific
    // pattern. It looks only for the quality cut!
    if (gEDBDEBUGLEVEL>3)  cout << "seg->W()* eCutp1[PatternAtNr] - eCutp0[PatternAtNr] = " << seg->W()* eCutp1[PatternAtNr] - eCutp0[PatternAtNr] << endl;
 
    // Constant BT density cut:
    if (seg->Chi2() >= seg->W()* eCutp1[PatternAtNr] - eCutp0[PatternAtNr]) return kFALSE;
 
    if (gEDBDEBUGLEVEL>3) cout <<"EdbShowPVRQuality::CheckSegmentQualityInPattern_ConstBTDens()   Segment " << seg << " has passed ConstBTDens cut!" << endl;
    return kTRUE;
}

CheckSegmentQualityInPattern_ConstQual()

Bool_t EdbShowPVRQuality::CheckSegmentQualityInPattern_ConstQual	(	EdbPVRec *	ali,
		Int_t	PatternAtNr,
		EdbSegP *	seg
	)

{
    // Core function to check if a basetrack of the specific pattern matches the expectations
    // for the desired quality cut (calculated on the estimations in CheckEdbPVRec(). ).
    // Implementation for the Cuttype 1: Constant Quality.
    // ---
    // See comments in CheckSegmentQualityInPattern_ConstBTDens
    // Constant BT quality cut:
    Float_t agreementChi2=TMath::Sqrt( ( (seg->Chi2()-eAgreementChi2CutMeanChi2)/eAgreementChi2CutRMSChi2)*((seg->Chi2()-eAgreementChi2CutMeanChi2)/eAgreementChi2CutRMSChi2)  +   ((seg->W()-eAgreementChi2CutMeanW)/eAgreementChi2CutRMSW)*((seg->W()-eAgreementChi2CutMeanW)/eAgreementChi2CutRMSW) );
    if (agreementChi2>eAgreementChi2WDistCut[PatternAtNr]) return kFALSE;
 
    if (gEDBDEBUGLEVEL>3) cout <<"EdbShowPVRQuality::CheckSegmentQualityInPattern_ConstQual()   Segment " << seg << " has passed ConstQual cut!" << endl;
    return kTRUE;
}

ClassDef()

EdbShowPVRQuality::ClassDef	(	EdbShowPVRQuality	,
		1
	)

CreateEdbPVRec()

void EdbShowPVRQuality::CreateEdbPVRec

(

)

{
    // Creates the cleaned EdbPVRec object, containing only segments that
    // satisfied the cutcriteria.
    // Attention: the couples structure and the tracking structure will be lost,
    // this EdbPVRec object is only useful for the list of segments (shower reco).
    // DO NOT USE THIS ROUTINE FOR GENERAL I/O and/or GENERAL EdbPVRec operations!
 
    cout << "-----     void EdbShowPVRQuality::CreateEdbPVRec()     -----" << endl;
    if (gEDBDEBUGLEVEL>2) {
        cout << "-----     " << endl;
        cout << "-----     This function makes out of the original eAli" << endl;
        cout << "-----     a new EdbPVRec object having only those seg-" << endl;
        cout << "-----     ments in it which satisfy the cutcriteria " << endl;
        cout << "-----     determined in Execute_ConstantBTDensity, Execute_ConstantQuality" << endl;
        cout << "-----     " << endl;
        cout << "-----     WARNING: the couples structure and the tracking structure" << endl;
        cout << "-----     will be lost, this PVR object is only useful for the" << endl;
        cout << "-----     list of Segments (==ShowReco...) ... " << endl;
        cout << "-----     DO NOT USE THIS ROUTINE FOR GENERAL I/O and/or EdbPVRec operations!" << endl;
        cout << "-----     " << endl;
        cout << "CreateEdbPVRec()  Mode 0:" << eCutMethodIsDone[0] << endl;
        cout << "CreateEdbPVRec()  Mode 1:" << eCutMethodIsDone[1] << endl;
        cout << "-----     " << endl;
        cout << "-----     ----------------------------------------------" << endl;
    }
 
    // Checks
    if (NULL==eAli_orig || eIsSource==kFALSE) {
        cout << "WARNING!   EdbShowPVRQuality::CreateEdbPVRec    NULL==eAli_orig   || eIsSource==kFALSE    return."<<endl;
        return;
    }
    // Checks
    if (eCutMethodIsDone[0]==kFALSE && eCutMethodIsDone[1]==kFALSE) {
        cout << "WARNING!   EdbShowPVRQuality::CreateEdbPVRec    eCutMethodIsDone[0]==kFALSE && eCutMethodIsDone[1]==kFALSE   return."<<endl;
        return;
    }
 
    // Make a new PVRec object anyway
    eAli_modified = new EdbPVRec();
 
    // These two lines dont compile yet ... (???) ...
    //  EdbScanCond* scancond = eAli_orig->GetScanCond();
    //  eAli_modified->SetScanCond(*scancond);
 
    Float_t agreementChi2;
 
    // This makes pointer copies of patterns with segments list.
    // wARNING: the couples structure and the tracking structure
    // will be lost, this PVR object is only useful for the
    // list of Segments (==ShowReco...) ...
 
    // Priority has the Execute_ConstantQuality cut.
    // If this was done we take this...:
 
    // Loop over patterns
    for (int i = 0; i <eAli_orig->Npatterns(); i++ ) {
        EdbPattern* pat = eAli_orig->GetPattern(i);
        EdbPattern* pt= new EdbPattern();
        // SetPattern Values to the parent patterns:
        pt->SetID(pat->ID());
        pt->SetPID(pat->PID());
        pt->SetZ(pat->Z());
 
        // Loop over segments
        for (int j = 0; j <pat->N(); j++ ) {
            EdbSegP* seg = pat->GetSegment(j);
 
            // Put here the cut condition ...
            if (eCutMethodIsDone[0]==kTRUE && eCutMethodIsDone[1]==kFALSE) {
                // Constant BT density cut:
                if (seg->Chi2() >= seg->W()* eCutp1[i] - eCutp0[i]) continue;
            }
            else if (eCutMethodIsDone[1]==kTRUE) {
                // Constant Quality cut:
                agreementChi2=TMath::Sqrt( ( (seg->Chi2()-eAgreementChi2CutMeanChi2)/eAgreementChi2CutRMSChi2)*((seg->Chi2()-eAgreementChi2CutMeanChi2)/eAgreementChi2CutRMSChi2)  +   ((seg->W()-eAgreementChi2CutMeanW)/eAgreementChi2CutRMSW)*((seg->W()-eAgreementChi2CutMeanW)/eAgreementChi2CutRMSW) );
                if (agreementChi2>eAgreementChi2WDistCut[i]) continue;
            }
            else {
                // do nothing;
            }
 
            // Add segment:
            pt->AddSegment(*seg);
        }
        eAli_modified->AddPattern(pt);
    }
 
    eIsTarget=kTRUE;
 
    //---------------------
    cout << "-----     void EdbShowPVRQuality::CreateEdbPVRec()...Created new EdbPVR object at address " << eAli_modified << endl;
    //---------------------
    cout << "-----     void EdbShowPVRQuality::CreateEdbPVRec()...done." << endl;
    return;
}

Execute_ConstantBTDensity()

void EdbShowPVRQuality::Execute_ConstantBTDensity

(

)

{
    // Execute the modified cut routines to achieve the basetrack density level,
    // after application the specific cut on the segments of the specific plate (pattern).
    // The Constant BT Density is defined by the number of BT/mm2 in the histogram.
 
    for (int i=0; i<80; ++i) cout << "-";
    cout << endl;
    cout << "-";
    cout << endl;
    cout << "EdbShowPVRQuality::Execute_ConstantBTDensity " << endl;
    cout << "-";
    cout << endl;
    for (int i=0; i<80; ++i) cout << "-";
    cout << endl;
 
    if (!eIsSource) return;
    if (NULL==eAli_orig) return;
 
    // Check the patterns of the EdbPVRec:
    eAli_maxNpatterns= eAli_orig->Npatterns();
    if (eAli_maxNpatterns>57) cout << " This tells us not yet if we do have one/two brick reconstruction done. A possibility could also be that the dataset was read with microtracks. Further investigation is needed! (On todo list)." << endl;
    if (eAli_maxNpatterns>114) {
        cout << "ERROR!  EdbShowPVRQuality::Execute_ConstantBTDensity  eAli_orig->Npatterns()=  " << eAli_maxNpatterns << " is greater than possible basetrack data two bricks. This class does (not yet) work with this large number of patterns. Set maximum patterns to 114!!!." << endl;
        eAli_maxNpatterns=114;
    }
 
    TH1F* histPatternBTDensity = new TH1F("histPatternBTDensity","histPatternBTDensity",200,0,200);
 
    // Loop over the patterns:
    for (int i=0; i<eAli_maxNpatterns; i++) {
        if (i>56) {
            cout << "WARNING     EdbShowPVRQuality::Execute_ConstantBTDensity()    Your EdbPVRec object has more than 57 patterns! " << endl;
            cout << "WARNING     EdbShowPVRQuality::Execute_ConstantBTDensity()    Check it! " << endl;
        }
        if (gEDBDEBUGLEVEL>2) cout << "Execute_ConstantBTDensity   Doing Pattern " << i << endl;
 
        // Now the condition loop:
        // Loop over 20 steps a 0.15,0.145,0.14 ...  down to 0.07
        for (int l=0; l<20; l++) {
            eHistYX->Reset(); // important to clean the histogram
            eHistChi2W->Reset(); // important to clean the histogram
            histPatternBTDensity->Reset(); // important to clean the histogram
 
            EdbPattern* pat = (EdbPattern*)eAli_orig->GetPattern(i);
            Int_t npat=pat->N();
            EdbSegP* seg=0;
            // Loop over the segments.
            for (int j=0; j<npat; j++) {
                seg=(EdbSegP*)pat->At(j);
                // Very important:
                // For the data case, we assume the following:
                // Data (MCEvt<0) will     be taken for BTdensity calculation
                // Sim (MCEvt>0)      will NOT be taken for BTdensity calculation
                // We take it ONLY and ONLY into account if it is especially wished
                // by the user!
                // Therefore (s)he needs to know how many Gauge Coupling Parameters
                // in the Standard Model exist (at least)...
                Bool_t result=kTRUE;
                if (seg->MCEvt()>0) {
                    if (eBTDensityLevelCalcMethodMCConfirmationNumber==18&&eBTDensityLevelCalcMethodMC==kTRUE) {
                        result = kTRUE;
                        // cout << "result = kTRUE !! " << endl;
                    }
                    else {
                        result = kFALSE;
                    }
                }
 
                if (gEDBDEBUGLEVEL>4)  cout << "Doing segment " << j << " result for bool query is: " << result << endl;
 
                // Main decision for segment to be kept or not (seg is of MC or data type).
                if ( kFALSE == result ) continue;
                // Constant BT density cut:
                if (seg->Chi2() >= seg->W()* eCutp1[i] - eCutp0[i]) continue;
 
                eHistYX->Fill(seg->Y(),seg->X());
                eHistChi2W->Fill(seg->W(),seg->Chi2());
            }
 
            if (gEDBDEBUGLEVEL>2) cout << "I have filled the eHistYX Histogram. Entries = " << eHistYX->GetEntries() << endl;
 
            int nbins=eHistYX->GetNbinsX()*eHistYX->GetNbinsY();
            for (int k=1; k<nbins-1; k++) {
                if (eHistYX->GetBinContent(k)==0) continue;
                histPatternBTDensity->Fill(eHistYX->GetBinContent(k));
            }
 
            ePatternBTDensity_modified[i]=histPatternBTDensity->GetMean();
            if (gEDBDEBUGLEVEL>2) cout <<"Execute_ConstantBTDensity      Loop l= " << l << ":  for the eCutp1[i] : " << eCutp1[i] <<   "  we have a dens: "  << ePatternBTDensity_modified[i] << endl;
 
            // Now the condition check:
            if (ePatternBTDensity_modified[i]<=eBTDensityLevel) {
                if (gEDBDEBUGLEVEL>2)  cout << "Execute_ConstantBTDensity      We reached the loop end due to good BT density level ... and break loop." << endl;
                break;
            }
            else {
                // Next step, tighten cut:
                eCutp1[i] += -0.005;
            }
 
        } // of condition loop...
 
        // Fill target profile histogram:
        Float_t bincontentXY=histPatternBTDensity->GetMean();
        eProfileBTdens_vs_PID_target->Fill(i,bincontentXY);
 
    } // of Npattern loops..
 
    eCutMethodIsDone[0]=kTRUE;
 
    // Check if modified or original PVRec object should be returned:
    if (eProfileBTdens_vs_PID_source_meanY>eBTDensityLevel) {
        eNeedModified=kTRUE;
        cout << "----------void EdbShowPVRQuality::Execute_ConstantBTDensity() Check if modified or original PVRec object should be returned: " << endl;
        cout << "----------void EdbShowPVRQuality::Execute_ConstantBTDensity() " << eNeedModified << endl;
    }
 
 
    // This will be commented when using in batch mode...
    // For now its there for clarity reasons.
    TCanvas* c1 = new TCanvas();
    c1->Divide(2,2);
    c1->cd(1);
    eHistYX->DrawCopy("colz");
    c1->cd(2);
    eHistChi2W->DrawCopy("colz");
    c1->cd(3);
    histPatternBTDensity->DrawCopy("");
    c1->cd(4);
    eProfileBTdens_vs_PID_source->Draw("profileZ");
    eProfileBTdens_vs_PID_source->GetXaxis()->SetRangeUser(0,eAli_maxNpatterns+2);
    eProfileBTdens_vs_PID_target->SetLineStyle(2);
    eProfileBTdens_vs_PID_target->Draw("profileZsame");
    c1->cd();
    histPatternBTDensity->Reset();
    eHistYX->Reset();
    eHistChi2W->Reset();
 
 
    eProfileBTdens_vs_PID_source_meanX=eProfileBTdens_vs_PID_source->GetMean(1);
    eProfileBTdens_vs_PID_source_meanY=eProfileBTdens_vs_PID_source->GetMean(2);
    eProfileBTdens_vs_PID_source_rmsX=eProfileBTdens_vs_PID_source->GetRMS(1);
    eProfileBTdens_vs_PID_source_rmsY=eProfileBTdens_vs_PID_source->GetRMS(2);
 
    eProfileBTdens_vs_PID_target_meanX=eProfileBTdens_vs_PID_target->GetMean(1);
    eProfileBTdens_vs_PID_target_meanY=eProfileBTdens_vs_PID_target->GetMean(2);
    eProfileBTdens_vs_PID_target_rmsX=eProfileBTdens_vs_PID_target->GetRMS(1);
    eProfileBTdens_vs_PID_target_rmsY=eProfileBTdens_vs_PID_target->GetRMS(2);
 
    delete histPatternBTDensity;
 
    cout << "----------void EdbShowPVRQuality::Execute_ConstantBTDensity() Cuts are done and saved to obtain desired BT density. " << endl;
    cout << "----------void EdbShowPVRQuality::Execute_ConstantBTDensity() If you want to apply the cuts now, run the  CreateEdbPVRec()  function now.  " <<   endl;
    cout << "----------void EdbShowPVRQuality::Execute_ConstantBTDensity()....done." << endl;
    return;
}

Execute_ConstantQuality()

void EdbShowPVRQuality::Execute_ConstantQuality

(

)

______ now same code as in the function Execute_ConstantBTDensity ___________________

{
    // The cut method that compares the passing muon tracks with all scanned segments.
    // This may help to improve, since it takes into account the actual segment quality,
    // which varies from scan to scan anyway.
    // Works for tracks passing the volume to extract their mean chi2/W. Then a distance
    // measurement Sqrt{0.5*[ ((chi2-chi2mean)/chi2rms)^2 + ((W-Wmean)/Wrms)^2 )] } is
    // build and the value of single basetracks is compared to this variable. Starting from
    // dist_max = 3 going down until desired target level is achieved.
    // If eAli->Tracks is there we take them from there.
    // If not, we try if there is a file linked_tracks.root, and we take tracks from there.
    // If that doesnt work either, nothing is done.
 
    if (!eIsSource) return;
 
    for (int i=0; i<80; ++i) cout << "-";
    cout << endl;
    cout << "-";
    cout << endl;
    cout << "EdbShowPVRQuality::Execute_ConstantBTDensity " << endl;
    cout << "-";
    cout << endl;
    for (int i=0; i<80; ++i) cout << "-";
    cout << endl;
    cout << "----------    Works for tracks passing the volume to extract their mean chi2/W " << endl;
    cout << "----------    If eAli->Tracks is there we take them from there." << endl;
    cout << "----------    If not, we try if there is a file linked_tracks.root " << endl;
    cout << "----------    we take tracks from there. " << endl;
    cout << "----------    If that doesnt work either, nothing is done." << endl;
    for (int i=0; i<80; ++i) cout << "-";
    cout << endl;
    cout << "-";
    cout << endl;
 
    cout << "EdbShowPVRQuality::Execute_ConstantQuality()   eAli_orig.eTracks :" << eAli_orig->eTracks << endl;
 
    Float_t meanChi2=0.5;
    Float_t rmsChi2=0.2;
    Float_t meanW=22;
    Float_t rmsW=4;
    Float_t agreementChi2=100;
 
    // No  eAli.Tracks ? Look for tracks in linked_track.root:
    if (NULL == eAli_orig->eTracks) {
        cout << "EdbShowPVRQuality::Execute_ConstantQuality()   No eAli.Tracks. Look for tracks in linked_track.root" << endl;
        TFile* trackfile = new TFile("linked_tracks.root");
        trackfile->ls();
        TTree* tracks = (TTree*)trackfile->Get("tracks");
        if (NULL == tracks) {
            cout << "EdbShowPVRQuality::Execute_ConstantQuality()   No tracks in linked_track.root file. Return, leave eAli_orig unchanged and dont do any cleaning. You might try Execute_ConstantBTDensity instead. " << endl;
            return;
        }
        //      TH1F* h1; Attention the usage of  _npl_ is disvavoured
        //      because this can cause problems.
        tracks->Draw("nseg>>h(60,0,60)","","");
        TH1F *h1 = (TH1F*)gPad->GetPrimitive("h");
 
        // Short implementation of getting the last bin filled:
        int lastfilledbin=0;
        for (int k=1; k<h1->GetNbinsX()-1; k++) {
            if (h1->GetBinContent(k)>0) lastfilledbin=k;
        }
 
        Float_t travellenghtpercentage=Float_t(lastfilledbin)/Float_t(eAli_maxNpatterns);
 
        TString cutstring = TString(Form("nseg>=%d",int(h1->GetBinCenter(lastfilledbin-3)) ));
        tracks->Draw("s.eChi2>>hChi2(100,0,2)",cutstring);
        TH1F *hChi2 = (TH1F*)gPad->GetPrimitive("hChi2");
 
        cout << "EdbShowPVRQuality::Execute_ConstantQuality()   Found tracks in the file! Good. Info:" << endl;
        cout << "EdbShowPVRQuality::Execute_ConstantQuality()   The long tracks in the file have an average percentage length of:" << travellenghtpercentage  <<  endl;
        cout << "EdbShowPVRQuality::Execute_ConstantQuality()   Mean(RMS) of Chi2 distribution of passing tracks: " << hChi2->GetMean() << "+-"  << hChi2->GetRMS() << endl;
 
        TCanvas* c1 = new TCanvas();
        c1->cd();
        tracks->Draw("s.eW>>hW(50,0,50)",cutstring);
        TH1F *hW = (TH1F*)gPad->GetPrimitive("hW");
        cout << "EdbShowPVRQuality::Execute_ConstantQuality()   Mean(RMS) of W distribution of passing tracks: " << hW->GetMean() << "+-"  << hW->GetRMS() << endl;
 
        meanChi2=hChi2->GetMean();
        rmsChi2=hChi2->GetRMS();
        meanW=hW->GetMean();
        rmsW=hW->GetRMS();
 
        // Quick check if these values are reasonable:
        if (TMath::Abs(meanChi2-0.5)>1 || TMath::Abs(meanW-22)>6) {
            cout << "WARNING   EdbShowPVRQuality::Execute_ConstantQuality()   The tracks might have a strange distribution. You are urgently requested to check these manually!!!"<<endl;
        }
        else {
            cout << "EdbShowPVRQuality::Execute_ConstantQuality()   The tracks have reasonable values."<<endl;
        }
 
 
        // since the values for the passing tracks are
        // calculated for the whole volume, we assume that the cutvalues
        // are valid for all plates.
        eAgreementChi2CutMeanChi2=meanChi2;
        eAgreementChi2CutRMSChi2=rmsChi2;
        eAgreementChi2CutMeanW=meanW;
        eAgreementChi2CutRMSW=rmsW;
    }
 
 
    // Check the patterns of the EdbPVRec:
    eAli_maxNpatterns= eAli_orig->Npatterns();
    cout << "EdbShowPVRQuality::Execute_ConstantQuality  eAli_orig->Npatterns()=  " << eAli_maxNpatterns << endl;
    if (eAli_maxNpatterns>57) cout << " This tells us not yet if we do have one/two brick reconstruction done. A possibility could also be that the dataset was read with microtracks. Further investigation is needed! (On todo list)." << endl;
    if (eAli_maxNpatterns>114) {
        cout << "ERROR!  EdbShowPVRQuality::Execute_ConstantQuality  eAli_orig->Npatterns()=  " << eAli_maxNpatterns << " is greater than possible basetrack data two bricks. This class does (not yet) work with this large number of patterns. Set maximum patterns to 114!!!." << endl;
        eAli_maxNpatterns=114;
    }
 
    TH1F* histPatternBTDensity = new TH1F("histPatternBTDensity","histPatternBTDensity",100,0,200);
    TH1F* histagreementChi2 = new TH1F("histagreementChi2","histagreementChi2",100,0,5);
 
 
    cout << "Execute_ConstantQuality   Loop over the patterns..." << endl;
    for (int i=0; i<eAli_maxNpatterns; i++) {
        if (i>56) {
            cout << "ERROR     EdbShowPVRQuality::Execute_ConstantQuality() Your EdbPVRec object has more than 57 patterns! " << endl;
            cout << "ERROR     EdbShowPVRQuality::Execute_ConstantQuality() Check it! " << endl;
        }
 
        if (gEDBDEBUGLEVEL>2) cout << "Execute_ConstantQuality   Doing Pattern " << i << endl;
 
        // Now the condition loop:
        // Loop over 30 steps agreementChi2 step 0.05
        for (int l=0; l<30; l++) {
 
            if (gEDBDEBUGLEVEL>2) cout << "Execute_ConstantQuality   Doing condition loop = " << l << endl;
 
            eHistYX->Reset(); // important to clean the histogram
            eHistChi2W->Reset(); // important to clean the histogram
            histPatternBTDensity->Reset(); // important to clean the histogram
 
            EdbPattern* pat = (EdbPattern*)eAli_orig->GetPattern(i);
            Int_t npat=pat->N();
            EdbSegP* seg=0;
 
            if (gEDBDEBUGLEVEL>2) cout << "Execute_ConstantQuality   Loop over segments of pattern " << i << ",Number of segments= " << npat <<  endl;
            for (int j=0; j<npat; j++) {
                seg=(EdbSegP*)pat->At(j);
 
                if (gEDBDEBUGLEVEL>4) cout << "Execute_ConstantQuality   Doing segment= " << j << endl;
                // Very important:
                // For the data case, we assume the following:
                // Data (MCEvt<0) will     be taken for BTdensity calculation
                // Sim (MCEvt>0)      will NOT be taken for BTdensity calculation
                // We take it ONLY and ONLY into account if it is especially wished
                // by the user!
                // Therefore (s)he needs to know how many Gauge Coupling Parameters
                // in the Standard Model exist (at least)...
                Bool_t result=kTRUE;
                if (seg->MCEvt()>0) {
                    if (eBTDensityLevelCalcMethodMCConfirmationNumber==18&&eBTDensityLevelCalcMethodMC==kTRUE) {
                        result = kTRUE;
                        // cout << "result = kTRUE !! " << endl;
                    }
                    else {
                        result = kFALSE;
                    }
                }
 
                if (gEDBDEBUGLEVEL>4)  cout << "Doing segment " << j << " result for bool query is: " << result << endl;
 
                // Main decision for segment to be kept or not (seg is of MC or data type).
                if ( kFALSE == result ) continue;
 
                // Change here to the quality with values obtained from the tracks.
                // Constant BT quality cut:
                agreementChi2=TMath::Sqrt(0.5 *( (seg->Chi2()-meanChi2)/rmsChi2)*((seg->Chi2()-meanChi2)/rmsChi2)  +   ((seg->W()-meanW)/rmsW)*((seg->W()-meanW)/rmsW) );
                histagreementChi2->Fill(agreementChi2);
 
                // Constant BT quality cut:
                if (agreementChi2>eAgreementChi2WDistCut[i]) continue;
 
                eHistYX->Fill(seg->Y(),seg->X());
                eHistChi2W->Fill(seg->W(),seg->Chi2());
            }
 
            if (gEDBDEBUGLEVEL>2) cout << "I have filled the eHistYX Histogram. Entries = " << eHistYX->GetEntries() << endl;
 
            Int_t nbins=eHistYX->GetNbinsX()*eHistYX->GetNbinsY();
            Int_t nemptybinsXY=0;
            for (int k=1; k<nbins-1; k++) {
                if (eHistYX->GetBinContent(k)==0) {
                    ++nemptybinsXY;
                    continue;
                }
                histPatternBTDensity->Fill(eHistYX->GetBinContent(k));
            }
 
            // failsafe warning in case that there are many bins with zero content.
            // for now we print a error message:
            CheckFilledXYSize();
 
            ePatternBTDensity_modified[i]=histPatternBTDensity->GetMean();
 
            if (gEDBDEBUGLEVEL>2) cout <<"Execute_ConstantQuality      Loop l= " << l << ":  for the eAgreementChi2WDistCut : " << eAgreementChi2WDistCut[i] <<   "  we have a dens: "  << ePatternBTDensity_modified[i] << endl;
 
            // Now the condition check:
            if (ePatternBTDensity_modified[i]<=eBTDensityLevel) {
                if (gEDBDEBUGLEVEL>2)  cout << "Execute_ConstantQuality      We reached the loop end due to good BT density level ... and break loop." << endl;
                // But dont forget to set values:
                eCutDistChi2[i]=meanChi2;
                eCutDistW[i]=meanW;
                eAgreementChi2CutMeanChi2=meanChi2;
                eAgreementChi2CutRMSChi2=rmsChi2;
                eAgreementChi2CutMeanW=meanW;
                eAgreementChi2CutRMSW=rmsW;
                break;
            }
            else {
                // Next step, tighten cut:
                eAgreementChi2WDistCut[i]+=  -0.05;
            }
 
        } // of condition loop...
 
        // Fill target profile histogram:
        Float_t bincontentXY=histPatternBTDensity->GetMean();
        cout << "Execute_ConstantQuality   histPatternBTDensity->GetMean()." << histPatternBTDensity->GetMean() <<  endl;
        eProfileBTdens_vs_PID_target->Fill(i,bincontentXY);
 
    } // of Npattern loops..
    cout << "Execute_ConstantQuality   Loop over the patterns...done." << endl;
 
    // Check if modified or original PVRec object should be returned:
    if (eProfileBTdens_vs_PID_source_meanY>eBTDensityLevel) {
        eNeedModified=kTRUE;
        cout << "----------void EdbShowPVRQuality::Execute_ConstantBTDensity() Check if modified or original PVRec object should be returned: " << endl;
        cout << "----------void EdbShowPVRQuality::Execute_ConstantBTDensity() " << eNeedModified << endl;
    }
 
    eCutMethodIsDone[1]=kTRUE;
 
    // This will be commented when using in batch mode...
    // For now its there for clarity reasons.
    TCanvas* c1 = new TCanvas();
    c1->Divide(2,2);
    c1->cd(1);
    eHistYX->DrawCopy("colz");
    c1->cd(2);
    eHistChi2W->DrawCopy("colz");
    c1->cd(3);
    histPatternBTDensity->DrawCopy("");
    c1->cd(4);
    eProfileBTdens_vs_PID_source->Draw("profileZ");
    eProfileBTdens_vs_PID_source->GetXaxis()->SetRangeUser(0,eAli_maxNpatterns+2);
    eProfileBTdens_vs_PID_target->SetLineStyle(2);
    eProfileBTdens_vs_PID_target->Draw("profileZsame");
    c1->cd();
 
    eProfileBTdens_vs_PID_source_meanX=eProfileBTdens_vs_PID_source->GetMean(1);
    eProfileBTdens_vs_PID_source_meanY=eProfileBTdens_vs_PID_source->GetMean(2);
    eProfileBTdens_vs_PID_source_rmsX=eProfileBTdens_vs_PID_source->GetRMS(1);
    eProfileBTdens_vs_PID_source_rmsY=eProfileBTdens_vs_PID_source->GetRMS(2);
 
    eProfileBTdens_vs_PID_target_meanX=eProfileBTdens_vs_PID_target->GetMean(1);
    eProfileBTdens_vs_PID_target_meanY=eProfileBTdens_vs_PID_target->GetMean(2);
    eProfileBTdens_vs_PID_target_rmsX=eProfileBTdens_vs_PID_target->GetRMS(1);
    eProfileBTdens_vs_PID_target_rmsY=eProfileBTdens_vs_PID_target->GetRMS(2);
 
    histPatternBTDensity->Reset();
    eHistYX->Reset();
    eHistChi2W->Reset();
    delete histPatternBTDensity;
 
 
    cout << "----------void EdbShowPVRQuality::Execute_ConstantQuality() Cuts are done and saved to obtain desired BT density. " << endl;
    cout << "----------void EdbShowPVRQuality::Execute_ConstantQuality() If you want to apply the cuts now, run the  CreateEdbPVRec()  function now.  " <<   endl;
 
    cout << "----------void EdbShowPVRQuality::Execute_ConstantQuality()....done." << endl;
    return;
}

FindFirstBinAbove()

Int_t EdbShowPVRQuality::FindFirstBinAbove	(	TH1 *	hist,
		Double_t	threshold,
		Int_t	axis
	)

                                                                                    {
    // The TH1 function  FindFirstBinAbove  is only implemented in
    // root version >= 5.24. But since many scanning labs use old root
    // versions persistently, I had to copy the TH1 function as a
    // memberfunction of EdbShowPVRQuality
    // Code taken from
    // http://root.cern.ch/root/html/src/TH1.cxx.html#biA7FC
    TAxis* ax=0;
    if (axis==1) ax = hist->GetXaxis();
    if (axis==2) ax = hist->GetYaxis();
    if (axis==3) ax = hist->GetZaxis();
    int nb = ax->GetNbins();
    for (Int_t i=0; i<=nb; i++) {
        if (hist->GetBinContent(i)>threshold) return i;
    }
    return 0;
}

FindLastBinAbove()

Int_t EdbShowPVRQuality::FindLastBinAbove	(	TH1 *	hist,
		Double_t	threshold,
		Int_t	axis
	)

                                                                                   {
    // The TH1 function  FindFirstBinAbove  is only implemented in
    // root version >= 5.24. But since many scanning labs use old root
    // versions persistently, I had to copy the TH1 function as a
    // memberfunction of EdbShowPVRQuality
    // Code taken from
    // http://root.cern.ch/root/html/src/TH1.cxx.html#biA7FC
    TAxis* ax=0;
    if (axis==1) ax = hist->GetXaxis();
    if (axis==2) ax = hist->GetYaxis();
    if (axis==3) ax = hist->GetZaxis();
    int nb = ax->GetNbins();
    for (Int_t i=nb; i>=1; i--) {
        if (hist->GetBinContent(i)>threshold) return i;
    }
    return 0;
}

GetagreementChi2Cut() [1/2]

Float_t * EdbShowPVRQuality::GetagreementChi2Cut ( )

inline

                                            {
        return eAgreementChi2WDistCut;
    }

GetagreementChi2Cut() [2/2]

Float_t EdbShowPVRQuality::GetagreementChi2Cut ( Int_t  patNR )

inline

                                                       {
        return eAgreementChi2WDistCut[patNR];
    }

GetagreementChi2CutMeanChi2()

Float_t EdbShowPVRQuality::GetagreementChi2CutMeanChi2 ( )

inline

                                                    {
        return eAgreementChi2CutMeanChi2;
    }

GetagreementChi2CutMeanW()

Float_t EdbShowPVRQuality::GetagreementChi2CutMeanW ( )

inline

                                                 {
        return eAgreementChi2CutMeanW;
    }

GetagreementChi2CutRMSChi2()

Float_t EdbShowPVRQuality::GetagreementChi2CutRMSChi2 ( )

inline

                                                   {
        return eAgreementChi2CutRMSChi2;
    }

GetagreementChi2CutRMSW()

Float_t EdbShowPVRQuality::GetagreementChi2CutRMSW ( )

inline

                                                {
        return eAgreementChi2CutRMSW;
    }

GetBTDensityLevel()

Float_t EdbShowPVRQuality::GetBTDensityLevel ( )

inline

                                         {
        return eBTDensityLevel;
    }

GetBTDensityLevelCalcMethodMC()

Bool_t EdbShowPVRQuality::GetBTDensityLevelCalcMethodMC ( )

inline

                                                  {
        return eBTDensityLevelCalcMethodMC;
    }

GetBTDensityLevelCalcMethodMCConfirmation()

Int_t EdbShowPVRQuality::GetBTDensityLevelCalcMethodMCConfirmation ( )

inline

                                                             {
        return eBTDensityLevelCalcMethodMCConfirmationNumber;
    }

GetCutMethod()

Int_t EdbShowPVRQuality::GetCutMethod ( )

inline

                                    {
        return eCutMethod;
    }

GetCutMethodIsDone()

Bool_t EdbShowPVRQuality::GetCutMethodIsDone ( Int_t  type )

inline

                                                    {
        if (type>1) return 0;
        return eCutMethodIsDone[type];
    }

GetCutp0() [1/2]

Float_t * EdbShowPVRQuality::GetCutp0 ( )

inline

                                 {
        return eCutp0;
    }

GetCutp0() [2/2]

Float_t EdbShowPVRQuality::GetCutp0 ( Int_t  patNR )

inline

                                            {
        return eCutp0[patNR];
    }

GetCutp1() [1/2]

Float_t * EdbShowPVRQuality::GetCutp1 ( )

inline

                                 {
        return eCutp1;
    }

GetCutp1() [2/2]

Float_t EdbShowPVRQuality::GetCutp1 ( Int_t  patNR )

inline

                                            {
        return eCutp1[patNR];
    }

GetEdbPVRec() [1/3]

EdbPVRec * EdbShowPVRQuality::GetEdbPVRec ( )

inline

                                   {
        return GetEdbPVRec(eNeedModified);
    }

GetEdbPVRec() [2/3]

EdbPVRec * EdbShowPVRQuality::GetEdbPVRec ( Bool_t  NeedModified )

inline

                                                      {
        cout << "Inline EdbPVRecType= " <<  NeedModified << endl;
        if (NeedModified==1) {
            return eAli_modified;
        }
        else {
            return eAli_orig;
        }
    }

GetEdbPVRec() [3/3]

EdbPVRec * EdbShowPVRQuality::GetEdbPVRec ( Int_t  EdbPVRecType )

inline

                                                     {
        cout << "Inline EdbPVRecType= " <<  EdbPVRecType << endl;
        if (EdbPVRecType==1) {
            return eAli_modified;
        }
        else {
            return eAli_orig;
        }
    }

GetEdbPVRec_modified()

EdbPVRec * EdbShowPVRQuality::GetEdbPVRec_modified ( )

inline

                                            {
        return GetEdbPVRec(1);
    }

GetEdbPVRec_orig()

EdbPVRec * EdbShowPVRQuality::GetEdbPVRec_orig ( )

inline

                                        {
        return GetEdbPVRec(0);
    }

GetHistChi2W()

TH2F * EdbShowPVRQuality::GetHistChi2W ( )

inline

                                  {
        return eHistChi2W;
    }

GetHistYX()

TH2F * EdbShowPVRQuality::GetHistYX ( )

inline

                               {
        return eHistYX;
    }

GetTracksFromLinkedTracksRootFile()

TObjArray * EdbShowPVRQuality::GetTracksFromLinkedTracksRootFile

(

)

{
    cout << "EdbShowPVRQuality::GetTracksFromLinkedTracksRootFile()" << endl;
// cout << "___TODO___    STEP  1 :  Check if linked_tracks.root exists   " << endl;
// cout << "___TODO___    STEP  2 :  Open linked_tracks   " << endl;
// cout << "___TODO___    STEP  3 :  Check if there are tracks inside   " << endl;
// cout << "___TODO___    STEP  4 :  Get Tracks as this TObjArray  " << endl;
 
    // Array containing EdbTrackP objects:
    TObjArray* trackarray = new TObjArray();
 
    TFile*  file = new TFile("linked_tracks.root");
    TTree* tr= (TTree*)file->Get("tracks");
    TClonesArray *seg= new TClonesArray("EdbSegP",60);
    EdbSegP* segment=0;
    EdbSegP*    s2=0;
    EdbSegP*    t=0;
 
    int nentr = int(tr->GetEntries());
    //   check if tracks has entries: if so then ok, otherwise return directly:
    if (nentr>0) {
        cout << "EdbShowPVRQuality::GetTracksFromLinkedTracksRootFile()   " << nentr << "  entries Total"<<endl;
    }
    else {
        cout << "EdbShowPVRQuality::GetTracksFromLinkedTracksRootFile()   linked_tracks.root file has  NO  entries...Return NULL." << endl;
        return NULL;
    }
 
    int nseg,n0,npl;
    tr->SetBranchAddress("t."  , &t  );
    tr->SetBranchAddress("s"  , &seg  );
    tr->SetBranchAddress("nseg"  , &nseg  );
    tr->SetBranchAddress("n0"  , &n0  );
    tr->SetBranchAddress("npl"  , &npl  );
 
    for (int i=0; i<nentr; i++ ) {
        tr->GetEntry(i);
        EdbTrackP* realTrack = new EdbTrackP(t);
        for (int j=0; j<nseg; j++ ) {
            s2=(EdbSegP*) seg->At(j);
            segment=(EdbSegP*)s2->Clone();
            realTrack->AddSegment(segment);
        }
        //realTrack->PrintNice();
        realTrack->SetNpl();    // Necessary otherwise these variables are not filled.
        realTrack->SetN0();     // Necessary otherwise these variables are not filled.
        trackarray->Add(realTrack);
    }
    delete tr;
    file->Close();
    delete file;
    cout << "EdbShowPVRQuality::GetTracksFromLinkedTracksRootFile()...done." << endl;
    return trackarray;
}

Help()

void EdbShowPVRQuality::Help

(

)

{
    // Basic Help Function.
 
    cout << "----------------------------------------------" << endl;
    cout << "-----     void EdbShowPVRQuality::Help()     -----" << endl;
    cout << "-----" << endl;
    cout << "-----     This Class helps you to determine the Quality Cut Plate by Plate" << endl;
    cout << "-----     for suited BG level for shower reco." << endl;
    cout << "-----     You use it like this:" << endl;
    cout << "-----        EdbShowPVRQuality* QualityClass = new EdbShowPVRQuality(gAli)" << endl;
    cout << "-----     where  gAli is an EdbPVRec object pointer (usually the one" << endl;
    cout << "-----     from the cp files). Then you can get the CutValues with" << endl;
    cout << "-----     GetCutp0() (for each//all plate) back." << endl;
    cout << "-----" << endl;
    cout << "-----     On TODO list: to be interfaced with the libShower reconstruction mode." << endl;
    cout << "-----" << endl;
    cout << "-----" << endl;
    cout << "-----" << endl;
    cout << "-----  " << endl;
 
    cout << "-----     void EdbShowPVRQuality::Help()     -----" << endl;
    cout << "----------------------------------------------" << endl;
}

Init()

void EdbShowPVRQuality::Init ( void  )

protected

TEMPORARY, later the Standard Geometry should default be OPERA. For now we use a very large histogram that can contain both case x-y ranges. This takes slightly more memory but it shouldnt matter.

{
    // Basic Init function that creates objects in the memory which have
    // to be created only ONE time per class instance.
 
    Log(2,"EdbShowPVRQuality::Init","EdbShowPVRQuality::Init");
 
    eProfileBTdens_vs_PID_source = new TProfile("eProfileBTdens_vs_PID_source","eProfileBTdens_vs_PID_source",114,-0.5,113.5,0,200);
    eProfileBTdens_vs_PID_target = new TProfile("eProfileBTdens_vs_PID_target","eProfileBTdens_vs_PID_target",114,-0.5,113.5,0,200);
 
    eHistChi2W = new TH2F("eHistChi2W","eHistChi2W",40,0,40,90,0,3);
    eHistYX = new TH2F("eHistYX","eHistYX",100,0,1,100,0,1);
 
//     if (eHistGeometry==0) SetHistGeometry_OPERA();
//     cout << "EdbShowPVRQuality::Init()   /// TEMPORARY  SetHistGeometry_OPERAandMC " <<  endl;
    SetHistGeometry_OPERAandMC();
 
 
    Log(2,"EdbShowPVRQuality::Init","EdbShowPVRQuality::Init...done.");
    return;
}

Print()

void EdbShowPVRQuality::Print

(

)

{
    Log(2,"EdbShowPVRQuality::Print","EdbShowPVRQuality::Print");
    for (int i=0; i<80; ++i) cout << "-";
    cout << endl;
    cout << "-";
    cout << endl;
    cout << "EdbShowPVRQuality::Print()" << endl;
    cout << "-";
    cout << endl;
    for (int i=0; i<80; ++i) cout << "-";
    cout << endl;
    cout << "-" << endl;
    cout << "EdbShowPVRQuality::Print() --- GENERAL SETTINGS: ---" << endl;
    cout << "EdbShowPVRQuality::Print() --- " << setw(40) << "eCutMethod = " << eCutMethod << endl;
    cout << "EdbShowPVRQuality::Print() --- " << setw(40) << "eCutMethodIsDone[0] = " << eCutMethodIsDone[0] << endl;
    cout << "EdbShowPVRQuality::Print() --- " << setw(40) << "eCutMethodIsDone[1] = " << eCutMethodIsDone[1] << endl;
    cout << "EdbShowPVRQuality::Print() --- " << setw(40) << "eBTDensityLevel = " << eBTDensityLevel << endl;
    cout << "EdbShowPVRQuality::Print() --- " << setw(40) << "eBTDensityLevelCalcMethodMC = " << eBTDensityLevelCalcMethodMC << endl;
    cout << "EdbShowPVRQuality::Print() --- " << setw(40) << "eBTDensityLevelCalcMethodMCConfirmationNumber = " << eBTDensityLevelCalcMethodMCConfirmationNumber << endl;
    cout << "EdbShowPVRQuality::Print() --- " << setw(40) << "eHistGeometry = " << eHistGeometry <<     endl;
    cout << "-" << endl;
    cout << "EdbShowPVRQuality::Print() --- SETTINGS: Input data:" << endl;
    cout << "EdbShowPVRQuality::Print() --- " << setw(40) << "eAli_orig = " << eAli_orig << endl;
    cout << "EdbShowPVRQuality::Print() --- " << setw(40) << "eAli_maxNpatterns = " << eAli_maxNpatterns << endl;
    cout << "EdbShowPVRQuality::Print() --- " << setw(40) << "eIsSource = " << eIsSource << endl;
    cout << "EdbShowPVRQuality::Print() --- " << setw(40) << "eProfileBTdens_vs_PID_source_meanX = " << eProfileBTdens_vs_PID_source_meanX << endl;
    cout << "EdbShowPVRQuality::Print() --- " << setw(40) << "eProfileBTdens_vs_PID_source_meanY = " << eProfileBTdens_vs_PID_source_meanY << endl;
    cout << "EdbShowPVRQuality::Print() --- " << setw(40) << "eProfileBTdens_vs_PID_source_rmsX = " << eProfileBTdens_vs_PID_source_rmsX << endl;
    cout << "EdbShowPVRQuality::Print() --- " << setw(40) << "eProfileBTdens_vs_PID_source_rmsY = " << eProfileBTdens_vs_PID_source_rmsY << endl;
    cout << "-" << endl;
    cout << "EdbShowPVRQuality::Print() --- SETTINGS: Output data:" << endl;
    cout << "EdbShowPVRQuality::Print() --- " << setw(40) << "eAli_modified = " << eAli_modified << endl;
    cout << "EdbShowPVRQuality::Print() --- " << setw(40) << "eAli_maxNpatterns = " << eAli_maxNpatterns << endl;
    cout << "EdbShowPVRQuality::Print() --- " << setw(40) << "eIsTarget = " << eIsTarget << endl;
    cout << "EdbShowPVRQuality::Print() --- " << setw(40) << "eProfileBTdens_vs_PID_target_meanX = " << eProfileBTdens_vs_PID_target_meanX << endl;
    cout << "EdbShowPVRQuality::Print() --- " << setw(40) << "eProfileBTdens_vs_PID_target_meanY = " << eProfileBTdens_vs_PID_target_meanY << endl;
    cout << "EdbShowPVRQuality::Print() --- " << setw(40) << "eProfileBTdens_vs_PID_target_rmsX = " << eProfileBTdens_vs_PID_target_rmsX << endl;
    cout << "EdbShowPVRQuality::Print() --- " << setw(40) << "eProfileBTdens_vs_PID_target_rmsY = " << eProfileBTdens_vs_PID_target_rmsY << endl;
    cout << "EdbShowPVRQuality::Print() --- " << setw(40) << "eNeedModified = " << eNeedModified << endl;
    cout << "-" << endl;
    cout << "EdbShowPVRQuality::Print()....done." << endl;
    for (int i=0; i<80; ++i) cout << "-";
    cout << endl;
 
    Log(2,"EdbShowPVRQuality::Print","EdbShowPVRQuality::Print...done.");
    return;
}

PrintCutType()

void EdbShowPVRQuality::PrintCutType

(

)

{
    // Print PrintCutType data for this class.
    cout << "EdbShowPVRQuality::Print()   " << endl;
    cout << "EdbShowPVRQuality::Print()   eCutMethodIsDone[0] " << eCutMethodIsDone[0] << endl;
    cout << "EdbShowPVRQuality::Print()   eCutMethodIsDone[1] " << eCutMethodIsDone[1] << endl;
    if (eCutMethodIsDone[0]) PrintCutType0();
    if (eCutMethodIsDone[1]) PrintCutType1();
    return;
}

PrintCutType0()

void EdbShowPVRQuality::PrintCutType0

(

)

{
    //     Prints quality cut values for each plate of the original EdbPVRec object that were
    //     applied to achieve the basetrack density level, after application the specific cut
    //     on the segments of the specific plate (=pattern).
 
    if (!eIsSource) return;
 
    if (!eCutMethodIsDone[0]) {
        cout << "---   EdbShowPVRQuality::PrintCutType0()----------" << endl;
        cout << "---      This CutType has not been done. Run  Execute_ConstantBTDensity() ---" << endl;
        cout << "---      to see the results (return now). ---" << endl;
        return;
    }
 
    cout << "----------void EdbShowPVRQuality::PrintCutType0()----------" << endl;
    cout << "Pattern || Z() || Nseg || eCutp0[i] || eCutp1[i] || BTDensity_orig ||...|| Nseg_modified || BTDensity_modified ||"<< endl;
 
    if (NULL==eAli_modified) {
        eAli_modified=eAli_orig;
        cout << "WARNING eAli_modified==NULL  ==>>  Take eAli_orig instead. To build eAli_modified please run CreateEdbPVRec(eAli_orig)" << endl;
    }
 
    int Npat_orig = eAli_orig->Npatterns();
    for (int i=0; i<Npat_orig; i++) {
        EdbPattern* pat_orig = (EdbPattern*)eAli_orig->GetPattern(i);
        Int_t npatO=pat_orig->N();
        EdbPattern* pat_modified = (EdbPattern*)eAli_modified->GetPattern(i);
        Int_t npatM=pat_modified->N();
        //
        cout << i;
        cout << "   ";
        printf("%.1f  %d  %.3f  %.3f  %.2f",pat_orig->Z(),npatO, eCutp0[i], eCutp1[i] , ePatternBTDensity_orig[i]);
        cout << "   ... ";
        printf("%.1f  %d  %.3f  %.3f  %.2f",pat_modified->Z(),npatM,  eCutp0[i] ,eCutp1[i],  ePatternBTDensity_modified[i]);
        cout << endl;
    }
 
    return;
}

PrintCutType1()

void EdbShowPVRQuality::PrintCutType1

(

)

{
    //     Prints quality cut values for each plate of the original EdbPVRec object that were
    //     applied to achieve the basetrack density level, after application the specific cut
    //     on the segments of the specific plate (=pattern).
 
    if (!eIsSource) return;
    if (!eCutMethodIsDone[1]) {
        cout << "---   EdbShowPVRQuality::PrintCutType1()----------" << endl;
        cout << "---      This CutType has not been done. Run  Execute_ConstantQuality() ---" << endl;
        cout << "---      to see the results (return now). ---" << endl;
        return;
    }
 
    cout << "----------void EdbShowPVRQuality::PrintCutType1()----------" << endl;
 
    cout << "Pattern || Z() || Nseg || BTDensity_orig || Chi2CutMeanChi2 || Chi2CutRMSChi2 || Chi2CutMeanW || Chi2CutRMSW || Chi2Cut[i] || BTDensity_modified ||"<< endl;
 
    if (NULL==eAli_modified) {
        eAli_modified=eAli_orig;
        cout << "WARNING eAli_modified==NULL  ==>>  Take eAli_orig instead. To build eAli_modified please run CreateEdbPVRec(eAli_orig)" << endl;
    }
 
    int Npat_orig = eAli_orig->Npatterns();
    for (int i=0; i<Npat_orig; i++) {
        EdbPattern* pat_orig = (EdbPattern*)eAli_orig->GetPattern(i);
        Int_t npatO=pat_orig->N();
        EdbPattern* pat_modified = (EdbPattern*)eAli_modified->GetPattern(i);
        Int_t npatM=pat_modified->N();
        cout << i;
        cout << "   ";
        printf("%.1f  %d  %.2f  %.2f  %.2f  %.2f  %.2f  %.2f",pat_orig->Z(),npatO, eAgreementChi2CutMeanChi2 , eAgreementChi2CutRMSChi2,  eAgreementChi2CutMeanW , eAgreementChi2CutRMSW, eAgreementChi2WDistCut[i], ePatternBTDensity_orig[i]);
        cout << "   ... ";
        printf("%.1f  %d  %.2f  %.2f  %.2f  %.2f  %.2f  %.2f",pat_modified->Z(),npatM, eAgreementChi2CutMeanChi2 , eAgreementChi2CutRMSChi2,  eAgreementChi2CutMeanW , eAgreementChi2CutRMSW, eAgreementChi2WDistCut[i], ePatternBTDensity_modified[i]);
        cout << endl;
    }
    return;
}

Remove_DoubleBT()

EdbPVRec * EdbShowPVRQuality::Remove_DoubleBT

(

EdbPVRec *

aliSource

)

{
    // Remove double counted basetracks. This is a scanning induced error, at the
    // edges of the objective field, corrections are different, thus one
    // basetrack is seen differently in two different views.
    // In Bern data, we do see clearly this double bump, that is characterized by
    // _very_ close values of two BTs in X,Y,TX,TY (and also Chi2 and W).
    // Separation threshold values are 2microns in postion and 10 mrad in angle.
    // (again obtained from Bern data).
 
    // Quick and Dirty implementation !
    cout << "EdbShowPVRQuality::Remove_DoubleBT()" << endl;
    cout << "EdbShowPVRQuality::Remove_DoubleBT()  aliSource = " << aliSource << endl;
    cout << "EdbShowPVRQuality::Remove_DoubleBT()  eAli_orig = " << eAli_orig << endl;
 
    if (NULL==aliSource) {
        cout << "WARNING!----EdbShowPVRQuality::Remove_DoubleBT()  Source EdbPVRec is NULL. Try to change to object eAli_orig: " << eAli_orig << endl;
        if (NULL==eAli_orig) {
            cout << "WARNING!----EdbShowPVRQuality::Remove_DoubleBT() Also eAli_orig EdbPVRec is NULL. Do nothing and return NULL pointer!" << endl;
            return NULL;
        }
        else {
            aliSource=eAli_orig;
        }
    }
 
    // Make a new PVRec object anyway
    EdbPVRec* aliTarget = new EdbPVRec();
    aliTarget->Print();
 
    Bool_t seg_seg_close=kFALSE;
    EdbSegP* seg=0;
    EdbSegP* seg1=0;
    Int_t NdoubleFoundSeg=0;
 
    //gEDBDEBUGLEVEL=3;
    cout << "EdbShowPVRQuality::Remove_DoubleBT() Start looping now. Attention, this might take a while!" << endl;
 
    for (int i = 0; i <aliSource->Npatterns(); i++ ) {
        if (gEDBDEBUGLEVEL==2) cout << "." << flush;
        if (gEDBDEBUGLEVEL>1) cout << "Looping over Ali_source->Pat()=" << i <<  ". Until now double Candidates found: " << NdoubleFoundSeg << endl;
 
        // Create Target Pattern:
        EdbPattern* pat = aliSource->GetPattern(i);
        EdbPattern* pt= new EdbPattern();
        // SetPattern Values to the parent patterns:
        pt->SetID(pat->ID());
        pt->SetPID(pat->PID());
        pt->SetZ(pat->Z());
        // Helper Variable for cout purpose
        Int_t nPat=pat->N();
 
        for (int j = 0; j < nPat-1; j++ ) {
            if (gEDBDEBUGLEVEL>2) if (j%(nPat/10)==0) cout << "10%more of loop done." << flush;
 
            seg = pat->GetSegment(j);
            seg_seg_close=kFALSE;
 
            for (int k = j+1; k <pat->N(); k++ ) {
                if (seg_seg_close) continue;
 
                if (gEDBDEBUGLEVEL>3) cout << "Looping over pattern for segment pair nr=" << j << "," << k << endl;
                seg1 = pat->GetSegment(k);
 
                // Here decide f.e. which segments to check...
                Bool_t toContinue=kFALSE;
 
                // Skip already large distances
                if (TMath::Abs(seg->X()-seg1->X())>20.1) continue;
                if (TMath::Abs(seg->Y()-seg1->Y())>20.1) continue;
                if (TMath::Abs(seg->TX()-seg1->TX())>0.1) continue;
                if (TMath::Abs(seg->TY()-seg1->TY())>0.1) continue;
 
                //
                // a) all BT parings within a square DeltaR and Delta Theta:
                //
                // b) all BT pairings within dT<->dR line:
                //
                Float_t dR_allowed=0;
                Float_t deltaR=0;
                Float_t deltaTheta=0;
                Float_t dR_allowed_m_deltaR=0;
 
                deltaTheta=TMath::Sqrt( TMath::Power(seg->TX()-seg1->TX(),2)+TMath::Power(seg->TY()-seg1->TY(),2) );
                if (deltaTheta>0.1) continue;
 
                deltaR=TMath::Sqrt( TMath::Power(seg->X()-seg1->X(),2)+TMath::Power(seg->Y()-seg1->Y(),2) );
 
                // line defined by experimental values:
                // dR_allowed = 11(micron)/0.1(rad) * dTheta
                dR_allowed = 10.75/0.1*deltaTheta;
                dR_allowed_m_deltaR=TMath::Abs(deltaR-dR_allowed);
 
                if (deltaR>3.0||deltaTheta>0.005) toContinue=kTRUE;
                if (dR_allowed_m_deltaR<1.0) toContinue=kFALSE;
                if (toContinue) continue;
 
                if (gEDBDEBUGLEVEL>3) cout << "EdbPVRQuality::Remove_DoubleBT()   Found incompatible segment for dR_dT line condition: " << endl;
 
                // if (gEDBDEBUGLEVEL>3) cout << "EdbPVRQuality::Remove_DoubleBT()   Do last check if both are MCEvt segments of different event number! " << endl;
                // if ((seg->MCEvt()!=seg1->MCEvt())&&seg->MCEvt()<0) continue;
                // I have doubts if I shall take MC events also out, but I guess Yes, because if
                // in data these small pairings appear, then they will fall also under the
                // category "fake doublet" and then be removed, even if they are real double BT
                // from a signal.
 
 
                ++NdoubleFoundSeg;
                seg_seg_close=kTRUE;
                //                 if (seg_seg_close) break;
            }
            if (seg_seg_close) continue;
 
            // Add segment:
            if (gEDBDEBUGLEVEL>3) cout << "// Add segment:" << endl;
            pt->AddSegment(*seg);
        }
        if (gEDBDEBUGLEVEL>2) cout << "// Add Pattern:" << endl;
        aliTarget->AddPattern(pt);
    }
 
    if (gEDBDEBUGLEVEL>1) aliSource->Print();
    if (gEDBDEBUGLEVEL>1) aliTarget->Print();
 
    cout << "EdbShowPVRQuality::Remove_DoubleBT() Statistics: " << endl;
    cout << "EdbShowPVRQuality::Remove_DoubleBT() We found " << NdoubleFoundSeg  << " double segments too close to each other to be different BTracks." << endl;
    cout << "EdbShowPVRQuality::Remove_DoubleBT()...done." << endl;
    return aliTarget;
}

Remove_Passing()

EdbPVRec * EdbShowPVRQuality::Remove_Passing

(

EdbPVRec *

aliSource

)

{
    // Removes Passing Tracks from the EdbPVRec source object.
    // Unfortunately, there does Not Exist an easy function like
    // ->RemoveSegment from EdbPVRec object. That makes implementation complicated.
 
    cout << "EdbShowPVRQuality::Remove_Passing()." << endl;
    cout << "EdbShowPVRQuality::Remove_Passing()  aliSource = " << aliSource << endl;
    cout << "EdbShowPVRQuality::Remove_Passing()  eAli_orig = " << eAli_orig << endl;
 
    if (NULL==aliSource) {
        cout << "WARNING!----EdbShowPVRQuality::Remove_Passing()  Source EdbPVRec is NULL. Try to change to object eAli_orig: " << eAli_orig << endl;
        if (NULL==eAli_orig) {
            cout << "WARNING!----EdbShowPVRQuality::Remove_Passing() Also eAli_orig EdbPVRec is NULL. Do nothing and return NULL pointer!" << endl;
            return NULL;
        }
        else {
            aliSource=eAli_orig;
        }
    }
 
    // Get the tracks from the source object:
    TObjArray* Tracks=aliSource->eTracks;
 
    // If eAli_source has no tracks, we return here and stop.
    if (NULL == Tracks) {
        cout << "WARNING!----EdbShowPVRQuality::Remove_Passing() NULL == eTracks." << endl;
        cout << "EdbShowPVRQuality::Remove_Passing() Read tracks from a linked_tracks.root file if there is any." << endl;
        Tracks=GetTracksFromLinkedTracksRootFile();
 
        if (NULL == Tracks) {
            cout << "WARNING!----EdbShowPVRQuality::GetTracksFromLinkedTracksRootFile() NULL == eTracks." << endl;
            cout << "WARNING!----EdbShowPVRQuality::Remove_Passing() Reading tracks from a linked_tracks.root file failed." << endl;
            cout << "WARNING!----EdbShowPVRQuality::Remove_Passing() Return the original object back." << endl;
            return aliSource;
        }
    }
 
    // Now the object array  Tracks  should be filled anyway, either from the source EdbPVR
    // or from the linked_tracks.root file:
    Int_t TracksN=Tracks->GetEntries();
    EdbTrackP* track;
    EdbSegP* trackseg;
 
    if (gEDBDEBUGLEVEL>1) {
        cout << "EdbShowPVRQuality::Remove_Passing()  aliSource/linked_tracks.root has Tracks N=" << TracksN << endl;
    }
 
    // Make a new PVRec object anyway
    EdbPVRec* aliTarget = new EdbPVRec();
    Bool_t seg_in_eTracks=kFALSE;
    Int_t totallyRemovedSegments=0;
    Int_t totallyAddedSegments=0;
    Int_t aliSourceNPat=aliSource->Npatterns();
 
    for (int i = 0; i < aliSourceNPat; i++ ) {
        if (gEDBDEBUGLEVEL>2) cout << "Looping over aliSource->Pat()=" << i << endl;
 
        EdbPattern* pat = aliSource->GetPattern(i);
        EdbPattern* pt= new EdbPattern();
        // SetPattern Values to the parent patterns:
        pt->SetID(pat->ID());
        pt->SetPID(pat->PID());
        pt->SetZ(pat->Z());
 
        for (int j = 0; j <pat->N(); j++ ) {
            EdbSegP* seg = pat->GetSegment(j);
            seg_in_eTracks=kFALSE;
 
            if (gEDBDEBUGLEVEL>3) cout << "Checking segment j= " << j << " and loop over tracks to check correspondance." <<endl;
 
            for (int ll = 0; ll<TracksN; ll++ ) {
                track=(EdbTrackP*)Tracks->At(ll);
                Int_t TrackSegN=track->N();
 
                // Compare if this track is regarded as "passing"
                // when its less than5  plates than the original pvr object
                if (track->Npl()<aliSourceNPat-5) continue;
//                                  cout << track->Npl() << " " << aliSourceNPat << endl;
 
                // Here decide f.e. which tracks to check...
                // On cosmics it would be nice that f.e. Npl()>=Npatterns-4 ...
                // if ....()....
                // Since addresses of objects can vary, its better to compare them
                // by absolute positions.
 
                for (int kk = 0; kk<TrackSegN; kk++ ) {
                    if (seg_in_eTracks) continue;
                    trackseg=track->GetSegment(kk);
                    if (TMath::Abs(seg->Z()-trackseg->Z())>10.1) continue;
                    if (TMath::Abs(seg->X()-trackseg->X())>5.1) continue;
                    if (TMath::Abs(seg->Y()-trackseg->Y())>5.1) continue;
                    if (TMath::Abs(seg->TX()-trackseg->TX())>0.05) continue;
                    if (TMath::Abs(seg->TY()-trackseg->TY())>0.05) continue;
                    //cout << "Found compatible segment!! " << endl;
                    totallyRemovedSegments++;
                    seg_in_eTracks=kTRUE;
                }
                if (seg_in_eTracks) continue;
            }
            if (seg_in_eTracks) continue;
            // Arrived here then the segmen has no equivalent in any segments
            // of the tracks array.
 
            // Add segment:
            if (gEDBDEBUGLEVEL>3) cout << "EdbShowPVRQuality::Remove_Passing() Add segment:" << endl;
            pt->AddSegment(*seg);
            totallyAddedSegments++;
        }
        if (gEDBDEBUGLEVEL>2) cout << "EdbShowPVRQuality::Remove_Passing()  TotallyAddedSegments (up to now) = " << totallyAddedSegments << " // Add Pattern:" << endl;
        aliTarget->AddPattern(pt);
    }
 
    if (gEDBDEBUGLEVEL>1) aliSource->Print();
    if (gEDBDEBUGLEVEL>1) aliTarget->Print();
 
 
 
    cout << "EdbShowPVRQuality::Remove_Passing() Totally removed segments= " << totallyRemovedSegments << endl;
    cout << "EdbShowPVRQuality::Remove_Passing()...done." << endl;
    return aliTarget;
}

Remove_Segment()

EdbPVRec * EdbShowPVRQuality::Remove_Segment	(	EdbSegP *	seg,
		EdbPVRec *	aliSource = NULL,
		Int_t	Option = 0
	)

                                                                                           {
    // Remove one segment (EdbSegP* seg type).
    // Return a new EdbPVRec* with previous element removed.
    cout << "EdbShowPVRQuality::Remove_Segment():" << endl;
    TObjArray* segArray= new TObjArray();
    segArray->Add(seg);
    Remove_SegmentArray(segArray,aliSource,Option);
    delete segArray;
    cout << "EdbShowPVRQuality::Remove_Segment()...done." << endl;
    return NULL;
}

Remove_SegmentArray()

EdbPVRec * EdbShowPVRQuality::Remove_SegmentArray	(	TObjArray *	segarray,
		EdbPVRec *	aliSource = NULL,
		Int_t	Option = 0
	)

{
    // Every "Remove_XXY()" function will call Remove_SegmentArray() since the removal of
    // a segment array is easiest to implement. As said earlier, there is no easy function
    // to remove tracks from a EdbPVRec object directly. So instead we have to create a
    // new EdbPVRec object, fill with all basetracks from the old one, except where the
    // segments of the SegmentArray coincide with those which are in the source EdbPVRec.
 
    cout << "EdbShowPVRQuality::Remove_SegmentArray():" << endl;
 
    if (NULL==aliSource) {
        cout << "EdbShowPVRQuality::Remove_SegmentArray()   ERROR   aliSource=NULL pointer. I cannot" << endl;
        cout << "EdbShowPVRQuality::Remove_SegmentArray()   ERROR   do something here and I will return a NULL pointer now. Stop." << endl;
        return NULL;
    }
 
    if (gEDBDEBUGLEVEL>2) {
        cout << "EdbShowPVRQuality::Remove_SegmentArray()   INFO   " << endl;
        cout << "EdbShowPVRQuality::Remove_SegmentArray()   segarray =    " << segarray << " with " << segarray->GetEntries() << " entries." << endl;
        cout << "EdbShowPVRQuality::Remove_SegmentArray()   aliSource =    " << aliSource << " with " << aliSource->Npatterns() << " patterns." << endl;
        cout << "EdbShowPVRQuality::Remove_SegmentArray()   Option =    " << Option << "." << endl;
    }
 
 
    cout << "EdbShowPVRQuality::Remove_SegmentArray()...done." << endl;
    return NULL;
}

Remove_Track()

EdbPVRec * EdbShowPVRQuality::Remove_Track	(	EdbTrackP *	track,
		EdbPVRec *	aliSource = NULL,
		Int_t	Option = 0
	)

                                                                                             {
    // Remove one track (EdbTrackP* track type).
    // Return a new EdbPVRec* with previous element removed.
    cout << "EdbShowPVRQuality::Remove_Track():" << endl;
    TObjArray* segArray= new TObjArray();
    segArray=  TrackToSegmentArray(track);
    Remove_SegmentArray(segArray,aliSource,Option);
    delete segArray;
    cout << "EdbShowPVRQuality::Remove_Track()...done." << endl;
    return NULL;
}

Remove_TrackArray()

EdbPVRec * EdbShowPVRQuality::Remove_TrackArray	(	TObjArray *	trackArray,
		EdbPVRec *	aliSource = NULL,
		Int_t	Option = 0
	)

                                                                                                       {
    // Remove a whole track array (Array of EdbTrackP type).
    // Return a new EdbPVRec* with previous element removed.
    cout << "EdbShowPVRQuality::Remove_TrackArray():" << endl;
    TObjArray* segArray= new TObjArray();
    segArray=  TrackArrayToSegmentArray(trackArray);
    Remove_SegmentArray(segArray,aliSource,Option);
    delete segArray;
    cout << "EdbShowPVRQuality::Remove_TrackArray()...done." << endl;
    return NULL;
}

Set0()

void EdbShowPVRQuality::Set0 ( )

protected

{
    // Reset Values
 
    Log(2,"EdbShowPVRQuality::Set0","EdbShowPVRQuality::Set0");
 
    eAli_orig=NULL;
    eAli_modified=NULL;
    eNeedModified=kFALSE;
    eIsSource=kFALSE;
    eIsTarget=kFALSE;
    eAli_maxNpatterns=0;
    for (int i=0; i<2; i++) eCutMethodIsDone[i]=kFALSE;
 
    // Default BT density level for which the standard cutroutine
    // will be put:
    eBTDensityLevel=20; // #BT/mm2
 
    // Default BT density level will use only segments for
    // data calculation .
    eBTDensityLevelCalcMethodMC=kFALSE;
    eBTDensityLevelCalcMethodMCConfirmationNumber=0;
 
    // Reset Default Geometry: 0 OperaGeometry, 1: MC Geometry
    eHistGeometry=0;
 
 
 
    eHistYX->Reset();
    eHistChi2W->Reset();
 
    for (int i=0; i<114; i++) {
        ePatternBTDensity_orig[i]=0;
        ePatternBTDensity_modified[i]=0;
        eCutp1[i]=0.15;
        eCutp0[i]=1.0; // Maximum Cut Value for const, BT dens
        eAgreementChi2WDistCut[i]=3.0;  // Maximum Cut Value for const, BT dens
    }
 
 
    eProfileBTdens_vs_PID_source_meanX=0;
    eProfileBTdens_vs_PID_source_meanY=0;
    eProfileBTdens_vs_PID_source_rmsX=0;
    eProfileBTdens_vs_PID_source_rmsY=0;
    eProfileBTdens_vs_PID_target_meanX=0;
    eProfileBTdens_vs_PID_target_meanY=0;
    eProfileBTdens_vs_PID_target_rmsX=0;
    eProfileBTdens_vs_PID_target_rmsY=0;
 
    // Default values for cosmics, taken from a brick data:
    // (I cant remember which one, I hope it doesnt matter).
    eAgreementChi2CutMeanChi2=1.0;
    eAgreementChi2CutRMSChi2=0.3;
    eAgreementChi2CutMeanW=23;
    eAgreementChi2CutRMSW=3;
 
    Log(2,"EdbShowPVRQuality::Set0","EdbShowPVRQuality::Set0...done.");
    return;
}

SetBTDensityLevel()

void EdbShowPVRQuality::SetBTDensityLevel ( Float_t  BTDensityLevel )

inline

                                                          {
        eBTDensityLevel=BTDensityLevel;
    }

SetBTDensityLevelCalcMethodMC()

void EdbShowPVRQuality::SetBTDensityLevelCalcMethodMC ( Bool_t  BTDensityLevelCalcMethodMC )

inline

                                                                                 {
        eBTDensityLevelCalcMethodMC=BTDensityLevelCalcMethodMC;
    }

SetBTDensityLevelCalcMethodMCConfirmation()

void EdbShowPVRQuality::SetBTDensityLevelCalcMethodMCConfirmation ( Int_t  BTDensityLevelCalcMethodMCConfirmationNumber )

inline

                                                                                                              {
        eBTDensityLevelCalcMethodMCConfirmationNumber=BTDensityLevelCalcMethodMCConfirmationNumber;
    }

SetCutMethod()

void EdbShowPVRQuality::SetCutMethod

(

Int_t

CutMethod

)

{
    // Set Cut Method of the background reduction:
    // 0: (default): Do cut based on the linear relation: seg.Chi2()<seg.eW*a-b
    // 1: (testing): Do cut based on a chi2-variable that compares with passing tracks
    // (if available), i.e. cosmics.
 
    Log(2,"EdbShowPVRQuality::SetCutMethod","EdbShowPVRQuality::SetCutMethod");
 
    eCutMethod=CutMethod;
    cout << "EdbShowPVRQuality::SetCutMethod  eCutMethod=  " << eCutMethod << endl;
    if (CutMethod>1) {
        eCutMethod=0;
        cout << "WARNING   EdbShowPVRQuality::SetCutMethod  eCutMethod invalid, Set back to default eCutMethod=  " << eCutMethod << endl;
    }
 
    Log(2,"EdbShowPVRQuality::SetCutMethod","EdbShowPVRQuality::SetCutMethod...done.");
    return;
}

SetEdbPVRec()

void EdbShowPVRQuality::SetEdbPVRec ( EdbPVRec *  Ali_orig )

inline

                                                            {
        eAli_orig=Ali_orig;
        eIsSource=kTRUE;
        eAli_maxNpatterns=Ali_orig->Npatterns();
    }

SetHistGeometry_MC()

void EdbShowPVRQuality::SetHistGeometry_MC

(

)

{
    // Set the geometry of the basetrack density evaluation using simulation case,
    // MC ORFEO size conventions: x=-xmax..0..+xmax;y=-ymax..0..ymax).
    // BinArea is 1mmx1mm
    eHistYX->Reset();
    eHistYX->SetBins(100,-50000,50000,100,-50000,50000);
    //cout << "SetHistGeometry_MC :binwidth (micron)= " << eHistYX->GetBinWidth(1) << endl;
    return;
}

SetHistGeometry_OPERA()

void EdbShowPVRQuality::SetHistGeometry_OPERA

(

)

{
    // Set the geometry of the basetrack density evaluation using OPERA case,
    // European Scanning System size conventions: x=0..125000;y=0..100000).
    // BinArea is 1mmx1mm.
    eHistYX->Reset();
    eHistYX->SetBins(100,0,100000,120,0,120000);
    //cout << "SetHistGeometry_OPERA :binwidth (micron)= " << eHistYX->GetBinWidth(1) << endl;
    return;
}

SetHistGeometry_OPERAandMC()

void EdbShowPVRQuality::SetHistGeometry_OPERAandMC

(

)

{
    // Set the geometry of the basetrack density evaluation covering MC and DATA case,
    // size conventions: x=-125000..0..+125000;y=-125000..0..125000).
    // BinArea is 1mmx1mm
 
    Log(3,"EdbShowPVRQuality::SetHistGeometry_OPERAandMC","EdbShowPVRQuality::SetHistGeometry_OPERAandMC");
    eHistYX->Reset();
    eHistYX->SetBins(250,-125000,125000,250,-125000,125000);
    if (gEDBDEBUGLEVEL>2) cout << "SetHistGeometry_OPERAandMC::binwidth (micron)= " << eHistYX->GetBinWidth(1) << endl;
    Log(3,"EdbShowPVRQuality::SetHistGeometry_OPERAandMC","EdbShowPVRQuality::SetHistGeometry_OPERAandMC...done");
    return;
}

TrackArrayToSegmentArray()

TObjArray * EdbShowPVRQuality::TrackArrayToSegmentArray

(

TObjArray *

trackArray

)

                                                                            {
    // Convert segments of tracks of a track array into an TObjArray containing segments.
    TObjArray* segArray= new TObjArray();
    Int_t nSegTotal=0;
    for (int j=0; j<trackArray->GetEntries(); j++) {
        EdbTrackP* track = (EdbTrackP*)trackArray->At(j);
        for (int i=0; i<track->N(); i++) {
            EdbSegP* seg = track->GetSegment(i);
            segArray->Add(seg);
            ++nSegTotal;
        }
    }
    if (gEDBDEBUGLEVEL>3) cout << "EdbShowPVRQuality::TrackArrayToSegmentArray() Totally added: " << nSegTotal << " segments from the track array."<<endl;
    return segArray;
}

TrackToSegmentArray()

TObjArray * EdbShowPVRQuality::TrackToSegmentArray

(

EdbTrackP *

track

)

                                                                  {
    // Convert segments of a track into an TObjArray containing segments.
    TObjArray* segArray= new TObjArray();
    Int_t nSegTotal=0;
    for (int i=0; i<track->N(); i++) {
        EdbSegP* seg = track->GetSegment(i);
        segArray->Add(seg);
        ++nSegTotal;
    }
    if (gEDBDEBUGLEVEL>3) cout << "EdbShowPVRQuality::TrackToSegmentArray() Totally added: " << nSegTotal << " segments from the track."<<endl;
    return segArray;
}

Member Data Documentation

eAgreementChi2CutMeanChi2

Float_t EdbShowPVRQuality::eAgreementChi2CutMeanChi2

private

eAgreementChi2CutMeanW

Float_t EdbShowPVRQuality::eAgreementChi2CutMeanW

private

eAgreementChi2CutRMSChi2

Float_t EdbShowPVRQuality::eAgreementChi2CutRMSChi2

private

eAgreementChi2CutRMSW

Float_t EdbShowPVRQuality::eAgreementChi2CutRMSW

private

eAgreementChi2WDistCut

Float_t EdbShowPVRQuality::eAgreementChi2WDistCut[114]

private

eAli_maxNpatterns

Int_t EdbShowPVRQuality::eAli_maxNpatterns

private

eAli_modified

EdbPVRec* EdbShowPVRQuality::eAli_modified

private

eAli_orig

EdbPVRec* EdbShowPVRQuality::eAli_orig

private

eBTDensityLevel

Float_t EdbShowPVRQuality::eBTDensityLevel

private

eBTDensityLevelCalcMethodMC

Bool_t EdbShowPVRQuality::eBTDensityLevelCalcMethodMC

private

eBTDensityLevelCalcMethodMCConfirmationNumber

Int_t EdbShowPVRQuality::eBTDensityLevelCalcMethodMCConfirmationNumber

private

eCutDistChi2

Float_t EdbShowPVRQuality::eCutDistChi2[114]

private

eCutDistW

Float_t EdbShowPVRQuality::eCutDistW[114]

private

eCutMethod

Int_t EdbShowPVRQuality::eCutMethod

private

eCutMethodIsDone

Bool_t EdbShowPVRQuality::eCutMethodIsDone[2]

private

eCutp0

Float_t EdbShowPVRQuality::eCutp0[114]

private

eCutp1

Float_t EdbShowPVRQuality::eCutp1[114]

private

eHistChi2W

TH2F* EdbShowPVRQuality::eHistChi2W

private

eHistGeometry

Int_t EdbShowPVRQuality::eHistGeometry

private

eHistYX

TH2F* EdbShowPVRQuality::eHistYX

private

eIsSource

Bool_t EdbShowPVRQuality::eIsSource

private

eIsTarget

Bool_t EdbShowPVRQuality::eIsTarget

private

eNeedModified

Bool_t EdbShowPVRQuality::eNeedModified

private

ePatternBTDensity_modified

Float_t EdbShowPVRQuality::ePatternBTDensity_modified[114]

private

ePatternBTDensity_orig

Float_t EdbShowPVRQuality::ePatternBTDensity_orig[114]

private

eProfileBTdens_vs_PID_source

TProfile* EdbShowPVRQuality::eProfileBTdens_vs_PID_source

private

eProfileBTdens_vs_PID_source_meanX

Float_t EdbShowPVRQuality::eProfileBTdens_vs_PID_source_meanX

private

eProfileBTdens_vs_PID_source_meanY

Float_t EdbShowPVRQuality::eProfileBTdens_vs_PID_source_meanY

private

eProfileBTdens_vs_PID_source_rmsX

Float_t EdbShowPVRQuality::eProfileBTdens_vs_PID_source_rmsX

private

eProfileBTdens_vs_PID_source_rmsY

Float_t EdbShowPVRQuality::eProfileBTdens_vs_PID_source_rmsY

private

eProfileBTdens_vs_PID_target

TProfile* EdbShowPVRQuality::eProfileBTdens_vs_PID_target

private

eProfileBTdens_vs_PID_target_meanX

Float_t EdbShowPVRQuality::eProfileBTdens_vs_PID_target_meanX

private

eProfileBTdens_vs_PID_target_meanY

Float_t EdbShowPVRQuality::eProfileBTdens_vs_PID_target_meanY

private

eProfileBTdens_vs_PID_target_rmsX

Float_t EdbShowPVRQuality::eProfileBTdens_vs_PID_target_rmsX

private

eProfileBTdens_vs_PID_target_rmsY

Float_t EdbShowPVRQuality::eProfileBTdens_vs_PID_target_rmsY

private

maxX

Float_t EdbShowPVRQuality::maxX

private

maxY

Float_t EdbShowPVRQuality::maxY

private

minX

Float_t EdbShowPVRQuality::minX

private

minY

Float_t EdbShowPVRQuality::minY

private

NbinsX

Int_t EdbShowPVRQuality::NbinsX

private

NbinsY

Int_t EdbShowPVRQuality::NbinsY

private

---

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

• /home/antonio/fedra_doxygen/src/libShowRec/EdbShowPVRQuality.h
• /home/antonio/fedra_doxygen/src/libShowRec/EdbShowPVRQuality.cxx