AngClusterFitAna.cxx File Reference

#include <iostream>
#include <algorithm>
#include <TCanvas.h>
#include <TRotation.h>
#include <TFile.h>
#include <TH2F.h>
#include <TH1F.h>
#include <TH3F.h>
#include <TPolyMarker3D.h>
#include <TPolyLine3D.h>
#include "StandardNtuple/NtpStRecord.h"
#include "CandNtupleSR/NtpSRRecord.h"
#include "CandNtupleSR/NtpSREvent.h"
#include "CandNtupleSR/NtpSRTrack.h"
#include "CandNtupleSR/NtpSRStrip.h"
#include "MessageService/MsgService.h"
#include "NueAna/AngClusterFitAna.h"
#include "AnalysisNtuples/ANtpDefaultValue.h"

Go to the source code of this file.

Functions

 CVSID ("$Id: AngClusterFitAna.cxx,v 1.21 2007/03/01 16:38:49 rhatcher Exp $")
static Double_t shwLongFunc (Double_t *x, Double_t *par)
static Double_t shwTransFunc (Double_t *x, Double_t *par)

Variables

const Float_t kDetectorPitch = 0.05949
const Float_t kStripWidth = 0.041
const Float_t kRadLengthToPlane = 1.54
const Float_t kMolRadDet = 0.0391
const Float_t kMeuToGeV = 1.0

Function Documentation

CVSID ( "$Id: AngClusterFitAna.  cxx,
v 1.21 2007/03/01 16:38:49 rhatcher Exp $"   
)
Id
AngClusterFitAna.cxx,v 1.21 2007/03/01 16:38:49 rhatcher Exp

class AngClusterFitAna

NueAna package

Purpose: Fit a EM shower profile to the 3D hits contained in the primary cluster determined by AngClusterAna..

Author: Alex Sousa <asousa@minos.phy.tufts.edu>

Created on: Fri May 06 2005

static Double_t shwLongFunc ( Double_t *  x,
Double_t *  par 
) [static]

Definition at line 135 of file AngClusterFitAna.cxx.

References kRadLengthToPlane.

Referenced by AngClusterFitAna::FitShower().

00136 {
00137 
00138     // Longitudinal profile to be fitted to em showers 
00139     //xx=Distance from the vertex in units of X0 
00140     //par[0]=a, par[1]=b, par[2]=E0
00141    
00142     Float_t xx=x[0];
00143     
00144     Double_t lnf = TMath::Log(kRadLengthToPlane*par[2]*par[1])+(par[0]-1)
00145         *TMath::Log(xx*par[1])-par[1]*xx-TMath::LnGamma(par[0]);
00146     
00147     Double_t f = TMath::Exp(lnf);
00148 
00149     return f;
00150 }

static Double_t shwTransFunc ( Double_t *  x,
Double_t *  par 
) [static]

Definition at line 152 of file AngClusterFitAna.cxx.

Referenced by AngClusterFitAna::FitShower().

00153 {
00154     //Transverse fit profile
00155     //x[0]=Radial distance from shower axis in units of Moliere radius
00156     //par[0]=Lambda1 attenuation length for high energy core
00157     //par[1]=Lambda2 attenuation length for shower halo
00158     //par[2]=Relative weight between lambda1 and lambda2
00159     //par[3]=E0 Energy normalization
00160 
00161     Float_t xx=x[0];
00162     Double_t f=par[3]*(TMath::Exp(-TMath::Sqrt(xx/par[0]))
00163                        +par[2]*TMath::Exp(-xx/par[1]));    
00164     return f;
00165 }


Variable Documentation

const Float_t kDetectorPitch = 0.05949

Definition at line 40 of file AngClusterFitAna.cxx.

const Float_t kMeuToGeV = 1.0

Definition at line 49 of file AngClusterFitAna.cxx.

const Float_t kMolRadDet = 0.0391

Definition at line 46 of file AngClusterFitAna.cxx.

Referenced by AngClusterFitAna::Analyze().

const Float_t kRadLengthToPlane = 1.54

Definition at line 44 of file AngClusterFitAna.cxx.

Referenced by AngClusterFitAna::Analyze(), and shwLongFunc().

const Float_t kStripWidth = 0.041

Definition at line 42 of file AngClusterFitAna.cxx.


Generated on 15 Jul 2018 for loon by  doxygen 1.6.1